1. Introduction
A flexible robot is characteristically lightweight, with high flexibility and a large load [
1]. Through the structure design, the robot itself can be further optimized by considering the flexibility factor [
2]. Since the robot design needs to meet the actual working conditions, designing its control method poses additional difficulties [
3]. In practice, the stiffness change is absolute; that is, the elastic deformation is absolute. The dynamics of flexible robots are considered based on the comprehensive factors of deformation and control [
4]. Therefore, dynamic modeling and characteristic analysis of flexible robots are essential.
Flexible manipulators have incomparable advantages over rigid manipulators, so current research is more inclined to focus on flexible manipulators [
5]. But flexible manipulators have many limitations, including control complexity, uncertainty, MIMO and nonlinear systems, control overflow, and observed overflow [
6]. The main reason for the complexity is the selection of dynamic models [
7]. In the past three decades, people have developed many discrete methods for flexible robots. These methods include the finite element method (FEM), assumed mode method (AMM), and lumped parameter method [
8]. Dynamics modeling methods are usually established based on AMM and Lagrangian methods. These methods have the advantages of high computational efficiency and flexible selection of boundary conditions [
9]. Fukuda et al. established the dynamic characteristics of the two-link flexible manipulator under the influence of payload and gravity. He only analyzed and studied the effect of bending vibration on the manipulator [
10]. Buffinton et al. [
11] took the translation member of the standard arm as the elastic beam and established the motion equation containing the translation movement of the elastic member of the standard arm. Using the finite variable mode amplitude and space characteristic mode function, he further established the dynamic equation. Rahimi and Nazemizadeh [
12] present a new method of selecting linkage boundary conditions and their modal eigenfunctions. They further used a dynamic model to discuss its deflection in terms of configuration vibration modes. Yi and Chen [
13] presented dynamic modeling of the combined payload using assumed mode method (AMM). Moreover, he conducted numerical and experimental verification. Scaglioni and Ferretti [
14] proposed a closed type dynamics model of a highly flexible robot based on Neutons-Euler. They used Matlab/Simulink, Modelica/Dymola, MSC/Adams and other software to verify the model. Khairudin et al. [
15] built a two-link flexible robot to carry out the dynamic model based on comprehensively considering the influence of inertia and load on the damping hub of the structure, and they verified the effectiveness of the model through numerical simulation. Pradhan and Subudhi [
16] used recursive extended Least squares (RELS) to achieve real-time identification of the NARMAX model. The model was verified by considering different loads and the experimental results were compared with the NARMAX model. Khairudin et al. [
17] proposed a dynamic modeling and characterization of plane two-link flexible manipulator TLFM with payload. Using the Euler–Lagrange model and the AMM model, they established the dynamic model of the manipulator, which was verified by numerical simulation and experiment. Based on the above studies, the Lagrangian method was used to study the dynamic modeling of flexible robots from the perspective of its boundary conditions and load characteristics. In this paper, the boundary between the central rigid body and the load characteristics is considered comprehensively, and the dynamics are further constructed.
The energy source of robot movement mainly depends on the driving moment of joint, in which the energy loss of the robot must be considered. Many scholars mainly focus on the study of joint friction loss. At present, the most common friction models are the Stribeck friction model, the coulomb-Viscous friction model, the mane model, and the coulomb friction model [
18]. Alinalbu et al. [
19] used both coulomb friction and static friction models to describe the friction phenomenon of joint motion, and established a dynamic model considering both joint flexibility and friction. They analyzed the influence of friction on the dynamic characteristics of flexible joints. Moreno J et al. [
20] took the two-degree-of-freedom manipulator as the research object and established the dynamic model of the viscous frictional manipulator. They further analyzed the influence of friction at the joints on the dynamic characteristics and end position errors of the manipulator. He et al. [
21] used the LuGre friction model to establish the joint dynamics model containing friction, and analyzed the energy consumption caused by friction at robot joints. Wu et al. [
22] proposed an extended joint friction model and analyzed the joint friction, finding that the low-speed motion performance of the mechanical arm deteriorates, and the operation accuracy is greatly affected. Taking a space manipulator with two degrees of freedom (2-DOF) as the research object, Liu et al. [
23] solved the friction force based on the Coulomb friction model, and they explored and analyzed the influence of joint friction on the driving force of a space manipulator joint. The results show that gravity has a great influence on the dynamic characteristics of the manipulator. By integrating Coulomb’s friction theory with the Lagrangian multiplier algorithm, Dou and Yue [
24] derived a friction model with sliding and viscous states. They analyzed the energy loss caused by friction in the collision system. Therefore, the coulomb friction model is used to solve the friction moment at the joint of the manipulator. The force of friction is proportional to the normal load. It is opposite to the direction of relative motion and does not depend on the friction theory of contact area. This coulomb friction model is relatively intuitive and easy to calculate and apply [
25].
Many control methods of flexible robots have been developed and improved based on classical control methods of traditional rigid robots. The proportional differential/proportional differential integral (PD/PID) controller with position feedback and state feedback control with tracking and calculating torque was modified for flexible manipulator operation. Moreover, further control studies for different payloads were conducted. Mohamed et al. [
26] used a linear matrix inequality PD (LMI PD) controller to control the position and end vibration suppression of the two-link flexible manipulator under the condition of changing the load effect. The robustness of the control technology is verified by the comparison between numerical simulation and experimental results. Walsh et al. [
27] proposed the modeling and control of a scalable, flexible robot. The terminal position control of manipulator using a PD controller was realized. Matsuno and Yamamoto [
28] used singular perturbation technology to divide the dynamic system into a slow subsystem and a fast subsystem. They further adopted state feedback control to control and suppress vibration at the end of the flexible two-link. Sawada and Itamiya [
29] used feedback linearized extension (calculated torque method) to control the position of the flexible robot. In the case of uncertainty and interference, the proposed method was numerically simulated. Vandini et al. [
30] proposed a flexible two-link motion control based on a visual sensor in which feedback signals are generated by visual sensors and PID controllers are used to control the movement of flexible robots. Chu and Cui [
31] propose a method to actively suppress the internal vibration of the two-link flexible manipulator to adapt to the variation of model parameters. The model parameters are composed of input shaper and multi-mode adaptive positive-position feedback.
In this paper, coupling characteristics between joint friction, terminal load, and flexible machinery are considered comprehensively. Taking the structure of the central rigid body-rotating flexible Euler–Bernouli beam-terminal load as the research object, the dynamics of the flexible manipulator are established by using assumed modes method (AMM) and the Euler-Lagrange method. Finally, the dynamics of the flexible manipulator are verified and its characteristics are analyzed by simulation. The vibration characteristics of the manipulator under different loads and friction conditions are analyzed respectively. The coupling dynamics of the flexible robot are further decomposed and the track tracking and vibration suppression are carried out by using multiple feedback control methods.
2. Dynamic Modeling of the Flexible Robot
The structure of the flexible robot is mainly composed of motors, long links, terminal grasping mechanisms, and loads. The concrete structure is simplified as a central rigid body, flexible manipulator, and terminal load. The structure is shown in
Figure 1.
In
Figure 1, the manipulator structure motor drives the rotor, profile, flexible beam, and end load. The drive motor shaft is located at the fixed end with a rotor radius
R. This coordinate system is the inertial coordinate system
O-
XYZ. The connection point between the contour and the beam is the following coordinate system is
o-
xyz. And the rotation between the two frames is
θ, so it is a rigid rotation. Its rotational moment is
τθ. The length of the manipulator is
L, and the density of the mechanical manipulator
ρ. The load mass of the free end is
ML. The position and velocity of any point
p on the flexible robot and the position and velocity of the end load are shown in Equation (1).
Its length L is substituted into Equation (1), and the terminal position and velocity of the flexible manipulator are obtained.
The position vector of flexible deformation of point
p on the flexible manipulator is
w(
x,
t). The
N-dimensional modal functions are selected by the modal truncation method. The displacement vector of the flexible deformation is expressed as Equation (2).
where
ϕi(
x) is the modal function,
εi(
t) is the generalized coordinates.
According to the existing research on the cantilever beam, the deformation of flexible manipulator was analyzed. Based on the flexural elasticity of material mechanics, when L/hz > 10, L/hb > 10, the flexible manipulator is described by the Euler–Bernoulli beam. The axial deformation and shear deformation were not taken into account. The parameters hz and hb were the height and width of the flexible manipulator.
According to the bending vibration principle, the forced bending equation of the flexible manipulator is expressed as Equation (3).
where
E is Modulus of elasticity,
I is Moment of inertia section,
S is Cross-sectional area,
F(
x,
t) is the external force.
To solve its general solution, the free vibration Equation (4) is obtained.
Using the separate variables method, the bending variables are decomposed spatially and temporally.
Equation (5) can be equivalent to a constant χ2. The equation is decomposed into Equations (6) and (7).
where
β is
EI/
ρS. Equation (6) is the variable function of the flexible manipulator space. Equation (7) is the time function of the flexible manipulator.
The solution of Equation (6) requires four boundary conditions. Firstly, four boundary conditions are determined by mechanical analysis of the bar. As the end point of the beam is fixed, the bending displacement and slope of the fixed end of the flexible manipulator are 0.
At the free end point, the bending moment and shear force of the beam are obtained under load.
where
ML is the mass of the end load. Equation (11) is the shear force under load, and the right formula of Equation (11) is the product of the acceleration and mass at the end of the link under bending changes.
The parameter
βχ2 is equivalent to
α4.
The general solution of Equation (6) is Equation (13).
In the general solution formula,
C1,
C2,
C3, and
C4 can be solved by the above equations.
Substituting Equations (8)–(11) into Equations (14)–(16), the frequency equation and mode function of the single flexible manipulator can be obtained.
Using the above boundary conditions, the vibration frequency equation can be obtained.
where
αi(
i = l, 2, …) is the eigenvalue of the system, and its corresponding natural frequency of vibration is Equation (19).
By Equations (18) and (19), each vibration natural frequency
β corresponds to a particular vibration mode function
ϕ(
x) and time function coefficient tau
τ(
t). By using the principle of linear superposition, the solution of Equation (20) can be obtained.
The kinetic energy of the flexible manipulator is obtained by using the Lagrange method.
where
is mass per unit length.
The potential energy of the flexible manipulator is obtained.
In Equation (22), the first term is the energy generated by the deformation of the flexible manipulator. The second and third items are the gravitational potential energy of the flexible manipulator and the end-point load respectively. According to the Lagrange equation of the second kind, the dynamic equation of the flexible manipulator can be obtained.
where
qi = [
θ,
q1,
q2,……
qn]
Tis the generalized coordinates of flexible manipulator,
Qi = [
τθ,0,….0]
Tis the generalized force of the flexible manipulator.
The dynamics of a flexible manipulator is expressed as an equation of state.
where
where the state matrix of
mij (
i,
j = 1, 2) is the mass inertia matrix, which is symmetric and bounded.
C is the coupling matrix of Coriolis force and centrifugal force.
K is the stiffness matrix,
G is the gravity matrix,
Q is the matrix of the input torque. The specific equation of state implementation is given in the
Appendix A.
3. Improved Friction Model
Based on the analysis of the instantaneous static equilibrium principle, the friction torque of articulation at the joint is obtained based on the Coulomb friction force. It is shown in
Figure 2. The improved friction model takes into account the influence of the flexible characteristics of the manipulator and the lumped mass of the end. Therefore, the friction torque at the joint is recalculated [
20].
According to the velocity vector in Equation (1), the acceleration vector at any point on the manipulator can be obtained through further derivation.
As x = L is substituted into Equation (27), the acceleration vector of the terminal load is obtained. The force at any point on the mechanical arm is decomposed and projected on the XOY plane, and Fx and Fy can be obtained. Since the flexible manipulator is studied by distributed parameters, the projection force can be obtained by using the fixed-base integral method and the force balance principle.
According to D’alembert’s principle and Newton’s second law, the forces on each microunit of the manipulator are obtained. The force is integrated over the length range of the manipulator arm [0,
L] and the force components on the mechanical arm members in
x and
y directions are obtained.
The force on the terminal load can be further calculated.
The contact forces for the joints of the flexible robot are obtained.
According to the Coulomb friction model, the friction moment of the joint is obtained.
where
μ is the friction coefficient.
The friction factor takes into account the effect of the flexibility characteristics on the mass of the end load. Combining Equations (24) and (31), a dynamic model with flexibility characteristics and terminal load is established.
For convenience, Equation (32) is shortened to Equation (33).
4. The Controller Design for Multiple Closed-Loop Adaptive Computing Torques
To better control terminal load and the coupling dynamics of flexible and friction characteristics, the dynamics of the flexible manipulator are first decomposed and multiple closed-loop controllers are applied to the manipulator controller. Secondly, considering that the two subsystems as the load at the end of the flexible deformation of its kinetic energy is smaller, we are therefore able to ignore the flexible deformation of the load in the end as well as the influence of the friction torque at the same time. The influence of the flexible and terminal load specific controller design is shown in
Figure 3.
In
Figure 3,
θd is the rigid joint trajectory, and
qd is the specified trajectory of the flexible vibration. However, the flexible trajectory is a virtual quantity. According to the need of the vibration suppression target, the specified trajectory should be as small as possible or guaranteed to be within a small range. According to Equation (32), rigid motion and flexible vibration are interactive. The dynamics module of the flexible manipulator is represented. Based on the rigid feedback of the manipulator, the flexible feedback is designed to track the trajectory of the whole manipulator system. Compared with the previous singular perturbation method and intelligent method, this method is realized by decomposing the subsystem of the coupling system. The difference between them is that the flexible system constrains the control error by the virtual control objective, and it is equivalent in the way of the black box. The realization of the control method can also be identified by the intelligent control, but this will increase the amount of control tasks. According to Equations (27)–(31), both rigid and flexible variables are taken into account for joint friction. Therefore, rigid variables and flexible variables are introduced into the friction compensation design.
As the systems
M,
C,
G, and
K are deterministic systems, and they are known, Equation (33) is further derived.
The second term in Equation (34) is substituted into the first term, and Equation (34) is further rewritten into Equation (35).
where
,
,
,
As the expected trajectories
and
are given, the flexible vibration feedback control is carried out based on rigid motion. Equation (35) will be expressed as the dynamic equation of the error space.
Let
tends to
,
,
and
is minimal or 0.
According to the motion relation of unit mass system, the control closed-loop system satisfies the following conditions.
Substituting Equation (37) into Equation (39), the joint driving torque is obtained.
The control law is expressed in Equation (40).
where
According to the motion equation of the error space above, it can be seen that appropriate
kv,
kp,
kqv, and
kqp are selected. Substituting Equation (40) into Equation (34), Equation (41) can be obtained.
Finally, the stability of the control system is proved.
6. Discussion
Considering the effects of a central rigid body and a curved flexible beam and loads, the dynamics of a complex coupled flexible manipulator were established. For further precise dynamics modeling, the effect of friction on the joints were considered. Compared with the previous friction model, the friction model not only considers the influence of gravity and rigid motion, but also considers the influence of flexible factors. Considering the friction effect, the dynamic model was further modified. The model is based on the combination of flexibility, terminal load, and improved friction effect. To control the flexible manipulator more conveniently and effectively, the coupling dynamics of the manipulator were decomposed. Flexible reaction targets were limited to a very small range and the flexible robot was controlled by a multi-loop feedback torque method. Compared with the previous feedforward control or complex control methods, it was simple and easy to implement. As the effective flexible control target approaches 0, the flexible manipulator can ensure that the control will be in the effective range.
Through the simulation analysis of the experiment, it can be seen that when the friction coefficient of the joint increases, the vibration of the mechanical arm presents an aperiodic vibration, and its vibration frequency decreases, and the amplitude of vibration significantly decreases. This is due to the great loss of driving energy caused by friction, and it is superimposed with the original vibration. As the terminal load increases in a certain range, the amplitude of vibration increases, and the vibration frequency almost has no obvious change. Under the multivariable feedback torque controller, the control method can effectively track the control trajectory, and the vibration of the internal flexible manipulator can be well suppressed within the effective time (t < 4 s). Compared with the rigid joint control, this control method has obvious good control effect.
7. Conclusions
In this paper, the center rigid body–flexible beam–load structure was selected as the research object. Considering the friction effect, the interaction effect between the flexible characteristics, and the terminal load, the dynamic model of the flexible robot with coupling characteristics was established. To realize track tracking and vibration suppression, the multiple closed-loop controllers were designed according to their rigid and flexible effects. The results show that the larger the terminal load is, the smaller the flexible vibration frequency is. but the larger the terminal load is, the larger the amplitude is. As the friction coefficient increases, the amplitude and vibration frequency decreases, and the vibration is aperiodic. Trajectory tracking and vibration suppression control were carried out for the nominal model with multiple feedback control. When the target of the flexible feedback approached the minimum value, and considering the interaction between the friction and the terminal load and the flexible vibration, the joint trajectory tracking and modal vibration suppression can be achieved within 4s. Finally, compared with the conventional PID rigid joint controller, it can be seen that the controller is effective and better.
In this paper, the dynamics of the manipulator with multiple factors such as flexibility, friction and terminal load were studied and analyzed. However, many factors such as joint flexibility, temperature effect, and air environment damping were not taken into account. The study of complex control methods under the influence of multiple factors needs further study. To meet the practical need, the dynamics of precision robots and more complex control research are the next steps.