Modélisation B ============== Caractéristiques du maillage ---------------------------- On discrétise la structure en 5 mailles de type HEXA8. 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. .. image:: images/10000000000000C3000001C07FFA573042BEA7F2.png :width: 1.5236in :height: 3.5in .. _RefImage_10000000000000C3000001C07FFA573042BEA7F2.png: **Figure** 3.1-a **: maillage avec 5** HEXA8 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:`\mathrm{N1}` , :math:`\mathrm{N6}` , :math:`\mathrm{N11}` , :math:`\mathrm{N16}` ): :math:`\mathrm{DX}=0` , :math:`\mathrm{DY}=0` et :math:`\mathrm{DZ}=0` Face supérieure (Nœuds :math:`\mathrm{N21}` , :math:`\mathrm{N22}` , :math:`\mathrm{N23}` , :math:`\mathrm{N24}` ): :math:`\mathrm{DX}=0` , :math:`\mathrm{DY}=0` et :math:`\mathrm{DZ}=\mathrm{uz}` Ceci constitue le 1er cas de chargement. En fait, on prend la liberté de déplacer la partie supérieure de la structure suivant les trois directions, on choisira donc comme 2ème cas de chargement: Face inférieure (Nœuds :math:`\mathrm{N1}` , :math:`\mathrm{N6}` , :math:`\mathrm{N11}` , :math:`\mathrm{N16}` ): :math:`\mathrm{DX}=0` , :math:`\mathrm{DY}=0` et :math:`\mathrm{DZ}=0` Face supérieure: :math:`\mathrm{DX}=\mathrm{ux}` , :math:`\mathrm{DY}=\mathrm{uy}` et :math:`\mathrm{DZ}=\mathrm{uz}` .. csv-table:: ":math:`\mathrm{ux}={1.10}^{-6}`" ":math:`\mathrm{uy}=2.{10}^{-6}`" ":math:`\mathrm{uz}=3.{10}^{-6}`" .. _Ref76289108: Résolution analytique --------------------- La solution d’un tel problème est bien sûr encore évidente. Tous les déplacements suivant :math:`x` et :math:`y` sont nuls, tous les déplacements suivant :math:`x` , :math:`y` et :math:`z` en dessous de la level set sont nuls et tous les déplacements suivant :math:`x` , :math:`y` et :math:`z` au dessus de la level set sont égaux au déplacement imposé :math:`{u}_{x}` , :math:`{u}_{y}` et :math:`{u}_{z}` 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:`3.3 ` ] pour les 2 cas de chargements. On obtient le tableau suivant pour le 1er cas de chargement. .. csv-table:: "**Identification**", "**Référence**", "**Tolérance**" ":math:`\mathit{DX}` pour tous les nœuds juste en dessous de l'interface", "0.00", "1.0E-16" ":math:`\mathit{DY}` pour tous les nœuds juste en dessous de l'interface", "0.00", "1.0E-16" ":math:`\mathit{DZ}` pour tous les nœuds juste en dessous de l'interface", "0.00", "1.0E-16" ":math:`\mathit{DX}` pour tous les nœuds juste au dessus de l'interface", "0.00", "1.0E-16" ":math:`\mathit{DY}` pour tous les nœuds juste au dessus de l'interface", "0.00", "1.0E-16" ":math:`\mathit{DZ}` pour tous les nœuds juste au dessus de l'interface", "1.0E-6", "1.0E-9%" On obtient le tableau suivant pour le 2ème cas de chargement. .. csv-table:: "**Identification**", "**Référence**", "**Tolérance**" ":math:`\mathit{DX}` pour tous les nœuds juste en dessous de l'interface", "0.00", "1.0E-16" ":math:`\mathit{DY}` pour tous les nœuds juste en dessous de l'interface", "0.00", "1.0E-16" ":math:`\mathit{DZ}` pour tous les nœuds juste en dessous de l'interface", "0.00", "1.0E-16" ":math:`\mathit{DX}` pour tous les nœuds juste au dessus de l'interface", "1.0E-6", "1.0E-9%" ":math:`\mathit{DY}` pour tous les nœuds juste au dessus de l'interface", "2.0E-6", "1.0E-9%" ":math:`\mathit{DZ}` pour tous les nœuds juste au dessus de l'interface", "3.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. 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.