Timestepping schemes for nonsmooth dynamics based on discontinuous Galerkin methods: Definition and outlook

Abstract

International audienceThe contribution deals with timestepping schemes for nonsmooth dynamical systems. Traditionally, these schemes are locally of integration order one, both in non-impulsive and impulsive periods. This is inefficient for applications with infinitely many events but large non-impulsive phases like circuit breakers, valve trains or slider-crank mechanisms. To improve the behaviour during non-impulsive episodes, we start activities twofold. First, we include the classic schemes in time discontinuous Galerkin methods. Second, we split non-impulsive and impulsive force propagation. The correct mathematical setting is established with mollifier functions, Clenshaw-Curtis quadrature rules and an appropriate impact representation. The result is a Petrov-Galerkin distributional differential inclusion. It defines two Runge-Kutta collocation families and enables higher integration order during non-impulsive transition phases. As the framework contains the classic Moreau-Jean timestepping schemes for constant ansatz and test functions on velocity level, it can be considered as a consistent enhancement. An experimental convergence analysis with the bouncing ball example illustrates the capabilities

Nonsmooth Contact Dynamics for the large-scale simulation of granular material

“…The DVI (6) is a variational formulation in a function space setting and resembles the weak form of a partial differential equation where the spacial descretization is achieved by the Galerkin projection. A similar approach is used here in time rather than space, based on [18]. As starting point, consider a one dimensional problem where q : [t 0 , t 1 ] → R is absolutely continuous and…”

A conical interior point method for nonsmooth rigid body dynamics

“…A stable, semi-implicit time stepping scheme is derived from the equations of motion (4) using a Petrov-Galerkin method [8]. We identify finite dimensional ansatz functions for v, λ ∈ BV(R) and finite dimensional test functions δ q ∈ AC(R) and satisfy (4) in these subspaces.…”

Time finite element based Moreau‐type integrators

“…While the algorithm of Fetecau et al belongs to the class of event‐driven algorithms, Johnson et al proposed an approach allowing for discontinuous trajectories. The same discontinuous behavior of the trajectories was assumed in the work of Schindler and Acary, where the discretization of a variational formulation of the equality of measures with a discontinuous Galerkin approach in time was discussed. This leads then to Moreau's time‐stepping scheme together with some higher‐order schemes.…”