Vibration observation for a translational flexiblelink manipulator based on improved Luenberger observer
Jinyong Ju^{1} , Wei Li^{2} , Yuqiao Wang^{3} , Mengbao Fan^{4} , Xuefeng Yang^{5} , Yufei Liu^{6}
^{1, 2, 3, 4, 5, 6}School of Mechatronic Engineering, China University of Mining and Technology, Xuzhou 221116, P. R. China
^{2}Corresponding author
Journal of Vibroengineering, Vol. 18, Issue 1, 2016, p. 238249.
Received 10 July 2015; received in revised form 26 August 2015; accepted 6 October 2015; published 15 February 2016
JVE Conferences
The residual vibration is a very universal problem in flexible manipulators which are widely used in robot technology. This paper focuses on the soft measurement of the vibration signals for a translational flexiblelink manipulator (TFLM) system. A vibration observer based on the improved Luenberger observer, which only requires the practical measurement values of the boundary positions, is designed to obtain the vibration signals of the TFLM. The main contribution of the vibration observer is its ability to simplify system structure and get the vibration signals of any point of the TFLM which is unrealistic by infinite sensors in practice. Furthermore, the improved part of the Luenberger observer is the added feedback coefficients for the tip vibration signals which can correct the observed mode and reduce the observation error markedly. And according to the stable conditions of observer, the added feedback coefficients are designed by Lyapunov technique and multiple population genetic algorithms (MPGA). Finally, the efficiency of the designed vibration observer is verified by combinedsimulation.
Keywords: TFLM, vibration observer, observation error, tip vibration, Lyapunov technique, MPGA.
1. Introduction
With the highprecision and highspeed development of robot’s technology, the requirements for the mechanical arm are increasingly. Compared with the rigid manipulators, the flexiblelink manipulator has numerous advantages, such as smaller damping, higher speed and larger load. As a consequence, it has increasingly applied in many engineering fields such as precision assembly [1], modern manufacturing [2, 3], medical equipment [4], space equipment [5]. However, the flexiblelink manipulator exhibits long time residual vibrations which are aroused during the positioning and seriously impacts the positioning speed and accurate tracking ability. Because the flexiblelink manipulator is essentially an infinitedimensional distributedparameter system and the actuators are fewer than its degrees of freedom which means that the flexiblelink manipulator system is underactuated, the residual vibration control for the flexiblelink manipulator is rather challenging [68].
In order to constitute effective feedback control to suppress the residual vibrations of the translational flexiblelink manipulator (TFLM), the vibration signals should be obtained firstly. Currently, the researchers often measure the vibration signals of TFLM by sensors such as accelerometers [9], strain gauges [10], positionsensitive detector (PSD) [11] and so on. Qiu [12] presented a characteristic model based nonlinear golden section adaptive control algorithm for the TFLM system and significantly suppressed the larger amplitude vibration and the lower amplitude vibration. Malki and Misir [13] applied a fuzzy PID controller to the vibration control of the flexiblejoint manipulator and verified the effectiveness and robustness of this control method. In [14], an active vibration controller based on fuzzy logic and neural networks for the vibration suppression of the singlelink flexible manipulator was proposed and the experimental results showed its feasibility. Unfortunately, because the low stiffness and large deflection of the TFLM, the addition of sensors is bound to affect dynamic characteristics of the system [15]. Thus, many researchers investigate the location configuration, structure optimization of sensors and control mechanism optimization for reducing the influence of sensors on the system. Sun [16] described the validity of piezoelectric materials (PZT) for the vibration suppression of a singlelink flexible manipulator and the PSD sensor, which is a kind of disconnection displacement sensor, is applied to obtained the vibration signals of the flexible manipulator. In [17], a novel observerbased robust dynamic feedback tracking controller without velocity measurements was developed to the motion tracking control of flexiblejoint robotic manipulators and the number of sensors can be greatly reduced. An optimization problem with genetic algorithms (GA) is formulated in [18], which is aim to optimize the location and size of the sensors on the flexible manipulators. However, in those methods, the influences of sensors on the dynamic characteristics and structural optimization still exist while the sensors are used to obtain the vibration signals. Besides, the residual vibration of the TFLM is often suppressed by adding additional actuators [19, 20]. Due to the position limitation, the sensors can only obtain the residual vibration signals near the actuators which seriously affects the control effect. Furthermore, because the actuators positions should be optimized firstly [11, 16], the vibration signals of any point of the TFLM is indispensable which is unrealistic by infinite sensors.
Fortunately, the Luenberger observer can estimate the internal state variables by the measurable input and output variables which is a kind of method for indirect measurement [21]. Based on this, the vibration signals of the TFLM can be effectively estimated by rational Luenberger observer. In [22], an observer based on partial differential equation (PDE) model was presented to achieve the estimate of the infinite dimensional states for a flexiblelink manipulator. Mosayebi [23] developed a nonlinear high gain observer to estimate the elastic degrees of freedom and their time derivatives for flexiblelink manipulators. But those observer is only based on the feedback values of the joint angles which is an openloop observation method for the elastic degrees of freedom and the observation error cannot be solved. Above all, there are two major methods in observer design of the TFLM: one is based on the ordinary differential equation (ODE) models [2325] and another is based on PDE models [22]. Through comparing these two methods, although the PDE models can describe the TFLM system more precisely, the complicated calculations seriously affect the effect of realtime observation. Then, many researchers apply ODE models to construct observer in practical application. On the other side, because observers are based on the theoretical models, the intrinsic deviations of ODE models, which include the imprecise part caused by mode truncation and the nonlinear part, obviously affect the observation accuracy. Thus, the solutions of observation error should be further analyzed.
Therefore, based on the ODE dynamic model of the TFLM, a vibration observer is constructed in this paper. In this way, the vibration signals of any point on the TFLM can be got with limited sensors. In view of the modal truncation errors of ODE model, the tip vibration signals are fed back to improve the observation precision of the other positions of the TFLM. The design process of the vibration observer gains deals with two parts. The basic part is based on the thought of no deviation Kalman estimator and the part of observation error compensation is designed by Lyapunov technique and multiple population genetic algorithms (MPGA).
The remainder of this paper is organized as follows. The dynamic modeling of TFLM is given in Section 2. The main contribution of this paper is described in Section 3, including the design and analysis of the observer. Section 4 provides the union simulations to illustrate the effectiveness of the proposed control strategy. Finally, conclusions are drawn in Section 5.
2. Dynamic modelling of the TFLM
Schematic diagram of a TFLM system is shown in Fig. 1. During the construction of the dynamic model, assumptions are made as follows: (a) The TFLM is simplified as a EulerBernoulli beam and the impact of shear and axial deformation can be neglected. (b) The connection between the TFLM and the slider is pure rigid. (c) The influence of gravity can be neglected.
Fig. 1. Schematic diagram of the TFLM system
In accordance with the assumed mode methods [24, 25], the axial absolute coordinates of $H$ which is a random point on the TFLM can be expressed as:
where $Z\left(t\right)$ is the slider displacement, $\omega ({x}_{l},t)$ denotes the elastic deformation of $H$, ${\varphi}_{i}\left({x}_{l}\right)$ is the $i$th modal shape, ${q}_{i}\left(t\right)$ is the modal coordinates. The boundary conditions, which are one end fixed and another end free, formulate the modal shape as follows:
where ${\beta}_{i}$ is constant coefficient which is also determined by the boundary conditions.
Through dynamic analysis of the TFLM, the kinetic energy of the system can be written as:
where ${m}_{b}$ is mass of the slider, ${m}_{t}$ is mass of the payload. ${\rho}_{b}$, ${A}_{b}$ and ${L}_{b}$ are the density, the sectional area and the length of the TFLM, respectively.
Potential energy of the system mainly consider the elastic potential energy caused by the elastic deformation of the TFLM and can be expressed as:
where ${E}_{b}$ and ${I}_{b}$ are the elastic modulus and moment of inertia of the TFLM, respectively.
The external forces of the TFLM primarily include the driving force and the friction between the slider and the linear guide. The driving force $F\left(t\right)$, as the control input of TFLM system, is generated by the motor through a reducer and a ball screw. So the virtual work [26] of the TFLM system is represented by:
where ${\upsilon}_{c}$ is the friction coefficient between the slider and the linear guide and ${\mu}_{s}$ is the structural damping of the TFLM. According to the second Lagrange equation [27] as follows:
where ${Q}_{j}^{\text{'}}$ denotes generalized force which corresponds to generalized coordinates ${q}_{j}^{\text{'}}$. $L={T}_{e}U$ is Lagrange multiplier and can be expressed as:
By substituting Eq. (1) into Eq. (7) and Lagrange approach, the dynamic model of the TFLM can be deduced which is showed as:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=F\left(t\right){v}_{c}\dot{Z}\left(t\right),$
where ${\lambda}_{i}$ is the $i$th natural frequency of the TFLM. According to the orthogonality of the modal shape, Eq. (9) can be further simplified as:
By defining ${m}_{i}={\rho}_{b}{A}_{b}\underset{0}{\overset{{L}_{b}}{\int}}{\varphi}_{i}\left({x}_{l}\right)d{x}_{l}$, Eq. (8) and Eq. (10) can be converted as:
Taking the state variables as $x=[\begin{array}{llllllllll}Z& {q}_{1}& {q}_{2}& ...& {q}_{i}& \dot{Z}& {\dot{q}}_{1}& {\dot{q}}_{2}& ...& {\dot{q}}_{i}\end{array}{]}^{T}$, the system dynamic equation can be transformed into the state space equation whose form is:
where:
3. Observer design
Through the analysis of the observability [28], the TFLM system is completely observable, which meet the existence conditions of a Luenberger observer. By using the dynamic model of the TFLM in Section 2, the vibration observer based on Luenberger observer is constructed, as illustrated in Fig. 2. Where $\widehat{\omega}({L}_{b},t)$ expresses the estimates of tip vibration signals of the TFLM. $\widehat{Z}\left(t\right)$ and $\widehat{\omega}({x}_{l},t)$ are denoted as the estimates of $Z\left(t\right)$ and $\omega ({x}_{l},t)$, where ${x}_{l}\in \left[0,{L}_{b}\right]$. Then, the estimate of $Y({x}_{l},t)$ can be expressed as:
Fig. 2. Principle diagram of vibration observer
Because the number of modal shapes only affect the complexity of the calculation, the following illustrates the design process of vibration observer by only considering the firstorder mode. Vibration observer equation can be written as:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\mathbf{A}\widehat{x}+\mathbf{B}u+{K}_{e}\stackrel{~}{Z}+{K}_{p}\stackrel{~}{\omega}\left({L}_{b},t\right),$
where ${K}_{e}$ and ${K}_{p}$ are the designed observer gains. $\stackrel{~}{Z}\left(t\right)$ and $\stackrel{~}{\omega}({L}_{b},t)$ express the observation error of the slider and the tip vibration displacement, respectively.
Firstly, $Z\left(t\right)$ and $\omega ({L}_{b},t)$ are directly measured through sensors. Subtracting Eq. (14) from Eq. (12), the error model of vibration observer can be obtained as:
where ${\mathbf{C}}_{1}=\left[\begin{array}{llll}1& 0& 0& 0\end{array}\right]$ and $\stackrel{~}{x}$ represents the observation error of the state variables ($x$) and $\dot{\stackrel{~}{x}}\to \dot{x}$. The designed observer gains of vibration observer compose of two parts: ${K}_{e}$ is based on the thought of no deviation Kalman estimator and ${K}_{p}$ is designed by the purpose to guarantee the stability of the observer system and minimize observation error of the vibration signals. The form of ${K}_{e}$ appears as follow:
where $\mathbf{P}$ is the positive definite solution of Ricatti equation by defining appropriate $\mathbf{Q}$ and $\mathbf{R}$ which are positive definite diagonal matrix. The form of Ricatti equation can be written as:
In order to guarantee the stability of the observer system, the Lyapunov function of the error model is defined as:
Differentiating $V$ with respect to time yields:
${\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\left[\right(\mathbf{A}{K}_{e}{\mathbf{C}}_{1})\stackrel{~}{\mathbf{x}}{K}_{p}\stackrel{~}{\omega}({L}_{b},t\left)\right]}^{T}{\mathbf{P}}^{1}\stackrel{~}{\mathbf{x}}+{\stackrel{~}{\mathbf{x}}}^{T}{\mathbf{P}}^{1}\left[\right(\mathbf{A}{K}_{e}{\mathbf{C}}_{1})\stackrel{~}{\mathbf{x}}{K}_{p}\stackrel{~}{\omega}({L}_{b},t\left)\right]$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\left[\right(\mathbf{A}{K}_{e}{\mathbf{C}}_{1})\stackrel{~}{\mathbf{x}}{]}^{T}{\mathbf{P}}^{1}\stackrel{~}{\mathbf{x}}+{\stackrel{~}{\mathbf{x}}}^{T}{\mathbf{P}}^{1}(\mathbf{A}{K}_{e}{\mathbf{C}}_{1})\stackrel{~}{\mathbf{x}}{\stackrel{~}{\mathbf{x}}}^{T}[{\mathbf{C}}_{2}^{T}{\varphi}_{1}\left({L}_{b}\right){K}_{p}^{T}{\mathbf{P}}^{1}$
$+{\mathbf{P}}^{1}{K}_{p}{\varphi}_{1}\left({L}_{b}\right){\mathbf{C}}_{2}]\stackrel{~}{\mathbf{x}}={\stackrel{~}{\mathbf{x}}}^{T}{\mathbf{A}}^{T}{\mathbf{P}}^{1}\stackrel{~}{\mathbf{x}}{\stackrel{~}{\mathbf{x}}}^{T}{\mathbf{C}}_{1}^{T}{K}_{e}^{T}{\mathbf{P}}^{1}\stackrel{~}{\mathbf{x}}+{\stackrel{~}{\mathbf{x}}}^{T}{\mathbf{P}}^{1}\mathbf{A}\stackrel{~}{\mathbf{x}}{\stackrel{~}{\mathbf{x}}}^{T}{\mathbf{P}}^{1}{K}_{e}{\mathbf{C}}_{1}\stackrel{~}{\mathbf{x}}$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}{\stackrel{~}{\mathbf{x}}}^{T}{\mathbf{P}}^{1}\left[\mathbf{P}{\mathbf{C}}_{2}^{T}{\varphi}_{1}\right({L}_{b}){K}_{p}^{T}+{K}_{p}{\varphi}_{1}({L}_{b}\left){\mathit{C}}_{2}\mathit{P}\right]{\mathit{P}}^{1}\stackrel{~}{\mathbf{x}}$
$={\stackrel{~}{\mathit{x}}}^{T}({\mathbf{A}}^{T}{\mathbf{P}}^{1}+{\mathbf{P}}^{1}\mathbf{A}{\mathbf{C}}_{1}^{T}{K}_{e}^{T}{\mathbf{P}}^{1}{\mathbf{P}}^{1}{\mathrm{K}}_{\mathrm{e}}{\mathbf{C}}_{1})\stackrel{~}{\mathbf{x}}{\dot{V}}_{2}$
$={\stackrel{~}{\mathbf{x}}}^{T}{\mathbf{P}}^{1}\left(\mathbf{P}{\mathbf{A}}^{T}+\mathbf{A}\mathbf{P}\mathbf{P}{\mathbf{C}}_{1}^{T}{K}_{e}^{T}{K}_{e}{\mathbf{C}}_{1}\mathbf{P}\right){\mathbf{P}}^{1}\stackrel{~}{\mathbf{x}}{\dot{V}}_{2}.$
where:
By substituting Eq. (18) into Eq. (19) and combining with Eq. (17), Eq. (19) can be further simplified as:
$={\stackrel{~}{\mathbf{x}}}^{T}{\mathit{P}}^{1}\left(\mathit{Q}+{\mathit{P}}^{T}{\mathbf{C}}_{1}^{T}{\mathit{R}}^{T}{\mathit{C}}_{1}\mathit{P}\right){\mathit{P}}^{1}\stackrel{~}{\mathbf{x}}{\dot{V}}_{2}.$
Because $\mathbf{Q}$ and $\mathbf{R}$ are symmetric positive definite matrix and $\mathbf{P}$ is positive definite matrix too, the first item of $\dot{V}$ is obviously less than zero which is shown as:
It can be got form Eq. (21) that ${\dot{V}}_{1}=0$ only if $\stackrel{~}{\mathit{x}}=0$. Then in order to pledge $\dot{V}\le 0$, Eq. (22) which the coefficients of ${K}_{p}$ are satisfied can be obtained:
With the above analysis, the designed observer system is asymptotically stable when the observer gains satisfy Eq. (16) and Eq. (22). Because Eq. (22) only determines the range of ${K}_{p}$, the MPGA is adopted to achieve the optimization of ${K}_{p}$. To implement the evolution steps of MPGA, the conventional strategies of selection, crossover, mutation and the unique strategies of immigration, artificial selection are necessary. In the optimization process of ${K}_{p}$, MPGA can search the optimum and change the best individuals in multipopulations by immigration. Then the purpose of cooperative coevolution can be achieved. Artificial selection is added to avoid the outstanding individuals being destroyed by crossover and mutation operation in the population. The procedure of optimizing ${K}_{p}$ is shown in Fig. 3.
In order to achieve the purpose of minimizing the observation error, the fitness function ($J$) is defined as:
where ${t}_{f}$ is the terminal time of simulation. The optimized target is minimize $J$. Termination condition is the retentive generations of the optimal individual which can take full advantages of the knowledge accumulation ability of GA and guarantee ${K}_{p}$ convergence towards global optimum.
Fig. 3. Flow chart of optimizing ${K}_{p}$
4. Simulation experiment of the TFLM
To illustrate the effectiveness of the designed vibration observer, the physical model of TFLM, which is constructed in ADAMS, is regarded as the controlled object. The material of the TFLM is stainless and relevant physical parameters are listed in Table 1.
The physical model of TFLM consists of two parts: one part is the slider and the support, the other is the flexible manipulator. The two parts are concatenated together by fixed constraints. The vibration signals of any point of the TFLM can be figured out by ADAMS/Flex. The input is defined as the driving force of slider and the outputs are set to the displacements of slider and the desired points of TFLM. Then, the physical model is converted into Simulink subsystem which is used to establish combinedsimulation platform.
Table 1. Physical parameters of the TFLM
Project

Length
${L}_{b}$ / m

Density
${\rho}_{b}$ / kg.m^{3}

Sectional area ${A}_{b}$ / m^{2}

Modulus of elasticity ${E}_{b}$ / GPa

Slider quality ${m}_{b}$_{}/ kg

Payload quality ${m}_{t}$ / kg

Value

0.635

7850

2.83e5

197

0.5

0.2

Fig. 4. Five times evolutionary process by MPGA
As indicated in [29], the low order modes play a leading role of the vibrations and the higher modes can be neglected in the assumed mode methods. Thus, in order to simplify the computation, only the firstorder mode is taken into consideration in the process of optimizing ${K}_{p}$. The total simulation time is 2.5 s and the retentive generations of the optimal individual is 4 which is deemed as the termination condition. The optimization results of ${K}_{p}$ is shown in Fig. 4.
In Fig. 4, the optimal values are consistent in five times evolutionary process which means that the evolutionary process of ${K}_{p}$ is sufficient by MPGA. Through the above analysis, the optimal observer gains of vibration observer are ${K}_{e}={\left[\begin{array}{llll}\text{90.0}& \u2013\text{4.3}& \text{271.2}& \u2013\text{112.9}\end{array}\right]}^{T}$ and ${K}_{p}={\left[\begin{array}{llll}\u2013\text{2000.0}& \u2013\text{1412.1}& \text{200.0}& \text{199.7}\end{array}\right]}^{T}$.
Afterwards, the tracking effect of observer is examined, under the initial state of $Z\left(0\right)=$0.3 m which is seen as the initial disturbances. The tracking effects of slider displacement ($Z\left(t\right)$) and slider speed ($\dot{Z}\left(t\right)$) are shown in Figs. 5 and 6.
It is obtained from Figs. 5, 6 that the designed observer can quickly track the slider displacement and speed whether in the uniform acceleration or constant velocity phase. The tracking accuracy of slider displacement and speed is high and the initial disturbances can be quickly eliminated.
Fig. 5. Tracking effect of slider displacement
In order to verify the tracking effect of the designed observer for the vibration signals of TFLM, the vibration signals of ${x}_{l}=\mathrm{}$0.4 m and ${x}_{l}=\mathrm{}$0.5 m are selected to be observered by compared with conventional Luenberger observer. The tracking effects are shown in Figs. 7, 8, where VS indicates the actual vibration signals, CO means the estimation of conventional Luenberger observer and ILO indicates the estimation of improved Luenberger observer.
Fig. 6. Tracking effect of slider speed
Fig. 7. Tracking effect of vibration signals (${x}_{l}=$ 0.4 m, ${x}_{l}=\mathrm{}$0.5 m)
Fig. 8. Tracking error of vibration signals (${x}_{l}=\mathrm{}$0.4 m, ${x}_{l}=\mathrm{}$0.5 m)
The results shown in Figs. 7, 8 indicate, as expected, that the designed observer is effective. Tracking trend of the vibrations is consistent with the output of ADAMS in a very short period of time whether the designed observer or the conventional Luenberger observer. Although the designed observer is based on the ODE model of the TFLM which is identical with the conventional Luenberger observer, the tracking error of amplitude is well solved after correcting the observed mode by feeding back the tip vibration signals.
5. Conclusions
In this work, a vibration observer based on improved Luenberger observer is proposed to obtain the vibration signals of the TFLM. Then a physical model of the TFLM is established by ADAMS software to verify the effectiveness of the vibration observer. The results demonstrate that the proposed vibration observer has a good observation effect of the vibration signals, meanwhile reduces the observational vibration error remarkably by compared with the conventional Luenberger observer. Through comparing observation results with the ADAMS outputs, the practicability of the designed observer is proved. The results obtained in this paper can be applied for the robot manipulators in such as precision assembly, modern manufacturing, medical equipment, space equipment and has positive effects on motion planning, optimization of actuators positions and residual vibrations control for flexiblelink manipulator in future research. Moreover, it has positive significance for the structural optimization and cost reduction in practical application by replacing the sensors.
Acknowledgements
This research work was partially supported by the National Natural Science Foundation of China (Nos. 51305444 and 51307172), the Doctoral Fund of Ministry of Education under Grants (No. 20120095120013), the Scientific and Technological Project of Jiangsu Province (BY201402806), the Jiangsu Provincial Natural Science Foundation of China (BK2012567), the Postgraduate Cultivation and Innovation Project of Jiangsu Province (KYLX_1375), and the Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD). The authors would like to thank the editor and the reviewers for their comments and suggestions which have improved the presentation.
References
 Dwivedy S. K., Eberhard P. Dynamic analysis of flexible manipulators, a literature review. Mechanism and Machine Theory, Vol. 41, Issue 7, 2006, p. 749777. [Search CrossRef]
 Abe A. Trajectory planning for residual vibration suppression of a twolink rigidflexible manipulator considering large deformation. Mechanism and Machine Theory, Vol. 44, Issue 9, 2009, p. 16271639. [Search CrossRef]
 Katsura S., Suzuki J., Ohnishi K. Pushing operation by flexible manipulator taking environmental information into account. IEEE Transactions on Industrial Electronics, Vol. 53, Issue 5, 2006, p. 16881697. [Search CrossRef]
 Riviere C. N., Ang W. T., Khosla P. K. Toward active tremor canceling in handheld microsurgical instruments. IEEE Transactions on Robotics and Automation, Vol. 19, Issue 5, 2003, p. 793800. [Search CrossRef]
 Gouliaev V. I., Zavrazhina T. V. Dynamics of a flexible multilink cosmic robot manipulator. Journal of Sound and Vibration, Vol. 243, Issue 4, 2001, p. 641657. [Search CrossRef]
 Rahimi H. N., Nazemizadeh M. Dynamic analysis and intelligent control techniques for flexible manipulators: a review. Advanced Robotics, Vol. 28, 2014, p. 6376. [Search CrossRef]
 Shan J, Liu H. T., Sun D. Modified input shaping for a rotating singlelink flexible manipulator. Journal of Sound and Vibration, Vol. 285, Issue 1, 2005, p. 187207. [Search CrossRef]
 Hassan M., Dubay R., Li C., Wang R. Active vibration control of a flexible onelink manipulator using a multivariable predictive controller. Mechatronics, Vol. 7, Issue 6, 2007, p. 311323. [Search CrossRef]
 Staufer P., Gattringer H. State estimation on flexible robots using accelerometers and angular rate sensors. Mechatronics, Vol. 22, Issue 8, 2012, p. 10431049. [Search CrossRef]
 Ramos F., Feliu V. New online payload identification for flexible robots. Application to adaptive control. Journal of Sound and Vibration, Vol. 315, Issues 12, 2008, p. 3457. [Search CrossRef]
 Shan J. J., Liu H. T., Sun D. Slewing and vibration control of a singlelink flexible manipulator by positive position feedback (PPF). Mechatronics, Vol. 15, Issue 4, 2005, p. 487503. [Search CrossRef]
 Qiu Z. C. Adaptive nonlinear vibration control of a Cartesian flexible manipulator driven by a ball screw mechanism. Journal of Mechanical Systems, Vol. 30, 2012, p. 248266. [Search CrossRef]
 Malki H. A., Misir D. Fuzzy PID control of a flexiblejoint robot arm with uncertainties from imevarying loads. IEEE Transactions on Control Systems Technology, Vol. 5, 1997, p. 371378. [Search CrossRef]
 Jnifene A., Andrews W. Experimental study on active vibration control of a singlelink flexible manipulator using tools of fuzzy logic and neural networks. IEEE Transactions on Instrumentation and Measurement, Vol. 54, Issue 3, 2005, p. 12001208. [Search CrossRef]
 João C. P. R., José S. D. C. Motion planning and actuator specialization in the control of activeflexible link robots. Journal of Sound and Vibration, Vol. 331, 2012, p. 32553270. [Search CrossRef]
 Sun D., Mills J. K., Shan J. J., Tso S. K. A PZT actuator control of a singlelink flexible manipulator based on linear velocity feedback and actuator placement. Mechatronics, Vol. 14, Issue 4, 2004, p. 381401. [Search CrossRef]
 Chang Y. C., Yen H. M. Robust tracking control for a class of electrically driven flexiblejoint robots without velocity measurements. International Journal of Control, Vol. 85, Issue 2, 2012, p. 194212. [Search CrossRef]
 Molter A., Da Silveira O. A. A., Fonseca J. S. O., Bottega V. Simultaneous piezoelectric actuator and sensor placement optimization and control design of manipulators with flexible links using SDRE method. Mathematical Problems in Engineering, Vol. 2010, 2010, p. 123. [Search CrossRef]
 Dadfarnia M., Jalili N., Xian B., Dawson D. M. A Lyapunovbased piezoelectric controller for flexible Cartesian robot manipulators. Journal of Dynamic Systems Measurement and Control – Transactions of the ASME, Vol. 126, Issue 2, 2004, p. 347358. [Search CrossRef]
 Qiu Z. C., Han J. D., Zhang X. M., Wang Y. C., Wu Z. W. Active vibration control of a flexible beam using a noncollocated acceleration sensor and piezoelectric patch actuator. Journal of Sound and Vibration, Vol. 326, Issues 35, 2009, p. 438455. [Search CrossRef]
 Kolhe J. P., Shaheed M., Chandar T. S., Talole S. E. Robust control of robot manipulators based on uncertainty and disturbance estimation. International Journal of Robust and Nonlinear Control, Vol. 23, Issue 1, 2013, p. 104122. [Search CrossRef]
 Yang H. J., Liu J. K., Lan X. Observer design for a flexiblelink manipulator with PDE model. Journal of Sound and Vibration, Vol. 341, 2015, p. 237245. [Search CrossRef]
 Mosayebi M., Ghayour M., Sadigh M. J. A nonlinear high gain observer based inputoutput control of flexible link manipulator. Mechanics Research Communications, Vol. 45, 2012, p. 3441. [Search CrossRef]
 Green A., Sasiadek J. Z. Dynamics and trajectory tracking control of a twolink robot manipulator. Journal of Vibration and Control, Vol. 10, Issue 10, 2004, p. 14151440. [Search CrossRef]
 Yang T. W., Xu W.L., Tso S. K. Dynamic modeling based on realtime deflection measurement and compensation control for flexible multilink manipulators. Dynamics and Control, Vol. 11, Issue 1, 2001, p. 524. [Search CrossRef]
 Vidoni R., Gasparetto A., Giovagnoni M. A method for modeling threedimensional flexible mechanisms based on an equivalent rigidlink system. Journal of Vibration and Control, Vol. 20, 2014, p. 483500. [Search CrossRef]
 Rognant M., Courteille E., Maurine P. A systematic procedure for the elastodynamic modeling and identification of robot manipulators. IEEE Transactions on Robotics, Vol. 26, Issue 6, 2010, p. 10851093. [Search CrossRef]
 Liu K., Sun X. System identification and model reduction for a singlelink flexible manipulator. Journal of Sound and Vibration, Vol. 242, Issue 5, 2001, p. 867891. [Search CrossRef]
 Kiang C. T., Spowage A., Yoong C. K. Review of control and sensor system of flexible manipulator. Journal of Intelligent and Robotic Systems, Vol. 77, Issue 1, 2015, p. 187213. [Search CrossRef]