v6.03.121 SSNP121 – Intégration des termes de contact en 2D et 3D#
Résumé:
Ce problème correspond à une analyse quasi-statique d’un problème de mécanique avec contact sans frottement. On s’intéresse tout particulièrement ici à l’intégration des termes de contact ( patch-test ). Il s’agit d’étudier deux blocs identiques maillés différemment et soumis à des déplacements imposés symétriques.
Ce test comporte six modélisations en \(\mathrm{2D}\) (éléments linéaires SEG2):
modélisation A: METHODE=‘CONTINUE’. La méthode d’intégration par sous-éléments proposée aussi par Bathe [bib1] est utilisée avec trois sous-éléments,
modélisation B: METHODE=‘CONTRAINTE’,
modélisation C: METHODE=‘PENALISATION’,
modélisation G: METHODE=‘CONTINUE’. Une variante de la modélisation Aoù l’on teste l’activation initiale du statut de contact.
modélisation L: METHODE=‘GCP’,
modélisation P: METHODE=‘CONTINUE’. Une variante de la modélisation Goù l’on teste l’algorithme de pénalisation du contact.
une modélisation en \(\mathrm{2D}\) (éléments linéaires SEG2 en vis-à-vis d’éléments quadratiques SEG3):
modélisation H: METHODE=‘CONTINUE’. Les surfaces de contact sont composées d’éléments SEG2 en vis-à-vis d’éléments SEG3,
six modélisations en \(\mathrm{3D}\) (éléments quadratiques):
modélisation D: METHODE=‘CONTINUE’. Les surfaces de contact sont composées d’éléments QUAD8 en vis-à-vis d’éléments QUAD8,
modélisation E: METHODE=‘CONTINUE’. Les surfaces de contact sont composées d’éléments TRIA6 en vis-à-vis d’éléments TRIA6,
modélisation F: METHODE=‘CONTINUE’. Les surfaces de contact sont composées d’éléments TRIA6 en vis-à-vis d’éléments QUAD8.
modélisation Q: METHODE=‘CONTINUE’. Les surfaces de contact sont composées d’éléments QUAD9 en vis-à-vis d’éléments QUAD9, les mailles \(\mathrm{3D}\) sont des HEXA27,
modélisation R: METHODE=‘CONTINUE’. Les surfaces de contact sont composées d’éléments TRIA6 en vis-à-vis d’éléments TRIA6, les mailles \(\mathrm{3D}\) sont des PENTA18,
et quatre modélisations en \(\mathrm{3D}\) (éléments linéaires en vis-à-vis d’éléments quadratiques):
modélisation I: METHODE=‘CONTINUE’. Les surfaces de contact sont composées d’éléments QUAD4 en vis-à-vis d’éléments QUAD8,
modélisation J: METHODE=‘CONTINUE’. Les surfaces de contact sont composées d’éléments TRIA3 en vis-à-vis d’éléments TRIA6,
modélisation K: METHODE=‘CONTINUE’. Les surfaces de contact sont composées d’éléments TRIA3 en vis-à-vis d’éléments QUAD8.
Solution de référence#
Méthode de calcul utilisée pour la solution de référence#
La solution de référence, analytique, peut être déduite d’un calcul très simple. Le coefficient de Poisson étant nul, la déformation uniaxiale s’écrit \(\varepsilon ={U}_{0}/\mathit{HD}=–0.1/2.\) La pression vaut alors \(E\varepsilon =–{10}^{5}\mathit{Pa}\) .
Résultats de référence#
La pression de contact est constante et égale à \(–{10}^{5}\mathit{Pa}\) sur toute la surface de contact. De même le déplacement vertical (selon \(y\) ) est constant sur la surface de contact et égal à \({U}_{0}/2=–0.05\text{m}\) .
Référence bibliographique#
EL-ABBASI and K.J. BATHE: « Stability and Patch Test Performance of Contact Discretizations and a New Solution Algorithm « , Computers & Structures, 79, 1473-1486, 2001
Modélisation A#
Caractéristiques de la modélisation#
On utilise une modélisation D_PLAN pour les éléments solides avec la méthode CONTINUE pour le traitement du CONTACT avec intégration de type ‘SIMPSON2’.
Caractéristiques du maillage#
Nombre de nœuds: |
313 |
Nombre de mailles et types: |
265 QUAD4et 132 SEG2 |
Les maillages des 2 surfaces de contact sont incompatibles. 12 éléments finis SEG2 sont disposés sur la surface de contact initiale de la plaque 1 et seulement 11 sur la surface de contact de l’autre surface.
En activant le mot-clé SIMPSON2 trois sous-éléments sont utilisés pour l’intégration des termes de contact.
Grandeurs testées et résultats#
Identification |
Valeur de référence |
Type de référence |
Tolérance |
LAGS_Cau noeud \(\mathit{N14}\) |
–1.E+5 |
“ANALYTIQUE” |
1 % |
LAGS_Cau point \(A\) |
–1.E+5 |
“ANALYTIQUE” |
1 % |
LAGS_Cau point \(B\) |
–1.E+5 |
“ANALYTIQUE” |
1 % |
DYau noeud \(\mathit{N14}\) |
-0.05 |
“ANALYTIQUE” |
1 % |
DYau point \(A\) |
-0.05 |
“ANALYTIQUE” |
1 % |
DYau point \(B\) |
-0.05 |
“ANALYTIQUE” |
1 % |
Modélisation B#
Caractéristiques de la modélisation#
On utilise une modélisation D_PLAN pour les éléments solides avec la méthode CONTRAINTE pour le traitement du CONTACT.
Caractéristiques du maillage#
On utilise le même maillage que pour la modélisation précédente.
Nombre de nœuds: |
313 |
Nombre de mailles et types: |
265 QUAD4et 132 SEG2 |
Grandeurs testées et résultats#
Identification |
Valeur de référence |
Type de référence |
Tolérance |
SIYYau noeud \(\mathit{N14}\) |
–1.E+5 |
“ANALYTIQUE” |
1 % |
SIYYau point \(A\) |
–1.E+5 |
“ANALYTIQUE” |
6.1% |
SIYYau point \(B\) |
–1.E+5 |
“ANALYTIQUE” |
6.1% |
DYau noeud \(\mathit{N14}\) |
-0.05 |
“ANALYTIQUE” |
1 % |
DYau point \(A\) |
-0.05 |
“ANALYTIQUE” |
1,67 % |
DYau point \(B\) |
-0.05 |
“ANALYTIQUE” |
1,67 % |
Remarques#
Les formulations discrètes ne vérifient pas le patch-test contrairement à la formulation continue, c’est pourquoi les tolérances sont plus lâches.
Modélisation C#
Caractéristiques de la modélisation#
On utilise une modélisation D_PLAN pour les éléments solides avec la méthode PENALISATION avec un coefficient de pénalisation de 1.E7 pour le traitement du contact.
Caractéristiques du maillage#
On utilise le même maillage que pour la modélisation précédente.
Nombre de nœuds: |
313 |
Nombre de mailles et types: |
265 QUAD4et 132 SEG2 |
Grandeurs testées et résultats#
Identification |
Valeur de référence |
Type de référence |
Tolérance |
SIYYau noeud \(\mathit{N14}\) |
–1.E+5 |
“ANALYTIQUE” |
2,13 % |
SIYYau point \(A\) |
–1.E+5 |
“ANALYTIQUE” |
7,5% |
SIYYau point \(B\) |
–1.E+5 |
“ANALYTIQUE” |
7,5% |
DYau noeud \(\mathit{N14}\) |
-0.05 |
“ANALYTIQUE” |
2,13 % |
DYau point \(A\) |
-0.05 |
“ANALYTIQUE” |
0,7 % |
DYau point \(B\) |
-0.05 |
“ANALYTIQUE” |
0,7 % |
Remarques#
Les formulations discrètes ne vérifient pas le patch-test contrairement à la formulation continue, c’est pourquoi les tolérances sont plus lâches.
Modélisation D#
Caractéristiques de la modélisation#
On utilise une modélisation 3D pour les éléments solides avec la méthode CONTINUE pour le traitement du CONTACT. Les surfaces de contact en vis-à-vis sont constituées d’éléments QUAD8.
Caractéristiques du maillage#
Nombre de nœuds: |
850 |
Nombre de mailles et types: |
128 HEXA20et 64 QUAD8 |
Grandeurs testées et résultats#
Identification |
Valeur de référence |
Type de référence |
Tolérance |
LAGS_Cau noeud \(\mathit{A21}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A22}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A23}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A24}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
MAX des LAGS_Csur la surface |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
MIN des LAGS_Csur la surface |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A21}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A22}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A23}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A24}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
Modélisation E#
Caractéristiques de la modélisation#
On utilise une modélisation 3D pour les éléments solides avec la méthode CONTINUE pour le traitement du CONTACT. Les surfaces de contact en vis-à-vis sont constituées d’éléments TRIA6.
Caractéristiques du maillage#
Nombre de nœuds: |
1010 |
Nombre de mailles et types: |
256 PENTA15et 128 TRIA6 |
Grandeurs testées et résultats#
Identification |
Valeur de référence |
Type de référence |
Tolérance |
LAGS_Cau noeud \(\mathit{A21}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A22}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A23}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A24}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A21}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A22}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A23}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A24}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
MAX des LAGS_Csur la surface |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
MIN des LAGS_Csur la surface |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
Modélisation F#
Caractéristiques de la modélisation#
On utilise une modélisation 3D pour les éléments solides avec la méthode CONTINUE pour le traitement du CONTACT. Les surfaces de contact en vis-à-vis sont constituées d’éléments TRIA6 en vis-à-vis de mailles QUAD8.
Caractéristiques du maillage#
Nombre de nœuds: |
930 |
Nombre de mailles et types: |
64 HEXA20, 128 PENTA15, 32 QUAD8et 64 TRIA6 |
Grandeurs testées et résultats#
Identification |
Valeur de référence |
Type de référence |
Tolérance |
LAGS_Cau noeud \(\mathit{A21}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A22}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A23}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A24}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A21}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A22}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A23}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A24}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
MAX des LAGS_Csur la surface |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
MIN des LAGS_Csur la surface |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
Remarques#
La surface de contact esclave est maillée en TRIA6 et la surface de contact maître est maillée en QUAD8. Une utilisation avec des éléments QUAD8 pour la surface esclave et des TRIA6 pour la surface maître ne satisfait pas les conditions de compatibilité nécessaires à la bonne intégration des termes de contact [R5.03.52]. De manière générale, si l’on n’appréhende pas correctement cette notion de compatibilité, on conseille à l’utilisateur non averti d’utiliser les mêmes éléments pour le maillage des surfaces maître et esclave.
Modélisation G#
Caractéristiques de la modélisation#
On utilise une modélisation D_PLAN pour les éléments solides avec la méthode CONTINUE pour le traitement du CONTACT avec intégration de type SIMPSON2.
Caractéristiques du maillage#
Nombre de nœuds: |
313 |
Nombre de mailles et types: |
265 QUAD4et 132 SEG2 |
Grandeurs testées et résultats#
Identification |
Valeur de référence |
Type de référence |
Tolérance |
LAGS_Cau noeud \(\mathit{N14}\) |
–1.E+5 |
“ANALYTIQUE” |
1 % |
LAGS_Cau point \(A\) |
–1.E+5 |
“ANALYTIQUE” |
1 % |
LAGS_Cau point \(B\) |
–1.E+5 |
“ANALYTIQUE” |
1 % |
DYau noeud \(\mathit{N14}\) |
-0.05 |
“ANALYTIQUE” |
1 % |
DYau point \(A\) |
-0.05 |
“ANALYTIQUE” |
1 % |
DYau point \(B\) |
-0.05 |
“ANALYTIQUE” |
1 % |
Remarques#
Dans cette modélisation, le chargement est une pression imposée sur la plaque supérieure (et non un déplacement). Les 2 approches sont parfaitement équivalentes, cependant l’utilisation d’une pression impose ici de bloquer les mouvements de corps rigide (suivant \(\mathrm{DZ}\) ), cela est fait en déclarant un contact initial (CONTACT_INIT).
Modélisation H#
Caractéristiques de la modélisation#
On utilise une modélisation D_PLAN pour les éléments solides avec la méthode CONTINUE pour le traitement du CONTACT entre éléments mixtes linéaires et éléments quadratiques.
Caractéristiques du maillage#
Nombre de nœuds: |
721 |
Nombre de mailles et types: |
144 QUAD4, 121 QUAD8, 48 SEG2et 44 SEG3 |
Les maillages sont incompatibles. 12 éléments finis SEG2 sont disposés sur la surface de contact initiale de la plaque esclave et seulement 11 éléments finis SEG3 sur la surface de contact maître.
En activant le mot-clé NCOTES2, un schéma de type Newton-Cotes couplé à une technique de subdivision en sous-éléments a été utilisé pour l’intégration des termes de contact.
Grandeurs testées et résultats#
Identification |
Valeur de référence |
Type de référence |
Tolérance |
LAGS_Cau noeud \(\mathit{N14}\) |
–1.E+5 |
“ANALYTIQUE” |
0,5 % |
LAGS_Cau point \(A\) |
–1.E+5 |
“ANALYTIQUE” |
0,5 % |
LAGS_Cau point \(B\) |
–1.E+5 |
“ANALYTIQUE” |
0,5 % |
DYau noeud \(\mathit{N14}\) |
-0.05 |
“ANALYTIQUE” |
0,5 % |
DYau point \(A\) |
-0.05 |
“ANALYTIQUE” |
0,5 % |
DYau point \(B\) |
-0.05 |
“ANALYTIQUE” |
0,5 % |
Modélisation I#
Caractéristiques de la modélisation#
On utilise une modélisation 3D pour les éléments solides avec la méthode CONTINUE pour le traitement du CONTACT entre éléments mixtes linéaires et éléments quadratiques.
Caractéristiques du maillage#
Nombre de nœuds: |
227 |
Nombre de mailles et types: |
16 HEXA20, 25 HEXA8, 32 QUAD8et 50 QUAD4 |
Les maillages sont incompatibles. Les surfaces de contact sont constituées de 5 éléments QUAD4 en vis‑à-vis de 4 éléments QUAD8. Le schéma d’intégration utilisé est de type Newton-Cotes couplé à une technique de subdivision en sous-éléments.
Grandeurs testées et résultats#
Identification |
Valeur de référence |
Type de référence |
Tolérance |
LAGS_Cau noeud \(\mathit{A21}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A22}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A23}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A24}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
MAX des LAGS_Csur la surface |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
MIN des LAGS_Csur la surface |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A21}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A22}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A23}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A24}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
Modélisation J#
Caractéristiques de la modélisation#
On utilise une modélisation 3D pour les éléments solides avec la méthode CONTINUE pour le traitement du CONTACT entre éléments mixtes linéaires et éléments quadratiques.
Caractéristiques du maillage#
Nombre de nœuds: |
630 |
Nombre de mailles et types: |
128 PENTA6, 128 PENTA15, 64 TRIA3et 64 TRIA6 |
Les maillages sont compatibles. Les surfaces de contact sont constituées d’éléments TRIA3 en vis-à-vis d’éléments TRIA6. Le schéma d’intégration utilisé est de type Newton-Cotes.
Grandeurs testées et résultats#
Identification |
Valeur de référence |
Type de référence |
Tolérance |
LAGS_Cau noeud \(\mathit{A21}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A22}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A23}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A24}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A21}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A22}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A23}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A24}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
MAX des LAGS_Csur la surface |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
MIN des LAGS_Csur la surface |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
Modélisation K#
Caractéristiques de la modélisation#
On utilise une modélisation 3D pour les éléments solides avec la méthode CONTINUE pour le traitement du CONTACT pour les éléments mixtes linéaires/quadratiques .
Caractéristiques du maillage#
Nombre de nœuds: |
550 |
Nombre de mailles et types: |
64 HEXA20, 128 PENTA6, 32 QUAD8et 64 TRIA3 |
Les maillages sont compatibles. Les surfaces de contact sont constituées d’éléments TRIA3 en vis-à-vis d’éléments QUAD8. Le schéma d’intégration utilisé est de type Newton-Cotes.
Grandeurs testées et résultats#
Identification |
Valeur de référence |
Type de référence |
Tolérance |
LAGS_Cau noeud \(\mathit{A21}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A22}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A23}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A24}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A21}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A22}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A23}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A24}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
MAX des LAGS_Csur la surface |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
MIN des LAGS_Csur la surface |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
Remarques#
La surface de contact esclave est maillée en TRIA3 et la surface de contact maître est maillée en QUAD8. Une utilisation avec des éléments QUAD8 pour la surface esclave et des TRIA3 pour la surface maître ne satisfait pas les conditions de compatibilité nécessaires à la bonne intégration des termes de contact [R5.03.52].
Modélisation L#
Caractéristiques de la modélisation#
On utilise une modélisation D_PLAN pour les éléments solides avec la méthode GCP pour le traitement du CONTACT.
Caractéristiques du maillage#
On utilise le même maillage que pour la modélisation A.
Nombre de nœuds: |
313 |
Nombre de mailles et types: |
265 QUAD4et 132 SEG2 |
Grandeurs testées et résultats#
Identification |
Valeur de référence |
Type de référence |
Tolérance |
SIYYau noeud \(\mathit{N14}\) |
–1.E+5 |
“ANALYTIQUE” |
2,13 % |
SIYYau point \(A\) |
–1.E+5 |
“ANALYTIQUE” |
7,5% |
SIYYau point \(B\) |
–1.E+5 |
“ANALYTIQUE” |
7,5% |
DYau noeud \(\mathit{N14}\) |
-0.05 |
“ANALYTIQUE” |
2,13% |
DYau point \(A\) |
-0.05 |
“ANALYTIQUE” |
0,7 % |
DYau point \(B\) |
-0.05 |
“ANALYTIQUE” |
0,7 % |
Remarques#
Les formulations discrètes ne vérifient pas le patch-test contrairement à la formulation continue, c’est pourquoi les tolérances sont plus lâches.
Modélisation P#
Caractéristiques de la modélisation#
On utilise une modélisation D_PLAN pour les éléments solides avec la méthode CONTINUE pour le traitement du CONTACT avec intégration de type SIMPSON2. La formulation du contact est pénalisée, le coefficient de pénalisation est de 1.0E16.
Caractéristiques du maillage#
Nombre de nœuds: |
313 |
Nombre de mailles et types: |
265 QUAD4et 132 SEG2 |
Grandeurs testées et résultats#
Identification |
Valeur de référence |
Type de référence |
Tolérance |
LAGS_Cau noeud \(\mathit{N14}\) |
–1.E+5 |
“ANALYTIQUE” |
1 % |
LAGS_Cau point \(A\) |
–1.E+5 |
“ANALYTIQUE” |
1 % |
LAGS_Cau point \(B\) |
–1.E+5 |
“ANALYTIQUE” |
1 % |
DYau noeud \(\mathit{N14}\) |
-0.05 |
“ANALYTIQUE” |
1 % |
DYau point \(A\) |
-0.05 |
“ANALYTIQUE” |
1 % |
DYau point \(B\) |
-0.05 |
“ANALYTIQUE” |
1 % |
Remarques#
Dans cette modélisation, le chargement est une pression imposée sur la plaque supérieure (et non un déplacement). Les 2 approches sont parfaitement équivalentes, cependant l’utilisation d’une pression impose ici de bloquer les mouvements de corps rigide (suivant \(\mathrm{DZ}\) ), cela est fait en déclarant un contact initial (CONTACT_INIT).
Modélisation Q#
Caractéristiques de la modélisation#
On utilise une modélisation 3D pour les éléments solides avec la méthode CONTINUE pour le traitement du CONTACT. Les surfaces de contact en vis-à-vis sont constituées d’éléments QUAD9.
Caractéristiques du maillage#
Nombre de nœuds: |
850 |
Nombre de mailles et types: |
128 HEXA27et 64 QUAD9 |
Grandeurs testées et résultats#
Identification |
Valeur de référence |
Type de référence |
Tolérance |
LAGS_Cau noeud \(\mathit{A21}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A22}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A23}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A24}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
MAX des LAGS_Csur la surface |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
MIN des LAGS_Csur la surface |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A21}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A22}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A23}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A24}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
Remarques#
Dans cette modélisation, on utilise des mailles HEXA27 qui ont été transformées par l’opérateur CREA_MAILLAGE.
Modélisation R#
Caractéristiques de la modélisation#
On utilise une modélisation 3D pour les éléments solides avec la méthode CONTINUE pour le traitement du CONTACT. Les surfaces de contact en vis-à-vis sont constituées d’éléments TRIA6.
Caractéristiques du maillage#
Nombre de nœuds: |
1458 |
Nombre de mailles et types: |
256 PENTA18, 128 TRIA6 |
Grandeurs testées et résultats#
Identification |
Valeur de référence |
Type de référence |
Tolérance |
LAGS_Cau noeud \(\mathit{A21}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A22}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A23}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
LAGS_Cau point \(\mathit{A24}\) |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
MAX des LAGS_Csur la surface |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
MIN des LAGS_Csur la surface |
–1.E+5 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A21}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A22}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A23}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
DZau point \(\mathit{A24}\) |
-0.05 |
“ANALYTIQUE” |
0,1 % |
Remarques#
Dans cette modélisation, on utilise des mailles PENTA18 qui ont été transformées par l’opérateur CREA_MAILLAGE.
Modélisation S#
Caractéristiques de la modélisation#
On utilise une modélisation D_PLAN pour les éléments solides avec la méthode CONTINUE pour le traitement du CONTACT. On teste la méthode de PENALISATION avec le Newton généralisé.
Caractéristiques du maillage#
Le maillage est le même que celui de la modélisation A.
Grandeurs testées et résultats#
Identification |
Valeur de référence |
Type de référence |
Tolérance |
LAGS_Cau noeud \(\mathit{N14}\) |
–1.E+5 |
“ANALYTIQUE” |
1 % |
LAGS_Cau point \(A\) |
–1.E+5 |
“ANALYTIQUE” |
1 % |
LAGS_Cau point \(B\) |
–1.E+5 |
“ANALYTIQUE” |
1 % |
DYau noeud \(\mathit{N14}\) |
-0.05 |
“ANALYTIQUE” |
1 % |
DYau point \(A\) |
-0.05 |
“ANALYTIQUE” |
1 % |
DYau point \(B\) |
-0.05 |
“ANALYTIQUE” |
1 % |
Remarques#
Dans cette modélisation, on utilise des mailles PENTA18 qui ont été transformées par l’opérateur CREA_MAILLAGE.
Synthèse des résultats#
On cherche sur cet exemple très simple (connu sous le nom de patch-test ) à tester une nouvelle technique d’intégration des termes de contact fondée sur la subdivision par sous-éléments disponible uniquement pour la formulation ‘CONTINUE’. Cette technique a pour objectif d’atténuer l’amplitude d’oscillation de la pression de contact. Dans le cas étudié ici, la pression est constante sur toute la surface de contact (remarquons que le coefficient de Poisson est nul).
On constate ainsi qu’avec les formulations discrètes ‘CONTRAINTE’ et ‘PENALISATION’ la solution présente des oscillations non physiques de l’ordre de 6 à 7 %. En utilisant la méthode ‘CONTINUE’, avec un schéma d’intégration adapté les oscillations disparaissent presque totalement et les résultats obtenus sont très proches de la solution de référence (<0,5%).