65K15 Numerical methods for variational inequalities and related problems
We present a globally convergent method for the solution of frictionless large deformation contact problems involving hyperelastic materials.
For the discretisation we apply the dual mortar method which is known to be more stable than node-to-segment approaches. The resulting non-convex constrained minimisation problems are solved using a filter–trust-region scheme.
This method combines several techniques from non-linear optimisation to achieve global convergence towards first-order optimal points. To speed up the method inexact linearisations of the non-penetration constraint are used
whenever the current iterate is far away from a critical point. A monotone multigrid method is applied for the fast solution of the constrained Newton
problems.
We introduce and analyze nonsmooth Schur-Newton methods for a class of nonsmooth saddle point problems. The method is able to solve problems where the primal energy decomposes into a convex smooth part and a convex separable but nonsmooth part. The method is based on nonsmooth Newton techniques for an equivalent unconstrained dual problem. Using this we show that it is globally convergent even for inexact evaluation of the linear subproblems.
Convergence Analysis of Smoothing Methods for Optimal Control of Stationary Variational Inequalities
(2012)
In the article an optimal control problem subject to a stationary variational inequality is investigated. The optimal control problem is complemented with pointwise
control constraints. The convergence of a smoothing scheme is analyzed. There, the variational inequality is replaced by a semilinear elliptic equation. It is shown that solutions of the regularized optimal control problem converge to solutions of the original one. Passing to
the limit in the optimality system of the regularized problem allows to prove C-stationarity of local solutions of the original problem. Moreover, convergence rates with respect to the regularization parameter for the error in the control are obtained. These rates coincide with
rates obtained by numerical experiments, which are included in the paper.
A capillary surface in a negative gravitational field describes the shape of the surface of a hanging drop in a capillary tube with wetting material on the bottom. Mathematical modeling leads to the volume- and obstacle-constrained minimization of a nonconvex nonlinear energy functional of mean curvature type which is unbounded from below. In 1984 Huisken proved the existence and regularity of local minimizers of this energy under the condition on gravitation being sufficiently weak. We prove convergence of a first order finite element approximation of these minimizers. Numerical results demonstrating the theoretic convergence order are given.
The aim of this paper is to devise an adaptive timestep control in the contact--stabilized Newmark method (CONTACX) for dynamical contact problems between two viscoelastic bodies in the framework of Signorini's condition. In order to construct a comparative scheme of higher order accuracy, we extend extrapolation techniques. This approach demands a subtle theoretical investigation of an asymptotic error expansion of the contact--stabilized Newmark scheme. On the basis of theoretical insight and numerical observations, we suggest an error estimator and a timestep selection which also cover the presence of contact. Finally, we give a numerical example.