A Unified Framework for Dynamics and Lyapunov Stability of Holonomically Constrained Rigid Bodies

Author(s):  
Khoder Melhem ◽  
◽  
Zhaoheng Liu ◽  
Antonio Loría ◽  
◽  
...  

A new dynamic model for interconnected rigid bodies is proposed here. The model formulation makes it possible to treat any physical system with finite number of degrees of freedom in a unified framework. This new model is a nonminimal realization of the system dynamics since it contains more state variables than is needed. A useful discussion shows how the dimension of the state of this model can be reduced by eliminating the redundancy in the equations of motion, thus obtaining the minimal realization of the system dynamics. With this formulation, we can for the first time explicitly determine the equations of the constraints between the elements of the mechanical system corresponding to the interconnected rigid bodies in question. One of the advantages coming with this model is that we can use it to demonstrate that Lyapunov stability and control structure for the constrained system can be deducted by projection in the submanifold of movement from appropriate Lyapunov stability and stabilizing control of the corresponding unconstrained system. This procedure is tested by some simulations using the model of two-link planar robot.

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.


2019 ◽  
Vol 2019 (2) ◽  
pp. 59-75
Author(s):  
Jarosław Stanisławski

Summary The paper presents simulation method and results of calculations determining behavior of helicopter and landing site loads which are generated during phase of the helicopter take-off and landing. For helicopter with whirling rotor standing on ground or touching it, the loads of landing gear depend on the parameters of helicopter movement, occurrence of wind gusts and control of pitch angle of the rotor blades. The considered model of helicopter consists of the fuselage and main transmission treated as rigid bodies connected with elastic elements. The fuselage is supported by landing gear modeled by units of spring and damping elements. The rotor blades are modeled as elastic axes with sets of lumped masses of blade segments distributed along them. The Runge-Kutta method was used to solve the equations of motion of the helicopter model. According to the Galerkin method, it was assumed that the parameters of the elastic blade motion can be treated as a combination of its bending and torsion eigen modes. For calculations, data of a hypothetical light helicopter were applied. Simulation results were presented for the cases of landing helicopter touching ground with different vertical speed and for phase of take-off including influence of rotor speed changes, wind gust and control of blade pitch. The simulation method may help to define the limits of helicopter safe operation on the landing surfaces.


SIMULATION ◽  
1964 ◽  
Vol 2 (2) ◽  
pp. R-9-R-23
Author(s):  
Edward E. Markson ◽  
John L. Stricker

Space mission simulator programs may be divided into two broad categories: (1) training tools (quali tative devices often simulating a continuous mission), and (2) laboratory tools (quantitative devices treating the mission in phases, each phase being programmed separately to obtain optimum scaling). This paper describes the development of an analog program capable of continuously simulating an entire lunar mission in six degrees of freedom with high resolu tion throughout. The reported work logically traces the program development through the equations of motion, the guidance and control equations, and the analog mechanization. The translation equations are de veloped using a modified form of Encke's method; two reference origins are utilized at the two points of primary interest—the landing site and the target vehicle—such that the displacements are approach ing a minimum in the regions where the highest reso lution is required. The variables are rescaled as this region is approached to obtain maximum accuracy. Relays, stepping switches and diode gates are used for rescaling and to re-reference origins. A particular Euler angle sequence is selected based on matrix validity criteria applied to the mission. A previously reported guidance technique is shown to be appli cable to all phases of the mission. It is concluded that the method demonstrated in this paper leads to minimum computer loading for simulating a manned space mission without program discontinuities. Supporting data include an analog- computed trajectory representative of a long-dura tion mission, which is compared in detail with a digital solution.


Author(s):  
B. J. Gilmore ◽  
R. J. Cipra

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.


2020 ◽  
Vol 39 (10-11) ◽  
pp. 1239-1258
Author(s):  
Shameek Ganguly ◽  
Oussama Khatib

Multi-surface interactions occur frequently in articulated-rigid-body systems such as robotic manipulators. Real-time prediction of contact-interaction forces is challenging for systems with many degrees of freedom (DOFs) because joint and contact constraints must be enforced simultaneously. While several contact models exist for systems of free rigid bodies, fewer models are available for articulated-body systems. In this paper, we extend the method of Ruspini and Khatib and develop the contact-space resolution (CSR) model by applying the operational space theory of robot manipulation. Through a proper choice of contact-space coordinates, the projected dynamics of the system in the contact space is obtained. We show that the projection into the dynamically consistent null space preserves linear and angular momentum in a subspace of the system dynamics complementary to the joint and contact constraints. Furthermore, we illustrate that a simultaneous collision event between two articulated bodies can be resolved as an equivalent simultaneous collision between two non-articulated rigid bodies through the projected contact-space dynamics. Solving this reduced-dimensional problem is computationally efficient, but determining its accuracy requires physical experimentation. To gain further insights into the theoretical model predictions, we devised an apparatus consisting of colliding 1-, 2-, and 3-DOF articulated bodies where joint motion is recorded with high precision. Results validate that the CSR model accurately predicts the post-collision system state. Moreover, for the first time, we show that the projection of system dynamics into the mutually complementary contact space and null space is a physically verifiable phenomenon in articulated-rigid-body systems.


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.


Author(s):  
Tamer Wasfy

A new technique for modeling rigid bodies undergoing spatial motion using an explicit time-integration finite element code is presented. The key elements of the technique are: (a) use of the total rotation matrix relative to the inertial frame to measure the rotation of the rigid bodies; (b) time-integration of the rotational equations of motion in a body fixed (material) frame, with the resulting incremental rotations added to the total rotation matrix; (c) penalty formulation for creating connection points (virtual nodes which do not add extra degrees of freedom) on the rigid-body where joints can be placed. The use of the rotation matrix along with incremental rotation updates circumvents the problem of singularities associated with other types of three and four parameter rotation measures. Benchmark rigid multibody dynamics problems are solved to demonstrate the accuracy of the present technique.


Author(s):  
S. M. Mehdi Ansarey M. ◽  
M. J. Mahjoob

In this paper, the dynamics and control of an automated guided vehicle (AGV) is described. The objective is to control the vehicle direction and location with respect to a prescribed trajectory. This is accomplished based on an optimum control strategy using vehicle state variables. A four-wheel vehicle with three degrees of freedom including longitudinal, lateral and yaw motion is considered. The nonlinearity of the tire and steering mechanism is also included. The control system design for circular, straight forward and composite path is presented based on feedback linearization. Some trajectory simulation for discrete curvatures is carried out. The controller was implemented within MATLAB environment. The design was also evaluated using ADAMS full vehicle assembly. The results demonstrated the accuracy of the model and the effectiveness of the developed control system.


Author(s):  
Ganesh Venkatesan ◽  
William Clark

The application of computational fluid dynamics method to the submarine maneuvering simulations of ONR Body 1 is presented. ONR Body 1 is an unclassified submarine radio controlled model with propeller and control surfaces. Unsteady Reynolds-averaged Naviers-Stokes equations of fluid flow is coupled to the six degrees-of-freedom equations of motion of a rigid body via user coding to predict the instantaneous position and body orientation. Propeller and control surface motions are accounted for by using the moving mesh feature integrated into the solution procedure which allows sliding interfaces between different mesh blocks of the computational domain (for propeller rotation), as well as mesh distortion (for control surface deflection). This offers the flexibility of using a single computational grid for the entire simulation period. The maneuvers simulated include a constant depth and heading run as well as a horizontal overshoot maneuver using conditions consistent with the experiment. Predicted results show favorable agreement with experimental measurements.


2005 ◽  
Vol 2005 (4) ◽  
pp. 365-382 ◽  
Author(s):  
Hazem Ali Attia

A dynamic model for multi-rigid-body systems which consists of interconnected rigid bodies based on particle dynamics and a recursive approach is presented. The method uses the concepts of linear and angular momentums to generate the rigid body equations of motion in terms of the Cartesian coordinates of a dynamically equivalent constrained system of particles, without introducing any rotational coordinates and the corresponding rotational transformation matrix. For the open-chain system, the equations of motion are generated recursively along the serial chains. A closed-chain system is transformed to open-chain by cutting suitable kinematical joints and introducing cut-joint constraints. An example is chosen to demonstrate the generality and simplicity of the developed formulation.


Sign in / Sign up

Export Citation Format

Share Document