Full Content is available to subscribers

Subscribe/Learn More  >

Efficient and Robust Implementation of the TLISMNI Method

[+] Author Affiliations
Ahmed Aboubakr, Ahmed A. Shabana

University of Illinois at Chicago, Chicago, IL

Paper No. DETC2015-48105, pp. V006T10A039; 16 pages
  • ASME 2015 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference
  • Volume 6: 11th International Conference on Multibody Systems, Nonlinear Dynamics, and Control
  • Boston, Massachusetts, USA, August 2–5, 2015
  • Conference Sponsors: Design Engineering Division, Computers and Information in Engineering Division
  • ISBN: 978-0-7918-5716-8
  • Copyright © 2015 by ASME


The dynamics of large scale and complex multibody systems (MBS) that include flexible bodies and contact/impact pairs is governed by stiff equations. Because explicit integration methods can be very inefficient and often fail in the case of stiff problems, the use of implicit numerical integration methods is recommended in this case. This paper presents a new and efficient implementation of the two-loop implicit sparse matrix numerical integration (TLISMNI) method proposed for the solution of constrained rigid and flexible MBS differential and algebraic equations. The TLISMNI method has desirable features that include avoiding numerical differentiation of the forces, allowing for an efficient sparse matrix implementation, and ensuring that the kinematic constraint equations are satisfied at the position, velocity and acceleration levels. In this method, a sparse Lagrangian augmented form of the equations of motion that ensures that the constraints are satisfied at the acceleration level is first used to solve for all the accelerations and Lagrange multipliers. The generalized coordinate partitioning or recursive methods can be used to satisfy the constraint equations at the position and velocity levels. In order to improve the efficiency and robustness of the TLISMNI method, the simple iteration and the Jacobian-Free Newton-Krylov approaches are used in this investigation. The new implementation is tested using several low order formulas that include Hilber–Hughes–Taylor (HHT), L-stable Park, A-stable Trapezoidal, and A-stable BDF methods. The HHT method allow for including numerical damping. Discussion on which method is more appropriate to use for a certain application is provided. The paper also discusses TLISMNI implementation issues including the step size selection, the convergence criteria, the error control, and the effect of the numerical damping. The use of the computer algorithm described in this paper is demonstrated by solving complex rigid and flexible tracked vehicle models, railroad vehicle models, and very stiff structure problems. The results, obtained using these low order formulas, are compared with the results obtained using the explicit Adams-Bashforth predictor-corrector method. Using the TLISMNI method, which does not require numerical differentiation of the forces and allows for an efficient sparse matrix implementation, for solving complex and stiff structure problems leads to significant computational cost saving as demonstrated in this paper. In some problems, it was found that the new TLISMNI implementation is 35 times faster than the explicit Adams-Bashforth method.

Copyright © 2015 by ASME



Interactive Graphics


Country-Specific Mortality and Growth Failure in Infancy and Yound Children and Association With Material Stature

Use interactive graphics and maps to view and sort country-specific infant and early dhildhood mortality and growth failure data and their association with maternal

Citing articles are presented as examples only. In non-demo SCM6 implementation, integration with CrossRef’s "Cited By" API will populate this tab (http://www.crossref.org/citedby.html).

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In