No abstract

No abstract

SUMMARYA common approach for the numerical simulation of non-linear multi-body contact problems is the use of Lagrange multipliers to model the contact conditions. The stability of standard algorithms is improved by introducing a modified mass matrix which assigns no mass to the potential contact nodes. By this, the spurious algorithmic oscillations in the multiplier do not occur any more, which facilitates the application of the primal-dual active set strategy to dynamical contact problems. The new mass matrix is calculated via a modified quadrature formula that needs no extra computational cost. In addition the conservation properties of the underlying algorithm are transferred to the modified mass version. Different numerical examples for frictional two-body contact problems illustrate the improvement in the results for the contact stresses.

In this paper, a new space-time discretization is proposed which is based on a modified mass matrix. The mass associated with a surface layer of elements is redistributed such that the inertia at the boundary is removed. This approach is motivated by the observation that standard space-time discretization schemes applied to dynamic contact problems yield spurious oscillations. A widely used approach for the numerical simulation of these problems is based on Lagrange multipliers which represent the contact stresses. But the algebraic contact conditions in combination with the inertia volume terms often yield nonphysical results for the contact stresses, and the stability of the algorithm can be lost. Our modified matrix is calculated via nonstandard quadrature formulas that require no extra computational effort. In addition, the conservation properties of the underlying algorithm are carried over to the modified method, and the standard optimal a priori estimates are still satisfied. Numerical examples confirm the optimality of the approach and its stabilization effect applied to contact problems.

No abstract

International audienceFrictional dynamic contact problems with complex geometries are a challenging task from the compu tational as well as from the analytical point of view since they generally involve space and time multi scale aspects. To be able to reduce the complexity of this kind of contact problem, we employ a non conforming domain decomposition method in space, consisting of a coarse global mesh not resolving the local struc ture and an overlapping fine patch for the contact computation. This leads to several benefits: First, we resolve the details of the surface only where it is needed, i.e., in the vicinity of the actual contact zone. Second, the subproblems can be discretized independently of each other which enables us to choose a much finer time scale on the contact zone than on the coarse domain. Here, we propose a set of interface conditions that yield optimal a priori error estimates on the fine meshed subdomain without any artificial dissipation. Further, we develop an efficient iterative solution scheme for the coupled problem that is robust with respect to jumps in the material parameters. Several complex numerical examples illustrate the performance of the new scheme

Inequality constraints occur in many different fields of application, e.g., in structural mechanics, flow processes in porous media or mathematical finance. In this paper, we make use of the mathematical structure of these conditions in order to obtain an abstract computational framework for problems with inequality conditions. The constraints are enforced locally by means of Lagrange multipliers which are defined with respect to dual basis functions. The reformulation of the inequality conditions in terms of nonlinear complementarity functions leads to a system of semismooth nonlinear equations that is solved by a generalized version of Newton's method for semismooth problems. By this, both nonlinearities in the pde model and inequality constraints are treated within a single Newton iteration which converges locally superlinear. The scheme can efficiently be implemented in terms of an active set strategy with local static condensation of the non-essential variables. Numerical examples from different fields of application illustrate the generality and the robustness of the method.

In order to obtain a fast solution scheme, the trajectory piecewise linear (TPWL) method is applied to the transient elastohydrodynamic (EHD) line contact problem for the first time. TPWL approximates the nonlinearity of a dynamical system by a weighted superposition of reduced linearized systems along specified trajectories. The method is compared to another reduced order model (ROM), based on Galerkin projection, Newton–Raphson scheme and an approximation of the nonlinear reduced system functions. The TPWL model provides further speed-up compared to the Newton–Raphson based method at a high accuracy.