Advances in Nonlinear VibrationView this Special Issue
Research Article | Open Access
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, in which the motion is considered to be “frozen.” The aim of this paper is to present some characteristic properties of these systems.
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, periods of time, when the coefficients can be approximate as constant, are being considered. It is difficult to estimate the length of this interval 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] where : inertia matrix, symmetric; : matrix of Coriolis terms, skew symmetric; : classical rigidity matrix, symmetric; : modification of rigidity due to the angular acceleration of the mobile reference system; : modification of rigidity due to the angular speed of the mobile reference system; : rigidity due to the second-order effects; : the external forces concentrated and distributed; : inertial forces due to the angular acceleration of the mobile reference system; : inertial forces due to the angular speed of the mobile reference system; : inertial forces due to rotation only at bars; : 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 : inertia matrix, symmetric; : matrix of Coriolis terms, skew symmetric; : rigidity matrix, symmetric; : modification of rigidity due to the angular acceleration of the mechanical system components, skew symmetric; : modification of rigidity due to the angular speed of the mechanical system components, symmetric; : matrix of the total nodal loads.
To summarize, matrices ,, of the whole structure are symmetric and 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.
2. 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.
2.1. 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 , we get
Since and are skew-symmetric, we have and (2.4) becomes where (see Rayleigh ) Moreover, these values are real, since and are symmetric; thus there is no damping in the system. It results that the matrix , 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: We can enounce the following property.
2.2. P2. Matrix A Defined by Relation (2.11) Has Two Null Eigenvalues
The matrices and 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 and symmetric and skew-symmetric. With the classical substitution the first-order linear system is obtained or where and the unit matrix was named .
If is the system dimension the dimension of matrix will be .
Let be the eigenvalues for matrix and the corresponding eigenvectors. We write The matrix of eigenvectors is clear. By transforming , the system is reduced to where
The equations allow the solutions .
If in the condition we plug in , we get the condition det .
As a proof, Developing on rows, we get finally because is singular.
We will show in what follows that allows another eigenvalue 0. If an eigenvalue is known, matrix can be reduced to a matrix of dimension which has eigenvalues equal to the eigenvalues of where has been eliminated. It is known  that matrix has the same eigenvalues as , but the eigenvectors are Then, if a matrix has the form with then the problem of eigenvalues is to be written resuming to determine the eigenvalues of .
It is checked by direct computation that the transformation is as follows:
By simple computation, we get
It results in the possibility of avoiding the computation of the eigenvalues for matrix with det = 0, considering the fact that we know an eigenvalue , and, by the transformation presented, we will compute the other eigenvalues as eigenvalues of .
For , the eigenvector corresponding to matrix will be computed. We will have or The first rows give us .
Using these results, the other rows give us or, carrying out the multiplications of the matrix blocks, The result will be that the first components of the eigenvector are components of the vector corresponding to value for matrix which characterizes the system With these components, we build where . We compute .
Consider where is the matrix having the dimensions . It is obtained from having the dimension by eliminating the first column.
We have where the matrix columns have been named .
By developing the matrix determinant on rows, we obtain (we took into account that ).
The following relation has been considered: We have Taking into account the previous relation, we may write since with , the lines of the matrix have been named. It results in det , thus, the matrix has a null eigenvalue. Matrix will have two null eigenvalues.
To these two eigenvalues corresponds the nonharmonic solution which will represent the rigid motion of the multibody system in a first approximation ( and 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.
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 , 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 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.
This paper is supported by the Sectoral Operational Programme Human Resources Development (SOP HRD), financed from the European Social Fund, and by the Romanian Government under contract no. POSDRU/89/1.5/S/59323.
- S. Vlase, “Dynamical response of a multibody system with flexible elements with a general three-dimensional motion,” Romanian Journal of Physics, vol. 57, pp. 3676–4693, 2012 (Romanian).
- C. Bagci, “Elastodynamic response of mechanical systems using matrix exponential mode uncoupling and incremental forcing techniques with finite element method,” in Proceedings of the 6th Word Congress on Theory of Machines and Mechanisms, p. 472, New Delhi, India, December 1983.
- B. M. Bahgat and K. D. Willmert, “Finite element vibrational analysis of planar mechanisms,” Mechanism and Machine Theory, vol. 11, no. 1, pp. 47–71, 1976.
- W. L. Cleghorn, R. G. Fenton, and B. Tabarrok, “Finite element analysis of high-speed flexible mechanisms,” Mechanism and Machine Theory, vol. 16, no. 4, pp. 407–424, 1981.
- A. G. Erdman, G. N. Sandor, and A. Oakberg, “A general method for kineto-elastodynamic analysis and synthesis of mechanisms,” ASME Journal of Engineering for Industry, vol. 94, no. 4, p. 1193, 1972.
- P. K. Nath and A. Ghosh, “Kineto-elastodynamic analysis of mechanisms by finite element method,” Mechanism and Machine Theory, vol. 15, no. 3, pp. 179–197, 1980.
- S. Vlase, “A method of eliminating lagrangean multipliers from the equations of motion of interconnected mechanical systems,” Journal of Applied Mechanics, vol. 54, no. 1, pp. 235–237, 1987.
- S. Vlase, “Elimination of lagrangian multipliers,” Mechanics Research Communications, vol. 14, no. 1, pp. 17–22, 1987.
- S. Vlase, “Finite element analysis of the planar mechanisms: numerical aspects,” in Applied Mechanics, vol. 4, pp. 90–100, Elsevier, New York, NY, USA, 1992.
- P. Fanghella, C. Galletti, and G. Torre, “An explicit independent-coordinate formulation for the equations of motion of flexible multibody systems,” Mechanism and Machine Theory, vol. 38, no. 5, pp. 417–437, 2003.
- J. Mayo and J. Domínguez, “Geometrically non-linear formulation of flexible multibody systems in terms of beam elements: geometric stiffness,” Computers and Structures, vol. 59, no. 6, pp. 1039–1050, 1996.
- X. Zhang and A. G. Erdman, “Dynamic responses of flexible linkage mechanisms with viscoelastic constrained layer damping treatment,” Computers and Structures, vol. 79, no. 13, pp. 1265–1274, 2001.
- W. Hou and X. Zhang, “Dynamic analysis of flexible linkage mechanisms under uniform temperature change,” Journal of Sound and Vibration, vol. 319, no. 1-2, pp. 570–592, 2009.
- G. Piras, W. L. Cleghorn, and J. K. Mills, “Dynamic finite-element analysis of a planar high-speed, high-precision parallel manipulator with flexible links,” Mechanism and Machine Theory, vol. 40, no. 7, pp. 849–862, 2005.
- Y. M. Shi, Z. F. Li, H. X. Hua, Z. F. Fu, and T. X. Liu, “The modelling and vibration control of beams with active constrained layer damping,” Journal of Sound and Vibration, vol. 245, no. 5, pp. 785–800, 2001.
- B. S. Thompson and C. K. Sung, “A survey of finite element techniques for mechanism design,” Mechanism and Machine Theory, vol. 21, no. 4, pp. 351–359, 1986.
- Z. Xianmin, L. Jianwei, and S. Yunwen, “Simultaneous optimal structure and control design of flexible linkage mechanism for noise attenuation,” Journal of Sound and Vibration, vol. 299, no. 4-5, pp. 1124–1133, 2007.
- C. M. Harris and C. E. Crede, Shock and Vibration Handbook Volume 1, McGraw-Hill, New York, NY, USA, 1961.
- C. H. Massonet et al., Computer Structures Calculus, Tehnică, Bucuresti, Romania, 1974.
Copyright © 2012 Maria Luminiţa Scutaru and Sorin Vlase. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.