6. Eléments de structure#

6.1. Généralités#

  • Calcul d’un champ EGRU d’efforts généralisés en repère global (#33252/33584 - 17.1)

    Dans la commande MODI_REPERE, un nouveau mot-clé NOM_CHAM_RESU a été ajouté et ne concerne que le champ EFGE_ELNO. Si NOM_CHAM_RESU=EGRU_ELNO est indiqué, alors un nouveau champ nommé « EGRU_ELNO » est créé et contient le champ EFGE_ELNO calculé dans le nouveau repère. La sd resultat du MODI_REPERE contiendra également EFGE_ELNO qui reste inchangé (LE CHAMP INITIAL). Si NOM_CHAM_RESU est omit, alors le comportement initial de MODI_REPERE est conservé et le champ renseigné dans NOM_CHAM se retrouve modifié (écrasé). Ces champs sont utilisables avec l’opérateur COMB_SISM_MODAL.

  • Blocage de la normale dans AFFE_CHAR_MECA (#33329 - 17.1)

    On introduit un nouveau mot-clé DRNOR dans AFFE_CHAR_MECA / FACE_IMPO afin d’imposer des déplacements suivant la normale locale.

  • Amélioration performance de l’opérateur AFFE_CARA_ELEM (#34561 - 17.4)

    On modifie le mot-clé VERIF =[OUI | NON]. ( defaut = OUI). Les vérifications sont maintenant centralisées et le mot-clé VERIF permet de les désactiver si besoin.

6.2. Eléments plaques, coques et membranes#

  • Nouveau mot-clé FORCE_COQUE_FO dans AFFE_CHAR_MECA (#33522 - 17.1)

    Ce nouveau mot-clé FORCE_COQUE_FO permet de fournit un chargement dépendants de X, Y, Z , EP (epaisseur), et ou RHO. On évalue ces fonctions dans l’opérateur pour chaque maille. Cette évaluation n’est plus nécessaire à chaque pas de temps comme c’est le cas avec AFFE_CHAR_MECA_F/FORCE_COQUE.

  • Vérification planéité des plaques (#33517 - 17.1)

    La fiche propose de déplacer la vérification de la planéité des plaques (QUAD4) dans AFFE_MODELE au lieu de l’émettre à chaque itération. L’option VERI_PLAN est créée. Elle utilise le même critère de planéité que l’ancien (dans le te). On ajoute le mot-clef VERI_PLAN dans AFFE_MODELE.

  • Définition d’un repère cylindrique pour AFFE_CARA_ELEM / GRILLE_MEMBRANE (#33303 - 17.1)

    Pour déclarer des repères locaux dans AFFE_CARA_ELEM / GRILLE, l’utilisation de VECT_1 ou VECT_2 est pratique et simple. Mais il existe des configurations géométriques qui ne permettent l’usage ni de VECT_1, ni de VECT_2 (repère polaire dans un plan horizontal). Avec le nouveau choix, REPERE='CYLINDRIQUE', ORIGINE et AXE_Z, les vecteurs VECT_1 / VECT_2 peuvent maintenant être exprimés dans le repère cylindrique.

  • Gestion des matrices non-symétriques pour les DKT (#31432 - 17.2)

    On ajoute le traitement du cas des matrices non symétriques pour les éléments DKT.

  • Nouvel élément plaque via la bibiothèque Fenics (#34936 - 17.4)

    Cette fiche traite de l’ajout de la modélisation PLAQ_MITC sur un support QUAD9 (POC sur l’implémentation d’un élément fini avec la bibliothèque Fenics). La modélisation adoptée pour les plaques repose sur la théorie de Reissner–Mindlin, enrichie par une interpolation mixte de type MITC (Mixed Interpolation of Tensorial Components), dans la variante proposée par Durán et Liberman. Cette approche permet de corriger efficacement le verrouillage numérique lié au cisaillement transverse, assurant ainsi une meilleure précision dans la reproduction des déformations, même pour des plaques minces.

6.3. Eléments poutres et tuyaux#

  • Gestion Poutres proches de l”axe Z avec GROT_GDEP (#33398 - 17.2)

    Les calculs de poutres en GROT_GDEP, quand celles-ci sont orientées selon l’axe Z ou très proche de Z, ont parfois des difficultés à converger. L’actualisation du repère local des poutres en grands déplacements est en cause. Une solution plus robuste est mise en place pour améliorer la convergence.

  • Ajout de possibilités de definition de l”orientation des sections dans AFFE_CARA_ELEM (#31109 - 17.3).

    Lorsque l’utilisateur utilise MACR_CARA_POUTRE, en sortie il obtient une table avec les caractéristiques mécaniques nécessaires à AFFE_CARA_ELEM. Les caractéristiques nécessaires à AFFE_CARA_ELEM sont toujours celles qui sont calculées dans le repère principal d’inertie. On ajoute plusieurs manières de définir la section de poutre autour de son axe neutre ( possible avec « VECT_Y », « ANGL_VRIL »):

    • VECT_Z : qui est similaire à VECT_Y mais pour définir l’axe z3.

    • VECT_MAIL_Y/VECT_MAIL_Z : qui est similaire à VECT_Y/Z mais qui permet de définir la projection des axes Y/Z du maillage donné en entrée de MACR_CARA_POUTRE. On a donc besoin de l’angle “Alpha” qui est dans la table et qui donne l’orientation du repère principal d’inertie par rapport au repère du maillage.

  • Evolution des données en sortie de MACR_CARA_POUTRE (#31038 - 17.3)

    Ce développement consiste à

    • Modifier les propriétés pour que le ALPHA soit toujours dans le premier cadran.

    • Ajouter les coordonnées du centre de torsion dans le repère du maillage.

  • Ajout d’un chargement PRE_EPSI en repère global pour POU_D_T (#34241 - 17.3)

    On ajoute les mots-clé KN et VECT_N dans AFFE_CHAR_MECA et AFFE_CHAR_MECA_F pour PRE_EPSI. Ce chargement doit s’appliquer aux éléments POU_D_E et POU_D_T. Les fonctions dépendent de l’abscisse curviligne.

  • Ajout d’un post-traitement pour les éléments poutres du champ EFGE_ELNO (#34544 - 17.3).

    Étant donné un RESU de type DYNA, on calcule un champ maximal du même type que EFGE_ELNO. Pour les composantes FX, FY, FZ, on prend la valeur absolue maximale sur l’ensemble des instants du résultat. Pour les composantes MX, MY, MZ, on sélectionne l’instant où la quantité MX² + MY² + MZ² est maximale, puis on considère les valeurs absolues des composantes correspondantes.

  • Amélioration de l’intégration des éléments Tuyaux (#32525 / #34450 / #34557 /#30681 - 17.3) via la sous-intégration des intégrales.

    Parmi les autres développements, on note l’interdiction de MODI_METRIQUE=OUI pour les tuyaux (#34450 17.3), la sous-intégration de la matrice de masse (#34557 17.3) et la reprise des raccords Coque/Tuyau et 3D/Tuyau sont mal écrits dans le cas des grands déplacements. (#30681 -17.3)

  • Refactoring éléments de Tuyau : séparation formulation poutre et coque (#34415 - 17.4)

    Ce développement traite de la reprise de la programmation des éléments tuyaux en essayant de séparer informatiquement la partie poutre (courbe) et la décomposition de Fourier (coque) et en anonymisant au maximum des données d’entrée.

6.4. Eléments discrets#

  • Dilatation d’un discret soumis à un champ de température (#33437 - 17.1)

    On souhaite prendre en compte la dilatation thermique des discrets. Les développements concernent seulement les discrets K_TR_D_L (en remplacement des poutres, mêmes DDL). Les déplacements aux noeuds dus à la température pour les poutres est du type α*(T-Tref)*L, et dépend donc de la longueur de l’élément. Les déplacements aux noeuds pour les discrets seront de la même forme et donc dépendront de la longueur du discret. Le coefficient de dilatation « alpha » aura donc la même signification pour les poutres et les discrets.

  • Développement de la méthode des grilles (#34612 - 17.4)

    Cette fiche traite de l’implémentation de la méthode RIGI_GRILLE. Il s’agit d’un nouveau mot-clé dans l’opérateur AFFE_CARA_ELEM. Cette fonctionnalité correspond à une méthodologie utilisée pour déterminer les caractéristiques d’éléments discrets à appliquer au radier à partir de raideurs globales. Cette option est disponible en 3D. Le radier sera modélisé par une surface, les discrets sont des DIS_TR. Il faut distinguer un groupe de mailles pour le radier, à déclarer derrière le mot clé GROUP_MA et les 5 groupes de mailles de type SEG2 à déclarer derrière le mot clé GROUP_MA_SEG2. Ces mailles SEG2 doivent être affecté par des DIS_TR dans AFFE_MODELE. Ces mailles possède chacune un noeud en commun avec le radier, et un noeud sur lequel une LIAISON_SOLIDE sera appliquée par ailleurs.