Refine
Year of publication
- 2012 (7) (remove)
Language
- English (7)
Keywords
- coupling conditions (2)
- energy minimization (2)
- optimal control (2)
- path-following (2)
- C-stationarity (1)
- Cosserat rod (1)
- Cosserat rods (1)
- PDE constrained optimization (1)
- Variational inequalities (1)
- contact problem (1)
Project
- A17 (7) (remove)
Application Area
- A (7)
We formulate the static mechanical coupling of a geometrically exact Cosserat rod
to a nonlinearly elastic continuum. In this setting, appropriate coupling conditions have
to connect a one-dimensional model with director variables to a three-dimensional
model without directors.
Two alternative coupling conditions are proposed,
which correspond to two different configuration trace spaces.
For both we show existence of solutions of the coupled problems, using the direct
method of the calculus of variations. From the first-order optimality conditions
we also derive the corresponding conditions for the dual variables. These are
then interpreted in mechanical terms.
We present a time-dependent finite element model of the human knee joint of full 3D geometric complexity together with advanced numerical algorithms needed for its simulation. The model comprises bones, cartilage and the major ligaments, while patella and menisci are still missing. Bones are modeled by linear elastic materials, cartilage by linear viscoelastic materials, and ligaments by one-dimensional nonlinear Cosserat rods. In order to capture the dynamical contact problems correctly, we solve the full PDEs of elasticity with strict contact inequalities. The spatio--temporal discretization follows a time layers approach (first time, then space discretization). For the time discretization of the elastic and viscoelastic parts we use a new contact-stabilized Newmark method, while for the Cosserat rods we choose an energy--momentum method. For the space discretization, we use linear finite elements for the elastic and viscoelastic parts and novel geodesic finite elements for the Cosserat rods. The coupled system is solved by a Dirichlet--Neumann method. The large algebraic systems of the bone--cartilage contact problems are solved efficiently by the truncated non-smooth Newton multigrid method.
We derive a-priori estimates on the length of the primal-dual path that results from a Moreau-Yosida approximation of the feasible set for state constrained optimal control problems. These bounds depend on the regularity of the state and the dimension of the problem. Comparison with numerical results indicates that these bounds are sharp and are attained for the case of a single active point.
We formulate the static mechanical coupling of a geometrically exact Cosserat rod to an elastic continuum. The coupling conditions accommodate for the difference in dimension between the two models. Also, the Cosserat rod model incorporates director variables, which are not present in the elastic continuum model. Two alternative coupling conditions are proposed, which correspond to two different configuration trace spaces. For both we show existence of solutions of the coupled problems. We also derive the corresponding conditions for the dual variables and interpret them in mechanical terms.
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.
We consider a shape implant design problem that arises in the context of facial surgery. We introduce a reformulation as an optimal control problem, where the control acts as a boundary force. The state is modelled as a minimizer of a polyconvex hyperelastic energy functional. We show existence of optimal solutions and derive - on a formal level - first order optimality conditions. Finally, preliminary numerical results are presented.
We present a discretization for dynamic large deformation contact problems without friction. Our model is based on Hamilton’s principle, which avoids the explicit appearance of the contact forces. The resulting differential inclusion is discretized in time using a modified midpoint rule. This modification, which concerns the evaluation of the generalized gradient, allows to achieve energy dissipativity. For the space discretization we use a dual-basis mortar method. The resulting spatial algebraic
problems are nonconvex minimization problems with nonconvex inequality constraints. These can be solved efficiently using a trust-region SQP framework with a monotone multigrid inner solver.