1. Reference problem#
1.1. Theoretical framework#
The unknowns of the problem are the degrees of freedom of movement and nodal damage. It is then a question of minimizing an energy of the form:
\(\phi (u,\alpha )\mathrm{=}\frac{1}{2}A(d)E{\epsilon }^{2}+\psi (d)+\frac{c}{2}\mathrm{\nabla }\alpha \mathrm{.}\mathrm{\nabla }\alpha\)
Where \(E\) is the Young’s modulus of the material, \(A(d)\) the stiffness function, the stiffness function, \(\psi (d)\) the dissipation, and \(c\) the non-local coefficient.
In the case of law ENDO_CARRE:
\(A(d)={(1-d)}^{2}\) and \(\psi (d)=\frac{{\sigma }_{y}^{2}}{E}d\)
The criterion corresponding to law ENDO_CARRE, for a homogeneous solution (\(\nabla \alpha =0\)), is therefore written:
\(d=1-(\frac{{W}_{y}}{{W}_{\mathrm{el}}})\)
Where \({W}_{\mathrm{el}}\) is the elastic deformation energy and:
\({W}_{y}=\frac{{\sigma }_{y}^{2}}{2E}\)
1.2. Geometry#
Consider a square with side \(L=1\text{m}\).
Dx
Figure 1: Representation of the problem
1.3. Material properties#
1.3.1. Law of damage: material ENDO_CARRE#
Elastic characteristics:
\(E=1\text{Pa}\)
\(\nu =0.\)
Characteristics related to the law of damage:
Elasticity limit:
\(\mathrm{SY}=0.01\text{Pa}\)
Non-local characteristics:
\(c=1.0\text{N}\)
1.4. Boundary conditions and loads#
Embed: Imposed movements zero DY = 0 \(m\) on the bottom and top horizontal edges (\(y=0.\) and \(y=1.\)) and DX = 0 \(m\) on the left edge (\(x=0.\)). See Figure 1.
Loading 1: Imposed displacement \({U}_{1}\) on the right vertical edge:
At moment \({t}_{1}\): DX = \(0.01\text{m}\)
At moment \({t}_{2}\): DX = \(0.0125\text{m}\)
At moment \({t}_{3}\): DX = \(0.02\text{m}\)