#### Abstract

Considering the complex dynamic modeling of multi-DOF planar flexible manipulators, a general-purpose method for the rigid-flexible coupling dynamic modeling of *N*-DOF flexible manipulators is proposed in this paper, and symbolic calculation software is developed. The modeling method is based on the Lagrange equation and assumed mode method (AMM). First, the *N*-DOF flexible manipulator is divided into two parts, which are assumed to be rigid and flexible. On this basis, the rigid part and the flexible part are coupled, and the calculation process of the model is further simplified. Then, the simplest general symbolic expression of the dynamic model of the *N*-DOF flexible manipulator is obtained with the induction method. According to the modeling method, “symbolic expression computation software for dynamic equations of *N*-DOF flexible manipulators” is developed using the symbolic calculation software Mathematica. Finally, the dynamic modeling method and the symbolic calculation software are verified by a trajectory tracking experiment with a PD control applied to a 2-DOF flexible manipulator. Compared with the traditional modeling method, the calculation time can be reduced by more than using the modeling method proposed in this paper, which significantly reduces the complexity of the modeling process.

#### 1. Introduction

Robotic manipulators are widely used in environments that are hazardous to humans and in situations where there are substantial amounts of repetitive work to reduce the burden on humans and to assist in accomplishing special tasks. To ensure sufficient positioning accuracy at the end of a manipulator, the links of the manipulator can generally be stiffened by a bulky structure to prevent the links from vibrating and to avoid elastic deflection during operation [1–4]. However, heavy links limit the speed of the manipulator’s movement and consume excessive amounts of energy. In some cases, thickening the structure of the links to ensure the positioning accuracy of the end of the manipulator may not be practical. For example, the weight of a manipulator operating on the space station must be as light as possible because of launch costs [5]. In addition, to ensure that the manipulator can perform in a greater workspace, the boom of a concrete pump truck structure should be more slender [6]. Robotic manipulators are also used in nuclear maintenance [7, 8], for spraying large aircraft [9], etc. Under these circumstances, these robotic manipulators have to be designed as flexible-link manipulators (FLMs) with a slender structure.

FLMs have been found in extensive applications because of their advantages in terms of being lightweight and having high load ratios (the ratio of load mass to manipulator mass), high movement speeds, and low energy consumption [10–12]. Moreover, in human-robot interactions, a lightweight FLM is safer to humans and has smaller workload requirements on joint drivers, sensors [13, 14], etc. However, such an FLM experiences vibration and elastic deformation during movement, which will greatly affect the stability of the entire system and the positioning accuracy of the end [15, 16]. Therefore, it is important to study the effects of the dynamic characteristics of FLMs, in terms of the vibration suppression and elastic deformation compensation of the flexible links. For example, the dynamic model of a single-joint flexible manipulator has been developed to suppress the vibration of the flexible manipulator [3, 13].

Because multijoint FLMs are highly nonlinear, strongly coupled, multi-degree-of-freedom (DOF) systems, the dynamic modeling is complex. In addition, considering the payload mass at the end of the manipulator and the lumped mass (for example, the mass of the actuators) on each link, multi-DOF rigid-flexible coupling dynamic modeling is more complicated. These challenging topics have attracted a large number of excellent researchers, who have conducted a considerable amount of research on the dynamic modeling of FLMs. In Giorgio and Del Vescovo’s work [17], the lumped parameter method is used to model the dynamics of flexible links, and although this modeling method is simple, the number of elements must be as large as possible to obtain an adequate model accuracy. The modeling and control of a special class of single-link, an arbitrary number of lumped masses, flexible arms is presented in [18], in which the arm consists of massless flexible structures. Subudhi and Morris adopted the Lagrange equation and the assumed mode method (AMM) to establish the dynamic model of a manipulator with flexible links and joints [19]. A dynamic model of 3-PRR parallel manipulator with flexible links based on AMM was proposed in [20]. Wei et al. provided a global mode method for dynamic modeling of a flexible arm flexible joint manipulator with a tip mass [21]. Some other modeling methods, such as the finite element method (FEM), are also commonly used. For example, the dynamic modeling of the flexible links using FEM is presented in [22]. The FEM is used to derive the dynamic equations for finding the maximum allowable dynamic load of flexible-link mobile manipulators [23]. In addition, the Newton–Euler method is also used to dynamically model the flexible arms [24]. Nonlinear dynamics of motion are developed for flexible manipulator links consisting of rotary joints that connect pairs of flexible links [25]. Bascetta et al. developed a closed-form dynamic model of flexible manipulators based on the Newton–Euler equations [26]. However, the dynamic modeling in the above studies only considers flexible arms with a certain number of DOFs or flexible links without payload mass, and the expressions of the dynamic model are still very computationally expensive in practice.

Considering the complexity of the dynamic modeling of multi-DOF FLMs with a lumped mass on the links, a general-purpose dynamic modeling method for *N*-DOF FLMs based on the Lagrange equation is presented in this paper, and the symbolic expression computation software for the dynamic equations of *N*-DOF flexible manipulators is developed. In a more general situation, where the lumped mass is located not only on the end of the links but also at other positions of the links (for example, the flexible arms composed of functionally graded materials (FGM) in the fields of spacecraft, nuclear industries, robotics, etc., the flexibility of the arms can be increased or changed, or changing the vibration frequency of the flexible arm structure by adding concentrated mass to the intermediate position of flexible arms, see [27]) is also considered in this paper. The modeling method is as follows. First, the flexible arms are divided into imaginary rigid and flexible parts and modeled separately. Then, the rigid and flexible parts are coupled, and the dynamic model is further simplified. Finally, the simplest general expression of the dynamic model of the flexible arm with *N*-DOFs is obtained by the induction method, and “symbolic expression computation software for dynamic equations of *N*-DOF flexible manipulators” is developed using the symbolic calculation software Mathematica. The test results indicate that compared with the traditional modeling method, the proposed modeling method can reduce the computation time by more than and the software can quickly calculate and generate the dynamic model of the flexible manipulator, thereby greatly reducing the difficulty of the dynamic modeling of the flexible manipulator. Finally, the position tracking control experiment of a 2-DOF flexible manipulator is taken as an example and the dynamic model is generated automatically by the program. Meanwhile, a physical model is also constructed. In an experiment on the end effector of the manipulator tracking a trajectory of a circle, the numerical simulation results of the dynamic model are compared with the experimental results. The dynamic model and the dynamic modeling method are both verified to be correct.

The remainder of this paper is organized as follows. In Section 2, the motion of *N*-DOF FLM systems is described. Then, in Section 3, the dynamic equations of *N*-DOF FLMs are derived and established. “The symbolic expression computation software for dynamic equations of *N*-DOF flexible manipulators” is developed in Section 4, in which the dynamic equations of two-DOF FLMs are automatically generated by the program, and the experiments are conducted. Finally, conclusions are drawn in Section 5.

#### 2. Motion of Flexible Links

To study the general situation of dynamic modeling for serial multijoint FLMs and to simplify the derivation process, fixed-base and small-deformation flexible manipulators moving in the vertical plane are the focus of this paper. Before we derive and establish the dynamic model of the multi-DOF FLM, the appropriate float frame for each flexible link should be assigned, as shown in Figure 1.

The fixed base of the manipulator is typically treated as the global inertial reference frame and denoted , and the *i*-th flexible link is associated with the float frame . The *x*-axis of the float frame passes through the start and end points of the *i*-th flexible link, and its origin coincides with the end point of the *i*-th flexible link. Assuming that the length of the *i*-th flexible link without deformation is , the angle between the *i*-th and *i* − 1-th links along the *x*-axis is . Thus, according to reference [28], the homogeneous transformation matrix from frame to the coordinate system is derived as

In equation (1), the superscript of the homogeneous transformation matrix represents the reference frame and the subscript *i* represents the current frame . Thus, the position vector of an arbitrary point on link *i* can be expressed in the global inertial frame aswhere , is the position vector of point on link *i* in float frame , and .

According to the aforementioned assumption, the flexible link here is considered to be an Euler beam with a constant cross section. For a small deformation of the link, because the deflection of the flexible links in the length direction is very small, we can ignore it and only consider the transverse deformation. Because the exact mathematical description and calculation of the flexible link deformation are complex given that the flexible link deformation is infinite dimensional, a finite-dimensional deformation is generally introduced to approximate the link deformation. The assumed mode method (AMM) is adopted in this paper. In the AMM, the elastic displacements are generally described by a truncated modal series function and time-dependent mode amplitudes (for example, see [29]) aswhere is a *k*-th-order mode shape function, are generalized coordinates with respect to the of link *i*, and *s* is the number of modes used to describe the deflection of link *i*.

Generally, this method has sufficient approximation accuracy for small deformations of flexible links when we truncate at . Thus, in the following, we only consider the first second-order mode shape function and its generalized coordinates. For the Euler beam with clamped boundary conditions at both ends of the link, we choose the mode shape function . Thus, equation (3) can be rewritten as

The velocity vector of in frame can be expressed as

In general, the elements in matrix are a function of and the elements in vector are a function of .

Note that the partial derivative of with respect to can be calculated aswhere, for a revolute joint, the matrix is defined as

The advantage of using is that we can transform the differential operation of into a constant matrix multiplying the matrix, which can still be used; then, we can obtain

To simplify the derivation of dynamic equations of multi-DOF FLMs, we can divide the flexible link system into two parts: an imaginary rigid part and an imaginary flexible part. The expressions include the variables , representing the rigid part, and the variables or , representing the flexible part. Now, we imagine that the *x*-axis of frame from the start point to the end point of flexible link *i* is a virtual rigid link *i*, which can be regarded as a rigid link of the flexible link *i* without deformation; thus, equation (8) can be interpreted as the effect of the motion of joint *j* on all the points on virtual rigid link *i*. To simplify the notations, we define

Thus, equation (5) can be rewritten as

For multi-DOF manipulators, the dynamic equations are a strongly coupled system; thus, considering the interaction effects between virtual rigid links, we can obtain

Equation (11) can be interpreted as the interaction effects of the motions of joint *j* and joint *k* on all points on virtual link *i*.

#### 3. Derivation of Dynamic Equations for Flexible Links with Lumped Mass

##### 3.1. Calculation of Kinetic Energy

After obtaining the velocity of points on link *i*, we can calculate the kinetic energy of link *i*. Assuming that the mass of flexible link *i* is , let be the kinetic energy of link *i* and let be the kinetic energy of a particle with differential mass on link *i*; then, we can obtain thatwhere Tr is a trace operator and . Now, let us compute the kinetic energy of lumped mass on flexible link *i*. We assume that the lumped mass can locate at any position on flexible link *i*. Similarly, we use the same method as equations (10) and (12) and obtainwhere is the lumped mass and and are the position vector and velocity vector of the lumped mass in the global inertial reference frame , respectively. and are the position vector and velocity vector of lumped mass in float frame , respectively. Thus, we let , as shown in Figure 2.

Because we ignore the deformation of the flexible link in the longitudinal direction, the length of flexible link *i* is always . To simplify the derivation, we define , . When , lumped mass locates on the end of flexible link *i* if represents at the middle of flexible link *i*. Thus, we can obtain , and according to equation (4), we can obtain

To simplify the notation, we define , . The total kinetic energy of flexible link *i* and lumped mass can be computed aswhere , , and . Additionally, is the pseudoinertia matrix of flexible link *i*, which, for a homogeneous Euler beam, can be calculated as

As shown, the pseudoinertia matrix is a function of generalized coordinates . The variation in the generalized coordinates can reflect the deformation of the flexible link. Similarly, is calculated as

Note that , and we can easily obtain and :

and are computed in the same manner:

##### 3.2. Calculation of Potential Energy

Let the gravitational potential energy of flexible link *i* with lumped mass be ; thus,where is a gravity vector expressed in the inertial reference frame and is the gravitational constant. Because of the small deformation of the flexible link, the variation in the gravitational potential energy of the flexible link caused by the flexible link’s deformation can be neglected and is small compared with the total gravitational potential energy of the flexible link. Let and be the mass center vector of flexible link *i* without deformation in inertial frame and float frame , respectively.

Meanwhile, the elastic potential energy of flexible link *i* caused by deformation can be expressed aswhere and are Young’s modulus and the second moment of the area of flexible link *i*, respectively. Now, from equations (12) to (21), we define the Lagrangian function of FLMs with *N*-DOFs as

We now apply the Lagrange–Euler formulation to the Lagrangian function of the *N*-DOF FLM equation (22) to obtain the necessary generalized torque to drive the manipulators:where and are a vector of generalized coordinates and its derivative with respect to time, respectively. We define generalized coordinates , , and , and is a generalized torque vector.

Before we apply the Lagrange–Euler formulation to the Lagrangian function of equation (22), we should note that the transformation matrix and the matrix are only functions of generalized coordinates ; thus, the derivation of or with respect to will be zero. The matrices , , and are only functions of the generalized coordinates and generalized velocity ; thus, the derivatives of , , and with respect to will be zero.

First, we only consider the generalized coordinate ; thus, take as and substitute equation (22) into equation (23) to obtain the rigid generalized torque vector :where . Then, we derive the flexible generalized torque vector by substituting equation (22) into equation (23) when considering as

For , where is a partial derivative of with respect to the generalized coordinate and .

Equations (24) and (25) are the dynamic equations of FLMs with *N*-DOFs. Because the entire process of deriving the dynamic equations from equations (23) to (25) is very complex, we divide the equations into two parts: a rigid part with respect to the generalized coordinate and a flexible part with respect to the generalized coordinate . Therefore, the complexity is greatly reduced. Note that in equations (24) and (25), we should calculate the derivative of the matrices , , and with respect to time *t* becausewhere and because the is canceled during the calculation, and thus, we define

It is easy to calculate and according to the expressions of and for :

From equations (24) and (25), we find that there are not only rigid motion terms with but also flexible motion terms with , in addition to the coupling terms between rigid motion and rigid motion and flexible link deformation . Notably, there are no flexible deformation coupling terms, which should have been . This makes solving the dynamic equations easier, and this is an important reason for establishing such float frames to describe the FLMs. In essence, the two equations also reflect the vibration of the FLMs during movement.

##### 3.3. Relationship between the Generalized Forces and Actual Driving Forces

We continue to discuss the relationship between the generalized forces corresponding to the generalized coordinates and the actual forces applied to the manipulator’s joints. As shown in Figure 3, is the actual force applied to joint *i*, where the first subscript *i* denotes the link that the force acts on and the second subscript denotes the link that the force reacts on. According to the principle of virtual work, the virtual work acting on link *i* iswhere and ; thus, the total virtual work on all joints iswhere and , and the generalized forces corresponding to the generalized coordinates are

The generalized forces corresponding to the generalized coordinates arewhere .

##### 3.4. Matrix Form of Dynamic Equations for FLM

To make the dynamic equations easy to understand, we rewrite the equations from equations (24) to (31) in matrix form aswhere , , , and are mass matrices and , , and are Coriolis and centrifugal matrices. and are gravity vectors. is a diagonal stiffness matrix. The detailed expressions are as follows:where andwhere . Because we only take the first two mode shape functions and generalized coordinates, we can regard every two-column vector as one group in the matrix , and “*λ*” means the *λ*-th column vector in the *j*-th group. Hence, represents the *i*-th row, *j*-th group, and *λ*-th element in the matrix :where . From the expression, we find that the matrix is the transpose of :where , and is a diagonal matrix.where , withwhere . The computation of the matrix is the same as that of :where . The calculation of the matrix is the same as that of the matrix :where .where . is a diagonal matrix, and the element is also a diagonal matrix.where . Here, we also regard every two-row vector as one group, similar to equation (35). Thus, represents the *i*-th group, *λ*-th element of the matrix .

We now rewrite equation (33) in a more compact and general matrix form aswhere is the mass matrix, the matrix consists of the centrifugal force and Coriolis force, the column vector denotes the gravitational and elastic force, and is a generalized force.

##### 3.5. Symbolic Expression Calculation Procedure

The complete symbolic expression of the rigid-flexible coupling dynamic equation of the *N*-DOF flexible manipulator with the lumped mass on an arbitrary location of links has been obtained. According to the whole derivation process of this paper and from equations (33) to (44), software that can automatically generate the rigid-flexible coupling dynamic equation of an *N*-DOF flexible manipulator is developed by using Wolfram Mathematica symbolic calculation software. The operating interface of the software is shown in Figure 4.

In the program interface, only the number of *N*-DOFs, the arm length, the section inertia moment, the elastic modulus, the mass of the flexible arm, the symbolic parameters of the end mass, and the location of the lumped mass are input. According to the program, the symbolic expression of the dynamic equation of the flexible manipulator with N-DOFs can be obtained conveniently and quickly. Here, we rearrange the symbolic expression computation procedures as follows:(i)*Step 0*. Initialization of the algorithm The number of DOFs and the parameters of the FLM, such as , are predefined Computation of are predefined and from (1); from (2) from (9); from (11)(ii)*Step 1*. Computation of the global mass matrix Computation of from (34); from (35) from (36); from (37)(iii)*Step 2*. Computation of the global matrices of the centrifugal and Coriolis forces Computation of from (38); from (39) from (40)(iv)*Step 3*. Computation of the global vector of the gravity, stiffness matrix, and generalized force Computation of from (41); from (43) from (42) from (31); from (32)(v)*Step 4.* Postprocessing of the result Simplicity of the dynamic model: some terms can be ignored for special situations Output of the symbolic expression for the model.

The flowchart of the symbolic expression computation procedures of the dynamic equation of the flexible manipulator with N-DOFs is shown in Figure 5.

##### 3.6. Comparison with Traditional Modeling Methods

The traditional dynamic modeling methods of flexible manipulators include the Lagrange method, Newton–Euler method, and Kane method, which have their advantages and disadvantages. In this paper, the traditional modeling method is compared with the proposed method, which is also based on the Lagrange equation. The traditional modeling method is as shown in [19].

The traditional modeling method is relatively simple, but it is quite complex in terms of the actual calculation process for solving the dynamic model expression and there is no general formula for flexible manipulators with different DOFs. Therefore, it often takes substantial time to calculate the dynamic model of various types of flexible manipulators in practice.

Here, the program developed in Section 3.5 is tested. Additionally, according to the modeling method in [19], the calculation programs for flexible manipulators with different DOFs developed using Mathematica are also tested. The time required for the two modeling methods to run the computational program (on a PC with a i3-7100 3.90 GHz CPU with 8 GB of RAM) was tested when *N* = 1, 2, …, 5. The comparison results are shown in Table 1.

As shown in Table 1, the modeling method proposed in this paper can reduce the modeling time by more than compared with the traditional modeling method. This result is mainly because in the traditional modeling process, the Lagrange energy function of the manipulator system requires a large number of integral operations; then, in the derivation of the generalized coordinates, there are many differential operations. Moreover, there are many repeated intermediate calculations in the differential calculations. The method proposed in this paper eliminates the repetitive intermediate calculation process theoretically and directly derives the simplest expression of the dynamic model equation.

Note that to explain the modeling method proposed in this paper, the simplest planar flexible manipulator is selected as the research object. However, from the entire modeling process shown above, the flexibility characteristics of rigid-flexible coupling manipulators are related only to the matrices , , and and their derivatives. Thus, for the dynamic modeling of complex spatial flexible manipulators, large-deformation flexible arms, or even irregularly shaped flexible links, we only need to modify the matrices , , and and their derivatives. The structure of the symbolic expressions remains unchanged; thus, this dynamic modeling approach also has good portability. It can be predicted that to satisfy the increasingly complex working conditions, the structure of the flexible arm will be more complex, and it will be a difficult problem to establish a sufficiently accurate dynamic model. Therefore, the modeling method proposed in this paper will have substantial advantages in terms of modeling convenience and computational efficiency.

#### 4. Experimental Study

To verify the correctness of the modeling method proposed in this paper, the dynamic model of a 2-DOF flexible manipulator is automatically generated using the software. Simultaneously, a flexible manipulator experimental platform is set up, as shown in Figure 6. Figure 7 is a schematic diagram of the assembly structure of the flexible arm device.

The expressions of the dynamic equations for the 2-DOF flexible manipulator are shown in the appendix. The symbols in the formula are the same as those in the derivation process.

In this experiment, the two arms of the flexible manipulator are controlled to follow the given trajectories. The torque output from the controller and the angles of the joints of the flexible manipulator measured by the encoders are recorded. Simultaneously, the angles are transmitted to the dynamics simulation model, and the estimations of the control torques are calculated for comparison with the actual situation. The whole experimental scheme and recorded data are shown in Figure 8. Figure 9 is a schematic diagram of the experimental signal acquisition and control.

According to the scheme shown in Figure 8, the Simulink control model is established and compiled into a “.so file,” the input and output (IO) interface between the control model and the real manipulator is configured using VeriStand, and finally, all the related control model files are downloaded to the real-time cRIO-9035 controller. The output torque calculated by the controller is converted into a PWM signal by dividing the motor torque constant and the gear ratio of the reducer. A 40 MHz FPGA chip is mounted onto the back of the control slot to acquire the high-speed pulse signal produced by the encoders through the digital IO interface. The motors, gears, and drivers are manufactured by Times Brilliant in China. The type of brushless DC motor 1 with the power of 150 W and brushless DC motor 2 with the power of 100 W is 57BL95S15-230TF9B and 57BL75S10-230TF9, respectively. The type of driver 1 and driver 2 is ZM-6610M. Reduction ratios of PX series planetary gears matched with motors are 64 and 32, respectively. The other specific parameters of the components in the experimental device are listed in Table 2.

##### 4.1. Case 1: Free Vibration of the End

In the experiment, the initial working positions of the flexible manipulator are selected as and . At the beginning, the external support makes the manipulator as free as possible from gravity and the second flexible manipulator keeps the horizontal position. Then, the flexible manipulator system is allowed to vibrate freely under the force of gravity, and the acceleration of the end of the flexible manipulator under free vibration is observed and recorded. In addition, the accelerometer type is WT901C485, the acquisition time and sampling time are 10 s and 20 ms, respectively, and the simulation interval of Simulink is 0.001 s. The measured acceleration graph is compared with the simulation graph, as shown in Figure 10(a). It is apparent that there are very good correlations between the measured values and simulated values and the amplitude of simulation is smaller, which can be due to model approximation and initial attitude error between the simulation model and the physical model. The spectrum of the end acceleration curve of the simulation model is obtained by fast Fourier transform, as shown in Figure 10(b). It can be seen that the first-order vibration frequency is *ω*_{1} = 12.246 rad/s and the second-order vibration frequency is *ω*_{2} = 40.82 rad/s. From the experiment, it is obvious that the simulation curves are basically consistent with the experimental curves, which verifies the correctness of the proposed modeling method to some extent. The similar results of vibration characteristics of flexible manipulators can be referred to [30], which can be used to compare with these experimental results.

**(a)**

**(b)**

##### 4.2. Case 2: End Effector Tracking a Trajectory of Circle

In this experiment, the initial working positions of the flexible manipulator are selected as and . A circle with a radius of 0.1 m in the XOY plane is selected as the end track of the flexible arm, as shown in Figure 11. The equations of the circle are , , whereas the center coordinate of the circle is , , . The end of the flexible arm is substituted into the inverse equation of rigid kinematics, and the desired joint angles are obtained by real-time operation. In addition, the simple PD control is adopted in the cRIO-9035 controller to make the joints track the desired angles, in which and were debugged through repeated simulation experiments and finally determined to be and , and the calculation interval is 0.001 s. Figures 12(a) and 12(b) are the angle curves of joint 1 and joint 2. It can be seen from the figures that when the flexible arm starts to move, the experimental curve and the simulation curve have some errors, and then the errors gradually decreases and the curves tracking effect is better. The input torque curves of the joint 1 and the joint 2 are as shown in Figures 13(a) and 13(b), and the vertical acceleration curves of the flexible arm 1 and the flexible arm 2 are as shown in Figures 14(a) and 14(b). It can be seen that the torque simulation curves and the experimental curves follow the same trend, and there are small tracking errors, due to the inevitable vibration of the flexible mechanical arm during the tracking process, the acceleration will cause frequent fluctuations, but the simulation curve and the experimental curve can also achieve a better tracking effect on the whole. The main cause of some vibration of the manipulator during movement is that, to simplify the control model, the distance between the output shafts of two motors is regarded as the arm length of the flexible arms, the mass of the motor reducers is treated as the lumped mass, and the unavoidable clearance of the motor reducers makes the simulation model and the actual model have certain errors, but these errors are within acceptable ranges.

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

#### 5. Conclusions

A general modeling method for the dynamic equations of an *N*-DOF flexible manipulator based on the Lagrange equation and the small deformation assumption of an Euler beam is proposed in this paper. According to this modeling method, symbolic calculation software for the dynamic equations of *N*-DOF flexible manipulators is developed, greatly reducing the computational complexity. Compared with the traditional modeling method, the modeling time can be reduced by more than . Finally, the correctness of the modeling method is verified by experiments.

The modeling method proposed in this paper takes the planar flexible manipulator as the research object, but according to the modeling process, it can be applied to the more complex rigid-flexible coupling dynamic modeling of flexible manipulators only through modifying the corresponding matrices of the flexible part. This is also future work that needs to be studied. Simultaneously, the modeling method proposed in this paper is also helpful for the dynamic modeling of *N*-DOF spatial flexible manipulators, large deformation flexible manipulators, etc.

#### Appendix

The dynamic equations for the 2-DOF flexible manipulator are shown as follows:

The matrix is omitted in the above expressions because , which yields

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

This work was supported by the Development and Demonstration of Intelligent Slag Feeder Robot for Continuous Casting Crystallizer, Hubei Science and Technology Support Plan (grant no. 2015BAA023).