Modélisation D =============== Cette modélisation est basée sur la modélisation A. Le type d’élément choisi pour le maillage est la seule différence entre ces deux modélisations. Caractéristiques du maillage ---------------------------- On discrétise la structure en 6 éléments finis TETRA4 L'interface est présente au sein de ces 6 éléments par le biais des level sets. .. image:: images/1000000000000239000004238C5733C18C337071.jpg :width: 2.9634in :height: 4.302in .. _RefImage_1000000000000239000004238C5733C18C337071.jpg: **Figure** 5.1-a **: Maillage** Conditions aux limites ---------------------- Les conditions aux limites sont celles de 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. Résolution analytique --------------------- La solution analytique est celle présentée dans la modélisation A [§ :ref:`2.3 ` ] : tous les degrés de liberté suivant :math:`x` et :math:`y` sont nuls et tous les degrés de liberté suivant :math:`z` valent :math:`\mathrm{uz}/2` , où :math:`\mathrm{uz}={10}^{-6}` Grandeurs testées et résultats ------------------------------ On teste les valeurs du déplacement juste en dessous et au dessus de l'interface 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:`2.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" "DZ 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", "0.00", "1.0E-16" "DZ 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. 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.