v6.04.219 SSNV219 – Méthode des solutions manufacturées en contact 3D et grandes déformations#

Résumé:

L’objectif de ce test est de vérifier la modélisation du contact 3D en grandes déformations grâce à la méthode des solutions manufacturées [bib1].

Solution de référence#

Méthode de calcul#

La solution de référence analytique est donnée par :

(4805)#\[\begin{split}\begin{array}{ccc}\mathit{Ux}& =& 0.2\times {Z}^{2}\times X\times Y\\ \mathit{Uy}& =& 0.2\times {Z}^{2}\times X\times Y\\ \mathit{Uz}& =& -0.2\times (1.+{X}^{2}+{Y}^{2})\times (1+0.01\times Z)-0.01\mathrm{\ast }Z-0.3\end{array}\end{split}\]

Les conditions de Dirichlet, de Neumann et le terme source sont obtenus par la méthode des solutions manufacturées [bib1].

On commence par déterminer le gradient de la transformation \(\underline{\underline{F}}\) :

(4806)#\[\underline{\underline{F}}=\underline{\nabla}\underline{U}+\underline{\underline{\mathrm{Id}}}\]

Connaissant la normale \(\underline{N}={[0,0,-1]}^{T}\) à la surface \(\text{ESCLAVE}\) dans la configuration non-déformée, on obtient son expression dans la configuration déformée par la formule de Nanson :

(4807)#\[\underline{n}=\frac{{\underline{\underline{F}}}^{-T}\underline{N}}{\parallel {\underline{\underline{F}}}^{-T}\underline{N}\parallel }\]

Connaissant le tenseur de Hooke \(\underline{\underline{\underline{\underline{A}}}}\) et le tenseur de Green-Lagrange \(\underline{\underline{E}}\) , on calcule le second tenseur de Piola-Kirchhoff \(\underline{\underline{S}}\) :

(4808)#\[\underline{\underline{E}}=\frac{1}{2}({\underline{\underline{F}}}^{T}\cdot \underline{\underline{F}}-\underline{\underline{\mathrm{Id}}})\]
(4809)#\[\underline{\underline{S}}=\underline{\underline{\underline{\underline{A}}}}:\underline{\underline{E}}\]

On rappelle que le second tenseur de Piola-Kirchhoff \(\underline{\underline{S}}\) permet d’obtenir des efforts en configuration non déformée par unité de surface non déformée :

(4810)#\[\frac{d\underline{{f}_{0}}}{\mathrm{dA}}=\underline{\underline{S}}\cdot \underline{N}\]

Comme nous cherchons à déterminer des efforts en configuration déformée, nous allons déterminer le premier tenseur de Piola-Kirchhoff \(\underline{\underline{\Pi}}\)

(4811)#\[\underline{\underline{\Pi}}=\underline{\underline{F}}\cdot \underline{\underline{S}}\]

On peut ainsi déterminer les forces de volume \({\underline{f}}_{\mathrm{vol}}\) :

(4812)#\[{\underline{f}}_{\mathrm{vol}}=\text{-div}\underline{\underline{\Pi}}\]

Connaissant la normale en configuration initiale sur les différentes faces et le premier tenseur de Piola-Kirchhoff \(\underline{\underline{\Pi}}\) , on peut calculer les efforts de surface en configuration déformée :

(4813)#\[{\underline{f}}_{\text{surf}}=\underline{\underline{\Pi}}\cdot \underline{N}\]

Sur la surface \(\text{BAS}\) qui est en contact, il faut un traitement particulier. En effet, les efforts normaux y sont pris en compte par le contact :

(4814)#\[\begin{split}\begin{array}{ccccc}{\underline{f}}_{\text{surf}}^{\text{BAS}}& =& {\underline{f}}_{{\text{surf}}_{n}}^{\text{BAS}}& +& {\underline{f}}_{{\text{surf}}_{t}}^{\text{BAS}}\\ & =& {\underline{f}}_{\text{contact}}& +& {\underline{f}}_{{\text{surf}}_{t}}^{\text{BAS}}\\ & =& p\ast \underline{n}& +& {\underline{f}}_{{\text{surf}}_{t}}^{\text{BAS}}\end{array}\end{split}\]

Où \(p\) désigne la pression de contact. Elle peut être déterminée par l’expression :

(4815)#\[p=(\underline{\underline{\Pi}}\cdot \underline{N})\cdot \underline{n}\]

Il ne faut donc y appliquer que les efforts tangentiels. On les calcule par l’expression :

(4816)#\[\begin{split}\begin{array}{ccc}{\underline{f}}_{{\text{surf}}_{t}}^{\text{BAS}}& =& {\underline{f}}_{\text{surf}}^{\text{BAS}}-{\underline{f}}_{{\text{surf}}_{n}}^{\text{BAS}}\\ & =& {\underline{f}}_{\text{surf}}^{\text{BAS}}-({\underline{f}}_{{\text{surf}}_{n}}^{\text{BAS}}\cdot \underline{n})\underline{n}\end{array}\end{split}\]

Concernant les efforts de contact, il est absolument indispensable de construire la solution manufacturée de manière à ce qu’ils vérifient les équations du contact [bib2], à savoir :

(4817)#\[\begin{split}\begin{array}{ccc}\text{gap}(\underline{U})& \text{}\ge \text{}& 0\\ p& \text{}\le \text{}& 0\\ p\cdot \text{gap}(\underline{U})& =& 0\end{array}\end{split}\]

Cette vérification se fait après avoir calculé de façon analytique la pression et le saut de déplacement associés à la solution manufacturée, en général avec un outil de calcul formel (en l’occurence, il s’agit du module Python sympy ). On doit alors les visualiser, afin de vérifier rétrospectivement que la solution que l’on a contruite vérifie bien (5010). Dans le cas de ce test, nous avons représenté pression et saut de déplacement analytiques en figures Fig. 663 et Fig. 664. On remarque qu’ils vérifient \(p<0\) et \(\text{gap}(\underline{U})=0\) , ce qui est caractéristique d’une surface entièrement contactante, et conforme à (5010).

../../../../_images/100000000000029A000002888BE59863C946E150.png

Fig. 663 Validité de la solution manufacturée : pression p#

../../../../_images/10000000000002640000028D95A7401FA089DB08.png

Fig. 664 Validité de la solution manufacturée : jeu Gap(U)#

Grandeurs et résultats de référence#

La valeur de l’écart entre solutions analytique et calculée sur le maillage : \(\sum^{\text{noeuds}n}\mid {\underline{U}}_{\text{n}}^{\text{calc}}-{\underline{U}}_{\text{n}}^{\text{ref}}\mid\) et \(\sum^{\text{noeuds}n}\mid {p}_{\text{n}}^{\text{calc}}-{p}_{\text{n}}^{\text{ref}}\mid\) .

Dans le cas des modélisations qui réalisent une analyse de convergence avec la finesse du maillage, la vitesse de convergence avec la finesse du maillage de la solution calculée vers la solution analytique en norme \({L}_{2}\) :

  • le plus grand réel \({\alpha}_{U}>0\) tel que \({\parallel {\underline{U}}^{\text{calc}}-{\underline{U}}^{\text{ref}}\parallel }_{0,\Omega }<{C}_{U}\times {h}^{{\alpha}_{U}}\) où \({C}_{U}\) est indépendant de \(h\) pour le déplacement;

  • le plus grand réel \({\alpha}_{p}>0\) tel que \({\parallel {p}^{\text{calc}}-{p}^{\text{ref}}\parallel }_{0,{\Gamma}_{C}}<{C}_{p}\times {h}^{{\alpha}_{p}}\) où \({C}_{p}\) est indépendant de \(h\) pour la pression de contact.

Incertitudes sur la solution#

Aucune

Références bibliographiques#

[bib1] (1,2)

Document U2.08.08, Utilisation de la Méthode des Solutions Manufacturées pour la validation logicielle, Documentation U2 de Code_Aster

[bib2] (1,2)

Document R5.03.50, Formulation discrète du contact-frottement, Documentation R de Code_Aster

Modélisation A#

Caractéristiques de la modélisation#

On utilise une modélisation 3D et la méthode CONTINUE de traitement du contact.

Caractéristiques du maillage#

Le maillage contient 1 élément de type QUAD8, 768 éléments de type TRIA6 et 3072 éléments de type TETRA10.

Grandeurs testées et résultats#

On teste la somme des valeurs absolues de l’écart entre la solution calculée et la solution analytique.

Identification

Type de référence

Valeur de référence

\(\sum^{\text{noeuds}n}\mid {\underline{U}}_{\text{n}}^{\text{calc}}-{\underline{U}}_{\text{n}}^{\text{ref}}\mid\)

“NON_REGRESSION”

0.0410211809958

Modélisation B#

Caractéristiques de la modélisation#

On utilise une modélisation 3D et la méthode CONTINUE de traitement du contact.

Caractéristiques du maillage#

On réalise une étude de convergence avec la finesse du maillage de la solution calculée vers la solution analytique. Une suite de maillages obtenus par raffinement uniforme à l’aide de la commande MACR_ADAP_MAIL est utilisée :

  • maillage 0 : 1 élément de type QUAD8, 12 éléments de type TRIA6 et 6 éléments de type TETRA10

  • maillage 1 : 1 élément de type QUAD8, 48 éléments de type TRIA6 et 48 éléments de type TETRA10

  • maillage 2 : 1 élément de type QUAD8, 192 éléments de type TRIA6 et 384 éléments de type TETRA10

  • maillage 3 : 1 élément de type QUAD8, 768 éléments de type TRIA6 et 3072 éléments de type TETRA10

Grandeurs testées et résultats#

On teste la vitesse de convergence avec la finesse du maillage de la solution calculée vers la solution analytique en norme \({L}_{2}\) :

  • le plus grand réel \({\alpha}_{U}>0\) tel que \({\Vert {\underline{U}}^{\text{calc}}-{\underline{U}}^{\text{ref}}\Vert }_{0,\Omega }<{C}_{U}\times {h}^{{\alpha}_{U}}\) où \({C}_{U}\) est indépendant de \(h\) pour le déplacement;

  • le plus grand réel \({\alpha}_{p}>0\) tel que \({\Vert {p}^{\text{calc}}-{p}^{\text{ref}}\Vert }_{0,{\Gamma}_{C}}<{C}_{p}\times {h}^{{\alpha}_{p}}\) où \({C}_{p}\) est indépendant de \(h\) pour la pression de contact.

On teste aussi la somme des valeurs absolues de l’écart entre la solution calculée et la solution analytique pour le déplacement.

Identification

Type de référence

Valeur de référence

\(\sum^{\text{noeuds}n}|{\underline{U}}_{\text{n}}^{\text{calc}}-{\underline{U}}_{\text{n}}^{\text{ref}}|\)

“NON_REGRESSION”

4.13935026178E-05

\({\alpha}_{U}\)

“ANALYTIQUE”

3.0

\({\alpha}_{p}\)

“NON_REGRESSION”

2.534025066720

\({\alpha}_{p}\)

“ANALYTIQUE”

2.5

Modélisation C#

Caractéristiques de la modélisation#

On utilise une modélisation 3D et la méthode DISCRETE du gradient conjugué projeté (GCP) de traitement du contact.

Caractéristiques du maillage#

On réalise une étude de convergence avec la finesse du maillage de la solution calculée vers la solution analytique. Une suite de maillages obtenus par raffinement uniforme à l’aide de la commande MACR_ADAP_MAIL est utilisée :

  • maillage 0 : 14 éléments de type TRIA6 et 6 éléments de type TETRA10

  • maillage 1 : 50 éléments de type TRIA6 et 48 éléments de type TETRA10

  • maillage 2 : 194 éléments de type TRIA6 et 384 éléments de type TETRA10

  • maillage 3 : 770 éléments de type TRIA6 et 3072 éléments de type TETRA10

On note que, par rapport aux modélisations A et B, la base est maillée avec 2 TRIA6 au lieu d’un QUAD8. En effet, les méthodes discrètes ne sont pas adaptées à l’usage de ces éléments (voir [bib2]).

Grandeurs testées et résultats#

On teste la vitesse de convergence avec la finesse du maillage de la solution calculée vers la solution analytique en norme \({L}_{2}\) :

  • le plus grand réel \({\alpha}_{U}>0\) tel que \({\Vert {\underline{U}}^{\text{calc}}-{\underline{U}}^{\text{ref}}\Vert }_{0,\Omega }<{C}_{U}\times {h}^{{\alpha}_{U}}\) où \({C}_{U}\) est indépendant de \(h\) pour le déplacement dans le domaine \(\Omega\) ;

  • le plus grand réel \({\alpha}_{s}>0\) tel que \({\mathrm{\parallel }{U}^{\text{calc}}-{U}^{\text{ref}}\mathrm{\parallel }}_{0,{\Gamma}_{C}}<{C}_{s}\times {h}^{{\alpha}_{s}}\) où \({C}_{s}\) est indépendant de \(h\) pour le déplacement sur la surface \({\Gamma}_{C}\) .

On teste aussi la somme des valeurs absolues de l’écart entre la solution calculée et la solution analytique pour le déplacement.

Identification

Type de référence

Valeur de référence

\(\sum^{\text{noeuds}n}\mid {\underline{U}}_{\text{n}}^{\text{calc}}-{\underline{U}}_{\text{n}}^{\text{ref}}\mid\)

“NON_REGRESSION”

4.25881911029E-05

\({\alpha}_{U}\)

“ANALYTIQUE”

3.0

\({\alpha}_{p}\)

“ANALYTIQUE”

3.5

Modélisation D#

Caractéristiques de la modélisation#

On utilise une modélisation 3D et la méthode DISCRETE du gradient conjugué projeté (GCP) de traitement du contact.

Caractéristiques du maillage#

On réalise une étude de convergence avec la finesse du maillage de la solution calculée vers la solution analytique. Une suite de maillages obtenus par raffinement uniforme à l’aide de la commande MACR_ADAP_MAIL est utilisée :

  • maillage 0 : 2 éléments de type TRIA6, 6 éléments de type QUAD9 et 1 élément de type HEXA27

  • maillage 1 : 2 éléments de type TRIA6, 24 éléments de type QUAD9 et 8 éléments de type HEXA27

  • maillage 2 : 2 éléments de type TRIA6, 96 éléments de type QUAD9 et 64 éléments de type HEXA27

  • maillage 3 : 2 éléments de type TRIA6, 384 éléments de type QUAD9 et 512 éléments de type HEXA27

On note que, par rapport aux modélisations A et B, la base est maillée avec 2 TRIA6 au lieu d’un QUAD8. En effet, les méthodes discrètes ne sont pas adaptées à l’usage de ces éléments (voir r5.03.50).

Grandeurs testées et résultats#

On teste la somme des valeurs absolues de l’écart entre la solution calculée et la solution analytique pour le déplacement. On ne teste pas le taux de convergence car l’élément HEXA27 fournit la solution exacte avec un unique élément.

Identification

Type de référence

Valeur de référence

\(\sum^{\text{noeuds}n}\mid {\underline{U}}_{\text{n}}^{\text{calc}}-{\underline{U}}_{\text{n}}^{\text{ref}}\mid\)

“NON_REGRESSION »

0

Modélisation E#

Caractéristiques de la modélisation#

On utilise une modélisation 3D avec le type d’élément TETRA4 et la méthode LAC de traitement du contact.

Caractéristiques du maillage#

On réalise une étude de convergence avec la finesse du maillage de la solution calculée vers la solution analytique. Une suite de maillages obtenus par raffinement uniforme à l’aide de la commande MACR_ADAP_MAIL est utilisée :

  • maillage 0 : 16 éléments de type TRIA3 et 10 éléments de type TETRA4

  • maillage 1 : 64 éléments de type TRIA3 et 64 éléments de type TETRA4

  • maillage 2 : 256 éléments de type TRIA3 et 448 éléments de type TETRA4

  • maillage 3 : 1024 éléments de type TRIA3 et 3328 éléments de type TETRA4

  • maillage 4 : 4096 éléments de type TRIA3 et 25600 éléments de type TETRA4

La surface \(\text{MAITRE}\) courbée est maillée avec un unique QUAD8.

Grandeurs testées et résultats#

On teste la vitesse de convergence avec la finesse du maillage de la solution calculée vers la solution analytique en norme \({L}_{2}\) :

  • le plus grand réel \({\alpha}_{U}>0\) tel que \({\mathrm{\parallel }{\underline{U}}^{\text{calc}}-{\underline{U}}^{\text{ref}}\mathrm{\parallel }}_{0,\Omega }<{C}_{U}\times {h}^{{\alpha}_{U}}\) où \({C}_{U}\) est indépendant de \(h\) pour le déplacement;

  • le plus grand réel \({\alpha}_{p}>0\) tel que \({\Vert {p}^{\text{calc}}-{p}^{\text{ref}}\Vert }_{0,{\Gamma}_{C}}<{C}_{p}\times {h}^{{\alpha}_{p}}\) où \({C}_{p}\) est indépendant de \(h\) pour la pression de contact.

Identification

Type de référence

Valeur de référence

\({\alpha}_{U}\)

“ANALYTIQUE”

2.0

\({\alpha}_{p}\)

“ANALYTIQUE”

0.5

Modélisation F#

Caractéristiques de la modélisation#

On utilise une modélisation 3D avec le type d’élément TETRA10 et la méthode LAC de traitement du contact.

Caractéristiques du maillage#

On réalise une étude de convergence avec la finesse du maillage de la solution calculée vers la solution analytique. Une suite de maillages obtenus par raffinement uniforme à l’aide de la commande MACR_ADAP_MAIL est utilisée :

  • maillage 0 : 16 éléments de type TRIA6 et 10 éléments de type TETRA10

  • maillage 1 : 64 éléments de type TRIA6 et 64 éléments de type TETRA10

  • maillage 2 : 256 éléments de type TRIA6 et 448 éléments de type TETRA10

  • maillage 3 : 1024 éléments de type TRIA6 et 3328 éléments de type TETRA10

La surface \(\text{MAITRE}\) courbée est maillée avec un unique QUAD8.

Grandeurs testées et résultats#

On teste la vitesse de convergence avec la finesse du maillage de la solution calculée vers la solution analytique en norme \({L}_{2}\) :

  • le plus grand réel \({\alpha}_{U}>0\) tel que \({\mathrm{\parallel }{\underline{U}}^{\text{calc}}-{\underline{U}}^{\text{ref}}\mathrm{\parallel }}_{0,\Omega }<{C}_{U}\times {h}^{{\alpha}_{U}}\) où \({C}_{U}\) est indépendant de \(h\) pour le déplacement;

  • le plus grand réel \({\alpha}_{p}>0\) tel que \({\Vert {p}^{\text{calc}}-{p}^{\text{ref}}\Vert }_{0,{\Gamma}_{C}}<{C}_{p}\times {h}^{{\alpha}_{p}}\) où \({C}_{p}\) est indépendant de \(h\) pour la pression de contact.

Identification

Type de référence

Valeur de référence

\({\alpha}_{U}\)

“ANALYTIQUE”

2.5

\({\alpha}_{p}\)

“ANALYTIQUE”

1.0

Modélisation G#

Caractéristiques de la modélisation#

On utilise une modélisation 3D avec le type d’élément HEXA8 et la méthode LAC de traitement du contact.

Caractéristiques du maillage#

On réalise une étude de convergence avec la finesse du maillage de la solution calculée vers la solution analytique. Une suite de maillages obtenus par raffinement uniforme à l’aide de la commande MACR_ADAP_MAIL est utilisée :

  • maillage 0 : 10 éléments de type QUAD4 et 6 éléments de type HEXA8

  • maillage 1 : 40 éléments de type QUAD4 et 28 éléments de type HEXA8

  • maillage 2 : 160 éléments de type QUAD4 et 144 éléments de type HEXA8

  • maillage 3 : 640 éléments de type QUAD4 et 832 éléments de type HEXA8

  • maillage 4 : 2560 éléments de type QUAD4 et 5376 éléments de type HEXA8

On note que, par rapport aux modélisations E et F, la surface \(\text{MAITRE}\) courbée est maillée avec 2 TRIA6 au lieu d’un unique QUAD8.

Grandeurs testées et résultats#

On teste la vitesse de convergence avec la finesse du maillage de la solution calculée vers la solution analytique en norme \({L}_{2}\) :

  • le plus grand réel \({\alpha}_{U}>0\) tel que \({\mathrm{\parallel }{\underline{U}}^{\text{calc}}-{\underline{U}}^{\text{ref}}\mathrm{\parallel }}_{0,\Omega }<{C}_{U}\times {h}^{{\alpha}_{U}}\) où \({C}_{U}\) est indépendant de \(h\) pour le déplacement;

  • le plus grand réel \({\alpha}_{p}>0\) tel que \({\Vert {p}^{\text{calc}}-{p}^{\text{ref}}\Vert }_{0,{\Gamma}_{C}}<{C}_{p}\times {h}^{{\alpha}_{p}}\) où \({C}_{p}\) est indépendant de \(h\) pour la pression de contact.

Identification

Type de référence

Valeur de référence

\({\alpha}_{U}\)

“ANALYTIQUE”

2.0

\({\alpha}_{p}\)

“ANALYTIQUE”

0.5

Modélisation H#

Caractéristiques de la modélisation#

On utilise une modélisation 3D avec le type d’élément HEXA20 et la méthode LAC de traitement du contact.

Caractéristiques du maillage#

On réalise une étude de convergence avec la finesse du maillage de la solution calculée vers la solution analytique. Une suite de maillages obtenus par raffinement uniforme à l’aide de la commande MACR_ADAP_MAIL est utilisée :

  • maillage 0 : 10 éléments de type QUAD8 et 6 éléments de type HEXA20

  • maillage 1 : 40 éléments de type QUAD8 et 28 éléments de type HEXA20

  • maillage 2 : 160 éléments de type QUAD8 et 144 éléments de type HEXA20

  • maillage 3 : 640 éléments de type QUAD8 et 832 éléments de type HEXA20

On note que, par rapport aux modélisations E et F, la surface \(\text{MAITRE}\) courbée est maillée avec 2 TRIA6 au lieu d’un unique QUAD8.

Grandeurs testées et résultats#

On teste la vitesse de convergence avec la finesse du maillage de la solution calculée vers la solution analytique en norme \({L}_{2}\) :

  • le plus grand réel \({\alpha}_{U}>0\) tel que \({\mathrm{\parallel }{\underline{U}}^{\text{calc}}-{\underline{U}}^{\text{ref}}\mathrm{\parallel }}_{0,\Omega }<{C}_{U}\times {h}^{{\alpha}_{U}}\) où \({C}_{U}\) est indépendant de \(h\) pour le déplacement;

  • le plus grand réel \({\alpha}_{p}>0\) tel que \({\Vert {p}^{\text{calc}}-{p}^{\text{ref}}\Vert }_{0,{\Gamma}_{C}}<{C}_{p}\times {h}^{{\alpha}_{p}}\) où \({C}_{p}\) est indépendant de \(h\) pour la pression de contact.

Identification

Type de référence

Valeur de référence

\({\alpha}_{U}\)

“ANALYTIQUE”

2.5

\({\alpha}_{p}\)

“ANALYTIQUE”

1.0

Modélisation I#

Caractéristiques de la modélisation#

On utilise une modélisation 3D avec le type d’élément HEXA27 et la méthode LAC de traitement du contact.

Caractéristiques du maillage#

On réalise une étude de convergence avec la finesse du maillage de la solution calculée vers la solution analytique. Une suite de maillages obtenus par raffinement uniforme à l’aide de la commande MACR_ADAP_MAIL est utilisée :

  • maillage 0 : 6 éléments de type QUAD9 et 1 éléments de type HEXA27

  • maillage 1 : 24 éléments de type QUAD9 et 8 éléments de type HEXA27

  • maillage 2 : 96 éléments de type QUAD9 et 64 éléments de type HEXA27

  • maillage 3 : 384 éléments de type QUAD9 et 512 éléments de type HEXA27

On note que, par rapport aux modélisations E et F, la surface \(\text{MAITRE}\) courbée est maillée avec 2 TRIA6 au lieu d’un unique QUAD8.

Grandeurs testées et résultats#

On teste la vitesse de convergence avec la finesse du maillage de la solution calculée vers la solution analytique en norme \({L}_{2}\) :

  • le plus grand réel \({\alpha}_{U}>0\) tel que \({\mathrm{\parallel }{\underline{U}}^{\text{calc}}-{\underline{U}}^{\text{ref}}\mathrm{\parallel }}_{0,\Omega }<{C}_{U}\times {h}^{{\alpha}_{U}}\) où \({C}_{U}\) est indépendant de \(h\) pour le déplacement;

  • le plus grand réel \({\alpha}_{p}>0\) tel que \({\Vert {p}^{\text{calc}}-{p}^{\text{ref}}\Vert }_{0,{\Gamma}_{C}}<{C}_{p}\times {h}^{{\alpha}_{p}}\) où \({C}_{p}\) est indépendant de \(h\) pour la pression de contact.

Identification

Type de référence

Valeur de référence

\({\alpha}_{U}\)

“ANALYTIQUE”

2.5

\({\alpha}_{p}\)

“ANALYTIQUE”

1.0

Modélisation J#

Caractéristiques de la modélisation#

On utilise une modélisation 3D avec le type d’élément PYRAM5 et la méthode LAC de traitement du contact.

Caractéristiques du maillage#

On réalise une étude de convergence avec la finesse du maillage de la solution calculée vers la solution analytique. Une suite de maillages obtenus par raffinement uniforme à l’aide de la commande MACR_ADAP_MAIL est utilisée comme dans la modélisation G pour les éléments HEXA8.

Cependant, en utilisant l’option DECOUPE_HEXA= ‘PYRA’ dans la commande CREA_MAILLAGE, on peut traiter les mailles esclaves avec le type d’élément PYRAM5 pour le traitement du contactpar la méthode LAC :

  • maillage 0 : 4 éléments de type TRIA3, 5 éléments de type QUAD4, 5 éléments de type PYRAM5 et 0 élément de type HEXA8

  • maillage 1 : 16 éléments de type TRIA3, 20 éléments de type QUAD4, 20 éléments de type PYRAM5 et 4 éléments de type HEXA8

  • maillage 2 : 64 éléments de type TRIA3, 80 éléments de type QUAD4, 80 éléments de type PYRAM5 et 48 éléments de type HEXA8

  • maillage 3 : 256 éléments de type TRIA3, 320 éléments de type QUAD4, 320 éléments de type PYRAM5 et 448 éléments de type HEXA8

  • maillage 4 : 1024 éléments de type TRIA3, 1280 éléments de type QUAD4, 1280 éléments de type PYRAM5 et 3840 éléments de type HEXA8

On note que, par rapport aux modélisations E et F, la surface \(\text{MAITRE}\) courbée est maillée avec 2 TRIA6 au lieu d’un unique QUAD8.

Grandeurs testées et résultats#

On teste la vitesse de convergence avec la finesse du maillage de la solution calculée vers la solution analytique en norme \({L}_{2}\) :

  • le plus grand réel \({\alpha}_{U}>0\) tel que \({\mathrm{\parallel }{\underline{U}}^{\text{calc}}-{\underline{U}}^{\text{ref}}\mathrm{\parallel }}_{0,\Omega }<{C}_{U}\times {h}^{{\alpha}_{U}}\) où \({C}_{U}\) est indépendant de \(h\) pour le déplacement;

  • le plus grand réel \({\alpha}_{p}>0\) tel que \({\Vert {p}^{\text{calc}}-{p}^{\text{ref}}\Vert }_{0,{\Gamma}_{C}}<{C}_{p}\times {h}^{{\alpha}_{p}}\) où \({C}_{p}\) est indépendant de \(h\) pour la pression de contact.

Identification

Type de référence

Valeur de référence

\({\alpha}_{U}\)

“ANALYTIQUE”

2.0

\({\alpha}_{p}\)

“ANALYTIQUE”

0.5

Modélisation K#

Caractéristiques de la modélisation#

On utilise une modélisation 3D avec le type d’élément PYRAM13 et la méthode LAC de traitement du contact.

Caractéristiques du maillage#

On réalise une étude de convergence avec la finesse du maillage de la solution calculée vers la solution analytique. Une suite de maillages obtenus par raffinement uniforme à l’aide de la commande MACR_ADAP_MAIL est utilisée comme dans la modélisation G pour les éléments HEXA20.

Cependant, en utilisant l’option DECOUPE_HEXA= ‘PYRA’ dans la commande CREA_MAILLAGE, on peut traiter les mailles esclaves avec le type d’élément PYRAM13 pour le traitement du contactpar la méthode LAC :

  • maillage 0 : 4 éléments de type TRIA6, 5 éléments de type QUAD8, 5 éléments de type PYRAM13 et 0 élément de type HEXA20

  • maillage 1 : 16 éléments de type TRIA6, 20 éléments de type QUAD8, 20 éléments de type PYRAM13 et 4 éléments de type HEXA20

  • maillage 2 : 64 éléments de type TRIA6, 80 éléments de type QUAD8, 80 éléments de type PYRAM13 et 48 éléments de type HEXA20

  • maillage 3 : 256 éléments de type TRIA6, 320 éléments de type QUAD8, 320 éléments de type PYRAM13 et 448 éléments de type HEXA20

On note que, par rapport aux modélisations E et F, la surface \(\text{MAITRE}\) courbée est maillée avec 2 TRIA6 au lieu d’un unique QUAD8.

Grandeurs testées et résultats#

On teste la vitesse de convergence avec la finesse du maillage de la solution calculée vers la solution analytique en norme \({L}_{2}\) :

  • le plus grand réel \({\alpha}_{U}>0\) tel que \({\mathrm{\parallel }{\underline{U}}^{\text{calc}}-{\underline{U}}^{\text{ref}}\mathrm{\parallel }}_{0,\Omega }<{C}_{U}\times {h}^{{\alpha}_{U}}\) où \({C}_{U}\) est indépendant de \(h\) pour le déplacement;

  • le plus grand réel \({\alpha}_{p}>0\) tel que \({\Vert {p}^{\text{calc}}-{p}^{\text{ref}}\Vert }_{0,{\Gamma}_{C}}<{C}_{p}\times {h}^{{\alpha}_{p}}\) où \({C}_{p}\) est indépendant de \(h\) pour la pression de contact.

Identification

Type de référence

Valeur de référence

\({\alpha}_{U}\)

“ANALYTIQUE”

2.5

\({\alpha}_{p}\)

“ANALYTIQUE”

1.0

Modélisation L#

Caractéristiques de la modélisation#

On utilise une modélisation 3D avec le type d’élément PENTA6 et la méthode LAC de traitement du contact.

Caractéristiques du maillage#

On réalise une étude de convergence avec la finesse du maillage de la solution calculée vers la solution analytique. Une suite de maillages obtenus par raffinement uniforme à l’aide de la commande MACR_ADAP_MAIL est utilisée pour le type d’élément PENTA6.

Cependant, en utilisant l’option DECOUPE_HEXA= ‘PYRA’ dans la commande CREA_MAILLAGE, les mailles esclaves (PENTA6) sont découpés en TETRA4 et PYRAM5 pour le traitement du contactpar la méthode LAC :

  • maillage 0 : 8 éléments de type TRIA3, 4 éléments de type QUAD4, 2 éléments de type TETRA4, 6 éléments de type PYRAM5 et 0 élément de type PENTA6

  • maillage 1 : 32 éléments de type TRIA3, 16 éléments de type QUAD4, 8 éléments de type TETRA4, 24 éléments de type PYRAM5 et 8 éléments de type PENTA6

  • maillage 2 : 128 éléments de type TRIA3, 64 éléments de type QUAD4, 32 éléments de type TETRA4, 96 éléments de type PYRAM5 et 96 éléments de type PENTA6

  • maillage 3 : 512 éléments de type TRIA3, 256 éléments de type QUAD4, 128 éléments de type TETRA4, 384 éléments de type PYRAM5 et 896 éléments de type PENTA6

  • maillage 4 : 2048 éléments de type TRIA3, 1024 éléments de type QUAD4, 512 éléments de type TETRA4, 1536 éléments de type PYRAM5 et 7680 éléments de type PENTA6

On note que, par rapport aux modélisations E et F, la surface \(\text{MAITRE}\) courbée est maillée avec 2 TRIA6 au lieu d’un unique QUAD8.

Grandeurs testées et résultats#

On teste la vitesse de convergence avec la finesse du maillage de la solution calculée vers la solution analytique en norme \({L}_{2}\) :

  • le plus grand réel \({\alpha}_{U}>0\) tel que \({\mathrm{\parallel }{\underline{U}}^{\text{calc}}-{\underline{U}}^{\text{ref}}\mathrm{\parallel }}_{0,\Omega }<{C}_{U}\times {h}^{{\alpha}_{U}}\) où \({C}_{U}\) est indépendant de \(h\) pour le déplacement;

  • le plus grand réel \({\alpha}_{p}>0\) tel que \({\Vert {p}^{\text{calc}}-{p}^{\text{ref}}\Vert }_{0,{\Gamma}_{C}}<{C}_{p}\times {h}^{{\alpha}_{p}}\) où \({C}_{p}\) est indépendant de \(h\) pour la pression de contact.

Identification

Type de référence

Valeur de référence

\({\alpha}_{U}\)

“ANALYTIQUE”

2.0

\({\alpha}_{p}\)

“ANALYTIQUE”

0.5

Modélisation M#

Caractéristiques de la modélisation#

On utilise une modélisation 3D avec le type d’élément PENTA15 et la méthode LAC de traitement du contact.

Caractéristiques du maillage#

On réalise une étude de convergence avec la finesse du maillage de la solution calculée vers la solution analytique. Une suite de maillages obtenus par raffinement uniforme à l’aide de la commande MACR_ADAP_MAIL est utilisée pour le type d’élément PENTA15.

Cependant, en utilisant l’option DECOUPE_HEXA= ‘PYRA’ dans la commande CREA_MAILLAGE, les mailles esclaves (PENTA15) sont découpés en TETRA10 et PYRAM13 pour le traitement du contactpar la méthode LAC :

  • maillage 0 : 8 éléments de type TRIA6, 4 éléments de type QUAD8, 2 éléments de type TETRA10, 6 éléments de type PYRAM13 et 0 élément de type PENTA15

  • maillage 1 : 32 éléments de type TRIA6, 16 éléments de type QUAD8, 8 éléments de type TETRA10, 24 éléments de type PYRAM13 et 8 éléments de type PENTA15

  • maillage 2 : 128 éléments de type TRIA6, 64 éléments de type QUAD8, 32 éléments de type TETRA10, 96 éléments de type PYRAM13 et 96 éléments de type PENTA15

  • maillage 3 : 512 éléments de type TRIA6, 256 éléments de type QUAD8, 128 éléments de type TETRA10, 384 éléments de type PYRAM13 et 896 éléments de type PENTA15

On note que, par rapport aux modélisations E et F, la surface \(\text{MAITRE}\) courbée est maillée avec 2 TRIA6 au lieu d’un unique QUAD8.

Grandeurs testées et résultats#

On teste la vitesse de convergence avec la finesse du maillage de la solution calculée vers la solution analytique en norme \({L}_{2}\) :

  • le plus grand réel \({\alpha}_{U}>0\) tel que \({\mathrm{\parallel }{\underline{U}}^{\text{calc}}-{\underline{U}}^{\text{ref}}\mathrm{\parallel }}_{0,\Omega }<{C}_{U}\times {h}^{{\alpha}_{U}}\) où \({C}_{U}\) est indépendant de \(h\) pour le déplacement;

  • le plus grand réel \({\alpha}_{p}>0\) tel que \({\Vert {p}^{\text{calc}}-{p}^{\text{ref}}\Vert }_{0,{\Gamma}_{C}}<{C}_{p}\times {h}^{{\alpha}_{p}}\) où \({C}_{p}\) est indépendant de \(h\) pour la pression de contact.

Identification

Type de référence

Valeur de référence

\({\alpha}_{U}\)

“ANALYTIQUE”

2.5

\({\alpha}_{p}\)

“ANALYTIQUE”

1.0

Synthèse des résultats#

Les résultats sont en très bon accord avec la théorie.

Sur cette page