Some Properties of Motion Equations Describing the Nonlinear Dynamical Response of a Multibody System with Flexible Elements

The industrial applications use instruments and machines operating at high speeds, developing high forces, low temperatures, corrosive environments, extreme pressures, and so forth. Under these conditions, the elasticity of elements such a machine is built of cannot be ignored anymore, and models are needed to more accurately “grasp” the mechanical phenomena accompanying the operation. The vibrations and the loss of stability are the main effects occurring under these conditions. For the study on this kind of systems with rigid motion and elastic elements, numerous models have been elaborated, the main idea being the discretization of the elements and the use of finite element method. Finally, second-order differential equations with variable coefficients are obtained; these equations are strong nonlinear ones due to the time-dependent values of angular speed and acceleration, and they can be linearized considering a very short period of time, inwhich the motion is considered to be “frozen.” The aim of this paper is to present some characteristic properties of these systems.


Introduction
A mechanical system, a machine or instrument, is made up of elastic elements, the elasticity manifesting itself more or less.The rigid elements assumption generally made when studying such technical systems represents a first approximation leading to rapid results closer to reality.Depending on the given instrument operation conditions, this assumption may lead to correct results or to results considerably deviating from the real situation.If the instrument or the machine works with low operation speeds or if it is subjected to lower loads, then the model built based on the rigid elements assumption may lead to excellent results.The elasticity becomes a significant element if the loads occurred are high and/or the operation speeds are high Figure 1 .In this case, the deformations of the machine element will influence, usually in a negative way, the correct operation of the system.The resonance and the loss of stability represent the main forms of manifestation of elasticity.They will occur in the case of an inadequate design leading to a fast machine damage.The main method of approach of such a problem is the method of finite elements, a method used in a lot of works for elaborating models describing the behavior of machines containing elastic elements.
As a result of modeling a set of differential equations having a complex nonlinear form is obtained.It contains a series of additional terms due to the relative movements of nodal coordinates and the general motion, and the operation motion is called "rigid motion" of the system.As the geometrical configuration of a mechanical system with rigid motion changes from one moment to the other, the motion equations are valid for that very moment, in which the system motion is to be considered to be "frozen".In this case, Δt periods of time, when the coefficients can be approximate as constant, are being considered.It is difficult to estimate the length of this interval Δt that depends on the topology and geometry of the multibody system.We consider that it is short enough in order to regard the motion as "frozen."The analysis made is an incremental analysis.The motion equations written in the local coordinate system may be expressed under the form 1-9 m e δe 2c e ω δe k e k e ε k e ω 2 k G e δ e q e q * e − q i e ε − q i e ω 2 − m i Ee Iε L − m i oe R T roG , where m e : inertia matrix, symmetric; c e ω : matrix of Coriolis terms, skew symmetric; k e : classical rigidity matrix, symmetric; k e ε : modification of rigidity due to the angular acceleration of the mobile reference system; k e ω 2 : modification of rigidity due to the angular speed of the mobile reference system; k G e : rigidity due to the second-order effects; q e q * e : the external forces concentrated and distributed; q i e ε : inertial forces due to the angular acceleration of the mobile reference system; q i e ω 2 : inertial forces due to the angular speed of the mobile reference system; m i Ee Iε L : inertial forces due to rotation only at bars; m i oe R T roG : inertial forces due to the shifting of the mobile reference system.After passing to a common global coordinate system neglecting the second-order effects , finally, the second-order differential equations system with variable coefficients is obtained 10-12 where M: inertia matrix, symmetric; C: matrix of Coriolis terms, skew symmetric; K: rigidity matrix, symmetric; K ε : modification of rigidity due to the angular acceleration of the mechanical system components, skew symmetric; K ω 2 : modification of rigidity due to the angular speed of the mechanical system components, symmetric; Q: matrix of the total nodal loads.
To summarize, matrices M, K, K ω 2 of the whole structure are symmetric and C, K ε are skew symmetric.New researches increase the complexity of the problem of the study of the multibody systems with flexible elements 13-17 , but there are not many useful results concerning solving the equations.Some properties of such a system will be presented later.

Properties of Motion Equations of the Mechanical System
The motion equations of such system have properties allowing an easier solution of the equation system obtained but also a qualitative interpretation of the dynamic response of the system.We are presenting these properties as follows.

P1. In the Rayleigh Quotient, the Eigenvalues Do Not Depend Directly on the Damping Matrix
We consider a very short period of time in which the variation of the matrix coefficients of motion equations is of no importance very low .The system solution of the following form is selected: By derivation, the following is being successively obtained: Introducing into the homogeneous system associated to 1.2 , we get By premultiplying 2.3 by A, we get

2.4
Since C and K ε are skew-symmetric, we have and 2.4 becomes where see Rayleigh 18 Moreover, these values are real, since K K ω 2 and M are symmetric; thus there is no damping in the system.It results that the matrix C, having in equations the significance of a viscous damping, does not introduce a damping term in the obtained solutions.Its skew symmetry denotes the system property not dissipating the energy to the exterior.

Let us now write the relation 2.1 as follows:
Δ A 1 cos ωt A 2 sin ωt.

2.1
We can enounce the following property.

P2. Matrix A Defined by Relation 2.11 Has Two Null Eigenvalues
The matrices K ε and K ω 2 are neglected; they have generally a small contribution to the rigidity matrix.We propose ourselves to solve the homogeneous system associated to 1.2 which becomes under this assumption with M and K symmetric and C skew-symmetric.With the classical substitution the first-order linear system is obtained where and the unit matrix was named E.
If n is the system dimension the dimension of matrix A will be 2n × 2n.

2.12
The matrix of eigenvectors is clear.By transforming x Qq, the system is reduced to q λq, 2.13 where

2.14
The equations qi λ i q i allow the solutions q i C i exp λ i t .
If in the condition det A − λE 0 we plug in λ 0, we get the condition det A 0. As a proof, Developing on rows, we get finally because K is singular.We will show in what follows that A allows another eigenvalue 0. If an eigenvalue λ 1 is known, matrix A can be reduced to a matrix of dimension 2n − 1 × 2n − 1 which has 2n − 1 eigenvalues equal to the eigenvalues of A where λ 1 has been eliminated.It is known 19 that matrix B T −1 AT has the same eigenvalues as A, but the eigenvectors are Then, if a matrix B has the form

2.20
By simple computation, we get

2.21
It results in the possibility of avoiding the computation of the eigenvalues for matrix A with det A 0, considering the fact that we know an eigenvalue λ 1 0, and, by the transformation presented, we will compute the other eigenvalues as eigenvalues of B 22 .
For λ 1 0, the eigenvector V 1 corresponding to matrix A will be computed.We will have . . .

2.22
The n first rows give us e n 1 0; e n 2 0, . . ., e 2n 0. Using these results, the other n rows give us . .
or, carrying out the multiplications of the matrix blocks, . .

2.23
The result will be that the first n components of the eigenvector are components of the vector corresponding to value λ 1 0 for matrix −M −1 K which characterizes the system M Δ KΔ 0.

2.24
With these components, we build where e 1 1.We compute B T −1 AT.
where α * is the matrix having the dimensions 2n× 2n−1 .It is obtained from −M −1 K having the dimension 2n × 2n by eliminating the first column.We have

2.29
We have

2.30
Taking into account the previous relation, we may write     To these two eigenvalues corresponds the nonharmonic solution q C 1 t C 2 which will represent the rigid motion of the multibody system in a first approximation C 1 and C 2 are integration constants .The other values are different from null for a multibody system with only one degree of mobility generally being different from each other.

Conclusions
The paper presents a few properties of motion equations in the case of mechanical systems having elastic elements.The properties are due to the existence of the skew-symmetric matrix C, by which the relative motion of nodal coordinates is manifested in case of applying the finite element method by the Coriolis effects.These properties allow a qualitative analysis of the obtained motion equations.Thus, the Coriolis effects due to the relative motions will determine a modification generally small of the systems eigenvalues.The Coriolis damping is not dissipative, meaning that the systems energy is not being influenced by the terms in which the skew-symmetric matrix C occurs.The last conclusion is that in case of modeling with the finite element method, the rigid motion, considered in a first approximation an a uniform motion of the system, can be eliminated from the motion equations written.
This fact also suggests the incremental solution of the problem on small periods of time in which the motion may be considered as "frozen" or uniform.
written λ − λ 1 det B 22 − λE 0 resuming to determine the eigenvalues of B 22 .It is checked by direct computation that the transformation T is as follows: