v6.04.303 SSNV303 – Poutre encastrée non-linéaire#

Résumé

Ce test permet de valider la prise en compte des effets non-linéaires (cinématique et matériau) pour les éléments COQUE_SOLIDE. Les résultats sont comparés par rapport à une modélisation 3D très fine (cas-test ssnv302).

  • Modélisation A: cinématique PETIT en élasticité

  • Modélisation B: cinématique GDEF_LOG en élasticité

  • Modélisation C: plasticité de type VMIS_ISOT_LINE

  • Modélisation D: test du raccord entre modélisations avec LIAISON_MAILLE

Solution de référence#

Solution de référence#

La solution de référence est prise sur un calcul 3D quadratique fin (test ssnv302): tous les tests sont donc de type AUTRE_ASTER pour les déplacements et les contraintes.

En flexion, on teste:

  • Déplacement DZ au nœud I

  • Contraintes de von Mises SIEQ_NOEU, composante VMIS aux nœuds C, E et G

En traction, on teste:

  • Déplacement DX au nœud I

  • Contraintes de von Mises SIEQ_NOEU, composante VMIS aux nœuds C, E et G

Résultats de référence#

Élasticité petites déformations, en flexion (modélisation A), cas-test ssnv302a:

Lieu

NOM_CHAM

NOM_CMP

Valeur (test ssnv302)

Nœud I

DEPL

DZ

15.84 mm

Nœud C

SIEQ_NOEU

VMIS

180018 Pa

SIEF_NOEU

SIXX

-180032.Pa

Nœud E

SIEQ_NOEU

VMIS

120015 Pa

SIEF_NOEU

SIXX

-120000 Pa

Nœud G

SIEQ_NOEU

VMIS

60032 Pa

SIEF_NOEU

SIXX

-60000 Pa

Élasticité, grandes déformations, en flexion (modélisation B), cas-test ssnv302b:

Lieu

NOM_CHAM

NOM_CMP

Valeur (test ssnv302)

Nœud I

DEPL

DZ

6.984 mm

Nœud C

SIEQ_NOEU

VMIS

96137 Pa

Nœud E

SIEQ_NOEU

VMIS

53027Pa

Nœud G

SIEQ_NOEU

VMIS

22960Pa

Plasticité, petites déformations, en traction (modélisationC), cas-test ssnv302c:

Lieu

NOM_CHAM

NOM_CMP

Valeur (test ssnv302)

Nœud I

DEPL

DZ

0.035mm

Nœud C

SIEQ_NOEU

VMIS

10000092 Pa

SIEF_NOEU

SIXX

10000151 Pa

Nœud E

SIEQ_NOEU

VMIS

10000000Pa

SIEF_NOEU

SIXX

10000000Pa

Nœud G

SIEQ_NOEU

VMIS

10000000Pa

SIEF_NOEU

SIXX

10000000Pa

Élasticité petites déformations, en flexion (modélisation D) avec raccord 3D, cas-test ssnv303a:

Lieu

NOM_CHAM

NOM_CMP

Valeur (test ssnv302a)

Nœud I

DEPL

DZ

15.84 mm

Nœud E

SIEQ_NOEU

VMIS

120015 Pa

Nœud G

SIEQ_NOEU

VMIS

60032 Pa

Incertitude sur la solution#

Solutions numériques par rapport au 3D avec maillage overkill (tests ssnv302 dans validation)

Modélisation A#

Caractéristiques de la modélisation#

On utilise une modélisation COQUE_SOLIDE. On est en flexion, élasticité isotrope linéaire, petites déformations (ELAS / PETIT).

Caractéristiques du maillage#

Le maillage est le suivant:

  • Nombre de nœuds : 294

  • Nombre de mailles et types : 120 HEXA9

Résultats#

Les résultats sont comparés à ssnv302a (calcul 3D overkill), référence AUTRE_ASTER :

Lieu

NOM_CHAM

NOM_CMP

PRECISION

Nœud I

DEPL

DZ

0, 11%

Nœud C

SIEQ_NOEU

VMIS

0, 042%

SIEF_NOEU

SIXX

0,025%

Nœud E

SIEQ_NOEU

VMIS

0, 0038%

SIEF_NOEU

SIXX

0,020%

Nœud G

SIEQ_NOEU

VMIS

0, 015%

SIEF_NOEU

SIXX

0,0007%

Modélisation B#

Caractéristiques de la modélisation#

On utilise une modélisation COQUE_SOLIDE. On est en flexion, élasticité isotrope linéaire, grandes déformations (ELAS / GDEF_LOG).

Caractéristiques du maillage#

Le maillage est le suivant:

  • Nombre de nœuds : 294

  • Nombre de mailles et types : 120 HEXA9

Résultats#

Les résultats sont comparés à ssnv302b (calcul 3D overkill), référence AUTRE_ASTER :

Lieu

NOM_CHAM

NOM_CMP

PRECISION

Nœud I

DEPL

DZ

0, 061%

Nœud C

SIEQ_NOEU

VMIS

0, 74%

Nœud E

SIEQ_NOEU

VMIS

0, 65%

Nœud G

SIEQ_NOEU

VMIS

0, 55%

Modélisation C#

Caractéristiques de la modélisation#

On utilise une modélisation COQUE_SOLIDE. On est en traction, plasticité de von Mises à écrouissage isotrope linéaire, petites déformations (VMIS_ISOT_LINE / PETIT).

Caractéristiques du maillage#

Le maillage est le suivant:

  • Nombre de nœuds : 294

  • Nombre de mailles et types : 120 HEXA9

Résultats#

Les résultats sont comparés à ssnv302c (calcul 3D overkill), référence AUTRE_ASTER :

Lieu

NOM_CHAM

NOM_CMP

PRECISION

Nœud I

DEPL

DX

0, 51%

Nœud C

SIEQ_NOEU

VMIS

0, 10%

SIEF_NOEU

SIXX

0,13%

Nœud E

SIEQ_NOEU

VMIS

0, 0023%

SIEF_NOEU

SIXX

0,02%

Nœud G

SIEQ_NOEU

VMIS

0, 00052%

SIEF_NOEU

SIXX

0,0053%

Modélisation D#

Caractéristiques de la modélisation#

On utilise des modélisations3D_SI et COQUE_SOLIDE, liés par un LIAISON_MAILLE. On est en flexion, élasticité isotrope linéaire, petites déformations (ELAS / PETIT).

Caractéristiques du maillage#

La partie en modélisation 3D_SI représente le quart de la longueur de la poutre, près de l’encastrement. Le reste de la poutre est réalisé en modélisation COQUE_SOLIDE. Le maillage n’est pas coïncident à la jonction et celle-ci est réalisée par un LIAISON_MAILLE.

../../../../_images/100000000000046800000246CA313C40DE1553E4.png

Le maillage est le suivant:

  • Nombre de nœuds : 3816

  • Nombre de mailles et types : 720 HEXA9et 420 HEXA20

Résultats#

Les résultats sont comparés à ssnv302a (calcul 3D overkill), référence AUTRE_ASTER :

Lieu

NOM_CHAM

NOM_CMP

PRECISION

Nœud I

DEPL

DZ

4,81%

Nœud C

SIEQ_NOEU

VMIS

0, 0015%

Nœud G

SIEQ_NOEU

VMIS

0, 0046 %

Synthèse des résultats#

Les résultats sont excellents en déplacements et en contraintes par rapport à un calcul très fin en quadratique et en modélisation 3D.