The flexible Cartesian robotic manipulator (FCRM) is coming into widespread application in industry. Because of the feeble rigidity and heavy deflection, the dynamic characteristics of the FCRM are easily influenced by external disturbances which mainly concentrate in the driving end and the load end. Thus, with the influence of driving base disturbance and terminal load considered, the motion differential equations of the FCRM under the plane motion of the base are constructed, which contain the forced and non-linear parametric excitations originated from the disturbances of base lateral and axial motion respectively. Considering the relationship between the coefficients of the motion differential equations and the mode shapes of the flexible manipulator, the analytic expressions of the mode shapes with terminal load are deduced. Then, based on multiple scales method and rectangular coordinate transformation, the average equations of the FCRM are derived to analyze the influence mechanism of base disturbance and terminal load on the system parametric vibration stability. The results show that terminal load mainly affects the node locations of mode shapes and mode frequencies of the FCRM, and the axial motion disturbance of the driving base introduces parametric excitation while the lateral motion disturbance generates forced excitation for the transverse vibration model of the FCRM. Furthermore, with the increase of the base excitation acceleration and terminal load, the parametric vibration instability region of the FCRM increases significantly. This study will be helpful for the dynamic characteristics analysis and vibration control of the FCRM.
Cartesian coordinate robot is a common structure of industrial robots, which has been widely applied in handling operation, surface mount technology (SMT), welding operation and CNC machine (Park et al., 2013; Wu and Zhang, 2014; Srinivasan et al., 2015; Qian et al., 2016). With the development of modern machinery equipment to integration and lightweight, the flexible Cartesian robotic manipulator (FCRM) has received increasing attention. Compared with rigid structure, the FCRM has many advantages such as low energy consumption and heavy load. However, due to the low mode frequency and structure damping, the FCRM easily exhibits elastic vibration during the operation process, which has a heavy influence on the positioning accuracy and service life of the whole system, especially in high-speed operations (Neto et al., 2013). Thus, the dynamic characteristics and vibration control of the FCRM should be further analyzed (Feliu et al., 2014; Kiang et al., 2015).
The FCRM system is composed of driving unit, flexible manipulator, and end effector, which is a typical multi-coupling system. And the flexible manipulator can be assumed as a flexible beam to describe its dynamic characteristics (Mann et al., 2014). The calculated results by Zhang et al. (1995) showed that the results of the traditional modelling method were divergent when the base rotation frequency of the slender beam reached or exceeded the fundamental frequency of the beam. The origin of this phenomenon was that the rigid-flexible coupling effect of the flexible beam caused the total stiffness of the traditional zero-order approximation system model less than zero under the high-speed motion of the base (Feng and Hu, 2002; Duc et al., 2016). Then, with the modern robot technology developing to high-speed, the high-order coupling terms between the large overall motion of the base and the elastic vibration of the flexible manipulator for the modelling of the FCRM should be taken into consideration. On this basis, it is obvious that the base motion has a significant impact on the dynamic response of the FCRM.
Unfortunately, most of the existing literatures mainly analyze the parameters of the FCRM itself on its dynamic character by assuming that the output speed or torque of the driving unit are constant. Qiu (2012) studied the nonlinear vibration control for a Cartesian flexible manipulator by ignoring the torque ripple of the driving motor. The trajectory planning of the base was applied for the vibration control of a ?exible Cartesian robot manipulator with the assumption that the base motion can completely in accordance with the specified manner (Abe, 2011). The eigenspectrum and eigenfunctions of a flexible manipulator are investigated by ignoring the disturbance of the base movement (Coleman and Mcsweeney, 2004). However, considering the coupling fact of the flexible manipulator and the driving base, there are some inadequacies for the dynamic analysis of the flexible manipulator in existing research. Taking the motor drives as an example, it is difficult for the driving base to move with uniform velocity or constant axial force under the influence of power harmonics, motor rotor eccentricity and mechanism gap.
With the fluctuation components of the transmission velocity or axial force of the driving base considered, the FCRM may appear large amplitude vibration of parametric excitation and present complex dynamical characteristics. Atsumi et al. (2005) analyzed the influence of the transmission error on the servo accuracy and the results indicated that the high frequency components had a significant effect. Liu et al. (2015) derived the analytical solutions of a harmonic-disturbance exciting system for a translational flexible-link manipulator. Accordingly, the influence of base disturbance and other external factors on the dynamic characteristics of the flexible manipulator is obvious. With the development of nonlinear dynamics theory, the vibration mechanism of the flexible manipulator has been deeply understood. Based on Lie group theory, Özkaya and Pakdemirli (2002) analyzed the vibration stability of an Euler beam with time-varying velocity motion. By multiple scales method, Chen and Yang (2005) investigated the dynamic stability problems of an axially accelerating beam under the conditions of combination and principal parametric resonance. Similarly, with an added mass placed at an intermediate position, Pratiher and Bhowmick (2011) investigated the nonlinear dynamic characteristics of a Cartesian manipulator. De Luca and Siciliano (1991) analyzed the influences of the payload on the first-order mode shapes of the link-1 and link-2 for a planar multilink lightweight robot, but this research did not further study the impact of the payload on the system dynamics stability. Through the summary of the existing literature, it is obvious that the model parameters of the flexible manipulator are closely related to its mode shapes. However, there is very little research about considering the influence of terminal load on the mode shapes, especially the high-order modes, of the flexible manipulator before the dynamics stability analysis, which may cause dynamic instability phenomenon and unstable response. The FCRM is designed to operate the appropriate objects. Due to the structure characteristics, the terminal load, especially the variable load conditions, can significantly affect the mode frequency and vibration mode of the flexible manipulator. Thus, during the nonlinear dynamics and stability analysis of the FCRM, the impact of terminal load should be taken into consideration firstly.
The objectives of this paper are to construct the nonlinear dynamics model of the FCRM with combining the high-order coupling terms between the large overall motion of the driving base and the elastic vibration of the flexible manipulator. Considering the influence of terminal load on the boundary conditions, the mode shapes of the flexible manipulator are deduced. Based on multiple scales method, the impact mechanisms of terminal load and driving base disturbance on the nonlinear dynamics and vibration stability of the FCRM are analyzed. The structure of this paper is organized as follows. The nonlinear dynamic modelling of the FCRM is given in Sect. 2. Section 3 analyzes the influence of terminal load on the mode shapes of the flexible manipulator. Section 4 provides the nonlinear response analysis of the FCRM, including the amplitude-frequency characteristics and parametric vibration stability of the FCRM. Finally, conclusions are drawn in Sect. 5.
Schematic diagram of the FCRM system under large overall motion.
The structural diagram of the FCRM is shown in Fig. 1, where
Only the tensile and bending deformations of the FCRM are considered. Then,
according to deformation geometry condition, one can obtain:
Spreading Eq. (1) by Taylor series and keeping to
O(
The axial force and lateral bending moment of the flexible manipulator can be
expressed as:
Then, the elastic potential energy of the FCRM system can be further shown
as:
Under large overall plane motion, the radius vector after deformation of
Differentiating Eq. (6) with respect to time yield, the speed vector of
Moreover, the radius vector of the terminal operated object can be written
as:
Differentiating Eq. (8) with respect to time yield, the speed vector of the
terminal operated object can be shown as:
Then, the kinetic energy of the FCRM system can be expressed as:
Based on Rayleigh-Ritz method (Huang and Li, 2010), the longitudinal and
transverse vibration of the FCRM can be obtained as:
The second Lagrange equation expression can be written as:
By substituting Eqs. (1) and (10)–(12) into Eq. (13), the coupling dynamic
model of the longitudinal and transverse vibration for the FCRM system can be
deduced as:
One may note that the coupling dynamic model of the FCRM system contains
inertia
Because there is small deformation in the longitudinal vibration compared
with the transverse vibration, the effect caused by the longitudinal
vibration is not taken into consideration. Compared with the forced
vibration, the parametric excitation is more complex and has a greater
impact on the FCRM system (Yan et al., 2012). Thus, the parametric
excitation, caused by the axial movement disturbance of the base, of the
FCRM system is investigated. To the FCRM system, the excitations of the
unstable motion of the base, caused by the power harmonics, motor rotor
eccentricity, mechanism gap and so on of the motor driving system, have
obvious periodicity. Considering the periodic disturbance signals can be
transformed into trigonometric functions by Fourier transform, for
simplifying the analysis, the disturbance acceleration of the axial
excitation is defined as
The following dimensionless parameters are introduced for simplifying the
coupling dynamic model:
Through substituting Eq. (16) into Eq. (15) and introducing the linear
viscous damping and a book keeping parameter
Contrast curves of the first three order mode shapes for the flexible manipulator with terminal load.
From Eq. (17), we can obtain that the associated parameters of the dynamic
model for the transverse vibration of the FCRM system are mainly determined
by the mode shapes of the flexible manipulator. Thus, the mode shapes under
the influence of terminal load should be analyzed firstly (Chu et al., 2013).
The structural parameters of the FCRM system are defined as:
Because the deflection and the angle of the fixed end for the flexible
manipulator are equal to zero, the fixed end boundary conditions of the
flexible manipulator can be expressed as:
Considering the influence of terminal load, the bending moment and shear
force of the free end for the flexible manipulator are equal to zero and the
inertia force of terminal load, respectively. Then, the free end boundary
conditions of the flexible manipulator can be specified as:
Substituting Eq. (12) into Eqs. (19) and (20), the mode shapes of the
flexible manipulator satisfy the following expressions which are shown as:
Substitution of Eq. (21) into Eq. (18) yields:
Contrast curves of the first three order mode shapes for the flexible manipulator with terminal load by FEM.
By substituting Eq. (22) into Eq. (18), the result combined with Eq. (23)
yields:
In order to guarantee the existence of nontrivial solution for
The relationship between
If
Equations (27) and (28) are the frequency equations and mode shapes of the flexible manipulator without terminal load which is consistent with the expressions of a standard cantilever beam concluded by Singiresu (2004). Then, the validity of Eqs. (25) and (26) are verified.
Mode frequencies of the first three order modes with terminal load by FEM.
Figure 2 shows the first three order mode shapes of the flexible manipulator with terminal load. It is obtained that the terminal load has an apparent effect on the mode shapes and there may cause error if the terminal load is ignored. With the increase of terminal load, it is seen that the amplitude of the first-order mode shape decreases obviously, and the second node of the second-order mode shape as well as the third node of the third-order mode shape move away with the fixed end of the flexible manipulator. Additionally, because the second node of the second mode shape and the third node of the third mode shape are already close to the free end of the flexible manipulator, the influence degree of terminal load on the second mode shape and third mode shape of the flexible manipulator will decrease with the increase of terminal load. Then, the relative vibration amplitude near the free end of the flexible manipulator reduces observably. In order to further verify the correctness of the deduced mode shapes functions, the modal analysis of the flexible manipulator with terminal load is implemented in ANSYS. Additionally, the analysis results by finite element method (FEM) in ANSYS are exported to describe more clearly and the results are shown in Fig. 3. Through comparing Figs. 2 and 3, it is obvious that the variation tendency of the first three mode shapes, obtained by theoretical analysis, is consistent with the results by ANSYS. Thus, the correctness of the deduced mode shape functions is verified.
Table 1 shows the mode frequencies of the first three order modes for the flexible manipulator with terminal load by FEM. It is known that the increase of terminal load reduces the mode frequencies of the flexible manipulator. Owing to the relationship between the mode frequency and the mode shape is same as the relationship between eigenvalue and eigenvector, the change of the mode frequency inevitably lead to the change of the mode shape. Moreover, by comparing the mode frequencies and mode shapes of the first two order modes, it is seen that the closer the node to the fixed end, the higher the mode frequency of the flexible manipulator. From this relationship, we can conclude that the existence of terminal load causes the decrease of the mode frequencies, and then leads to the change of the nodes locations of the mode shapes for the flexible manipulator. Thus, in the nonlinear response analysis of the FCRM system, the influence of terminal load on the mode shapes should be taken into consideration.
Because Eq. (17) contains multiple nonlinear terms, the accurate analytical
solution is hard to solve. Thus, based on multiple scales method, the
approximate analytical solution for Eq. (17) is obtained which can be written
as:
Substituting Eq. (29) into Eq. (17) and equating the coefficients of
The general solution of Eq. (30) can be represented as:
Steady-state response curves of the principal parameter resonance
for the first-order mode of the FCRM with different terminal load:
Steady-state response curves of the principal parameter resonance
for the second-order mode of the FCRM with different terminal load:
Substituting Eq. (32) into Eq. (31) yields:
In order to analyze the principal parametric resonance of the
Substituting Eq. (34) into Eq. (33) and eliminating the secular terms leads
to:
Defining the form of
Through substituting Eq. (36) into Eq. (35) and separating of real and
imaginary parts, one may obtain:
Equation (37) is the average equation of the original system, which can fully
reflect the system dynamics behavior near the balance point. For steady state
response (
Equation (38) is a resonance curve equation of the response amplitude
(
Steady-state response curves of the principal parameter resonance
for the first two order modes of the FCRM with different base excitation:
The stability boundaries for the principal parametric resonance of
the first-order mode of the FCRM with different terminal load in
The stability boundaries for the principal parametric resonance of
the second-order mode of the FCRM with different terminal load in
Figures 4 and 5 show the steady-state response curves of the principal parameter resonance for the first two order modes of the FCRM with different terminal load, where the solid lines indicate the amplitude-frequency response curves of the FCRM under the impact of different terminal load and the dash-dotted lines show the spine lines of the corresponding amplitude-frequency response curves in Figs. 4a and 5a. It is obtained that the amplitude-frequency characteristics of the first two order modes are hard character and soft character respectively, which is identical with the results concluded by Anderson et al. (1996). Then, the correctness of the constructed dynamic model and nonlinear response analysis of the FCRM is verified.
It is seen from Figs. 4a and 5a that the increase of terminal load causes the increase of the bending degree and the resonance field of the amplitude-frequency characteristic curve for the principal parameter resonance of the FCRM system. To the first-order mode, with the increase of terminal load, the bending degree to the right of the spine line becomes larger, which means the hard character of the amplitude frequency characteristic of the first-order mode becomes stronger, as well as the resonance field of the principal parameter resonance for the first-order mode of the FCRM becomes wider. To the second-order mode, with the increase of terminal load, the soft character of the amplitude frequency characteristic becomes stronger and the resonance field becomes wider too.
Figures 4b and 5b show the relationship between the axial excitation
amplitude of the driving base and the response amplitude of the first two
order modes of the FCRM. It is obtained that the excitation acceleration
value of saddle-node bifurcation becomes smaller with the increase of
terminal load for both the first two order modes of the FCRM
(
Figure 6 shows the steady-state response curves of the principal parameter resonance for the first two order modes of the FCRM under different axial excitation of the driving base. It is seen that the spine line has no changes with the increase of the axial excitation acceleration, which means that the axial excitation of the driving base has no influence on the hard or soft character of the amplitude frequency characteristic of the first two order modes. However, with the increase of the axial excitation acceleration, the response amplitude and the resonance field of the principal parameter resonance becomes lager, which indicates that the instability of the FCRM is intensified dramatically.
In order to intuitively express the influence of terminal load on the parametric vibration stability of the FCRM system, the stability boundaries for the principal parametric resonance of the first two order modes of the FCRM are analyzed.
Through defining
The trivial solution of Eq. (39) corresponds to the periodic solution of
Eq. (35). It is obvious that the equilibrium point of Eq. (39) is the origin
of coordinates
The characteristic equation of Eq. (40) is:
The existence of the positive real parts for the eigenvalues of the Jacobin
matrix leads to the instability of the trivial solution of Eq. (39).
According to the Routh-Hurwitz criterion, the sufficient and necessary
conditions for all the roots of Eq. (41) have negative real parts are that
all the Routh-Hurwitz determinants are greater than zero, which is shown as:
Thus, the condition for the stability of the trivial solution of Eq. (39) is:
The linear viscous damping of the first two order modes is defined as 0.1 and
the range of the detuning parameter
This paper mainly investigates the influence of base disturbance and terminal load on the nonlinear response and dynamic stability of the FCRM system. The nonlinear dynamic model of the transverse vibration for the FCRM system, considering the high-order rigid-flexible coupling terms between the large overall motion of the driving base and the elastic vibration of the flexible manipulator, is constructed. Based on the established nonlinear dynamic model and multiple scales method, the amplitude-frequency characteristics and parametric vibration stability of the FCRM are studied. The results demonstrate that the amplitude of the first-order mode shape decreases obviously, and the second node of the second-order mode shape as well as the third node of the third-order mode shape move away with the fixed end of the flexible manipulator, with the increase of terminal load. Furthermore, the excitation acceleration value of saddle-node bifurcation becomes smaller with the increase of terminal load for both the first two order modes of the FCRM. And the instability of the FCRM is intensified with the increase of the base excitation acceleration and terminal load, which is revealed in the increase of the resonance fields and instability areas of the principal parameter resonance of the first two order modes for the FCRM. Thus, the influence of base disturbance and terminal load should be taken into consideration during the kinetic analysis of the FCRM.
All datasets used in the manuscript can be requested from the corresponding author Wei Li (cmeecumt512@yahoo.com).
WL and JYJ put forward the initial ideas; YQW guided the use of multiple scales method; MB Fan and JYJ developed the model code and performed the simulations; XFY contributed the analysis tools; JYJ prepared the manuscript with contributions from all co-authors.
The authors declare that they have no conflict of interest.
This work was partially supported by the National Natural Science Foundation of China (no. U1610111 and no. 51307172), the Six Talent Peaks Project in Jiangsu Province (no. ZBZZ-041), the Postgraduate Cultivation and Innovation Project of Jiangsu Province (KYZZ16-0213) and the Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD). Edited by: Ali Konuralp Reviewed by: two anonymous referees