Generalized Vector-Network Formulation for the Dynamic Simulation of Multibody Systems

1986 ◽  
Vol 108 (4) ◽  
pp. 322-329 ◽  
Author(s):  
M. J. Richard ◽  
R. Anderson ◽  
G. C. Andrews

This paper describes the vector-network approach which is a comprehensive mathematical model for the systematic formulation of the nonlinear equations of motion of dynamic three-dimensional constrained multi-body systems. The entire procedure is a basic application of concepts of graph theory in which laws of vector dynamics have been combined. The main concepts of the method have been explained in previous publications but the work described herein is an appreciable extension of this relatively new approach. The method casts simultaneously the three-dimensional inertia equations associated with each rigid body and the geometrical expressions corresponding to the kinematic restrictions into a symmetrical format yielding the differential equations governing the motion of the system. The algorithm is eminently well suited for the computer-aided simulation of arbitrary interconnected rigid bodies; it serves as the basis for a “self-formulating” computer program which can simulate the response of a dynamic system, given only the system description.

2011 ◽  
Vol 27 (4) ◽  
pp. 479-491 ◽  
Author(s):  
W. Jiang ◽  
X. D. Chen ◽  
X. Luo ◽  
Y. T. Hu ◽  
H. P. Hu

ABSTRACTMany kinds of mechanical systems can be modeled as spatial rigid multibody systems (SR-MBS), which consist of a set of rigid bodies interconnected by joints, springs and dampers. Vibration calculation of SR-MBS is conventionally conducted by approximately linearizing the nonlinear equations of motion and constraint, which is very complicated and inconvenient for sensitivity analysis. A new algorithm based on constraint-topology transformation is presented to derive the oscillatory differential equations in three steps, that is, vibration equations for free SR-MBS are derived using Lagrangian method at first; then, an open-loop constraint matrix is derived to obtain the vibration equations for open-loop SR-MBS via quadric transformation; finally, a cut-joint constraint matrix is derived to obtain the vibration equations for closed-loop SR-MBS via quadric transformation. Through mentioned above, the vibration calculation can be significantly simplified and the sensitivity analysis can be conducted conveniently. The correctness of the proposed method has been verified by numerical experiments in comparison with the traditional approaches.


Author(s):  
Quan Gu ◽  
Jinghao Pan ◽  
Yongdou Liu

Consistent tangent stiffness plays a crucial role in delivering a quadratic rate of convergence when using Newton’s method in solving nonlinear equations of motion. In this paper, consistent tangent stiffness is derived for a three-dimensional (3D) wheel–rail interaction element (WRI element for short) originally developed by the authors and co-workers. The algorithm has been implemented in finite element (FE) software framework (OpenSees in this paper) and proven to be effective. Application examples of wheelset and light rail vehicle are provided to validate the consistent tangent stiffness. The quadratic convergence rate is verified. The speeds of calculation are compared between the use of consistent tangent stiffness and the tangent by perturbation method. The results demonstrate the improved computational efficiency of WRI element when consistent tangent stiffness is used.


Author(s):  
Cagkan Yildiz ◽  
Tamer M. Wasfy ◽  
Hatem M. Wasfy ◽  
Jeanne M. Peters

In order to accurately predict the fatigue life and wear life of a belt, the various stresses that the belt is subjected to and the belt slip over the pulleys must be accurately calculated. In this paper, the effect of material and geometric parameters on the steady-state stresses (including normal, tangential and axial stresses), average belt slip for a flat belt, and belt-drive energy efficiency is studied using a high-fidelity flexible multibody dynamics model of the belt-drive. The belt’s rubber matrix is modeled using three-dimensional brick elements and the belt’s reinforcements are modeled using one dimensional truss elements. Friction between the belt and the pulleys is modeled using an asperity-based Coulomb friction model. The pulleys are modeled as cylindrical rigid bodies. The equations of motion are integrated using a time-accurate explicit solution procedure. The material parameters studied are the belt-pulley friction coefficient and the belt axial stiffness and damping. The geometric parameters studied are the belt thickness and the pulleys’ centers distance.


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.


Author(s):  
Tamer M. Wasfy ◽  
Hatem M. Wasfy

Abstract Belt-drives are used to transmit power between rotational machine elements in many mechanical systems such as industrial machines, home appliances, and internal combustion engines. The belt cross-section typically consists of axially stiff tension cords (made of steel or polyester strands) embedded in a rubber matrix. The rubber matrix provides the friction interface between the belt and the pulleys through which mechanical torque is transmitted. In this paper, the effect of the rubber’s Young’s modulus and Poisson’s ratio on the steady-state belt normal, tangential and axial stresses, average belt slip, and belt-drive energy efficiency is studied using a high-fidelity flexible multibody dynamics model of a flat belt-drive. The belt’s rubber matrix is modeled using three-dimensional brick elements and the belt’s cords are modeled using one dimensional truss elements. Friction between the belt and the pulleys is modeled using an asperity-based Coulomb friction model. The pulleys are modeled as rigid bodies with a cylindrical contact surface. The equations of motion are integrated using a time-accurate explicit solution procedure.


Author(s):  
Thorstein R. Rykkje ◽  
Eystein Gulbrandsen ◽  
Andreas Fosså Hettervik ◽  
Morten Kvalvik ◽  
Daniel Gangstad ◽  
...  

Abstract This paper extends research into flexible robotics through a collaborative, interdisciplinary senior design project. This paper deploys the Moving Frame Method (MFM) to analyze the motion of a relatively high multi-link system, driven by internal servo engines. The MFM describes the dynamics of the system and enables the construction of a general algorithm for the equations of motion. Lie group theory and Cartan’s moving frames are the foundation of this new approach to engineering dynamics. This, together with a restriction on the variation of the angular velocity used in Hamilton’s principle, enables an effective way of extracting the equations of motion. The result is a dynamic 3D analytical model for the motion of a snake-like robotic system, that can take the physical sizes of the system and return the dynamic behavior. Furthermore, this project builds a snake-like robot driven by internal servo engines. The multi-linked robot will have a servo in each joint, enabling a three-dimensional movement. Finally, a test is performed to compare if the theory and the measurable real-time results match.


Author(s):  
Matteo Filippi ◽  
Alfonso Pagani ◽  
Erasmo Carrera

This work explores the effects of geometrical nonlinearities in the vibration analysis of rotating structures and helicopter blades. Structures are modelled via higher-order beam theories with variable kinematics. These theories fall in the domain of the Carrera Unified Formulation (CUF), according to which the nonlinear equations of motion of rotating blades can be written in terms of fundamental nuclei, whose formalism is an invariant of the theory approximation. The inherent three-dimensional nature of CUF enables one to include all Green-Lagrange strain components as well as all coupling effects due to the geometrical features and the three-dimensional constitutive law. Numerical solutions are considered and opportunely discussed. Also, linearized and full nonlinear solutions for vibrating rotating blades are compared both in case of small amplitudes and in the large deflections/rotations regime.


1999 ◽  
Vol 66 (4) ◽  
pp. 986-996 ◽  
Author(s):  
S. K. Saha

Constrained dynamic equations of motion of serial multibody systems consisting of rigid bodies in a serial kinematic chain are derived in this paper. First, the Newton-Euler equations of motion of the decoupled rigid bodies of the system at hand are written. Then, with the aid of the decoupled natural orthogonal complement (DeNOC) matrices associated with the velocity constraints of the connecting bodies, the Euler-Lagrange independent equations of motion are derived. The De NOC is essentially the decoupled form of the natural orthogonal complement (NOC) matrix, introduced elsewhere. Whereas the use of the latter provides recursive order n—n being the degrees-of-freedom of the system at hand—inverse dynamics and order n3 forward dynamics algorithms, respectively, the former leads to recursive order n algorithms for both the cases. The order n algorithms are desirable not only for their computational efficiency but also for their numerical stability, particularly, in forward dynamics and simulation, where the system’s accelerations are solved from the dynamic equations of motion and subsequently integrated numerically. The algorithms are illustrated with a three-link three-degrees-of-freedom planar manipulator and a six-degrees-of-freedom Stanford arm.


Sign in / Sign up

Export Citation Format

Share Document