Mechatronics and Robotics
Mechatronics and Robotics
Edited by:
Zoran Gacovski
ARCLER
P r e s s
www.arclerpress.com
Mechatronics and Robotics
Zoran Gacovski
Arcler Press
224 Shoreacres Road
Burlington, ON L7L 2H2
Canada
www.arclerpress.com
Email: [email protected]
This book contains information obtained from highly regarded resources. Reprinted
material sources are indicated. Copyright for individual articles remains with the au-
thors as indicated and published under Creative Commons License. A Wide variety of
references are listed. Reasonable efforts have been made to publish reliable data and
views articulated in the chapters are those of the individual contributors, and not neces-
sarily those of the editors or publishers. Editors or publishers are not responsible for
the accuracy of the information in the published chapters or consequences of their use.
The publisher assumes no responsibility for any damage or grievance to the persons or
property arising out of the use of any materials, instructions, methods or thoughts in the
book. The editors and the publisher have attempted to trace the copyright holders of all
material reproduced in this publication and apologize to copyright holders if permission
has not been obtained. If any copyright holder has not been acknowledged, please write
to us so we may rectify.
Notice: Registered trademark of products or corporate names are used only for explana-
tion and identification without intent of infringement.
Arcler Press publishes wide variety of books and eBooks. For more information about
Arcler Press and its products, visit our website at www.arclerpress.com
DECLARATION
Some content or chapters in this book are open access copyright free
published research work, which is published under Creative Commons
License and are indicated with the citation. We are thankful to the
publishers and authors of the content and chapters as without them this
book wouldn’t have been possible.
ABOUT THE EDITOR
Dr. Zoran Gacovski has earned his PhD degree at Faculty of Electrical
engineering, Skopje. His research interests include Intelligent systems and
Software engineering, fuzzy systems, graphical models (Petri, Neural and
Bayesian networks), and IT security. He has published over 50 journal and
conference papers, and he has been reviewer of renowned Journals. Currently,
he is a professor in Computer Engineering at European University, Skopje,
Macedonia.
TABLE OF CONTENTS
List of Contributors........................................................................................xv
List of Abbreviations..................................................................................... xxi
Preface................................................................................................... ....xxiii
x
Dynamic Equations................................................................................. 104
Control Design....................................................................................... 106
Simulation Results.................................................................................. 111
Conclusion............................................................................................. 115
References.............................................................................................. 116
xi
Carrying Manipulator.................................................................... 171
The Kinematics Simulation Algorithm of Manipulator............................. 176
The Kinematics Simulation of Manipulator.............................................. 177
Conclusions............................................................................................ 181
References.............................................................................................. 182
xii
Signals For Vehicle Detection And Tracking............................................. 234
Algorithms For Vehicle Detection And Tracking....................................... 235
Analysis Of Algorithms........................................................................... 244
Summary Of Algorithms......................................................................... 248
Conclusions And Future Directions......................................................... 248
Acknowledgment.................................................................................... 249
References.............................................................................................. 250
xiii
References.............................................................................................. 342
Chapter 16 Theory and Simulation Analysis of the Mode Shape and Normal
Shape Actuators and Sensors.................................................................. 345
Abstract.................................................................................................. 345
Introduction............................................................................................ 346
Actuator Equation................................................................................... 347
Sensor Equation...................................................................................... 356
Case Study: Theory And Simulation Analysis........................................... 358
Results And Discussion........................................................................... 360
Conclusion............................................................................................. 362
Acknowledgements................................................................................ 367
References.............................................................................................. 368
Index...................................................................................................... 403
xiv
LIST OF CONTRIBUTORS
A. Gasparetto
Polytechnic Department of Engineering and Architecture, University of Udine, Udine,
Italy
L. Scalera
Polytechnic Department of Engineering and Architecture, University of Udine, Udine,
Italy
Theodoros Foradis
Electrical and Computer Engineering, University of Patras, Patras, Greece
Kleanthis Thramboulidis
Electrical and Computer Engineering, University of Patras, Patras, Greece
Kleanthis Thramboulidis
Visiting Professor Helsinki University of Technology, Electrical & Computer
Engineering, University of Patras, Patras, Greece.
Syazili Roslan
Faculty of Bioresources & Food Industry, Universiti Sultan Zainal Abidin, Tembila
Campus, Besut, Terengganu Darul Iman, Malaysia
Hayan Basit
Faculty of Bioresources & Food Industry, Universiti Sultan Zainal Abidin, Tembila
Campus, Besut, Terengganu Darul Iman, Malaysia
Ashraf Elfasakhany
Department of Mechanical Engineering, Faculty of Engineering, Taif University, Al-
xv
Haweiah, Saudi Arabia
Tecnológico de Monterrey, Campus Ciudad Juárez, Ciudad Juarez, Mexico
Eduardo Yanez
Tecnológico de Monterrey, Campus Ciudad Juárez, Ciudad Juarez, Mexico
Karen Baylon
Tecnológico de Monterrey, Campus Ciudad Juárez, Ciudad Juarez, Mexico
Ricardo Salgado
Tecnológico de Monterrey, Campus Ciudad Juárez, Ciudad Juarez, Mexico
Pedram Bagheri
Department of Electrical and Computer Engineering, University of Calgary, Calgary,
Alberta, Canada
Laleh Behjat
Department of Electrical and Computer Engineering, University of Calgary, Calgary,
Alberta, Canada
Qiao Sun
Department of Mechanical and Manufacturing Engineering, University of Calgary,
Calgary, Alberta, Canada
Firas A. Raheem
Control and Systems Department, University of Technology, Baghdad, Iraq
Umniah I. Hameed
Control and Systems Department, University of Technology, Baghdad, Iraq
Valder Steffen, Jr
School of Mechanical Engineering, Federal University of Uberlândia, Uberlândia,
Brazil
xvi
Chen Liu
Department of Transportation and Environmental Engineering, Hiroshima University,
Higashi-Hiroshima, Japan
Yukio Fujimoto
Department of Transportation and Environmental Engineering, Hiroshima University,
Higashi-Hiroshima, Japan
Yoshikazu Tanaka
Department of Transportation and Environmental Engineering, Hiroshima University,
Higashi-Hiroshima, Japan
Shuncai Yao
Department of Electric Engineering, North University of China, Taiyuan, China
Jindong Tan
Department of Electrical Computer Engineering, Michigan Technological University,
Houghton, USA
Hongxia Pan
Department of Electric Engineering, North University of China, Taiyuan, China
G. Padmavathi
Department of Computer Science, Avinashilingam University for Women, Coimbatore,
Tamil Nadu, India
D. Shanmugapriya
Department of Information Technology, Avinashilingam University for Women,
Coimbatore, Tamil Nadu, India
M. Kalaivani
Department of Computer Science, Avinashilingam University for Women, Coimbatore,
Tamil Nadu, India
Mohammed A. Al-Moniee
Saudi Arabian Oil Company (Saudi Aramco), Dhahran, Saudi Arabia
Xiangyang Zhu
Saudi Arabian Oil Company (Saudi Aramco), Dhahran, Saudi Arabia
Rikke Markfoged
Danish Technological Institute (DTI), Aarhus, Denmark
xvii
Aabdullah H. Al-Wadei
Saudi Arabian Oil Company (Saudi Aramco), Dhahran, Saudi Arabia
Poul L. Pedersen
Danish Technological Institute (DTI), Aarhus, Denmark
Anders K. Tuxen
Danish Technological Institute (DTI), Aarhus, Denmark
Fuad I. Al-Nuwaiser
Saudi Arabian Oil Company (Saudi Aramco), Dhahran, Saudi Arabia
Lone Tang
Danish Technological Institute (DTI), Aarhus, Denmark
Thomas Lundgaard
Danish Technological Institute (DTI), Aarhus, Denmark
Jingyang Peng
Department of Mechanical Engineering, University of Saskatchewan, Saskatoon,
Canada
Xiongbiao Chen
Department of Mechanical Engineering, University of Saskatchewan, Saskatoon,
Canada
Mosfequr Rahman
Georgia Southern University, Statesboro, GA, USA
Mahbub Ahmed
Southern Arkansas University, Magnolia, AR, USA
Masud Nawaz
Georgia Southern University, Statesboro, GA, USA
Gustavo Molina
Georgia Southern University, Statesboro, GA, USA
Abdur Rahman
Edinboro University, Edinboro, PA, USA
xviii
Jwo Ming Jou
Department of Mechanical Engineering, Cheng Shiu University, Kaohsiung City,
Taiwan
Minoru Yaga
Department of Mechanical Systems, Faculty of Engineering, University of the Ryukyus,
Okinawa, Japan
Yusuke Uechi
Mie Metal Industry Co. Ltd., Mie, Japan
Hiroaki Ozono
Graduate School of Engineering and Science, University of the Ryukyus, Okinawa,
Japan
Masaaki Ishikawa
Department of Mechanical Systems, Faculty of Engineering, University of the Ryukyus,
Okinawa, Japan
Isao Teruya
Department of Mechanical Systems, Faculty of Engineering, University of the Ryukyus,
Okinawa, Japan
Y. Cao, X. B. Chen
Department of Mechanical Engineering, University of Saskatchewan, Saskatoon,
Canada
xix
LIST OF ABBREVIATIONS
xxii
PREFACE
1
A Brief History of Industrial Robotics
in the 20th Century
A. Gasparetto, L. Scalera
Polytechnic Department of Engineering and Architecture, University of
Udine, Udine, Italy
ABSTRACT
Industrial robotics is a branch of robotics that gained paramount importance
in the last century. The presence of robots totally revolutionized the
industrial environment in just a few decades. In this paper, a brief history
of industrial robotics in the 20th century will be presented, and a proposal
for classifying the evolution of industrial robots into four generations is
set forward. The characteristics of the robots belonging to each generation
are mentioned, and the evolution of their features is described. The most
INTRODUCTION
The idea to design and build some sort of beings, or devices that could
carry out repetitive or heavy tasks, thus relieving men from this burden,
dates back to ancient times. Since the Greek-Hellenistic age some of these
devices, which were named automata, have been designed and created by
ingenious inventors, belonging to several different civilizations through the
centuries. The term “automata” mainly refers to human-like devices, while
the term “robot” has a more general meaning.
The origin of the term “robot” is placed in more recent times: namely,
it comes from the Czech word “robota”, meaning “heavy work” or “forced
labor”. The introduction of this term is due to the Czech writer Karel Čapek
(1890-1938), who used it for the first time in 1920 in his novel “R.U.R.:
Rossum’s Universal Robots” (Čapek, 2004) .
On the other hand, the word “Robotics” was employed for the first time
by Isaac Asimov (1920-1992) in his novel “Runaround” (1942), contained in
the famous series “I, Robot”. In that novel he defined three rules concerning
the behavior of robots and the interaction with humans: these rules would
later be named the three Laws of Robotics (Asimov, 1942).
In the literature, there are not so many works dealing with Robotics
from a historical perspective. For instance, Ceccarelli dealt with this
topic in (Ceccarelli, 2001) and (Ceccarelli, 2004), while Gasparetto
(2016) presents a historical outline of robotics from ancient times until the
Industrial Revolution.
Industrial applications of Robotics gained a paramount importance in the
last century. The beginning of “Industrial Robotics”, as we currently define
it, can be dated back to the 1950’s, although some kinds of automatization in
the industrial environment started to appear since the times of the Industrial
Revolution.
In this paper, the main milestones of the history of industrial robotics,
from its beginning (in the 1950’s and even earlier) to the end of the 20th
A Brief History of Industrial Robotics in the 20th Century 3
Figure 1. Joseph Engelberger and George Devol (left); the Unimate robot
(right).
The PUMA (an acronym for Programmable Universal Machine for Assembly)
was considered for many decades the archetype of the anthropomorphic
robots, and its kinematics is still taken as an example in several robotic
books at the undergraduate as well at the graduate level.
In the meanwhile, other companies developed and manufactured other
types of industrial robots. For instance, KUKA developed in 1973 the
Famulus robot, whose name in Latin means “servant”. In 1974 Cincinnati
Milacron, a major machine tool manufacturer, developed a robot called T3
(an acronym for “The Tomorrow Tool”), which was installed in several
automotive plants, and especially in the Volvo plants in Sweden. T3 was
the first commercially available minicomputer-controlled industrial robot
(Figure 5).
In 1974, the Swedish company ASEA (now ABB) started the production
of the robots of the famous and successful IRB series, well known worldwide
also for their typical orange color. The first robot of this series, that was
issued for more than 20 years, was the IRB-6, which was largely employed
in productive sites for complex tasks (machining, arc-welding), for its ability
to move smoothly along continuous paths (Figure 6).
In the same year, the Japanese company Hitachi developed the robot HI-
T-HAND Expert, that is relevant in the history of industrial robotics for the
precision it reached in insertion operations (it could insert mechanical parts
with a clearance of about 10 micrometers). It was also provided with a force
feed-back control system and a flexible wrist mechanism.
on-line (the operator could use a teach box with a keyboard) or off-line, being
connected to a PLC or a PC, which allowed to use a high-level language for
motion programming and enabled the robots to be interfaced with a CAD
or a database. The possibility of high-level, off-line programming enlarged
the operational potential of the robots: for in-stance, they could elaborate
data from sensor reading, in order to adjust the robot movements taking into
account changes in the environment (e.g. changes in position and orientation
of the work pieces). Moreover, the diagnostic capabilities could be greatly
enhanced: these robots could produce not only an indication of failure
detection, but also a report on the location and on the type of the failure.
In addition, some sort of “intelligence” was present in the robots of the
third generation, with some (although limited) adaptive capabilities. These
capabilities could be employed in some more complex tasks (such as tactile
inspection, assembly operations, arc welding), by using the data coming
from vision or perception systems to locate the objects and the work pieces
and guide the joint movements according to the task to be performed, taking
into account the possibility of small changes in the position of the objects.
Between the end of the 1970’s and the beginning of the 1980’s, other
scientific and technical improvements contributed to the diffusion of robots.
In 1978, a novel kinematic structure was proposed by the Japanese
scientist Hiroshi Makino from Yamanashi University. Such a structure was
made of three revolute joints with parallel axes and a prismatic joint lying at
the end of the kinematic chain (Figure 7). The robot with this structure was
named SCARA (an acronym from “Selective Compliance Assembly Robot
Arm”), since its compliance in the horizontal direction resulted lower than
the compliance in the vertical direction. For this reason, as well as for the
lightness of the kinematic chain (that allowed a simpler and faster controller),
this robot was suitable to be employed in tasks such as the assembly of small
objects (Makino & Furuya, 1980).
Another relevant technical improvement in industrial robotics was the
appearance of direct drive actuated robot. The first prototype of this kind was
the CMU Direct Drive Arm (Asada & Kanade, 1983), developed in 1981 by
Kanade and Asaka at Carnegie Mellon University (Pittsburgh, USA). This
kind of robot featured higher accuracy and faster operations because the
motors connected directly to the arms eliminating the need for intermediate
gear or chain systems.
Both the aforementioned findings were employed in the AdeptOne (Figure
8), the first commercially available direct-driven SCARA robot (1984).
A Brief History of Industrial Robotics in the 20th Century 11
lines, especially in Japan. For this reason, Japan became in the 1980’s not
only the world leader in robot manufacturing (with more than 40 companies
that built and sold robots worldwide), but also in the production of electronic
consumer goods.
Robotics in the 1980’s was a rising star, not only in Japan but in all the
developed countries. It appeared as a promising field that drew the interest
of journalist, scientists, policy makers and also common people.
From the scientific and technical point of view, this is the time when the
robots became even more versatile, by exploiting important improvements
both with respect to the hardware and the software. Concerning the former
aspect, robots started to be provided with advanced sensors (e.g. cameras,
force sensors, laser scanners); concerning the latter aspects, the control
software became more “intelligent” by introducing come techniques related
to Artificial Intelligence. Both these aspects increased the versatility and
the flexibility of the robots, which could be employed in more and more
complex tasks.
Despite the significant progress undergone in the 1980’s, the need for
robots that could carry out task at high speed pushed the scientific research
to design innovative kinematic structures. The idea of employing parallel
kinematic chains instead of the classical serial kinematic chains was put
forward and led to a type of lightweight robot featuring the capability of
moving at high speed. The archetype of this kind of robots was the Delta
robot (that appeared in 1992), conceived by the Swiss scientist Reymond
Clavel at the Ecole Poly-technique Fédérale de Lausanne (EPFL). This
type of robot, designed by Clavel in his PhD thesis, had three translational
DOFs and one rotational DOF (Clavel, 1991). With respect to serial robots,
parallel robots featured a smaller work-space, but the capability of operating
at much higher speed. The kinematic architecture of the Delta robot was
copied in many parallel manipulators, devoted to high speed pick-and-place
operations.
The first application of Delta robots was developed by the Swiss
company Demaurex in 1992: six Delta robots were operating inside a work
cell for loading pretzels into trays (Figure 9). Some years afterwards (1998)
ABB developed the Flex-Picker, the world’s fastest picking robot, based on
the structure of the Delta robot (Figure 10).
In this decade, the first attempts to synchronize robots took place:
in 1994 Motoman developed MRC, the first robot control system which
provided the synchronized control of two robots (up to 21 axes).
A Brief History of Industrial Robotics in the 20th Century 13
The end of the third generation is conventionally set to the end of the century;
beginning from the year 2000, the industrial robots are considered to belong
to the fourth generation (which extends up to the current days). Such robots
feature high-level “intelligent” capabilities (such as performing advanced
computations, logical reasoning, deep learning, complex strategies,
collaborative behavior).
Figure 9. Delta robots in the Demaurex work cells for the packaging of pretzels.
CONCLUSION
In this paper, a brief history of industrial robotics in the 20th century
was presented. The evolution of the industrial robots was conventionally
categorized into four generations, of which the first three cover the timespan
from the 1950’s to the end of the century.
14 Mechatronics and Robotics
In this historical sketch, not only the scientific and technical evolution
was taken into account, but some considerations about the economic and
geopolitical issues that determined the diffusion of industrial robots, were
also done.
The evolution of industrial robotics is not over, but is still developing in
the current days: innovative ideas and novel hardware devices, together with
some new programming techniques connected with Artificial Intelligence,
are revolutionizing the concept of industrial automation and giving a new
youth to the factory environment.
CONFLICTS OF INTEREST
The authors declare no conflicts of interest regarding the publication of this
paper.
A Brief History of Industrial Robotics in the 20th Century 15
REFERENCES
1. Asada, H., & Kanade, T. (1983). Design of Direct-Drive Mechanical
Arms. Journal of Vibration, Acoustics, Stress, and Reliability, 105,
312-316. https://doi.org/10.1115/1.3269106
2. Asimov, I. (1942). I, Robot. In Runaround. New York: Spectra Books.
3. Birnie, J. V. (1974). Practical Implications of Programmable
Manipulators. Industrial Robot: An International Journal, 1, 122-128.
https://doi.org/10.1108/eb004720
4. Capek, K. (2004). R.U.R. Rossum’s Universal Robots. Fayetteville:
Penguin Ed.
5. Ceccarelli, M. (2001). A Historical Perspective of Robotics toward the
Future. Fuji International Journal of Robotics and Mechatronics, 13,
299-313. https://doi.org/10.20965/jrm.2001.p0299
6. Ceccarelli, M. (2004). Fundamentals of Mechanics of Robotic
Manipulation. Dordrecht: Kluwer/Springer. https://doi.
org/10.1007/978-1-4020-2110-7
7. Clavel, R., (1991). Conception d’un robot parallèle rapide à 4 degrés
de liberté. Ph.D. Thesis, Lausanne, Switzerland: EPFL.
8. Devol, G. (1954). Programmable Article Transfer. https://patents.
google.com/patent/US2988237
9. Gasparetto, A. (2016). Robots in History: Legends and Prototypes
from Ancient Times to the Industrial Revolution. In Explorations in
the History of Machines and Mechanisms. History of Mechanism and
Machine Science (vol. 32, pp. 39-49). Berlin/Heidelberg: Springer.
https://doi.org/10.1007/978-3-319-31184-5_5
10. Karlsson, J. M. (1991). A Decade of Robotics; Analysis of the
Diffusion of Industrial Robots in the 1980s by Countries, Application
Areas, Industrial Branches and Types of Robots. Stockholm, Sweden:
Mekanforbundets Forlag.
11. Koetsier, T. (2019). Towards the Global Intelligent Machine. In The
Ascent of GIM, the Global Intelligent Machine. History of Mechanism
and Machine Science (vol 36, pp. 293-317). Berlin/Heidelberg:
Springer. https://doi.org/10.1007/978-3-319-96547-5
12. Lagerberg, T., & Jonson, J. (2016). Une bio du robot. Petite histoire de
la robotique industrielle. ABB Review 3/16. http://new.abb.com/about/
technology/abb-review
16 Mechatronics and Robotics
2
From Mechatronic Components to
Industrial Automation Things: An
IoT Model for Cyber-Physical
Manufacturing Systems
ABSTRACT
IoT is considered as one of the key enabling technologies for the fourth
industrial revolution that is known as Industry 4.0. In this paper, we consider
the mechatronic component as the lowest level in the system composition
hierarchy that tightly integrates mechanics with the electronics and software
required to convert the mechanics to intelligent (smart) object offering well
defined services to its environment. For this mechatronic component to be
INTRODUCTION
Based on one of the most commonly used definitions, the term Mechatronics
emphasizes on the synergistic integration of the three discipline areas, i.e.,
mechanical engineering, electronics and intelligent computer control, in the
design and manufacture of products and processes [1] , i.e., it emphasizes
on synergy. What is not clear by this definition is the level at which this
integration should be performed, i.e., at the system level, which is the
traditional approach, or at the subsystem or even at the mechanical unit
(component) level.
The latter is proposed in Model Integrated Mechatronics [2] and
refined with the 3+1SysML-view model [3] [4]. This approach defines
the Mechatronic component as the main building block that abstracts the
mechanical object to the software level, and transforms it to a smart object
by adding additional functionality to the one offered by the mechanical part.
The so constructed mechatronic components are integrated with cyber
components and humans to construct the industrial automation system. This
approach slightly finds its road to production in the context of Industry 4.0,
e.g., [5], since it greatly reduces the coupling between the system components
compared to the traditional one, which considers the integration of the three
disciplines at the system integration level.
From Mechatronic Components to Industrial Automation ... 19
Contiki based Industrial Automation Thing is presented along with the case
study. The model-to-model transformer for the C language and the Contiki
operating system is presented in section 5 and the paper is concluded in the
last section.
RELATED WORK
IoT offers new levels of connectivity in the industrial domain that may
lead to higher efficiency, flexibility, and interoperability among industries
[20]. However, not only many definitions exist for the IoT but also several
models. These models, e.g., ETSI, IETF, SENSEI, have been developed to
capture the key concepts of the domain and provide the infrastructure to
develop frameworks and architectures for the systems based on IoT. Fei
et al. [21] claim that even though IoT has been used in various application
domains there is still no clear and uniform definition and architecture
about it. Therefore, the IoT domain suffers from an inconsistent usage and
understanding of the meaning of several key terms, as also claimed in [22].
The definition of an IoT domain model is a prerequisite for defining IoT
Architectures. A detailed discussion on the IoT models can be found in
[23], where the IoT-A reference model for IoT is presented and validated.
Authors in [24] present and discuss applications of IoT in Manufacturing
and describe a five-layer architecture for manufacturing based on IoT.
Authors in [25] focus on the device nature of Thing and consider sensors,
actuators and controllers as IoT devices, i.e., things. They focus on the data
field structures and evaluate the benefit of using an IP smart gateway as the
decentralized peripheral to integrated sensors, actuators and the controller
and claim that this may improve the performance of IoT devices. We do not
agree with the use of sensors and actuators as first-class model elements
in the high-level design specification of the system. Sensors and actuators
are just technology artefacts used to integrate the physical with the cyber
world so they have no place in the high-level design spec of the system.
In our approach, the IAT, which encapsulates sensors, actuators and the
controller, plays the role of the Thing and represents the key construct in the
IoT manufacturing environment.
MDD becomes more and more popular in the development of embedded
software systems and various reports refer efficiency gains, from up to 50%,
for example, in the development in the car industry [26], with high error
reductions and a rapid increase of the maturity level of developed products.
MDD is considered as a promising solution to address the complexity of
22 Mechatronics and Robotics
[30]. Authors do not define any DSML but they construct class diagrams
and state charts using only primitive UML model elements. Authors in [31]
describe an approach to define a visual DSML for the IoT based on UML.
They model the Thing, which they consider as key construct for building
an IoT system, using the UML component construct and its interface using
provided and required interfaces. Authors do not address the mapping of the
conventional object- oriented (OO) interfaces of the Thing with the ones of
the REST paradigm.
Yingfeng Zhang et al. [32] present a real-time information capturing
and integration architecture of the internet of manufacturing things (IoMT)
to provide a new paradigm by extending the techniques of IoT to the
manufacturing field.
They use the term manufacturing thing but they do not focus on its
structure and its development process. Instead, they focus on the overall
architecture for a manufacturing system and describe a framework with
focus on real-time tracking and tracing for the dynamic monitoring of the
manufacturing process.
To the best of our knowledge there is no other work that focuses on the
automation of the transformation process of the conventional mechatronic
component to an IoT compliant one, i.e., to an Industrial Automation Thing
ready to be integrated into the IoT-based industrial environment.
Figure 1. Example key concepts and interaction in the IoT-A model (De et al.
2011).
Thus, the heating Completed and mixing Completed signals sent by the
Smart Silo will trigger an asynchronous without a reply reaction to the
Smart Silo client, e.g., the type A liqueur generation process, through the
corresponding Receptions captured in the Smart Silo User If.
or call back functions, imply a tight coupling among the smartSilo and the
components that use its behaviour, in the sense that these interfaces should
be known in advance for the development of the component’s clients.
A model-to-model transformer is required to automate the process of
generating the IAT. This transformer will use as input the properly annotated
with the DSML conventional mechatronic component. One approach is to
mark the UML design specification of the mechatronic component with the
stereotypes defined by the UML4IoT profile. If a UML design is not available
then the source code of the cyber part of the mechatronic component is
properly annotated with specific annotations that have been defined based on
the UML profile. An example of annotated code with Java-like annotations
is given in Listing 2, where objects and object types of the cyber-physical
component as well as their properties that should be exported to the IoT
environment are properly annotated.
Figure 7. The silo industrial automation thing based on a hardware silo simula-
tor.
and used on demand based on requirements assuming that the physical part
supports the requested new behavior. Resources, Resource Instances, Objects,
Object Instances which are exposed by the IAT as well as their attributes,
are accessed by the clients of the IAT through the device management and
service enablement interface (DM&SE If). IPSO smart objects [38] have
been adopted in this work to satisfy the requirement for interoperability.
The Guile module that implements the c-parser scheme procedures was
developed as a high-level recursive descent parser, based on the guile-parser-
combinators module (https://git.dthompson.us/guile-parser-combinators.
git). GNU Guile is an implementation of the Scheme programming language
(https://www.gnu.org/software/guile/). Autogen collaborates with the Guile
VM through guile procedure calls embedded in the template file. The output
of the transformation process is the source code of the IAT in a format ready
to be compiled with the Contiki operating system and then deployed.
Transformation Rules
The following rules that apply for each object of the lwm2m client have been defined.
Rule 1:Create wrapper functions for annotated behaviors.
For each function with the BehaviorResource annotation create a
wrapper function with input parameters:
lwm2m_context_t *ctx, const uint8_t *arg, size_targsize, uint8_t
*outbuf, size_t outsize
E.g., Source:static int fill(void);
Target: static intfill(lwm2m_context_t *ctx, const uint8_t *arg, size_
targsize, uint8_t *outbuf, size_t outsize);
Rule 2:Create getters/setters functions for each property annotated with
the PrimitiveRes annotation.
For each attribute annotated with the PrimitiveRes annotation create the
corresponding read and write function depending on the applied operations
on the attribute defined in the annotation.
E.g., for the silo_state property
staticintget_silo_state(lwm2m_context_t *ctx, uint8_t *outbuf, size_t
outsize) {
char *value;
value = get_silo_state_inString(silo->state);
returnctx->writer->write_string(ctx, outbuf, outsize, value,
strlen(value));}
Rule 3: Construct the Resource model.
For each annotated attribute or function create an entry using the
LWM2M_RESOURCE_CALLBACK macro of the lwm2m implementation
which is integrated into the Contiki OS. E.g., for the silo_state
36 Mechatronics and Robotics
CONCLUSION
In this paper, we consider the tight integration of the physical world with the
cyber one at the mechatronic component level. A mechatronic component
offers its functionality through well-defined mechanical, electrical and
software interfaces. In this sense the industrial automation system is a
composition of mechatronic components along with cyber components and
humans. IoT is adopted for the integration of these components to exploit
the benefit of this technology and UML4IoT is utilized to automatically
transform the conventional mechatronic component into an IoT compliant
cyber-physical one, i.e., to an Industrial Automation Thing. The IoT model
used in the UML4IoT approach is extended towards a complete IoT model
for the manufacturing domain. The transformation rules required for the
development of the model-to-model transformer have been developed and
validated through a prototype implementation of the liqueur production
laboratory system. The prototype implementation of the silo industrial
automation thing based on a Contiki enabled embedded board and the C
language is used to demonstrate the applicability and the effectiveness of
the proposed approach.
38 Mechatronics and Robotics
REFERENCES
1. UNESCO Chair (2016) Home Page. On Mechatronics and Mechatronics
Research and Application Center. http://mecha.ee.boun.edu.tr/
2. Thramboulidis, K. (2005) Model Integrated Mechatronics—Towards
a new Paradigm in the Development of Manufacturing Systems.
IEEE Transactions on Industrial Informatics, 1, 54-61. https://doi.
org/10.1109/TII.2005.844427
3. Thramboulidis, K. and Scholz, S. (2010) Integrating the 3+1 SysML
View Model with Safety Engineering. 2010 IEEE Conference on
Emerging Technologies and Factory Automation (ETFA), Bilbao, Spain,
13-16 September 2010, 1-8.https://doi.org/10.1109/etfa.2010.5641353
4. Thramboulidis, K. (2012) Overcoming Mechatronic Design
Challenges: The 3+1 SysML-View Model. The Computing Science
and Technology International Journal, 2, 6-14.
5. FESTO (2016) Industry 4.0: Efficient Engineering Processes
with “OPAK”. https://www.festo.com/net/en_corp/SupportPortal/
MobilePressDetails.aspx?documentId=368146&q=
6. OASIS (2016) Devices Profile for Web Services (DPWS) Specification.
http://docs.oasis-open.org/ws-dd/ns/dpws/2009/01
7. Angulo, P., Guzmán, C.C., Jiménez, G. and Romero, D. (2016) A
Service-Oriented Architecture and Its ICT-Infrastructure to Support
Eco-Efficiency Performance Monitoring in Manufacturing Enterprises.
International Journal of Computer Integrated Manufacturing, 30, 202-
214. https://doi.org/10.1080/0951192X.2016.1145810
8. Morgan, J. and O’Donnell, G.E. (2015) Enabling a Ubiquitous
and Cloud Manufacturing Foundation with Field-Level Service-
Oriented Architecture. International Journal of Computer Integrated
Manufacturing, 30, 442-458.
9. He, W. and Xu, L. (2015) A State-of-the-Art Survey of Cloud
Manufacturing. International Journal of Computer Integrated
Manufacturing, 28, 239-250.https://doi.org/10.1080/095119
2X.2013.874595
10. Ghimire, S., Luis-Ferreira, F., Nodehi, T. and Jardim-Goncalves,
R. (2016) IoT Based Situational Awareness Framework for Real-
Time Project Management. International Journal of Computer
Integrated Manufacturing, 30, 74-83.https://doi.org/10.1080/095119
2X.2015.1130242
From Mechatronic Components to Industrial Automation ... 39
11. Ren, L., Zhang, L., Wang, L., Tao, F. and Chai, X. (2014) Cloud
Manufacturing: Key Characteristics and Applications. International
Journal of Computer Integrated Manufacturing, 30, 501-515. https://
doi.org/10.1080/0951192X.2014.902105
12. Bi, Z., Xu, L.D. and Wang, C. (2014) Internet of Things for Enterprise
Systems of Modern Manufacturing. IEEE Transactions on Industrial
Informatics, 10, 1537-1546. https://doi.org/10.1109/TII.2014.2300338
13. Bradley, D., Russell, D., Ferguson, I., Isaacs, J., MacLeod, A. and
White, R. (2015) The Internet of Things—The Future or the End of
Mechatronics. Mechatronics, 27, 57-74.
14. Thramboulidis, K. and Christoulakis, F. (2016) UML4IoT—A UML
Profile to Exploit IoT in Cyber-Physical Manufacturing Systems.
Computers in Industry, 82, 259-272.
15. Dunkels, A., Gronvall, B. and Voigt, T. (2004) Contiki—A Lightweight
and Flexible Operating System for Tiny Networked Sensors. 29th
Annual IEEE International Conference on Local Computer Networks,
16-18 November 2004, 455-462. https://doi.org/10.1109/lcn.2004.38
16. Thramboulidis, K. (2015) A Cyber-Physical System-Based Approach
for Industrial Automation Systems. Computers in Industry, 72, 92-102.
17. Eclipse (2016) Papyrus for IoT—A Modeling Solution for IoT. https://
www.eclipse.org/community/eclipse_newsletter/2016/april/article3.
php
18. Open Mobile Alliance (OMA) (2015) Lightweight Machine to
Machine Technical Specification. OMA-TS-LightweightM2M-V1_0-
20151214-C, Candidate Version 1.0. 14 December 2015.
19. Kovatsch, M., Duquennoy, S. and Dunkels, A. (2011) A Low-Power
CoAP for Contiki. 8th IEEE International Conference on Mobile Ad-
Hoc and Sensor Systems, Valencia, Spain, 17-21 October 2011, 855-
860. https://doi.org/10.1109/MASS.2011.100
20. Khaleel, H., Conzon, D., Kasinathan, P., Brizzi, P., Pastrone, C., et al.
(2015) Heterogeneous Applications, Tools, and Methodologies in the
Car Manufacturing Industry through an IoT Approach. IEEE Systems
Journal, PP, 1. https://doi.org/10.1109/JSYST.2015.2469681
21. Tao, F., Cheng, Y., Xu, L.D., Zhang, L. and Li, B.H. (2014) CCIoT-
CMfg: Cloud Computing and Internet of Things-Based Cloud
Manufacturing Service System. IEEE Transactions on Industrial
Informatics, 10, 1435-1442.https://doi.org/10.1109/TII.2014.2306383
40 Mechatronics and Robotics
22. Bauer, M., Bui, N., De Loof, J., Magerkurth, C., Nettstrter, A., Stefa, J.
and Walewski, J. (2013) Iot Reference Model. In: Bassi, A., et al., Eds.,
Enabling Things to Talk, Springer, Berlin, Heidelberg, 113-162.https://
doi.org/10.1007/978-3-642-40403-0_7
23. IoT-A (2016) Final Architectural Reference Model for the IoT. http://
www.meet-iot.eu/deliverables-IOTA/D1_5.pdf
24. Tao, F., Zuo, Y., Xu, L.D. and Zhang, L. (2014) IoT-Based Intelligent
Perception and Access of Manufacturing Resource toward Cloud
Manufacturing. IEEE Transactions on Industrial Informatics, 10, 1547-
1557.https://doi.org/10.1109/TII.2014.2306397
25. Diaz-Cacho, M., Delgado, E., Falcon, P. and Barreiro, A. (2015) IoT
Integration on Industrial Environments. 2015 IEEE World Conference
on Factory Communication Systems (WFCS), Palma de Mallorca, 27-
29 May 2015, 1-7. https://doi.org/10.1109/wfcs.2015.7160553
26. Broy, M., Kirstan, S., Krcmar, H., Schatz, B. and Zimmermann, J.
(2013) What Is the Benefit of a Model-Based Design of Embedded
Software Systems in the Car Industry? In: Software Design and
Development: Concepts, Methodologies, Tools, and Applications, IGI
Global, 310-334.
27. Nguyen, X.T., Tran, T., Baraki, H. and Geihs, K. (2015) FRASAD:
A Framework for Model-Driven IoT Application Development.
2015 IEEE 2nd World Forum on Internet of Things (WF-IoT),
Milan, 14-16 December 2015, 387-392. https://doi.org/10.1109/WF-
IoT.2015.7389085
28. Malavolta, I. and Muccini, H. (2014) A Study on MDE Approaches
for Engineering Wireless Sensor Networks. 40th EUROMICRO
Conference on Software Engineering and Advanced Applications
(SEAA), Verona, 27-29 August 2014, 149-157.https://doi.org/10.1109/
seaa.2014.61
29. Conzon, D., Brizzi, P., Kasinathan, P., Pastrone, C., Pramudianto,
F. and Cultrona, P. (2015) Industrial Application Development
Exploiting IoT Vision and Model Driven Programming. 2015 18th
International Conference on Intelligence in Next Generation Networks
(ICIN), Paris, 17-19 February 2015, 168-175. https://doi.org/10.1109/
icin.2015.7073828
30. Prehofer, C. (2015) Models at REST or Modelling RESTful Interfaces
for the Internet of Things. IEEE 2nd World Forum on Internet of Things
(WF-IoT), Milan, 14-16 December 2015, 251-255.
From Mechatronic Components to Industrial Automation ... 41
31. Eterovic, T., Kaljic, E., Donko, D., Salihbegovic, A. and Ribic, S. (2015)
An Internet of Things Visual Domain Specific Modeling Language
Based on UML. 2015 XXV International Conference on Information,
Communication and Automation Technologies (ICAT), Sarajevo, 29-
31 October 2015, 1-5. https://doi.org/10.1109/icat.2015.7340537
32. Zhang, Y., Zhang, G., Wang, J., Sun, S., Si, S. and Yang, T. (2015)
Real-Time Information Capturing and Integration Framework of the
Internet of Manufacturing Things. International Journal of Computer
Integrated Manufacturing, 28, 811-822. https://doi.org/10.1080/09511
92X.2014.900874
33. De, S., Barnaghi, P., Bauer, M. and Meissner, S. (2011) Service
Modelling for the Internet of Things. 2011 Federated Conference on
Computer Science and Information Systems (FedCSIS), Szczecin, 18-
21 September 2011, 949-955.
34. Haller, S. (2010) The Things in the Internet of Things. Internet of
Things Conference 2010, Tokyo, 29 November-1 December 2010.
35. Ferry, P., et al. (2014) IoTLink: An Internet of Things Prototyping
Toolkit. IEEE International Conference on Ubiquitous Intelligence and
Computing, Bali, 9-12 December 2014, 1-9.
36. Basile, F., Chiacchio, P. and Gerbasio, D. (2013) On the Implementation
of Industrial Automation Systems Based on PLC. IEEE Transactions
on Automation Science and Engineering, 10, 990-1003. https://doi.
org/10.1109/TASE.2012.2226578
37. OMG (2015) Object Management Group, OMG Unified Modeling
Language (OMG UML), Version 2.5, OMG Document Number
Formal/2015-03-01. http://www.omg.org/spec/UML/2.5
38. Internet Protocol for Smart Objects (IPSO) Alliance (2014) IPSO
Smart Object Guideline. IPSO Smart Object Committee, 21 September.
https://www.ipso-alliance.org/ipso-community/resources/smart-
objects-interoperability/
CHAPTER
3
The 3+1 SysML View-Model
in Model Integrated
Mechatronics
Kleanthis Thramboulidis
Visiting Professor Helsinki University of Technology, Electrical & Computer
Engineering, University of Patras, Patras, Greece.
ABSTRACT
Software is becoming the driving force in today’s mechatronic systems. It
does not only realize a significant part of their functionality but it is also
used to realize their most competitive advantages. However, the traditional
development process is wholly inappropriate for the development of these
systems that impose a tighter coupling of software with electronics and
INTRODUCTION
Software does not only implement a significant part of the functionality of
today’s mechatronic systems, but it is also used to realize their most competitive
advantages. It is the evolving driver for innovations in many mechatronic
systems and in general it is considered as the driving force in improving this
kind of systems. However, the traditional development process is wholly
inappropriate for the development of systems characterized by complexity,
dynamics and uncertainty as is the case with today’s mechatronic systems
[1]. According to the traditional development process the constituent parts
of the mechatronic system, i.e. the mechanical, electronic and software, that
constitute the system are developed independently and then are integrated to
compose the final system.
The software development starts when the development of electronic
and mechanical is already at a stage where any change in these parts is
expensive and time consuming. This is why the mechanical and electronic
properties impose several constrains and narrow the solution space for
software development. Moreover, as claimed in [2] “the actual cooperation
during the construction is less developed. There is no joint development
process, no joint tool usage, no joint modeling formalism and no joint
analysis. Every discipline has its own approaches”. As a result of this, the
current process that is traditionally divided into software, electronics and
mechanics, emphasizes on domain-spanning design methods and tools and
is unable to address the demand for synergetic mechatronic dependability
predictions; this is why many products suffer from severe dependability
problems [3].
The 3+1 SysML View-Model in Model Integrated Mechatronics 45
Related Work
Several researchers are already working in the direction of improving the
effectiveness of the development process of mechatronic systems. Schafer
and Wehrheim [2] survey on current developments in mechatronics and
present the architecture of their mechatronic rail system that seems to provide
an excellent platform for studding and analyzing future developments and
research chalenges in mechatronic systems. They identify the need of an
integrated framework for the construction of mechatronic systems and
they discuss future trends in mechatronics especially from the software
engineering point of view. Habib [10] argues on the urgent need for theories,
models, and tools that should facilitate modeling, analysis, synthesis,
simulation, and prototyping of mechatronic systems. He emphasizes the
argument that the approach based on optimization within each domain
separately will not result in the optimum system design, and he proposes
a data and model integration approach to address the integration problem.
Burmester et al. [11] claim that in today’s mechatronic systems most of
the control and reconfiguration functionality is realized in software. They
present mechatronic UML to exploit the Model Driven Architecture
approach for the design of hybrid mechatronic real-time systems that have
to fulfill safety-critical requirements. “Mechatronic UML” is defined as an
extension of UML to built platform independent models for mechatronic
systems. Various UML models have been extended to cover the requirements
of modeling the structural view as well as the behavioral view of the system.
However, the proposed extension is used to model only the software part of
the mechatronic system. Authors in [12] briefly refer to a process model of
Robert Bosch GmbH for the development of mechatronic systems in Motor
Vehicles to support aspects such as reuse, exchangeability, scalability and
distributed development. They use the concept of mechatronic component,
even though not well defined, as the basic construct of their process. They
argue: a) on the need of a clear specification of component interfaces;
and b) the great contribution of re-use to increase the quality properties
of mechatronic systems and decrease development time. Nordmann [13]
is using the concept of mechatronic component and presents an example
of using Active Magnetic Bearings to increase performance, reliability,
reusability and longer lifetime. Authors in [14] propose for the development
of multidisciplinary systems, such as mechatronics, the integration of the
various domain-specific tools. They mainly focus on the integration of used
data and models and not on a process level integration. Moreover, none
of these approaches provide a high level architecture for an integrated,
50 Mechatronics and Robotics
its perspective so it has to provide specific interfaces to the tools of the other
views, in order to implement the interactions of its own part to the other
parts of the primitive MTC. The arrows that cross the boundaries between
the three views in Figure 3 represent the interactions of the corresponding
models and have to be implemented by specific interfaces of the model
refinement, analysis and execution tools of the three disciplines. The AP233
or more formally the ISO 10303-233 standard for systems engineering [17],
that provides a data exchange format for the reliable interchange of data
between software tools may be exploited to effectively implement these
interactions. The execution and analysis of the primitive MTC is obtained
through a collaboration of the corresponding tools of the three views. It is
clear that the contribution of the three views’ specific tools is restricted at
the primitive MTC internal level while the execution and analysis of MTC
models is done at the MTS level with the coordination and synchronization
between MTCs carried by this level. This makes the tool integration a major
challenge in the domain of mechatronic systems.
The 3+1 SysML view-model when used with the MTS V-Model that is
described in one of the following sections, promote the synergistic integration
of the three constituent parts of the mechatronic system and emphasizes the
importance of a common model for the system. However, this model can
also be used with the traditional development process that is based on the
independent development of constituent parts of the mechatronic system and
their subsequent integration. Even in this case the existence of a common
model for the system greatly improves the effectiveness of the development
process.
It should be noted that in each view the corresponding discipline’s
specific architectures and tools may be exploited, as for example the 4+1
architectural view [18] that may be exploited in the context of the e-view by
the software engineer or the MTC developer/integrator.
the MTC. In this section the modeling of the MTS and MTC levels using
SysML is considered. It is evident that specific interfaces have to be defined
for the integration of the different views and these interfaces have to be
realized by the tools used in the various disciplines to create a completely
integrated tool chain to support the MTS development process. The SysML
to AP233 mapping [19] is towards this direction.
combination of the above. It is assumed that the MTC developer has already
performed a QoS analysis for the MTC. All this information comprises the
offered QoS characteristics of the MTC [4]. It should also be stated that the
MTC developer does not know during the MTC’s development time all the
systems that this MTC will be used in the future.
According to the top-down approach, for every system-level
responsibility that has not been assigned to a single MTC, a set of abstract
MTCs is defined and the required MTC responsibilities are specified
along with the required collaboration scheme. Required system-level QoS
characteristics are decomposed to derive the component-level required QoS
characteristics.
This process results to the definition of the required QoS characteristics
at the level of constituent MTCs. The process of deriving MTC-level QoS
characteristics from system-level ones is a complex process and has to be
defined. At this time the engineer has well defined required specifications
(functional and nonfunctional, including QoS characteristics) for every
abstract MTC. Using these required QoS characteristics the engineer is able
to select from the market or his components repository the ones that their
offered QoS characteristics meet the required ones [8]. If such MTC’s do not
exist they have to be further analyzed in order to be developed.
Advantages and disadvantages for the above approaches that result in
the definition of the architecture of the MTS of its composite MTCs may
be identified but it is expected that in the real MTS development process
a combination of both approaches will be used resulting to a hybrid more
efficient approach.
The above process, bottom-up or top-down, is again applied to every
composite MTC that has to be developed. It is applied iteratively down to
the primitive-MTC level; the identification of primitive MTCs signals the
end of this iteration. For each composite MTC the system modeling process
as defined by the left-hand side part of the V-Model is followed. Analysis
is applied and its architecture is defined in terms of constituent components
(composite and/or primitive). Sequence diagrams are defined to realize use
cases of the MTC and identify the activities that are involved in the specific
use case.
This is not the case for primitive MTCs that have to bypass the system
process and follow a synergistic integration of the three constituent parts,
i.e. mechanic, electronic and software (MTC synthesis) as shown in the
bottom of the V-model in Figure 7.
The 3+1 SysML View-Model in Model Integrated Mechatronics 59
CONCLUDING REMARKS
The traditional approach in the development of mechatronic systems is unable
to address the needs of today’s complex mechatronic systems. An integrated
framework for the construction of mechatronic systems is missing. The work
presented in this paper attempts to contribute to this direction by: a) using
SysML to define the artifacts of the MIM paradigm; b) proposing the 3+1
SysML view-model imposed by the MIM architecture; and c) extending the
well accepted and widely used in the software domain V-model to address
the demands of the mechatronic system development process. However,
the challenges for a fully automated MTS development process crosses the
boundaries of the three disciplines of mechatronic systems and impose a
joint effort and collaboration between computer science, electronics and
mechanics. The current status of discipline isolation imposed in many cases
by the existing structure of engineering degree programs makes the task
even more complicated.
ACKNOWLEDGMENTS
Part of this work has been funded by TIKOSU, a project belonging to
the Digital Product Process -program of the Finnish Funding Agency for
Technology and Innovation (TEKES). The author wishes to thank the
partners and especially Jarmo Alanen, Kari Koskinen, and Seppo Sierla for
fruitful comments on these ideas. Thanks are also due to Piero Larizza, Eric
Coatanea and Jussi Suomela for discussions on these concepts.
62 Mechatronics and Robotics
REFERENCES
1. G. Rzevski, “On conceptual design of intelligent mecha- tronic
systems,” Mechatronics, 2003.
2. W. Schafer and H. Wehrheim, “The challenges of build- ing advanced
mechatronic systems,” Future of Software Engineering, International
Conference on Software Engi- neering, IEEE Computer Society, 2007.
3. Philipp Limbourg, “Dependability modelling under un- certainty: An
imprecise probabilistic approach,” Springer, 2008.
4. K. Thramboulidis, “Model integrated mechatronics: To- wards a new
paradigm in the development of manufac- turing systems,” IEEE
Transactions on Industrial Infor- matics, Vol. 1, No. 1, February 2005.
5. OMG, “OMG Systems Modeling Language (OMG SysML™),” V1.0,
September 2007.
6. GD250, “Lifecycle process model ‘V-Model’,” Available online: http://
www.informatik.uni-bremen.de/gd-pa/vmodel/ vm1. htm#application
7. K. Thramboulidis, G. Doukas, and G. Koumoutsos, “A SOA-
based embedded systems development environment for industrial
automation,” EURASIP Journal on Embed- ded Systems, Article ID
312671, pp. 15, 2008.
8. K. Thramboulidis, “Challenges in the development of mechatronic
systems: The mechatronic component,” 13th IEEE International
Conference on Emerging Technolo- gies and Factory Automation
(ETFA), Hamburg, Ger- many, September 2008.
9. OMG, “Unified modeling language: Superstructure,” version 2.1.1,
formal/2007-02-03.
10. M. Habib, “Mechatronics,” IEEE Industrial Electronics Magazine,
Vol. 1, No. 2, Summer 2007.
11. S. Burmester, H. Giese, and M. Tichy, “Model-driven development of
reconfigurable mechatronic systems with mechatronic ‘UML’ in model
driven architecture,” Springer Berlin/Heidelberg, Vol. 3599, 2005.
12. K. Knorr, A. Lapp, P. Torre Flores, J. Schirmer, D. Kraft
13. J. Petersen, M. Bourhaleb, and T. Bertram “A process model for
distributed development of networked mecha- tronic components
in motor vehicles,” Proceedings of the IEEE Joint International
Conference on Requirements Engineering (RE’02), 2002
14. R. Nordmann, “Use of mechatronic components in rotat- ing machinery,”
The 3+1 SysML View-Model in Model Integrated Mechatronics 63
4
Design and Development of
Mechatronic Application in Agricultural
Irrigation Device
ABSTRACT
In order to help the small-scale farmer, an automatic irrigation control
system was proposed. This system will provide an irrigation system that
will ease the burden of the citizen to take care of the plant. This system will
run automatically by referring to the time set by the user. As the name itself
Citation:- Razali, M. , Roslan, S. , Halim, A. and Basit, H. (2016), “Design and Devel-
opment of Mechatronic Application in Agricultural Irrigation Device”. World Journal
of Engineering and Technology, 4, 450-459. doi: 10.4236/wjet.2016.43045.
Copyright: © 2016 by authors and Scientific Research Publishing Inc. This work is li-
censed under the Creative Commons Attribution International License (CC BY). http://
creativecommons.org/licenses/by/4.0
66 Mechatronics and Robotics
is a water control system, this system will only start irrigating when the time
set triggered the water control level for the plant to grow healthily. It will
automatically stop when the timer is off (1 hour). The brain of the system is
the PLC (Programmable Logic Controller). This is the place where all the
activities are done. The irrigation will be provided by a pump that is also
connected to the microcontroller. The pump will be activated until the timer
has reached its time set. This system will continue running until the user
presses the OFF button.
Keywords: Irrigation System, Water Control System, PLC
(Programmable Logic Controller)
INTRODUCTION
As technology undergoes rapid advancement, complicated task such as
control system is accomplished with a highly automated control system.
Example of the system may be in the form of Programmable Logic Controller
(PLC) and possibly a host computer and accompanied with signal interfacing
to the field devices such as opera- tor panel, motors, sensors, switches, and
solenoid valves. Network communication is capable of providing a large
scale implementation and process coordination besides enabling greater
versatility in understanding distributed control system. Every single integral
component module in a control system plays a major role regardless of size.
Generally, PLC can be explained as a digital electronic device that uses a
programmable memory to store direction and to accomplish functions such
as counting, logic sequencing, arithmetic and timing in order to control
machine and processes. The term logic was used because the programming
is basically concerned with implementing logic and switching operations
[1]. Among the most important implementation for PLCs in Southeast Asia
is the oil & gas industry. For example, Malaysia is world famous for its
palm oil, accounting for a large portion of total global production. Here,
many mid and high-end PLCs are used in the production of palm oil, fed
mostly by government funding. Meanwhile, large PLCs are used in the oil
and natural gas industry in Indonesia and Malaysia. Another important for
PLCs implementations in Southeast Asia is the automotive manufacturing
industry. Thailand is world-famous for its production of automotive
electronics; with industry requirements for PLCs multiply after the country
was hit by critical flooding in 2011. Malaysia has a big automotive industry
markets in the area among the locals, where a very large population will
soon requesting more favorable modes of transport. A third application is
Design and Development of Mechatronic Application in ... 67
Problem Statement
The world’s population is estimated to reach 9.2 billion by 2050.The UN
Food and Agriculture Organization (FAO) has projected that farmers will
need to increase their production by 70% more food than in 2006 to meet
this demand.
Inability to do so will causing in food scarcity and poorer healthcare
in developing countries, with damaging impact for development and
the potential for conflict within and between nations. While agriculture
productivity has been increasing, production capacity is stagnant, and food
security remains a major issue in many nations due to hiking prices as well
as availability. Continuing increase costs of input of agriculture, commodity
speculation and competition with other uses for crops. In poorer nations, this
causing the number of people considered to be ‘food insecure’ increasing
because they are not able to afford enough food even if it were available
to buy. Both productivity and absolute production need to increase if this
issue is to be addressed. In Malaysia, small farms remain at the Centre of
agriculture and rural development. However, one of the main causes for the
low agricultural productivity in most developing countries in the region is
the lack of appropriate machineries that cater to and suit the requirements
of small-scale farms. For this reason, many small farms are deemed as
unproductive and inefficient. Farm mechanization plays a significant role
in every nation’s economy. However, it is often misconstrued to mean
modernization, beneficial only to industrialized countries with highly
mechanized agriculture. Developing countries often have to rely on a
variety of imported farm machines, which are seldom appropriate for small
farms [2]. One of the issues that need to be address is the lack of usage of
68 Mechatronics and Robotics
Significant of Study
Automation and robotic can play a large role in society to meet the
agricultural production needs. For decades automation and robots have
played a core role in multiplying the efficiency and lowering the cost of
industrial production and products. In the past several years, a similar shift
has started to take place in agriculture industry, with GPS- and vision-based
such as self-guided tractors and harvesters already being available in the
market commercially.
More recently, developers have started to test with autonomous systems
that combine operations such as thinning, pruning and harvesting, as well
as sowing, spraying, mowing and weed removal. In the fruit production
industry, for example, robotic platforms rode by workers have shown to be
doubling the efficiency compared to workers using ladders. Advancement
in sensors technology and control systems allow for optimum resource and
integrated pest and disease management.
This will be a revolution in the way that food nowadays is grown,
tended, and harvested. The level of input of engineering technologies into
agriculture is generally still low [3]. Automation component such as PLC is
an important module in agricultural automation.
It can be found used in irrigation systems, food processing, building
environmental control, grain drying, aquaculture production, and tractor and
machinery systems. Using automation such as PLC enables us to increase the
number of production by reducing error such as human and environmental
factor (damaged produce).
Objective of Study
• To develop the hardware and software by using PLC
(Programmable Logic Controller) as a main controller.
• To interface PLC module with input and output component.
• To design prototype system for motion control.
Design and Development of Mechatronic Application in ... 69
LITERATURE REVIEW
Input Device
Intelligence of an automated system is greatly depending on the ability of
a PLC to read in the signal from various types of automatic sensing and
manually input field devices. Toggle switches, push buttons, and keypad,
which from the basic man machine interface, are types of manual input
device. On the other hand, for detection of workplace, monitoring of moving
mechanism, checking on pressure and or liquid level and many others, the
Output Devices
An automatic system is incomplete and the PLC system is virtually paralyzed
without means of interface to the field output devices. Some of the most
commonly controlled device is relay indicator, motors, solenoids, buzzers
Design and Development of Mechatronic Application in ... 71
and so on (Figure 3). Through activation of motors and solenoids the PLC
can control from simple pick and place system to a much complex servo
positioning system. These type of output devices are the mechanism of an
automated system and so its direct effect on the system performance.
Ladder Logic
To ease the use of PLC’s we programmed it using ladder logic format. Ladder
logic is a visual representation of a set of inputs and outputs. Ladder Logic
format resembles familiar hardware systems and it doesn’t require extra
training for technicians, engineers and also those who without backgrounds
in these hardware systems. The name is derived from the fact that the diagram
upon completion resembles a ladder. A Ladder Logic diagram consists of a
vertical line on the left hand side, known as the hot rail, and a vertical line
on the right hand side, known as the neutral rail [6] . They are connected
by lines, known as rungs; with several different symbols each represents an
input or output. The logic is following the path and determining if the rung
is true. For the rung to be true it is necessary to be able to flow across closed
contacts to the opposite rail. If the rung is true the output is then considered
true. Figure 4 is sample ladder logic diagram.
Relay
Relay is one of the switch types which are electrically operated. Most
of the relay uses an electromagnet to operate a switching mechanism
mechanically. It is used to control a circuit by a low-power signal (with
complete electrical isolation between control and controlled circuits), where
several circuits must be controlled with one signal. A contactor is one of
the relay types that can handle a high power that required controlling an
electric motor electrical circuit from overload or faults; in modern electric
power systems these functions are performed by digital instruments still
called “protective relays”. In general, components in relay are inductor coil,
a spring (not shown in figure), swing terminal, and two high power contacts
named as normally closed (NC) and normally open (NO). Relay uses an
Electromagnet to move swing terminal between two contacts (NO and NC).
When there is no power applied to the inductor coil (Relay is OFF), the
spring holds the swing terminal is attached to NC contact. Figure 5 shows
the pin diagram for relay:
72 Mechatronics and Robotics
The relay will start operating when the required power is applied to the
inductor coil will generates a magnetic field which will move the swing
terminal from normally close contact to normally open contact. When the
power is OFF, the spring will move back to normally close contact.
Materials
1) PLC training kit
• ・ PLC microcontroller (CP1E-N30DR-D)
• ・ CX one PLC programmer
• ・ Stepper motor
• ・ Touch-screen teach pendant
74 Mechatronics and Robotics
Methodology
Figure 6 showed the methodology or work flow of the project that has been
used as the guideline in order to do the project. In started with investigate
the topic and objectives with supervisor. After doing literature review,
equipment that needed was investigated like relay, water pump motor, PLC
and AC/DC converter.
Pilot Experiment
Pilot experiment was done first so we will learn first on how to write a
ladder logic system and install it into PLC.
This experiment was done using PLC training kit as a prototype that
already had been assembled so it can be used as a learning tool before actual
experiment.
The programming that we are going to use is CX one programmer that
can be used to write ladder logic system. First, we need to connect PLC to a
computer so we will be able to install the programing into our PLC. We put
the logic:
When push button is ON, turn the stepper motor ON. After that, we
write this logic into a ladder logic system (CX one).
Then we load this program into the PLC. We connect the sensor input to
the PLC, the one we specified in our program (push button).
Lastly, we connect the PLC External Output Terminal (specified by our
program) to the stepper motor. Now, we execute the logic program on the
PLC.
Design and Development of Mechatronic Application in ... 75
Construction of Model
Firstly, we need to connect our PLC to a computer using an adapter. Then
we run the program (CX One programmer shown in Figure 7) and write the
intended logic sequence for our experiment that is:
1) When switch is ON, turn the Timer A for 100 seconds.
2) During Timer A ON, the water pump motor is ON until Timer A
is OFF.
3) When Timer A is OFF, turn ON Timer B (100 seconds).
4) When Timer B is OFF, switch ON the water pump motor and
Timer C (100 seconds).
5) When Timer c is OFF, turn ON Timer D (100 seconds).
6) When Timer D is OFF, the counter resets the main switch and the
cycle continue.
Next, we load the program that has been written into the PLC. After
that, connect the push switch to the PLC by using connecting wire, the one
we specified in our programming earlier. Then, connect the PLC External
Output Terminal to the water pump motor (output). Now, we execute the
logic program on the PLC. After the model has been completed, we do
several test run for our model to make sure the PLC programming work as it
intended. Any failure will be observed and documented so we will be able to
do troubleshooting for our model. Figure 8 shows the model of construction
software for irrigation control.
76 Mechatronics and Robotics
RESULT
Hardware Development
Figure 9 shows an overview of work table in Farm Mechanization laboratory
in Universiti Sultan Zainal Abidin campus. Figure 10 shows constructed
model for irrigation control.
Simulation
Figure 11 shows an algorithm for automatic irrigation control.
DISCUSSION
Introduction
Both hardware and software part are developed and built successfully. The
system is complete when both of them are connected together. This chapter
explains all the results related to both part.
Hardware
The microcontroller is connected directly to AC to DC adaptor that is
functioning as power supply to switch on the control unit. The moisture
sensor will start measuring and the data obtain will be sent to control unit
in voltage form that will be processed by PIC microcontroller. Since the
microcontroller has the capability of 10-bit analog to digital converter, it
will convert and manipulate the data obtained to get the actual value of
humidity. After that, these values are displayed on the LCD. The control unit
also will instruct the pump to start working.
CONCLUSION
The project involved the design activity on both software and hardware
development of PLC in agriculture mini hydroponic system. The design on
circuitry system is involved for connecting between input and output system
including push button device, microcontroller device, 240 volts power
supply and water motor pump. The software development is created from
basic concept on timer and counter application. By using this technique, this
project was successful to automate basic irrigation system for hydroponic.
The system will automatically circulate the fertilizer mix water throughout
the system on selected period of morning and afternoon daily with less
intervention of the user.
Design and Development of Mechatronic Application in ... 79
REFERENCES
1. Selvaraj, K. (2010) Development of a “Programmable Logic Controller
Circuitry” for Optimal Power Distribution in a Manufacturing
Industries. Journal of Computer Science, 6, 250-252.http://dx.doi.
org/10.3844/jcssp.2010.250.252
2. (2005) Small Farm Mechanization Systems Development, Adoption
and Utilization. Report, FFTC Annual Parts, Sri Lanka, 13-17.
3. Kamaruddin, R., Rukunuddin, I.H. and Seng, O.H. (2007) Research and
Development of Agricultural Engineering in Malaysia. Paper, Country
Asian, United Nations Centre, Pacific Engineering, Agricultural, 1-11.
4. Halim, M.Z.B. (2007) The Temperature Control System Using PLC.
5. Ismail, M.H. (2008) Rain Water PLC Based Detector and Valve
Switcher.
6. Schumann, D., Fietsam, J., Rochel, E., Dixon, K., Imbayan, M. and
Ibayan, M. (n.d.) Literature Review for the Design of a SCADA
System.
7. Shah, C. (2004) Sensor Less Control of Stepper Motor Using Kalman
Filter. Department of Electrical and Computer Engineering, Cleveland
State University, Master of Science in Electrical Engineering.
CHAPTER
5
Design and Development of a
Competitive Low-Cost Robot Arm with
Four Degrees of Freedom
ABSTRACT
The main focus of this work was to design, develop and implementation of
competitively robot arm with enhanced control and stumpy cost. The robot
arm was designed with four degrees of freedom and talented to accomplish
INTRODUCTION
The term robotics is practically defined as the study, design and use of robot
systems for manufacturing [1]. Robots are generally used to perform unsafe,
hazardous, highly repetitive, and unpleasant tasks. They have many different
functions such as material handling, assembly, arc welding, resistance
welding, machine tool load and unload functions, painting, spraying, etc.
There are mainly two different kinds of robots: a service robot and
an industrial robotic. Service robot is a robot that operates semi or fully
autonomously to perform services useful to the well-being of humans and
equipment, excluding manufacturing operations [2]. Industrial robot, on
the other hand, is officially defined by ISO as an automatically controlled
and multipurpose manipulator programmable in three or more axis [1].
Industrial robots are designed to move material, parts, tools, or specialized
devices through variable programmed motions to perform a variety of tasks.
An industrial robot system includes not only industrial robots but also any
devices and/or sensors required for the robot to perform its tasks as well as
sequencing or monitoring communication interfaces.
In 2007 the world market grew by 3% with approximately 114,000 new
installed industrial robots. At the end of 2007 there were around one million
industrial robots in use, compared with an estimated 50,000 service robots
for industrial use [3].
Due to increase using of industrial robot arms, an evolution to that topic
began trying to imitate human movements in a detail mode. For example a
group of students in Korea made a design of innovations that robotic arm
Design and Development of a Competitive Low-Cost Robot Arm with ... 83
MECHANICAL DESIGN
The mechanical design of the robot arm is based on a robot manipulator
with similar functions to a human arm [6-8]. The links of such a manipulator
are connected by joints allowing rotational motion and the links of the
manipulator is considered to form a kinematic chain. The business end of
the kinematic chain of the manipulator is called the end effector or end-of-
arm-tooling and it is analogous to the human hand. Figure 1 shows the Free
Body Diagram for mechanical design of the robotic arm. As shown, the
end effector is not included in the design because a commercially available
gripper is used.
84 Mechatronics and Robotics
(1)
(2)
(3)
(4)
(5)
(6)
(7)
(8)
The servo motor that was selected, based on the calculations, is the
Hextronik HX12K, which has a torque of 280 oz/in. This motor was
recommended because it is much cheaper than any other motor with same
specifications. Since we need more torque at joint B, see Equation (8),
we used two motors at point B to comply with the torque requirements;
however, one motor is enough for the other joints. Using two motors at joint
B is much cheaper than using one big motor with 560 oz/in. Other relevant
Design and Development of a Competitive Low-Cost Robot Arm with ... 87
characteristics of the motors, which can be shown in Figure 5, are that they
can turn 60 degrees in 130 milliseconds and they have a weight of 47.9
grams each.
Once the initial dimensions for the robot arm and the motor were defined,
the design were carried out using the SolidWorks platform; design should
carefully take into account the thickness of the acrylic sheet and the way that
the pieces would be attached to each other. The acrylic sheet used to make
the robot is 1/8 thickness and that thin sheet was chosen because it easier for
machining and less weight with a good resistance.
There were critical points that did not resist the fastening and, in turn, may
break down; hence, reinforcements in these points were considered. The
final result of the robot arm is shown in Figure 7.
positions the final element in the y axis, the motors B and C positions the
final element in the x and z axis.
The problem was simplified by using the xz plane, as shown in Figure 9.
In which the following known values were defined [9]:
LAB: the forearm length.
LBC: the arm length.
z: the position in the z axis.
x: the position in the x axis. y: the position in the y axis.
Using trigonometry relations, as shown in Figure 9, the motor angles θ2
and θ1 are obtained, as seen in Equations (9) and (10).
` (9)
(10)
(11)
The motor B is going to use θ1 and the motor C is going to use θ2. The
angle for the motor A is calculated as seen in Equation (11).
Figure 9. xz Plane.
With these calculations, the angles of servomotors are obtained and in
turn they take the action to move the whole structure to the specific position.
END-EFFECTOR SELECTION
The end effector is probably one of the most important and most complex
parts of the system. Wisely, it is much easier and economical to use a
commercial one than build it. The end effector varies mainly according to
the application and the task that the robot arm accomplishes for; it can be
pneumatic, electric or hydraulic. Since our robot arm is based in an electric
system, we may choose electric basis of end effector. Besides, the main
application of our system is handling, accordingly, the recommended type of
our end effector is a gripper, as shown in Figure 10. Please note that the end
effector is controlled by a servo motor and, in turn, the total servo motors
used for our robot arm will be 5 motors that move the structure.
Manual Control
This type of control is an extra option for our system that useful in specific
positions. In case of mandatory positions that the inverse kinematics mode
cannot calculate their valid angles, we may use the manual control instead.
Basically, manual control consists of a series of analog inputs, such as
potentiometers, that are connected with the microcontroller which will
interpret the values and send a command to the servo driver. In order to
implement this, a control board, as shown in Figure 16, should be built
to work as an interface with the user. Possible implementation includes a
teaching feature where the microcontroller stores positions in memory and
by a keypad or a series of switches we may recall these positions.
Current Consumption
The current consumption depends on the load and the type of motion
of the robotic arm. In the current study, there are 4 levels of current
consumptions:
• Low (from 0 to 200 mA). This consumption takes place when the
robot is at rest (not motion case).
• Normal (from 200 to 500 mA). This happened when the robot
arm is moving with capability to go to the target without needs
of great torque.
• High (from 500 mA to 900 mA). This range is reached at the
beginning of carrying loads. By overcoming the initial moment
of inertia for loads, the normal range takes a place.
Over current (more than 900 mA). The load is too heavy and the motor
cannot move at all. For being under this condition for more than one minute,
the motor will burn, i.e. it is not possible to be used any more.
Maximum Load
These results were obtained using different weights; a bag of corn was used
with a scale to determine bag’ weight. Results carried out by using the robot
arm to pick up the bag and move it to specific positions. Table 2 presents
the current consumption at different weights of bag of corn. From Table 2,
it can be seen that the robot can move without problems at loads lower than
50 grams. At loads 60 grams, the robot arm start having difficulties and
after passing 80 grams severe condition occurred where irreversible damage
could be happened in motors.
98 Mechatronics and Robotics
Final Position
Results show the precision of the robot arm to move different weight (<50
grams) is presented in Table 3. As shown, the robot arm is able to perform the
movement to the position specified. However, this movement is not smooth
and sometimes the motors do not have enough force, especially when the
load is heavy. In addition, some problems may appear due to synchronizing
the two bottom motors. The steps of the two motors were not coincidental
and that causes tension in the acrylic parts, which in case of being too much
will break the parts.
CONCLUSIONS
This paper presents the design, development and implementation of robot
arm, which has the talent to accomplish simple tasks, such as light material
handling. The robot arm was designed and built from acrylic material where
servo motors were used to perform links between arms and execute arm
movements. The servo motors include encoder so that no controller was
implemented; however, the rotation range of the motor is less than 180º
span, which greatly decreases the region reached by the arm and the possible
positions. The design of the robot arm was limited to four degrees of freedom
since this design allows most of the necessary movements and keeps the
Design and Development of a Competitive Low-Cost Robot Arm with ... 99
costs and the complexity of the robot competitively. The end effector is not
included in the design because a commercially available gripper is used
since it is much easier and economical to use a commercial one than build it.
During design, we faced some difficulties due to the way of joining thin
acrylic parts strongly. A mechanical junction based on screws and nuts is
used and in order to accomplish that, a small feature was designed which
allowed fastening the bolts with the nuts without having to screw in the thin
acrylic layer.
To control the robot arm, three approaches are implemented: a
microcontroller, a driver, and a computer-based user interface. This system
has unique characteristics that allow flexibility in programming and
controlling method, which was implemented using inverse kinematics;
besides it could also be implemented in a full manual mode. This robotic
arm is contrast with others as being much cheaper than available robot
arms, also it can be controlled all of its movements from a computer, using
a Labview interface.
Several tests were carried out to validate the robot arm where the testes
covered both the particular elements and the overall system; results at
different operating conditions show trustful of the robot arm presented.
100 Mechatronics and Robotics
REFERENCES
1. Manipulating Industrial Robots—Vocabulary, International
Organization for Standardization Standard 8373, 1994.
2. Industrial and Service Robots, IFR International Federation of
Robotics, 2010. http://www.ifr.org/home
3. Case Studies and Profitability of Robot Investment, The IFR Statistical
Department, 2008. http://www.ifrstat.org/downloads/2008_Pressinfo_
english .pdf
4. R. J. Wang, J. W. Zhang, et al., “The Multiple-Function Intelligent
Robotic Arms,” FUZZ-IEEE Journal, Korea, 20-24 August 2009, pp.
1995-2000.
5. L. B. Duc, M. Syaifuddin, et al., “Designing 8 Degrees of Freedom
Humanoid Robotic Arm,” International Conference on Intelligent and
Advanced Systems, Kuala Lumpur, 25-28 November 2007, pp. 1069-
1074.
6. C. R. Carignan, G. G. Gefke and B. J. Roberts, “Intro to Space Mission
Design: Space Robotics,” Seminar of Space Robotics, University of
Maryland, Baltimore, 26 March 2002.
7. Occupational Safety and Health Administration Technical Manual,
OSHA 3167, United States Department of Labor, 1970.
8. B. Siciliano, L. Sciavicco, L. Villani and G. Oriolo, “Robotics,
Modelling, Planning and Control,” Springer, London, 2009.
9. M. P. Groover and M. Weiss, “Robotica Industrial, Technologia,
Programacion y Aplicaciones,” Mc-Graw Hills, Mexico D.F., 1989.
CHAPTER
6
An Adaptive Robust Approach to
Modeling and Control of Flexible
Arm Robots
ABSTRACT
In this paper, a novel adaptive robust approach to modeling and control of
a class of flexible-arm robots subject to actuators unmodeled dynamics is
proposed. It is shown how real-time signals measured from a dynamical
Citation:- Bagheri, P., Behjat, L. and Sun, Q. (2018), “An Adaptive Robust Approach
to Modeling and Control of Flexible Arm Robots”. World Journal of Engineering and
Technology, 6, 282-295. doi: 10.4236/wjet.2018.62017.
Copyright: © 2018 by authors and Scientific Research Publishing Inc. This work is li-
censed under the Creative Commons Attribution International License (CC BY). http://
creativecommons.org/licenses/by/4.0
102 Mechatronics and Robotics
INTRODUCTION
The area of flexible-arm robots has attracted much attention during the last
few decades [1] [2]. This interest is due to the advantages that flexible arms
offer compared to their rigid counterparts. Weight reduction, lower energy
consumption, and faster system response are among several benefits utilized
in their numerous applications such as space missions [3].
There has been a great number of studies coping with controlling
flexible-arm robots, many of which investigate both theoretical and
experimental aspects in this field [4] [5]. On grounds of the flexibility of
the arms, along with a trajectory-tracking control problem, vibration control
should be also considered for such systems to improve the control system
performance [6]. As a result of the vibrations caused by the flexibility of
the arms, designing controllers for such systems becomes a challenging
task. There exist several research works in the literature addressing the
flexibility of the arm. Passive control methods are one way to deal with the
vibrations of elastic arms which require modification in physical parameters
of the system structure [7]. Due to these structural modifications, absorption
properties of the structure can be employed to increase damping properties
of the arms. Active control approaches have also been widely used to control
An Adaptive Robust Approach to Modeling and Control of ... 103
Feedback linearization methods are wildly used for nonlinear systems. Nev-
ertheless, As a result of the underactuated nature of flexible arms, the exact
feedback linearization method [21] cannot be employed for such mechanical
systems. Hence, the partial feedback linearization method [22] is considered
as one the most suitable way to cope with underactuated systems. That is,
the entire system can be either linearized with respect to the active degrees
of freedom (collocated problem) or with respect to the passive degrees of
freedom (non-collocated problem). Be that as it may, such a control strategy
fails to handle uncertainties as control inputs to feedback-linearized systems
depend on the governing equations.
In this paper, a novel adaptive robust nonlinear controller is designed for
a class of multi-link flexible arms subject to uncertainties and unmodeled
dynamics of the actuators. The contribution of the current research work
is to utilize the system’s real-time responses to improve the accuracy of
the available mathematical model. Additionally, despite uncertainties and
underactuated nature of the system, the proposed controller is able to
track desired trajectories asymptotically. That is, despite the presence of
uncertainties and disturbance sources, the tracking error converges to zero.
To do such, first, an adaptive robust controller is designed for the active
degrees of freedom. Meanwhile, adaptive laws are proposed to estimate the
system uncertainties as well as approximating actuators uncertain dynamics.
Furthermore, another adaptive robust controller is designed to control the
passive degrees of freedom where the stabilities of the overall closed-loop
system, in both cases, in the presence of uncertainties are established. In
addition, due to the flexibility of the arms, the extraction of the system
dynamics is demanding and, for the most part, should be done analytically.
In the present approach, two adaptive signals are introduced and synthesized
such that the cumbersome analytical part of the model extraction can be
done numerically.
DYNAMIC EQUATIONS
In the current study, a multi-link flexible-arm robot is considered for
modeling and control purposes. The investigated system is assumed to
include m links containing l elastic modes. Therefore, the total degrees of
freedom come to n=m+l.
Equations governing a flexible arm using the Euler-Bernoulli beam can
be written as follows
An Adaptive Robust Approach to Modeling and Control of ... 105
(1)
(2)
(3)
where U denotes potential energy of the system. Extracting the dynamic
equations of flexible robot arms is time consuming and can be analytically
complicated. Due to its analytical nature, there exists a high chance of
computational errors as the equations are derived. The extraction of the
dynamic equations can be divided into two phases: 1) Extracting the portions
of the dynamics that are straightforward and demand less calculations such
as kinetic and potential energy. It is to be noted that the inertia matrix
M(q) can be arrived at directly from the rearrangement of the kinetic energy
of the system. 2) Obtaining the portions of the dynamics which require
lengthy analytical calculations such as G(q) and . It is worth noting
that calculating G(q) and require calculus of variations which can be
cumbersome for systems with a large degrees of freedom. The idea of this
paper is to propose a way to procure G(q) and numerically using the
system responses. For this purpose, (1) is rewritten as
(4)
where
(5)
(6)
106 Mechatronics and Robotics
and d(t) is an unknown vector added to the equations to account for uncertain
dynamics of the actuators. It is also assumed that where γ is an
unknown constant.
In order to complete the dynamic equations given in (4), and
should be computed. As it follows from (5) and (6), and
are expressed in terms of the matrix M, system states, and its potential
energy. Therefore, one can use the measurement of the states along with the
kinetic and potential energy of the system and use numerical computations to
obtain h1 and h2. Hence, the system dynamics is computed as it runs without
going through analytical calculations. However, such dynamic equations are
prone to numerical and measurement errors. Thus, a special measure should
be taken to improve the accuracy of such computations. To compensate for
such errors, two unknown auxiliary signals P and S are employed as follows
(7)
(8)
where are two unknown auxiliary signals
introduced to account for numerical and measurement errors;
are two diagonal matrices constructed as
and denotes the
jth element of the vector hi where i=1,2 . Given the physical and geometric
properties of flexible arm robots, it is not restrictive to assume that
both H1 and H2, as defined above, have full ranks. Therefore, their column
space can be utilized as a basis for the vectors h1 and h2 at any given time.
Thus, for any exact value of h1 and h2 at any given time, there exist two
unique vectors P and S such that h1=H1P and h2=H2S. It is clear that the
vectors P and S change in values over time. Hence, their values need to be
updated according to the system response. Therefore, an adaptive procedure
is developed for each of these vectors along with designing a controller in
Section 3. Consequently, as the system runs, its response is used to update
the values of P and S.
CONTROL DESIGN
In this section, an adaptive controller is designed for the considered system
in the presence of uncertainties. Alongside the adaptive controller, two
adaptive laws are derived for updating the values of P and S which are used
An Adaptive Robust Approach to Modeling and Control of ... 107
(9)
Two matrices N and R are defined as .
Consequently, (9) comes to
(10)
As N is positive definite, (10) can be expressed as follows
(11)
where , and . Tracking error and its dynamics
can be defined as , respectively; q2d is the desired
trajectory for the active degrees of freedom. A measure of the tracking error
is also defined as
(12)
where is a positive definite matrix. The estimation errors of the
uncertainty vectors are shown by which are defined as and
are the estimations of the uncertainty vectors. Using the
definitions of (11) is rewritten as follows
(13)
Considering (12), ˙ is calculated as follows
(14)
To proceed with the control design, the following Lyapunov function is
proposed
108 Mechatronics and Robotics
(15)
where is the Euclidean norm of the estimation error of the
uncertain dynamics of the actuators. Taking time-derivative of (15) gives
(16)
where is decomposed into . Given the structure of (16), the
following control and adaptive laws are proposed
(17)
(18)
(19)
where is a positive definite matrix; (17) is the proposed control
law, and (18)-(19) are the adaptive laws to update the two auxiliary signals
introduced in (7)-(8). Substituting (17)-(19) into (16) yields
(20)
(20) can be rewritten as follows
(21)
where ∥ . ∥ stands for the Euclidean norm. Since N is a bounded positive
definite matrix, the adaptive law to update is proposed as
(22)
Employing (22), (21) becomes
(23)
To establish the asymptotic stability of the closed-loop system, the
second part of the control law is introduced as follows in (24). It is to be
noted that the second part of the control law compensates for the actuators
An Adaptive Robust Approach to Modeling and Control of ... 109
uncertain dynamics.
(24)
where is an arbitrary signal where . By substituting (24)
into (23), the following is obtained
(25)
Since are positive, the following inequality is concluded
(26)
Hence, (25) is simplified as follows
(27)
Taking integral from both sides of (27) gives
(28)
Considering the fact that , the following is achieved
(29)
where a is a positive constant. It is observed that the left side of (29) is
non-negative, thus, it follows that (V V (0) − ∞( )) is bounded. As a result, it
is deduced that the right side of (29) remains bounded as well. Consequently,
QL∈2. It is also easy to check that . Therefore, it is concluded from
Barbalat’s lemma that which implies . Thus, it was proven
that adopting the extracted control laws (17), (24), and adaptive laws (18),
(19)) and (22) guarantees the asymptotic stability of the closed-loop system.
That is, the system is able to track desired trajectories asymptotically in the
presence of model and actuators uncertainties.
110 Mechatronics and Robotics
(30)
Substituting (30) into (7) yields
(31)
where . It follows from the structure of M that
M22 is nonsingular. In addition, it can be easily shown that F is a positive
definite matrix as well. Therefore (31) can be expressed as follows
(32)
where , and J is a control signal yet to
be defined. The actual control input can be retrieved from J as follows
(33)
It is concluded from (33) that G should be nonsingular. Thus, according
to the definition of G, M12 needs to be invertible as well. Hence, in order to
extract a non-singular control law for the passive degrees of freedom, the
following assumption is made:
Assumption 1 It is assumed that the number of passive degrees of
freedom is less or equal than the number of active degrees of freedom.
The tracking error for the passive degrees of freedom and its dynamics
are defined as where q1d is the desired trajectory
for the passive degrees of freedom. Moreover, a measure of the tracking
error is defined as is a positive definite matrix.
Using the definitions of which are defined in the previous section, (32)
is expressed as follows
(34)
Based on the structure of (34), the following Lyapunov function is introduced
(35)
The time-derivative of (35) is calculated as
An Adaptive Robust Approach to Modeling and Control of ... 111
(36)
where . Given (36), the following control and adaptive laws
are proposed
(37)
(38)
(39)
Where is a positive definite matrix? Substituting (37)-(39) into
(36) and taking the Euclidean norm form both sides lead to the following
(40)
To make (40) negative semi-definite, the following control and adaptive
laws are proposed
(41)
(42)
where σ(t) is a positive arbitrary signal belonging to L1 . Substituting
(40)-(42) into (40) gives
(43)
Similar to the stability analysis presented in the previous section, it is
concluded from (43) that , which implies that . Hence, it is
guaranteed that the tracking error converges to zero asymptotically in the
presence of uncertainties.
SIMULATION RESULTS
In this section, a simulation is conducted to validate the analytical results.
The simulated case study is a one-link flexible-arm robot with a payload at
the tip of the arm. It is assumed that the system is subject to uncertainties
112 Mechatronics and Robotics
and unmodeled dynamics in the actuators. The strategy is, first, to control the
active degrees of freedom to their desired values. Then, the control law and
adaptive laws switch to the non-collocated control to dampen the vibrations
of the flexible arm. The physical properties of the investigated arm are as
follows: Arm s Length = 1 m; Density = 7850 kg/m3; Beams section area
= 2e−4 m2; Elasticity = 2.07e11; Load mass = 0.1 kg. The actual values
uncertainty signals are S=8.2 and P=5.5.
It is shown in Figure 1 that the flexible arms hub angle is successfully
controlled to the desired angle. It is also observed from Figure 1 that the
hubs angle remains in its desired angle after the controller switches.
In Figure 2, the vibrations of the flexible arm are displayed, and its
zoomed-in version is depicted in Figure 3. It follows from Figure 3 that the
arm residual vibrations are dampened as the controller switches.
Figure 4 and Figure 5 correspond to the estimation of the introduced
auxiliary signals to compensate for uncertainties. As noticed from the
figures, the estimation of the uncertain parameters have converged to their
limits after about two seconds through the simulation.
Finally, the required actuator torque to control the arm is shown in Figure
6. It follows from the simulations results that the arm tip can be controlled
to its desired position more accurately. The reason stems from the fact that
the non-collocated control helps decrease the response time. It is because
that it would dramatically take a longer time for the residual vibrations to
dampen naturally.It was shown both analytically, in Section 3, and through
simulation that the proposed controller can achieve asymptotic tracking in
the presence of unknown and bounded disturbance sources.
a desired trajectory, the output can stay in the vicinity of the trajectory.
However, the convergence of the output to the desired trajectory is not
guaranteed.
CONCLUSION
As presented in detail, an adaptive robust nonlinear controller was designed
for a class of flexible-arm robots based on the Lyapunov stability theory.
A new numerical technique was presented to facilitate the extraction
of the governing equations of flexible-arm robots. It was shown how
computational errors caused by numerical operations could be compensated
and approximated as uncertainty signals. In the control design process,
first, an adaptive controller was designed for the active degrees of freedom
to enable the system to follow desired trajectories in the presence of
uncertainties. Further, it was shown that under some feasible conditions,
another adaptive robust controller could be designed for the passive degrees
of freedom. Therefore, as the system was being controlled, its responses
were utilized to improve the accuracy of the mathematical model. Thus, the
presented method leads to improvements in the controller performance and
increases tracking accuracy. Furthermore, due to a numerical approach to
extracting the dynamics, the amount of analytical computations is reduced
dramatically.
116 Mechatronics and Robotics
REFERENCES
1. Giorgio, I., Della Corte, A. and Del Vescovo, D. (2017) Modelling Exible
Multi-Link Robots for Vibration Control: Numerical Simulations and
Real-Time Experiments. Mathematics and Mechanics of Solids, 21,
1234-1246.
2. Kiang, C.T., Spowage, A. and Yoong, C.K. (2015) Review of Control
and Sensor System of Exible Manipulator. Journal of Intelligent &
Robotic Systems, 77, 187-213. https://doi.org/10.1007/s10846-014-
0071-4
3. Dwivedy, S.K. and Eberhard, P. (2006) Dynamic Analysis of Exible
Manipulators, a Literature Review. Mechanism and Machine Theory,
41, 749-777. https://doi.org/10.1016/j.mechmachtheory.2006.01.014
4. He, W., Ouyang, Y. and Hong, J. (2017) Vibration Control of a
Exible Robotic Manipulator in the Presence of Input Deadzone.
IEEE Transactions on Industrial Informatics, 13, 48-59. https://doi.
org/10.1109/TII.2016.2608739
5. Kumar, A., Pathak, P.M. and Sukavanam, N. (2011) Reduced Model
Based Control of Two Link Exible Space Robot. Intelligent Control
and Automation, 2, 112-120. https://doi.org/10.4236/ica.2011.22013
6. Nani, V.-M. and Cires, I. (2014) Vibration Control for Two Air
Compressor Units Used in Tandem. World Journal of Engineering and
Technology, 2, 314-321. https://doi.org/10.4236/wjet.2014.24033
7. Sales, T., Rade, D. and De Souza, L. (2013) Passive Vibration Control of
Exible Spacecraft Using Shunted Piezoelectric Transducers. Aerospace
Science and Technology, 29, 403-412. https://doi.org/10.1016/j.
ast.2013.05.001
8. Bai, K., Luo, M., Li, T., Wu, J., Yang, L., Liu, M. and Jiang, G. (2016)
Active Vibration Adaptive Fuzzy Backstepping Control of a 7-Dof
Dual-Arm of Humanoid Robot with Input Saturation. Journal of
Intelligent & Fuzzy Systems, 31, 2949-2957. https://doi.org/10.3233/
JIFS-169179
9. Zhang, Q., Li, C., Zhang, J. and Zhang, J. (2017) Smooth Adaptive
Sliding Mode Vibration Control of a Exible Parallel Manipulator
with Multiple Smart Linkages in Modal Space. Journal of Sound and
Vibration, 411, 1-19. https://doi.org/10.1016/j.jsv.2017.08.052
10. Mottershead, J.E. and Ram, Y.M. (2006) Inverse Eigenvalue Problems
in Vibration Absorption: Passive Modification and Active Control.
An Adaptive Robust Approach to Modeling and Control of ... 117
7
Interactive Heuristic D* Path Planning
Solution Based on PSO for Two-Link
Robotic Arm in Dynamic Environment
ABSTRACT
This paper is devoted to find an intelligent and safe path for two-link robotic
arm in dynamic environment. This paper focuses on computational part of
motion planning in completely changing dynamic environment at every
motion sample domains, since the local minima and sharp edges are the most
common problems in all path planning algorithms. In addition, finding a path
solution in a dynamic environment represents a challenge for the robotics
INTRODUCTION
The words “path planning” refer to find a path solution in robotics where
the process is totally or sub automated. It is used to indicate the type of
computational process for moving a robot from one place to another with
respect to obstacles. In other words, it represents the search for initial
feasible path which is considered as the first step to solve the path planning
problem and includes finding a path for a robot that must move from initial
given start point to the goal point which is also given as the destination
position. The environment usually contains a fully de-fined set of obstacles.
Therefore, the robot does not collide with any of the obstacles [1]. Generally,
obstacles cannot always be static. When the obstacle changes its location
over time, it called dynamic. Dynamic environment represents a revolution
and important side in modern automation. It redounds in many important
diverse applications like air and sea traffic control negotiating, freeway
traffic, intelligent vehicles, automated assembly, and automated wheelchairs.
Because of existence the moving objects, the robot needs to make decision
quickly for avoiding possible collisions [2] [3]. There are three types of
dynamic environment. The first type is called a known dynamic environment
when all the information regarding the obstacles motions, sizes, shapes, and
locations is completely known. The second type is called a partially known
dynamic environment when not all information about the obstacles exists
at the planning time. In this state, it needs to calculate the robot motion
Interactive Heuristic D* Path Planning Solution Based on ... 121
position and time which does not used before. These modification
cases not included in the original D* theory.
3) Finding a path solution utilizing the proposed method in this
paper that use the modified heuristic D* algorithm based on the
PSO optimization technique in known dynamic environment.
4) Applying the proposed method to simulate the path solution for
the two-link robot arm as case study, which is more difficult
than the mass-point case study usually used in path planning
researches.
5) The resulting path solution from this proposed method can solve
the problem of finding an optimal path in dangerous, harsh, toxic
environments and inaccessible areas.
(1)
(2)
124 Mechatronics and Robotics
(3)
Since sin θ2 is
(4)
To solve θ2, it has to use the cos−1 function. Although it is more benefical
to use tan−1tan−1 for numerical solution, in the software implementation,
we will use the function tan−1 (b/c) which is ATAN2 (b, c). This function
has a uniform accuracy over the range of its arguments, returns a unique
value for the angle depending on the signs of b and c, and gives the correct
solution if b and/or c is zero [11].
Subsequently, the two conceivable solutions for θ2 can be written as:
(5)
At last, θ1θ1 can be solved by:
(6)
where l1 and l2 are lengths of the first and the second links, respectively.
Interactive Heuristic D* Path Planning Solution Based on ... 125
Figure 2. Proposed planning for the two-Link robotic arm in dynamic environ-
ment.
The other nodes will go to the open list with taking into consideration
the obstacles having a huge cost function and sorting the initial locations
in the closed list. With each motion samples, the closed list will be updated
to catch the new obstacle location. This means, the closed list is changing
dynamically over motion samples to save the obstacle nodes in the current
motion sample only, and the process complete like that until reaching the
start node. When the search is completed the robot will move from the start
Interactive Heuristic D* Path Planning Solution Based on ... 127
(8)
The position of each particle (the point that is chosen from D* path) is
updated by using the new velocity for that particle, according to equation
below:
(9)
where c1 and c2 are the cognitive coefficients (c1+c2≤4 ), and r1,j and r2,j
are random real numbers between [0, 1], the inertia weight w controls the
particle momentum.
128 Mechatronics and Robotics
SIMULATION RESULT
First Environment
Here we have applied the two-Link robot arm as a case study. The arm length
is equal to the half of the environment with 50 cm for each link and the arm
has 360 degrees. The obstacle motion behavior is shown in Figure 4. Since
Interactive Heuristic D* Path Planning Solution Based on ... 129
the multi FCS is applied as clarified in Figure 5, where the column (5(a))
represents the elbow down FCS analysis at some selected motion samples;
1, 4, 5, 7, and 8, while in the column (5(b)) and (5(c)) represents the elbow
up FCS analysis and the gathering between elbow up and elbow down FCS
analysis respectively at the same motion samples. In each FCS figure the
colored area (blue, red or green) is represented the free and safe area for
moving and the white area is represented the unreachable area. Therefore, it
has calculated the elbow up and elbow down FCS to enable the robot from
making decision which free space is the suitable for moving.
Figure 3. Proposed D* via Point flowchart for position and time consideration
dynamic environment.
on. Figure 6shows the two-link robot arm D* path environment where part
(a) represents the beginning of the planning while the end-effector is on the
start point, part (b) is the arm location after 8 motion sample from planning
process. At part (c) the arm is return backward at motion sample 15 to avoid
collision, part (d) and
Figure 5. FCS for the D* path and PSO path where (a) represents elbow down,
(b) represents elbow up and (c) represents gathering between elbow up and
elbow down.
Interactive Heuristic D* Path Planning Solution Based on ... 131
Second Environment
For more verification, a second more complicated environment with
four moving obstacles in such different way in comparing with the first
environment and more close to the robot links. Two from these obstacles
132 Mechatronics and Robotics
Figure 7. (a) Final D* and PSO path and (b) the path cost function.
Figure 8. Two-link robot arm result for the D* based on PSO path.
In each FCS figures the colored area (blue, red or green) indicates the
free and safe area for moving robot links while the white area indicates the
unreachable area.
Interactive Heuristic D* Path Planning Solution Based on ... 133
Figure 9. The two-link robot arm joint angles for the optimized path.
Figure 10. The change in joint angles for the optimized path.
Environment 1 2
Obstacles num- 3 4
ber
Population 50 50
iteration 100 100
D* path 72.4264 cm 78.2132 cm
PSO path 69.6865 cm 73.8408 cm
Therefore, the robot can take a proper decision for selecting the next
points in a free and intelligent path inside the permitted motion area. Figure
12 shows the two-link robot arm D* path environment where part (a)
represents the beginning of the planning while the end-effector is on the
start point, part (b) shows the arm location after 11 motion samples, part (c)
indicates that the arm is returning backward case at the motion sample 21 to
avoid collision, part (d) and
134 Mechatronics and Robotics
Figure 11. FCS for the D* path and PSO path for the second environment where
(a) represents elbow down, (b) represents elbow up and (c) represents gathering
between elbow up and elbow down.
Figure 12. The two-link robot arm result for the D* path for the second envi-
ronment.
Interactive Heuristic D* Path Planning Solution Based on ... 135
part (e) show the planning after 33 and 42 motion samples respectively,
while after 46 motion samples, the arm is returned backward again to evade
the closest obstacle before reaching the goal point safely as shown in part
(g) and (h).Because of the known position and time motion behavior of the
obstacles, the process of the path planning can be made off-line, therefore the
PSO optimization technique can plays an important role for finalizing this
process and finding the best possible shortest and optimal path to get rid of the
heuristic sharp path edges. Both of the D* and PSO paths are shown in Figure
13(a), where the blue line represents the D* path and the black line represents
the PSO path. Also, Figure 13(b) represents the optimization cost function
progress and proves that the final D* based on PSO path length was shorter by
4.3724 cm than the original path. Figure 14 confirms the path safety validation
by showing clearly all the robot configurations during motion according to
the optimized D* based on PSO path from the initial point as in part (a) to
the goal point as in part (f). In addition, Figure 15 and Figure 16 represent
the two link robot arm joint angles and the changes in these angles during
the optimized path respectively. Also, Table 2 clarifies the environment
information.In comparison with the nearest research literature [4] that try to
solve the same problem for mobile robot using the D* heuristic method, our
proposed method has different results because of our original modification
that applied on D* algorithm which guarantees the interactivity and safety of
the path and because of finalizing the path using PSO optimization technique.
CONCLUSION
This paper proposed an optimized heuristic approach for solving the
problem of two-link robot arm path planning in known position and time
changes dynamic environment. This approach utilizes the use of heuristic
D* algorithm for finding the best path and the use of PSO optimization
technique for getting the final optimal path.
Figure 13. (a) Final D* and PSO path and (b) the path cost function for the
second environment.
136 Mechatronics and Robotics
Figure 14. The two-link robot arm result for the D* based on PSO path for the
second environment.
Figure 15. The two-link robot arm joint angles for the optimized path of the
second environment.
Interactive Heuristic D* Path Planning Solution Based on ... 137
Figure 16. The change in joint angles for the optimized path of the second
environment
Also, a modification on the D* algorithm has been done by adding a
stop case and return backward case to ensure the path solution interactivity
response to the obstacles position and time continuous changes during robot
motion. From this approach, it was possible to get an intelligent path in
spite of the environmental difficulties where a completely collision free path
could be found at the expense of optimality. Therefore, the PSO optimization
technique has been used to enhance the final path by removing the sharp
edges until reaching path optimality. The result clearly proves the advantage
of using this proposed methodology since the robotic arm has preserved the
intelligence interactive behavior even after applying PSO Technique. The
length of the final path was smart, smooth, safe, and optimal, and traditional
D* algorithm could not find such a path. Our proposed method has a
new approach of solution that leads to different results since it keeps the
motion behavior of the robot to be intelligent and interactive with dynamic
environment changes. The future directions for developing this proposed
method include utilizing other intelligent methods and theories.]
Conflicts of Interest
The authors declare no conflicts of interest regarding the publication of this
paper.
138 Mechatronics and Robotics
REFERENCES
1. Miao, H., Tian, Y.-C. and Feng ,Y.M. (2013) Robot Path Planning in
Dynamic Environments using a Simulated Annealing Based Approach.
Faculty of Science and Technology, Queensland University of
Technology, Brisbane.
2. Raheem, F.A. and Bader, M.M. (2017) Development of Modified Path
Planning Algorithm Using Artificial Potential Field (APF) Based on
PSO for Factors Optimization. American Scientific Research Journal
for Engineering, Technology, and Science, 37, 316-328.
3. Fujimura, K. (1991) Motion Planning in Dynamic Environments.
Springer-Verlag, Tokyo. https://doi.org/10.1007/978-4-431-68165-6
4. Stentz, A. (1995) The Focused D* Algorithm for Real-Time
Replanning. Proceedings of the International Joint Conference on
Artificial Intelligence, Montreal, 20-25 August 1995, 1652-1659.
5. Shen, Y.S. and Huper, K. (2004) Optimal Joint Trajectory Planning for
Manipulator Robot Performing Constrained Motion Tasks. Department
of Information Engineering, Australian National University, National
ICT Australia Limited, Canberra.
6. Khelchandra, T. and Huang, J. (2015) Robot Motion Planning
with Neuro-Genetic- Fuzzy Approach in Dynamic Environment.
International Journal of Information Technology, 21, 1-23.
7. Ni, J.J., Wu, W.B., Shen, J.R. and Fan, X.N. (2014) An Improved
VFF Approach for Robot Path Planning in Unknown and Dynamic
Environments. Mathematical Problems in Engineering, 2014, Article
ID: 461237. https://doi.org/10.1155/2014/461237
8. Asadi, B. (2015) Single and Dual Arm Manipulator Motion Planning
Library. IEEE/RSJ International Conference on Intelligent Robots and
Systems (IROS), Hamburg, 28 September-2 October 2015.
9. Vemula, A. (2017) Safe and Efficient Navigation in Dynamic
Environments. Carnegie Mellon University, Pittsburgh.
10. Spong, M.W., Hutchinson, S. and Vidyasagar, M. (2005) Robot
Modeling and Control. 1st Edition, Wiley, Hoboken.
11. Craig, J.J. (2005) Introduction to Robotics Mechanics and Control. 3rd
Edition, Pearson Prentice Hall Pearson Education, Inc., Upper Saddle
River, NJ.
12. Sadiq, A.T., Raheem, F.A. and Abbas, N.A.F. (2017) Optimal Trajectory
Planning of 2-DOF Robot Arm Using the Integration of PSO Based on
Interactive Heuristic D* Path Planning Solution Based on ... 139
8
Optimal Task Placement of a Serial
Robot Manipulator for Manipulability
and Mechanical Power Optimization
Uberlândia, Brazil
Faculty of Mathematics, Federal University of Uberlândia, Uberlândia,
2
Brazil
ABSTRACT
Power consumption and accuracy are main aspects to be taken into account
in the movement executed by high performance robots. The first aspect is
important from the economical point of view, while the second is requested to
satisfy technical specifications. Aiming at increasing the robot performance,
a strategy that maximizes the manipulator accuracy and minimizes the
mechanical power consumption is considered in this work. The end-effector
is constrained to follow a predefined path during the optimal task positioning.
The proposed strategy defines a relation between mechanical power and
manipulability as a key element of the manipulator analysis, establishing a
performance index for a rigid body transformation. This transformation is
used to compute the optimal task positioning through the optimization of a
multicriteria objective function. Numerical simulations regarding a serial
robot manipulator demonstrate the viability of the proposed methodology.
Keywords: Optimal task placement, Optimal robot path planning,
Multicriteria optimization
INTRODUCTION
Minimization of production costs and maximization of productivity are
some of the major objectives of industrial automation. In this scenario, serial
robot manipulators have been proved to be very useful tools.
Due to the augmenting use of serial robot manipulators to perform a
number of tasks in industry, requirements concerning higher precision,
improved productivity, reduced costs, and better manufacturing quality
become very important. To effectively explore all the potential of the robot,
its path planning is a subject of major concern.
There are several works in the literature dedicated to different approaches
concerning this subject. Regarding the off-line path planning approach, in
Saramago [1] a solution of moving a robot manipulator with minimum cost
along a specified geometric path in the presence of moving obstacles is
presented. The optimal traveling time and the minimum mechanical energy
of the actuators are considered together to form a multiobjective function to
be minimized along the process.
In Santos [2] a strategy to determine the trajectory for a defined movement
taking into account the requirements of torque, velocity, operation time and
robot positioning along the movement is proposed. The analysis is performed
as an off-line path planning through spline interpolation techniques in the
manipulator joint space. In the paper of Chiddarwar [3] an off-line path
planning for coordinated manipulation is proposed. The swept sphere
volume technique is used to model multiple robots and static obstacles. In
Optimal Task Placement of a Serial Robot Manipulator for ... 143
In the present paper both kinematics and dynamics features are considered
in the proposed task placement strategy. From the best of the authors’
knowledge, this problem was not previously addressed in the literature,
despite its importance and applicability.
The position of the end-effector is fixed according to a Cartesian
reference system, and described as a sequence of sucessive positions. The
optimization process applies a translational and rotational transformation
matrix with respect to a given reference frame. As a result, an optimal
sequence of kinematical robot positions is determined.
The outline of the paper is the following. In Section 2, a review about
manipulability is presented. In Section 3, the mechanical power concept is
discussed.
A strategy to describe the task and its positioning in the design space
is proposed in Section 4. Additionally in this section the design variables
are defined. In Section 5 the multicriteria programming formulation is
presented. In Section 6, the optimization strategy is established and the
numerical experiments are performed in Section 7. Finally, the conclusions
are drawn in Section 8.
MANIPULABILITY
Some factors should be taken into account when the posture of the
manipulator in the workspace is determined for performing a given task
during operation. An important factor is the ease of arbitrarily changing the
position and orientation of the end-effector at the tip of the manipulator.
As an approach for evaluating quantitatively the ability of manipulators
from the viewpoint of kinematics, the concepts of manipulability ellipsoid
and manipulability measure (Yoshikawa [23]) are presented.
Consider a manipulator with n degrees of freedom. The joint variables are
denoted by an n-dimensional vector, q = [q1, q2, ..., qn]T. An m-dimensional
vector r = [r1, r2, ..., rm]T, (m ≤ n) describes the position and orientation of
the end-effector.
The relation between the velocity vector v corresponding to r and the
joint velocity is given by:
(1)
where J(q) is the Jacobian matrix, computed from the kinematics
description of the robot.
146 Mechatronics and Robotics
The set of all end-effector velocities v that are realizable by joint velocities
such that the Euclidean norm of satisfies defines an ellipsoid in the
m-dimensional Euclidean space. In the direction of the major axis of the
ellipsoid, the end-effector can move at high speed, and in the direction of the
minor axis, the end-effector can move only at low speed. Also, the larger the
ellipsoid is, the faster the end-effector can move.
Since this ellipsoid represents an ability of manipulation, it is called the
manipulability ellipsoid.
The principal axes of the ellipsoid can be found by making use of the
singular-value decomposition of J(q).
One of the representative measures for the ability of manipulation
derived from the manipulability ellipsoid is the volume of the ellipsoid.
Since this volume is proportional to the eigenvalues of the Jacobian, it can
be seen as a representative measure.
Therefore, the manipulability measure, f1, for the manipulator
configuration, q, is defined as
(2)
Generally f1 ≥ 0 holds, and f1 = 0 if and only if rank J(q) < m.
There is a direct relation between singular configuration and
manipulability (through the Jacobian concept). According to this relationship,
it can be shown that the larger the manipulability measure, the larger the
ability of avoiding singular configurations.
MECHANICAL POWER
The consideration of the dynamics behavior of a serial robot manipulator is
of great importance for its path planning. This information allows a detailed
analysis and consequently the development of a precise control specification.
Many efficient schemes have been proposed to model the dynamics of
rigid multibody mechanical systems (Richard [24], Mata [25], Mata [26]).
The dynamics model can be obtained explicitly through algebraically
computation, or numerically through iterative computation.
The techniques based on the Newton-Euler method starts from the
dynamics of all individual parts of the system. They look at the instantaneous
or infinitesimal aspects of the motion, using vector quantities such as
Cartesian velocities and forces.
Optimal Task Placement of a Serial Robot Manipulator for ... 147
Alternatively, the Euler-Lagrange based methods starts from the kinetic and
potential energy of the entire system, by considering the states of the system
during a finite time interval. This approach works with scalar quantities,
namely the energies.
Independently of the approach used, at the end of the process the
generalized forces are determined:
(3)
where and are the joint position, velocity and acceleration,
respectively. M(q) is the joint space mass matrix and is the vector
of Coriolis and centrifugal forces. Vector g(q) is the vector of gravitational
forces and τ is the generalized force vector.
The energy that is necessary to move the robot is an important design
issue, because in real applications energy supply is limited and any energy
reduction leads to smaller operational costs. Due to the relationship that
exists between energy and force, the minimal energy can be estimated from
the generalized force, , that is associated to each joint i at the time instant
t0 ≥ t ≥ tf.
In the present contribution, for a manipulator with n degrees of freedom,
the interval between the initial time t0 and the final time tf is represented by
a set of N points and the computation is numerically carried out through
a recursive Newton-Euler formulation (Craig [27]). Then, the mechanical
power is used for design purposes as defined by:
(4)
This expression is representative of the phenomenon under study
because it considers both the kinematics and the dynamic aspects of the
trajectory, simultaneously (Saramago [1]).
TASK SPECIFICATION
The task is specified as a set S of N Cartesian points (x, y, z) and the
respective time step. Hence,
(5)
148 Mechatronics and Robotics
(6)
where β is a rotation angle in the z axis (given in rad) and δx , δy and
δz are translations along the reference axes (given in meters), respectively.
If desired, it is possible to extend this formulation to consider the rotation
with respect to the reference axes.
The physical meaning of considering a rotation angle with respect to x or
y is that the robot’s base plane will not be parallel to the ground’s reference
plane. In other words, the robot’s z axis (its height) is not normal to the
ground’s reference plane (xy plane). As positioning two planes according to
a precise angle may not be possible in a number of practical situations (even
when it is possible, it is not an easy procedure to be done, in general terms),
the rotation along the z axis only is considered, maintaining the robot’s base
plane parallel to the task’s base plane.
After defining the set S that describes the task, the inverse kinematics
with respect to each Cartesian point is computed.
The inverse kinematics computation can be achieved algebraically, by
considering the analytical model of the manipulator and its geometry, or
numerically. With the aim of presenting a general procedure, in the present
contribution the inverse kinematics problem is solved numerically through
Optimal Task Placement of a Serial Robot Manipulator for ... 149
(7)
where Pbase, Pend and are the robot base point, the end-effector
point, and the kinematics representation of the manipulator, respectively.
The penalty parameter used as a constraint in the following optimization
formulation is given by
(8)
and represents the sum of all end-effector positioning errors during the
movement.
Equation (7) means that the end-effector error is considered null when
the accuracy is better than the one specified by ε (given in meters). Otherwise,
this error corresponds to the distance from the end-effector specified position
to the position obtained from the direct kinematics calculation.
In some robotic implementations it is necessary to plan the layout of
the workspace, i.e., it is required to locate the robot base in such a way that
dexterity is maximized at (or around) given targets. It is worth mentioning
that such cases are also covered by the present formulation. Since the optimal
design ((β, δx, δy, δz) in Equation (6)) means a coordinate transformation
of the end-effector path regarding the fixed robot base, it is also possible
to adjust the location of the robot base regarding the fixed end-effector
reference through the relation (-β, -δx, -δy, -δz).
150 Mechatronics and Robotics
MULTICRITERIA PROGRAMMING
In problems with multiple criteria, one deals with a design variable vector x,
which satisfies all constraints and makes the scalar performance index that
is calculated from the m components of an objective function vector f(x)
as small as possible. This goal can be achieved by the vector optimization
problem:
(9)
where is the function domain (the design space).
An important feature of such multiple criteria optimization problem is
that the optimizer has to deal with objective conflicts (Eschenauer [29]).
Other authors discuss the same topic by defining compromise programming,
because there is no unique solution to the problem (Vanderplaats [30]).
(10)
where f1 and f2 are given by Equations (2) and (4), respectively. The
reference values are obtained from the computation of the
manipulability and mechanical power sum at the initial task positioning.
Therefore, in the ideal case the objective function value is lower than
2, where (the manipulability was increased) and (the
mechanical power was decreased).
The specification of a reference value to both the manipulability measure
and the mechanical power are of paramount importance in the construction
of the multicriteria objective function.
The inverse kinematics error is included in the problem as an equality
constraint in the optimization formulation, through the expression sE = 0
(Equation (8)).
It should be pointed out that the proposed formulation is able to handle
the case in which the task is outside the robot workspace. This configuration
reflects in the inverse kinematics error, which is taken into account as a
constraint (Equation (8)) of the optimization formulation.
Optimal Task Placement of a Serial Robot Manipulator for ... 151
(11)
(12)
for each joint j in the initial and final times t0 and tf, respectively. To
represent a rest-to-rest motion, zero initial and final velocities are considered
above.
OPTIMIZATION STRATEGY
An important feature of the proposed analysis is the computation of highly
nonlinear equations. Despite the good performance of classical nonlinear
programming techniques, the global optimal design is hardly reached due to
the existence of local minima in the design space.
With the aim of improving the robustness of the optimization process, the
use of a global methodology is considered. The tunneling algorithm (Levy
[31,32]) is a heuristic methodology designed to find the global minimum of
a function. It is composed of a sequence of cycles, each cycle consists of two
phases: a minimization phase having the purpose of lowering the current
function value, and a tunneling phase that is devoted to finding a new initial
point (other than the last minimum found) for the next minimization phase.
This algorithm was first introduced by Levy [31], and the name derives from
its graphic interpretation.
In summary, the computation evolves according to the following phases:
a) Minimization phase: Given an initial point q0, the optimization
procedure computes a local minimum q* of f(q). At the end, it is
considered that a local minimum is found.
b) Tunneling phase: Given q* found above, since it is a local
minimum, there exists at least one , such that
(13)
In other words, there exists q0 Z = {q Ω – {q*} | f(q) ≤ f(q*) }. To
move from q* to q along the tunneling phase, it is defined a new initial point
0
(14)
152 Mechatronics and Robotics
NUMERICAL SIMULATIONS
The numerical simulations performed to illustrate the presented methodology
use two robot models, namely the Scara and the Puma 560 robot manipulators.
These two types of serial robot manipulators are classical configurations in
the field. The case studies have been proposed due to their complexity as
compared to those normally found in the literature. The following results
considers only optimal profiles whose sum of the end-effector error is zero
(Equation (8)), given ε = 0.001 m (Equation (7)). This means that the inverse
kinematics solution is found for all positions with the required precision.
Different tolerance values may result from different optimal profiles.
As different values for the weighting factor k lead to different
optimization problems (Equation (10)), the value of the global minimum is
normalized according to the following equation
(15)
Optimal Task Placement of a Serial Robot Manipulator for ... 153
First Experiment
In this experiment, a subset of 7 reference points is defined: (0, 0.4, 0, –0.3),
(1, 0.4, –0.1, –0.3), (2, 0.2, –0.1, –0.3), (3, 0.2, –0.2, –0.3), (4, 0.4, –0.25,
–0.3), (5, 0.4, –0.2, –0.3) and (6, 0.2, –0.25, –0.3). Comparing these values
with Equation (5), the initial task is inside a rectangle defined by 0.2 ≤ x
≤ 0.4, –0.3 ≤ y ≤ –0.1 and z = –0.3. The corresponding traveling time is
defined by 0 ≤ t ≤ 6. These data are computed by interpolating cubic splines
for t, x, y and z, respectively. The interpolated data encompass a set of N
= 27 points (Equation (5)) that describes the path-following problem, as
presented in Figure 1.
At the original position, the values of the design variables are: β= 0 rad, δx =
0 m, δy = 0 m, and δz = 0 m. These values are used as an initial guess for the
optimization process (Equation (10)).
154 Mechatronics and Robotics
Second Experiment
In this experiment a subset of 8 reference points is defined: (0, 0.2, 0,
–0.3), (1, 0.4, –0.1, –0.3), (2, 0.4, –0.05, –0.3), (3, 0.2, –0.1, –0.3), (4,
0.25, –0.15, –0.3), (5, 0.3, –0.17, –0.3), (6, 0.27, –0.13, –0.3) and (7, 0.26,
–0.2, –0.3). By comparing these values with Equation (10), the initial task
is inside a rectangle defined by 0.2 ≤ x ≤ 0.4, –0.2 ≤ y ≤ 0 and z = –0.3.
The corresponding traveling time is defined by 0 ≤ t ≤ 7. These data are
computed by interpolating cubic splines for t, x, y and z, respectively. The
interpolated data encompass a set of N = 54 points
Figure 5. Proposed (---) and optimal (––) paths from two different perspectives.
(Equation (5)) that describes the path-following problem, as presented
in Figure 7” target=”_self”> Figure 7.
Optimal Task Placement of a Serial Robot Manipulator for ... 157
Third Experiment
In this experiment, a subset of 7 reference points is defined: (0, 0.4, 0.1, 0), (1,
0.4, 0, 0), (2, 0.4, 0, –0.4), (3, 0.4, –0.3, –0.4), (4, 0.4, –0.4, 0), (5, 0.4, –0.3,
0), and (6, 0.4, –0.4, –0.4). By comparing these values with Equation (5), the
initial task is inside a rectangle defined by x = 0.4, –0.4 ≤ y ≤ 0.1, and –0.4
160 Mechatronics and Robotics
Figure 11. Proposed (---) and optimal (––) paths from two different perspec-
tives.
Figure 12. Joint coordinates and the corresponding manipulability index (Scara
manipulator).
Optimal Task Placement of a Serial Robot Manipulator for ... 161
Figure 13. The proposed path-following problem from two different perspec-
tives.
At the original position, the values of the design variables are: β = 0 rad,
δx = 0 m, δy = 0 m, and δz = 0 m. These values are used as an initial guess for
the optimization process (Equation (10)).
The optimal values of both the manipulability and the mechanical power
indexes demonstrate that there is no significant difference among most of
the results presented. The best values are provided by k = 0.4 and 0.5.
Details about the minimal performance index, k = 0.5, are presented
in Table 3” target=”_self”> Table 3. The overall performance was improved
40%, while the minimal manipulability and the mechanical power were
improved 64% and 41%, respectively.
Some of the values obtained by using other weighting factors may
result in further improvement of either the manipulability or the mechanical
power, separately. However, the above configuration is a global minimum
of the performance index f4.
The performance provided by the tunneling process is shown in Figure
14.
After the 5th iteration of the Tunneling process, no further improvement
was achieved for the next five conescutive iterations. This behavior defines
the stop criteria for the global optimization search. In this case, it is possible
to see that good results are obtained in a few iterations.
162 Mechatronics and Robotics
Figure 15 shows the robot performing the path following task and Figure
16 presents the optimal positioning, according to the optimal design variable
values given in Table 3.
Figure 14. Performance index along the tunneling iterations (for k = 0.5).
Figure 16. Proposed (---) and optimal (––) paths from two different perspec-
tives.
Optimal Task Placement of a Serial Robot Manipulator for ... 163
Figure 17. Joint coordinates and the corresponding manipulability index (Puma
560 manipulator).
The joint coordinates at each time step are obtained by inverse kinematics
computation. As the task is considered as a set of consecutive Cartesian
points with continuous transition, the corresponding joint coordinates
presents a smooth behavior. Figure 17 shows the initial and optimal joint
coordinates for the first three links and the corresponding manipulability
index.
The previous example considered a task initially defined as being
parallel to the xz plane, i.e., the path is parallel to a vertical plane. It can be
noticed that the optimal design corresponds to rotating the task about 0.391
rad anticlockwise with respect to the reference frame.
CONCLUSIONS
In this work the path-following problem of robot manipulators was
addressed. An approach to increase the manipulability while decreasing
dynamics requirements was proposed by using optimization techniques.
The concept of manipulability was revisited. The basic idea behind
manipulability analysis consists in describing directions either in the task
space or the joint space that optimize the ratio between a measure of effort in
the joint space and a measure of performance in the task space. In this work
a classical formulation to describe the relation between the two spaces (task
and joint spaces) was considered.
Since the end-effector presents various performances due to a number of
factors (control specification, joint flexibility, joint friction, etc.), a general
methodology to reduce high torques and to increase manipulability was
proposed. The corresponding equation considers the overall performance
of the entire manipulator, as the determinant of the Jacobian is proportional
164 Mechatronics and Robotics
to the product of its singular values, and the mechanical power index is
obtained by a sum involving all the joints. The methodology also includes the
manipulator positioning error analysis, through the definition of constraint
functions in the optimization model.
The serial manipulator usually performs differently in different zones of
the workspace. Considering also the cases for which different joint velocity
and torque profiles may result different performance levels, a better position
into the workspace (in the sense that the performance indexes are optimized)
for performing a specified task was determined.The optimization problem
was based in a global search heuristic, improving the robustness of the process
with respect to the initial guess.Numerical experiments showed that the task
repositioning successfully increases the manipulability and reduces the total
mechanical power, simultaneously. As changes in the task positioning may
result different associated manipulability indexes and different joint velocity
and torque profiles, the methodology searches for an optimal way of using
the robot possibilities (as described by the performance indexes) by finding
an optimal position into the workspace to perform the task.
In some practical cases, the manipulated objects cannot be displaced
because of physical constraints such as weight or geometry. Even in such
cases the proposed strategy is suitable to be applied by changing the position
of the manipulator itself. The corresponding changes in the robot’s base are
given by the optimal profile by using the opposite sign of the optimal values,
i.e., (-β, -δx, -δy, -δz). It is also possible to adjust the analysis to consider
different objectives, as the total torque instead of the mechanical power, for
example. The steps to be followed remain the same.
In summary, the main contributions of this paper are the simultaneous
analysis of both the required mechanical power and the manipulability
(dynamics and kinematics aspects), together with a global deterministic
heuristic to achieve the optimal task placement.
As the proposed methodology is suitable for serial manipulators of
arbitrary complexity, the authors believe that the present methodology
represents a contribution concerning general optimal robot path planning
and design.
ACKNOWLEDGEMENTS
The authors are thankful to the agencies FAPEMIG and CNPq (INCT-EIE)
for partially funding this work.
Optimal Task Placement of a Serial Robot Manipulator for ... 165
REFERENCES
1. S. F. P. Saramago and V. Steffen Jr., “Optimization of the Trajectory
Planning of Robot Manipulators Taking Into-Account the Dynamics
of the System,” Mechanism and Machine Theory, Vol. 33, 1998, pp.
883-894.
2. R. R. D. Santos, V. Steffen Jr. and S. F. P. Saramago, “Robot Path
Planning: Avoiding Obstacles,” In: 18th International Congress of
Mechanical Engineering, Ouro Preto-MG, Brazil, 2005.
3. S. S. Chiddarwar and N. R. Babu, “Offline Decoupled Path Planning
Approach for Effective Coordination of Multiple Robots,” Robotica,
Vol. 28, No. 4, 2010, pp. 477- 491.
4. R. R. D. Santos, V. Steffen Jr. and S. F. P. Saramago, “Optimal Path
Planning and Task Adjustment for Cooperative Flexible Manipulators,”
ABCM Symposium Series in Mechatronics, Associação Brasileira de
Engenharia e Ciências Mecânicas, ABCM, Vol. 3, 2008, pp. 236-245.
5. O. Von Stryk and M. Schlemmer, “Optimal Control of the Industrial
Robot Manutec r3,” Computational Optimal Control, International
Series of Numerical Mathematics, Vol. 115, 1994, pp. 367-382.
6. R. R. D. Santos, V. Steffen Jr. and S. F. P. Saramago, “Robot Path Planning
in a Constrained Workspace by Using Optimal Control Techniques,”
In: III European Conference on Computational Mechanics, Lisbon,
Portugal, 2006, pp. 159-177.
7. J. E. Bobrow, S. Dubowsky and J. S. Gibson, “Time-Optimal Control
of Robotic Manipulators along Specified Paths,” The Intemational
Journal of Robotics Research, Vol. 4, 1995, pp. 3-17.
8. J. E. Slotine and H. S. Yang, “Improving the Efficiency of Time-Optimal
Path-Following,” IEEE Transaction on Robotics and Automation, Vol.
5, No. 1, 1989, pp. 118- 124.
9. [9] D. Constantinescu and E. A. Croft, “Smooth and TimeOptimal
Trajectory Planning for Industrial Manipulators along Specified Paths,”
Journal of Robotic Systems, Vol. 17, 2000, pp. 233-249.
10. H. Seraji, “Reachability Analysis for Base Placement in Mobile
Manipulator,” Journal of Robotic Systems, Vol. 12, 1995, pp. 29-43.
11. S. Zeghloul and J. A. Pamanes, “Multi-Criteria Optimal Placement of
Robots in Constrained Environments,” Robotica, Vol. 11, No. 2, 1993,
pp. 105-110.
12. J. T. Feddema, “Kinematically Optimal Robot Placement for Minimum
166 Mechatronics and Robotics
9
Kinematical Analysis and Simulation of
High-Speed Plate Carrying Manipulator
Based on Matlab
ABSTRACT
In order to construct the more effective kinematics method for industry, by
taking a high-speed plate handing robot as an example, the structure and
parameters of the robot linkages are analyzed, and the standard Denavit-
Hartenberg method is applied to establish the coordinates and the kinematic
equation of the linkages. Depending on the graphics and matrix calculation
INTRODUCTION
As a typical representative and a main technical means of information
technology and advanced manufacturing technology, the complete set of
automatic stamping processing line has become a high technology to which
the developed countries have paid much attention. And its development
level has become one of the most important standards to measure a nation’s
technical development. It has been extensively applied in the industry of
mechanical manufacturing, nuclear, aerospace, energy and transportation,
petroleum chemistry, building, electronic and etc. [1]. But, one key aspect of
the automatic stamping processing line is to develop a mechanical manipulator
with the characteristic of the high-speed and dynamic transmission.
At present, the sheet metal forming production line is operated at the
high productive rate of 15 - 25 SPM. Removing the time spent by the press
stamping operations, there is only about 2 s left to the robot manipulator. In
such a short period of time, the manipulator has to operate at a high speed
of 200 - 250 m/min, so that the operations such as loading, moving and
unloading can be fulfilled. And this has put forward higher requirements
to the mechanical structure, material friction characteristics, and structural
dynamic characteristics of the robot manipulator.
Taking a high-speed plate carrying manipulator as the research object,
this paper firstly analyzes the structure and connecting rod parameters; then
adopts the standard D-H method [2] to establish the kinematics equation;
and then discusses the positive and inverse kinematics algorithm, and the
trajectory planning problems; finally in the environment of Matlab, the
kinematics model is built to take kinematics simulation by using of Robotics
Toolbox [3]. In simulation process, we can not only directly observe the
Kinematical Analysis and Simulation of High-Speed Plate Carrying ... 171
robot motion, but also get the required data in the graphic form. Therefore,
the virtual performance of the product can be tested in the conceptual design
stage, so as to improve the design performance, reduce design cost and
decrease product development time.
D-H Transformation
In order to describe the translational and rotational relationships between
adjacent rods, Denavit and Hartenberg (1955) have proposed a matrix
method to establish the possessed coordinate system for each rod in the
linkage chain. This method is to establish a homogeneous transformation
matrix for each of the join bar coordinate, which represents the relationship
of the previous coordinate system of the bar, and the principle [4,5] is
detailed as follows:
OXYZ: A fixed reference coordinate system of the fixed coordinates, is
called the coordinate system.
OiXiYiZi: Fixed connected with the member bar of number I of the robot,
the origin of coordinate is at the center point of the joint of the I + 1th.
174 Mechatronics and Robotics
(1)
Now, putting the D-H parameters and joint variables of manipulator
substitution of the Table 1 in (1), we get the homogeneous transformation
matrix of the i coordinate system relative to the position and posture of the
base coordinate system, expressed as:
(2)
Especially, when , can be obtained, it determines the position
and posture of the end of the manipulator relative to the base coordinate
system, the matrix of T can be expressed as:
(3)
Note:
the variables of each joint of the manipulator, here in our case the variables
are and .
The target point of the movement of the manipulator center is
. By using the transitional joint, the manipulator’s center
point can coincide with the target point, and the target point can be expressed
as: , then we can establish the following equation:
(4)
Combining (3) and (4) and solving the kinematics equation, we can get
the following equation.
(5)
(6)
The error is very small through the comparative analysis of the actual
value and set point, which explains that the kinematics equation and the
model are reliable and consistent.
Figure 5. The three-dimensional diagram and slide block control chart of the
high-speed plate carrying manipulator.
each joint’s displacement curve, as shown in Figure 7(a) and (b). (Note:
q represents the displacement, i represents the joint Numbers). Also we can
draw speed curve as shown in Figures 7(c)-(d), acceleration curve as shown
in Figures 7(e)-(f), through the corresponding functions:
CONCLUSIONS
Ke Wang, Jiping Zhou The paper carries out a kinematical analysis of a
highspeed plate carrying manipulator which is designed by using the module
function of Robotics Toolbox in the Matlab. Based on this, we have achieved
the following conclusions in three aspects: 1) the structure and parameters
of the plate carrying manipulator have been designed; 2) the kinematics
model has been established with the standard D-H method, and its positive
and inverse kinematics have been analyzed; 3) by using Robotics Toolbox
in Matlab, the kinematics model has been set up in order to make the study
of the motions of manipulator easier and more apparent, and to conduct a
kinematics simulation for the manipulator for obtaining the displacement,
velocity and acceleration curve of each manipulator joint. Thus, this proves
the validity of the design of manipulator structure parameters, as well
as provides the reliable basis for the theoretical analysis in terms of the
design, dynamics analysis, and trajectory planning of control system of the
manipulator.
182 Mechatronics and Robotics
REFERENCES
1. Q. Wang, X. J. Hu and L. Z. Li, “The Kinematic Analysis and Simulation
of the Humanoid Welding Manipulator Based on the MATLAB,”
Materials Science and Engineering College of Hefei University of
Technology, 2011.
2. J. J. Craig, “Introduction to Robotics,” Mechanical Industry Press,
Beijing, 2006.
3. P. I. Corke, “A Robotics Toolbox for MATLAB,” IEEE
Robotics and Automation Magazine, Vol. 3, No. 1, 1996, pp. 24-
32. doi:10.1109/100.486658
4. X. S. Jiang, “Introduction to Robotics,” Liaoning Science and
Technology Press, Shenyang, 1994.
5. Z. X. Cai, “The Robotics,” Tsinghua University Press, Beijing, 2000.
6. J. X. Luo and G. Q. Hu, “The Kinematic Simulation Research of Robot
Based on the MATLAB,” Journal of Xiamen University, JCR Science
Edition, Vol. 44, No. 5, 2005, pp. 640-644.
7. J. Han and L. Hao, “Trajectory Planning and Simulation of Robot in
Joint Coordinate System,” Journal of Nanjing University of Science
and Technology, Vol. 24, No. 6, 2000, pp. 540-543.
8. Z. X. Wang and W. X. Fan, “The Kinematics Analysis and Simulation
of Industrial Robot Based on the Matlab,” Mechanical and Electrical
Engineering, Vol. 29, No. 1, 2012, pp. 33-37.
9. The MathWorks Inc., “Matlab the Language of Technical Computing
Version 6,” The MathWorks Inc., 2002.
10. P. I. Corke, “Robotics Toolbox for Matlab (Release 7.1) [EB/OL],”
2002. http://www.cat.csiro.au/cmst/staff/pic/robot
CHAPTER
10
Flexible Impact Force
Sensor
ABSTRACT
When human body is punched by boxing glove, both body surface and
boxing glove deform in a complex shape. The purpose of this study is to
develop a flexible sensor that can be used in such interface. Firstly, several
mechanical phenomena, which are the cause of the error signal of the sensor,
are discussed. These are the influences of out-of-plane bending deformation,
Citation:- Liu, C. , Fujimoto, Y. and Tanaka, Y. (2014), “Flexible Impact Force Sen-
sor”. Journal of Sensor Technology, 4, 66-80. doi: 10.4236/jst.2014.42008.
Copyright: © 2011 by authors and Scientific Research Publishing Inc. This work is li-
censed under the Creative Commons Attribution International License (CC BY). http://
creativecommons.org/licenses/by/4.0
184 Mechatronics and Robotics
shear force caused by rubbing, shear force caused by the Poisson’s effect
of contact material, and the lateral compressive force caused by the
overhanging deformation of flexible material. Then, a prototype sensor
that can eliminate the error factors of these is developed. The sensor is
a distribution type impact sensor in which sixteen sensor elements are
arranged in a 4 × 4 matrix. Punching experiments using a boxing glove are
carried out by installing the sensor on the load cell, on the concrete wall and
on the sandbag. From the experiment, it is found that the impact force can
be measured with good accuracy by using the sensor. Despite the fact that
the sensor has inadequate distribution number of sensor elements, the sensor
structure meets mechanical requirements for the flexible impact sensor.
Keywords: Flexible Sensor, Impact Force, Boxing, Sandbag,
Distribution Type Sensor, Human Body
INTRODUCTION
When a boxer punches the human body, both body surface and boxing
glove deform in a complex shape. At present, the accurate measurement of
impact force in such interface is said to be difficult. The reason is because
the complex deformation and the complex force distribution are induced in
the flexible interface and they affect the sensor output. Up to now, many
kinds of tactile sensors such as resistance type [1], piezo-resistance type [2]
, capacitance type [3] and light intensity type [4] have been developed. Most
of them are assuming to mount the sensor on the robot members. However,
it has not been studied much about the influences of the deformation and
the associated force distribution at the flexible interface on the accuracy
of the sensor output. The reason of this is most of the tactile sensors are
intended to use under a small force level and a small deformation level on
the installation surface [5] [6].
By the way, it has been reported that peak impact force of about 4 kN
is generated when a strong boxer punches the face of crash test dummy [7].
When the hitting target is sandbag, the interface of glove and sandbag will
be deformed more greatly. A suitable sensor which can measure the impact
force in such condition is not seen at present. The purpose of this study
is to provide a sensor that can measure the impact force accurately at the
interface of large deformation and large force.
In this study, firstly, difficulty of force measurement in flexible interface
is discussed based on the experiments using a sheet type sensor (Chapter
2). Then, some mechanical factors necessary to achieve the flexible sensor
Flexible Impact Force Sensor 185
are discussed (Chapter 3). These factors are the out-of-plane bending
deformation of sensor, the shear force caused by rubbing on sensor surface,
the shear force generated by contact friction, and the transverse compression
caused by the overhanging deformation of the flexible material.
In Chapter 4, as a sensor that can overcome the mechanical problems in
flexible interface, a flexible sensor of distribution type is developed. In Chapter
5, boxing punching tests are carried out using a prototype sensor by installing
the sensor on the sponge rubber sheet which is placed on the load cell or the
concrete wall. Further, the punching test is carried out by installing the sensor on
a sandbag. In particular, the present study is focused to achieve a sensor structure
that is hard to cause trouble under a big impact force like boxing punching.
fixed on the top of a load cell. Sponge rubber plate with 40 mm thickness was
placed on the acryl plate. The sensor was put on it and pressed by spherical
body with 50, 100 or 150 mm diameter. The experiment was carried out by
changing the compressive force.
Figure 4 compares the relationship between the peak force of the load
cell and the sensor output. A good linear relationship is seen in the figure.
Despite the change in the diameter of spherical body, the output gradient does
not change much. The dashed line in the figure shows the output gradient
under the condition that the sensor is placed directly on the flat acryl plate.
The sensor detects the normal force on the curved surface, but the load cell
detects the vertical component. The difference in gradient of the two straight
lines is considered to be related with the detection characteristics.
Figure 2. Sheet type flexible sensor. (a) Main body of sensor with slits in outer
edges; (b) Sensor covered by stretch cloth.
Figure 4. Sensor output vs. load cell output. Dashed line is the output gradient
of the sensor on flat plane surface..
188 Mechatronics and Robotics
Figure 5. Fist press, palm press and fist punch of the sensor placedon the sponge
rubber plate.
Figure 6. Relationship between the output peak of sensor and the peak force of
load cell.
In the experiment of Section 2.1, the loading side was a hard material
with simple contact surface and the receiving side was a flexible material.
On the other hand, the contact surface of fist or palm is more complex and
it changes in every trial. Further, rubbing force and inclined force are also
working. The cause of the scatter of the output gradient can be considered as
the influences of complex force distribution and the excessive out-of-plane
deformation occurring in the sensor. This result means that a single sheet-
type sensor is not suitable to achieve the measurement of impact force on the
flexible material surface.
Figure 7 classifies the condition of impact force measurement into
four levels, in view of the difficulties of measurement. Figure 7(a) is the
condition that a hard object with simple contact geometry collides. Figure
7(b) is the condition that a soft object with simple contact geometry
collides. In accordance with the deformation of soft object, in this condition,
unexpected shear forces (rubbing force and contact friction) work at the
Flexible Impact Force Sensor 189
interface. Figure 7(c) is the condition that a hard object with complex contact
geometry collides. In this case, the out of-plane bending deformation of the
sensor and the distributions of shear forces at the interface become more
complex. Figure 7(d) is the condition that a soft object with complex contact
geometry collides. In this case, stress distribution and deformation becomes
most complicated. Boxing punching on sandbag surface is close to the
condition of Figure 7(d).
Figure 9. Example of the tensile strain (red arrow) generated in the sensor ele-
ment caused by rubbing force. F is compressive force, S is rubbing force and
RF and RS are the corresponding reaction forces.
Figure 11. Variability of sensor signal caused by the difference of contact mate-
rial.
1 mm thickness. Sponge rubber sheet (T = 5 mm, 10 mm thick) or natural
rubber sheet (T = 1 mm) or cork plate (T = 2 mm) was sandwiched between
the sensor element and the aluminum plate. As the contact material becomes
softer, the sensor output increases as shown in the figure. In the case of thick
sponge rubber (T = 10 mm), the magnitude becomes several times compared
with the case of the aluminum plate.
192 Mechatronics and Robotics
Figure 13. Impact force waveform in which vertical and lateral compression is
mixed.
Figure 16. Relationship between sensor output and load cell output obtained by
the hammer punching test.
Figure 17. Prototype sensor. (a) Flexible sensor consisting of 4 × 4 sensor ele-
ments; (b) Conductive cloth bag; (c) Hold the sensor bend.
Even though the sensor element itself is a rigid body, the connecting
portions of the sensor elements can keep the flexibility.
Figure 18 shows the cross section view of the flexible sensor. The
upper aluminum plates of all sensor elements are in contact with the upper
conductive cloth. Further, the lower aluminum plates of all sensor elements
are in contact with the lower conductive cloth. In order to prevent short-
circuiting of the upper and lower conductive cloth, bilayer deer skin is
sandwiched in the thickness center of the sensor as the insulation sheet.
After sewing the intermediate portion of the sensor element and element,
each sensor element is in a state put into a pocket.
The electrical wires are connected to the lower and upper conductive
cloth respectively. In other words, the output signal of the sensor element
is obtained from the conductive cloth via upper and lower aluminum plate.
Because there is no electrical wiring to connect each sensor element, there
is no risk of electric disconnection.
196 Mechatronics and Robotics
It is obvious that the flexibility of the sensor increases as the increase of the
number of sensor elements distributed. If the sensor element is assumed as
“a plate placed on the elastic floor”, the bending moment reduces along with
the sensor dimension becomes smaller. Therefore, it is possible to reduce the
sensor thickness as the size of sensor element is small.
EXPERIMENT
Figure 19. Impact test of flexible sensor on the load cell. (a) Sponge rubber (t =
40 mm) put on the load cell; (b) Sensor placed on the sponge rubber; (c) Punch-
ing test by boxing glove.
Flexible Impact Force Sensor 197
Also, similar experiment was conducted under the condition that the flexible
sensor was placed directly on the acryl plate.
Figure 20 shows an example of impact force waveform measured in the
above punching test. Solid and dashed curves indicate the experiments on
the sponge rubber and the acryl plate respectively. The red curve is sensor
output and the black curve is load cell output. The impact force waveform
of the sensor is in good agreement with that of load cell.
Figure 21 shows the relationship between output peak of the sensor and
output peak of the load cell. Regardless of the presence or absence of sponge
rubber, the sensor output shows almost same gradient. This means that the
sensor can measure the impact force in the flexible interface.
Figure 22 shows sponge sheet with 40 mm thickness overlaid on the
sponge rubber. The total thickness of the bilayer is 80 mm. The sensor was
placed on the sponge sheet and punched by boxing glove. In the figure, two
examples of the impact force waveform are compared with that of the load
cell output. As the sponge is so soft, the sensor was greatly deformed in a
concave shape. There is somewhat difference in the waveform of sensor and
load cell.
Figure 23 shows the sensor placed in bent condition on the wood
semicircular column with diameter 150 mm. Some positions of the sensor
were fixed with gum tape. Then, the sensor was punched by boxing glove.
The relationship between the load cell output and the sensor output are
indicated by black diamond. The red circle in the figure is the output gradient
of the experiment shown in Figure 22. These two gradients are substantially
coincident. Moreover, this output gradient is in good agreement with that
in Figure 21.
Figure 20. Impact force waveforms of the sensor and the load cell. Red curve
is sensor, black curve is load cell.
Figure 21. Comparison of the peak output between sensor and load cell.
Figure 22. Impact force waveform of sensor and load cell. Sponge is over-
laid on top of the sponge rubber.
Flexible Impact Force Sensor 199
Figure 23. Relationship between sensor output and load cell output. Black dia-
mond: semicircular column. Red circle: sponge & sponge rubber sheet.
Figure 24. Boxing punching of the sensor installed on the 40 mm thick sponge
rubber pasted on concrete wall.
Four impact force waveforms weakly punched by student “U” and two
waveforms strongly punched by student “F” are shown in the figure. The
peak impact force by the strong punching is around 2 kN. The contact time
is about 30 - 40 msec.
Therefore, the forward speed of the glove persists for a while even after
passing the peak load. Since the pull-back of the glove is delayed, the contact
between glove and sandbag is continued for a while. The difference of the
gradient of unloading segment is clear when compared to the waveform of
wall punching in Figure 24.
Figure 26. Impact force waveform of the sensor installed on the soft sandbag
surface and punched by boxing glove.
Flexible Impact Force Sensor 201
Figure 28(a) shows the sensor fixed to the instep of soccer shoes with cords
and rubber bands. Figure 28(b) shows the impact force waveforms of instep
kick of soccer ball. The strength of the kick is the extent to path to nearby
players. The peak impact force is about 1 kN and the impact time is about 12
- 15 msec. This contact time is close to the value shown in the reference [11].
From the experience, it seems possible to miniaturize the element size to
15 mm square and 5 mm thickness. When the element becomes smaller than
that, however, reconsideration of the manufacturing process is necessary. In
spite of the rough distribution density of the sensor elements, the sandbag
sensor was not much discomfort at the time of use. On the other hand, in
the soccer shoe sensor, the sensor weight and the lack of flexibility were
felt uncomfortable. “Flexibility” and “light weight” are thought to be more
important for the sensors installed on human body side.
CONCLUSIONS
This study is aimed to develop a flexible sensor which can be used in the
measurement of impact force in flexible interface. The boxing punching is
chosen as the measurement object. The results obtained are as follows.
Figure 27. (a) Sensor elements of 25 mm square and 8 mm thickness; (b) Flex-
ible sensor made fromsix sensor elements.
Figure 28. (a) Set up of the sensor on soccer shoe; (b) Impact force wave-
forms of instep kick of soccer ball.
202 Mechatronics and Robotics
ACKNOWLEDGEMENTS
Chen Liu, Yukio Fujimoto, Yoshikazu Tanaka The authors wish to thank
Tadashi Uesugi, Takahiro Namikawa, Atsushi Maeda and Takuya Okada
of the Graduate School of Engineering, Hiroshima University for their
assistance in experiments. This work was partially supported by the Venture
Business Laboratory of Hiroshima University.
Flexible Impact Force Sensor 203
REFERENCES
1. Ishikawa, M. and Shimojo, M. (1982) A Method for Measuring the
Center Position of a Two Dimensional Distributed Load Using
Pressure-Conductive Rubber. Transactions of the Society of Instrument
and Control Engineers, 18, 730-735.
2. Mei, T., Li, W.J., Ge, Y., Chen, Y., Ni, L. and Chan, M.H. (2000) An
Integrated MEMS Three-Dimensional Tactile Sensor with Large Force
Range. Sensors and Actuators, 80, 155-162. http://dx.doi.org/10.1016/
S0924-4247(99)00261-7
3. Lee, H.K., Chang, S.I. and Yoon, E. (2006) A Flexible Polymer
Tactile Sensor: Fabrication and Modular Expandability for Large Area
Deployment. Journal of Micro Electromechanical Systems, 15, 1681-
1686. http://dx.doi.org/10.1109/JMEMS.2006.886021
4. Donati, M., Vitiello, N., De Rossi, S.M.M., Lenzi, T., Crea, S.,
Persichetti, A., Giovacchini, F., Koopman, B., Porobnik, J., Munih, M.
and Carrozza, M.C. (2013) A Flexible Technology for the Distributed
Measurement of Interaction Pressure. Sensors, 13, 1021-1045. http://
dx.doi.org/10.3390/s130101021
5. Lee, M.H. and Nicholls, H.R. (1999) Review Article Tactile sensing for
Mechatronics—A State of the Art Survey. Mechatronics, 9, 1-31. http://
dx.doi.org/10.1016/S0957-4158(98)00045-2
6. Lumelsky, V.J., Shur, M.S. and Wagner, S. (2001) Sensitive Skin. IEEE
Sensors Journal, 1, 41-51. http://dx.doi.org/10.1109/JSEN.2001.923586
7. Walilko, T.J., Viano, D.C. and Bir, C.A. (2005) Biomechanics of the
Head for Olympic Boxer Punches to the Face. British Journal of Sports
Medicine, 39, 710-719. http://dx.doi.org/10.1136/bjsm.2004.014126
8. Fujimoto, Y. and Setyanto, A.T. (2007) Sheet Type Impact Force Sensor
by the Use of Piezoelectric Film. Transactions of the Japan Society
of Mechanical Engineers, Series C, 73, 184-191. (in Japanese) http://
dx.doi.org/10.1299/kikaic.73.184
9. Fujimoto, Y., Shintaku, E., Tanaka, Y. and Fujiyoshi, J. (2012) Pad
Type Compressive Force Sensor Suitable for High-Speed Impact
Force Measurement. Transactions of the Japan Society of Mechanical
Engineers, Series C, 78, 2438-2449. (in Japanese)
10. Fujimoto, Y., Liu, C., Uesugi, T., Tanaka, Y. and Shintaku, E. (2013)
Pipe Surface Sensor for Impulsive Force Measurement. Transactions
of the Japan Society of Mechanical Engineers, Series C, 79, 113-123.
204 Mechatronics and Robotics
11
A Sensing and Robot Navigation of
Hybrid Sensor Network
China
Department of Electrical Computer Engineering, Michigan Technological
2
ABSTRACT
Traditional sensor network and robot navigation are based on the map of
detecting fields available in advance. The optimal algorithms are explored to
solve the energy saving, shortest path problems, etc. However, in practical
Citation:- S. Yao, J. Tan and H. Pan, “A Sensing and Robot Navigation of Hybrid Sen-
sor Network,” Wireless Sensor Network, Vol. 2, No. 4, 2010, pp. 267-273. doi: 10.4236/
wsn.2010.24037.
Copyright: © 2010 by authors and Scientific Research Publishing Inc. This work is li-
censed under the Creative Commons Attribution International License (CC BY). http://
creativecommons.org/licenses/by/4.0
206 Mechatronics and Robotics
environment, there are many fields, whose map is difficult to get, and need to
detect. This paper explores a kind of ad-hoc navigation algorithm based on
the hybrid sensor network without the prior map. The system of navigation
is composed of static nodes and mobile nodes. The static nodes monitor
events occurring and broadcast. In the system, a kind of cluster broadcast
method is adopted to determine the robot localization. The mobile nodes
detect the adversary or dangerous fields and broadcast warning message.
Robot gets the message and follows ad-hoc routine to arrive the events
occurring place. In the whole process, energy saving has taken into account.
The algorithms of nodes and robot are given in this paper. The simu- late and
practical results are available as well.
Keywords: Robot Navigation, Hybrid Sensor Network, Routine Search
INTRODUCTION
Sensor network is an advance technology in recent years. There have been
many papers in this field. Most of these are concerned in static sensor
network or mobile sensor network separately. The communication and
related technology have been studied. The static sensor network consists of
the fixed nodes and adopts the ad-hoc communication. This kind of system
can monitor the environment and send the messages to the base station.
Energy would be saved, but the system can hardly do any- thing better if
the whole condition is changed or the adversary available. In mobile sensor
network, the nodes can adjust their places according the algorithms, which
are based on the changing condition. However, at the same time, the energy
would be cost much. It is well know that is a serious problem for sensor
networks.
In navigation, most of these two systems always adopt a method that
is based on mapping. That is to say, the whole environment map is known
in advance. In spite that the optimal algorithms have been adopted in these
systems, it is not so suitable for the genuine conditions yet.
In this paper, we explore a hybrid sensor network. It is composed of both
static sensor nodes and mobile sensor nodes. The static sensor nodes are
deployed regularly, and detect the events of environment. While, the mobile
sensor nodes detect the adversary or dangerous field. When these fields are
found the mobile nodes will stop at the edge of adversary or dangerous fields
and send a message to robot. Finally, there will be more nodes at the edge
of adversary or dangerous fields to protect the robot not to approach. The
method is a compromise between the energy saving and flexible of whole
A Sensing and Robot Navigation of Hybrid Sensor ... 207
RELATED WORK
In a static sensor network, the quantity of nodes is an important problem.
The more nodes means the more coverage of whole area and more accurate
of navigation. Lots of research has focused on the sensor localization [1],
structural monitoring [2] or protocols [3]. Recently, researchers have begun
to explore the method, which low down the quantity. And, it is taken into
account that the adversary and dangerous fields available in environment.
Since these methods adopt the skeleton graphs and optimal algorithm,
the quantity of nodes and the price of communication have been slowed
down.
The quantity of nodes is not such a pivotal problem in mobile sensor
network because its flexibility. The problem of energy costs is tried to low
down by sorts of algorithm based on the passive algorithm [4, 5] or statistic-
learning algorithm. On the other hand, these algorithms are complicated.
They are not so easy to come into realities and the use more capacity of
nodes memory.
The sensor network of this paper develops the hybrid sensor network
and adopts a kind of cluster navigation. At first, the static sensor nodes are
deployed as regular as possible in the whole environment. That is to say, the
localization of static sensor nodes should be arranged at the correct place,
but we need not care it so much. These static sensor nodes may local in the
adversary or dangerous fields, but it does not matter. We can take measure to
solve the problem. The sensing program can be downloaded in these sensor
nodes. They can sense the change of environment area that we care. The
mobile sensor nodes, actually they are the small mobile robots, but they just
do the dangerous fields detecting work and nothing else. Navigated robots
are not these nodes. Actually, there is another robot being navigated. In
Figure 1, the static sensor nodes are deployed in the whole area. They have
sensing and communication abilities. Actually, these nodes are tmote Sky
[6] made by Moteiv Corporation. They are ultra-low power wireless module
for use in sensor networks, monitoring applications, and rapid application
prototyping. Tmote Sky leverages industry standards like USB and IEEE
802.15.4 to interoperate seamlessly with other devices.
208 Mechatronics and Robotics
The mobile sensor nodes can gather around the edge of adversary or
dangerous. The robot is in a waiting status. The mobile nodes are four-wheel
cart, which carry the Mica nodes. The Mica wireless platform serves as a
foundation for the emerging possibilities. With sensing, communication,
and I/O capabilities, Mica can simultaneously act as a data router, sensor
interface, and control point. Mica’s flexible design serves as a building
block for creating efficient application-specific protocols. And Mica does
not require use of predefined protocols [7].
We have done many works both on the communication of static sensors
and moving of mobile sensor nodes. This is the base that we can explore the
scheme.
points of the detecting area. They are the working nodes. That does not mean
we do not need large number of static nodes. If the quantity of static nodes
is not enough, the navigation can become very difficult, we discuss this
problem at navigation. The nodes, which are not on the key points, fall in
sleep. That means these sleeping nodes are working at the low power state.
When the system needs them, they will be taken up. These sleeping nodes
are in working state, and afford robot navigation. That is to say, there is large
number of nodes in system, but, the working nodes are not so many and
energy is saved sequentially. On the other hand, the flexibility and accuracy
of system can be maintained. Those nodes, which are not in adversary fields
and not on the point, sleep unless the robot navigation needs them (We will
discuss in Section 5.).
Except these non-key point nodes, there are still some nodes falls into
sleep. They are the nodes lie in the adversary or dangerous fields. They
can keep the sleeping status through the all process. Actually, they have
quitted events monitor. Therefore, the street map and skeleton graph can be
described as Figure 2.
Given such preliminary works, we modify the static nodes algorithm
to arrange the nodes places, monitor events and broadcast the message as
shown in Algorithm 1.
Algorithm 1. Deploy of static sensor nodes
1: while TRUE do
2: detect the location
3: if location! = dangerous fields & location! = key point then
4: send navigation requirement
5: else
6: node keeps sleep & quit monitor
7: end if
8: if navigation response! = TRUE then
9: node keeps sleep
10: else
11: broadcast navigation message
12: end if
13: end while
210 Mechatronics and Robotics
Figure 2. The diagram of on point nodes and the sleeping nodes. These nodes
which is in the dangerous fields are put into sleeping status.
Once there is an event occurs, the monitoring message can be sent by ad-
hoc module to find the robot. In order to save the energy, we set a threshold
for the event. Therefore, the working static nodes are not always at “working
status”.
(1)
Figure 3. The mobile nodes (a) Raster fashion moving; (b) Mobile nodes loca-
tion.
(2)
where l is the path length, ~ means the order of polynomial, and d is the
dimensions of function. And the time can be estimated as:
(3)
After the algorithms of static sensor nodes and mobile nodes are
designed, we can go into the next step. That is the navigation of robot.
NAVIGATION OF ROBOT
In navigation, the robot is Pioneer Ⅱ-AT four-wheel drive. The robot
also has a variety of expansion power and I/O ports for attachment and close
integration of additional sensors and other accessories. Expansion includes
an addressable I/O bus for up to 16 devices, two RS-232 serial ports, eight
digital I/O ports, five A/D ports, PSU controllers and more—all accessible
through a common application interface to the robot server software, P2OS
[13]. Therefore, we use the Mica as a communicate de- vice and serial ports
send dictate to the robot.
Most navigation is based on the detect fields map, which is available in
advance. So, the location of robot can be decided by coordinate of the map
A Sensing and Robot Navigation of Hybrid Sensor ... 213
[14]. And the optimal navigation algorithm can be explored. There have
been many algorithms in these fields. But they can hardly adapt the non-
mapping condition. That is to say, if there is a novel environment or a region
that its condition is difficult to get to know, what these algorithms do is
detect the whole region and acquire the whole information of the region.
Only after the information of region is known, can the optimal algorithm
can be adopted. Therefore, the energy may be cost in acquiring information.
In order to solve the problem, we explore a kind of non-mapping navigation
algorithm to guide the robot to arrive the event occurring place.
Generally speaking, navigation can be described as 3 questions. That
is where I am, where I should go, and how to arrive [15]. We adopt a kind
of dynamic cluster navigation algorithm. It can decide the location of robot
and use both static nodes and mobile nodes to guide the robot to arrive the
event place.
1). Where I am (Decide the location of robot)
Since what we discuss is no prior field map in this paper; the first thing
is the location of robot should be decided. It is the important step.
In order to acquire the location of robot, we adopt a cluster-navigating
algorithm. It depends on the ad-hoc relay communication. The static nodes
work as event heads is 1, and so do as other nodes in the next cluster until
the robot receive the message [16]. Robot receives the hierarchy of the static
nodes. Therefore, the robot can decide approximately how long the distance
between event and itself not precise. It is enough. And the robot should
distinct the message of event or relay. It can decide itself whether arriving
the place or not. At the same time, robot should also distinct the static or
mobile nodes, so it can decide it is close to danger fields or not.
If the static nodes are too sparse to decide the robot in a field, it is difficult
for robot to make a decision its direction to move. The robot cannot decide
where to go. The system must take measure to make up. Those nodes, which
are not on the key point, should be taken up the sleeping nodes around them
and form a new cluster in order to inform the robot.
Figure 4 shows the process of detect the robot location. The detail
algorithm of finding robot location is showed as follow:
Algorithm 3 The Transmission of Event 1: while robot! = arrived do
2: select working static nodes as cluster
3: if the nodes is too sparse then
214 Mechatronics and Robotics
Finally, the robot will arrive the place where events occur when it receive
the event’s source message.
The distance from robot to dangerous fields can be calculated by
potential [17]:
(4)
where R is the Euclidean distance from the edge to the mobile nods. In
this paper, it can be shown as:
(5)
(xe, ye) is the location of dangerous fields. The energy of mobile nodes
can be defined as follow:
(6)
Therefore, we can conclude the algorithm of robot navigation. There are
still some handicaps (they are the points not area) in the robot routine, but
there is sonar in the robot. The robot can avoid these handicap points. It is
the algorithm 4 shown as follow:
Algorithm 4 Navigation of Robot
1: While event message received = TURE Do
2: if arrived! = TURE then
3: if mobile warning message! = TURE then
4: use sonar (or communicate with the static nodes) follow navigation
routine
5: else
6: adopt the smooth routine 7: end if
8: else
9: send arrive message
10: inform the event node turn off requirement 11: end if
12: end while
Since there is no map in advance, what the robot does is receive the
message and do probing motion. Thus, the time to arrive the event is not
shortest and navigation routine is not the optimal routine. It is the cost of
216 Mechatronics and Robotics
none prior map available and the lower cost of both static and mobile nodes.
On the other hand, there are no so many working nodes in this system, the
energy cost of nodes can be saved.
4). How to wake up the sleeping nodes
When the static nodes, which are on the key point, is too sparse to
accomplish the navigation task, it is necessary to wake up the sleeping nodes.
At first, the cluster head, the static nodes on the key point, acquires the
message and sends the message to next closed cluster head. If the closed
cluster head can receive the message, it will send back a message to the
sending cluster head. That means that the statistic nodes are not too sparse
there. The message can relay to the next static nodes. If the closed cluster
head cannot receive the message, it cannot send back the message to the
sending cluster head. It means that the statistic nodes are too sparse there.
That is to say, the sleeping nodes in the cluster have to be taken up in order
to relay the message
The construction of the information deliver is shown in Figure 5
The static nodes in this system adopt tmote Sky [6] made by Moteiv
Corporation. They are ultra-low power wireless module for use in sensor
networks, monitoring applications, and rapid application prototyping. In an
open field, the communication distance of tmote nodes can reach 100 meters.
Consequently, message can be relayed on. If the sending cluster head does
not receive the backing message of the next cluster head, the sending cluster
head will send a signal to the sleeping nodes in its cluster, in order to wake
up these sleeping nodes. That is to say, the low power costs’ nodes (sleeping
nodes) will transfer from the low power costs’ state (sleeping) to working
state (waken up). Once the sleeping nodes are waken up, they can afford the
message delay events.
After the works have been accomplished, we simulate the sensing and
navigation of hybrid sensor network by using. Figure 6 and Figure 7 is result
of simulating experiments.
It can be seen that the energy saving of mobile work is better than the
static network. While, the energy saving of hybrid network is much better
than the mono-static and mono-mobile network.
The practical experiments are performed with the setup shown in Figure
8. A sequence process of navigation is shown from (a) to (d). An emergent
event occurs by a node (the temperature is over the threshold), the robot
(Pioneer Ⅱ-AT) should find the node and arrive the place in time. Part (a)
shows the initial experiments situation. There are 100 tmote nodes (static
nodes) and 8 small robots with the tmote nodes on them (mobile nodes) in
a 25m2 (5m×5m) filed. The lawn simulates the dangerous fields, and the
scope can be detected by moisture sensor on tmote nodes.
Despite the difficulties with the dangerous fields’ definition, the navigation
is successfully performed many times. Although arriving the final location
of the robot is not so quickly (about 3 minutes), the robot avoids the danger
fields and being guided by this hybrid sensor network system indeed. We
plan to conduct further re- search to get higher speed of this navigation
system.
CONCLUSIONS
This implementation demonstrated robot guidance by a hybrid sensor
network. The energy saving in sensor net- work is an important problem.
We take it as the primary problem in this system. The skeleton of adversary
or dangerous fields can be outlined by using the mobile sensor nodes and the
robot navigation comes into reali- ties on the ad-hoc method. In spite of the
less optimal routine, it is an effective scheme when there is no map of the
whole detecting fields in advance.
However, the precision of the navigation system greatly depends on the
robot hardware. Navigation can be problematic if environmental factors are
not so ideal, such as electrical cables, higher power electrical equipment and
steel structures exist. Although it is possible to compensate for these effects
in a known environment, a search situation may not permit elaborate offline
calibration. In addition, the response speed of whole system should be raised
higher. We are going to do further re- search on them.
ACKNOWLEDGMENT
We would like to give special thanks to Dr. Jindong Tan, the professor of
Michigan Technological University in
U.S.A. for his providing of experimental equipment and environment.
The same thanks go to Qi An, who helped
me with studies of the project and sensor network,
TinyOS.
220 Mechatronics and Robotics
REFERENCES
1. D. Moore, J. Leonard, D. Rus and S. Teller, “Robust Distributed
Networklocalization with Noisy Range Meas- urements,” Proceedings
of the 2nd International Confer- ence on Embedded networked sensor
systems, Baltimore, 2004, pp. 50-61.
2. R. Szewczyk, A. Mainwaring, J. Polastre and D. Culler, “An Analysis
of a Large Scale Habitat Monitoring Ap- plication,” Proceedings of
the 2nd International Confer- ence on Embedded networked sensor
systems, Baltimore, 2004, pp. 214-226.
3. J. Polastre, J. Hill and D. Culler, “Versatile Low Power Media Access
for Wireless Sensor Networks,” Proceed- ings of the 2nd International
Conference on Embedded networked sensor systems, Baltimore, 2004,
pp. 95-107.
4. M. Batalin, M. Rahimi, Y. Yu, D. Liu, A. Kansal, G. Sukhatme, W.
Kaiser, M. Hansen, G. Pottie, M. Srivastava and D. Estrin, “Call and
Response: Experiments in Sampling the Environment,” Proceedings of
ACM Sensor Systems, Los Angeles, 2004.
5. M. Rahimi, R. Pon, W. J. Kaiser, G. S. Sukhatme, D. Estrin and
M. Srivastava, “Adaptive Sampling for Environmental Robotics,”
Proceedings of IEEE International Conference on Robotics and
Automation, Los Angeles, 2004.
6. http://www.moteiv.com
7. http://www.crossbow.com
8. R. Willett, A. Martin and R. Nowak, “Backcasting: Adaptive Sampling
for Sensor Networks,” Proceedings of IEEE International Conference
on Information process- ing in sensor networks, New York, 2004, pp.
124-133.
9. R. Castro, R. Willett and R. Nowak, “Faster Rates in Re- gression via
Active Learning,” Proceedings of Neural In- formation Processing
Systems, 2005. http://homepages.cae. wisc.edu/ rcastro/ ECE-05-3.pdf
10. M. Horstein, “Sequential Decoding Using Noiseless Feed- back,” IEEE
Transactions Information Theory, Vol. 9, No. 3, 1963, pp. 136-143.
11. J. Borenstein and H. R. Everett, “Navigating Mobile Ro- bots: Sensors
and Techniques,” John Wiley & Sons, 1992.
12. M. Burnashev and K. S. Zigangirov, “An Interval Esti- mation Problem
for Controlled Observations,” Problems in Information Transmission,
Vol. 10, 1974, pp. 223-231.
A Sensing and Robot Navigation of Hybrid Sensor ... 221
12
A Study on Vehicle Detection and
Tracking Using Wireless Sensor
Networks
ABSTRACT
Wireless Sensor network (WSN) is an emerging technology and has great
potential to be employed in critical situations. The development of wireless
sensor networks was originally motivated by military applications like
battlefield surveillance. However, Wireless Sensor Networks are also used
in many areas such as Industrial, Civilian, Health, Habitat Monitoring,
INTRODUCTION
The Wireless Sensor Networks comprise of relatively inexpensive sensor
nodes capable of collecting, processing, storing and transferring information
from one node to another. These nodes are able to autonomously form a
network through which sensor readings can be propagated. Since the sensor
nodes have some intelligence, data can be processed as it flows through
the network. Sensing devices will be able to monitor a wide variety of
ambient conditions: temperature, pressure, humidity, soil makeup, vehicular
movement, noise levels, lighting conditions, the presence or absence of
certain kinds of objects, mechanical stress levels on attached objects and
so on. These devices will also be equipped with significant processing,
memory and wireless communication capabilities. Emerging low-level
and low-power wireless communication protocols can be used to networks
the sensors. This capability will add a new dimension to the capabilities of
sensors. Sensors will be able to coordinate among themselves on a higher-
level sensing task. The sensors can be deployed in any facility or area,
A Study on Vehicle Detection and Tracking Using Wireless ... 225
Characteristics of WSN
Some of the unique characteristics of a WSN include:
• Limited power they can harvest or store.
• Ability to withstand harsh environmental conditions
• Ability to cope with node failures
• Mobility of nodes
• Dynamic network topology
226 Mechatronics and Robotics
• Communication failures
• Heterogeneity of nodes large scale of deployment
• Unattended operation
Applications
The various sensor network applications include,
1) Military Applications
• Monitoring Friendly Forces, Equipments and Ammunition
• Battlefield Surveillance
• Targeting
• Battle Damage Assessment
• Nuclear, Biological and Chemical Attack Detection
2) Environmental Applications
• Forest Fire Detection
• Flood Detection
• Monitoring through Internet
• Monitoring Biodiversity
3) Habitat Monitoring Applications
• Habitat Monitoring
• Great Duck Island System
4) Health Applications
• Tele-monitoring of Physical Data
• Tracking and Monitoring Doctors and Patients in a Hospital
• Drug Administration in Hospitals
5) Home and office Applications
• Smart Homes
• Managing Inventory Control
• Home automation
• Environmental control in office buildings
6) Other Applications
• Monitoring Nuclear Reactor.
• Target Tracking
• Suspicious Individual detection: Interactive museums
A Study on Vehicle Detection and Tracking Using Wireless ... 227
Magnetometer Less susceptible than loops to stresses of traffic. Installation requires pavement cut.
(two-axis fluxgate Insensitive to inclement weather such as snow, rain, Improper installation decreases pavement life.
magnetometer) and fog. Installation and maintenance require lane closure.
Some models transmit data over wireless radio fre- Models with small detection zones require multiple units for
quency (RF) link. full lane detection.
Magnetic (induction Can be used where loops are not feasible (e.g., bridge Installation requires pavement cut or boring under roadway.
or search coil mag- decks). Cannot detect stopped vehicles unless special sensor layouts
netometer) Some models are installed under roadway without and signal processing software are used.
need for pavement cuts. However, boring under
roadway is required.
Insensitive to inclement weather such as snow, rain,
and fog.
Less susceptible than loops to stresses of traffic.
Microwave radar Typically insensitive to inclement weather at the Continuous wave (CW) Doppler sensors cannot detect
relatively short ranges encountered in traffic man- stopped vehicles
agement applications.
Direct measurement of speed.
Multiple lane operation available.
Active infrared (laser Transmits multiple beams for accurate measurement Operation may be affected by fog when visibility is less than
radar) of vehicle position, speed, and class. 20 feet (ft) (6 m) or blowing snow is present.
Multiple lane operation available. Installation and maintenance, including periodic lens clean-
ing, require lane closure
Passive infrared Multizone passive sensors measure speed. Passive sensor may have reduced vehicle sensitivity in
heavy rain, snow and dense fog.
Some models not recommended for presence detection.
Ultrasonic Multiple lane operation available Environmental conditions such as temperature change and
Capable of over height vehicle detection. extreme air turbulence can affect performance. Temperature
Large Japanese experience base. compensation is built into some models.
Large pulse repetition periods may degrade occupancy
measurement on freeways with vehicles travelling at moder-
ate to high speeds.
Acoustic Passive detection. Cold temperatures may affect vehicle count accuracy.
Insensitive to precipitation. Specific models are not recommended with slow-moving
Multiple lane operation available in some models. vehicles in stop-and-go traffic.
Video image proces- Monitors multiple lanes and multiple detection Installation and maintenance, including periodic lens clean-
sor zones/lane. ing, require lane closure when camera is mounted over
Easy to add and modify detection zones. roadway (lane closure may not be required when camera is
Rich array of data available. mounted at side of roadway).
Provides wide-area detection when information gath- Performance affected by inclement weather such as fog,
ered at one camera location can be linked to another. rain, and snow; vehicle shadows; vehicle projection into
adjacent lanes; day-to-night transition; vehicle/road con-
trast; and water, salt grime, icicles, and cobwebs on camera
lens.
Reliable night-time signal actuation requires street lighting.
Requires 30- to 50-ft (9- to 15-m) camera mounting height
(in a side-mounting configuration) for optimum presence
detection and speed measurement.
Some models are susceptible to camera motion caused by
strong winds or vibration of camera mounting structure.
Generally cost effective when many detection zones within
the camera field of view or specialized data are required.
Intrusive Sensors
Intrusive sensors are those that need to be installed under the pavement,
in saw-cuts or holes on the roads. Popular intrusive sensors include
inductive loops, magnetometers, micro loop probes, pneumatic road
tubes, piezoelectric cables and other weigh-in-motion sensors. The main
advantage of these sensors is their high accuracy for vehicle detection while
A Study on Vehicle Detection and Tracking Using Wireless ... 229
the drawbacks include the disruption of traffic for installation and repair,
resulting in high installation and maintenance cost.
Non-Intrusive Sensors
To overcome the disadvantage of intrusive sensors, non-intrusive sensors
are developed eg, above ground vehicle detection sensors. Above ground
sensors can be mounted above the lane of traffic or on the side of a roadway
where they can view multiple lanes of the traffic at angles perpendicular to or
at a slanting angle to the flow direction. Technologies used in aboveground
sensors include video image processing (VIP), microwave radar, laser radar,
passive infrared, ultrasonic, passive acoustic array, and combinations of
these sensor technologies. However, these non-intrusive sensors tend to be
large size and power hunger.
Acoustic Sensors
The acoustic sensor in the Smart Dust sensor node is a condenser type
microphone. The schematic for a typical condenser acoustic sensor is shown
in Figure 4. It includes a stretched metal diaphragm that forms one plate of a
capacitor. A metal disk placed close to the diaphragm acts as a backplate. A
stable DC voltage is applied to the plates through a high resistance to keep
electrical charges on the plates. When a sound field excites the diaphragm,
the capacitances between the two plates vary according to the variation in
the sound pressure. The change in the capacitance generates an AC output
proportional to the sound pressure, which shows the ultra low-frequency
pressure variation. A high-frequency voltage (carrier) is applied across the
plates and the acoustic sensor output signal is the modulated carrier. The
photo in the right of Figure 4 shows the Panasonic
Magnetic Sensors
Magnetic sensors differ from most other detectors in that they do not
directly measure the physical property of interest. Devices that monitor
properties such as temperature, pressure, strain, or flow provides an output
that directly reports the desired parameter (Figure 6). Magnetic sensors, on
the other hand, detect changes, or disturbances, in magnetic fields that have
been created or modified, and from them derive information on properties
such as direction, presence, rotation, angle, or electrical currents. The output
signal of these sensors requires some signal processing for translation into
the desired parameter. Although magnetic detectors are somewhat more
difficult to use, they do provide accurate and reliable data without physical
contact.
The Honeywell HMC1002 magnetometer on the MICA sensor board is
a magnetoresistive sensor. Honeywell magnetic sensors and magnetometers
offer complete magnetic field sensing solutions that are highly accurate, and
allow for easy integration for virtually any application. Figure 7 shows the
honeywell magnetic sensors.
The Anisotropic Magneto Resistive (AMR) sensor is one type that has a
wide Earth’s field sensing range and can sense both the strength and direction
of the Earth field. The AMR sensor is made of a nickel-iron (Permalloy)
thin film deposited on a silicon wafer and patterned as a resistive strip. The
strip resistance changes about 2%-3% when a magnetic field is applied.
Typically, four of these resistive strips are connected in a wheatstone bridge
configuration so that both magnitude and direction of a field along a single
axis can be measured. The key benefit of AMR sensors is that they can be bulk
manufactures on silicon wafers mounted in commercial integrated circuit
packages. Figure 8 shows a typical vehicle magnetic signal waveforms.
Magnetic Field
A magnetic field is a vector field that surrounds magnets and electric currents,
and is detected by the force it exerts on moving electric charges and on
magnetic materials. When placed in a magnetic field, magnetic dipoles tend
to align their axes parallel to the magnetic field. Magnetic fields also have
their own energy with an energy density proportional to the square of the
field intensity.
Magnetic Data
Magnetic data is the term used for data that is acquired from magnetic (as
opposed to optical) motion capture systems. A central magnet is used to
create a field in which sensors can determine their position and orientation
as they move about in the field. Raw magnetic data has no hierarchy
234 Mechatronics and Robotics
information; the sensors do not know where they are relative to the other
sensors, and know their position in the magnetic field. The file formats for
magnetic motion capture data reflect this. In this respect, they are somewhat
similar to the BVA file format.
The problem associated with magnetometer vehicle detection is similar
to the acoustic sensors but the magnetic signals are much cleaner than
acoustic signals. The magnetometers available today can sense magnetic
fields within the earth’s field-below 1 gauss.
They can be used for detecting the vehicles, which are ferrous objects
that disturb the earth’s field. The earth’s field provides a uniform magnetic
field over wide area in the scale of kilometers and even a ferrous object can
create a local disturbance in this field.
This local field disturbance can be sensed by the magnetometers for
vehicle detection. Figure 9 shows the disturbance of earth’s magnetic field
by a car. After presenting the basic principles of sensors, the characteristics
of the measured acoustic and magnetic signals will be studied and algorithms
are proposed for reliable low cost vehicle detection.
Acoustic Signals
The acoustic signature is made up of a number of individual elements.
These include:
A Study on Vehicle Detection and Tracking Using Wireless ... 235
Magnetic Signals
Magnetic Signals are changes that happen in the magnetic field caused by
the movement of ferrous objects.
Basic Fusion
This algorithm is presented first since it has an extremely simple design and
implementation and serves as an ideal algorithm to use as the foundation
for explaining other more complex detection algorithms. This algorithm
collects the captured data from each of the nodes in a network. Then, to
236 Mechatronics and Robotics
remove any invalid data provided by faulty nodes, the largest and smallest
data values are dropped. The remaining data is then averaged together. If this
average exceeds a predetermined threshold, then the system is notified that a
target has been detected. This algorithm has two versions in which the data
provided by each node differs. The data is either the original data collected
by the node or it is a detection decision already made by the individual
nodes using their original data [5].
Hybrid
This algorithm is similar to the basic fusion algorithm. In that a detection
decision is based on data collected by numerous nodes throughout the
network. However, this algorithm differs in that either a yes/no target present
bit decision is sent or the collected data is sent. In this algorithm, each node
is given two threshold values for a previously determined detection signal
type, such as sound. One of these threshold values determines that if the
data collected by the node exceeds this value then a target has been detected.
Correspondingly, the other threshold is a minimum in that if the data is below,
then no target is present. If either of these two thresholds is met then a bit
message indicating the presence or lack of a target is transmitted. However,
if neither of these thresholds is exceeded then the raw data collected by the
node is passed on to the base station where it is used with the data of other
nodes to determine if a target is present or not [18].
Predicate Clustering
In this algorithm nodes are divided into groups called clusters in which each
group has a cluster head that all sub nodes of the cluster report their data to.
Typically, one cluster is formed around a single target with multiple clusters
being in the network if multiple targets are present. Cluster heads are chosen
by all nodes within the area of a detected target “electing” to become the
head by broadcasting their intent to become a cluster head. The node with
the peak sensor reading is chosen as the cluster head [6]. Once cluster heads
are chosen, the remaining nodes choose group to join by running a decision
predicate on their captured sensor data and the data of their neighbours. One
particularly interesting feature of this algorithm is that a cluster head can
be responsible for multiple targets. When it is detected, multiple targets are
going to enter the same cluster, both clusters provide the target information
to the “merged” cluster head that is now responsible for tracking both targets.
A Study on Vehicle Detection and Tracking Using Wireless ... 237
Grids
In this approach the sensor network is broken up into “virtual grids” where
each grid contains points. For each grid point, it is determined which nodes’
sensing area covers the grid point and then these nodes create a schedule so
that they all collaboratively cover this grid point at alternating times. So at
any point of time‘t’ at least one node’s sensing area covers grid point ‘x’.
Since each node covers multiple grid points, each node will contain several
sensing schedules to keep in sync so that at any point of time all grid points
are covered by at least one node. In order to do this, all nodes first determine
the schedules for each grid point, then schedule all their individual point
schedules together so that each node is not constantly waking and sleeping
instead sleeps a little and then wakes for a while. Once each node determines
its own individual schedule, it must be correlated with those of its neighbours
so that as few as possible nodes are awake at any point of time but all grid
points are covered at all times. Once the schedules are completed for each
node, the nodes start to follow their schedule, going through states of being
asleep and awake in order to detect a new target [7].
Wave
The Wave protocol has the same goal in mind as that of the Grid algorithm,
detecting targets while conserving energy by not having all the nodes on
all the time. The basic idea of the wave algorithm is to wake up a series
of nodes at certain times so that only nodes that are awake at one point of
time are in one concentrated area. Then, over time, this section moves over
the network by putting nodes to sleep and waking new ones. This gives the
appearance of an awake section of nodes flowing over the entire network
like a wave [8].
This algorithm proposes three types of waves:
• Two straight lines of nodes are awake at the same time, with the
lines being parallel to each other. The two lines start on opposite
sides of the network from each other and gradually work towards
each other to meet in the centre of the network.
• A wide line of nodes are awake at the same time and this line
moves across the network from one edge to the opposite edge.
• Waves of nodes are awake and the wave traverses from one edge
of the network to the other. This approach differs from the wide
line in that the boundary shape is curvy or wave like instead of
238 Mechatronics and Robotics
Patrolling
Patrolling algorithms are similar to the wave algorithms in that only a series
of nodes are awake at one point of time. This algorithm just chooses which
nodes to wake differently. This can be done on an on-demand basis or a
coverage-orientated basis. The difference between the two is that on-demand
is some predetermined path that is currently interesting whereas coverage-
orientated just repeatedly watches an area [9].
The algorithm for both approaches is the same. First, a patrol is defined
as a set of information containing the patrol speed, duration, iteration period
and path. Once this information is collected, a Patrol Host node is selected,
which is usually the first node on the patrol path. The duty of the Patrol
Host is to periodically transmit the patrol data along with the current patrol
time. The idea is that over the time the patrol virtually traverses down the
specified path. As the patrol advances, nodes that are closest to the patrol
path wake up and collect information. Nodes determine if they are near the
patrol path by using the patrol information that the Patrol Hosts broadcast.
Additionally, as the patrol moves, new Patrol Hosts must be elected. As the
patrol virtually moves out of the current patrol host’s range, a new patrol
host is selected on the path which then takes over the patrol host duties and
the previous patrol host can now sleep. This process of waking nodes and
handing off patrol host duties continues until the end of the path is reached.
Once the path end is reached, the patrol is reversed in order to traverse the
path back to the original starting point. This process continues according to
the duration and iteration parameters are set in the original patrol parameters
[9].
Mobile Nodes
Until this point, all the target detection algorithms discussed so far have
assumed that the sensor nodes remain stationary. This section looks at a
A Study on Vehicle Detection and Tracking Using Wireless ... 239
few target detection algorithms in which the nodes are mobile within the
network. Even though these algorithms allow nodes to be mobile, their goal
remains the same as the previously mentioned target detection algorithms.
These algorithms are still trying to detect targets in a timely and efficient
manner while trying to cover as much of the network as possible using as
little energy as possible. The theory behind allowing nodes to be mobile is
that they may now detect a target that would have otherwise gone undetected
in a stationary node network [10].
Past Detections
Like other algorithms, the network in this algorithm is divided into cell
areas. Each cell is unique in that it maintains the current state and past
history of targets that have been in the cell. This includes the location in
which the target was first detected within in the cell. Using this data, the
system calculates the areas of the cell that are most likely to contain a target
based on the locations of previously detected targets [11].
Once these locations are determined, each node is assigned an area
to watch and a priority. These priorities are used in the “coordination
mechanism” part of the algorithm. The coordination part of the algorithm is
responsible for altering the path nodes so that their observation paths overlap
as little as possible. For example, the node that watches the area where
targets are the most likely to occur is given the highest priority of the system
and its path is not altered. A second node, whose priority is slightly less than
the first node, whose path also slightly intersects that of the first node has its
path altered so that the two paths will not cross. This path altering continues
throughout the system from highest priority down to lowest priority until all
node paths are altered to prevent too much overlapping of observation areas
[11].
Collaborative Coverage
This algorithm is little different in that the paths the nodes traverse are not
coordinated in any way. Each node is allowed to create and follow its own
arbitrary path. While the node is travelling on this path at different locations
and times, it collects data and saves. The sensed data are saved with time
and location details [12]. While the nodes move across their paths they share
their collected data with those nodes they come into radio communication
with. Data is usually shared on a request basis. Each node is provided a
coverage threshold and confidence threshold. Actual target detection in this
240 Mechatronics and Robotics
system is determined by examining the data taken at several points that are
close to each other. If the combined data collected by this group of close
points exceeds some predetermined threshold then it is determined that a
target is detected at this area. Note that the value of this detection threshold
will determine how sensitive/ insensitive the system is at detecting targets.
Simple Triangulation
This algorithm is an extremely simple design and implementation and
A Study on Vehicle Detection and Tracking Using Wireless ... 241
Clusters
The cluster target tracking algorithm has been widely discussed through
many research papers. In this section, the basic idea of target tracking using
clusters is discussed, followed by an overview of the variances in each of the
different cluster algorithm research papers.
Dynamic Clustering
Like the super node algorithm, this algorithm also assumes that cluster head
nodes have more power than normal sensor nodes. However, sensor nodes
are not assigned to clusters in this algorithm. Instead, they are invited to join
a cluster by the cluster head. The cluster head does this by broadcasting a
join message that includes the time and signature of the target the cluster
head detected [19]. Those sensor nodes that have stored data that matches
the data in the broadcast message respond to the broadcast by sending their
captured data to the cluster head node. Interestingly, the cluster head only
waits for a certain number of replies and when the required numbers of
replies are received, the cluster head calculates the target’s location. Unlike
the super node algorithm, this algorithm has one active cluster head at a
time. In other words, cluster heads do not work together [22].
Cell Collaboration
Like predicate clustering all nodes in this algorithm are of the same type.
Nodes are formed into clusters (called cells) in which all nodes within the
cells collaboratively decide when a target has entered their cell and if they
should track it [22]. It is interesting to note that cells can have different
sizes; the size of cells is determined by the observed velocity of the target.
So cell size will increase for faster moving targets.
A Study on Vehicle Detection and Tracking Using Wireless ... 243
Probabilistic Localization
This algorithm deviates even further from the previously discussed
algorithms. Cluster heads are special high powered nodes that know the
location of every node within their cluster [23]. This algorithm takes
advantage of the fact that the cluster head knows these locations. When
sensor nodes detect a target, they send a very small notification message to
their cluster head and store the target location data, time and other relevant
data in their local memory. Upon receiving several notifications, the cluster
head performs a probabilistic localization algorithm to determine which
sensor nodes to query saved data from. In other words, the cluster head
runs an algorithm that helps to predict which sensor nodes are closest to the
target. Therefore, the cluster head has the best data to calculate an accurate
estimation of the target location. When the cluster head determines which
nodes to query, they are asked for the data they saved in their local memory
and the cluster head uses this information to calculate the location of the
target.
Rooted Tree
Related to the clusters tracking algorithm is the idea of the rooted tree
tracking algorithm. This algorithm differs somewhat in that instead of
having multiple cluster heads, there is only one head node and it is referred
to as the root. The root node is the node that is closest to the target.
244 Mechatronics and Robotics
Particle Filtering
The basic idea behind a particle filtering tracking algorithm is that numerous
object state descriptions are saved that contain data necessary to calculate
the target position at a certain time. These state descriptions are referred to
as particles and each particle has its own weight. The weight of a particle
determines how much the data it contains will contribute to the location
estimation of an object. When new particles are created, the weights of
the pre-existing particles are adjusted and then all the particles are used to
calculate a new target location [26]. Eventually, particles with weights that
are below a certain threshold are eliminated as duplicate particles.
Some of the particles filtering algorithms differ a little. Certain algorithms
have all the particles stored in a central node and this node does all the target
location processing [28]. Other algorithms distribute the stored particles
across the network nodes.
ANALYSIS OF ALGORITHMS
This section analyses each target tracking algorithm mentioned with a focus
on how it could be improved by combining it with one or more of the ideas
presented in the target detection algorithms.
Simple Triangulation
The main problem with this algorithm lies in the fact that when a node
has detected a target and hears two other nodes broadcast the same
target detection, the nodes starts to perform the triangulation localization
A Study on Vehicle Detection and Tracking Using Wireless ... 245
calculation. This means that at the same point of time there will be multiple
nodes calculating the same target location. This is a waste of computation and
energy resources. This algorithm needs to be improved so that when a node
starts to calculate the location of the target, it first broadcasts its intention to
calculate the location to all the other nodes, similar to the data sharing that
is done in the collaborative coverage mobile detection algorithm. Then all
the other nodes that would have performed this same calculation will not
perform it once they hear the broadcasted intention of the first node.
Basic Cluster
The majority of the cluster algorithms had all the nodes on during their
target detection phase. It would be better if they instead adopt a technique
similar to that used in the wave detection algorithm to detect targets. Since
each area of the network is already divided into clusters, each node within
the cluster could take turns monitoring the area for the arrival of a target
using one of the wave detection methods mentioned. This would increase
the lifetime of the network.However, using one of the mobile detection
algorithms is not recommended since the nodes have to register with the
cluster head. Depending on the mobile detection algorithm used, the nodes
would constantly be moving in and out of clusters. This would cause an
increase in the amount of message communications with these nodes in
order keep track of the nodes within the cluster. Mobile detection algorithm
of this type can be applied to the situation only if the mobility of the nodes
is limited and the node always stay within the same cluster.Additionally, the
cluster algorithms could benefit from the basic fusion detection algorithm.
Once a target has been detected and all the information is passed to the
cluster head the cluster head could use the technique used in the basic fusion
algorithm to eliminate any data provide by faulty nodes. This would help
make the network less vulnerable to faulty nodes.
Dynamic Clustering
The biggest problem with this algorithm is that it will not work with mobile
nodes. The network calculates the sensor node’s locations at start-up. Since
these locations need to be known by the cluster head and are calculated
only once, this prevents any nodes from moving or even being added to
the network. This problem could be easily fixed by periodically refreshing
the locations of the sensor nodes, although this will decrease the network’s
lifetime. So it is better to use a detection algorithm and constantly refreshing
the nodes’ location isn’t an optimal solution is better in this case.
246 Mechatronics and Robotics
Cell Collaboration
This clustering algorithm is already perfectly setup to use the grid detection
algorithm for detection of targets. Since the nodes in this algorithm are
already divided into cells, the basic structure to use the grid algorithm is
already in place. All that needs to be done is to add the grid algorithm for
the detection part. Once, this is done the algorithm would be more efficient
since the nodes would be alternating sleep schedules with each other in order
to detect the target instead of all being awake at the same time.
Probabilistic Localization
This algorithm already sounds like that it has incorporated parts of the hybrid
detection algorithm in it. When a target is detected, this algorithm sends a
small message to the cluster head. Similarly, when a target is detected in the
hybrid algorithm a small yes/no bit message is sent to the decision maker in
the network, or cluster head.
Then the cluster head uses this data to determine which nodes to query
for further data. In determining which nodes to query, it would be a good
idea to use an averaging technique similar to that used in the basic fusion
detection algorithm in order to eliminate the impact of any data provided by
a faulty node. The similarities to the hybrid algorithm continue. When the
cluster head determines which nodes to query, the nodes respond back to the
cluster head with more detailed information.
However, this algorithm should not be combined with a mobile node
detection algorithm because of the way it decides which nodes to query.
In this algorithm, the cluster heads create probability tables in order to
determine which sensors to query [23]. The use of a table is acceptable until
nodes move outside or are added to the system. When this happens, the
probability table must be updated. Therefore, the use of a mobile detection
scheme would cause a lot of additional computational overhead since these
tables are to be constantly refreshed.
Rooted Tree
One of the basic problems with the rooted tree technique is that it doesn’t
provide any fault tolerance facts. In order to make this algorithm less
vulnerable to problems caused by faulty nodes, it should use an averaging
technique similar to that used in the basic fusion detection algorithm.
Additionally, this algorithm encounters problems when the velocity of
the target is very high. As the velocity of the target increases so does the
number of times in which the tree must be reconfigured [25]. This will cause
problems when a threshold velocity of the target is reached in which the
tree can no longer be generated quickly enough to keep up with the moving
target. This problem could be eliminated by altering this algorithm to work
more like the patrolling algorithm. Instead of constantly rebuilding the tree,
the tree is treated as a path that can have steps (or nodes) along the path
added or deleted from it. This way when the target moves, just add a new
node to the path and remove any one that is too far away to contribute any
longer. The basic nodes needed to collect data are still members of the path
and new ones are added. It saves an enormous amount of time and message
communication needed to periodically rebuild the tree structure each time.
Particle Filtering
In order to track targets, this algorithm saves an enormous amount of data
in the form of particles. Over time, these particles take up a large amount
of memory space. The amount of memory actually used could be decreased
if this algorithm is to adopt the idea behind the TBD algorithm. Bayesian
Estimation Track Before Detect (TBD) based on particle filtering gives high
performance and hybrid detection algorithm stops recording data once the
data it currently contains exceeds a threshold to know whether a target is
present or not. The particle filtering algorithm could use the same principle.
Once enough particle data is collected to exceed some threshold, the node
would stop collecting data. This would not only decrease the amount of
memory used but also decrease the amount of work to be done to save energy.
SUMMARY OF ALGORITHMS
It is observed that, one of the best target detection and tracking algorithms
is the combination of the Bayesian Estimation Track before Detect (TBD)
Estimator and particle filtering algorithm. The fact that the TBD covers the
network area in a short amount of time and enables nodes to sleep at times
saves a huge amount of energy. This enables a target to be detected, without
too much delay, in an energy efficient manner.
The particle filtering algorithm is also on the right track since the pre-
existing particles are adjusted when the new particles are created and hence
particles with weights that are below a certain threshold are eliminated as
duplicate particles, which also conserves energy. Additionally, the use of
several central localization calculating points instead of one helps to reduce
the vulnerability of the system to an attacker. It eliminates the one point
of failure problem. It also shares the burden of energy usage that one node
would face across several nodes.
Particle filtering approach allows all the nodes to share the energy burden.
Hence particle filtering algorithm combined with Bayesian Estimation Track
before Detect (TBD) Estimator algorithm can be used for target detection
and tracking.
ACKNOWLEDGMENT
The authors would like to thank the Armament Research Board (ARMREB-
DRDO) for supporting this Research project by funding.
250 Mechatronics and Robotics
REFERENCES
1. Y. C. Tseng, S. P. Kuo, H. W. Lee, and C. F. Huang, “Location tracking
in a wireless sensor network by mobile agents and its data fusion
strategies,” Lecture notes in computer science, Springer Link, pp. 2–3,
2003.
2. E. Y. Luz and A. Mimbela, “Summary of vehicle detection and
surveillance technologies used in intelligent transportation systems,”
The Vehicle Detector Clearinghouse, Southwest Technology
Development Institute (SWTDI) at New Mexico State University
(NMSU), Fall 2007. http://www.nmsu.edu/traffic/.
3. S. Coleri, M. Ergen, and T. J. Koo, “Lifetime analysis of a sensor
network with hybrid automata modelling,” Processings of ACM
International Workshop on Wireless Sensor Networks and Applications
(Atlanta, GA), pp. 98–104, 2002.
4. D. Jiagen, S. Y. Cheung, C. W. Tan, and V. Pravin, “Signal processing of
sensor node data for vehicle detection,” IEEE Intelligent Transportation
Systems, pp. 70–75, 2004.
5. C. Thomas, K. S. Kewel, and R. Parameswaran, “Fault tolerance in
collaborative sensor networks for target detection,” IEEE Transactions
on Computers, pp. 320–333, 2004. http://citeseer.ist.psu.edu/
clouqueur03fault.html.
6. F. Qing, Z. Feng, and L. Guibas, “Lightweight sensing and
communication protocols for target enumeration and aggregation,”
Proceedings of the 4th ACM international symposium on Mobile ad
hoc networking & computing, pp. 165–176, 2003.
7. Y. Ting, H. Tian, and A. S. John, “Differentiated surveillance for
sensor networks,” Proceedings of the 1st international conference on
Embedded networked sensor systems, pp. 51–62, 2003. http://www.
cs.virginia.edu/papers/ diff_surveillance_sn_p51-yan.pdf.
8. S. S. Ren, Q. Li, H. N. Wang, and X. D. Zhang, “Design and analysis of
wave sensing scheduling protocols for object-tracking applications,”
Lecture notes in computer science, SpringerLink, pp. 228–243, 2005.
http://www.cs. wm. edu/~liqun/paper/dcoss05.pdf.
9. C. Gui and P. Mohapatra, “Virtual patrol: a new power conservation
design for surveillance using sensor networks,” Information Processing
in Sensor Networks, pp. 246–253, 2005.
10. B. Y. Liu, P. Brass, O. Dousse, P. Nain, and D. Towsley, “Mobility
A Study on Vehicle Detection and Tracking Using Wireless ... 251
13
Deployment of Pre-Industrial
Autonomous Microbe Sensor in Saudi
Arabia’s Injection Seawater System
ABSTRACT
Microbial growth in water injection systems can lead to many problems,
including biofouling, water quality deterioration, injectivity loss, microbial
corrosion, and reservoir formation damage. Monitoring of microbial
INTRODUCTION
Saudi Arabian Oil Company (Saudi Aramco) has the largest seawater injection
system in the world. The growth of microorganisms and the formation
of biofilm on the pipeline’s inner surfaces and process equipment lead to
an array of challenges, including biofouling, water quality deterioration,
infectivity loss, reservoir plugging, and microbiologically influenced
corrosion (MIC) [1] [2] . A critical cornerstone in an effective mitigation
strategy is the routine surveillance of the system microbiology, enabling
operators to apply and adjust microbiological countermeasures accordingly.
The microbial number and activity are traditionally monitored with
conventional growth methods that require manual sampling and handling.
Monitoring the microbial activity with conventional methods in this large
system, especially at the remote locations, presents a long-standing challenge
to Saudi Aramco. In an attempt to ensure the injection water quality and
maintain the integrity of injection pipelines network, Saudi Aramco and the
Danish Technological Institute (DTI) embarked on a three-phase project
Deployment of Pre-Industrial Autonomous Microbe Sensor in ... 255
PERFORMANCE ASSESSMENT
Lab validation
After construction of two pre-industrial AMS units, the AMS was thoroughly
tested and validated in the laboratories at DTI premises. The repeatability
and functionality of the units were tested with respect to procedures as well
as software. The major improvements in lab validation and their effect on
AMS performance are summarized in Table 2.
After implementation of above improvements on AMS units, a correlation
of fluorescence signals from AMS and direct counting of microbial cell
number under microscope after DAPI staining was established (Table 3).
Improvement Effect
Decrease of Pump A rate to 3 Eliminate air from surroundings into sy-
mL・min−1 ringe A
New loading procedure for internal Improved mixing of internal control (IC)
DNA control DNA and Milli-Q water
Increased total volume of diluted in- Improved loading of IC DNA
ternal DNA control to sample volume
New dye loading procedure. Dead Accurate dye dosing
space in Pump B valve used as dye
dosing volume
Selection of incubation conditions of Improved staining of cells with SYBR
13 min at 60˚C Green I
No air is moved through flow cell Eliminate introduction of air bubbles in
flow cell
Flow direction of flow cell changed Remove potential air bubbles in flow cell
New flow rate through flow cell Optimal flow rate for maximum reading
Increased rate of dispensing and Minimize total AMS measurement cycle
aspirating with Pump A for all non- time to 75 min
critical steps
Incorporated data export function For automatic plot of raw and calibrated
readings on graph in Excel
Table 3. Laboratory calibration data and equation for AMS1 and AMS2.
AMS1 AMS2
Fluorescence Cells mL −1
Fluorescence Cells mL−1
count (AMS) count (AMS)
Sterile filtered sys- 154 ± 7, n = 3 98 ± 3, n = 3
tem water (blank)
LLOD 167 6.0 × 103 104 8.8 × 102
LLOQ 188 1.6 × 104 114 2.3 × 103
Regression equa- 473X − 72,788 146X − 14,265
tion
Regression, R2 0.9902 0.9490
Linear range 6.4 × 10 - 5.0 × 10
3 6
1.1 × 104 - 2.4 × 106
Improvement Effect
Volume of sample 4 ml
Volume of IC DNA 4 ml (0.1 ml concentrated IC discarded at first,
0.2 ml IC in 3.8 ml sterile H2O)
SYBR Green I Diluted 1:20 in Dimethyl sulfoxide (DMSO).
Volume equals to dead space in Pump B (20 μL)
Incubation of sample 13 min at 60˚C
Incubation of IC Less than one minute at ambient temperature
Flow rate and time 3.5 μL・s−1 of 180 sec and 10 sec
Measurement wavelength inter- 520 - 530 nm
val
Fluorescence count Average of readings every one sec for 180 and 10
sec
Background reading (auto-fluo- To be subtracted from sample reading. Fluorescence
rescence) reading of sample without dye, for correction of any
auto-fluorescence of the sample water. Background
reading is not subtracted from internal control readings.
Figure 3. Relation between cell number and fluorescence count of AMS1 (blue)
and AMS2 (red) using enrichment culture from the system, diluted in sterile
filtered system water. All AMS readings are performed with an AMS interior
temperature of 30˚C.
A linear regression between AMS fluorescence measurement and DAPI
count for converting AMS fluorescence count to cell numbers in the sample
is summarized in Table 5. Compared to the calibration performed in the
laboratory, these data show that the fluorescence count relative to the cell
number is lower in the field, and therefore the regression equation in the lab
at 25˚C is not recommended for field trial. The field regression equations for
each AMS unit were loaded in the specific AMS units, and used for the data
conversion from fluorescence count to cell number.
(1)
where the fitted parameters are A = 2.4 × 104 s−1, B = −0.10˚C−1, x0 =
4.16˚C, C = 2.49 × 103 s−1, and x = AMS inside temperature (˚C).
Deployment of Pre-Industrial Autonomous Microbe Sensor in ... 263
The function f(x) between temperature and fluorescence intensity was used
to temperature compensate internal control (IC) readings. The compensation
is done by multiplying the IC reading with the ratio f(30)/f(T), where 30 is
the reference temperature, and T is the AMS inside temperature of the given
data point.
Table 5. Field calibration data and equation for AMS1 and AMS2.
AMS1 AMS2
Fluorescence Cells mL −1
Fluorescence Cells mL−1
count (AMS) count (AMS)
Sterile filtered sys- 114 ± 12, n = 3 123 ± 7, n = 3
tem water (blank)
LLOD 136 1.9 × 104 137 1.2 × 104
LLOQ 173 5.0 × 104 160 3.2 × 104
Regression equa- 843X − 96,020 878X − 107,935
tion
Regression, R2 0.9939 0.9940
Melting of the double stranded (ds) Oligo DNA (internal control) would
result in a lower binding between the SYBR molecule and the DNA, as
SYBR Green I only binds to ds DNA [6] [16] . The oligo DNA as well as the
internal control buffer have been designed to withstand high temperatures,
and based on calculations from the Oligo DNA supplier (Integrated DNA
Technologies), the melting temperature (Tm, 50%) is 73˚C for the conditions
provided from the internal control buffer solution (Figure 7).
The spectrometer must be cooled for proper readings, and cooling has been
included in the design of the AMS. To test whether the AMS spectrometer is
affected by AMS inside temperature, fluorescein was injected manually into
the flow cell at high and low temperatures (40˚C and 25˚C, respectively).
Fluorescein is a dye that emits fluorescent light with excitation and emission
spectra comparable to those of SYBR Green when bound to DNA [17]. The
effect of AMS inside temperature on spectrometer readings of fluorescein at
two concentrations was not detected (data not shown).
Field measurements
The AMS1 measurement data, from October 2016 to February 2017, on the
fluorescence count of the IC DNA and microbial cell number in the seawater
are plotted in Figure 8 and Figure 9, respectively. Figure 8 shows that the
IC readings are relatively stable in this 4-month period, after temperature
compensation of fluorescence intensity, verifying the stability of both IC
DNA and SYBR Green dye, and the functionality of AMS in the field. Figure
9 shows that the majority of the measurements of seawater cell number are
above the detection limit. When system water fluorescence count has been
temperature compensated, the cell numbers from winter are in the range
comparable to fall cell numbers.
Temperature effect on fluorescence intensity is determined to be
significant when SYBR Green I dye is used to stain IC DNA and cell DNA in
seawater samples. As a result, temperature compensation for all fluorescence
counts of both IC DNA and microbial cells is required. The temperature
compensation factor was incorporated in the software for converting the raw
data to temperature compensated data for data exporting and plotting.
LESSONS LEARNED
During the field testing and validation, a number of technical challenges
have been either encountered or identified, and the solutions have been
implemented to address these challenges (Table 6). Among the challenges,
266 Mechatronics and Robotics
the two main threats are the pressure regulator clogging and leakage from
the fluid handling system, which are discussed in detail in the following
text, to illustrate how the technical challenges have been resolved.
The pre-industrial AMS is equipped with a restrictor to reduce the side-
stream pressure from 10 bars to less than 1 bar before water enters the AMS
system. The restrictor is equipped with a pressure gauge to monitor inlet
pressure. Since the seawater pipeline pressure is much higher, it is important
that another pressure regulator is installed on the side-stream to reduce the
pressure to about 10 bars for prevention of damage to AMS hardware.
Figure 7. Hybridization profile of the Oligo DNA used as internal control. The
total DNA is intact up to 70˚C, and the Tm 50% is 73˚C.
ACKNOWLEDGMENTS
The authors would like to acknowledge the Saudi Arabian Oil Company
(Saudi Aramco) and the Danish Technological Institute for granting
permission to publish this paper.
Deployment of Pre-Industrial Autonomous Microbe Sensor in ... 271
REFERENCES
1. Ren, H., Wang, W., Liu, Y., Liu, S., Lou, L., Cheng, D., He, X., Zhou,
X., Qiu, S., Fu, L., Liu, J. and Hu, B. (2015) Pyrosequencing Analysis
of Bacterial Communities in Biofilms from Different Pipe Materials
in a City Drinking Water Distribution System of East China. Applied
Microbiology and Biotechnology, 99, 10713-10724. https://doi.
org/10.1007/s00253-015-6885-6
2. Moniee, M.A., Juhler, S., Sorensen, K., Abeedi, F.N., Lundgaard, T.
and Sanders, P.F. (2014) A Review of Saudi Aramco’s Water Flooding
System and Methods for Monitoring Microbial Activity. Proceedings
15th Middle East Corrosion Conference, Bahrain, 2-5 February 2014,
1-19.
3. Moniee, M.A., Juhler, S., Sorensen, K., Zhu, X.Y., Lundgaard, T.,
Abeedi, F.N. and Sanders, P.F. (2016) Laboratory-Scale Evaluation
of Single Analyte Bacterial Monitoring Strategies in Water Injection
Systems. Journal of Sensor Technology, 6, 11-26. https://doi.
org/10.4236/jst.2016.62002
4. Moniee, M.A., Zhu, X.Y., Tang, L., Juhler, S., Nuwaiser, F.I., Sanders,
P.F. and Abeedi, F.N. (2016) Optimization of DNA Staining Technology
for Development of Autonomous Microbe Sensor for Injection
Seawater Systems. Journal of Sensor Technology, 6, 27-45. https://doi.
org/10.4236/jst.2016.63003
5. Moniee, M.A., Zhu, X.Y., Tang, L., Nuwaiser, F.I., Voigt, N.V., Sanders,
P.F., Abeedi, F.N. and Habboubi, H.H. (2016) Validation of Autonomous
Microbe Sensor Prototype for Monitoring of Microorganisms in
Injection Seawater Systems. Journal of Sensor Technology, 6, 81-100.
https://doi.org/10.4236/jst.2016.64007
6. Dragan, A.I., Pavlovic, R., McGivney, J.B., Casas-Finet, J.R., Bishop,
E.S., Strouse, R.J., Schenerman, M.A. and Geddes, C.D. (2012) SYBR
Green I: Fluorescence Properties and Interaction with DNA. Journal
of Fluorescence, 22, 1189-1199. https://doi.org/10.1007/s10895-012-
1059-8
7. Noble, R.T. and Fuhrman, J.A. (1998) Use of SYBR Green I for Rapid
Epifluorescence Counts of Marine Viruses and Bacteria. Aquatic
Microbial Ecology, 14, 113-118. https://doi.org/10.3354/ame014113
8. Shibata, A., Goto, Y., Saito, H., Kikuchi, T., Toda, T. and Taguchi,
S. (2006) Comparison of SYBR Green I and SYBR Gold Stains for
272 Mechatronics and Robotics
14
A Survey of Modeling and Control
of Piezoelectric Actuators
ABSTRACT
Piezoelectric actuators (PEAs) have been widely used in microand
nanopositioning applications due to their fine resolution, fast responses,
and large actuating forces. However, the existence of nonlinearities such
as hysteresis makes modeling and control of PEAs challenging. This paper
reviews the recent achievements in modeling and control of piezoelectric
actuators. Specifically, various methods for modeling linear and nonlinear
INTRODUCTION
Piezoelectric actuators (PEAs) and PEA-driven positioning systems
have been widely used in the fields of microand nanopositioning such as
atomic force microscopes [1-3], adaptive optics [4], computer components
[5], machine tools [6], aviation [7], internal combustion engines [8],
micromanipulators [9], and synchrotron-based imaging systems [10] due
to their high displacement resolution (sub nanometer) and large actuating
force (typically a few hundreds of N). PEA-driven positioning systems have
also been developed with various configurations, for example, 1-degree-
of-freedom (1-DOF) positioning systems with flexure hinge mechanisms
[11,12], stick-slip actuators [13], multiple PEAs-driven inchworms [14],
or walking actuators [15]; multi-DOF positioning systems with series
mechanism [16], parallel mechanism [17], or stick-and-clamping actuators
[18].
In all positioning control applications, the hysteresis and creep effects of
PEAs have shown to be able to significantly degrade the system performance
and even system stability [19]. For improvement, models are desirable to
represent these effects and on this basis, tracking controllers can be designed
and implemented to achieve desired positioning requirements.
This paper reviews the modeling and control methods of PEAs in
microand nanopositioning applications and some remaining issues to be
solved, and is arranged as follows.
Firstly, the working principles and behaviors of PEAs are briefly
introduced in Section 2. In Section 3 various models and modeling techniques
for PEAs are examined, along with their limitation; and in Section 4 various
control schemes for PEAs are reviewed with associated challenges. Section
5 concludes with emerging issues identified in modeling and control of
PEAs for future research.
A Survey of Modeling and Control of Piezoelectric ... 275
Figure 1. (a) Hysteresis of a PEA and (b) creep of a PEA subject to a 30 V step
input.
A Survey of Modeling and Control of Piezoelectric ... 277
Hysteresis Sub-Models
, with . Each operator has two saturation values: 0 and 1, and its
contribution to the model output is adjusted by , referred to as the
Preisach weighting function.
(1)
Figure 2. (a) Preisach hysteresis operator and (b) working principle of the Prei-
sach hysteresis model: the sum of two hysteresis operators of different switch-
ing values as shown in (1) can represent a simple but unsmooth hysteresis loop
as shown in (2); by using the sum of an infinite number of hysteresis operators
of different switching values, smooth hysteresis loops including both major and
minor loops can then be represented as shown in (3).
As such, the resulting Preisach hysteresis model is expressed as where
is the input and is the output of the hysteretic system [47]. To
reduce the computational effort due to the involvement of the double
integration, alternative expressions of the Preisach hysteresis model that
only involves arithmetic operations have also been derived from Equation
(1) and the properties of the Preisach hysteresis operator [68,69].
280 Mechatronics and Robotics
(2)
where is the model output; i denotes the index of the play operator;
and are the weight parameter and the threshold parameter of the i-th
play operator, respectively.
Similar to the PI hysteresis model, the MRC hysteresis model represents
the hysteresis by the combined effect of a finite number of elasto-slide
elements or operators as shown in Figure 4 [64]. Each elasto-slide element
or operator (Figure 4(a)(1)) is composed of 1) a mass sliding on a surface
with a Coulomb friction , where is the friction coefficient and N is
the normal force between the mass and the surface, and 2) a spring of
stiffness k with one end connected to the mass whilst the displacement of
the other end (the free end) can be freely assigned and use as input to the
elasto-slide element.
Figure 3. (a) Play operator and (b) working principle of the PI hysteresis mod-
el: the sum of two play operators of different parameters as shown in (1) form a
simple PI hysteresis model as shown in (2).
A Survey of Modeling and Control of Piezoelectric ... 281
Figure 4. (a) An elasto-slide element and its input-output relationship; (b) The
physical interpretation of the MRC hysteresis model [63].
As such, hysteresis relationship exists between the (input) displacement
of the free end of the spring and the resulting force in the spring , as
can be seen in Figure 4(a)(2).
Connecting n elasto-slide elements as in Figure 4(b), the resulting MRC
hysteresis model or the relationship between the input displacement of the
free ends of the springs and the total force experienced at the free ends
of the springs can be expressed as
(3)
where is the force in the spring of the i-th elasto-slide element
induced by .
Creep Sub-Models
(4)
where is the creep model output, is the creep output at after the
static input voltage is applied, and is a coefficient depending on the input
voltage. The dependence of on the input voltage and the unsoundness of
such a nonlinear creep model as and pose difficulties in its
implementation [19].
To avoid the problem of unboundness encountered in the nonlinear
creep model Equation (4) and to avoid the nonlinearity, a linear creep model
was proposed as.
A Survey of Modeling and Control of Piezoelectric ... 285
(5)
where the creep transfer function is represented by the combined
effects of spring damper systems [1,88]. Besides, an operator-
based linear creep model was also reported in [60].
However, these linear creep models have linear and memory free
equilibrium values, whist the actual equilibrium value of the creep effect of
PEA exhibits hysteresis nonlinearity [89].
To further solve the problems encountered in both the nonlinear and
linear creep model above, the PI hysteresis model design approach was
adopted in [63] to develop an operator-based nonlinear creep model, whose
output is a weighted sum of the values of a series of elementary creep
operators.
Model Structures
With different sub-models representing hysteresis, creep, and vibration
dynamics of PEAs defined, macroscopic PEA models are constructed by
connecting such sub models according to certain model structures. There are
two categories of model structures for macroscopic PEA models.
In the first category, the behavior of a PEA, when subject to a voltage
input and an external loading force, is decoupled into several effects such as
rate-independent hysteresis, vibration dynamics, and creep based on physical
laws. Such effects are then represented by individual sub-models and such
sub-models are connected in various manners to form a comprehensive
model of the PEA.
The most well-known model structure of PEAs belonging to this
category is the electromechanical model proposed in [64], which is further
refined in [22]. Figure 5 shows the schematic representation of this model.
A Survey of Modeling and Control of Piezoelectric ... 287
(6)
(7)
(8)
(9)
(10)
(11)
The hysteresis effect, H, and the vibration dynamics, M, can be modeled
by different kinds of sub-models. Specifically, His usually represented by
rate-independent hysteresis models, e.g. the MRC model was used in [64]
whist the first-order nonlinear differential equation model, which was firstly
proposed in [78], was employed in [22].
CONTROL OF PEAS
Over the years, a number of control schemes have been reported for positioning
control of PEAs. Due to the detrimental effects of the nonlinearities
(especially hysteresis) of PEAs on positioning control performances in
terms of positioning error and stability [19], the compensation or rejection
of such nonlinear effects has been a major consideration in all of these
control schemes. In the following, the typical control schemes for PEAs are
reviewed.
which has the same effect as the aforementioned unknown effects on the
PEA output, is subtracted from . However, the effectiveness of this
scheme was not examined through experiments. In [102], the inverse
hysteresis model is updated adaptively to account for the unknown effects.
However, the influence of the unknown effects on the PEA output, though
largely suppressed, can still be clearly observed in the experiment results.
As such, the best solution to this problem remains introducing feedback into
the control schemes, which is to be discussed in the following sections.
Feedback
Feedback control schemes as shown in Figure 8 lead to strong suppression
of the unknown effects including model errors, external loads, and changes
in the dynamics of the PEA on the position control performances, hence they
are widely used.
Specifically, in its basic form, the control signal generated by a sliding mode
controller can be divided into two parts. The first part is discontinuous,
referred to as the switching control, for compensating for the matched
uncertainties while driving the plant trajectory to a prescribed geometric
entity in the phase space (a space with each dimension corresponding to
one state variable in the plant model), referred to as the sliding surface.
The second part is continuous, referred to as the equivalent control, for
keeping the plant trajectory on the sliding surface after the latter is reached.
Since the sliding surface and hence the closed loop system dynamics after
reaching the sliding surface can be freely designed, the closedloop system
becomes insensitive to the matched uncertainties after the sliding surface is
reached [110]. The major problem of implementing the sliding mode control
techniques on PEAs is chattering and the need of the knowledge about the
bounds of the matched uncertainties for designing the discontinuous control
signal [110], which is to be discussed in more detail in Section 5.2.
Based on the model structure shown in Figure 6, hysteresis and creep
in a PEA can be treated as (matched) uncertainties in the input channels
of the vibration dynamics sub-model of the PEA (given that a linear creep
model is used). As such, a sliding mode controller can be designed using
the vibration dynamics sub-model as a nominal model of the PEA, and
the effects of hysteresis and creep on the PEA output displacement can be
strongly suppressed without knowing the exact expression of the hysteresis
sub-model and the creep sub-model, thus largely reducing the efforts
involved in hysteresis and creep modeling and compensation [95,111-114].
Various modifications of the basic sliding mode control technique have also
been reported for controlling PEAs for improved performance. For instance,
adaptive methods are combined with sliding mode control to remove the
requirement of model parameter estimation or lead to better compensation
for the nonlinearities as in [115-117]. And in [118] the nonlinearities of the
positioning mechanism driven by a PEA in additional to those of the PEA
itself were considered in sliding mode controller design.
Robust control techniques constitute another branch of widely used
advanced control techniques for positioning control of PEAs. They try to
find the control law via optimizing an objective function that incorporates
the robustness objective, for example, minimizing the H∞ norm [88,119-
122] or the H2 norm [76] of the transfer functions relating the disturbances
to the plant output, thus minimizing the effects of the disturbances on the
plant output and enhancing robustness. Robust control techniques can be
combined with sliding mode control to suppress unmatched uncertainties
A Survey of Modeling and Control of Piezoelectric ... 293
(which the sliding mode control techniques cannot compensate for), while
sliding mode control are considered to be superior in dealing with matched
uncertainties because robust control techniques can only minimize the
effects of the matched uncertainties on the plant output whilst sliding mode
control techniques can completely reject such effects in theory [110].
Other advanced control techniques such as state feedback [123], optimal
control [109,124], adaptive control [125-128], and neuron network methods
[129-131] have also been applied to control PEAs.
Inverse sub-models of the nonlinearities of PEA are sometimes inserted
between the feedback controller and the controller in Figure 8 to linearize
the PEA thus facilitating the use of linear feedback controllers. For example,
an inverse Preisach hysteresis model is used to cancel out the hysteresis
and then a linear H2 controller is designed based on the remaining linear
dynamics of a PEA in [76], while an inverse PI hysteresis model and a
sliding mode controller are used in [113].
One problem with the feedback control schemes for PEAs is that the
resulting closed-loop systems have low gain margins, which limit the use
of high gain controllers. Such low gain margin is largely due to the quick
phase loss near the first resonant peak in the frequency response of a PEA
[23]. This problem can be alleviated by using a notch filter to reduce the gain
around the first resonant peak [20,132,133].
Figure 10. A feedback controller is used for nonlinear suppression whist a feed
forward controller is the inverse model of the closed loop system for canceling
out the dynamics of the closed loop system, which can be treated as linear [19].
As such, to compensate for hysteresis and other nonlinear/uncertain
effects in the PEA, one only needs to subtract from before being
applied to the PEA. A control scheme with such a disturbance observer is
shown in Figure 11. This scheme has been shown to be effective in both step
[140] and high frequency tracking (200 Hz sinusoidal) [141] operations.
Modeling Issues
besides the hysteresis loop (Figures 3(a) and 4(a)(2)). However, to limit the
computational effort, the number of hysteresis operators used in PI or MRC
hysteresis models in practice is usually very small (typically around 10)
compared to the Preisach hysteresis model (infinity).
As a result, the PI and MRC models tend to be less accurate than the
Preisach hysteresis model in representing the hysteresis loops, especially
in the cases of the small minor loops. So, given a tractable amount of
computational effort, there is a contradiction between the capability of
representing all hysteresis phenomena, including one-sided hysteresis (the
PI and MRC hysteresis models), and the accuracy (Preisach hysteresis
models) in the existing rate-independent hysteresis models. This remains an
issue to be solved in future research.
Figure 11. Augmenting the feedback control loop with a disturbance observer.
Figure 12. Actuation sequence of a PEA-driven stick-slip actuator; (a) Start; (b)
Slow expansion and (c) Fast contraction.
Control Issues
Controller Design
As can be seen from Section 4.2, sliding mode control has been one of the most
promising techniques for PEA control applications due to their capability to
largely suppress hysteresis and creep in practice. However, there are two
remaining problems. One problem of implementing the sliding mode control
techniques on PEAs is chattering. This is a kind of high frequency vibration
in the plant output displacement induced by the discontinuous control signal
generated by a sliding mode controller, and it can induce wear, noises, and
even resonance and instability in the plant. Chattering is especially severe
in plants with fast responses, such as PEAs, because they cannot filter out
the high frequency component in the switching control signal. Another
problem is that the bounds of the matched uncertainties need to be known
for designing the discontinuous control signal, but such bounds are usually
difficult to determine in practice while conservatively setting the bounds to
large values only leads to more severe chattering [110]. Existing research
efforts have been focused on solving the chattering problem. For example,
298 Mechatronics and Robotics
State Estimation
Many advanced control techniques reviewed in Section 4, such as sliding
mode control [114], require state feedback. Since usually not all of the states
of a PEA model are measurable in application, state observers must be used.
However, the problem of estimating the states of PEA models has not been
well addressed in the literature, as discussed below.
Currently, many kinds of state observers are applicable to PEA control
and they can be classified into two categories: non-model-based and model-
based. Non-modelbased filters/differentiators, such as the low-pass filter
plus ideal differentiator, the α-β filter [145], the highgain differentiator [146],
the integral-chain differentiator [146], and the sliding mode differentiator
based on the super-twisting algorithm [147] (given that the states needed are
the derivatives of the measured plant output), usually generate large phase
lags if the desired level of noise suppression is enforced (for the first four
aforementioned methods) or excessive chattering in noisy systems (for the
last method). Compared to non-modelbased filters/differentiators, model-
based observers, e.g., the extended Kalman filter [148], the unscented
Kalman filter [149,150], and the high-gain observer [114], can generate
more accurate estimations if the imperfection of the model can be ignored.
However, with the presence of uncertainties such as the effects of hysteresis,
creep, and external load on a PEA, however, the performances of the model-
based observers degrade. For improvement, in the cases that the system
uncertainties can be treated as a lumped unknown input to the system model,
a kind of model-based observers should be used to estimate the system
states even with the presence of the unknown input, which is referred to as
unknown input observers (UIOs) in the literature. Many UIOs have been
reported in the past three decades, including the full-order UIO [151], the
reduced-order UIO [152], the UIO designed based on projection operator
[153], and the sliding mode-based observers (SMOs) [110,153-155].
A Survey of Modeling and Control of Piezoelectric ... 299
model [16].Many control schemes that have been applied to control single-
DOF PEA-driven positioning systems can also be extended to control
fully-actuated PEA-driven positioning systems, e.g. PID control [112,168],
H2 control [168], and sliding mode control [112] since such methods can be
easily extended to fully-actuated multi-in putmulti-output cases.
Compared to the fully-actuated systems, the over actuated systems
have certain advantages in the performance point of view, e.g. singularity
elimination, dexterity improvement, and better load carrying ability [169].
However, the over-actuated structure poses difficulties in both modeling and
control. Specifically, the PEAs in the mechanism need to cooperate well or
excessive internal forces will occur and may damage the mechanism, whilst the
involvement of PEA nonlinearities makes such cooperation difficult to model
and control. As a result, very few researches on the modeling and control of
over-actuated PEA-driven parallel mechanism have been reported compared
to its popularity in macro parallel mechanisms. The only example is [112], in
which the modeling and control of a xθyθz 3-DOF positioning system driven
by 4 PEAs are concerned. Specifically, in [112] each PEA in the mechanism is
modeled and controlled individually and uses only the output displacements
of the corresponding PEA as feedback, while the inverse kinematic model
of the mechanism is used to generate the desired displacement for each PEA
according to the desired position of the end-effector. The cross-axis couplings
are considered as disturbances. As such, the control problem is transformed
into several 1-DOF control problems, which are much simpler. However, in
this scheme the actual position of the end-effector is not closed-loop controlled,
and the requirement of being able to measure the displacement of each PEA
in the mechanism limits the generality of this scheme. As such, more general
schemes for modeling and controlling over-actuated multi DOF PEA-driven
positioning systems that can relax the aforementioned requirement are still to
be developed in the future.
CONCLUSIONS
Jingyang Peng, Xiongbiao Chen (chattering) and the corresponding state
estimators (the incompatibility between the observer matching condition
and the existing PEA models); finally, 6) the modeling and control of
some PEA-driven positioning systems with unconventional structures, e.g.
PEA-driven stick-slip actuators and over-actuated multiDOF PEA-driven
positioning systems, are largely untouched in the literature, thus remain to
be addressed in future research.
A Survey of Modeling and Control of Piezoelectric ... 301
REFERENCES
1. D. Croft, G. Shed and S. Devasia, “Creep, Hysteresis, and Vibration
Compensation for Piezoactuators: Atomic Force Microscopy
Application,” Journal of Dynamic Systems, Measurement, and Control,
Vol. 123, No. 1, 2001, pp. 35-43. doi:10.1115/1.1341197
2. Q. Zou, K. K. Leang, E. Sadoun, M. J. Reed and S. Devasia, “Control
Issues in High-Speed AFM for Biological Applications: Collagen
Imaging Example,” Asian Journal of Control, Vol. 6, No. 2, 2004, pp.
164-178. doi:10.1111/j.1934-6093.2004.tb00195.x
3. R. Kassies, K. O. Van der Werf, A. Lenferink, C. N. Hunter, J. D.
Olsen, V. Subramaniam and C. Otto, “Combined AFM and Confocal
Fluorescence Microscope for Applications in Bio-Nanotechnology,”
Journal of Microscopy, Vol. 217, No. 1, 2005, pp. 109-116. doi:10.1111/
j.0022-2720.2005.01428.x
4. H. Song, G. Vdovin, R. Fraanje, G. Schitter and M. Verhaegen,
“Extracting Hysteresis from Nonlinear Measurement of Wavefront-
Sensorless Adaptive Optics System,” Optics Letters, Vol. 34, No. 1,
2009, pp. 61-63. doi:10.1364/OL.34.000061
5. W. Yang, S.-Y. Lee and B.-J. You, “A Piezoelectric Actuator with
a Motion-Decoupling Amplifier for Optical Disk Drives,” Smart
Materials and Structures, Vol. 19, No. 6, 2010, pp. 065027.1-
065027.10. doi:10.1088/0964-1726/19/6/065027
6. G. Stöppler and S. Douglas, “Adaptronic Gantry Machine Tool with
Piezoelectric Actuator for Active Error Compensation of Structural
Oscillations at the Tool Centre Point,” Mechatronics, Vol. 18, No. 8,
2008, pp. 426-433. doi:10.1016/j.mechatronics.2008.03.002
7. S. R. Viswamurthy, A. K. Rao and R. Ganguli, “Dynamic Hysteresis
of Piezoceramic Stack Actuators Used in Helicopter Vibration Control:
Experiments and Simulations,” Smart Materials and Structures, Vol.
16, No. 4, 2007, pp. 1109-1119. doi:10.1088/0964-1726/16/4/020
8. M. S. Senousy, F. X. Li, D. Mumford, M. Gadala and R. K. N. D.
Rajapakse, “Thermo-Electro-Mechanical Performance of Piezoelectric
Stack Actuators for Fuel Injector Applications,” Journal of Intelligent
Material Systems and Structures, Vol. 20, No. 4, 2009, pp. 387-
399. doi:10.1177/1045389X08095030
9. J.-J. Wei, Z.-C. Qiu, J.-D. Han and Y.-C. Wang, “Experimental
Comparison Research on Active Vibration Control for Flexible
302 Mechatronics and Robotics
Smart Materials and Structures, Vol. 16, No. 6, 2007, pp. 2354-
2367. doi:10.1088/0964-1726/16/6/040
19. S. Devasia, E. Eleftheriou and S. O. R. Moheimani, “A Survey of
Control Issues in Nanopositioning,” IEEE Transactions on Control
Systems Technology, Vol. 15, No. 5, 2007. pp. 802-823. doi:10.1109/
TCST.2007.903345
20. P. Instrumente, “The World of Nanopositioning and Micropositioning
2005/2006,” Physik Instrumente, Karlsruhe, 2005.
21. I. Mayergoyz and G.Bertotti (Eds.), “The Science of Hysteresis,” Vol.
3, Elsevier, St. Louis, 2005.
22. H. J. M. T. S. Adriaens, W. L. De Koning and R. Banning, “Modeling
Piezoelectric Actuators,” IEEE/ASME Transactions on Mechatronics,
Vol. 5, No. 4, 2000, pp. 331- 341. doi:10.1109/3516.891044
23. J. A. Main and E. Garcia, “Piezoelectric Stack Actuators and Control
System Design: Strategies and Pitfalls,” Journal of Guidance, Control,
and Dynamics, Vol. 20, No. 3, 1997, pp. 479-485. doi:10.2514/2.4066
24. T. Fett and G. Thun, “Determination of Room-Temperature Tensile
Creep of PZT,” Journal of Materials Science Letters, Vol. 17, No. 22,
1998, pp. 1929-1931. doi:10.1023/A:1006608509876
25. W. T. Ang, F. A. Garmón, P. Khosla and C. N. Riviere, “Modeling
Rate-Dependent Hysteresis in Piezoelectric Actuator,” 2003 IEEE/
RSJ International Conference on Intelligent Robots and Systems, Las
Vegas, 27 October-1 November 2003, pp. 1975-1980. doi:10.1109/
IROS.2003.1248937
26. W. T. Ang, P. K. Khosla and C. N. Riviere, “Feedforward Controller
with Inverse Rate-Dependent Model for Piezoelectric Actuators
in Trajectory-Tracking Applications,” IEEE/ASME Transactions
on Mechatronics, Vol. 12, No. 2, 2007, pp. 134-142. doi:10.1109/
TMECH.2007.892824
27. R. S. Robinson, “Interactive Computer Correction of Piezoelectric
Creep in Scanning Tunneling Microscopy Images,” Journal of
Computer-Assisted Microscopy, Vol. 2, No. 1, 1996, pp. 53-58.
28. X. B. Chen, Q. Zhang, D. Kang and W. Zhang, “On the Dynamics of
Piezoactuated Positioning Systems,” Review of Scientific Instruments,
Vol. 79, No. 11, 2008, pp. 116101.1-116101.3. doi:10.1063/1.2982238
29. W. Chen and C. S. Lynch, “A Micro-Electro-Mechanical Model for
Polarization Switching of Ferroelectric Materials,” Acta Materialia, Vol.
304 Mechatronics and Robotics
212-217. doi:10.1109/20.43892
49. P. Hejda and T. Zelinka, “Generalized Preisach Model of Hysteresis—
Theory and Experiment,” Czechoslovak Journal of Physics, Vol. 40,
No. 1, 1990, pp. 57-68. doi:10.1007/BF01598355
50. M.-J. Jang, C.-L. Chen and J.-R. Lee, “Modeling and Control of a
Piezoelectric Actuator Driven System with Asymmetric Hysteresis,”
Journal of the Franklin Institute, Vol. 346, No. 1, 2009, pp. 17-
32. doi:10.1016/j.jfranklin.2008.06.005
51. P. N. Sreeram, G. Salvady and N. G. Naganathan, “Hysteresis
Prediction for a Piezoceramic Material System,” The ASME Winter
Annual Meeting, New Orleans, 28 November-3 December 3 1993, pp.
35-42.
52. D. C. Hughes and J. T. Wen, “Preisach Modeling and Compensation
for Smart Material Hysteresis,” Proceedings of SPIE, Vol. 2427, 1995,
pp. 50-64.
53. P. Ge and M. Jouaneh, “Generalized Preisach Model for Hysteresis
Nonlinearity of Piezoceramic Actuators,” Precision Engineering, Vol.
20, No. 2, 1997, pp. 99-111. do:10.1016/S0141-6359(97)00014-7
54. H. Hu and R. B. Mrad, “On the Classical Preisach Model for Hysteresis
in Piezoceramic Actuators,” Mechatronics, Vol. 13, No. 2, 2002, pp.
85-94. doi:10.1016/S0957-4158(01)00043-5
55. G. Robert, D. Damjanovic and N. Setter, “Preisach Modeling of
Piezoelectric Nonlinearity in Ferroelectric Ceramics,” Journal ofApplied
Physics, Vol. 89, No. 9, 2001, pp. 5067-5074. doi:10.1063/1.1359166
56. G. Song, J. Zhao, X. Zhou and J. A. De Abreu-Garcia, “Tracking Control
of a Piezoceramic Actuator with Hysteresis Compensation Using
Inverse Preisach Model,” IEEE/ASME Transactions on Mechatronics,
Vol. 10, No. 2, 2005, pp. 198-209. doi:10.1109/TMECH.2005.844708
57. X. Zhao and Y. Tan, “Neural Network based Identification of Preisach-
Type Hysteresis in Piezoelectric Actuator Using Hysteretic Operator,”
Sensors and Actuators A: Physical, Vol. 126, No. 2, 2006, pp. 306-
311. doi:10.1016/j.sna.2005.10.023
58. X. Zhao and Y. Tan, “Modeling Hysteresis in Piezo Actuator based on
Neural Networks,” Advances in Computation and Intelligence, Lecture
Notes in Computer Science, Vol. 5370, 2008, pp. 290-296.
59. K. Kuhnen and H. Janocha, “Adaptive Inverse Control of Piezoelectric
Actuators with Hysteresis Operators,” Proceedings of European
A Survey of Modeling and Control of Piezoelectric ... 307
897-919. doi:10.1016/0020-7225(86)90023-6
79. R. Bouc, “Forced Vibrations of a Mechanical System with Hysteresis,”
Proceedings of the 4th Conference on Non-linear Oscillations, Prague,
5-9 September 1967, pp. 315.
80. Y.-K. Wen, “Method of Random Vibration of Hysteretic Systems,”
Journal of the Engineering Mechanics Division, Vol. 102, No. 2, 1976,
pp. 249-263.
81. M. Ismail, F. Ikhouane and J. Rodellar, “The Hysteresis Bouc-Wen
Model, a Survey,” Archives of Computational Methods in Engineering,
Vol. 16, No. 2, 2009, pp. 161- 188. doi:10.1007/s11831-009-9031-8
82. O. Gomis-Bellmunt, F. Ikhouane and D. MontesinosMiracle, “Control
of Bouc-Wen Hysteretic Dystems: Application to a Piezoelectric
Actuator,” The 13th Power Electronics and Motion Control Conference,
2008. EPEPEMC 2008, Poznan, 1-3 September 2008, pp. 1670-
1675. doi:10.1109/EPEPEMC.2008.4635507
83. M. Rakotondrabe, “Bouc-Wen Modeling and Inverse Multiplicative
Structure to Compensate Hysteresis Nonlinearity in Piezoelectric
Actuators,” IEEE Transactions on Automation Science and Engineering,
Vol. 8, No. 2, 2011, pp. 428-431. doi:0.1109/TASE.2010.2081979
84. D. H. Wang, W. Zhu and Q. Yang, “Linearization of Stack Piezoelectric
Ceramic Actuators Based on BoucWen Model,” Journal of Intelligent
Material Systems and Structures, Vol. 22, No. 5, 2011, pp. 401-
413. doi:10.1177/1045389X10386132
85. L. Deng and Y. Tan, “Modeling Hysteresis in Piezoelectric Actuators
Using NARMAX Models,” Sensors and Actuators A: Physical, Vol.
149, No. 1, 2009, pp. 106- 112. doi:10.1016/j.sna.2008.09.022
86. A. Badel, J. Qiu and T. Nakano, “A New Simple Asymmetric Hysteresis
Operator and Its Application to Inverse Control of Piezoelectric
Actuators,” IEEE Transactions on Ultrasonics, Ferroelectrics and
Frequency Control, Vol. 55, No. 5, 2008, pp. 1086-1094. doi:10.1109/
TUFFC.2008.761
87. H. Jung, J. Y. Shim and D. Gweon, “New Open-Loop Actuating
Method of Piezoelectric Actuators for Removing Hysteresis and
Creep,” Review of Scientific Instruments, Vol. 71, No. 9, 2000, pp.
3436-3440. doi:10.1063/1.1287627
88. K. Kuhnen and H. Janocha, “Compensation of Creep and Hysteresis
Effects of Piezoelectric Actuators with Inverse Systems,” 6th
310 Mechatronics and Robotics
025105.11. doi:10.1063/1.2841804
168. T. W. Seo, D. S. Kang and J. Kim, “Synthesis and Comparison of
Fine Actuator Controllers for a 3-DOF Micro Parallel Positioning
Platform,” Proceedings of the 10th International Conference on
Control, Automation, Robotics and Vision, Hanoi, 17-20 December
2008, pp. 1606- 1611. doi:10.1109/ICARCV.2008.4795765
169. M. Grotjahn, B. Heimann and H. Abdellatif, “Identification of Friction
and Rigid-Body Dynamics of Parallel Kinematic Structures for Model-
Based Control,” Multibody System Dynamics, Vol. 11, No. 3, 2004,
pp. 273-294. doi:10.1023/B:MUBO.0000029426.05860.c2
CHAPTER
15
Experimental Investigation of
Photostrictive Materials for MEMS
Application
ABSTRACT
Micro-Electro-Mechanical-Systems (MEMS) is an emerging technology
that integrates micromachined mechanical structures with integrated circuits
(IC). During the last two decades, it has grown vastly as a prominent research
area among scientists and engineers. Usually the main components of MEMS
INTRODUCTION
Photostriction is a phenomenon in which strain is induced in photostrictive
materials by incident light. In principle, this effect arises from a superposition
of the photovoltaic effect, i.e., the generation of large voltage from the
irradiation of light, and the converse piezoelectric effect, i.e., expansion
or contraction under the voltage applied. A “bulk photovoltaic effect” is
noticeable in some ferroelectric (or non-centrosymmetric) materials, such
as ferroelectric single crystals or polarized ferroelectric ceramics. Along
with this photovoltage, the converse piezoelectric effect induces mechanical
strain [1]. The combinations of these two effects make these ferroelectric
materials candidates for wireless actuator applications, that is, actuators
excited by incident light. Photostrictive ceramics are one promising material
for micro-fabrication technologies, which have photovoltaic, piezoelectric,
Experimental Investigation of Photostrictive Materials for ... 321
Research Hypotheses
Transverse deflection of the PLZT optical actuator depends on the light
intensity. The deflection of the PLZT optical actuator will vary for the
continuous and pulses of lights. Fast response and significant deflection of
the PLZT optical actuator under the illumination of light could be used in
MEMS applications as a wireless control.
Research Objectives
Photostrictive actuators have the superiority of remote sensing capabilities and
freedom from external electric/magnetic disturbances, unlike conventional
actuators. Based on this inspiration from current research, a PLZT thin
film placed on the top of a silicon wafer is used to support a systematic
characterization of photostrictive properties to develop an actuator. The
main goal of the research is to investigate the potential of the PLZT material
as a wireless actuator to convert photonic energy to mechanical energy.
The objectives of the present research follow.
1) Investigate the photostrictive effect of a continuous high intensity
near-ultraviolet (UV) light on a PLZT optical actuator cantilever
beam, which includes:
• measuring transverse deflection of the photostrictive cantilever
beam for UV light focused at three different locations on the beam;
• for each focused location on the cantilever beam, measuring
the transverse deflection along the beam at three different light
intensities: 250 mW/cm2, 200 mW/cm2 and 100 mW/cm2; and
Experimental Investigation of Photostrictive Materials for ... 323
(1)
where Es is the saturated photovoltaic field; as = La/ba (length/width)
is the aspect ratio; α is optical actuator constant and β is a constant related
to voltage leakage. The light intensity I(t) is a function of time, and El(tj)
denotes the induced electric field at time tj. While I(t) remains constant
during the interval Δt, then the electric field at tj=tj−1+Δttj=tj−1+Δt is El(tj).
The actuator temperature can be defined in a difference equation [27]:
(2)
where P is the absorbed heat from the light illumination; H is the heat
capacity of the actuator; g is the heat transfer rate from the actuator to the
structure; θ(tj) is the actuator body temperature at tj; I(tj) is an illumination
intensity at tj; and Δθ is a temperature rise in a time interval Δt.
The temperature change can trigger the pyroelectric effect from which
an additional voltage is generated. The electric field Eq(t), contributed by
pyroelectric effect, can be determined from [27] :
Experimental Investigation of Photostrictive Materials for ... 327
(3)
where Pn is the pyroelectric constant.
The total induced electric field E(t) including both the photovoltaic
effect and the pyroelectric effect is [27] :
(4)
The magnitude of the light induced strains can be written as [18] :
(5)
where d33 is the piezoelectric-strain constant, λ is a thermal stress
coefficient, and Ya is the Young’s modulus of the optical actuator.
The electric potential of the actuator is
(6)
The objective of the current research is to investigate the influence of
continuous and pulses of light on the transverse deflection of PLZT optical
actuator. How the transverse deflection of a cantilever beam depends
on the light intensity (both continuous and pulses) will be investigated
experimentally. If the intensity of the light and the focused locations of
the beam affect the time response will also be investigated in the current
research.
METHODOLOGY
Experimentation has been performed mainly using a high intensity light
source, a good repeatability displacement measurement sensor, and a
photostrictive optical actuator.
The goal of the experiment was to observe the transverse deflection of a
PLZT thin film on a silicon wafer cantilever beam under the illumination of
continuous and pulses of high intensity light.
Experiment Method
The set-up consists of a high pressure short arc xenon lamp with lamp
housing, power supply with igniter, hot mirror, band pass filter, optical
chopper, photostrictive cantilever type optical actuator, and laser sensor
head with controller. Based on the design, an actual experimental set-up was
developed as shown in Figure 2.
A sample of PLZT (3/52/48) ceramic was doped with 0.5% WO3 of 1 µm
thin film deposited on top of a 550 µm silicon wafer as shown in Figure 3.
PLZT material was deposited on the silicon wafer at the MEMS Laboratory
of Georgia Institute of Technology using the spattering method. PLZT
(3/52/48) means Pb concentration of 97%, La concentration of 3% with Zr/
Ti ratio 52/48. (3/52/48) ceramic also can be represented as Pb0.97 La 0.03
(Zr0.52Ti 0.48) 0.99 O3.
A 150 watt xenon short arc lamp was provided high intensity light for
maximum photovoltaic effect on the optical actuator. The arc lamp was
placed inside the Science Tech 201 - 100 series lamp housing. A model 500
- 200 series AC power supply and a model 500-IG igniter were used to turn
on the short arc lamp. To produce pulses of high intensity light on an optical
actuator the Scitech 300CD model an optical chopper was used. A Keyence
LK-H087 laser sensor head with 0.1 µm repeatability was used to measure
the transverse deflection of actuator beam.
the experiment. PLZT thin film on silicon wafer sample was mounted as
a cantilever beam on a stand. An LK-H087 laser displacement sensor was
mounted on a linear actuator at the other side of the light source. The linear
actuator can travel a minimum of 0.1 mm/step in the horizontal direction,
thus allowing the laser displacement sensor head to move along the length
of the sample cantilever beam. To irradiate pulses of light on the sample
a variable frequency optical chopper was placed between the sample and
the lamp housing window lens. Source light beam has a maximum light
intensity of 250 mW/cm2 at its focal point through the lamp housing window
lens at 20 cm distance, which was measured by a UV meter. The UV meter
sensor head was 12 mm in diameter, which measured the intensity on the
average of sensor head area.
It was not experimentally possible for the whole cantilever beam to be
exposed with equal light intensity. The light intensity starts to decrease as
far as from the center of focal point on its perpendicular plane. Experiments
were performed by focusing the light on the targeted cantilever beam at
three different locations; these are the free end, middle, and fixed end as
shown in Figure 4. Using the UV meter the light intensities were measured,
at a distance of 20 cm, 25 cm, and 35 cm from the lamp housing window
lens, were 250 mW/cm2, 200 mW/cm2, and 100 mW/cm2 respectively.
Experimental Procedure
In conducting the experiment to produce photostriction effect and measuring
the transverse deflection of the actuator, the following steps have been
followed:
• A high-pressure short arc xenon lamp was placed properly inside
the lamp housing with right polarity of cathode (−) and anode (+).
This was used as a UV light source for the experiment.
• A power supply and a lamp igniter were connected to the lamp
housing; and the lamp housing was placed on a vibration isolation
table.
• A hot mirror lens and a 428 nm wavelength filter were placed after
the lamp inside the housing filter slot. The source light first passed
through the hot mirror then through the 428 nm wavelength filter.
• Output light was focused on the PLZT optical actuator cantilever
beam at different distances from the output lens to get various
light intensities on the targeted area.
Experimental Investigation of Photostrictive Materials for ... 331
Figure 4. UV light was focused (a) At the fixed end; (b) At the middle; and (c)
At the free end of the cantilever beam.
RESULTS
Figure 6. Controller output data collected in every 0.1 sec interval for 480 sec.
Experimental Investigation of Photostrictive Materials for ... 333
Figure 10. Transverse deflection at the free end of the beam with three differ-
ent light intensities for three different focused locations on the cantilever beam.
When the gap between the light source and the targeted area of the beam
increased the light beam diffused and allowed the cantilever beam to expose
more light to the actuator surface; this happened especially when the light
was focused at the middle position compare to the other two positions. At
100 mW/cm2, the deflections were the same when the focused position of the
light was at the middle and at the fixed end of the cantilever beam as shown
in Figure 13. More diffused light fell on the beam surface when the light
was focused at the middle of the beam compared to the other two positions.
In addition, the light focused at the fixed end position allowed the beam to
bend more near the fixed end, which ultimately resulted in these having the
same beam deflection.
Experimental Investigation of Photostrictive Materials for ... 335
Figure 11. Variation of transverse deflection along the length of the cantilever
beam for three different light focus positions on the beam at 250 mW/cm2.
Figure 12. Variation of transverse deflection along the length of the cantilever
beam for three different light focus positions on the beam at 200 mW/cm2.
Figure 13. Variation of transverse deflection along the length of the cantilever
beam for three different light focus positions on the beam at 100 mW/cm2.
336 Mechatronics and Robotics
Figure 14. Comparison of time response on transverse deflection for three dif-
ferent light intensities when light was focused at the fixed end.
Experimental Investigation of Photostrictive Materials for ... 337
Figure 15. Comparison of time response on transverse deflection for three dif-
ferent light intensities when light was focused at the middle.
Figure 16. Comparison of time response on transverse deflection for three dif-
ferent light intensities when light was focused at the free end.
Figure 17. Hysteresis of the cantilever beam deflection at the free end.
338 Mechatronics and Robotics
Transverse deflection increases along the length of the actuator beam with
the increase of distance from the fixed end as shown in Figure 18. It is
also been observed that the deflection of the cantilever beam increases with
increase of the light intensity as shown in Figure 19.
Figure 18. Transverse deflection along the length of the cantilever beam at
three different light intensities for pulses of light condition.
Figure 19. Transverse deflection depends on the light intensity in pulses of light
condition.
Rising time was about 20 sec to reach 63% of the steady state maximum
deflection. The deflection of the cantilever beam was approaching a steady
state after about 180 sec.
Experimental Investigation of Photostrictive Materials for ... 339
Figure 20. Comparison of time response on transverse deflection for three dif-
ferent light intensities in the pulses of light condition.
CONCLUSIONS
From this experimental investigation, it has been observed that the continuous
lights and the pulses of lights have a significant influence on the transverse
deflection of PLZT optical actuator; less deflection was observed for the
pulses of light compare to the continuous light. The PLZT optical actuator
could be used in MEMS applications via a wireless control but the time
response of the actuator observed was slow.
Experimental Investigation of Photostrictive Materials for ... 341
From the current study, analysis, and results of this research, the following
conclusions can be made:
1) Experimental transverse deflection of a cantilever beam increases
with the increase of light intensity for both continuous lights and
pulses of lights.
2) Maximum tip deflection at the free end of the cantilever beam for
pulses of light was about 1/3rd deflection of that for continuous
light.
3) The light intensity and focused locations of the cantilever beam
have no major influence on the time response. Rising time was
observed to be about 20 sec to reach 63% of the steady state the
maximum deflection for both continuous lights and pulses of light
conditions. Deflection of the cantilever beam was approaching a
steady state after about 180 sec for both continuous lights and
pulses of lights.
342 Mechatronics and Robotics
REFERENCES
1. Fridkin, V.M. (1979) Photo Ferroelectrics. Springer Series in Solid
State Sciences, Vol. 9, 85-113. https://doi.org/10.1007/978-3-642-
81351-1_6
2. Rahman, M. (2005) A Finite Element Approach to Model and Analyze
Photostrictive Optical Actuators. PhD Dissertation, University of
Alabama, Tuscaloosa.
3. Kawaguchi, K., ichiki, M., Morikawa, Y. and Nakada, T. (2002)
Electrical Properties of Photovoltaic PLZT and Its Application to
Optical Motor. Ferroelectric, 273, 41-46.
4. Tanimura, M. and Uchino, K. (1988) Effect of Impurity Doping on
Photostriction in Ferroelectric Ceramics. Sensors and Materials, 1, 47-
56.
5. Sada, T., Inoue, M. and Uchino, K. (1987) Photostriction in PLZT
Ceramics. Journal of the Ceramic Society of Japan, International
Edition, 95, 499-504.
6. Okuym, M. (1998) Micro Sensors and Micro Actuators using
Ferroelectric Thin Films. Proceedings of the International Symposium
on Micromechatronics and Human Science, Nagoya, 29-34.
7. Satol, K., Sugama, A., Akahoshi, T., Kondo, M., Aoki, S. and Kurihara,
K. (2007) Electro-Optic Effect of Lanthanum-Modified Lead Zirconate
Titanate and Application to Optical Device. 16th IEEE International
Symposium on Applications of Ferroelectrics, 827-828.
8. Huff, M.A. (2012) What Is Mems Technology? MEMS and
Nanotechnology Exchange. https://www.mems-exchange.org/MEMS/
what-is.html
9. Chynoweth, A.G. (1956) Surface Space-Charge Layers in Barium
Titanate. Physics Review, 102, 705-714. https://doi.org/10.1103/
PhysRev.102.705
10. Brody, P.S. (1975) High Voltage Photovoltaic Effect in Barium Titanate and
Lead Titanate-Lead Zirconate Ceramics. Solid State Chemistry, 12, 193-200.
11. Brody, P.S. (1983) Optomechnical Bimorph Actuator. Ferroelectrics,
50, 27-33. https://doi.org/10.1080/00150198308014428
12. Qin, M., Kui, Y. and Yung, C.L. (2007) Photo Induced Current in
(Pb0.97La0.03) (Zr0.52Ti0.48)O3 Thin Films of Different Thicknesses.
Integrated Ferroelectrics, 88, 58-67. https://doi.org/10.1080/10584580
601098647
Experimental Investigation of Photostrictive Materials for ... 343
16
Theory and Simulation Analysis of the
Mode Shape and Normal Shape
Actuators and Sensors
ABSTRACT
In this paper, we will try to find a universal theoretical model and approximate
solutions which can be applied to both mode shape and normal shape actuators
and sensors, and which can be predicted the gain of the first three modes of
the mode shape and normal shape actuators and sensors, finally through
computer simulation analysis to validate. In order to prove the feasibility
of the theory and as well as convenient to use on the electro-mechanical
Citation:- Jou, J. (2014), “Theory and Simulation Analysis of the Mode Shape and
Normal Shape Actuators and Sensors”. Open Journal of Acoustics, 4, 184-203. doi:
10.4236/oja.2014.44019.
Copyright: © 2014 by authors and Scientific Research Publishing Inc. This work is li-
censed under the Creative Commons Attribution International License (CC BY). http://
creativecommons.org/licenses/by/4.0
346 Mechatronics and Robotics
INTRODUCTION
Due to the piezoelectric materials (such as PZT or PVDF) have the
advantages of fast response and high actuating force, so they are very
suitable made for sensors and actuators to sense and control the vibration
of the flexible structures, especially to sense or control the low frequency
vibration, because the high-frequency vibration is very easy to be absorbed
by the structure itself. So that the most piezoelectric sensors and actuators
are used to sense and control the low frequency vibration of flexible
structures. And since the design concept and theory of the modal sensors
and actuators have been proposed by C.K., Lee [1] [2] , making them to
sense and control the low frequency vibration of flexible structures are more
rapid, sensitive, precise, accurate and efficient. In which the main design
concepts and methods of the modal sensors and actuators are based on the
coating of the effective electrode surface and changing of the polarization
direction with specific modal or mode shape, shown as Figure 1. As for the
theoretical model, it is entirely based on the mode shape functions or spatial
distribution methods. It is successfully pulled out the transverse vibration
displacement of the coupling multi-modal by the mode shape functions
or spatial distribution methods. Next, it uses the coating methods or the
functions of effective electrode surface to sense or control the low frequency
vibration of the flexible structures. Since the relevant theoretical analysis,
computer simulation and experimental eleven were uncovered or presented
[3] - [9]. In addition, the theoretical model of modal or mode shape actuator
is also applied to the different types of piezoelectric transformers, thereby to
enhance their step-up ratio and conversion efficiency [10] [11] .
However, the theoretical model is only applicable to modal sensors and
actuators, but it can’t be applied on the normal shape sensors and actuators.
If the theoretical model of modal sensors and actuators used in the normal
shape sensors and actuators, we could not find the correct approximate
Theory and Simulation Analysis of the Mode Shape and ... 347
solution. That is, so far, we still can’t find one kind of universal equation
to simultaneously describe or predict the behavioral of the mode shape and
normal shape sensors and actuators, or to compare the gain or to distinguish
the difference between the mode shape and normal shape sensors and
actuators from past studies.
Figure 1. (a) The side view of the polarization direction of the bimorph piezo-
electric cantilever beam; (b) The top view of the effective electrode surface of
mode shape 1; (c) The side view of the mode shape 1; (d) The top view of the ef-
fective electrode surface of mode shape 2; (e) The side view of the mode shape
2; (f) The top view of the effective electrode surface of mode shaspe 3; (g) The
side view of the mode shape 3.
In this paper, we will try from the original actuator and sensor equation,
and electromechanical boundary conditions to find one kind of universal
theoretical model and reasonable approximate solutions to compare the
gain and difference between the mode shape and normal shape sensors and
actuators. Furthermore we will also through the analysis results of computer
simulation to verify the correctness and feasibility of the universal theoretical
model and approximate solutions.
ACTUATOR EQUATION
In order to understand the differences between the mode shape and
normal shape actuators and sensors, we just start from the one-dimension
actuator equation of bimorph type piezoelectric cantilever beam. The
main contribution of this equation comes from the bending effect, so the
membrane effect can be completely ignored. Therefore, the equation can be
simplified to as follows [12]:
348 Mechatronics and Robotics
(1)
where the moment M1 can be divided into two terms of the mechanical and
electrical moment M 1 E as follow:
(2)
As for the constants in the above Equations (1) and (2) can be defined
as follow:
(3)
(4)
(5)
wherein above the symbols of c11, D11, h, w, t, x, z and are defined as the Young’s
modulus, bending stiffness constant per unit width, thickness, transverse displacement,
time, longitudinal coordinate, transverse coordinate and density.
Furthermore the electrical moment of Equation (2) can be divided into two
terms of time domain and space domain as follow:
(6)
And
(7)
wherein above the symbols of are defined as the
function of effective electrode surface of the nth mode shape, length, natural angular
velocity (natural frequency), mode shape eigenvalues, mode shape constant and mode
shape function.
Let Equations (2)-(7) substituted into Equation (1), and after finishing,
we can get a non-homogeneous partial differential equation of motion with
transverse displacement of mode shape actuator as follows:
Theory and Simulation Analysis of the Mode Shape and ... 349
(8)
When the surface is uniformly coated and normal polarization with the
bimorph type piezoelectric cantilever beam, that is, the function of effective
electrode surface of Equation (7) can be redefined as:
(9)
Then, the second derivative of electrical moment becomes zero.
(10)
So for the normal shape actuator, Equation (8) can be simplified as:
(11)
In order for the mode shape and normal shape actuator equations can be
applied to structures of different lengths, we can try to make Equation (8)
and Equation (11) become the dimensionless equations:
(12)
And
(13)
where the dimensionless mode shape constant can be redefined as:
(14)
where the solution form of Equation (12) and Equation (13) can be
expressed as:
(15)
And
350 Mechatronics and Robotics
(16)
Furthermore let the homogeneous solution of Equation (16) and Equation
(17) are assumed to be as follows:
(17)
As for the particular solution of Equation (15) is assumed to be as
follows:
(18)
where the symbols of are defined as the undetermined coefficient and
dimensionless homogeneous transverse displacement.
We can take Equation (17) into Equation (12) and Equation (13), then
after finishing, we can get a dimensionless homogeneous four-order ordinary
differential equation, as:
(19)
where the dimensionless natural eigenvalue of mode shape is defined as:
(20)
And the i natural resonance frequency of piezoelectric beam
th
(21)
We can further get a dimensionless homogeneous solution from Equation
(19)
(22)
In addition, we can take Equation (18) into Equation (12), so get the
undetermined coefficient as:
Theory and Simulation Analysis of the Mode Shape and ... 351
(23)
where the relationship of the fourth-order differential function and zero-
order function is:
(24)
So Equation (23) can be rewritten as:
(25)
Thus we can get the solution of the undetermined coefficient.
(26)
where in normal or natural state.
At this point, we can find the general solution of mode shape and normal
shape actuator, as:
(27)
And
(28)
Next, we can find the special solutions of the bimorph type piezoelectric
cantilever beam from the electro- mechanical boundary conditions as follow:
(29)
And
(30)
Or
(31)
352 Mechatronics and Robotics
And
(32)
where the dimensionless mode shape function of mode shape actuator
under clamped-free boundary conditions as:
(33)
(34)
And
(35)
According to the electro-mechanical and clamped-free boundary conditions, we
can determine the constants of Equation (27) and Equation (28), as follow:
(36)
And
(37)
where the constants of mode shape actuator are defined as
(38)
And
(39)
where the second-order and third-order derivative mode shape function on the free
boundary as follow:
(40)
Theory and Simulation Analysis of the Mode Shape and ... 353
And
(41)
As for the constants of normal shape actuator were defined as
(42)
And
(43)
which let for the function of the effective electrode surface of normal
shape actuator.
Until now, we can get a dimensionless general solution of the mode
shape and normal shape actuator as follow:
(44)
And
(45)
Due to the high frequency vibrations can easily be absorbed by the
structure itself, so we are only interested in low-frequency vibration. And
in order to understand the difference between the mode shape and normal
shape actuators, we analyzed only for the first three resonant modes of the
structures. And according to Equation (34) and Equation (35), we can get
the eigenvalues and parameters of the first three modes, shown as Table 1.
As the same time, we can get the dimensionless mode shape functions of the
first three resonant modes by Equation (33) and Table 1, shown as Figure 2.
And we can also get the dimensionless function of the effective electrode
surface of first three modes when let Equation (33) take into Equation
(7), shown as Figure 3. In addition, we can get a dimensionless transverse
displacement of the mode shape and normal shape piezoelectric stator under
354 Mechatronics and Robotics
Figure 2. The dimensionless mode shape function of the first three modes of
mode shape actuator relative to the unit length of structure when i = n.
Table 1. The eigenvalues and parameters of the first three modes of mode shape
actuator under clamped-free boundary conditions
Another according to Equation (40) and Equation (41), we can get the second and
third derivative of mode shape function of the first three modes of mode shape actuator
relative to the unit length of structure, shown as Figure 4 & Figure 5. And we can further
that of the first three modes of mode shape actuator at free end.
Therefore in, we can let Equation (44) further simplified as follows:
Theory and Simulation Analysis of the Mode Shape and ... 355
(46)
So far, Equation (46) of mode shape actuator appeared to be consistent
with the particular solution of previously papers [1] [10] [11] .
Furthermore, we can make Equation (45) of normal shape actuator
simplifies to
(47)
where the dimensionless mode shape function of normal shape actuator
is defined as
(48)
And the dimensionless parameters of normal shape actuator is defined as
(49)
According to Equation (46) and Equation (47), we can get the ratio of
dimensionless transverse displacement
Figure 4. The second derivative of mode shape function of the first three modes
of mode shape actuator relative to the unit length of structure when i = n.
356 Mechatronics and Robotics
Figure 5. The third derivative of mode shape function of the first three modes of
mode shape actuator relative to the unit length of structure when i = n.
of mode shape and normal shape actuator under conditions of steady
state, constant driving voltage and the same bending stiffness constant per
unit width as follows:
(50)
In order to facilitate understanding of the difference between mode shape
and normal shape actuator, we can set the above ratio at the free end, that is,
, shown as Figure 2. There fore Equation (50) can
again be expressed as
(51)
SENSOR EQUATION
For the bimorph type piezoelectric cantilever beam under condition of
constant electric potential, the sensor or current equation per unit length and
width of the mode shape and normal shape sensor can be expressed as [12] :
(52)
And we can further take Equation (46) and Equation (47) into Equation
(52) respectively, let Equation (52) be divided into two types of sensor or
current equations as follow:
Theory and Simulation Analysis of the Mode Shape and ... 357
(53)
And
(54)
where the first derivative of dimensionless mode shape function of the
first three modes can be obtained from clamped-free boundary conditions,
or known from Figure 6. Therefore the ratio of dimensionless current of
the mode shape and normal shape sensor dependent of eigenvalues under
conditions of clamped-free, steady state, the same driving voltage and
bending stiffness constant per unit width as follows:
(55)
where the first derivative of dimensionless mode shape function of the
mode shape and normal shape sensor can be expressed as follow:
Figure 6. The first derivative of dimensionless mode shape function of the first
three modes relative to the unit length of structure when i = n.
358 Mechatronics and Robotics
(56)
And
(57)
Since the current is proportional to voltage under condition of the same
load RL, So Equation (55) can also be expressed the ratio of dimensionless
voltage of the mode shape and normal shape sensor as follows:
(58)
Figure 7. The schematic diagram of the mode shape and normal shape actuators
and sensors.
360 Mechatronics and Robotics
Table 2. The size and physical properties of one-dimension bimorph type piezo-
electric cantilever beam.
Figure 8. Modeling and setting physical properties of the mode shape and nor-
mal shape actuators and sensors. (a) The normal shape actuator and sensor; (b)
The mode shape 1 actuator and sensor; (c) The mode shape 2 actuator and sen-
sor; (d) The mode shape 3 actuator and sensor.
Theory and Simulation Analysis of the Mode Shape and ... 361
Figure 9. Meshing the mode shape and normal shape actuators and sensors. (a)
The normal shape actuator and sensor; (b) The mode shape 1 actuator and sen-
sor; (c) The mode shape 2 actuator and sensor; (d) The mode shape 3 actuator
and sensor.
362 Mechatronics and Robotics
CONCLUSION
According to the results of the theory and simulation analysis, on the whole,
we found the first modal of the normal shape and mode shape actuators and
sensors are the same gain.
sensor; (b) The mode shape 1 actuator and sensor; (c) The mode shape 2 actuator and
sensor; (d) The mode shape 3 actuator and sensor.
Table 6. The gain or ratio of dimensionless transverse displacement and voltage of the
mode shape and normal shape actuators and sensors relative to the frequency under
condition of different frequency spacing (Δfi = 1 Hz, 10 Hz & 100 Hz).
F r e q u e n c y Actuators Sensors
Spacing MS1A. MS2A. M S 3 A . MS1S. MS2S. MS3S.
N S 1 A NS2A N S 3 A N S 1 S NS2S NS3S
(m/m) (m/m) (m/m) (V/V) (V/V) (V/V)
1 Hz 0.87 1.16 2.40 0.64 1.19 2.41
10 Hz 0.87 1.12 2.03 0.64 1.15 2.03
100 Hz 1.08 1.20 1.94 0.72 1.20 1.94
364 Mechatronics and Robotics
Figure 11. Post-processing the first three modes and the maximum deformation
of the mode shape and normal shape actuators and sensors. (a) The first three
modes of the normal shape actuator and sensor; (b) The first three modes of the
mode shape 1 actuator and sensor; (c) The first three modes of the mode shape 2
actuator and sensor; (d) The first three modes of the mode shape 3 actuator and
sensor (analysis type: modal).
Or in other words, the design concept of the first modal of mode shape
actuator or sensor is not necessarily better than of the first modal of normal
shape actuator or sensor, or even worse. However, the gain of the second
and third modal of the mode shape actuators and sensors are better than the
normal shape actuators and sensors under any operating states.
Figure 12. Post-processing the first three modes and the maximum electric po-
tential of the mode shape and normal shape actuators and sensors. (a) The first
Theory and Simulation Analysis of the Mode Shape and ... 365
three modes of the normal shape actuator and sensor; (b) The first three modes
of the mode shape 1 actuator and sensor; (c) The first three modes of the mode
shape 2 actuator and sensor; (d) The first three modes of the mode shape 3 ac-
tuator and sensor (analysis type: modal).
Most importantly, we have to find a universal theoretical model and
approximate solutions in this paper, which can predict the gain of the first
three modes of normal shape and mode shape actuators and sensors, and
through the analysis results of computer simulation to confirm.
Figure 13. Post-processing the electric potential output state of the mode shape
and normal shape actuators and sensors. (a) The normal shape actuator and sen-
sor; (b) The mode shape 1 actuator and sensor; (c) The mode shape 2 actuator
and sensor; (d) The mode shape 3 actuator and sensor (analysis type: steady
state).
Table 7. The resonance frequency and the ratio of resonance frequency of the
first three modes of the mode shape and normal shape actuators and sensors by
theory and simulation analysis (analysis type: modal).
Table 8. The maximum deformation and the ratio of maximum deformation of the
mode shape and normal shape actuators by simulation analysis (analysis type: static &
steady state).
Anal- NSA (m) MS1A (m) MS2A (m) MS3A (m) MS1A. MS2A. MS3A.
ysis NSA NSA NSA
Type (m/m) (m/m) (m/m)
Static 5.18E−14 −1.52E−13 −9.57E−11 −1.74E−10 2.93 1846.88 3353.80
&
Steady
State
Table 9. The maximum sensing or output voltage and the ratio of maximum
sensing or output voltage of the first three modes of the mode shape and normal
shape sensors by simulation analysis (analysis type: modal).
ACKNOWLEDGEMENTS
Jwo Ming Jou (Project No.: MOST103-2221-E-230-007).
368 Mechatronics and Robotics
REFERENCES
1. Lee, C.-K. (1987) Piezoelectric Laminates for Torsional and Bending
Modal Control: Theory and Experiment. Ph.D. Dissertation, Cornell
University, Ithaca.
2. Chiang, W.-W. and Lee, C.-K. (1989) Critical Active Damping Control
of a Flexible Slender Plate Using a Distributed Modal Actuator and
Sensor. American Control Conference, 700-706.
3. Vaz, A.F. and Farr, T.J. (1997) Experimental Validation of Analytical
Models of Flexible Structures with Bonded Piezoelectric Film Actuators
and Sensors. IEEE Instrumentation and Measurement Technology
Conference, 2, 887-892.
4. Vat, A.F. (1998) Composite Modeling of Flexible Structures with
Bonded Piezoelectric Film Actuators and Sensors. IEEE Transactions
on Instrumentation and Measurement, 47, 513-520. http://dx.doi.
org/10.1109/19.744200
5. Takagi, K., Nagase, K., Oshima, K. and Hayakawa, Y. (2001) Modal
Filtering for Simply Supported Plate Using Equally Segmented
Piezoelectric Film. American Control Conference, 2, 1593-1598.
6. Ballato, A. and Kim, Y.K. (2002) Displacements and Rotations of
Practical Vibrational Modes of Piezoelectric Bimorph Cantilever
Beams. Proceedings of IEEE Sensors, 2, 1294-1297.
7. Gou, X.K. and Tian, H.M. (2007) Active Vibration Control of a
Cantilever Beam Using Bonded Piezoelectric Sensors and Actuators.
8th International Conference on Electronic Measurement and
Instruments, 4, 85-88.
8. Rios-Gutierrez, M. and Silva-Navarro, G. (2009) Finite Element
and Modal Modeling of a Cantilever Beam with Piezoelectric Patch
Actuator for Vibration Absorption. 6th International Conference on
Electrical Engineering Computing Science and Automatic Control,
1-6.
9. Sanchez-Rojas, J.L., Hernando, J., Donoso, A., Bellido, J.C., Gutierrez,
L.M., Ababneh, A., Seidel, H. and Schmid, U. (2009) Piezoelectric
Modal Sensors/Actuators Based on Microplates Applying Surface
Electrode Patterning. Solid-State Sensors, Actuators and Microsystems
Conference, 2346-2349.
10. Hsu, Y.-H., Lee, C.-K. and Hsiao, W.-H. (2005) Electrical and
Mechanical Fully Coupled Theory and Experimental Verification
Theory and Simulation Analysis of the Mode Shape and ... 369
17
Experimental Study of the Response
of Transonic Diffuser Flow to a
Piezoceramic Actuator at
Diffuser Throat
Okinawa, Japan
ABSTRACT
An experimental study of the response of a piezoceramic actuator set at the
throat to a transonic diffuser is carried out by measuring wall static pressure
fluctuations and by visualizing the flow field using schlieren technique. The
visualized flow fields are captured with a digital still camera and a digital
high speed video camera. The piezo ceramic actuator is attached at the
throat of the diffuser and driven by sinusoidal amplified voltage signals. The
diffuser used in this experiment is circular arc half nozzle with the height
h* and width w of 3 mm and 25 mm, respectively. The blockage factor of
the piezoceramic actuator to the diffuser throat is 9.2% assuring the effect
of change in the throat area rather than the boundary layer disturbances. The
piezoceramic actuator is driven at the frequency of 100 Hz, 200 Hz, and 300
Hz and its amplitude is about 1 mm. It is found that the wall static pressure
fluctuations and the behavior of the shock wave clearly correspond to the
vibration of the piezo ceramic actuator for all the frequency ranges whereas
the averaged shock position remains almost unchanged. All the results
mentioned above suggest that driving the piezo ceramic actuator at the
diffuser throat can be one of the promising techniques to control unsteady
transonic diffuser flow.
Keywords: Compressible Flow; Shock Wave; Transonic Diffuser;
Piezoceramic Actuator; Throat
INTRODUCTION
The unsteady flow field in a transonic diffuser has attracted a great deal
of interest not only because of the practical industrial importance but
also because of the complexity of the flow itself. It is well known that a
shock wave in a transonic flow is basically unsteady due to the interaction
between a shock wave and other flow phenomena, i.e., the shock wave
oscillation is triggered and maintained by a local interaction between the
shock foot and the boundary later developed along the wall surface or
upstream and downstream propagating disturbances. Meier [1] reported that
shock-induced separation causes the large-amplitude, unstable oscillation
of a shock wave. In addition to the shock-induced separation, the effects
of propagating disturbances toward the shock waves on the oscillations
have been analytically, numerically and experimentally investigated [2,3],
including forced oscillations by some functions or rotating rods as a forcing
function to drive the shock waves. According to previous studies, several
factors affect unsteady shock behavior in transonic diffusers. Therefore, in
addition to revealing the causes of the shock oscillations, there have also
been attempts to eliminate or reduce these unfavorable unsteady shock
oscillations both for internal and external devices. It is quite difficult to deal
with these types of oscillations unless the disturbances around the shock wave
Experimental Study of the Response of Transonic Diffuser Flow ... 373
Piezoceramic Actuator
The piezoceramic actuator is set at the diffuser throat, as shown in Figure
1(b). In order to achieve sufficient displacement by the actuator, a bimorph
type piezoceramic actuator is adopted and driven by the circuit as shown
Figure 1(c).
Figure 2. Schlieren images for frequency 200 Hz: (a) p0/pb = 1.2; (b) p0/pb = 1.4;
(c) p0/pb = 1.6; (d) p0/pb = 1.8.
downstream with the increase in the pressure ratio. However, the gradient
of the shock displacement to the pressure ratio decreases with the pressure
ratio because the shock-induced separation causes displacement delay.
Considering the effect of the piezoceramic actuator, it is easily expected that
the location of shock waves corresponds to the throat height because the
Mach number just upstream of the shock wave depends on the ratio of the
local cross sectional area to the throat area.
From the viewpoint of small-scale observation, the shock wave oscillates
around its averaged position at the fixed frequency of the piezoceramic
actuator, which will be discussed in the following section based on images
captured by a high-speed video camera. The behavior of the shock wave
is critical to the generation of large pressure fluctuations under the various
pressure ratios.
Figure 4. Wall static pressure fluctuations at throat and x/h* = 10.0 for input
frequency 200Hz: (a) p0/pb = 1.45; (b) p0/pb= 1.50; (c) p0/pb = 1.55; (d) p0/pb=
1.55.
This is also confirmed in Figure 3, which showing xs/h* = 10.0 at a
pressure ratio of 1.55. Both the fluctuating and averaged pressure are much
lower due to the downstream displacement of the pressure measurement
position. These figures contain information about the shock position and the
shock wave behavior. In other words, Figure 4(a) shows the severe pressure
fluctuation at x/h* = 10.0, whereas the fluctuation at the throat does not. This
means that the shock wave stands between the throat and the measurement
position at x/h* = 10.0. As the pressure ratio increases, the fluctuation at x/h*
= 10.0 gradually decreases the pressure fluctuation, achieving the smallest
pressure fluctuation, as shown in Figure 4(c). The pressure fluctuation as
shown Figure 4(c) denotes the supersonic state at x/h* = 10.0.
However, the sudden change in the average pressure was observed as
shown in Figure 4(d) under the same condition, which is caused by the flow
separation downstream of the shock wave. The increase in the pressure
fluctuation indicates the shock wave displacement upstream due to the flow
separation. However, sometime after the increase in pressure, the average
value and fluctuation decrease, denoting the shock wave located downstream
of the position at x/h* = 10.0. The pressure fluctuations at the throat are
found to oscillate sinusoidally for all pressure ratios. Figure 4(d) also reveals
that the wall pressure fluctuation at a pressure ratio of 1.55 suddenly changes
its level at time t = 0.03 s, which suggests that the flow has just recovered
from an unexpected separation downstream of the shock wave.
Experimental Study of the Response of Transonic Diffuser Flow ... 379
Figure 5. RMS of wall static pressure fluctuations: (a) x/h* = 0.0; (b) x/h* =
6.67; (c) x/h* = 10.0; (d) x/h* = 13.33.
In Figure 5(b), the clear peaks at a pressure ratio p0/pb of approximately
1.2 show that a single and relatively weak shock wave appears and moves
downstream with the increase in the pressure ratio. The sudden drop in the
rms indicates the completion of the process from the subsonic-to-supersonic
transition due to the downstream displacement of the shock wave. This
process is independent of the input frequencies of the piezoceramic actuator.
Figure 5(c) shows that the rms at a pressure ratio of approximately
1.7 sudden creases, which is caused by the unexpected separation of a
boundary layer, as mentioned in the previous section. This also indicate the
shock wave existence around the position of x/h* = 10.0. However, for a
pressure ratio greater than 1.8, the rms for all input frequencies becomes
approximately zero, which suggests that the flow is completely supersonic.
Note that the source of the signal originates at the throat, and then even if
the flow measurement position is supersonic, the signal from the throat can
be detected. Figure 5(d) indicates the same variations that is all the rms for
each input frequency suddenly decrease at the pressure ratio at 2.0 denoting
the supersonic sta this postion.
Experimental Study of the Response of Transonic Diffuser Flow ... 381
Figure 6. FFT analysis of wall static pressure fluctuations at x/h* = 6.67: (a) fp
= 0 Hz; (b) fp = 100 Hz; (c) fp = 200 Hz; (d) fp = 300 Hz.
382 Mechatronics and Robotics
Figures 7(a)-(d) show the FFT analysis of the pressure fluctuation at x/h* =
10.0 for the same input frequency as that in Figure 6. Figure 7(a) shows no
dominant frequency but only a relatively low frequency, which has the same
tendency as that in Figure 6(a)
Figures 7(b)-(d) show cleominant frequencies. This indicates the
possibility of reducing the pressure fluctuation by controlling the shock
position, which is expected to cancels the pressure fluctuation.
In the figures, a wide range of frequencies are observed at a pressure
ratio of approximately 1.7 due to the unexpected shock-induced separation
of the boundary layer.
All of the figures show a sudden decrease in the power spectrum at
a pressure ratio of approximately 1.4. This implies that the shock wave
completely passes the measurement position, resulting in a change in state
from subsonic to supersonic.
Figure 8. FFT analysis of Shock positions: (a) fp = 0 Hz; (b) fp = 100 Hz; (c) fp
= 200 Hz; (d) fp = 300 Hz.
Figures 8(a)-(d) show the results of the FFT analyses of the shock
positions under the influence of the piezoceramic actuator as well as with no
control. In case of no control, there is no dominant frequency, as shown in
Figure 8(a). The relatively low-frequency spectrum is observed as long as
the pressure ratio is smaller than 1.4. This shows that the shock wave passes
the measurement position monotonically as the pressure ratio increases.
In contrast, for a shock wave frequency of 100 Hz, a dominant frequency
is observed until a pressure ratio of 1.5. The peak value of the dominant
frequency decreases as the pressure ratio increases. This decrease in the
peak value of the dominant frequency is explained by the behavior of the
shock wave, as shown in Figure 3. That is, the gradient of the shock wave
displacement decreases as the pressure ratio increase as shown in Figure 3
due to the expansion rate of the cross sectional area of the diffuser.
This implies that shock wave located downstream is insensitive to the
pressure fluctuation compared with that located upstream position. Then, the
oscillation of the shock wave gradually decreases as it moves downstream.
The diminishing of the dominant frequency for higher pressure ratios might
be due to the relative increase in other frequencies.
On the other hand, when the higher frequency is applied to the throat,
the clearly peaky dominant frequency, which is the exact same frequency as
that input to the piezoceramic actuator, can be observed. Note also that the
applied frequency remains, even for the higher pressure ratio, which suggests
that the shock wave is always under the influence of the piezoceramic
384 Mechatronics and Robotics
actuator. The tendencies observed in Figure 8(c) also appear in Figure 8(d).
In other words, the input frequency of 300 Hz is the dominant frequency for
all pressure ratios. Then, actuation at the throat is a promising method for
controlling shock wave behavior.
CONCLUSIONS
A piezoceramic actuator is applied to the throat of a circular arc diffuser with
various driving frequencies in order to clarify the response of the flow field
and shock wave behaviors to the piezoceramic actuator. The piezoceramic
actuator, which moves periodically as a reference signal, is considered to be
a driving force for oscillation phenomena. The conclusions are summarized
as follows:
1) The starting shock wave moves downstream monotonically with
the increase in the wind tunnel pressure ratio, regardless of the
input frequencies to the piezoceramic actuator.
2) The rms values of the wall pressure fluctuations decrease
suddenly just after the shock wave completely passes over the
measurement position for all driving frequencies.
3) The pressure fluctuations at the throat and downstream of the
throat correspond to the exact same frequencies of the input
frequencies to the piezoceramic actuator in both cases that the
position is supersonic and subsonic state.
4) Shock wave behaviors are also confirmed to correspond to the
piezoceramic actuator behaviors.
Experimental Study of the Response of Transonic Diffuser Flow ... 385
REFERENCES
1. G. E. A. Meier, “Shock Induced Flow oscillations in a Laval Nozzle,”
In: K. Oswatitsch and D. Rues, Eds., Symposium Transonicum II,
Springer, Berlin, 1976, pp. 252-261.
2. T. Hsieh, T. J. Bogar and T. J. Coakley, “Numerical Simulation and
Comparison with Experiment for SelfExcited Oscillations in a Diffuser
Flow,” AIAA Journal,
3. Vol. 25, No. 7, 1987, pp. 936-943. doi:10.2514/3.9725
4. R. Bur, R. Benay, A. Galli and P. Berthouze, “Experimental and
Numerical Study of Forced Shock-Wave Oscillations,” Aerospace
Science and Technology, Vol. 10, No. 4, 2006, pp. 265-278.
doi:10.1016/j.ast.2005.12.002
5. E. Benini, R. Biollo and R. Ponza, “Efficiency Enhancement in
Transonic Compressor Rotor Blades Using Synthetic Jets: A numerical
Investigation,” Applied Energy Vol. 88, No. 3, 2011, pp. 953-962,
6. doi:10.1016/j.apenergy.2010.08.006
7. P. B. Salunkhe, J. Joseph and A. M. Pradeep, “Active Feed-Back
Control of Stall in a Axial Flow Fan under Dynamic Inflow Distortion,”
Experimental Thermal and Fluid Science, Vol. 35, No. 6, 2011, pp.
1135-1142. doi:10.1016/j.expthermflusci.2011.03.008
8. M. Yaga, T. Haga and K. Oyakawa, “Study on Passive Control in a
Transonic Diffuser,” AIAA Paper 2000-0902, 38th Aerospace Sciences
Meetings & Exhibit, Reno, January 2000, pp. 1-11.
CHAPTER
18
An ARX-Based PID-Sliding Mode
Control on Velocity Tracking Control of
a Stick-Slip Piezoelectric-Driven
Actuator
Y. Cao, X. B. Chen
Department of Mechanical Engineering, University of Saskatchewan,
Saskatoon, Canada
ABSTRACT
Piezoelectric-driven stick slip actuators have been drawn more and more
attention in the nano- positioning application due to the high accuracy and
theoretical unlimited displacement. However, the hysteresis of piezoelectric
actuator (PEA) and the nonlinear friction force between the end- effector and
the stage make control of piezoelectric-driven stick slip actuator challenge.
This paper presents the development of an autoregressive exogenous
Citation:- Cao, Y. and Chen, X. (2015), “An ARX-Based PID-Sliding Mode Control
on Velocity Tracking Control of a Stick-Slip Pi-ezoelectric-Driven Actuator”. Modern
Mechanical Engineering, 5, 10-19. doi: 10.4236/mme.2015.51002.
Copyright: © 2015 by authors and Scientific Research Publishing Inc. This work is li-
censed under the Creative Commons Attribution International License (CC BY). http://
creativecommons.org/licenses/by/4.0
388 Mechatronics and Robotics
INTRODUCTION
Piezoelectric-driven stick-slip actuators play important roles in nano-
positioning applications due to their simple configuration, high accuracy and
theoretically unlimited displacement [1] - [6]. The end-effector is supported
and guided by a movable stage, which is connected to the base through
flexible hinges and driven by a PEA, as shown in Figure 1. The end-effector
moves together with the stage under slow expansion of the PEA. Sudden
contraction of the PEA leads to the sliding motion of the end-effector since
the force due to the friction limit between the end-effector and the movable
stage. As a result, the end-effector moves a step ΔS, with respect to its
original position. Such steps can be accumulated to achieve a theoretically
unlimited displacement (actually limited by the size of the moveable stage).
The hysteresis of PEA and the friction dynamics of the end-effector make
the control of stick-slip actuator challenge. An analog electronic circuit was
developed in [7] for driving stick-slip piezoelectric linear actuators. The
key design task for the amplifier is to provide a high-voltage asymmetric
sawtooth-like signal and feed it into a capacitive load.
To compensate the hysteresis effect, combination of a charge control
scheme with switch- ing is proposed and the experiments confirm the
superiority of the proposed method over other existing techniques.
However, the charge control requires expensive hardware, which increases
the difficulty and cost of its implementation.
An ARX-Based PID-Sliding Mode Control on Velocity Tracking ... 389
(1)
where are the parameters and the numerator
and denominator are relatively prime polynomials. With the use of the
system identification method, this plant model can be mathematically repre-
sented by an ARX model.
(2)
Equation (2) is identical to a discrete transfer function without zeros.
The neglect of zeros definitely leads to model errors, which is considered to
be disturbance in Equation (2).
The ARX model can be rewritten in a state space form
(3)
(5)
392 Mechatronics and Robotics
(6)
Substituting Equation (5) into Equation (4) yields
(7)
For the system described by Equation (3), one has
. Thus, Equation (7) is simplified as
(8)
For the sliding function that takes the following form of
(9)
The control action can be considered consisting two parts, i.e.,
(10)
where is the equivalent control and is the
switching control.
Substituting Equation (10) into Equations (8) and (9) yields
(11)
(12)
For the convenience of following discussion, Equations (9) and (11) are
rewritten as
(13)
An ARX-Based PID-Sliding Mode Control on Velocity Tracking ... 393
(14)
where P, I and D are parameters of the discrete PID-based SMC; T is the
sampling period.
Theorem 1: If the closed-loop system (13) is stable, the zero steady state
error can be achieved [17] .
Theorem 2: There exist some P, I and D such that the closed-loop control
system (13) is stable [16] .
EXPERIMENTS
To verify the effectiveness of the proposed method, experiments were
carried out to control the motion velocity of a stick-slip piezoelectric-
driven actuator (Figure 2). The PEA used to drive the stage generates
displacement in a range of 15 μm with a resolution of 0.5 nm. An inductive
sensor (SMU9000, Kaman, Windsor, CT) with a resolution of 10 nm was
used to measure the displacement of the end-effector. Both the actuator
and the sensor were connected to a host computer via an I/O board (PCI-
DAS1602/16, Measurement Computing Corporation) and controlled via
SIMULINK programs. The movable stage is driven by the PEA and the
end-effector on the stage can slide horizontally on the stage via a pair of V
grooves machined on them. Polished silicon wafers are attached to the V
grooves as frictional surfaces to guarantee the uniformity of the frictional
characteristics throughout the traveling range of the end-effector. All the
displacements presented in this study were measured by the inductive sensor
with a sampling interval of 0.05 ms.
Figure 3 shows the model configuration of the stick-slip piezoelectric-
driven stage developed in [18] .
394 Mechatronics and Robotics
Model Identification
To drive the end-effector, a saw-tooth generator is required, as seen in Figure
4. Each saw-tooth is 20 ms wide. The amplitude of the each saw-tooth is
defined by the reference input voltage u(t). If u(t) < 0, a reverse saw- tooth
signal is demanded to move the end-effector to the backward direction, as
shown in Figure 5.
Figure 4. Saw-tooth generator for the stick-slip motion and the velocity calcu-
lator.
(15)
where T is the tooth width, mod represents the modulus after division.
The displacement of the end-effector is measured by the inductive
sensor. The average velocity v is estimated by
(16)
where is the displacement of the end-effector in one period of
the saw-tooth signal.
The dynamics of the integral system can be regarded approximately
as a second order system according to our previous research [19] [20]. As
such, a second order ARX model was employed to describe the system
dynamics, while the hysteresis exhibited by the actuator and the friction
between the stage and the end-effector are considered as extra disturbances
to the dynamic model. For the experiments, a 10 V step input reference
voltage was provided to the PEA and the corresponding output velocity was
estimated by Equation (23). With the recorded data, the parameters were
identified by using the least square method, leading to
396 Mechatronics and Robotics
(17)
Step Tracking
The ARX-based PID-SMC developed in this study was implemented in
experiments to control the velocity of the piezoelectric-driven slick-slip
actuator. The sliding surface was defined according to Equation (9), where
and m = −0.1. To minimize the rising time, the PID parameters in
the PID regulator is adjusted to be P = 0.2, I = 3, D = 0. For comparison, the
traditional PID controller was also applied, with its parameters adjusted by
trails-and-errors to be P = 0.4, I = 4, D = 0.
Figure 6 shows the velocity step response of the piezoelectric-driven
stick-slip actuator controlled by the developed control scheme and the
traditional PID controller. Table 1 shows the comparison of the rising time
and the overshoot.
It can be seen that the tracking errors approach to zero through the
use of I component in the PID regulator. The velocity step response of the
piezoelectric-driven stick-slip actuator controlled by the proposed method
is faster than that controlled by the traditional PID controller. For example,
when a 10 μm/s step reference input was provided, the rising time of the
velocity step response controlled by the ARX-based PID-SMC is 0.13 s, 0.2
ms less than that controlled by the PID controller. Faster velocity response
is also observed for the 20 μm/s step reference input. However, since the
dynamics performance of the piezoelectric-driven stick-slip actuator varies
with the amplitude of the input voltage, the model error increased when the
20 μm/s step reference input was provided to the actuator. As a result, the
overshoot increased to 10% in this case.
Sinusoidal Tracking
To further show the effectiveness of the proposed control method, sinusoidal
tracking experiments with different frequencies were carried out on the
piezoelectric-driven stick-slip actuator.
398 Mechatronics and Robotics
0.5 - 0.75 s and 1.5 - 1.75 s is 50% of the amplitude if controlled by the ARX-
based PID-SMC. The maximum tracking error even reaches 100% of the
reference amplitude of the traditional PID controller is applied. This might
be due to the nonlinear friction force between the end-effector and moving
stage when the velocity changes from positive to negative. Obviously, the
ARX-based PID-SMC partially compensates the nonlinearity. Consideration
of the nonlinear friction model and its integration in the controller design are
required for further improvement on the velocity tracking performance of
the piezoelectric-driven stick-slip actuator.
CONCLUSION
This paper presents the development of an ARX-based PID-SMC for the
velocity control of the piezoelectric- driven stick-slip actuator. Specifically,
by applying ARX model, the output-tracking problem is defined as one of
state tracking, while the “bang-bang” switching control in SMC is replaced
with PID-based one. With the developed control scheme, chattering or the
state oscillation at a high frequency can be eliminated and also the zero
steady state error can be achieved.
To verify the effectiveness of the developed control scheme, experiments
were carried out on a piezoelectric-driven stick-slip actuator, whose
dynamics was identified by experiments.
The results of velocity tracking performance with the proposed control
scheme were compared to that with the traditional PID controller. It was
shown that both the step response and the sinusoidal tracking performance
can be further improved by applying the proposed method, although the
ARX model error exists.
ACKNOWLEDGEMENTS
Y. Cao, X. B. Chen (CSC) and the Natural Sciences and Engineering
Research Council (NSERC) of Canada is acknowledged.
An ARX-Based PID-Sliding Mode Control on Velocity Tracking ... 401
REFERENCES
1. Devasia, S. and Moheimani, S.O.R. (2007) A Survey of Control Issues in
Nano-Positioning. IEEE Transactions on Control System Technology,
15, 802-823. http://dx.doi.org/10.1109/TCST.2007.903345
2. Shiratori, H., Takizawa, M., Irie, Y., Hirata, S. and Aoyama, H. (2012)
Development of the Miniature Hemispherical Tilt Stage Driven by
Stick-Slip Motion Using Piezoelectric Actuators. Mechatronics-REM,
Paris, 21-23 November.
3. Zhang, Q.S., Chen, X.B. Yang, Q. and Zhang, W.J. (2012) Development
and Characterization of a Novel Piezoelectric-Driven Stick-Slip
Actuator with Anisotropic-Friction Surfaces. International Journal
of Advanced Manufacture Technology, 61, 1029-1034. http://dx.doi.
org/10.1007/s00170-011-3771-y
4. Eigoli, A.K. and Vossoughi, G. (2010) Dynamic Modeling of Stick-
Slip Motion in a Legged Piezoelectric Driven Microrobot. International
Journal of Advanced Robotic Systems, 7, 201-208.
5. Kang, D. (2012) Modeling of the Piezoelectric-Driven Stick-Slip
Actuators. Master Thesis, College of Graduate Studies and Research,
University of Saskatchewan, Saskatoon.
6. Bergander, A. and Breguet, J.M. (1998) Performance Improvement for
Stick-Slip Positioners. International Symposium on Micromechatronics
and Human Science, 59-66.
7. Spiller, M. and Hurak, Z. (2011) Hybrid Charge Control for Stick-Slip
Piezoelectric Actuators. Mechatronics, 21, 100- 108. http://dx.doi.
org/10.1016/j.mechatronics.2010.09.002
8. Rakotondrabe, M., Haddab, Y. and Lutz, P. (2008) Voltage/Frequency
Proportional Control of Stick-Slip Micropositioning Systems. IEEE
Transactions on Control System Technology, 16, 1316-1322. http://
dx.doi.org/10.1109/TCST.2008.917232
9. Rakotondrabe, M., Haddab, Y. and Lutz, P. (2009) Development of
Modeling and Control of Mircro/Nano-Positioning 2-DOF Stick-Slip
Device. IEEE Transactions on Mechatronics, 14.
402 Mechatronics and Robotics
10. Breguet, J.M. and Clavel, R. (1998) Stick and Slip Actuators: Design,
Control, Performances and Applications. International Symposium on
Micromechatronics and Human Sciences, Nagoya, 25-28 November
1998, 89-95.
11. Young, K.D., Utkin, V.I. and Özgüner, Ü. (1999) A Control Engineer’s
Guide to Sliding Mode Control. IEEE Transactions on Control System
Technology, 17, 328-342. http://dx.doi.org/10.1109/87.761053
12. Edwards, C. and Spurgeon, S.K. (1998) Sliding Mode Control. Taylor
& Francis Ltd, London.
13. Liawa, H.C., Shirinzadeh, B. and Smith, J. (2007) Enhanced Sliding
Mode Motion Tracking Control of Piezoelectric Actuators. Sensors and
Actuators, A138, 194-202. http://dx.doi.org/10.1016/j.sna.2007.04.062
14. Wang, S., Habibi, S., Burton, R. and Sampson, E. (2006) Sliding Mode
Control of an Electrohydraulic Actuator System with Discontinuous
Nonlinear Friction. Proceedings of the Institution of Mechanical
Engineers, 222, Part I:J, Systems and Control Engineering.
15. Peng, J.Y. and Chen, X.B. (2014) Integrated PID-Based Sliding Mode State
Estimation and Control for Piezoelectric Actuators. IEEE/ASME Transactions
on Mechatronics, 19, 88-99. http://dx.doi.org/10.1109/TMECH.2012.2222428
16. Cao, Y. and Chen, X.B. (2014) Disturbance Observer Based Sliding
Mode Control for a Three-DOF Nano-Positioning Stage. IEEE/ASME
Transactions on Mechatronics, 19, 924-931. http://dx.doi.org/10.1109/
TMECH.2013.2262802
17. Cao, Y. and Chen, X.B. (2012) Integrated Inversion-Feedforward
and PID-Based Sliding Mode Control for Piezoelectric Actuators.
American Control Conference, Montreal, 27-29 June 2012.
18. Peng, J.Y. and Chen, X.B. (2011) Modeling of Piezoelectric-Driven
Stick-Slip Actuators. IEEE/ASME Transactions on Mechatronics, 16,
394-399. http://dx.doi.org/10.1109/TMECH.2010.2043849
19. Cao, Y. and Chen, X.B. (2012) A Novel Discrete ARMA-Based
Model for Piezoelectric Actuator Hysteresis. IEEE/ ASME
Transactions on Mechatronics, 17, 737-744. http://dx.doi.org/10.1109/
TMECH.2011.2128339
20. Chen, X.B., Zhang, Q.S., Kang, D. and Zhang, W.J. (2008) On the
Dynamics of Piezoelectric Positioning Systems. Review of Scientific
Instrument, 79, Article ID: 116101. http://dx.doi.org/10.1063/1.2982238
INDEX
Symbols 389
Automated Monitoring 254
1-degree-of-freedom (1-DOF) 274 automatic irrigation control system
1-DOF control problems 300 65
A automatic system 70
Automation 68, 74
Acoustic and Magnetic Signals 224 Automation component 68
Actuators 320, 323, 342, 343, 344 automotive manufacturing industry
Adaptive Control 102 66
adaptive controller 106, 115 Autonomous Microbe Sensor 253,
ad-hoc relay communication. 213 254, 271
aerospace 170 Autoregresive Progress 388
agricultural production 68 auto-regressive (ARX) 390
Anisotropic Magneto Resistive autoregressive exogenous (ARX)
(AMR) 232 388
ANSYS Code 346
arc welding 82 B
Arduino Robot Control 82 bang-bang 389, 399
Artificial Intelligence 12, 14 Barbalat’s lemma 109
Atomic Force Microscope (AFM)
404 Mechatronics and Robotics
J 62
Mechatronic component 18, 30
Japanese Robotics Association 5
Mechatronics 18, 38, 39
K micro-assembly system 389
Microbial Sensor 254
kinematics model 170, 181 microbiologically influenced corro-
Known Dynamic Environment 120 sion (MIC) 254
L microcontroller unit 225
Micro-Electro-Mechanical Systems
Laboratory calibration data 260 321
Labview interface 96, 99 micro-electro mechanical systems
Labview Robot Control 82 (MEMS) 321
Lab view software 95 microprocessors 6
Ladder logic 71, 72 Military Applications 226
Ladder logic programming 72 mobile nodes 206, 208, 210, 211,
laser scanners 12 212, 213, 214, 215, 216, 217
linear dynamic system 294 mobile sensor network 206, 207
Liqueur production mechatronic mobile sensor nodes 206, 207, 208,
system 28 210, 211, 219
logical reasoning 3, 13 Model Driven Development 44
Low-Cost 81, 82 model-driven development of com-
low-frequency spectrum 383 plex mechatronic systems
Lyapunov-based methods 103 (MTSs) 45
Lyapunov function 107, 110 modeling flexible robots 103
Lyapunov stability theory 103, 114, Model Integrated Mechatronics
115 (MIM) 44, 45, 47
M Mode Shape 345, 346
Monitoring Biodiversity 226
Magnetic sensors 232, 233 Monitoring Nuclear Reactor 226
Managing Inventory Control 226 Monitoring through Internet 226
manipulator control system 170 mono-static and mono-mobile net-
Maxwell resistive capacitor (MRC) work 217
278 Monterrey University of Technol-
mechanical design 83, 84 ogy 83
mechanical engineering 18 motors 66, 70, 73
mechanical manufacturing 170 Multicriteria optimization 142
mechanical unit 18, 29
mechatronic component 45, 49, 51, N
Nano positioning 277
Index 407
navigation 205, 206, 207, 209, 212, Photostriction 320, 324, 326, 342
213, 215, 216, 217, 218, 219 Photostrictive effect of PLZT 325
Network communication 66 Photostrictive materials 321, 322,
non-centrosymmetric 320, 324 326
non-mapping navigation 207, 213 photostrictive optical actuator 326,
normally closed (NC) 71 327, 332
normally open (NO) 71 Physical Object and a cyber-physi-
Normal Shape 345, 346 cal IoT enabler 25
nuclear 170 PID-based SMC 389, 393
numerical simulations 152 PID control 291, 300
PID regulator 389, 396
O piezoceramic actuator 371, 373,
object- oriented (OO) 23 374, 375, 377, 379, 380, 381,
Open-loop control schemes 289 383, 384
opera- tor panel 66 Piezoelectric Actuators 273, 274,
optimal algorithm 207, 213 275, 303, 306, 307, 308, 309,
optimal design 149, 151, 154, 158, 310, 311, 315
162, 163 Piezoelectric actuators (PEAs) 273,
optimal path 121, 123, 135 274
Optimal robot path planning 142 piezoelectric cantilever 346, 347,
Optimal task placement 142 349, 351, 356, 358, 360
out-of-plane bending deformation Piezoelectric-driven stick-slip actua-
183, 185, 189 tors 388
piezoelectric film (PVDF) 185, 189
P piezo-resistance 184
painting 82 Pilot experiment 74
Particle Swarm Optimization (PSO) Platform Independent Model (PIM)
120 22
passive algorithm 207 PLC’s internal memory 69
path planning 119, 120, 121, 122, polynomial based linear mapping
123, 126, 131, 135 strategy 284
PEA-driven parallel mechanism 300 Power consumption 141
PEA hysteresis 277, 282 Power supply 69
Pediatric Research Clinical Institute practical industrial 372
xix PrandtlIshlinskii (PI) 278
penalized estimate 210 Preisach hysteresis model 278, 279,
Personalized Medicine xvi 281, 282, 283, 288, 290, 293,
petroleum chemistry 170 295
photo-ferroelectric (PFE) 321 Preisach hysteresis operator 278,
408 Mechatronics and Robotics