A Generalized Component Mode Synthesis Approach for Multibody System Dynamics Leading to Constant Mass and Stiffness Matrices

Author(s):  
Johannes Gerstmayr ◽  
Astrid Pechstein

A standard technique to reduce the system size of flexible multibody systems is the component mode synthesis. Selected mode shapes are used to approximate the flexible deformation of each single body numerically. Conventionally, the (small) flexible deformation is added relatively to a body-local reference frame, which results in the floating frame of reference formulation (FFRF). The coupling between large rigid body motion and small relative deformation is nonlinear, which leads to computationally expensive non-constant mass matrices and quadratic velocity vectors. In the present work, the total (absolute) displacements are directly approximated by means of mode shapes, without a splitting into rigid body motion and superimposed flexible deformation. As the main advantage of the proposed method, the mass matrix is constant, the quadratic velocity vector vanishes and the stiffness matrix is a co-rotated constant matrix. Numerical experiments show the equivalence of the proposed method to the FFRF approach.

Author(s):  
Astrid Pechstein ◽  
Daniel Reischl ◽  
Johannes Gerstmayr

A standard technique to reduce the system size of flexible multibody systems is the component mode synthesis. Selected mode shapes are used to approximate the flexible deformation of each single body numerically. Conventionally, the (small) flexible deformation is added relatively to a body-local reference frame which results in the floating frame of reference formulation (FFRF). The coupling between large rigid body motion and small relative deformation is nonlinear, which leads to computationally expensive nonconstant mass matrices and quadratic velocity vectors. In the present work, the total (absolute) displacements are directly approximated by means of global (inertial) mode shapes, without a splitting into rigid body motion and superimposed flexible deformation. As the main advantage of the proposed method, the mass matrix is constant, the quadratic velocity vector vanishes, and the stiffness matrix is a co-rotated constant matrix. Numerical experiments show the equivalence of the proposed method to the FFRF approach.


Author(s):  
Pierangelo Masarati ◽  
Fanny Darbas ◽  
Israël Wander

Abstract Substructuring, or component mode synthesis, requires components to share interface regions. When components modeled with rather different, often incompatible levels of refinement need to be connected, correctly defining the interfaces may be important. This work proposes the definition of the reduction of interface regions to the equivalent rigid-body motion which minimizes the strain energy in the structural component. The proposed formulation provides a natural and physically sound solution for the connection of detailed structural components within coarse, multi-rigid-body and 1D flexible models.


Author(s):  
Pascal Ziegler ◽  
Alexander Humer ◽  
Astrid Pechstein ◽  
Johannes Gerstmayr

In industrial practice, the floating frame of reference formulation (FFRF)—often combined with the component mode synthesis (CMS) in order to reduce the number of flexible degrees-of-freedom—is the common approach to describe arbitrarily shaped bodies in flexible multibody systems. Owed to the relative formulation of the flexible deformation with respect to the reference frame, the equations of motion show state-dependent nonconstant inertia terms. Such relative description, however, comes along with considerable numerical costs, since both the mass matrix and gyroscopic forces, i.e., the quadratic velocity vector, need to be evaluated in every integration step. The state dependency of the inertia terms can be avoided by employing an alternative formulation based on the mode shapes as in the classical CMS approach. In this approach, which is referred to as generalized component mode synthesis (GCMS), the total (absolute) displacements are approximated directly. Consequently, the mass matrix is constant, no quadratic velocity vector appears, and the stiffness matrix is a corotated but otherwise constant matrix. In order to represent the same flexible deformation as in the classical FFRF-based CMS, however, a comparatively large number of degrees-of-freedom is required. The approach described in the present paper makes use of the fact that a majority of components in technical systems are constrained to motions showing large rotations only about a single spatially fixed axis. For this reason, the GCMS is adapted for multibody systems that are subjected to small flexible deformations and undergo a rigid body motion showing large translations, large rotations about one axis, but small rotations otherwise. Thereby, the number of shape functions representing the flexible deformation is reduced, which further increases numerical efficiency compared to the original GCMS formulation for arbitrary rotations.


1985 ◽  
Vol 107 (4) ◽  
pp. 431-439 ◽  
Author(s):  
A. A. Shabana

This paper is concerned with modeling inertia properties of flexible components that undergo large angular rotations. Consistent, lumped and hybrid mass techniques are presented in detail and used to model the inertia properties of flexible bodies. The consistent formulation allows using the finite-element method as well as Rayleigh-Ritz method to describe the deformation of elastic components. Lumped mass techniques allow using shape vectors or experimentally identified data. In the hybrid mass formulation, the flexibility mass matrix is evaluated using a consistent mass formulation, while the inertia coupling between gross rigid body motion and elastic deformation is formulated using a lumped mass technique. Different mass formulations require the evaluation of similar sets of time-invariant matrices that represent the inertia coupling. Consequently, these matrices have to be evaluated only once in advance for the dynamic analysis. A unified mathematical model, and accordingly a unified computer program (DAMS: Dynamic Analysis of Multibody Systems), that deal with different formulations are developed. A comparative study is presented in order to study the effect of the mass formulation on the dynamic response of elastodynamic constrained systems. The validity of the linear theory that neglects the effect of small oscillations on large rigid body motion is also discussed.


Author(s):  
C.-P. Fritzen ◽  
P. Gans ◽  
D. Wiese ◽  
M. Breitwieser

Abstract The following paper deals with a model-based method for unbalance identification. In contrast to conventional methods there is no need for the properties of the oil-film bearings to be known. First, the rotor is modeled as a conservative system in rigid bearings and the corresponding real mode shapes are calculated. The real motion of the rotor is represented by means of these mode shapes and the superimposed rigid-body motion of the rotor. The rigid-body motion is completely determined by the measured deflections in the oil-film bearings. Modal unbalances of the rotor are identified. Initial positive test results are presented in this paper.


Author(s):  
Guangbo Hao ◽  
Xianwen Kong ◽  
Xiuyun He

A planar reconfigurable linear (also rectilinear) rigid-body motion linkage (RLRBML) with two operation modes, that is, linear rigid-body motion mode and lockup mode, is presented using only R (revolute) joints. The RLRBML does not require disassembly and external intervention to implement multi-task requirements. It is created via combining a Robert’s linkage and a double parallelogram linkage (with equal lengths of rocker links) arranged in parallel, which can convert a limited circular motion to a linear rigid-body motion without any reference guide way. This linear rigid-body motion is achieved since the double parallelogram linkage can guarantee the translation of the motion stage, and Robert’s linkage ensures the approximate straight line motion of its pivot joint connecting to the double parallelogram linkage. This novel RLRBML is under the linear rigid-body motion mode if the four rocker links in the double parallelogram linkage are not parallel. The motion stage is in the lockup mode if all of the four rocker links in the double parallelogram linkage are kept parallel in a tilted position (but the inner/outer two rocker links are still parallel). In the lockup mode, the motion stage of the RLRBML is prohibited from moving even under power off, but the double parallelogram linkage is still moveable for its own rotation application. It is noted that further RLRBMLs can be obtained from the above RLRBML by replacing Robert’s linkage with any other straight line motion linkage (such as Watt’s linkage). Additionally, a compact RLRBML and two single-mode linear rigid-body motion linkages are presented.


Author(s):  
T. D. Burton ◽  
C. P. Baker ◽  
J. Y. Lew

Abstract The maneuvering and motion control of large flexible structures are often performed hydraulically. The pressure dynamics of the hydraulic subsystem and the rigid body and vibrational dynamics of the structure are fully coupled. The hydraulic subsystem pressure dynamics are strongly nonlinear, with the servovalve opening x(t) providing a parametric excitation. The rigid body and/or flexible body motions may be nonlinear as well. In order to obtain accurate ODE models of the pressure dynamics, hydraulic fluid compressibility must generally be taken into account, and this results in system ODE models which can be very stiff (even if a low order Galerkin-vibration model is used). In addition, the dependence of the pressure derivatives on the square root of pressure results in a “faster than exponential” behavior as certain limiting pressure values are approached, and this may cause further problems in the numerics, including instability. The purpose of this paper is to present an efficient strategy for numerical simulation of the response of this type of system. The main results are the following: 1) If the system has no rigid body modes and is thus “self-centered,” that is, there exists an inherent stiffening effect which tends to push the motion to a stable static equilibrium, then linearized models of the pressure dynamics work well, even for relatively large pressure excursions. This result, enabling linear system theory to be used, appears of value for design and optimization work; 2) If the system possesses a rigid body mode and is thus “non-centered,” i.e., there is no stiffness element restraining rigid body motion, then typically linearization does not work. We have, however discovered an artifice which can be introduced into the ODE model to alleviate the stiffness/instability problems; 3) in some situations an incompressible model can be used effectively to simulate quasi-steady pressure fluctuations (with care!). In addition to the aforementioned simulation aspects, we will present comparisons of the theoretical behavior with experimental histories of pressures, rigid body motion, and vibrational motion measured for the Battelle dynamics/controls test bed system: a hydraulically actuated system consisting of a long flexible beam with end mass, mounted on a hub which is rotated hydraulically. The low order ODE models predict most aspects of behavior accurately.


1997 ◽  
Vol 53 (6) ◽  
pp. 953-960 ◽  
Author(s):  
F. Belaj

The asymmetric units of both ionic compounds [N-(chloroformimidoyl)phosphorimidic trichloridato]trichlorophosphorus hexachlorophosphate, [ClC(NPCl3)2]+PCl^{-}_{6} (1), and [N-(acetimidoyl)phosphorimidic trichloridato]trichlorophosphorus hexachloroantimonate, [CH3C(NPCl3)2]+SbCl^{-}_{6} (2), contain two formula units with the atoms located on general positions. All the cations show cis–trans conformations with respect to their X—C—N—P torsion angles [X = Cl for (1), C for (2)], but quite different conformations with respect to their C—N—P—Cl torsion angles. Therefore, the two NPCl3 groups of a cation are inequivalent, even though they are equivalent in solution. The very flexible C—N—P angles ranging from 120.6 (3) to 140.9 (3)° can be attributed to the intramolecular Cl...Cl and Cl...N contacts. A widening of the C—N—P angles correlates with a shortening of the P—N distances. The rigid-body motion analysis shows that the non-rigid intramolecular motions in the cations cannot be explained by allowance for intramolecular torsion of the three rigid subunits about specific bonds.


Sign in / Sign up

Export Citation Format

Share Document