Modeling A
==============

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

The mesh of the beam is composed of a SEG2 element.

The mesh of the section is composed of 4 quadrilateral elements.

Node group :math:`\mathit{ENC}` is composed of node :math:`\mathit{P1}`.

Node group :math:`\mathit{CHA}` is composed of node :math:`\mathit{P2}`.


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

This is a POU_D_EM model (elements with 2 Gauss points). The spin angle is equal to :math:`0°`.

Modeling:

MOPOU = AFFE_MODELE (MAILLAGE = MAPOU,

 AFFE =_F (TOUT =' OUI ', PHENOMENE =' MECANIQUE',

     MODELISATION =' POU_D_EM ',),

)

Boundary conditions:

BLOCAGE = AFFE_CHAR_MECA (MODELE =MO,

 DDL_IMPO =_F (GROUP_NO =' ENC ',

     DX=0.0, DY=0.0, DZ=0.0, DRX =0.0, =0.0, DRY =0.0, DRZ =0.0,),

)

Mechanical loading:

CHARGE = AFFE_CHAR_MECA (MODELE =MO,

 FORCE_NODALE =_F (GROUP_NO =' CHA ', FX=100, FY=100, FZ=100. ),

)

Assignment of the characteristics of the elements:

POUCA_0 = AFFE_CARA_ELEM (MODELE = MOPOU, GEOM_FIBRE =GF,

 POUTRE =_F (GROUP_MA =( 'POUTRE'),

     SECTION =' RECTANGLE ', CARA =(' HY ',' HZ '), VALE =( 0.2,0.1),

     PREC_AIRE =5. , PREC_INERTIE =10.,),

 ORIENTATION =_F (GROUP_MA =( 'POUTRE'), CARA =' ANGL_VRIL ', VALE =0.0,),

 MULTIFIBRE =_F (GROUP_MA =( 'POUTRE'), GROUP_FIBRE =( 'SBET',),),

)


Tested values and results
----------------------------

.. csv-table::

    "Coordinated
    mesh SG01 ", "Integration point", "Sub-point", "**Reference**"
    "COOR_X ", "2", "1", "1", "0.377088184"
    "COOR_X ", "2", "2", "2", "0.447798863"
    "COOR_X ", "2", "3", "3", "0.468211277"
    "COOR_X ", "2", "4", "4", "0.397500599"
    "COOR_Y ", "2", "1", "1", "0.447798863"
    "COOR_Y ", "2", "2", "2", "0.377088184"
    "COOR_Y ", "2", "3", "3", "0.397500599"
    "COOR_Y ", "2", "4", "4", "0.468211277"
    "COOR_Z ", "2", "1", "1", "0.443062145"
    "COOR_Z ", "2", "2", "2", "0.443062145"
    "COOR_Z ", "2", "3", "3", "0.402237316"
    "COOR_Z ", "2", "4", "4", "0.402237316"


All the tolerances are the ones by default: :math:`1.0E-3`