Journal
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING
Volume 330, Issue -, Pages 351-367Publisher
ELSEVIER SCIENCE SA
DOI: 10.1016/j.cma.2017.10.030
Keywords
Interior point; DVI; Rigid contact; Multibody
Ask authors/readers for more resources
Time-stepping methods for non-smooth dynamics are based on the solution of multiple complementarity problems: solving this class of problems represents a major numerical bottleneck, especially when dealing with near-singular and ill-posed systems such as those which feature large mass ratios between mechanical parts. To overcome these difficulties, we propose the adoption of a log-barrier method whose convergence properties outperform classical methods used in literature. This leads to a primal-dual interior point method, based on a custom Mehrotra predictor-corrector scheme, that is able to operate on the typical sparse matrices of multibody dynamics and that can accommodate tangent stiffness matrices arising from flexible bodies, if any. Although our original method addresses unilateral constraints only, it is possible to formulate it in order to include equality constraints, for example when adding joints. We implemented the algorithm in C++ in the open-source Project Chrono. Results show that the method converges to small residuals even in ill-posed scenarios that downgrade the convergence of other methods. Finally, for the sake of high efficiency, we tested different options for the linear solvers used in the inner iteration. (C) 2017 Elsevier B.V. All rights reserved.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available