9. Modélisation H#
Cette modélisation est basée sur la modélisation F.
Le type d’élément choisi pour le maillage est la seule différence entre ces deux modélisations.
9.1. Caractéristiques du maillage#
Chaque maille QUAD4 de la modélisation F est décomposée en 2 TRIA3 pour la modélisation H.
Ainsi la structure est discrétisée en 10 éléments finis TRIA3.
Les nœuds de part et d’autre de l’interface sont des nœuds enrichis, donc les triangles contenus dans les trois mailles centrales de la modélisation F possédant de tels nœuds sont eux aussi enrichis. Seuls les triangles contenus dans les deux quadrilatères extrêmes de la modélisation F sont des mailles classiques n’ayant que des nœuds classiques.
On pourra donc imposer des conditions aux limites sur les mailles extrêmes de la manière habituelle.
Figure 9.1-a : Maillage H
9.2. Conditions aux limites#
On encastre les nœuds de la face inférieure et on impose un déplacement des nœuds de la face supérieure:
Face inférieure (Nœuds \(\mathrm{N1}\) , \(\mathrm{N2}\) ): \(\mathrm{DX}=0\) , et \(\mathrm{DY}=0\)
Face supérieure (Nœuds \(\mathrm{N3}\) , \(\mathrm{N4}\) ): \(\mathrm{DX}=0\) et \(\mathrm{DY}=\mathrm{uy}\)
9.3. Résolution analytique#
La solution d’un tel problème est bien sûr encore évidente : tous les déplacements suivant \(x\) sont nuls, tous les déplacements suivant \(y\) en dessous de la level set sont nuls et tous les déplacements suivant \(y\) au dessus de la level set sont égaux au déplacement imposé \({u}_{y}\) au sommet de la structure.
9.4. Grandeurs testées et résultats#
On teste les valeurs du déplacement après convergence des itérations de l’opérateur STAT_NON_LINE. On vérifie que l’on retrouve bien les valeurs déterminées au [§9.3].
Identification |
Référence |
Tolérance |
\(\mathit{DX}\) pour tous les nœuds juste en dessous de l’interface |
0.00 |
1.0E-16 |
\(\mathit{DY}\) pour tous les nœuds juste en dessous de l’interface |
0.00 |
1.0E-16 |
\(\mathit{DX}\) pour tous les nœuds juste au dessus de l’interface |
0.00 |
1.0E-16 |
\(\mathit{DY}\) pour tous les nœuds juste au dessus de l’interface |
1.0E-6 |
1.0E-9% |
Pour tester tous les nœuds en une seule fois, on teste le MINIMUM et le MAXIMUM de la colonne.
9.5. Commentaires#
Pour tester tous les nœuds en une seule fois, on teste le MINIMUM et le MAXIMUM de la colonne.