Modélisation F ============== Caractéristiques du maillage ---------------------------- On discrétise la structure en 5 éléments finis QUAD4. Les nœuds de part et d’autre de l'interface sont des nœuds enrichis, donc les trois mailles centrales possédant de tels nœuds sont elles aussi enrichies. Seules les deux mailles extrêmes 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. .. image:: images/10000000000000C8000001F431D73E167571B76C.png :width: 1.1665in :height: 2.9165in .. _RefImage_10000000000000C8000001F431D73E167571B76C.png: **Figure** 7.1-a **: Maillage F** Conditions aux limites ---------------------- Les conditions aux limites appliquées représentent le même phénomène physique que pour la modélisation A. 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 :math:`\mathit{N1}` et :math:`\mathit{N2}` ): :math:`\mathrm{DX}=0` et :math:`\mathrm{DY}=0` Face supérieure (Nœuds :math:`\mathit{N3}` et :math:`\mathit{N4}` ): :math:`\mathrm{DX}=0` , :math:`\mathrm{DY}=\mathrm{uy}={10}^{-5}` . .. _Ref124332567: Résolution analytique --------------------- La solution d’un tel problème est bien sûr encore évidente : tous les déplacements suivant :math:`x` sont nuls, tous les déplacements suivant :math:`y` en dessous de la level set sont nuls et tous les déplacements suivant :math:`y` au dessus de la level set sont égaux au déplacement imposé :math:`{u}_{y}` au sommet de la structure. 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 [§ :ref:`7.3 ` ]. .. csv-table:: "**Identification**", "**Référence**", "**Tolérance**" "DX pour tous les nœuds juste en dessous de l'interface", "0.00", "1.0E-16" "DY pour tous les nœuds juste en dessous de l'interface", "0.00", "1.0E-16" "DX pour tous les nœuds juste au dessus de l'interface", "0.00", "1.0E-16" "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. On teste aussi les valeurs du déplacement issu de la commande POST_CHAM_XFEM. On teste en fait la valeur de la somme des valeurs absolues des déplacements des nœuds du maillage fissuré. C’est un test de non-régression par rapport aux valeurs obtenues avec la version 8.2.13 pour :math:`\mathit{DX}` et 9.0.21 pour :math:`\mathit{DY}` . .. csv-table:: "Identification", "Référence", "Différence" "SOMM_ABS(DX)", "0.000", "1.E-12" "SOMM_ABS(DY)", "1.3E-05", "1.0E-04%" Commentaires ------------ On remarque la discontinuité du champ de déplacements en traversant l'interface qui est possible grâce à l'enrichissement des éléments avec le degré de liberté Heaviside.