C modeling
==============

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_TGM model (elements with 3 Gauss points). The spin angle is :math:`0°`.

Modeling:

MOPOU = AFFE_MODELE (MAILLAGE = MAPOU,

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

     MODELISATION =' POU_D_TGM ',),

)

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 ", "Point
    integration", "Sub-point", "**Reference**"
    "COOR_X ", "2", "1", "1", "0.954438454"
    "COOR_X ", "2", "2", "2", "1.025149132"
    "COOR_X ", "2", "3", "3", "1.045561546"
    "COOR_X ", "2", "4", "4", "0.974850868"
    "COOR_Y ", "2", "1", "1", "1.025149132"
    "COOR_Y ", "2", "2", "2", "0.954438454"
    "COOR_Y ", "2", "3", "3", "0.974850868"
    "COOR_Y ", "2", "4", "4", "1.045561546"
    "COOR_Z ", "2", "1", "1", "1.020412415"
    "COOR_Z ", "2", "2", "2", "1.020412415"
    "COOR_Z ", "2", "3", "3", "0.979587585"
    "COOR_Z ", "2", "4", "4", "0.979587585"


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