A Hybrid Highly Parallelizable Algorithm for the Dynamics of Rigid Body Chain Systems

Author(s):  
Shanzhong Duan ◽  
Kurt S. Anderson

Abstract The paper presents a new hybrid parallelizable low order algorithm for modeling the dynamic behavior of multi-rigid-body chain systems. The method is based on cutting certain system interbody joints so that largely independent multibody subchain systems are formed. These subchains interact with one another through associated unknown constraint forces f¯c at the cut joints. The increased parallelism is obtainable through cutting the joints and the explicit determination of associated constraint loads combined with a sequential O(n) procedure. In other words, sequential O(n) procedures are performed to form and solve equations of motion within subchains and parallel strategies are used to form and solve constraint equations between subchains in parallel. The algorithm can easily accommodate the available number of processors while maintaining high efficiency. An O[(n+m)Np+m(1+γ)Np+mγlog2Np](0<γ<1) performance will be achieved with Np processors for a chain system with n degrees of freedom and m constraints due to cutting of interbody joints.

Author(s):  
Shanzhong (Shawn) Duan ◽  
Yogesh Patel

This paper presents a hybrid parallelizable algorithm for the computer-aided modeling of the dynamic behavior of multibody open tree systems. The method is based on cutting certain system interbody joints at branched bodies so that a system of largely independent multibody subchains is formed. These subchains interact with one another through associated unknown constraint forces fc at the cut joints. The increased parallelism is obtainable through cutting joints and the explicit determination of associated constraint forces combined with a sequential O(n) procedure. Consequently, the sequential O(n) procedure is carried out within each subchain to form and solve the equations of motion, while parallel strategies are performed between the subchains to form and solve constraint equations concurrently. Two extreme cutting procedures are further discussed to conduct a comparison of computational efficiency. One case is to cut the joints at branched bodies so that the longest lengths of subchains are obtained and the other is to cut the joints at branched bodies so that the shortest lengths of subchains are formed. The algorithm can easily accommodate the available number of processors while maintaining high efficiency. The algorithm will also be implemented on both parallel homogeneous computing (PHC) systems and network-distributed heterogeneous computing (DHC) environment. The implementation of the algorithm in a DHC environment will permit engineers and researchers to conduct distributed simulation of dynamic behaviors of large and complex multibody systems on ubiquitous network-distributed PC workstations in their workplace. In short, the exploration of the parallelizable algorithm for a multibody tree system will provide a deep understanding of the relationship between computational load balancing and optimal joint locations to be cut. The computational efficiency of the algorithm can be increased further.


2019 ◽  
Vol 24 (2) ◽  
pp. 175-180
Author(s):  
Vladimir Dragoş Tătaru ◽  
Mircea Bogdan Tătaru

Abstract The present paper approaches in an original manner the dynamic analysis of a wheel which climbs on an inclined plane under the action of a horizontal force. The wheel rolls and slides in the same time. The two movements, rolling and sliding are considered to be independent of each other. Therefore we are dealing with a solid rigid body with two degrees of freedom. The difficulty of approaching the problem lies in the fact that in the differential equations describing the motion of the solid rigid body are also present the constraint forces and these are unknown. For this reason they must be eliminated from the differential equations of motion. The paper presents as well an original method of the constraint forces elimination.


Author(s):  
Shanzhong Duan

Molecular dynamics is effective for a nano-scale phenomenon analysis. This paper presents a hybrid parallelizable algorithm for the computer simulation of the motion behavior of molecular chain and open-tree structure on parallel computing system. The algorithm is developed from an approach of rigid body dynamics, in which interbody constraints are exposed so that a system of largely independent multibody subchains is formed. The increased parallelism is obtainable through bringing interbody constraints to evidence and the explicit determination of the associated constraint forces combined with a sequential O(n) procedure. Each subchain then is assigned to a processor for parallel computing. The algorithm offers a sequential O(n) performance if there is only one processor available. The algorithm has O(log2n) computational efficiency if there are as many processors available as number for molecular bodies. For most common scenario, the algorithm will give a computational complexity between O(n) and O(log2n) if number of available processor is less than number of molecular bodies.


Author(s):  
Andreas Müller ◽  
Shivesh Kumar

AbstractDerivatives of equations of motion (EOM) describing the dynamics of rigid body systems are becoming increasingly relevant for the robotics community and find many applications in design and control of robotic systems. Controlling robots, and multibody systems comprising elastic components in particular, not only requires smooth trajectories but also the time derivatives of the control forces/torques, hence of the EOM. This paper presents the time derivatives of the EOM in closed form up to second-order as an alternative formulation to the existing recursive algorithms for this purpose, which provides a direct insight into the structure of the derivatives. The Lie group formulation for rigid body systems is used giving rise to very compact and easily parameterized equations.


Joint Rail ◽  
2004 ◽  
Author(s):  
Mohammad Durali ◽  
Mohammad Mehdi Jalili Bahabadi

In this article a train model is developed for studying train derailment in passing through bends. The model is three dimensional, nonlinear, and considers 43 degrees of freedom for each wagon. All nonlinear characteristics of suspension elements as well as flexibilities of wagon body and bogie frame, and the effect of coupler forces are included in the model. The equations of motion for the train are solved numerically for different train conditions. A neural network was constructed as an element in solution loop for determination of wheel-rail contact geometry. Derailment factor was calculated for each case. The results are presented and show the major role of coupler forces on possible train derailment.


2017 ◽  
Vol 2017 ◽  
pp. 1-13 ◽  
Author(s):  
T. S. Amer

In this paper, we will focus on the dynamical behavior of a rigid body suspended on an elastic spring as a pendulum model with three degrees of freedom. It is assumed that the body moves in a rotating vertical plane uniformly with an arbitrary angular velocity. The relative periodic motions of this model are considered. The governing equations of motion are obtained using Lagrange’s equations and represent a nonlinear system of second-order differential equations that can be solved in terms of generalized coordinates. The numerical solutions are investigated using the fourth-order Runge-Kutta algorithms through Matlab packages. These solutions are represented graphically in order to describe and discuss the behavior of the body at any instant for different values of the physical parameters of the body. The obtained results have been discussed and compared with some previous published works. Some concluding remarks have been presented at the end of this work. The importance of this work is due to its numerous applications in life such as the vibrations that occur in buildings and structures.


Author(s):  
Michael J. Sadowski ◽  
Kurt S. Anderson

This paper presents an algorithm for the efficient numerical analysis and simulation of a category of contact/impact problems in multi-rigid-body dynamic systems with tree topologies. The algorithm can accommodate the jumps in structure which occur in the equations of motion of general multi-rigid-body systems due to a contact/impact event between bodies, or due to the locking of joints as long as the resulting system is a tree topology. The presented method uses a generalized momentum balance approach to determine the velocity jumps which take place across impacts in such multibody dynamic systems where event constraint forces are of the “non-working” category. The presented method does not suffer from the performance (speed) penalty encountered by most other momentum balance methods given its O(n) overall cost, and exact direct embedded consideration of all the constraints. Due to these characteristics, the presented algorithm offers superior computing performance relative to other methods in situations involving both large n and potentially many unilateral constraints.


Author(s):  
Pratheek Bagivalu Prasanna ◽  
Ashok Midha ◽  
Sushrut G. Bapat

Abstract Understanding the kinematic properties of a compliant mechanism has always proved to be a challenge. A concept of compliance number offered earlier emphasized the development of terminology that aided in its determination. A method to evaluate the elastic degrees of freedom associated with the flexible segments/links of a compliant mechanism using the pseudo-rigid-body model (PRBM) concept is provided. In this process, two distinct classes of compliant mechanisms are developed involving: (i) Active Compliance and (ii) Passive Compliance. Furthermore, these also aid in a better characterization of the kinematic behavior of a compliant mechanism. A more lucid interpretation of the significance of compliance number is provided. Applications of this method to both active and passive compliant mechanisms are exemplified. Finally, an experimental procedure that aids in visualizing the degrees of freedom as calculated is presented.


1971 ◽  
Vol 93 (1) ◽  
pp. 191-195 ◽  
Author(s):  
Desideriu Maros ◽  
Nicolae Orlandea

This paper is a further development of the kinematic problem presented in our 1967 paper [1] in which we have obtained the transmission functions for different orders of plane systems with many degrees of freedom. This paper establishes the corresponding system of differential equations of motion beginning with these functions. The purpose of this paper is to facilitate computer programming. Our study is based on the work of R. Beyer [2, 3] and is the first original addition to his papers. A second original contribution to Beyer’s theories is the deductive method of solution, from general to particular, which we have, incorporated in our work. Beyer concluded that the cases having two or three degrees of freedom can be considered as particular solutions to the results obtained.


Author(s):  
Shanzhong Duan ◽  
Andrew Ries

This paper presents an efficient parallelizable algorithm for the computer-aided simulation and numerical analysis of motion behaviors of multibody systems with closed-loops. The method is based on cutting certain user-defined system interbody joints so that a system of independent multibody subchains is formed. These subchains interact with one another through associated unknown constraint forces fc at the cut joints. The increased parallelism is obtainable through cutting joints and the explicit determination of associated constraint forces combined with a sequential O(n) method. Consequently, the sequential O(n) procedure is carried out within each subchain to form and solve the equations of motion while parallel strategies are performed between the subchains to form and solve constraint equations concurrently. For multibody systems with closed-loops, joint separations play both a role of creation of parallelism for computing load distribution and a role of opening a closed-loop for use of the O(n) algorithm. Joint separation strategies provide the flexibility for use of the algorithm so that it can easily accommodate the available number of processors while maintaining high efficiency. The algorithm gives the best performance for the application scenarios for n>>1 and n>>m, where n and m are number of degree of freedom and number of constraints of a multibody system with closed-loops respectively. The algorithm can be applied to both distributed-memory parallel computing systems and shared-memory parallel computing systems.


Sign in / Sign up

Export Citation Format

Share Document