4. Modélisation C#

4.1. Caractéristiques du maillage et de l’interface#

On considère une structure de dimensions \(\mathrm{LX}=5m\) , \(\mathrm{LY}=5m\) et \(\mathrm{LZ}=25m\) . Cette structure est discrétisée avec 5 mailles HEXA8. On s’intéresse à une interface plane de normale

\(n=(\begin{array}{}-1\\ 1\\ 1\end{array})\)

passant par le point \(A\) de coordonnées \((5,5\delta ,5)\) . La [] montre un zoom du 2ème élément où la trace de l’interface est représentée en rouge.

5

N9

_images/10000000000002D6000001E66F238399F736C296.png

Figure 4.1-a : Maillage c et zoom

L’interface est caractérisée par la level set normale ayant pour équation cartésienne:

\(\mathrm{lsn}=-x+y+z-5\delta\)

Remarque:

Avec le nouvel enrichissement de saut de déplacement [:ref:`R7.02.12 <R7.02.12>`], nous ne constatons pas de baisse de précision sur la solution en déplacement au voisinage de l'interface et de problème de conditionnement, bien que l'interface soit rasante par rapport aux nœuds du maillage.

4.2. Conditions aux limites#

Les conditions aux limites sont les mêmes que celles de la modélisation B. On encastre les nœuds de la face inférieure et on impose un déplacement de traction aux nœuds de la face supérieure:

Face inférieure: \(\mathrm{DX}=0\) , \(\mathrm{DY}=0\) et \(\mathrm{DZ}=0\)

Face supérieure: \(\mathrm{DX}=0\) , \(\mathrm{DY}=0\) et \(\mathrm{DZ}={10}^{-6}\)

4.3. Grandeurs testées et résultats#

Le bon déroulement du calcul permet a priori de valider le cas. On teste donc les valeurs du déplacement juste au dessus de l’interface après convergence des itérations de l’opérateur STAT_NON_LINE.

Identification

Référence

Tolérance

\(\mathit{DZ}\) pour tous les nœuds juste au dessus de l’interface

\(1.0E-6\)

1.0E-3%

Pour tester tous les nœuds en une seule fois, on teste le MINIMUM et le MAXIMUM de la colonne.