Control Method of Flexible Manipulator Servo System Based on a Combination of RBF Neural Network and Pole Placement Strategy

: Gravity and ﬂexibility will cause ﬂuctuations of the rotation angle in the servo system for ﬂexible manipulators. The ﬂuctuation will seriously affect the motion accuracy of end-effectors. Therefore, this paper adopts a control method combining the RBF (Radial Basis Function) neural network and pole placement strategy to suppress the rotation angle ﬂuctuations. The RBF neural network is used to identify uncertain items caused by the manipulator’s ﬂexibility and the time-varying characteristics of dynamic parameters. Besides, the pole placement strategy is used to optimize the PD (Proportional Differential) controller’s parameters to improve the response speed and stability. Firstly, a dynamic model of ﬂexible manipulators considering gravity is established based on the assumed mode method and Lagrange’s principle. Then, the system’s control characteristics are analyzed, and the pole placement strategy optimizes the parameters of the PD controllers. Next, the control method based on the RBF neural network is proposed, and the Lyapunov stability theory demonstrates stability. Finally, numerical analysis and control experiments prove the effectiveness of the control method proposed in this paper. The means and standard deviations of rotation angle error are reduced by the control method. The results show that the control method can effectively reduce the rotation angle error and improve motion accuracy.


Introduction
Flexible manipulators are complex systems with multiple inputs and multiple outputs. They are widely used in in-depth space exploration, industrial assembly, and other fields [1,2]. Compared with the traditional rigid manipulators, flexible manipulators have a larger radius of rotation and lighter weight. Flexible manipulators are more likely to vibrate in a low-frequency area. In order to clarify the reasons for the vibration of flexible manipulators, it is necessary to carry out a dynamic modeling analysis. At present, the modeling methods of flexible manipulators are mainly divided into two categories. One is the finite element (FE) method, and the other is the assumed mode method (AMM). In literature [3][4][5], the dynamic equations of flexible manipulators with the rotation angle coupled with modal coordinates are established by using the AMM. In literature [6], the FE method is used to establish the dynamic equation of flexible manipulators with flexible joints. Compared with the AMM, the FE method is more suitable for flexible manipulators with irregular shapes [7]. Nowadays, in the dynamic modeling processes of flexible manipulators, it is usually assumed that flexible manipulators move in the horizontal plane, as shown in the literature [8,9]. According to this assumption, the influence of gravity on flexible manipulators can be ignored. However, in some practical applications of flexible manipulators, the gravity factor needs to be considered [2,10]. In literature [11], the gravity Section 5 carries out the numerical simulation and experiments of flexible manipulators to verify the proposed control strategy's effectiveness. Section 6 states the conclusion.

Dynamic Modeling of the Flexible Manipulator Servo System
The servo system for flexible manipulators can be considered to be composed of a motor side and a flexible load side, as shown in Figure 1. The motor side contains a gear reducer that reduces the speed of the motor side. The flexibility of a flexible load is equivalent to the flexibility of a flexible manipulator. According to the literature [25,26], the flexible load can be equal to the Euler-Bernoulli flexible beam model. where T m , T l , T a represent the electromagnetic motor torque, the external torque of the load side, the input torque of the flexible load; w(x, t) represents the deformation of the flexible load; θ m , θ represent the rotation angle of motor side, and the rotation angle of flexible load; XOY, x 0 Oy 0 represent the static coordinate and the moving coordinate; mg represents the gravity.
If the torsional stiffness is ignored, the rotation angle of the motor side is equal to the rotation angle of the flexible load. Wherein, the relationship between T m , T l , and T a is shown in Equation (1): According to the literature [27], the deformation of flexible loads is considered to be a two-dimensional function of the position of the load section x and time t. In the servo system's modeling processes for flexible manipulators, the above model can be equivalent to the system that the cantilever beam rotates around the center. When the cantilever beam rotates, the transverse deformation of the flexible load is obvious. The longitudinal deformation occurs at the same time as the transverse deformation. According to the literature [28], the longitudinal deformation of the flexible load can be ignored since the longitudinal deformation is small. Therefore, the flexible load in the servo system for flexible manipulators is equivalent to an Euler-Bernoulli beam.

Description of Flexible Load-Deformation
According to the vibration theory of flexible beams, the deformation of the flexible load is described by the AMM. It is assumed that a flexible beam's deformation is related to the modal function and the modal coordinate. The deformation of the flexible load is shown in Equation (2): where φ i (x), δ i (t) represent the modal function and the modal coordinate.
The flexible manipulator's external load is ignored, and the flexible load is equivalent to the flexible cantilever beam. Thus, the flexible load boundary conditions can be obtained, as shown in Equation (3): According to the vibration theory, the partial differential equation of the flexible load's lateral vibration can be obtained, as shown in Equation (4): where ρ represents the volume density of flexible load; A represents the cross-section area of flexible load; EI represents the elastic modulus of flexible load. Substituting Equation (2) into Equation (4), Equation (5) can be obtained: where ω i represents the modal frequency of the flexible load. Equation (5) can be sorted out to get Equation (6): Assume the expression of the modal function solution, as shown in Equation (7): According to Equations (6)-(9) can be obtained: where β i represents the characteristic root of the modal function. Its value is related to the flexible load length and boundary conditions. When the boundary conditions of the flexible load are consistent with the cantilever beam, the expression of the characteristic root value of the modal function is shown in Equation (10): According to Equation (8), the relationship between the modal frequency of the flexible load and the characteristic root of the modal function can be obtained, as shown in Equation (11): Thus, the characteristic equation and characteristic root of the modal function can be obtained according to Equation (9), as shown in Equation (12): According to Equation (12), the solution of the modal function can be obtained, as shown in Equation (13).
where D 1 , D 2 , D 3 , and D 4 represent the coefficients that can be obtained by the boundary conditions. According to literature [27], Equation (13) is sorted out to obtain Equation (14).
where the expression of ς i can be obtained in Equation (15).

Modeling of the Flexible Manipulator Servo System
If the flexible manipulator moves in the horizontal plane, it does not need to consider gravity's influence. In this case, the vector representation of any point on the flexible load is shown in Equation (16): So, the kinetic energy of the flexible load is shown in Equation (17): The elastic potential energy of the flexible manipulator is shown in Equation (18): According to the Lagrange Lagrange's principle, Equation (19) can be obtained: where q 1 , q 2 represent the rotation angle of flexible load and the modal coordinates; Q i represents the generalized force. Equation (19) can be obtained by Equation (20): Equation (20) can be finally sorted into the form of Equation (21): Equations (22) and (23) are defined: According to Equations (22) and (23), the servo system's dynamic equation for flexible manipulators can be summarized in the form shown in Equation (24): After removing the nonlinear terms, Equation (24) can be sorted out, as shown in Equation (25): Equation (25) represents the servo system's dynamic equation for flexible manipulators moving in the horizontal plane.

Modeling of Flexible Manipulator Servo System Considering the Gravity
If flexible manipulators move in the vertical plane, the influence of gravity needs to be considered. The flexible manipulator's kinetic energy is still shown in Equation (17), but its potential energy can be approximately expressed, as shown in Equation (26): According to Lagrange's principle, the dynamic Equation of the servo system for flexible manipulators considering the influence of gravity can be obtained, as shown in Equation (27): Similarly, after removing the nonlinear term, the simplified dynamic equation can be obtained, as shown in Equation (28): .
According to Equations (22) and (23), Equation (27) can be sorted into the form shown in Equation (29): According to Equation (29), Equation (30) can be obtained after removing the nonlinear term: Equation (29) represents the dynamic equation of the servo system for flexible manipulators considering gravity. Equation (31) represents the simplified form when nonlinear terms are not considered.
In order to study the influence of gravity on the input torque of flexible manipulators, a single link maneuver proposed by the literature [29] is also adopted. The relationship between the rotation angle and time is shown in Equation (31): ..
When gravity is taken into account, parameters such as the length and elastic modulus of flexible manipulators will affect the input torque. If only the first-order model is considered and the mass of the flexible manipulator is unchanged, the input torque can be obtained according to Equations (25) and (30), as shown in Figure 2.  According to Figure 2, the following conclusions can be drawn: (1) If the flexible manipulator moves on the vertical plane, gravity's influence on the input torque is crucial and cannot be ignored. If the factor of gravity is considered, the input torque of the flexible manipulator is less affected by the maximum rotation speed. According to the above conclusion, when the flexible manipulator moves in the vertical plane, the influence of gravity should not be ignored.

Transfer Function of the Servo System
When the flexible manipulator moves in the horizontal plane, gravity's influence does not need to be considered. Equation (25) can be used as the dynamic equation of the servo system for flexible manipulators. The dynamic equation shown in Equation (25) can be regarded as a single-input multiple-output system. The input torque of the flexible manipulator is taken as the input, and the rotation angle and the modal coordinates are taken as the output.
After the Laplace variation of Equation (25), Equation (32) can be obtained: Only the first-order model, the servo system's transfer function for flexible manipulators, can be obtained according to Equation (32), as shown in Equation (34): Equation (34) represents the transfer function from the input torque of flexible manipulators to the rotation angle. Equation (35) illustrates the transfer function from the input torque of flexible manipulators to the modal function.
If the effect of flexibility is not considered, the transfer function from the input torque of the manipulator to the rotation angle is shown in Equation (36): If we consider that flexible manipulators have a third-order mode, the transfer function from the input torque to the rotation angle is obtained, as shown in Equation (37): where the expression of M is shown in Equation (38): According to Equation (37), the open-loop Bode diagram from the input torque to the flexible manipulator's rotation angle can be obtained, as shown in Figure 3. According to Figure 3, the flexibility of the manipulator will cause oscillation frequency. When the oscillation frequency is close to the speed, it will cause more significant vibration. Therefore, it is necessary to consider the flexibility of the manipulator.  According to Figure 4, each mode's oscillation frequency decreases with the increase of the manipulator's length and increases with the growth of the elastic modulus. The manipulator will change the oscillation frequency in the telescopic movement due to the flexible factor, so it has strong robustness. This robustness will cause the motion error of the flexible manipulator and affect the accuracy. Therefore, the flexible factor should be considered for high precision manipulator: By comparing Equations (36) and (37), the flexibility of manipulators is mainly reflected in the flexible coupling part, and its expression is shown in Equation (39). Compared with rigid manipulators, flexible manipulators have a flexible coupling part in the transfer function:

Design of the PD Controller Based on Pole Placement Strategy
In the servo system for flexible manipulators, it is considered in the literature [30] that the control of the speed loop and current loop of the servo system has been completed. Therefore, the influence of the current loop and speed loop is ignored in this paper and considers position control. The position loop of the servo system adopts a PD control strategy, as shown in Figure 5. where K p , K d represent the proportion coefficient and the differential coefficient.
According to Figure 5, the closed-loop transfer function of position control of the flexible manipulator can be obtained, as shown in Equation (40): When only the first-mode is considered, the position control's closed-loop transfer function is shown in Equation (41): According to the literature [21], the pole placement strategy is adopted to configure the closed-loop transfer function's poles. The denominator polynomial of Equation (41) is written in the form shown in Equation (42): where ξ a1 , ξ b1 denote the damping coefficients of the poles; ω a1 , ω b1 denote the natural frequency coefficients of the poles. According to Equation (42), Equation (43) can be obtained: Therefore, Equations (44)-(47) can be obtained: The determination of PD controller parameters is the core of using the PD control method to control the servo system for flexible manipulators. According to Equations (44)-(47), pole placement strategy is used to calculate controller parameters. The PD controller parameters are shown in Equations (48) and (49): In this paper, the pole placement strategy with the same damping coefficient is used to determine the controller parameters. That means that the closed-loop transfer function of the system has the same damping coefficient. Substitute ξ a1 = ξ b1 = ξ 1 into the Equations (44)-(47), and Equations (50) and (51) can be obtained after sorting out: According to Equations (50) and (51), Equations (52) and (53) can be obtained: According to Equations (52) and (53), Equations (54) and (55) can be obtained: According to Equations (54) and (55), the control parameters can be obtained as shown in Equations (56) and (57): According to the control parameters, the servo system's closed-loop transfer function for flexible manipulators can be obtained, as shown in Equation (58): The parameters of the PD controller can be adjusted indirectly by adjusting the damping coefficient of the poles. The controller parameters can be calculated according to Equations (56) and (57) to achieve a stable and fast output of the system.

Control Method of Flexible Manipulators Based on RBF Neural Network
When flexible manipulators move in a vertical plane, the gravity factor cannot be ignored. Due to the nonlinear factors caused by gravity and flexibility, the difficulty of flexible manipulators control is increased. Therefore, the nonlinear control method of flexible manipulators is adopted in this paper. Nonlinear compensation is introduced into the inner control loop to offset the influence of gravity nonlinear term. The nonlinear part caused by flexibility is regarded as the disturbance term, which increases uncertain items in the system. During the compound motion of flexible manipulators, the parameters such as inertia have significant time-varying characteristics. Therefore, the servo drive system for flexible manipulators has strong robustness, bringing difficulty to precise control. In order to solve the uncertain items caused by flexibility and time-varying characteristics of the parameters, the RBF neural network is used to approximate the uncertain items. The accurate input torque of flexible manipulators can be obtained by introducing an internal control loop with nonlinear terms and identify the system's uncertain items with the RBF neural network. Then, the angle error and the vibration amplitude of flexible manipulators are reduced.
According to Equation (30), the servo system's dynamic equation for flexible manipulators considering gravity can be obtained. Equation (30) is simplified into a general form, as shown in Equation (59): where C(θ), G(θ), represent the damping coefficient and the gravity item respectively. If the damping effect is ignored, then C(θ) = 0. The symbol τ and d represent the input torque of the flexible load and the disturbance term respectively.

Internal Control Loop with Nonlinear Compensation Term
According to Equation (30) and Figure 2, the gravity factor cannot be ignored during the movement of flexible manipulators in a vertical plane. When calculating the input torque of flexible manipulators, whether the gravity factor is considered or not has a significant difference in the calculation results. In this paper, the inner control loop is added based on the PD control strategy that controls the servo system's position loop for flexible manipulators. The parameters of the PD controller are determined by pole placement strategy. A nonlinear gravity compensation term is added to the internal control loop. The input torque generated by the system can compensate for the influence of gravity in real-time. The schematic diagram of the servo system's internal control loop for flexible manipulators with nonlinear compensation is shown in Figure 6.
It is assumed that the uncertain items are not considered in the servo system for flexible manipulators. According to Figure 6, the control law of the PD control strategy with the inner control loop can be obtained, as shown in Equation (60):

Control Law Considering the Uncertain Items
Due to the time-varying characteristics of parameters and flexible nonlinear terms in the servo system for flexible manipulators, proper system parameters cannot be obtained. Therefore, only nominal models can be established. For the nominal system, the control law is shown in Equation (62): whereÎ a ,Ĉ(θ),Ĝ(θ) represent the nominal coefficients of the inertia, the damping coefficient, and the gravity item respectively; θ * represents the desired angular trajectory. By substituting the control law of the nominal system into Equation (59), Equation (63) can be obtained: After simplification of Equation (63), Equation (64) can be obtained: Equation (64) shows the closed-loop system equation with deterministic considerations, as shown in Equation (65): where ∆I a , ∆C(θ), and ∆G(θ) represent the difference between nominal coefficients and real coefficients of the inertia, the damping coefficient, and the gravity item, respectively. According to Equation (65), the uncertain items can be obtained, and its expression is shown in Equation (66): Based on Equations (65) and (66), the error state equation of the system can be obtained under the condition of considering system uncertain items, as shown in Equation (67): ..
Equation (67) can be written as shown in Equation (68): Assuming that the uncertain items are known, the control law of servo system for flexible manipulators is shown in Equation (69): By substituting Equation (68) into Equation (59), the stable closed-loop system equation can be obtained, as shown in Equation (61).
Because flexible manipulators' parameters are time-varying and the flexible item is nonlinear, the uncertain items cannot be obtained. Therefore, it is necessary to identify the uncertain items to ensure the precise input torque of flexible manipulators.

RBF Neural Network to Identify the Uncertain Items of the Model
In this paper, the RBF neural network is used to approximate the uncertain items of the servo system for flexible manipulators. The result after approximation by the neural network replaces the real uncertain items of the system.
In the RBF neural network, the most common activation function is Gaussian [31]. The Gaussian activation function has the following advantages: (1) The representation is simple. It does not add too much complexity to multi-variable input.
(3) The smoothness is good, and the derivative of any order exists.
In this paper, the Gaussian function is used as the activation function of the RBF neural network, and the output of the ith node of the neural network can be expressed by Equation (70): where u i (x) represents the output of the ith hidden node; σ i represents the normalizing constant of the ith hidden node; n represents the number of nodes of the hidden layer; c i represents the center vector of the ith hidden node of the Gaussian function; x represents the input samples.
The hidden layer of the RBF neural network to the output layer realizes linear mapping. The result of the output layer is shown in Equation (71): where y k represents the output of the kth node of the output layer; w k represents the weight coefficient vector; w ki represents weight coefficient. In this paper, there is only one output node, which represents the system uncertain items. According to Equations (70) and (71), it can be known that the system uncertain items are a function of the input samples and the weight coefficient vector. Use thef (x, w k ) to represent the uncertain items of the system identified by the neural network. The real uncertain items are a function of the state vector of the error. Use the f (x) to represent the uncertain items of the real system. This article makes the following assumptions: (1) The RBF neural network output continuously.
(2) The output of the RBF neural network approximates a continuous function f (x), and there is a very small positive number ε 0 that satisfies the condition of Equation (72): According to the above assumptions, the error state-space equation shown in Equation (68) can be written into the form shown in Equation (73): where w * k represents the neural network weight coefficient vector for the optimal identification of system uncertain items. Its expression is shown in Equation (74): This paper uses the η to represent the neural network error, and its expression is shown in Equation (75): According to Equation (75), the error state space-equation is written in the form shown in Equation (76): After training, the neural network's weight coefficient vector becomes w * k , then the output value of the neural network at this time is shown in Equation (77): According to Equation (66), the error state-space equation is written in the form shown in Equation (78): After applying the RBF neural network to identify the uncertain items, the servo drive system's control law for flexible manipulators is shown in Equation (79): whereŵ k represents the estimated value of w * k . After using the RBF neural network to identify the uncertain items, the control loop's schematic diagram is shown in Figure 7.

Proof of Stability
The RBF neural network is used to approach the uncertain items of the servo system for flexible manipulators. The control law is obtained, as shown in Equation (79). By substituting Equation (79) into Equation (59), Equation (80) can be obtained: ..
According to Equation (83), the error state-space equation can be obtained, as shown in Equation (84): where w k represents the difference between the optimal weight coefficients vector of neural network after training and the estimated weight coefficients vector; u(x) represents the output vector of the hidden node of the neural network. The Lyapunov function is defined as Equation (85): where γ is a positive number; P is a positive definite matrix and satisfies the Lyapunov function, as shown in Equation (86): where Q is a positive definite matrix. According to Equation (85), Equation (87) can be obtained: According to Equation (87), Equation (88) can be obtained: In Equation (88), the expression of M is shown in Equation (89): So, Equation (90) can be obtained: The adaptive law of the RBF neural network is shown in Equation (91): By substituting Equation (91) into Equation (90), the derivative of the Lyapunov function can be obtained as shown in Equation (92): By choosing the matrix Q and the small error η, the Lyapunov function's derivative can be less than 0. Under this condition, the servo system for flexible manipulators is stable.

Numerical Simulation Analysis and Experiment
In this paper, the control method combining the pole placement strategy and the RBF neural network is used to improve the motion accuracy of flexible manipulators. The PD controller parameters are optimized according to the pole placement strategy. Besides, the RBF neural network is used to identify the uncertain items of flexible manipulators. So, the input torque of flexible manipulators is more accurate. Therefore, the precision of the rotation angle of flexible manipulators is improved. In order to verify the effectiveness of the proposed control strategy, numerical simulation and physical control experiments of flexible manipulators are carried out in this paper. Firstly, the influence of the physical parameters of flexible manipulators on the input torque and resonance frequency are analyzed by numerical simulation experiments. Next, the influence of the damping coefficient of poles on the rotation angle is analyzed. Finally, the numerical simulation experiment and experiment of the servo control for flexible manipulator are carried out.

The influence of Physical Parameters of the Manipulator
The length and elastic modulus of flexible manipulators will affect the servo system's input torque and resonance frequency. Among them, the value of the input torque is closely related to the length of the manipulator. Besides, the input torque is closely associated with the gravity term. In this paper, according to Equation (25), Equation (30), and Equation (31), the relationship between input torque and time under different manipulators' lengths is respectively calculated in two conditions, as shown in Figure 8. One of the cases ignores the gravitational factor, and the other case considers the gravitational factor. According to Figure 8, the servo system's maximum input torque increases with the increase of the length of manipulators. Among them, the gravity factor has a more significant influence on the maximum input torque. Relative to the manipulator's length, the gravity factor is the dominant factor for the maximum input torque. Therefore, the internal control loop with nonlinear term compensation can effectively reduce the maximum input torque. Moreover, this control strategy can effectively reduce the fluctuation degree of the servo system's output torque.
The length and elastic modulus reflect the flexibility of manipulators. They also have an impact on the resonance frequency of the servo system. According to Equation (37), the influence of the flexible manipulator's physical parameters on the resonant frequency can be obtained, as shown in Figure 9. According to Figure 9, the flexible manipulator's resonance frequency decreases as the length increases and increases as the elastic modulus increases. The conclusion drawn in Figure 9 is the same as that in Figure 4. It can be seen that the stronger the flexibility of manipulators, the smaller the resonance frequency. Vibration is more likely to occur in the low-frequency phase.

The Influence of Coefficients of the Pole on the Rotation Angle
Before using the RBF neural network to compensate for the uncertain items, the PD controller parameters need to be determined. The value of the controller parameter directly affects the control accuracy of the servo system output. In this paper, the pole placement strategy is used to optimize the parameters of the PD controller. It provides a theoretical basis for the parameter selection of the PD controller. According to the pole placement strategy, the selection of controller parameters can be attributed to selecting the damping coefficient of poles. Different damping coefficients of poles will affect the output of the servo system. The unit step signal is used as the system's input signal to obtain the variation rules of rotation angle with time, as shown in Figure 10. Among them, the numerical simulation selected 4 different parameters, as shown in Table 1.  According to Figure 10, it can be seen that the variation rule of the rotation angle is not only related to the damping coefficient of poles but also closely related to the length and elastic modulus of flexible manipulators.
In this paper, three indexes of maximum overshoot, peak time, and adjustment time are selected to evaluate the output characteristics [21]. Among them, the maximum overshoot reflects the system's dynamic accuracy; the peak time and adjustment time reflect the rapidity of the system. If the maximum overshoot is smaller, the system is more stable. On the contrary, it shows that system stability is poor. If the peak time and adjustment time are shorter, the system rapidity is better. On the contrary, it indicates that the system reflects slowly. The influence law of the length, elastic modulus, and damping coefficient of poles on the evaluation index is shown in Figure 11.
According to Figure 11, as the flexible manipulator's length increases, the maximum overshoot of rotation angle does not change, but the adjustment time and peak time increase significantly. It can be seen that with the rise in the length of flexible manipulators, the flexibility of the manipulator rises, and it takes a long time to stabilize the system. With the increase of flexible manipulators' elastic modulators, the maximum rotation angle does not change, but the adjustment time and peak time decrease. It can be seen that with the increase of the elastic modulus, the flexibility of manipulators is weakened, and it takes a shorter time to stabilize the system. However, the damping coefficient of poles has a greater influence on the overshoot. With the increase of the damping coefficient of poles, the overshoot decreases gradually. The damping coefficient of poles has a weak effect on the adjustment time and peak time. According to Figure 11, it is appropriate to select poles' damping coefficient between 0.5 and 0.7. In this interval, the system can obtain a faster response and receive a smaller overshoot.

Control Method Based on the Combination of RBF Neural Network and Pole Placement Strategy
In order to verify the effect of the combination control method, the parameters in Table 1 are selected to carry out a numerical simulation experiment. The simulation experiment takes the sine function as the input signal and obtains a change in the rotation angle over time, as shown in Figure 12. The influence of the manipulator's length and elastic modulus on the rotation angle error under different control strategies is shown in Figure 13.
According to Figure 12, if the pole placement strategy is used alone, stable tracking can be achieved when manipulators' length is small. With the increase of the length of manipulators, pole placement's control strategy cannot satisfy the precise tracking. However, after using the neural network to compensate for the uncertain items, the combination control method can ensure the system's stability.
According to Figure 13, it can be seen that the control method using pole placement strategy alone will have a large error. Additionally, the error increases as the length of manipulators increases. The error decreases with the increase of elastic modulus. It is shown that the control method using pole placement strategy can only satisfy the condition of low flexibility. The control effect is not suitable for the situation of large flexibility. However, the combined control method is not affected by the flexibility of manipulators. The combined control method can obtain a better control effect in any cases.

Experiment
This paper builds an experimental platform to verify flexible manipulators' control strategy, as shown in Figure 14. The experimental platform is composed of a servo motor, a transmission tendon-sheath, and a flexible manipulator. The tendon-sheath transmits the input torque generated by the servo motor to the flexible manipulator. The input torque of the flexible manipulator can be accurately calculated through the tension sensor on the tendon-sheath. The control platform uses CRIO-9053 as the lower computer that inputs the collected feedback signals into LabVIEW's control program. Simultaneously, the control program's output signal is transmitted to the NI-9246 module to control the output torque. Magnetic encoders collect the feedback signal. The magnetic encoder is installed at the back of the flexible manipulator to collect the angle signal. The control principle of the experimental platform is shown in Figure 15.  The flexible beam is used to simulate the flexible manipulator in the control experiment platform. Flexible beams with lengths of 0.5, 1, and 1.5 m are used to simulate flexible manipulators of different lengths. Two different control methods are used to experiment with rotation angle control. One of the control methods is to use the pole placement strategy alone, and the other is to use the combined control method. The relationship between the rotation angle with a length of 0.5 m and its error are shown in Figure 16. Six groups of experimental data can be obtained through the experiment. The statistical results of the rotation angle errors of the experimental data groups are shown in Figure 17. What is more, the error indicators of six groups of experimental data can be obtained as shown in Table 2.   It can be seen from Figure 16 that the combined control method has a smaller error. According to Figure 17 and Table 2, it can be seen that the combined control method can make the flexible beams of different lengths get smaller tracking errors. However, the pole placement strategy alone cannot guarantee the flexible beam's tracking accuracy in the case of a long length. Therefore, the effectiveness of the control strategy proposed in this paper is verified through experiment.

Conclusions
In this paper, a dynamic model of flexible manipulators with gravity is established. The innovation of this paper is to consider the factor of gravity and the combined control strategy is used to improve the motion accuracy. The control method combining the pole placement strategy and the RBF neural network is applied to reduce the fluctuation of the rotation angle of flexible manipulators. Then, the motion precision of flexible manipulators is improved. Among them, the RBF neural network is used to distinguish the uncertain items of the system. The uncertain items include both the flexible factors and the time-varying characteristics of the dynamic parameters. Simulation analysis and experiments show that the proposed control method can effectively suppress the rotation angle's vibration and improve the motion accuracy of the end-effector. The specific conclusions are as follows: (1) Gravity factors will seriously affect the input torque of flexible manipulators. If flexible manipulators move in the vertical plane, the gravity factor should be taken into account. (2) Its length and elastic modulus determine the flexibility of manipulators. The PD controller can control the manipulator with less flexible. The pole placement strategy is used to optimize the PD controller's parameters to obtain a stable output of the rotation angle. However, for the manipulator with strong flexibility, the pole placement strategy cannot receive a stable output of the rotation Angle. Using the RBF neural network control strategy to identify the uncertain items containing the flexible factors can effectively reduce flexibility. (3) The simulation and experimental results show that the control method combined with the RBF neural network and the pole placement strategy can effectively reduce the error of the flexible manipulator's rotation angle. According to the Table 2, compared with pole placement strategy alone, the mean error is reduced by nearly 60%. Therefore, the combined control method can effectively reduce the angle error. Therefore, the control method proposed in this paper can effectively improve the control accuracy of the flexible manipulator.
This paper uses the control method to improve the motion accuracy of flexible manipulators. It is hoped that the vibration isolation device can be added between flexible manipulators and end-effectors in the future. In this way, the motion accuracy of flexible manipulators can be improved. Data Availability Statement: Some data, models, or code generated or used during the study are available from the corresponding author by request.