3. Modélisation A#

3.1. Caractéristiques de la modélisation#

_images/10000000000004AE0000031090491E2344399B1A.png

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

  2. L’élément ponctuel \(Q\) (modélisation DIS_TR) permet de représenter simplement un accéléromètre présent dans le modèle de l’article [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.

3.2. Écriture des conditions aux limites#

Le fond du réservoir ne peut se déplacer que suivant la direction \(x\) :

CONDLIM=AFFE_CHAR_MECA( MODELE=MODELE, DDL_IMPO=(_F( GROUP_NO=(“FONDS”,”FONDP”,), DY=0.0, DZ=0.0,),),);

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

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:

PESA=AFFE_CHAR_MECA(    MODELE=MODELE,

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

3.3. Caractéristiques du maillage#

Le maillage contient :

24575 TETRA4

5120 PENTA6

4096 TRIA3

3.4. Valeurs testées#

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

Identification

Référence

\(\mathit{DX}(\mathit{N145},t=0,8s)\)

5.1624169321991e-04

\(\mathrm{DX}(\mathrm{N145},t=1,4s)\)

1.4970110314375e-04

\(\mathrm{DX}(\mathrm{N145},t=2,0s)\)

-2.3927413131721e-04

\(\mathrm{DX}(\mathrm{N3119},t=1,0s)\)

-9.9736272860105e-04

\(\mathrm{DX}(\mathrm{N3119},t=1,6s)\)

-8.7855056121762e-04

\(\mathrm{DX}(\mathrm{N3119},t=2,0s)\)

-2.3929161952584e-04

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

Identification

Référence

CHAR_CRIT

-2.477264942149