Coupled-Nonlinear Elastic Structure: An Innovative Parameterization Scheme of the Motion Equations

In this paper, I applied the Euler-Lagrange equations in order to obtain the coupled-nonlinear motion equations for an elastic structure. The model is composed of six coupled and strongly nonlinear ordinary differential equations. The new contribution of this work arises from the fact that a convenient and innovative parameterization of the motion equations for the elastic system was developed with all mathematical nonlinearities taken into account, without the usage of any simplifying linearization procedure, as found in most of the works presented in the literature. The results can be used as a source for conducting experiments and can be useful for a better understanding and control of such nonlinear elastic systems.


Introduction
The use of the lightweight structural elements in space applications, underwater interventions as well as in robotic manipulators under the requirement for precise positioning, easier transportation, less power consumption, has increased the interest in having a precise model that closely represents such mechanical systems for all their conditions.The position and velocity vector obtained, after imposing the inextensibility conditions, are used in kinetic energy expression while the curvature is used in the potential energy.The Lagrangian dynamics [1] in conjunction with the assumed modes method [2] are utilized in order to obtain the non-linear [3] equations of motion that are treated with all non-linearities taken into account, without the usage of any simplifying linearization procedure, as found in most of the works found in the literature.The resulting non-linear model is com-posed of six coupled and strongly non-linear ordinary differential equations are discussed, simulated and some results of this simulation are presented.The way in which the motion equations were treated in this paper, allows the monitoring of each contributing factor for the system elasticity.The main goal of this work is to treat the motion equations according to a general approach, without the usage of any linearization procedure, and this way to permit assesses the system behavior through controlled simulations.

Problem Description
In this work the dynamic modeling is performed for a system that contains two elastic beams and two rotational joints.
A convenient parameterization of the terms of the motion equations, which makes it easier to compare the simulation results for the rigid and for the elastic structure [4] is also developed.The studied elastic system is assumed to have a planar movement.Even with this assumption, the complexity of the resulting dynamic equations for this system is large when compared to the equations for rigid structures [5].
The model is established basically by the superposition of the elastic movement with the movement of a hypothetical rigid body.A convenient mathematical set of equations is developed for this purpose.The elastic [6] movement of the beams is truncated in the second mode, that is, it is considered that the amplitudes of all higher order vibration modes [7] are much smaller than the amplitude of vibration of the first mode.In such a case, a non-linear model with six degrees of freedom is obtained.
I outline that, one of the tasks of this work is to treat the motion equations according to a general approach, without simplifying linearizations procedure.

The Model
The physical model for the elastic system is established in accordance to the schematic drawing shown in the Figure 1.
The following assumptions are made: the system has planar movement and the relative movement between the two links is resulting from the torques applied in each joint of the system.
At the terminal of the beam 1, a concentrated mass represents both the servo-motor and the joint masses.At the terminal of the beam 2, a discrete mass is used to represent the load to be handled between two points of the plane.
In order to describe the movement, three reference systems are defined:

[ ]
, , O X Y -inertial referential system with origin in joint 1; , , O x y -referential system with origin in O and the x 1 axis tangent to beam 1 at point O; , , O x y -referential system with origin in joint 2 and with the x 2 axis tangent to beam 2 at point O 2 .
Two angles are defined: ( ) is the angle between the x 1 and x axes; ( ) θ is the angle between the x 1 and x 2 axes.It is also defined a new system that is formed by the two segments OO flexible movement of the beams 1 and 2 with respect to this mobile system.

Kinematics Description
The position can be described by a convenient description of a set of coordinates.Any point P i may be specified if a new variable ( ) u x t is defined as being the coordinate of the flexible movement with respect to the ref- O x y .In order to simplify the notations, a matrix representation form of the reference frames can be introduced.Let them be: , is the unit vector of the frame OXY; { } , is the unit vector of the frame and where [ ] are the rotational transformation matrices and can be written: [ ] Ox y .The (generic) position vector i P (Figure 2) will be: where { } T i

U
means the transpose matrix.With the usage aformentioned equation, can be written the vector position of any point on the beam 1 of the following manner: where { }

C
means the transpose matrix.

Beam 2
In order to define the vector position of any point on the beam 2, it will be convenient to assume that the displacements of the elastic beams with respect to reference frames [ ] Ox y and [ ] O x y be small enough to consider the paths of points O 2 and O p as straight lines normal to the respective reference frames.Then, the position vector of any point P 2 on the beam 2 is established in accordance to the schematic drawing shown in Figure 3.It is also possible to note that, OO O O O P P P (7) If now, u 1E = elastic linear displacement at the end of the beam 1; L 1 , L 2 = lengths of the beam 1 and 2, then, Consequently, ) Similarly, one can obtain j  R and p  R , for a mass concentrated at the joint 2 and for the payload, respectively, and will be shown in the next section.

Kinetic Energy
By means of the velocity vectors previously mentioned, the total kinetic energy of the system may be expressed by the following equation: where, is the kinetic energy of the beam 1 and it is given by: is the kinetic energy of link 2 and it is given by: The same procedure can be applied to both: a mass concentrated at the joint 2 and a load with moment of inertia (J xp ) with respect to an axis normal to the plane of motion and through the center of gravity.By doing this for a mass in the joint 2, the expression (13) can be modified to and 1 2 is the kinetic energy related to the concentrated mass that represents the servo-motor and the joint.The mass is located at point O 2 , in joint 2; Note that: 1) x 1 was replaced by L 1 , in which case we are not at any point on the first link, but at its end; 2) u 1 has been replaced by u 1E , because u 1E is the linear displacement at the end of the elastic beam 1.
Similarly too, ( ) Note, yet 1) The load is located at the end of the beam 2 and therefore, 2 x was replaced by 2 L .2) 2 u and 2 u  were replaced, respectively by 2 E u and 2 E u  , emphasizing that the displacement considered due to the flexibility at the end of the beam 2.
and the moment of inertia of the load with respect to an axis through of the point O 2 defined by J p , then is the kinetic energy related to the mass of the load and the rotational kinetic energy of the load due to its movement around the axis that passes through point O 2 and that is perpendicular to the plane shown in Figure 1.
With these facts in mind, one can obtain the following final expression for the total kinetic energy of the system: m L q q c q q nq q nq q c nq q nq q q q s m L q φ θ θ θ θ nq q s nq q nq q q q c nq q nq q q q s m L m L q q m q q m φ θ φ   q q q q c m q q q q s J q q θ θ φ φ φ φ θ θ φ φ φ φ φ θ φ φ

Potential Energy
In the calculation of the total potential energy of the system it is assumed an energy associated to the rigid movement (gravitational potential energy), plus the elastic potential energy of the links.Ox is taken as reference and the potential energy of the system (assuming u 1 and u 2 sufficiently small) is given by: where: g is the gravitational acceleration constant; L 1 and L 2 are the lengths of the links 1 and 2, respectively; EI 1 and EI 2 are the rigidity of the links 1 and 2, respectively, which in this model are assumed to be constants.In fact, the elastic displacements u 1 and u 2 (Figure 2 and Figure 3), were considered sufficient small.However, I would like to outline that the assumption that the deformation is sufficient small do not implies, necessarily, that the angles are small.

Motion Equations
In order to write down the motion equations of the system it will be used the assumed modes method.For the elastic displacements of the beams 1 and 2, it is possible to assume that ( ) ( ) where the admissible functions O x y .Hence the system becomes represented by ( ) n + degrees of freedom.If it is assumed that the amplitudes of the higher order vibration modes are very small when compared to the first vibration mode, the system may be truncated with n equal to 2, resulting in a problem involving six degrees of freedom.
Besides this, if it is assumed that φ are eigenfunctions of the problem of a clamped free beam and due to the ortogonality of this functions we will have: With this equation, one can assess each one of the integrals in the equations for both the kinetic and the potential energy considering 1 2 11 12 21 22 , , , , , q q q q θ θ the generalized coordinates and 1 2 , Γ Γ the non-conservative torques acting at the joint of the system, it is possible to write the equations of the movement-with convenient parameterization of the terms and without the usage linearization procedure-using the Lagrange equations for non-conservative systems, ( ) where: Q r are the time dependent generalized non-conservative forces (or torques).
The equations assume the final form, where )
The coefficients of these equations are presented as follows: general coefficients are relative of the rigid system.Table 1 and Table 2 show the specific coefficients for the elastic system.Thus, General coefficients: ( ) ( ) ( ) ( ) ( ) ( )   ( ) and, q q q q q q q q nq q nq q nq q nq q q q q q nq q nq q This convenient parameterization of the terms of equations of motion allows us a comparison with pre-established methods and the contribution of each elastic term of the parameterized system.

Some Numerical Simulation Results
The simulations are performed considering sinusoidal excitation.Some results are presented, according to the following methodology: 1) The elastic system was simulated with all its contributions taken into account.Full elastic system, equations simulated (28) to (33) and Figure 4 to Figure 7.
2) After that, the effects are individually and cumulatively subtracted and the system behavior is analyzed.Elastic system, Equations (28) to (33) subtracted flexibility in the J-term; (Figure 8 and Figure 9) and after subtracted flexibility in the perturb T term (Figure 10 and Figure 11).
3) The effects are subtracted until the limit condition in which the elastic system is reduced (mathematically) to a rigid one by means of vanishing the flexibility related terms, and the system response converges-as expected-for the case of the rigid system modeled separately.Equations reduced (28) to (33) and Figure 12 and

Conclusions
The Lagrangian dynamic in conjunction with the assumed modes method were utilized in order to obtain the non-linear equations of motion.
A convenient parameterization of the terms of the motion equations, which makes it easier to compare the simulation results for the rigid and for the elastic system, was also developed.
About the comparison with other equations of motion in current literature, I outline that the approach here formulated is innovative.The works founded in the literature involving dynamic modelling of the elastic structures uses some kind of the linearization procedure and presents your equations of motion in most cases like a matrices formulation, i.e., mass matrix, damping matrix, stiffness matrix and elasticity matrix.In this paper, equations are treated with all non-linearities taken into account and the approach entails parameterization of the nonlinear dynamics without the usage of any simplifying linearization procedure.The elastic structure may be mathematically reduced to a rigid one by means of vanishing the flexibility related terms.The same procedure may be extended to the simulations, which makes it possible to find a mathematical frontier between both systems.
The effects of the flexibility are explored by comparing the resulting simulation results.Besides, the way in which the motion equations are treated in this paper allows the monitoring of each con-
Figure 2 shows
geometrical boundary conditions with respect to the representation of the links in the reference systems [ ]

Table 1 .
Coefficients of the first three equations of motion.

Table 2 .
Coefficients of the last three equations of motion.