An Ordinary Differential Equation Formulation for Multibody Dynamics: Holonomic Constraints

2016 ◽  
Vol 16 (2) ◽  
Author(s):  
Edward J. Haug
Keyword(s):  
Mass Matrix ◽  
Equations Of Motion ◽  
Computational Cost ◽  
State Variables ◽  
Generalized Coordinates ◽  
Constraint Manifold ◽  
Local Coordinates ◽  
Numerical Integrators ◽  
Well Conditioning ◽  
Spinning Top

A method is presented for formulating and numerically integrating ordinary differential equations (ODEs) of motion for holonomically constrained multibody systems. Tangent space coordinates are defined as independent generalized coordinates that serve as state variables in the formulation, yielding ODEs of motion. Orthogonal dependent coordinates are used to enforce kinematic constraints at position, velocity, and acceleration levels. Criteria that assure accuracy of constraint satisfaction and well conditioning of the reduced mass matrix in the equations of motion are used as the basis for redefining local coordinates on the constraint manifold, as needed, transparent to the user and at minimal computational cost. The formulation is developed for holonomically constrained multibody models that are based on essentially any form of generalized coordinates. A spinning top with Euler parameter orientation coordinates is used as a model problem to analytically reduce Euler's equations of motion to ODEs. Numerical results using a fourth-order Nystrom integrator verify that accurate results are obtained, satisfying position, velocity, and acceleration constraints to computer precision. A computational algorithm for implementing the approach with state-of-the-art explicit numerical integrators is presented and used in solution of three examples, one planar and two spatial. Performance of the method in satisfying all three forms of kinematic constraint, based on error tolerances embedded in the formulation, is verified.

An Ordinary Differential Equation Formulation for Multibody Dynamics: Nonholonomic Constraints

2016 ◽  
Vol 17 (1) ◽  
Author(s):  
Edward J. Haug
Keyword(s):  
Differential Equations ◽  
Ordinary Differential Equations ◽  
Multibody Systems ◽  
Mass Matrix ◽  
Equations Of Motion ◽  
Computational Cost ◽  
Nonholonomic Constraints ◽  
State Variables ◽  
Well Conditioning ◽  
Fifth Order

A method is presented for formulating and numerically integrating ordinary differential equations of motion for nonholonomically constrained multibody systems. Tangent space coordinates are defined in configuration and velocity spaces as independent generalized coordinates that serve as state variables in the formulation, yielding ordinary differential equations of motion. Orthogonal-dependent coordinates and velocities are used to enforce constraints at position, velocity, and acceleration levels. Criteria that assure accuracy of constraint satisfaction and well conditioning of the reduced mass matrix in the equations of motion are used as the basis for updating local coordinates on configuration and velocity constraint manifolds, transparent to the user and at minimal computational cost. The formulation is developed for multibody systems with nonlinear holonomic constraints and nonholonomic constraints that are linear in velocity coordinates and nonlinear in configuration coordinates. A computational algorithm for implementing the approach is presented and used in the solution of three examples: one planar and two spatial. Numerical results using a fifth-order Runge–Kutta–Fehlberg explicit integrator verify that accurate results are obtained, satisfying all the three forms of kinematic constraint, to within error tolerances that are embedded in the formulation.

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.

A Discrete Quasi-Coordinate Formulation for the Dynamics of Elastic Bodies

2006 ◽  
Vol 74 (2) ◽  
pp. 231-239 ◽  
Author(s):  
G. M. T. D’Eleuterio ◽  
T. D. Barfoot
Keyword(s):  
Equations Of Motion ◽  
Computational Cost ◽  
Dynamical Equations ◽  
Generalized Coordinates ◽  
Order Form ◽  
Elastic Bodies ◽  
Alternative Approach ◽  
Elastic Systems ◽  
Computationally Intensive ◽  
Derivatives Of

The discretized equations of motion for elastic systems are typically displayed in second-order form. That is, the elastic displacements are represented by a set of discretized (generalized) coordinates, such as those used in a finite-element method, and the elastic rates are simply taken to be the time-derivatives of these displacements. Unfortunately, this approach leads to unpleasant and computationally intensive inertial terms when rigid rotations of a body must be taken into account, as is so often the case in multibody dynamics. An alternative approach, presented here, assumes the elastic rates to be discretized independently of the elastic displacements. The resulting dynamical equations of motion are simplified in form, and the computational cost is correspondingly lessened. However, a slightly more complex kinematical relation between the rate coordinates and the displacement coordinates is required. This tack leads to what may be described as a discrete quasi-coordinate formulation.

A novel PID control parameters tuning approach for robot manipulators mounted on oscillatory bases

Robotica ◽  
2007 ◽  
Vol 25 (4) ◽  
pp. 467-477 ◽  
Author(s):  
J. Lin ◽  
Z.-Z. Huang
Keyword(s):  
Pid Control ◽  
Pid Controller ◽  
Low Cost ◽  
Equations Of Motion ◽  
Grey Theory ◽  
Control Parameters ◽  
State Variables ◽  
Relational Analysis ◽  
Parameters Tuning ◽  
Grey Relational

SUMMARYThis research focuses on the issue of dynamic modeling and controlling a robotic manipulator attached to a compliant base. Such a system is known under the name macro–micro system, characterized by the number of control actuators being less than the number of state variables. The equations of motion for a two-link planar elbow arm mounted on an oscillatory base has been presented in this investigation. In order to study the sensitivity of tuning the PID parameters to achieve the desired performance, the Grey relational analysis has first been proposed. Therefore, the aim of this work is to apply Grey theory to optimize parameters for partial states feedback of a PID controller for such a structure. The experimental results of the proposed methodology also show that it is technically and economically feasible to develop a low-cost, reliable, automatic, less time-consuming controller for robotics mounted on oscillatory bases.

Trajectory Tracking of Underactuated Multibody Systems

2011 ◽  
Author(s):  
Stefan Reichl ◽  
Wolfgang Steiner
Keyword(s):  
Trajectory Tracking ◽  
Multibody Systems ◽  
Degrees Of Freedom ◽  
Inverse Dynamics ◽  
Feedforward Control ◽  
Equations Of Motion ◽  
Three Dimensional ◽  
Differential Algebraic Equations ◽  
State Variables ◽  
Algebraic Equations

This work presents three different approaches in inverse dynamics for the solution of trajectory tracking problems in underactuated multibody systems. Such systems are characterized by less control inputs than degrees of freedom. The first approach uses an extension of the equations of motion by geometric and control constraints. This results in index-five differential-algebraic equations. A projection method is used to reduce the systems index and the resulting equations are solved numerically. The second method is a flatness-based feedforward control design. Input and state variables can be parameterized by the flat outputs and their time derivatives up to a certain order. The third approach uses an optimal control algorithm which is based on the minimization of a cost functional including system outputs and desired trajectory. It has to be distinguished between direct and indirect methods. These specific methods are applied to an underactuated planar crane and a three-dimensional rotary crane.

Multi-Parameter Optimization of Damped Linear Continuous Systems

1972 ◽  
Vol 94 (1) ◽  
pp. 1-7 ◽  
Author(s):  
O. B. Dale ◽  
R. Cohen
Keyword(s):  
Optimal Parameter ◽  
Equations Of Motion ◽  
Analytic Solutions ◽  
System Response ◽  
Continuous Systems ◽  
Frequency Interval ◽  
State Variables ◽  
Initial Value ◽  
Unknown Parameters ◽  
One Dimensional

A method is presented for obtaining and optimizing the frequency response of one-dimensional damped linear continuous systems. The systems considered are assumed to contain unknown constant parameters in the boundary conditions and equations of motion which the designer can vary to obtain a minimum resonant response in some selected frequency interval. The unknown parameters need not be strictly dissipative nor unconstrained. No analytic solutions, either exact or approximate, are required for the system response and only initial value numerical integrations of the state and adjoint differential equations are required to obtain the optimal parameter set. The combinations of state variables comprising the response and the response locations are arbitrary.

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

Toward a Minimal Dynamic Model of a 6-DOF Parallel Robot

1997 ◽  
Author(s):  
L. Beji ◽  
M. Pascal ◽  
P. Joli
Keyword(s):  
Dynamic Model ◽  
Degrees Of Freedom ◽  
Equations Of Motion ◽  
Computational Cost ◽  
Closed Form Solution ◽  
Parallel Robot ◽  
Form Solution ◽  
Lagrangian Formalism ◽  
Inverse Kinematic Problem ◽  
Geometrical Properties

Abstract In this paper, an architecture of a six degrees of freedom (dof) parallel robot and three limbs is described. The robot is called Space Manipulator (SM). In a first step, the inverse kinematic problem for the robot is solved in closed form solution. Further, we need to inverse only a 3 × 3 passive jacobian matrix to solve the direct kinematic problem. In a second step, the dynamic equations are derived by using the Lagrangian formalism where the coordinates are the passive and active joint coordinates. Based on geometrical properties of the robot, the equations of motion are derived in terms of only 9 coordinates related by 3 kinematic constraints. The computational cost of the obtained dynamic model is reduced by using a minimum set of base inertial parameters.

Simulation of Planar Dynamic Mechanical Systems With Changing Topologies: Part I — Characterization and Prediction of the Kinematic Constraint Changes

1987 ◽  
Author(s):  
B. J. Gilmore ◽  
R. J. Cipra
Keyword(s):  
Equations Of Motion ◽  
Mechanical Systems ◽  
Rigid Bodies ◽  
State Variables ◽  
Kinematic Constraint ◽  
Kinematic Constraints ◽  
Force Closure ◽  
Simulation Strategy ◽  
Reaction Forces ◽  
Discontinuous Equations

Abstract Due to changes in the kinematic constraints, many mechanical systems are described by discontinuous equations of motion. This paper addresses those changes in the kinematic constraints which are caused by planar bodies contacting and separating. A strategy to automatically predict and detect the kinematic constraint changes, which are functions of the system dynamics, is presented in Part I. The strategy employs the concepts of point to line contact kinematic constraints, force closure, and ray firing together with the information provided by the rigid bodies’ boundary descriptions, state variables, and reaction forces to characterize the kinematic constraint changes. Since the strategy automatically predicts and detects constraint changes, it is capable of simulating mechanical systems with unpredictable or unforeseen changes in topology. Part II presents the implementation of the characterizations into a simulation strategy and presents examples.

Dynamics Modeling of Elastic Multibody Systems Using Kane’s Method As Applied to a Flexible Robot

1997 ◽  
Author(s):  
Ali Meghdari ◽  
Farbod Fahimi
Keyword(s):  
Multibody Systems ◽  
Equations Of Motion ◽  
Component Mode Synthesis ◽  
Flexible Robot ◽  
Dynamics Modeling ◽  
Generalized Coordinates ◽  
Kane’S Method ◽  
Kane's Method ◽  
Elastic Multibody Systems ◽  
Elastic Form

Abstract Generalization of Kane’s equations of motion for elastic multibody systems is considered. Initially, finite element techniques are used to generate the elastic form of generalized coordinates. Then, the number of elastic coordinates are reduced by the component mode synthesis. Finally, Kane’s method is applied to obtain the equations of motion of such systems. Using this method, dynamic model of an elastic robot with one degree of freedom is presented.

Sign in / Sign up

Export Citation Format

Share Document