Reference problem
=====================


Geometry
---------

According to the modeling :math:`\mathrm{2D}` (plane deformation) or :math:`\mathrm{3D}`, we consider respectively a rectangle or a bar (see Figure 1.1-1).


.. image:: images/100002010000025C000001B9E3FEB8F31B7BFF89.png
    :width: 6.2047in
    :height: 4.2917in


.. _RefImage_100002010000025C000001B9E3FEB8F31B7BFF89.png:

**Figure** 1.1-a **: Geometry and Meshing of a rectangle (a) and a bar (b)**


Material properties
----------------------

Elasticity:


.. csv-table::

    ":math:`E=190000\text{MPa}` ", "Young's module"
    ":math:`\mathrm{\nu }=0.3` ", "Poisson's ratio"


Work hardening curve:

:math:`R(\mathrm{\kappa })=488.36+57.13(1-\mathrm{exp}(-8613\mathrm{\kappa }))+238.73(1-\mathrm{exp}(-10.39\mathrm{\kappa }))`

Ductile damage law GTN:


.. csv-table::

    ":math:`{q}_{1}=1.5` ", "Model parameter GTN"
    ":math:`{q}_{2}=1.07` ", "Model parameter GTN"
    ":math:`{f}_{0}=0.0002i` ", "Initial porosity"
    ":math:`{f}_{n}=0` ", "Germination Parameter"
    ":math:`{f}_{c}=0.05` ", "Coalescence porosity"
    ":math:`\mathrm{\delta }=3` ", "Coalescence coefficient related to coalescence"
    ":math:`c=1N` ", "Non-local parameter"
    ":math:`r=5000\mathit{MPa}` ", "Lagrange penalty parameter"


In particular, the distribution of the initial porosity is not homogeneous. It increases with altitude: :math:`{f}_{0}=0.0002i` where :math:`i` refers to the i-th element layer and :math:`1\le i\le 10`.

In DEFI_MATERIAU, the following information should be filled in:


.. csv-table::

    "**ELAS**", "**ECRO_** NL", "**GTN**", "**NON_LOCAL**", "****"
    "E = 190000", "R0= 488.361123569", "Q1 = 1.5", "C_ GRAD_VARI = 1"
    "NU = 0.3", "R1 = 57.1333673502", "Q2 = 1.07"," PENA_LAGR =5000"
    "", "GAMMA_1 = 8613"," PORO_INIT = 0.0002i", ""
    "", "R2 = 238.731127339"," COAL_PORO = 0.05", ""
    "", "GAMMA_2 =10.386585592"," COAL_ACCE = 3", ""



Boundary conditions and loads
-------------------------------------

For modeling :math:`2D` (plane deformation), the displacements along the :math:`X` axis of all the nodes are controlled: :math:`{u}_{x}=5x`, the displacements along the :math:`Y` axis of the nodes of the same altitude remain uniform, the node :math:`A` is blocked along the :math:`Y` axis (see Figure 1.1-1 (a) for geometry).

For modeling :math:`3D`, the movements along the :math:`X` axis of all the nodes are controlled: :math:`{u}_{x}=5x`, the movements along the :math:`Y` axis of all the nodes are blocked, the movements along the :math:`Z` axis of the nodes of the same altitude remain uniform, the node :math:`A` is blocked according to :math:`Z`. (see Figure 1.1-1 (b) for geometry).

Boundary conditions and loads are imposed this way so that the problem in :math:`2D` and the problem in :math:`3D` are the same.

The load is imposed using 1000 identical time steps. The pseudo-calculation time is 1.