B modeling
==============


Characteristics of modeling
-----------------------------------

3D modeling (linear elements) is used.


.. image:: images/Object_4.svg
    :width: 641
    :height: 362


.. _RefImage_Object_4.svg:

The discrete contact formulation is used (active stress method).


Characteristics of the mesh
----------------------------


.. image:: images/Object_5.png
    :width: 4.6102in
    :height: 2.2854in


.. _RefImage_Object_5.png:

Number of knots: 3292

Number of meshes: 2150 HEXA8, 260 PENTA6, 260, 1814 QUAD4 and 68 TRIA3

Number of nodes in contact: 210


Tested sizes and results
------------------------------

The following displacement :math:`Y` of the point :math:`K` of the piece of paper with respect to the surface :math:`\mathrm{ABCDEFG}` of the die is tested.

.. csv-table::

    "**Identification**", "**Reference type**", "**Reference value**", "**Tolerance**"
    "Point :math:`K`/Point :math:`B` - :math:`\mathrm{DY}` ", "'ANALYTIQUE'", "5.0000"," 0.1%"
    "Point :math:`K`/Point :math:`C` - :math:`\mathrm{DY}` ", "'ANALYTIQUE'", "20.8250"," 1.0%"
    "Point :math:`K`/Point :math:`D` - :math:`\mathrm{DY}` ", "'ANALYTIQUE'", "55.8800"," 1.0%"
    "Point :math:`K`/Point :math:`E` - :math:`\mathrm{DY}` ", "'ANALYTIQUE'", "78.6900"," 1.0%"
    "Point :math:`K`/Point :math:`F` - :math:`\mathrm{DY}` ", "'ANALYTIQUE'", "144.8950"," 1.0%"
    "Point :math:`K`/Point :math:`G` - :math:`\mathrm{DY}` ", "'ANALYTIQUE'", "155.0960"," 1.0%"


We test the number of iterations of Newton when the point :math:`K` of the plot is in front of the points :math:`B`, :math:`C`, :math:`D`, :math:`E`, :math:`F` and :math:`G` of the string.


.. csv-table::

    "**Identification**", "**Reference type**", "**Reference value**", "**Tolerance**"
    "ITER_GLOB/Point :math:`B` ", "'NON_REGRESSION'", "2"," 0.00%"
    "ITER_GLOB/Point :math:`C` ", "'NON_REGRESSION'", "2"," 0.00%"
    "ITER_GLOB/Point :math:`D` ", "'NON_REGRESSION'", "2"," 0.00%"
    "ITER_GLOB/Point :math:`E` ", "'NON_REGRESSION'", "2"," 0.00%"
    "ITER_GLOB/Point :math:`F` ", "'NON_REGRESSION'", "2"," 0.00%"
    "ITER_GLOB/Point :math:`G` ", "'NON_REGRESSION'", "2"," 0.00%"



notes
---------


* The calculation is carried out by imposing a displacement on the back side of the block :math:`\text{(MN)}`. The displacement is imposed as follows:


.. csv-table::

    "* from", ":math:`\mathrm{0.mm}` to", ":math:`\mathrm{5.mm}` in 5 steps"
    "* from", ":math:`\mathrm{5.mm}` to", ":math:`\mathrm{15.mm}` in 10 steps"
    "* from", ":math:`\mathrm{5.mm}` to", ":math:`\mathrm{20.mm}` in 5 steps"
    "* from", ":math:`\mathrm{20.mm}` to", ":math:`\mathrm{50.mm}` in 5 steps"
    "* from", ":math:`\mathrm{50.mm}` to", ":math:`\mathrm{70.mm}` in 10 steps"
    "* from", ":math:`\mathrm{70.mm}` to", ":math:`\mathrm{140.mm}` in 35 steps"
    "* from", ":math:`\mathrm{140.mm}` to", ":math:`\mathrm{155.mm}` in 15 steps"



* The calculations do not converge with the default keywords used for convergence in STAT_NON_LINE, because the default value of RESI_GLOB_RELA = 1.E-6 is too restrictive (the forces to which the piece is subjected are initially relatively weak). To overcome this problem, use the keyword RESI_GLOB_MAXI = 1.E-6.