A Fast Algorithm for Contact Dynamics of Multibody Systems Using the Box Friction Model

Author(s):  
Farnood Gholami ◽  
Mostafa Nasri ◽  
József Kövecses ◽  
Marek Teichmann

One of the major challenges in dynamics of multibody systems is to handle redundant constraints appropriately. The box friction model is one of the existing approaches to formulate the contact and friction phenomenon as a mixed linear complementarity problem (MLCP). In this setting, the contact redundancy can be handled by relaxing the constraints, but such a technique might suffer from certain drawbacks, specially in the case of large number of redundant constraints. Most of the common pivoting algorithms used to solve the resulting mixed complementarity problem might not converge when the relaxation terms are chosen as small as they should be. To overcome the aforementioned shortcoming, we propose a novel approach which takes advantage of the sparse structure of the formulated MLCP. This novel approach reduces the sensitivity of the solution of the problem to the relaxation terms and decreases the number of required pivots to obtain the solution, leading to shorter computational times. Furthermore, as a result of the proposed approach, much smaller relaxation terms can be used while the solution algorithms converge.

Author(s):  
Albert Peiret ◽  
József Kövecses ◽  
Josep M. Font-Llagunes

The dynamics of multibody systems with many contacts are frequently formulated as a Linear Complementarity Problem (LCP), for which several direct or iterative algorithms are available to solve it efficiently. These formulations rely on discretized friction models that approximate the friction cone of the Coulomb model to a pyramid. However, they produce rank-deficient LCPs even though the physical problem does not have constraint redundancy and has a unique solution. Here, a new discretized friction model is presented which results in an LCP formulation with a full-rank lead matrix. This model relies on an inertial term to couple the equations of the model, which behaves as close to the Coulomb model as the other discretized models. Moreover, it is shown through some simulations that some algorithms can be used with this formulation, which could not be used with the other rank-deficient LCP formulations.


Author(s):  
P. Flores ◽  
J. Ambro´sio ◽  
J. C. P. Claro ◽  
H. M. Lankarani

This work deals with a methodology to assess the influence of the spherical clearance joints in spatial multibody systems. The methodology is based on the Cartesian coordinates, being the dynamics of the joint elements modeled as impacting bodies and controlled by contact forces. The impacts and contacts are described by a continuous contact force model that accounts for geometric and mechanical characteristics of the contacting surfaces. The contact force is evaluated as function of the elastic pseudo-penetration between the impacting bodies, coupled with a nonlinear viscous-elastic factor representing the energy dissipation during the impact process. A spatial four bar mechanism is used as an illustrative example and some numerical results are presented, being the efficiency of the developed methodology discussed in the process of their presentation. The results obtained show that the inclusion of clearance joints in the modelization of spatial multibody systems significantly influences the prediction of components’ position and drastically increases the peaks in acceleration and reaction moments at the joints. Moreover, the system’s response clearly tends to be nonperiodic when a clearance joint is included in the simulation.


2021 ◽  
Author(s):  
Adwait Verulkar ◽  
Corina Sandu ◽  
Daniel Dopico ◽  
Adrian Sandu

Abstract Sensitivity analysis is one of the most prominent gradient based optimization techniques for mechanical systems. Model sensitivities are the derivatives of the generalized coordinates defining the motion of the system in time with respect to the system design parameters. These sensitivities can be calculated using finite differences, but the accuracy and computational inefficiency of this method limits its use. Hence, the methodologies of direct and adjoint sensitivity analysis have gained prominence. Recent research has presented computationally efficient methodologies for both direct and adjoint sensitivity analysis of complex multibody dynamic systems. The contribution of this article is in the development of the mathematical framework for conducting the direct sensitivity analysis of multibody dynamic systems with joint friction using the index-1 formulation. For modeling friction in multibody systems, the Brown and McPhee friction model has been used. This model incorporates the effects of both static and dynamic friction on the model dynamics. A case study has been conducted on a spatial slider-crank mechanism to illustrate the application of this methodology to real-world systems. Using computer models, with and without joint friction, effect of friction on the dynamics and model sensitivities has been demonstrated. The sensitivities of slider velocity have been computed with respect to the design parameters of crank length, rod length, and the parameters defining the friction model. Due to the highly non-linear nature of friction, the model dynamics are more sensitive during the transition phases, where the friction coefficient changes from static to dynamic and vice versa.


2018 ◽  
Vol 2018 ◽  
pp. 1-13
Author(s):  
Leiming Ning ◽  
Jichang Chen ◽  
Mingbo Tong

A high-fidelity cargo airdrop simulation requires the accurate modeling of the contact dynamics between an aircraft and its cargo. This paper presents a general and efficient contact-friction model for the simulation of aircraft-cargo coupling dynamics during an airdrop extraction phase. The proposed approach has the same essence as the finite element node-to-segment contact formulation, which leads to a flexible, straightforward, and efficient code implementation. The formulation is developed under an arbitrary moving frame with both aircraft and cargo treated as general six degrees-of-freedom rigid bodies, thus eliminating the restrictions of lateral symmetric assumptions in most existing methods. Moreover, the aircraft-cargo coupling algorithm is discussed in detail, and some practical implementation details are presented. The accuracy and capability of the present method are demonstrated through four numerical examples with increasing complexity and fidelity.


Sign in / Sign up

Export Citation Format

Share Document