Next Article in Journal
A Learning Analytics Framework to Analyze Corporal Postures in Students Presentations
Previous Article in Journal
Two-Capacitor Direct Interface Circuit for Resistive Sensor Measurements
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Study of Dynamic Modeling and Multivariable Feedback Control for Flexible Manipulators with Friction Effect and Terminal Load

School of Automation, Northwestern Polytechnical University, Xi’an 710072, China
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(4), 1522; https://doi.org/10.3390/s21041522
Submission received: 28 November 2020 / Revised: 4 February 2021 / Accepted: 19 February 2021 / Published: 22 February 2021
(This article belongs to the Section Sensors and Robotics)

Abstract

:
The flexible manipulato is widely used in the aerospace industry and various other special fields. Control accuracy is affected by the flexibility, joint friction, and terminal load. Therefore, this paper establishes a robot dynamics model under the coupling effect of flexibility, friction, and terminal load, and analyzes and studies its control. First of all, taking the structure of the central rigid body, the flexible beam, and load as the research object, the dynamic model of a flexible manipulator with terminal load is established by using the hypothesis mode and the Lagrange method. Based on the balance principle of the force and moment, the friction under the influence of flexibility and load is recalculated, and the dynamic model of the manipulator is further improved. Secondly, the coupled dynamic system is decomposed and the controller is designed by the multivariable feedback controller. Finally, using MATLAB as the simulation platform, the feasibility of dynamic simulation is verified through simulation comparison. The results show that the vibration amplitude can be reduced with the increase of friction coefficient. As the load increases, the vibration can increase further. The trajectory tracking and vibration suppression of the manipulator are effective under the control method of multi-feedback moment calculation. The research is of great significance to the control of flexible robots under the influence of multiple factors.

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).
r = ( R + x ) c o s θ w ( x , t ) s i n θ ( R + x ) s i n θ + w ( x , t ) c o s θ r ˙ = w ˙ ( x , t ) s i n θ R + x s i n θ + w ˙ ( x , t ) c o s θ θ ˙ w ˙ ( x , t ) cos θ + R + x c o s θ w ˙ ( x , t ) sin θ θ ˙
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).
w ( x , t ) = i = 1 N ϕ i x ε t
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).
E I 4 w x , t x 4 + ρ S 2 w x , t t 2 = F x , t
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.
E I 4 w x , t x 4 + ρ A 2 w x , t t 2 = 0
Using the separate variables method, the bending variables are decomposed spatially and temporally.
E I ϕ x 4 ϕ x x 4 = ρ A τ t 2 τ t t 2
Equation (5) can be equivalent to a constant χ2. The equation is decomposed into Equations (6) and (7).
4 ϕ x x 4 β χ 2 ϕ x = 0
2 τ t t 2 χ 2 τ t = 0
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.
w x , t x = 0 = 0
w x , t x x = 0 = 0
At the free end point, the bending moment and shear force of the beam are obtained under load.
2 w x , t x 2 x = L = 2 ϕ x x 2 x = L = 0
3 w x , t x 3 x = L = 3 ϕ x x 3 x = L = M L 2 w L , t t 2
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.
4 ϕ x x 4 α 2 ϕ x = 0
The general solution of Equation (6) is Equation (13).
ϕ x = C 1 cos α x + C 2 sin α x + C 3 cosh α x + C 4 sinh α x
In the general solution formula, C1, C2, C3, and C4 can be solved by the above equations.
ϕ x / x = α C 1 sin α x + C 2 cos α x + C 3 sinh α x + C 4 cosh α x
2 ϕ x / x 2 = α 2 C 1 cos α x C 2 sin α x + C 3 cosh α x + C 4 sinh α x
3 ϕ x / x 3 = α 3 C 1 sin α x C 2 cos α x + C 3 sinh α x + C 4 cosh α x
Substituting Equations (8)–(11) into Equations (14)–(16), the frequency equation and mode function of the single flexible manipulator can be obtained.
ϕ i x = cosh α i x cos α i x cosh α i L + cos α i L sinh α i L + sin α i L sinh α i x sin α i x
Using the above boundary conditions, the vibration frequency equation can be obtained.
1 + cos α i l cosh α i l = M L / ρ S L α i l sin α i l cosh α i l cos α i l sinh α i l
where αi(i = l, 2, …) is the eigenvalue of the system, and its corresponding natural frequency of vibration is Equation (19).
β = α l 2 E I / ρ l 4
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.
w x , t = i = 1 ϕ x τ t
The kinetic energy of the flexible manipulator is obtained by using the Lagrange method.
T = 1 2 J θ ˙ 2 + 1 2 m ¯ i = 1 N 0 L ϕ i 2 x τ i t d t 2 d x + 1 2 m ¯ θ ˙ i = 1 N 0 L ϕ i 2 x τ i t d t 2 d x + m ¯ i = 1 N 0 L R + x ϕ i 2 x τ i t d t 2 d x + m p L θ ˙ 2 + 1 + θ ˙ 2 + 2 L θ ˙ i = 1 N ϕ i L τ i t 2
where m ¯ = M / l is mass per unit length.
The potential energy of the flexible manipulator is obtained.
V = i = 1 n 0 L E I 2 ϕ i x x 2 2 q i 2 d x + 1 2 m g L cos θ + m p g L cos θ
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.
L = T V d d t L q ˙ i + L q i = Q i
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.
m 11 m 12 m 21 m 22 θ ¨ q ¨ 1 q ¨ n + C 0 θ , θ , q , q ˙ C 1 θ , θ , q , q ˙ C N θ , θ , q , q ˙ + 0 K 1 q 1 K N q N + G 0 0 = Q
where
m 11 = J + m L L 2 + i = 1 N ϕ i 2 L q i 2 + m ¯ i = 1 N 0 L ϕ i 2 x q i 2 d x m 12 = m 21 T = [ γ 1 , γ 2 , γ N ] γ n = m ¯ 0 L R + x ϕ i x d x + M L L ϕ i L n = 1 , 2 N m 22 = d i a g m 1 , m 2 , , m N m i = m ¯ 0 L ϕ i 2 x d x + m ϕ i 2 L
C 0 = 2 θ ˙ n = 1 N m ¯ q ˙ n q n 0 L ϕ n 2 x d x + 2 θ ˙ n = 1 N m q ˙ n q n ϕ n 2 L C n = m ¯ 0 L ϕ n 2 x θ ˙ 2 q n d x + m ϕ n 2 L θ ˙ 2 q n K n = E I 0 L d 2 ϕ i x d x 2 2 d x q n G = 1 2 m g L cos θ + m p g L cos θ Q = τ 0 0
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.
x ¨ = R + x sin θ + w ˙ x , t cos θ θ ¨ R + x sin θ + w ˙ x , t cos θ θ ˙ 2 w ˙ x , t cos θ + w ¨ x , t cos θ θ ˙ w ¨ sin θ y ¨ = R + x cos θ + w ˙ x , t sin θ θ ¨ R + x sin θ + w ˙ x , t cos θ θ ˙ 2 w ˙ x , t cos θ + w ¨ x , t cos θ θ ˙ + w ¨ cos θ
As x = L is substituted into Equation (27), the acceleration vector of the terminal load x ¨ L y ¨ L 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.
F x = ρ 1 S 1 R L x ¨ d x , F y = ρ 1 S 1 R L y ¨ d x
The force on the terminal load can be further calculated.
F M L x = M L x ¨ L , F M L y = M L y ¨ L
The contact forces for the joints of the flexible robot are obtained.
T x = ρ 1 S 1 0 L x ¨ d x + M L x ¨ L T y = ρ 1 S 1 0 L y ¨ d x + M L y ¨ L
According to the Coulomb friction model, the friction moment of the joint is obtained.
τ f = τ f x + τ f y = μ T x R s g n θ ˙ + μ T y R s g n θ ˙
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.
m 11 m 12 m 21 m 22 θ ¨ q ¨ 1 q ¨ n + C 0 θ , θ , q , q ˙ C 1 θ , θ , q , q ˙ C N θ , θ , q , q ˙ + 0 K 1 q 1 K N q N + G 0 0 + τ f 0 0 = Q
For convenience, Equation (32) is shortened to Equation (33).
M θ ¨ q ¨ + C θ , θ , q , q ˙ + K q + G + τ f = Q τ f = μ T x R s g n θ ˙ + μ T y R s g n θ ˙

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.
m 11 θ ¨ + m 12 q ¨ + C 11 θ ˙ + C 12 q ˙ + G θ + τ f θ ˙ , q ˙ = τ m 21 θ ¨ + m 22 q ¨ + C 21 θ ˙ + C 22 q ˙ + K q = 0 τ f = μ T x R s g n θ ˙ + μ T y R s g n θ ˙
The second term in Equation (34) is substituted into the first term, and Equation (34) is further rewritten into Equation (35).
M A θ ¨ + C A θ ˙ + G q + M B q ¨ + C B q ˙ + K q + τ θ ˙ , q ˙ = τ
where M A = m 11 + m 21 , M B = m 12 + m 22 , C A = C 11 + C 21 , C B = C 12 + C 22 As the expected trajectories θ d , θ ˙ d , θ ¨ d T and q d , q ˙ d , q ¨ d T 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.
M A θ ¨ θ ¨ d + C A θ ˙ θ ˙ d + M B q ¨ q ¨ d + C B q ˙ q ˙ d + K q + G θ + τ f θ ˙ , q ˙ = τ
Let lim t ε d , ε ˙ d , ε ¨ d tends to δ , δ ˙ , δ ¨ , δ , δ ˙ and δ ¨ is minimal or 0.
M A θ ¨ θ ¨ d + C A θ ˙ θ ˙ d + M B ε ¨ δ ¨ + C B ε ˙ δ ˙ + K q + G θ + τ f θ ˙ , q ˙ = τ
According to the motion relation of unit mass system, the control closed-loop system satisfies the following conditions.
θ ¨ = θ ¨ d k v θ ˙ θ ˙ d k p θ θ d q ¨ = δ ¨ k q v q ˙ δ ˙ k q p q δ
Substituting Equation (37) into Equation (39), the joint driving torque is obtained.
τ = M A θ ¨ d + C A M A k v θ ˙ θ ˙ d k p θ θ d + M B δ ¨ + C B M B k q v q ˙ δ ˙ k q p q δ + G θ + τ f θ ˙ , q ˙
The control law is expressed in Equation (40).
τ = τ ^ + τ e q + K q q + G θ + τ f θ ˙ , q ˙ .
where τ ^ = M A θ ¨ d + C A M A k v θ ˙ θ ˙ d k p θ θ d τ e q = M B q ¨ d + C B M B k q v q ˙ δ ˙ k q p q δ
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.
θ ¨ θ ¨ d k v θ ˙ θ ˙ d k p θ θ d k q v q ˙ δ ˙ k q p q δ = 0
Finally, the stability of the control system is proved.

5. Results

5.1. Characteristic Analysis of Friction and Terminal Load of the Flexible Manipulator

To verify the feasibility of the coupling dynamics, a series of simulations were carried out on a single degree of freedom flexible robot. The radius of the central rigid body is R = 0.02. The parameters of the manipulator are shown in Table 1.
The driving moment of the joint is 80 N × m. The initial value of the generalized coordinates is θ 0 θ ˙ 0 q q ˙ T = 0 0 0 0 T . The dynamics of flexible robots with different loads and friction coefficients are studied. As the friction coefficient is 0, the terminal load mass is 1 Kg, 2 Kg, and 3 Kg respectively. As the terminal load mass is 1 Kg, the friction coefficient at the joint is 0.01.0.02 and 0.03, respectively. The joint angles and vibration modes of the flexible manipulator were compared under different parameters. Figure 4, Figure 5 and Figure 6 respectively show the joint angle and vibration mode diagrams under different loads.
Figure 7, Figure 8 and Figure 9 respectively show the joint angle and vibration modal diagrams with different friction coefficients.
In Figure 4, under the same torque, when the load increases, the joint angle of the flexible manipulator becomes smaller. Figure 5 and Figure 6 show that when the load increases, the frequency of the modal of the first order and the modal of the second order does not change, and the vibration amplitude increases. As shown in Figure 7, when the load mass is ML = 1 Kg and the friction coefficient is small (the friction coefficient is 0.01), the joint angle is mainly affected by the opposite force of elastic vibration and end-load. This also reflects the low friction coefficient of energy consumption. The variation trend of the joint angle under the influence of flexibility is basically the same as that of a rigid joint angle. In Figure 7, the joint angles of the flexible manipulator under the influence of different friction coefficients are further compared, clearly showing that the joint angle generally decreases with the increase of friction coefficient, and the vibration amplitude increases due to the presence of flexibility. Figure 8 and Figure 9 show that the frequency of the first and second vibration modes is reduced, and as the friction coefficient is 0.01, they present a periodized motion. The amplitude of the modal decreases with the increase of the friction coefficient. This is because friction consumes the energy of the driving force. The influence of the flexible characteristics on friction is taken into account and the original periodic changes are superimposed and offset. This makes the vibration an aperiodic motion.

5.2. Simulation of Control Method

To verify the effectiveness of the control method, the Runge–Kutta method was used for simulation verification based on Matlab simulation platform. The mechanism parameters of the flexible manipulator are shown in Table 1. The terminal mass of the flexible manipulator is 2 Kg. The friction coefficient is 0.03. The desired objective function of the joint of the flexible manipulator is Equation (42). The generalized coordinate initial variable of the flexible manipulator is θ 0 θ ˙ 0 q q ˙ T = 0 0 0 0 T . The control parameters of the flexible manipulator are kv = 20, kp = 8, Kqv = [300 300]T, Kqp = [300 300]T. The joint angles and vibration modes of the flexible manipulator were compared with the rigid feedback and multiple feedback modes. Figure 10 and Figure 11 show the joint angle and velocity of the flexible robot.
Figure 12 and Figure 13 show the first two modes of the flexible robot.
Figure 14 shows the end position of the flexible robot.
To further show the effectiveness of the control method, Power spectral density PSD of the uncontrolled and controlled model was further analyzed and studied. The second derivative and derivative of the given desired trajectory are substituted into Equation (32), and the uncontrolled moment is obtained. Using the periodic diagram power spectral density estimation method, the spectral density of the uncontrolled vibration mode variable is obtained. It is compared with the spectral density of vibration modal variables under the control model in Figure 15 and Figure 16.
θ d = 0.2 × 6 t 5 / T 5 15 t 4 / T 4 + 10 t 3 / T 3 t < T , T = 1           0.2          t T , T = 1
In Figure 11, Figure 12, Figure 13 and Figure 14, the red line vibration is mainly caused by the elastic body of the flexible manipulator. In this paper, due to the joint trajectory tracking, the driving torque is time-varying. The effect of time-varying friction force was further considered. This makes the driving torque coupling and time-varying. The vibration of the flexible manipulator is inevitable under the action of time-varying driving moment and flexible characteristics.
Through comparing the changes in the curve in Figure 10, the joint angle variation under multiple variable feedback control can not only track the target quickly but also suppress the vibration. Compared with the feedback of rigid joints, it can still achieve better feedback control. The further analysis of the velocity shows that the variation amplitude of joint angular velocity under multivariable feedback is large in Figure 11, but its change time is relatively short. Based on the comparison and analysis of the curve changes in Figure 12 and Figure 13, the amplitude of modal changes under multiple variable feedback control is large, but it can achieve fast tracking suppression in a very short time. Modal vibration is affected by the coupling of friction, gravity, and flexible characteristics. The modal amplitude is reduced under the feedback control of the rigid joints, but the vibration cannot be completely suppressed for a limited time. Similarly, through the analysis of the curve changes in Figure 14, it can be seen that the vibration at the end position under the multiple feedback controller can be quickly suppressed and the target can be well tracked after 4 s. According to the curve analysis in Figure 15a and Figure 16a, when the low frequency changes, the power density spectrum of the vibration modal displacement changes greatly. This is consistent with the low-frequency modal-based control research that is the main factor of vibration offset. When the frequency increases, the variation of the power density spectrum is small. This is consistent with the fact that the contribution of the higher frequency to the vibration modal offset is small. These trends prove the effectiveness of the control model. In Figure 15 and Figure 16, the PSD of modal variables under the uncontrolled (a) and controlled model (b) are compared. When the system frequency is close to the natural frequency of the system, the PSD reaches the peak point at the frequency. However, the whole flexible manipulator system has not only external damping, but also internal damping. The PSD will decrease gradually under the influence of mixed damping. Compared with the PSD in the non-control model, the PSD in the control model gradually decreases on the whole, and its amplitude fluctuates little. Therefore, the effectiveness of the control model is further proved by this conclusion.

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.

Author Contributions

Conceptualization, Z.Y.; writing—original draft preparation, F.Z. Both authors have read and agreed to the published version of the manuscript.

Funding

This study is supported by National Natural Science Foundation of China (Grant No. 51505381, 51275423).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects in-volved in the study.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

As the modal interception order N of the flexible manipulator is generally 2, this paper takes N = 2 as an example for derivation. According to Equation (18), the following equations are further expressed.
A i = 0 L ϕ i 2 x d x , B i = 0 L R + x ϕ i x d x , E i = ϕ i L , η i = E I 0 L d 2 ϕ i x / d x 2 2 d x (i = 1,…,N)
when N = 2.
M = m 11 m 12 m 13 m 21 m 22 m 32 m 31 m 32 m 33 m 11 = J + m p L 2 + m / L A 1 q 1 2 + A 2 q 2 2 + m p E 1 2 q 1 2 + E 2 2 q 2 2
m 22 = m / L A 1 + m p E 1 2 m 33 = m / L A 2 + m p E 2 2 m 12 = m 21 = m / L B 1 + m p L E 1
m 13 = m 31 = m / L B 2 + m p L E 2 m 23 = m 32 = 0
C 0 = 2 θ ˙ m / L A 1 + m p E 1 2 q ˙ 1 + m / L A 2 + m p E 2 2 q ˙ 2 C 1 = m / L A 1 + m p E 1 2 θ ˙ 2 q 1
C 2 = m / L A 2 + m p E 2 2 θ ˙ 2 q 2 C = C 0 C 1 C 2 C 3 T
G = m g L cos θ / 2 + m p g L cos θ K = d i a g 0 η 1 2 q 1 η 2 2 q 2
A 1 = 0 L ϕ 1 2 x d x A 2 = 0 L ϕ 2 2 x d x B 1 = 0 L R + x ϕ 1 x d x B 2 = 0 L R + x ϕ 2 x d x
E 1 = ϕ 1 L E 2 = ϕ 2 L η 1 = E I 0 L d 2 ϕ 1 x / d x 2 2 d x η 2 = E I 0 L d 2 ϕ 2 x / d x 2 2 d x
where ϕ(x) is a shape function, its specific equation is given as Equation (17).

References

  1. Flores-Abad, A.; Ma, O.; Pham, K.; Ulrich, S. A review of space robotics technologies for on-orbit servicing. Prog. Aerosp. Sci. 2014, 68, 1–26. [Google Scholar] [CrossRef] [Green Version]
  2. Patane, S.; Joyce, E.R.; Snyder, M.P.; Shestople, P. In-Space Manufacturing and Assembly for Next-Generation Space Habitats. In Proceedings of the AIAA SPACE and Astronautics Forum and Exposition, Orlando, FL, USA, 12–14 September 2017; pp. 1–10. [Google Scholar]
  3. Kaouk, J.H.; Haber, G.P.; Autorino, R. A novel robotic system for single-port Urologic Surgery: First clinical investigation. Eur. Urol. 2014, 66, 1033–1043. [Google Scholar] [CrossRef] [PubMed]
  4. Xu, P.; Cheung, C.F.; Wang, C.; Zhao, C. Novel hybrid robot and its processes for precision polishing of freeform surfaces. Precis. Eng. 2020, 64, 53–62. [Google Scholar] [CrossRef]
  5. Cao, H.J.; Zhou, J.; Pei, J.; Hon, K.K.B. An integrated processing energy modeling and optimization of automated robotic polishing system. Robot. Comput. Integr. Manuf. 2020, 65, 101973. [Google Scholar] [CrossRef]
  6. Xin, X. Swing-up Control for a Two-link Underactuated Robot with a Flexible Elbow Joint: New Results Beyond the Passive Elbow Joint. In Proceedings of the 50th IEEE Conference on Decision and Control and European Control Conference, Orlando, FL, USA, 12–15 December 2011; pp. 2481–2486. [Google Scholar]
  7. Spector, V.A.; Flashner, H. Modeling and design implications of non-collocated control in flexible systems. J. Dyn. Sys. Meas. Control. 1990, 112, 186–193. [Google Scholar] [CrossRef]
  8. Zhang, N.; Feng, Y.; Yu, X. Optimization of Terminal Sliding Control for Two-link Flexible Manipulators. In Proceedings of the 30th Annual Conference of IEEE Industrial Electronics Society, Busan, Korea, 2–6 November 2004; pp. 1318–1322. [Google Scholar]
  9. Wang, F.; Gao, Y. Advanced Studies of Flexible Robotic Manipulators Modeling, Design, Control and Applications. In Series in Intelligent Control and Intelligent Automation; World Scientific: Tuck Link, Singapore, 2003; Volume 4. [Google Scholar]
  10. Fukuda, T.; Arakawa, A. Modeling and control characteristics for a two-degrees-of-freedom coupling system of flexible robotic arms. JSME Int. J. 1987, 30, 1458–1464. [Google Scholar] [CrossRef]
  11. Buffinton, K.W. Dynamics of Elastic Manipulators with Prismatic Joints. J. Dyn. Syst. Meas. Control. 1992, 114, 41–49. [Google Scholar] [CrossRef]
  12. Rahimi, H.; Nazemizadeh, M. Dynamic analysis and intelligent control techniques for flexible manipulators: A review. Adv. Robot. 2013, 28, 63–76. [Google Scholar] [CrossRef]
  13. Zhang, C.-Y.; Bai, C.-C. Extremum response surface method of reliability analysis on two-link flexible robot manipulator. J. Cent. South Univ. 2012, 19, 101–107. [Google Scholar] [CrossRef]
  14. Scaglioni, B.; Ferretti, G. Closed form model of manipulators with highly flexible links. IFAC-PapersOnLine 2015, 48, 653–654. [Google Scholar] [CrossRef]
  15. Khairudin, M.; Mohamed, Z.; Husain, A.R.; Ahmad, M.A. Dynamic Modelling and Characterisation of a Two-Link Flexible Robot Manipulator. J. Low Freq. Noise Vib. Act. Control. 2010, 29, 207–219. [Google Scholar] [CrossRef]
  16. Pradhan, S.K.; Subudhi, B. NARMAX Modeling of a Two-link Flexible Robot. In Proceedings of the Annual IEEE India Conference, Hyderabad, India, 16–18 December 2011; pp. 1–5. [Google Scholar]
  17. Khairudin, M.; Mohamed, Z.; Husain, A.; Mamat, R. Dynamic characterisation of a two-link flexible manipulator: Theory and experiments. Adv. Robot. Res. 2014, 1, 61–79. [Google Scholar] [CrossRef] [Green Version]
  18. Li, R.J.; Liu, H.Z.; Li, P.F. Dynamic modeling of planar two-degree-of-freedom five-bar mechanism with joint dissipation considered. J. Agric. Mach. 2008, 39, 157–161. [Google Scholar]
  19. Alinalbu, S.; Gerd, H. A globally stable state feedback controller for flexible joint robots. Adv. Robot. 2001, 15, 799–814. [Google Scholar]
  20. Moreno, J.; Kelly, R. Pose Regulation of Robot Manipulators with Dynamic Friction Compensation. In Proceedings of the 44th IEEE Conference on Decision and Control, Seville, Spain, 15 December 2005; pp. 4368–4372. [Google Scholar]
  21. He, Y.; Mei, J.G.; Fang, Z.W. Study on energy consumption of a single joint servo system for a robot. Manuf. Technol. Mach. Tools 2017, 9, 84–90. [Google Scholar]
  22. Wu, W.X.; Zhu, S.Q.; Wang, X.Y. Based on friction fuzzy modeling and compensating robot low speed control. J. Electr. Mach. Control. 2013, 17, 46–53. [Google Scholar]
  23. Liu, F.C.; Liu, C.; Feng, Q.L. Influence of friction under gravity on motion control of space manipulator. High Technol. Commun. 2014, 24, 65–71. [Google Scholar]
  24. Dou, B.; Yue, X.K. Collision dynamics analysis and control with consideration of friction in space capture process. China Space Sci. Technol. 2018, 38, 54–62. [Google Scholar]
  25. Mahamood, R.M.; Pedro, J.J. Hybrid PD/PID Controller Design for Two-Link Flexible Manipulators. In Proceedings of the 8th Asian Control Conference, Kaohsiung, Taiwan, 15–18 May 2011; pp. 1358–1363. [Google Scholar]
  26. Mohamed, Z.; Khairudin, M.; Husain, A.; Subudhi, B. Linear matrix inequality-based robust proportional derivative control of a two-link flexible manipulator. J. Vib. Control. 2014, 22, 1244–1256. [Google Scholar] [CrossRef]
  27. Walsh, A.; Richard, F.J. Modeling and control of flexible telescoping manipulators. IEEE Trans. Robot. 2015, 31, 936–947. [Google Scholar] [CrossRef]
  28. Matsuno, F.; Yamamoto, K. Dynamic hybrid force/position control of a two degree-of-freedom flexible manipulator. J. Robot. Syst. 2007, 11, 355–366. [Google Scholar] [CrossRef]
  29. Sawada, M.; Itamiya, K. A Position Control of 2 DOF Flexible Link Robot Arms Based on Computed Torque Method. In Proceedings of the 2012 International Conference on Advanced Mechatronic Systems, Tokyo, Japan, 18–21 September 2012; pp. 547–552. [Google Scholar]
  30. Vandini, A.; Salerno, A.; Payne, C.J.; Yang, G. Vision-based Motion Control of a Flexible Robot for Surgical Applications. In Proceedings of the 2014 IEEE International Conference on Robotics and Automation, Hong Kong, China, 31 May–7 June 2014; pp. 6205–6211. [Google Scholar]
  31. Chu, Z.; Cui, J. Experiment on vibration control of a two-link flexible manipulator using an input shaper and adaptive positive position feedback. Adv. Mech. Eng. 2015, 7, 1–13. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The deformation diagram of the central rigid body rotating flexible beam.
Figure 1. The deformation diagram of the central rigid body rotating flexible beam.
Sensors 21 01522 g001
Figure 2. The friction model of flexible manipulator.
Figure 2. The friction model of flexible manipulator.
Sensors 21 01522 g002
Figure 3. Design of flexible manipulator controller with multiple closed loops.
Figure 3. Design of flexible manipulator controller with multiple closed loops.
Sensors 21 01522 g003
Figure 4. Joint Angle changes of the manipulator under different loads.
Figure 4. Joint Angle changes of the manipulator under different loads.
Sensors 21 01522 g004
Figure 5. The first order modal changes of the manipulator under different loads, (a) 30 kg, (b) 20 Kg, (c) 10 Kg.
Figure 5. The first order modal changes of the manipulator under different loads, (a) 30 kg, (b) 20 Kg, (c) 10 Kg.
Sensors 21 01522 g005
Figure 6. The second order modal changes of the manipulator under different loads, (a) 30 kg, (b) 20 Kg, (c) 10 Kg.
Figure 6. The second order modal changes of the manipulator under different loads, (a) 30 kg, (b) 20 Kg, (c) 10 Kg.
Sensors 21 01522 g006
Figure 7. Joint angle changes of the manipulator under different friction coefficients.
Figure 7. Joint angle changes of the manipulator under different friction coefficients.
Sensors 21 01522 g007
Figure 8. The first-order modal changes of the manipulator with different friction coefficients, (a) 0.01, (b) 0.02, (c) 0.03.
Figure 8. The first-order modal changes of the manipulator with different friction coefficients, (a) 0.01, (b) 0.02, (c) 0.03.
Sensors 21 01522 g008
Figure 9. The second order modal changes of the manipulator with different friction coefficients, (a) 0.01, (b) 0.02, (c) 0.03.
Figure 9. The second order modal changes of the manipulator with different friction coefficients, (a) 0.01, (b) 0.02, (c) 0.03.
Sensors 21 01522 g009
Figure 10. Joint angles variation of different controls.
Figure 10. Joint angles variation of different controls.
Sensors 21 01522 g010
Figure 11. Joint angular velocity of different controls.
Figure 11. Joint angular velocity of different controls.
Sensors 21 01522 g011
Figure 12. The first-order mode of different controls.
Figure 12. The first-order mode of different controls.
Sensors 21 01522 g012
Figure 13. The second-order mode of different controls.
Figure 13. The second-order mode of different controls.
Sensors 21 01522 g013
Figure 14. The terminal position change of flexible manipulator.
Figure 14. The terminal position change of flexible manipulator.
Sensors 21 01522 g014
Figure 15. PSD of first-order modal variable of different controls, (a) controlled, (b) uncontrolled.
Figure 15. PSD of first-order modal variable of different controls, (a) controlled, (b) uncontrolled.
Sensors 21 01522 g015
Figure 16. PSD of second-order modal variables of different controls, (a) controlled, (b) uncontrolled.
Figure 16. PSD of second-order modal variables of different controls, (a) controlled, (b) uncontrolled.
Sensors 21 01522 g016
Table 1. The parameters of the flexible robot mechanism.
Table 1. The parameters of the flexible robot mechanism.
Length
L (m)
Density
ρ (Kg/m)
Width
b (m)
Height
h (m)
Rotational Inertia
J (Kg × m2)
Elasticity Modulus
EI (N × m)
Tip Mass
MTip (Kg)
55.92920.060.045.126.21 × 107--
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, F.; Yuan, Z. The Study of Dynamic Modeling and Multivariable Feedback Control for Flexible Manipulators with Friction Effect and Terminal Load. Sensors 2021, 21, 1522. https://doi.org/10.3390/s21041522

AMA Style

Zhang F, Yuan Z. The Study of Dynamic Modeling and Multivariable Feedback Control for Flexible Manipulators with Friction Effect and Terminal Load. Sensors. 2021; 21(4):1522. https://doi.org/10.3390/s21041522

Chicago/Turabian Style

Zhang, Fuli, and Zhaohui Yuan. 2021. "The Study of Dynamic Modeling and Multivariable Feedback Control for Flexible Manipulators with Friction Effect and Terminal Load" Sensors 21, no. 4: 1522. https://doi.org/10.3390/s21041522

APA Style

Zhang, F., & Yuan, Z. (2021). The Study of Dynamic Modeling and Multivariable Feedback Control for Flexible Manipulators with Friction Effect and Terminal Load. Sensors, 21(4), 1522. https://doi.org/10.3390/s21041522

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop