v6.03.154 SSNP154 – Benchmark NAFEMS de validation du contact 1 : cylinder roller contact#
Résumé:
Ce problème constitue le premier cas-test d’un benchmark NAFEMS de validation du contact-frottement. Les références du benchmark sont obtenues avec les codes Abaqus et MARC.
Ce test modélise un contact de Hertz avec et sans frottement entre un massif et un cylindre de raideurs différentes maillés avec des finesses différentes.
Une solution analytique existe pour le cas sans frottement.
Ce test comporte 6 modélisations permettant de tester:
les éléments linéaires et quadratiques,
les formulations de traitement du contact ‘CONTINUE’ sans et avec frottement.
les formulations de traitement du contact ‘DISCRETE’sans frottement
Solution de référence#
Méthode de calcul#
La solution de référence provient de résultats obtenus avec les codes Abaqus et MARC dans un benchmark NAFEMS de validation du contact-frottement [bib1].
Pour les modélisations sans frottement, une comparaison à la solution analytique de Hertz est réalisée.
Grandeurs et résultats de référence#
Déplacement vertical du point \(B\) (selon \(y\) ) (référence externe).
Pression de contact au point \(B\) (référence externe). La pression de contact relevée est celle extrapolée à partir des contraintes dans le volume.
Incertitudes sur la solution#
Importante (moyenne de codes).
Référence bibliographique#
[1] A. KONTER. «Advanced Finite Element Contact Benchmarks» . NAFEMS, 2006.
Modélisation A#
Caractéristiques de la modélisation#
La modélisation est D_PLAN, la formulation du contact est CONTINUE, le frottement est désactivé.
Caractéristiques du maillage#
Nombre de nœuds: 2625
Nombre de mailles et types: 2467 QUAD4 et 53 TRIA3.
Grandeurs testées et résultats#
Identification |
Type de référence |
Valeur de référence |
Tolérance |
\(\mathit{DY}\) au point \(B\) |
“SOURCE_EXTERNE” |
-1,3069040401077 |
2,6% |
\(\mathit{SIYY}\) au point \(B\) |
“SOURCE_EXTERNE” |
-3800,4949548283 |
0, 15% |
\(\text{PRES\_NOEU}\) au point \(B\) |
“SOURCE_EXTERNE” |
-3800,4949548283 |
0, 15% |
Remarques#
Les résultats obtenus en linéaire sans frottement avec la formulation continue sont en bon accord avec la solution analytique lorsque l’on relève la pression de contact extrapolée sur la surface maître.
On notera cependant la forte oscillation de la pression de contact lorsqu’elle est relevée à partir du degré de liberté LAGS_C spécifique à la formulation continue. L’amplitude des oscillations diminue lorsque l’on raffine le maillage. Les pressions extrapolées maître/esclave ont été obtenues avec la commande CALC_PRESSION.
Modélisation B#
Caractéristiques de la modélisation#
La modélisation est D_PLAN, la formulation du contact est CONTINUE, le frottement est désactivé.
Caractéristiques du maillage#
Nombre de nœuds: 7768
Nombre de mailles et types: 2467 QUAD8 et 53 TRIA6.
Grandeurs testées et résultats#
Identification |
Type de référence |
Valeur de référence |
Tolérance |
\(\mathit{DY}\) au point \(B\) |
“SOURCE_EXTERNE” |
-1,3069040401077 |
2,6% |
\(\mathit{SIYY}\) au point \(B\) |
“SOURCE_EXTERNE” |
-3800, 4949548282998 |
2,2% |
Remarques#
Les résultats obtenus en quadratique sans frottement avec la formulation continue sont en très bon accord avec la solution analytique quelle que soit la manière dont la pression est relevée.
On notera que les oscillations ont quasiment disparu par rapport à la modélisation A avec des éléments linéaires. Pour cela on a utilisé une intégration numérique des termes de contact de type Gauss (INTEGRATION=”GAUSS”). Les pressions extrapolées maître/esclave ont été obtenues avec la commande CALC_PRESSION.
Modélisation C#
Caractéristiques de la modélisation#
La modélisation est D_PLAN, la formulation du contact est CONTINUE, le frottement est pris en compte.
Caractéristiques du maillage#
Nombre de nœuds: 2625
Nombre de mailles et types: 2467 QUAD4 et 53 TRIA3.
Grandeurs testées et résultats#
Identification |
Type de référence |
Valeur de référence |
Tolérance |
\(\mathit{DY}\) au point \(B\) |
“SOURCE_EXTERNE” |
-1,3010795638714 |
2,5% |
\(\mathit{SIYY}\) au point \(B\) |
“SOURCE_EXTERNE” |
-3670,2561088662 |
1,6% |
\(\mathit{VX}\) au point \(B\) |
non-régression |
272.8972461688923 |
0,0001% |
\(\mathit{VY}\) au point \(B\) |
non-régression |
2.9322911927895867 |
0,0001% |
Remarques#
Les résultats obtenus en linéaire avec frottement avec la formulation continue sont en bon accord avec la solution analytique lorsque l’on relève la pression de contact extrapolée sur la surface maître.
On notera cependant la forte oscillation de la pression de contact lorsqu’elle est relevée à partir du degré de liberté LAGS_C spécifique à la formulation continue. L’amplitude des oscillations diminue lorsque l’on raffine le maillage. On utilise CALC_PRESSION pour les pressions extrapolées.
Les résultats obtenus par les codes commerciaux Abaqus et MARC dans le benchmark présentent eux aussi des oscillations pour le cas avec frottement.
Modélisation D#
Caractéristiques de la modélisation#
La modélisation est D_PLAN, la formulation du contact est CONTINUE, le frottement est pris en compte.
Caractéristiques du maillage#
Nombre de nœuds: 7768
Nombre de mailles et types: 2467 QUAD8 et 53 TRIA6.
Grandeurs testées et résultats#
Identification |
Type de référence |
Valeur de référence |
Tolérance |
\(\mathit{DY}\) au point \(B\) |
“SOURCE_EXTERNE” |
-1,3026356448389 |
2,5% |
\(\mathit{SIYY}\) au point \(B\) |
“SOURCE_EXTERNE” |
-3542,6878770218 |
2,9% |
\(\mathit{VX}\) au point \(B\) |
non-régression |
47.50792998023794 |
0,0001% |
\(\mathit{VY}\) au point \(B\) |
non-régression |
0.08685264239123569 |
0,0001% |
Remarques#
Les résultats obtenus en quadratique avec frottement avec la formulation continue sont en très bon accord avec ceux obtenus dans le benchmark (ici ceux du code MARC).
On notera que la pression de contact ainsi que la contrainte tangentielle relevées à partir des degrés de liberté LAGS_C et LAGS_F1 de la formulation continue collent parfaitement avec les grandeurs calculées par MARC. On utilise CALC_PRESSION pour les pressions extrapolées.
Comme c’était le cas sans frottement, les oscillations sont beaucoup moins importantes en quadratique, c’est d’ailleurs aussi le cas pour MARC.
Modélisation E#
Caractéristiques de la modélisation#
La modélisation est D_PLAN, la formulation du contact est DISCRETE, le frottement est désactivé. L’algorithme de contact utilisé est celui par défaut.
Caractéristiques du maillage#
Nombre de nœuds: 2625
Nombre de mailles et types: 2467 QUAD4 et 53 TRIA3.
Grandeurs testées et résultats#
Identification |
Type de référence |
Valeur de référence |
Tolérance |
\(\mathit{DY}\) au point \(B\) |
“SOURCE_EXTERNE” |
-1,3067596639633 |
2,6% |
\(\mathit{SIYY}\) au point \(B\) |
“SOURCE_EXTERNE” |
-3814,1811065942 |
0,18% |
Remarques#
Les résultats obtenus en linéaire sans frottement avec la formulation discrète sont identiques à ceux obtenus en formulation continue (modélisation A). On utilise CALC_PRESSION pour les pressions extrapolées.
En formulation discrète on est obligé de lier les deux nœuds du bloc et du cylindre en vis-à-vis initialement pour empêcher le mouvement de corps rigide suivant l’axe \(y\) .
Modélisation F#
Caractéristiques de la modélisation#
La modélisation est D_PLAN, la formulation du contact est DISCRETE, le frottement est désactivé.
L’algorithme de contact utilisé est celui par défaut.
Caractéristiques du maillage#
Nombre de nœuds: 7768
Nombre de mailles et types: 2467 QUAD8 et 53 TRIA6.
Grandeurs testées et résultats#
Identification |
Type de référence |
Valeur de référence |
Tolérance |
\(\mathit{DY}\) au point \(B\) |
“SOURCE_EXTERNE” |
-1,3080923611548 |
2,6% |
\(\mathit{SIYY}\) au point \(B\) |
“SOURCE_EXTERNE” |
-3502,7937782466 |
2,2% |
Remarques#
Les résultats obtenus en quadratique sans frottement avec la formulation discrète sont identiques à ceux obtenus en formulation continue (modélisation B). On utilise CALC_PRESSION pour les pressions extrapolées.
En formulation discrète on est obligé de lier les deux nœuds du bloc et du cylindre en vis-à-vis initialement pour empêcher le mouvement de corps rigide suivant l’axe \(y\) .
Synthèse des résultats#
Ce test permet de valider le contact-frottement par rapport à une solution analytique mais aussi par rapport à des références données par des codes de calcul commerciaux (Abaqus et MARC).
On observe un bon accord entre les résultats analytiques et ceux obtenus par Code_Aster ainsi qu’une bonne concordance avec les résultats de référence pour les calculs avec frottement.
On notera que les formulations continue et discrète donnent des résultats identiques avec cependant les restrictions suivantes:
le calcul en formulation continue est plus facile à mener car le blocage du mouvement de corps rigide vertical est réalisé automatiquement
le calcul avec frottement sur un maillage quadratique n’a pas pu être mené en formulation discrète
En quadratique, la pression de contact obtenue à partir du degré de liberté LAGS_C en formulation continue présente beaucoup moins d’oscillations qu’en linéaire.