35K85 Linear parabolic unilateral problems and linear parabolic variational inequalities [See also 35R35, 49J40]
Complete damage in linear elastic materials — Modeling, weak formulation and existence results
(2013)
In this work, we introduce a degenerating PDE system with a time-depending
domain for complete damage processes under time-varying
Dirichlet boundary conditions. The evolution of the system is
described by a doubly nonlinear differential inclusion for the damage
process and a degenerating quasi-static balance equation for the displacement field
which are strongly nonlinearly coupled.
In our proposed model, the material
may completely disintegrate which is indispensable for a realistic modeling of
damage processes in elastic materials. Complete damage theories
lead to several mathematical problems since, for instance, coercivity properties
of the free energy are lost and, therefore, several difficulties arise.
For the introduced complete damage model, we propose a classical
formulation and a corresponding suitable weak formulation in an
$SBV$-framework. The main aim is to prove existence of weak solutions
for the introduced degenerating model. In addition, we show that the classical
differential inclusion can be regained from the notion of weak solutions under
certain regularity assumptions which is a novelty in the theory of complete damage
models of this type.
For the existence results, we had to handle the following problem:
During the damage process it might occur that not completely damaged material regions are isolated
from the Dirichlet boundary. In this case, the
deformation field cannot be controlled in the transition from incomplete
to complete damage. To tackle this problem, we consider the evolution
process on a time-depending domain. In this context, two major challenges
arise:
Firstly, the time-dependent domain approach leads to jumps in the energy
which have to be accounted for in the energy inequality of the notion of
weak solutions. To handle this problem, several energy estimates are established
by $\Gamma$-convergence techniques. Secondly, the time-depending domain
might have bad smoothness properties such that Korn's inequality cannot be
applied. To this end, a covering result for such sets with smooth
compactly embedded domains has been shown.
A typical phase field approach for describing phase separation and
coarsening phenomena in alloys is the Cahn-Hilliard model. This
model has been generalized to the so-called Cahn-Larche system
by combining it with elasticity to capture non-neglecting deformation
phenomena,
which occur during phase separation and coarsening processes in the material.
In order to account for damage effects,
we extend the existing framework of Cahn-Hilliard and Cahn-Larche
systems by incorporating an internal damage variable of local character. This
damage variable allows to model the effect that damage of a material point is influenced by its local surrounding.
The damage process is described by a unidirectional rate-dependent evolution inclusion
for the internal variable.
For the introduced Cahn-Larche systems
coupled with rate-dependent
damage processes, we
establish a suitable notion of weak solutions and prove existence of weak solutions.
A model for the evolution of damage that allows for complete disintegration is addressed. Small strains and a linear response function are assumed. The ``flow rule'' for the damage parameter is rate-independent. The stored energy involves the gradient of the damage variable, which determines an internal length-scale. Quasi-static fully rate-independent evolution is considered as well as rate-dependent evolution including viscous/inertial effects. Illustrative 2-dimensional computer simulations are presented, too.
Quasistatic small-strain plasticity in the limit of small hardening and its numerical approximation
(2011)
The quasistatic rate-independent evolution of the Prager-Ziegler-type model of linearized plasticity with hardening is shown to converge to the rate-independent evolution of the Prandtl-Reuss elastic/perfectly plastic model. Based on the concept of energetic solutions we study the convergence of the solutions in the limit for hardening coefficients converging to 0 by using the abstract method of Gamma-convergence for rate-independent systems. An unconditionally convergent numerical scheme is devised and 2D and 3D numerical experiments are presented. A two-sided energy inequality is a posteriori verified to document experimental convergence rates.
A general abstract approximation scheme for rate-independent processes in the
energetic formulation is proposed and its convergence is proved under various
rather mild data qualifications. The abstract theory is illustrated on several
examples: plasticity with isotropic hardening, damage, debonding,
magnetostriction, and two models of martensitic transformation in
shape-memory alloys.