000015979 001__ 15979
000015979 005__ 20161115135335.0
000015979 04107 $$aeng
000015979 046__ $$k2013-06-12
000015979 100__ $$aEsefeld, B.
000015979 24500 $$aCoupling of Ode- and Dae-Integration With the Time-Stepping Integration Method for Multibody Systems With Impacts and Friction

000015979 24630 $$n34.$$pComputational Methods in Structural Dynamics and Earhquake Engineering
000015979 260__ $$bNational Technical University of Athens, 2013
000015979 506__ $$arestricted
000015979 520__ $$2eng$$aMultibody systems with unilateral constraints feature difficulties in the numerical integration because of their non-smooth characteristics. Depending on the number of unilateral contacts and their state transitions either time-stepping or event-driven schemes are commonly applied [1]. The hybrid method presented in [2] couples Moreau’s half-explicit time-stepping method on velocity level with a Livermore Solver for ordinary differential equations (LSODAR). It profits from the respective advantages and avoids many drawbacks. As an extension, the hybrid integration scheme [2] is complemented by an integration method for differential algebraic equations (DASKR). This optionally allows solving the non-impulsive phases as an index-2 DAE by replacing LSODAR with DASKR. Therefore, the non-impulsive contact conditions are formulated on velocity level. Hence, they are fulfilled on the same kinematic level as in the time-stepping mode reserved for impulsive transitions. In contrast to multibody systems with only bilateral constraints, the algebraic equations become more complicated due to possible state transitions. Complementarity conditions for forces and kinematic quantities are expressed via a proximal point formulation which transforms the inequality conditions into a system of non-smooth equations. This can be solved by specially designed root solving method like non-smooth variants of Newton’s method [3].For making the Newton iteration more robust, the corresponding iteration matrix can be provided partly analytically. Particularly, this includes analytic representations of the discontinuous Jacobians in the proximal point formulation. The behavior of the semi-analytical approach is compared to the fully numerical approximation of the iteration matrix provided automatically by DASKR. ODE- and DAE-type hybrid integration schemes are implemented in the multibody simulation software MBSim [4]. Test examples are used for a comparison. Therefore, the constraint conditions are regarded on different kinematic levels, as well as necessities to switch between time-stepping and non-impulsive integration modes, re-initialization procedures and numerical differences in the solutions are discussed.

000015979 540__ $$aText je chráněný podle autorského zákona č. 121/2000 Sb.
000015979 653__ $$a

000015979 7112_ $$aCOMPDYN 2013 - 4th International Thematic Conference$$cIsland of Kos (GR)$$d2013-06-12 / 2013-06-14$$gCOMPDYN2013
000015979 720__ $$aEsefeld, B.$$iSchindler, T.$$iUlbrich, H.
000015979 8560_ $$ffischerc@itam.cas.cz
000015979 8564_ $$s43570$$uhttps://invenio.itam.cas.cz/record/15979/files/2055.pdf$$yOriginal version of the author's contribution as presented on CD, section: SC-MS 15 COMPUTATIONAL CONTACT MECHANICS
.
000015979 962__ $$r15525
000015979 980__ $$aPAPER