scholarly journals Another Multibody Dynamics in Natural Coordinates through Automatic Differentiation and High-Index DAE Solving

2020 ◽  
Vol 24 (3) ◽  
pp. 315-341
Author(s):  
John D Pryce ◽  
Nedialko Nedialkov
Keyword(s):  
Automatic Differentiation ◽  
Mass Matrix ◽  
Equations Of Motion ◽  
3D Models ◽  
High Index ◽  
Natural Coordinates ◽  
Differential Algebraic Equation ◽  
Sparse Models ◽  
Algorithmic Differentiation ◽  
A Minor

The Natural Coordinates (NCs) method for Lagrangian modelling and simulation of multi-body systems is valued for giving simple, sparse models. We describe our version of it (NPNCs) and compare with the classical ap- proach of Jalón and Bayo (JBNCs). NPNCs use the high-index differential- algebraic equation solver DAETS. Algorithmic differentiation, not symbolic algebra, forms the equations of motion from the Lagrangian. NPNCs give significantly smaller equation systems than JBNCs, at the cost of a non- constant mass matrix for fully 3D models—a minor downside in the DAETS context. A 2D and a 3D example are presented, with numerical results.

Natural Coordinates in the Optimal Control of Multibody Systems

2011 ◽  
Author(s):  
Peter Betsch ◽  
Ralf Siebert ◽  
Nicolas Sa¨nger
Keyword(s):  
Optimal Control ◽  
Multibody Systems ◽  
Mass Matrix ◽  
Equations Of Motion ◽  
Differential Algebraic Equations ◽  
Open Loop ◽  
Algebraic Equations ◽  
Natural Coordinates ◽  
Joint Forces ◽  
New Energy

The formulation of multibody dynamics in terms of natural coordinates (NCs) leads to equations of motion in the form of differential-algebraic equations (DAEs). A characteristic feature of the natural coordinates approach is a constant mass matrix. The DAEs make possible (i) the systematic assembly of open-loop and closed-loop multibody systems, (ii) the design of state-of-the-art structure-preserving integrators such as energy-momentum or symplectic-momentum schemes, and (iii) the direct link to nonlinear finite element methods. However, the use of NCs in the optimal control of multibody systems presents two major challenges. First, the consistent application of actuating joint-forces becomes an issue since conjugate joint-coordinates are not directly available. Secondly, numerical methods for optimal control with index-3 DAEs are still in their infancy. The talk will address the two aforementioned issues. In particular, a new energy-momentum consistent method for the optimal control of multibody systems in terms of NCs will be presented.

Natural Coordinates in the Optimal Control of Multibody Systems

2011 ◽  
Vol 7 (1) ◽  
Author(s):  
Peter Betsch ◽  
Ralf Siebert ◽  
Nicolas Sänger
Keyword(s):  
Optimal Control ◽  
Multibody Systems ◽  
Mass Matrix ◽  
Equations Of Motion ◽  
Differential Algebraic Equations ◽  
Open Loop ◽  
Algebraic Equations ◽  
Natural Coordinates ◽  
Joint Forces ◽  
New Energy

The formulation of multibody dynamics in terms of natural coordinates (NCs) leads to equations of motion in the form of differential-algebraic equations (DAEs). A characteristic feature of the natural coordinates approach is a constant mass matrix. The DAEs make possible (i) the systematic assembly of open-loop and closed-loop multibody systems, (ii) the design of state-of-the-art structure-preserving integrators such as energy-momentum or symplectic-momentum schemes, and (iii) the direct link to nonlinear finite element methods. However, the use of NCs in the optimal control of multibody systems presents two major challenges. First, the consistent application of actuating joint-forces becomes an issue since conjugate joint-coordinates are not directly available. Second, numerical methods for optimal control with index-3 DAEs are still in their infancy. The talk will address the two aforementioned issues. In particular, a new energy-momentum consistent method for the optimal control of multibody systems in terms of NCs will be presented.

Flexible Multibody Dynamics Formulation by Hamilton’s Equations

2010 ◽  
Author(s):  
Martin M. Tong
Keyword(s):  
Multibody Dynamics ◽  
Multibody System ◽  
Mass Matrix ◽  
Equations Of Motion ◽  
Flexible Multibody Dynamics ◽  
Flexible Body ◽  
Flexible Multibody System ◽  
Generalized Coordinates ◽  
Flexible Multibody ◽  
The Matrix

Numerical solution of the dynamics equations of a flexible multibody system as represented by Hamilton’s canonical equations requires that its generalized velocities q˙ be solved from the generalized momenta p. The relation between them is p = J(q)q˙, where J is the system mass matrix and q is the generalized coordinates. This paper presents the dynamics equations for a generic flexible multibody system as represented by p˙ and gives emphasis to a systematic way of constructing the matrix J for solving q˙. The mass matrix is shown to be separable into four submatrices Jrr, Jrf, Jfr and Jff relating the joint momenta and flexible body mementa to the joint coordinate rates and the flexible body deformation coordinate rates. Explicit formulas are given for these submatrices. The equations of motion presented here lend insight to the structure of the flexible multibody dynamics equations. They are also a versatile alternative to the acceleration-based dynamics equations for modeling mechanical systems.

scholarly journals NON-RIGID KINEMATIC EXCITATION FOR MULTIPLY-SUPPORTED SYSTEM WITH HOMOGENEOUS DAMPING

2018 ◽  
Vol 14 (4) ◽  
pp. 152-157
Author(s):  
Alexander G. Tyapin
Keyword(s):  
Mass Matrix ◽  
Equations Of Motion ◽  
Linear Equations ◽  
Internal Damping ◽  
Static Response ◽  
Kinematic Excitation ◽  
Hand Part ◽  
Internal Forces ◽  
The Right ◽  
Damping Model

This paper continues the discussion of linear equations of motion. The author considers non-rigid kinematic excitation for multiply-supported system leading to the deformations in quasi-static response. It turns out that in the equation of motion written down for relative displacements (relative displacements are defined as absolute displacements minus quasi-static response) the contribution of the internal damping to the load in some cases may be zero (like it was for rigid kinematical excitation). For this effect the system under consideration must have homogeneous damping. It is the often case, though not always. Zero contribution of the internal damping to the load is different in origin for rigid and non-rigid kinematic excitation: in the former case nodal loads in the quasi-static response are zero for each element; in the latter case nodal loads in elements are non-zero, but in each node they are balanced giving zero resulting nodal loads. Thus, damping in the quasi-static response does not impact relative motion, but impacts the resulting internal forces. The implementation of the Rayleigh damping model for the right-hand part of the equation leads to the error (like for rigid kinematic excitation), as damping in the Rayleigh model is not really “internal”: due to the participation of mass matrix it works on rigid displacements, which is impossible for internal damping

Conditions for the planes of symmetry of an elastically supported rigid body

2009 ◽  
Vol 223 (8) ◽  
pp. 1755-1766 ◽  
Author(s):  
S J Jang ◽  
Y J Choi
Keyword(s):  
Rigid Body ◽  
Eigenvalue Problems ◽  
Mass Matrix ◽  
Equations Of Motion ◽  
Mass Matrices ◽  
Out Of Plane ◽  
Modes Of Vibration ◽  
Compliance Matrices ◽  
Special Points ◽  
Elastically Supported

Introducing the planes of symmetry into an oscillating rigid body suspended by springs simplifies the complexity of the equations of motion and decouples the modes of vibration into in-plane and out-of-plane modes. There have been some research results from the investigation into the conditions for planes of symmetry in which prior conditions for the simplification of the equations of motion are required. In this article, the conditions for the planes of symmetry that do not need prior conditions for simplification are presented. The conditions are derived from direct expansions of eigenvalue problems for stiffness and mass matrices that are expressed in terms of in-plane and out-of-plane modes and the orthogonality condition with respect to the mass matrix. Two special points, the planar couple point and the perpendicular translation point are identified, where the expressions for stiffness and compliance matrices can be greatly simplified. The simplified expressions are utilized to obtain the analytical expressions for the axes of vibration of a vibration system with planes of symmetry.

Three-Dimensional Solid Brick Element Using Slopes in the Absolute Nodal Coordinate Formulation

2013 ◽  
Vol 9 (2) ◽  
Author(s):  
Alexander Olshevskiy ◽  
Oleg Dmitrochenko ◽  
Chang-Wan Kim
Keyword(s):  
Finite Element ◽  
Finite Elements ◽  
Absolute Nodal Coordinate Formulation ◽  
Mass Matrix ◽  
Equations Of Motion ◽  
Test Problems ◽  
Absolute Nodal Coordinate ◽  
Highly Nonlinear ◽  
Brick Element ◽  
The Absolute

The present paper contributes to the field of flexible multibody systems dynamics. Two new solid finite elements employing the absolute nodal coordinate formulation are presented. In this formulation, the equations of motion contain a constant mass matrix and a vector of generalized gravity forces, but the vector of elastic forces is highly nonlinear. The proposed solid eight node brick element with 96 degrees of freedom uses translations of nodes and finite slopes as sets of nodal coordinates. The displacement field is interpolated using incomplete cubic polynomials providing the absence of shear locking effect. The use of finite slopes describes the deformed shape of the finite element more exactly and, therefore, minimizes the number of finite elements required for accurate simulations. Accuracy and convergence of the finite element is demonstrated in nonlinear test problems of statics and dynamics.

A Recursive Implementation Method With Implicit Integrator for Multibody Dynamics

2001 ◽  
Author(s):  
H. J. Cho ◽  
H. S. Ryu ◽  
D. S. Bae ◽  
J. H. Choi ◽  
B. Ross
Keyword(s):  
Equations Of Motion ◽  
Flexible Body ◽  
Differential Algebraic Equation ◽  
Dynamic Correction ◽  
Body Dynamics ◽  
Mode Method ◽  
Numerical Integrators ◽  
Flexible Body Dynamics ◽  
Body Concept ◽  
Implementation Method

Abstract A recursive implementation method for the equations of motion and kinematics is presented. Computational structure of the kinematic and dynamic equations is exploited to systematically implement a dynamic analysis program RecurDyn. A differential algebraic equation solution method with implicit numerical integrators is discussed. Virtual body concept is introduced for the flexible body dynamics. The accuracy of the flexible body solutions is estimated by an error measure and is improved by the dynamic correction mode method. Several examples are solved to demonstrate the efficiency of the proposed methods.

Modeling and simulation of planar flexible link manipulators with rigid tip connections to revolute joints

Robotica ◽  
2004 ◽  
Vol 22 (3) ◽  
pp. 285-300 ◽  
Author(s):  
S. M. Megahed ◽  
K. T. Hamza
Keyword(s):  
Mass Matrix ◽  
Equations Of Motion ◽  
Flexible Link ◽  
Link Length ◽  
Finite Segment ◽  
Revolute Joints ◽  
Body Dynamics ◽  
Lumped Masses ◽  
Multi Body ◽  
Rigid Zones

This paper presents the basis of a mathematical model for simulation of planar flexible-link manipulators, taking into consideration the effect of higher stiffness zones at the link tips. The proposed formulation is a variation of the finite segment multi-body dynamics approach. The formulation employs a consistent mass matrix in order to provide better approximation than the traditional lumped masses often encountered in the finite segment approach. The formulation is implemented into a computational code and tested through three examples; cantilever beam, rotating beam and three-link manipulator. In these examples, the length of the rigid tips at both sides of each link ranges from 0% to 6.25% of the whole link length. The zones of higher stiffness at the link tips are treated as short rigid zones. The effect of the rigid zones is averaged along with some portions of the flexible links, thereby allowing further simplification of the dynamic equations of motion. The simulation results demonstrate the effectiveness of the proposed modeling technique and show the importance of not ignoring the effect of the rigid tips.

ASPECTS OF Q-MATTER STABILITY

1992 ◽  
Vol 07 (28) ◽  
pp. 7169-7184 ◽  
Author(s):  
MINOS AXENIDES
Keyword(s):  
Order Parameter ◽  
Mass Matrix ◽  
Equations Of Motion ◽  
Finite Energy ◽  
Field Theories ◽  
Strong Interactions ◽  
Mass Relation ◽  
Bosonic Field ◽  
Meson Octet ◽  
Continuous Symmetries

Relativistic bosonic field theories in 3+1 dimensions with exact global continuous symmetries and conserved charges Q may admit stable, finite energy, time dependent configurations (Q-balls) as solutions to their equations of motion. Previous work established their existence for both Abelian and non-Abelian symmetries. In the present work we elaborate on some more issues of stability and uniqueness that arise in the SO(3) and SU(3) renormalizable models. We consider the effect of explicit symmetry breaking in the spectrum of the SU(3) model, by identifying its order parameter with the meson octet and by choosing a mass matrix consistent with the Gell-Mann-Okubo mass relation. We demonstrate the existence of “isospin” and “strange” balls whose stability is due to the presence of residual global symmetries which are identified with the exact symmetries of isospin and strangeness of strong interactions.

Sign in / Sign up

Export Citation Format

Share Document