A Modal Analysis Solution Technique to the Equations of Motion for Elastic Mechanism Systems Including the Rigid-Body and Elastic Motion Coupling Terms

1993 ◽  
Vol 115 (2) ◽  
pp. 314-323 ◽  
Author(s):  
A. G. Jablokow ◽  
S. Nagarajan ◽  
D. A. Turcic

This paper presents a modal analysis solution technique for the matrix equations of motion of elastic mechanism systems including the rigid-body elastic motion coupling terms and general damping. In many cases, researchers have neglected these terms because they complicate and diminish the efficiency of the solution process. This has been justified by assuming that the terms are small and do not affect the system response. The results obtained using the techniques developed within show this is true for some, but not all mechanisms. The solution technique adds the rigid-body elastic motion coupling terms to the system mass, damping, and stiffness matrices thus allowing them to affect the system response appropriately. The resulting nonsymmetric matrices are then rewritten in first order form allowing general damping to be included in the analysis. Modal analysis techniques are utilized to solve for the steady-state response of the elastic mechanism system. Thus the methods developed in this work provide a technique for including the rigid-body elastic motion coupling terms and general damping in the equations of motion while maintaining the advantage of using efficient modal analysis techniques of finding the response of the system. A number of examples are presented that establish the validity of this approach to the solution of the matrix equations of motion for elastic mechanism systems. The results show that the rigid-body elastic motion coupling terms can become significant at higher operating speeds.

Author(s):  
Simon A. Neild ◽  
Andrea Cammarano ◽  
David J. Wagg

In this paper we discuss a theoretical technique for decomposing multi-degree-of-freedom weakly nonlinear systems into a simpler form — an approach which has parallels with the well know method for linear modal analysis. The key outcome is that the system resonances, both linear and nonlinear are revealed by the transformation process. For each resonance, parameters can be obtained which characterise the backbone curves, and higher harmonic components of the response. The underlying mathematical technique is based on a near identity normal form transformation. This is an established technique for analysing weakly nonlinear vibrating systems, but in this approach we use a variation of the method for systems of equations written in second-order form. This is a much more natural approach for structural dynamics where the governing equations of motion are written in this form as standard practice. In fact the first step in the method is to carry out a linear modal transformation using linear modes as would typically done for a linear system. The near identity transform is then applied as a second step in the process and one which identifies the nonlinear resonances in the system being considered. For an example system with cubic nonlinearities, we show how the resulting transformed equations can be used to obtain a time independent representation of the system response. We will discuss how the analysis can be carried out with applied forcing, and how the approximations about response frequencies, made during the near-identity transformation, affect the accuracy of the technique. In fact we show that the second-order normal form approach can actually improve the predictions of sub- and super-harmonic responses. Finally we comment on how this theoretical technique could be used as part of a modal testing approach in future work.


1990 ◽  
Vol 112 (3) ◽  
pp. 307-314 ◽  
Author(s):  
Ye-Chen Pan ◽  
R. A. Scott ◽  
A. Galip Ulsoy

A dynamic model for flexible manipulators with prismatic joints is presented in Part I of this study. Floating frames following a nominal rigid body motion are introduced to describe the kinematics of the flexible links. A Lagrangian approach is used in deriving the equations of motion. The work done by the rigid body axial force through the axial shortening of the link due to transverse deformations is included in the Lagrangian function. Kinematic constraint equations are used to describe the compatibility conditions associated with revolute joints and prismatic joints, and incorporated into the equations of motion by Lagrange multipliers. The small displacements due to the flexibility of the links are then discretized by a displacement based finite element method. Equations of motion are derived for the cases of prescribed rigid body motion as well as prescribed joint torques/forces through application of Lagrange’s equations. The equations of motion and the constraint equations result in a set of differential algebraic equations. A numerical procedure combining a constraint stabilization method and a Newmark direct integration scheme is then applied to obtain the system response. An example, previously treated in the literature, is presented to validate the modeling and solution methods used in this study.


1990 ◽  
Vol 112 (1) ◽  
pp. 30-35 ◽  
Author(s):  
B. Samanta

A procedure is presented to study the dynamics of interconnected flexible systems using bond graphs. The concept of Lagrange multipliers, which are commonly used in analysis of constrained systems, is introduced in the procedure. The overall motions of each of the component bodies are described in terms of large rigid body motions and small elastic vibrations. Bond graphs are used to represent both rigid body and flexible dynamics of each body in a unified manner. Bond graphs of such sub-systems are coupled to one another satisfying the kinematic constraints at the interfaces to get the overall system model. Constraint reactions are introduced in the form of Lagrange multipliers at the interfaces without disturbing the integral causality in the subsystem models, which leads to easy derivation of system equations. The equations of motion and higher derivatives of the constraint relations are integrated to obtain the constraint reactions and the system response. The procedure is illustrated by an example system and results are in good agreement with those presented earlier.


1990 ◽  
Vol 112 (2) ◽  
pp. 215-224 ◽  
Author(s):  
S. Nagarajan ◽  
D. A. Turcic

The first step in the derivation of the equations of motion for general elastic mechanism systems was described in Part I of this work. The equations were derived at the elemental level using Lagrange’s equation and the generalized coordinates were both the rigid body degrees of freedom, and the elastic degrees of freedom of element ‘e’. Each rigid body degree of freedom gave rise to a scalar equation of motion, and the elastic degrees of freedom of element e gave rise to a vector equation of motion. Since both the rigid body degrees of freedom and elastic degrees of freedom are considered as generalized coordinates, the equations derived take into account the mutual dependence between the rigid body and elastic motions. This is important for mechanisms that are built using lightweight and flexible members and which operate at high speeds. A schematic diagram of how the equations of motion are obtained in this work is shown in Fig. 1 in Part I. The transformation step in the figure refers to the rotational transformation of the nodal elastic displacements (which were measured in the element coordinate system), so that they are measured in terms of the reference coordinate system. This transformation is necessary in order to ensure compatibility of the displacement, velocity and acceleration of the degrees of freedom that are common to two or more links during the assembly of the equations of motion. This final set of equations after assembly are obtained in closed form, and, given external torques and forces, can be solved for the rigid body and elastic response simultaneously taking into account the mutual dependence between the two responses.


1965 ◽  
Vol 7 (2) ◽  
pp. 185-192 ◽  
Author(s):  
P. Grootenhuis ◽  
D. J. Ewins

The equations of motion for a rigid body supported on four springs are derived for the general case of the centre-of-gravity being anywhere within the body and allowing for the sideways as well as the longitudinal stiffnesses of the springs. This constitutes a six-degrees-of-freedom case with three degrees of asymmetry. Coupling between motions in all directions occurs even when the centre-of-gravity is at the geometric centre with the exception then of vertical oscillations and rotation about the vertical axis. Any number of additional springs can be allowed for by adding terms to the expression for the potential energy stored in the springs. Allowance is made in the expression for kinetic energy for the products of inertia which arise with an offset centre-of-gravity. The real case is simulated for purposes of analysis by replacing the rigid body by a rectangular box with a light framework and all the mass concentrated at the eight corners. The matrix solution is changed into dimensionless parameters and the effect of an offset centre-of-gravity upon the eigenvalues and eigenvectors studied. Only the proportions of the box and the stiffness ratio between sideways to longitudinal stiffness of the springs remain as factors. The numerical example given is for proportions of height to width to length of 3/4/5 and for a stiffness ratio of 5. Small amounts of offset of the centre-of-gravity from the geometric centre do not alter the dynamic behaviour of the system much but displacing the total mass towards either a lower or an upper corner has marked effects. Some of the natural frequencies associated with motion in rotation when the system is symmetric become less than the frequencies connected with motion in translation for the centre-of-gravity being close to a corner connected to a spring. A large region free from any natural frequency arises when the centre-of-gravity is moved towards a corner furthest removed from the plane containing the springs. The asymptotic conditions for the position of the centre-of-gravity are also considered.


Author(s):  
Andreas Müller ◽  
Shivesh Kumar

AbstractDerivatives of equations of motion (EOM) describing the dynamics of rigid body systems are becoming increasingly relevant for the robotics community and find many applications in design and control of robotic systems. Controlling robots, and multibody systems comprising elastic components in particular, not only requires smooth trajectories but also the time derivatives of the control forces/torques, hence of the EOM. This paper presents the time derivatives of the EOM in closed form up to second-order as an alternative formulation to the existing recursive algorithms for this purpose, which provides a direct insight into the structure of the derivatives. The Lie group formulation for rigid body systems is used giving rise to very compact and easily parameterized equations.


1995 ◽  
Vol 62 (3) ◽  
pp. 685-691 ◽  
Author(s):  
F. Ma ◽  
T. K. Caughey

The coefficients of a linear nonconservative system are arbitrary matrices lacking the usual properties of symmetry and definiteness. Classical modal analysis is extended in this paper so as to apply to systems with nonsymmetric coefficients. The extension utilizes equivalence transformations and does not require conversion of the equations of motion to first-order forms. Compared with the state-space approach, the generalized modal analysis can offer substantial reduction in computational effort and ample physical insight.


2004 ◽  
Vol 126 (1) ◽  
pp. 175-183 ◽  
Author(s):  
E. P. Petrov

An effective method for analysis of periodic forced response of nonlinear cyclically symmetric structures has been developed. The method allows multiharmonic forced response to be calculated for a whole bladed disk using a periodic sector model without any loss of accuracy in calculations and modeling. A rigorous proof of the validity of the reduction of the whole nonlinear structure to a sector is provided. Types of bladed disk forcing for which the method may be applied are formulated. A multiharmonic formulation and a solution technique for equations of motion have been derived for two cases of description for a linear part of the bladed disk model: (i) using sector finite element matrices and (ii) using sector mode shapes and frequencies. Calculations validating the developed method and a numerical investigation of a realistic high-pressure turbine bladed disk with shrouds have demonstrated the high efficiency of the method.


2012 ◽  
Vol 2012 ◽  
pp. 1-6
Author(s):  
Xuefeng Duan ◽  
Chunmei Li

Based on the alternating projection algorithm, which was proposed by Von Neumann to treat the problem of finding the projection of a given point onto the intersection of two closed subspaces, we propose a new iterative algorithm to solve the matrix nearness problem associated with the matrix equations AXB=E, CXD=F, which arises frequently in experimental design. If we choose the initial iterative matrix X0=0, the least Frobenius norm solution of these matrix equations is obtained. Numerical examples show that the new algorithm is feasible and effective.


Author(s):  
Shanzhong Duan ◽  
Kurt S. Anderson

Abstract The paper presents a new hybrid parallelizable low order algorithm for modeling the dynamic behavior of multi-rigid-body chain systems. The method is based on cutting certain system interbody joints so that largely independent multibody subchain systems are formed. These subchains interact with one another through associated unknown constraint forces f¯c at the cut joints. The increased parallelism is obtainable through cutting the joints and the explicit determination of associated constraint loads combined with a sequential O(n) procedure. In other words, sequential O(n) procedures are performed to form and solve equations of motion within subchains and parallel strategies are used to form and solve constraint equations between subchains in parallel. The algorithm can easily accommodate the available number of processors while maintaining high efficiency. An O[(n+m)Np+m(1+γ)Np+mγlog2Np](0<γ<1) performance will be achieved with Np processors for a chain system with n degrees of freedom and m constraints due to cutting of interbody joints.


Sign in / Sign up

Export Citation Format

Share Document