Modélisation A
==============


Caractéristiques de la modélisation
-----------------------------------


.. image:: images/10000000000004AE0000031090491E2344399B1A.png
    :width: 6.889in
    :height: 4.5075in


.. _RefImage_10000000000004AE0000031090491E2344399B1A.png:


1.  La maillage total comporte 8163 nœuds, soit environ 125000 degrés de liberté,


2.  L’élément ponctuel :math:`Q` (modélisation DIS_TR) permet de représenter simplement un accéléromètre présent dans le modèle de l’article [:ref:`bib1 <bib1>`],


3.  La plaque déformable est modélisée par 5120 éléments de solide massif (modélisation 3D) pentaèdriques à 6 nœuds (10 couches dans l’épaisseur pour une bonne approximation du comportement en flexion malgré la linéarité des éléments),


4.  la surface libre est modélisée par 512 éléments MEFP_FACE3 (modélisation 2D_FLUI_PESA) triangles à 3 nœuds,


5.  le volume fluide est modélisé par 24576 éléments de fluide (modélisation 3D_FLUIDE) tétraèdriques à 4 nœuds.


Écriture des conditions aux limites
-----------------------------------

Le fond du réservoir ne peut se déplacer que suivant la direction :math:`x` :


.. code-block:: text
CONDLIM=AFFE_CHAR_MECA(    MODELE=MODELE,
DDL_IMPO=(_F(
GROUP_NO=('FONDS','FONDP',),
DY=0.0, DZ=0.0,),),);

Suivant cette direction :math:`x` , on impose au fond du réservoir un déplacement sinusoïdal en temps, de fréquence :math:`1,7704\mathrm{Hz}` et d’amplitude :math:`0,001m` :


.. code-block:: text
FREQ = 1.7704;
LFONC=DEFI_LIST_REEL(DEBUT=0.0, INTERVALLE=_F( JUSQU_A=10.0,
PAS=0.01,),);
FONC = FORMULE(    REEL = '''(REEL:INST) =
(0.001)*SIN(2*PI*FREQ*INST)''');
DEPLX=CALC_FONC_INTERP(        FONCTION=FONC,
NOM_PARA='INST',
LIST_PARA=LFONC,);
CHARG_SE=AFFE_CHAR_MECA_F(    MODELE=MODELE,
DDL_IMPO=_F(
GROUP_NO=('FONDS','FONDP',),DX=DEPLX,),);

Le chargement volumique de pesanteur est défini ainsi:


.. code-block:: text


    PESA=AFFE_CHAR_MECA(    MODELE=MODELE,
PESANTEUR=_F(    GRAVITE=9.81,
DIRECTION=(0.,0.,-1.,),),);


Caractéristiques du maillage
----------------------------


.. csv-table::

    "Le maillage contient :", "24575 TETRA4"
    "", "5120 PENTA6"
    "", "4096 TRIA3"



Valeurs testées
---------------

Les tests se font sur la valeur du déplacement suivant :math:`x` (notée :math:`\mathit{DX}` ) pour différents instants et pour les nœuds :math:`\mathit{N145}` et :math:`\mathit{N3119}` .


.. csv-table::

    "**Identification**", "**Référence**"
    ":math:`\mathit{DX}(\mathit{N145},t=0,8s)`", "5.1624169321991e-04"
    ":math:`\mathrm{DX}(\mathrm{N145},t=1,4s)`", "1.4970110314375e-04"
    ":math:`\mathrm{DX}(\mathrm{N145},t=2,0s)`", "-2.3927413131721e-04"
    ":math:`\mathrm{DX}(\mathrm{N3119},t=1,0s)`", "-9.9736272860105e-04"
    ":math:`\mathrm{DX}(\mathrm{N3119},t=1,6s)`", "-8.7855056121762e-04"
    ":math:`\mathrm{DX}(\mathrm{N3119},t=2,0s)`", "-2.3929161952584e-04"


On teste aussi la valeur propre critique calculée avec CRIT_STAB.


.. csv-table::

    "**Identification**", "**Référence**"
    "CHAR_CRIT", "-2.477264942149"