Automating the Derivation of the Equations of Motion of a Multibody Dynamic System With Uncertainty Using Polynomial Chaos Theory and Variational Work

2019 ◽  
Vol 15 (1) ◽  
Author(s):  
Paul S. Ryan ◽  
Sarah C. Baxter ◽  
Philip A. Voglewede

Abstract Understanding how variation impacts a multibody dynamic (MBD) system's response is important to ensure the robustness of a system. However, how the variation propagates into the MBD system is complicated because MBD systems are typically governed by a system of large differential algebraic equations. This paper presents a novel process, variational work, along with the polynomial chaos multibody dynamics (PCMBoD) automation process for utilizing polynomial chaos theory (PCT) in the analysis of uncertainties in an MBD system. Variational work allows the complexity of the traditional PCT approach to be reduced. With variational work and the constrained Lagrangian formulation, the equations of motion of an MBD PCT system can be constructed using the PCMBoD automated process. To demonstrate the PCMBoD process, two examples, a mass-spring-damper and a two link slider–crank mechanism, are shown.

Author(s):  
Paul S. Ryan ◽  
Sarah C. Baxter ◽  
Philip A. Voglewede

Variation occurs in many closed loop multi-body dynamic (MBD) systems in the geometry, mass, or forces. Understanding how MBD systems respond to variation is imperative for the design of a robust system. However, simulation of how variation propagates into the solution is complicated as most MBD systems cannot be simplified into to a system of ordinary differential equations (ODE). This paper investigates polynomial chaos theory (PCT) as a means of quantifying the effects of uncertainty in a closed loop MBD system governed by differential algebraic equations (DAE). To demonstrate how PCT could be used, the motion of a two link slider-crank mechanism is simulated with variation in the link lengths. To validate and show the advantages and disadvantages of PCT in closed loop MBD systems, the PCT approach is compared to Monte Carlo simulations.


Author(s):  
William Prescott

This paper will examine the importance of applying scaling to the equations of motion for multibody dynamic systems when applied to industrial applications. If a Cartesian formulation is used to formulate the equations of motion of a multibody dynamic system the resulting equations are a set of differential algebraic equations (DAEs). The algebraic components of the DAEs arise from appending the joint equations used to model revolute, cylindrical, translational and other joints to the Newton-Euler dynamic equations of motion. Stability issues can arise in an ill-conditioned Jacobian matrix of the integration method this will result in poor convergence of the implicit integrator’s Newton method. The repeated failures of the Newton’s method will require a small step size and therefore simulations that require long run times to complete. Recent advances in rescaling the equations of motion have been proposed to address this problem. This paper will see if these methods or a variant addresses not only stability concerns, but also efficiency. The scaling techniques are applied to the Gear-Gupta-Leimkuhler (GGL) formulation for multibody problems by embedding them into the commercial multibody code (MBS) Virtual. Lab Motion and then use them to solve an industrial sized automotive example to see if performance is improved.


Author(s):  
Edward J. Haug ◽  
Mirela Iancu ◽  
Dan Negrut

Abstract An implicit numerical integration approach, based on generalized coordinate partitioning of the descriptor form of the differential-algebraic equations of motion of multibody dynamics, is presented. This approach is illustrated for simulation of stiff mechanical systems using the well known Newmark integration method from structural dynamics. Second order Newmark integration formulas are used to define independent generalized coordinates and their first time derivative as functions of independent accelerations. The latter are determined as the solution of discretized equations obtained using the descriptor form of the equations of motion. Dependent variables in the formulation, including Lagrange multipliers, are determined to satisfy all the kinematic and kinetic equations of multibody dynamics. The approach is illustrated by solving the constrained equations of motion for mechanical systems that exhibit stiff behavior. Results show that the approach is robust and has the capability to integrate differential-algebraic equations of motion for stiff multibody dynamic systems.


Author(s):  
J. P. Meijaard ◽  
V. van der Wijk

Some thoughts about different ways of formulating the equations of motion of a four-bar mechanism are communicated. Four analytic methods to derive the equations of motion are compared. In the first method, Lagrange’s equations in the traditional form are used, and in a second method, the principle of virtual work is used, which leads to equivalent equations. In the third method, the loop is opened, principal points and a principal vector linkage are introduced, and the equations are formulated in terms of these principal vectors, which leads, with the introduced reaction forces, to a system of differential-algebraic equations. In the fourth method, equivalent masses are introduced, which leads to a simpler system of principal points and principal vectors. By considering the links as pseudorigid bodies that can have a uniform planar dilatation, a compact form of the equations of motion is obtained. The conditions for dynamic force balance become almost trivial. Also the equations for the resulting reaction moment are considered for all four methods.


Author(s):  
Francisco González ◽  
Pierangelo Masarati ◽  
Javier Cuadrado ◽  
Miguel A. Naya

Formulating the dynamics equations of a mechanical system following a multibody dynamics approach often leads to a set of highly nonlinear differential-algebraic equations (DAEs). While this form of the equations of motion is suitable for a wide range of practical applications, in some cases it is necessary to have access to the linearized system dynamics. This is the case when stability and modal analyses are to be carried out; the definition of plant and system models for certain control algorithms and state estimators also requires a linear expression of the dynamics. A number of methods for the linearization of multibody dynamics can be found in the literature. They differ in both the approach that they follow to handle the equations of motion and the way in which they deliver their results, which in turn are determined by the selection of the generalized coordinates used to describe the mechanical system. This selection is closely related to the way in which the kinematic constraints of the system are treated. Three major approaches can be distinguished and used to categorize most of the linearization methods published so far. In this work, we demonstrate the properties of each approach in the linearization of systems in static equilibrium, illustrating them with the study of two representative examples.


Author(s):  
Radu Serban ◽  
Jeffrey S. Freeman

Abstract Methods for formulating the first-order design sensitivity of multibody systems by direct differentiation are presented. These types of systems, when formulated by Euler-Lagrange techniques, are representable using differential-algebraic equations (DAE). The sensitivity analysis methods presented also result in systems of DAE’s which can be solved using standard techniques. Problems with previous direct differentiation sensitivity analysis derivations are highlighted, since they do not result in valid systems of DAE’s. This is shown using the simple pendulum example, which can be analyzed in both ODE and DAE form. Finally, a slider-crank example is used to show application of the method to mechanism analysis.


Author(s):  
Stefan Reichl ◽  
Wolfgang Steiner

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.


2003 ◽  
Vol 25 (3) ◽  
pp. 170-185
Author(s):  
Dinh Van Phong

The article deals with the problem of consistent initial values of the system of equations of motion which has the form of the system of differential-algebraic equations. Direct treating the equations of mechanical systems with particular properties enables to study the system of DAE in a more flexible approach. Algorithms and examples are shown in order to illustrate the considered technique.


Sign in / Sign up

Export Citation Format

Share Document