3. Modélisation A#
3.1. Caractéristiques de la modélisation#
C’est une modélisation PLAQ_MITC fait via fenicsx.
Le maillage de \(\Omega=[0,1]^2\) est directement construit par:
mesh = CA.Mesh.buildSquare(refine=FINE_REFINEMENT)
Conditions limites:
sur les bords \({0, 1}\times[0, 1] \cup [0, 1]\times{0, 1}\):
MECA_IMPO=_F(DZ=0, DRX=0, DRY=0, GROUP_MA=("LEFT", "RIGHT", "TOP", "BOTTOM"))
- Remarque:
A présent, PLAQ_MITC ne considère que la mouvement de type COQUE et on n’impose donc pas le déplacements DX, DY, DRZ.
Chargement:
PRES_REP=_F(GROUP_MA="SURFACE", PRES=PRES)
où PRES est une FORMULE qui représente l’équation dans le problème de référence.
3.2. Caractéristiques du maillage#
Nombre de nœuds : 16641
Nombre de mailles et types : 4352, 256 SEG3, 4096 QUAD9
3.3. Valeurs testées#
Localisation |
Épaisseur |
Type de valeur |
Référence |
Aster |
% différence |
Point \((0.5, 0.5)\) |
0.1 |
DZ |
\(9.254092261904761\cdot{10}^{-5}\) |
\(9.225572268464905\cdot{10}^{-5}\) |
< 1 |
Point \((0.5, 0.5)\) |
0.001 |
DZ |
\(8.13813244047619\cdot{10}^{-5}\) |
\(8.113126374037738\cdot{10}^{-5}\) |
< 1 |