Welcome to the forums. Please post in English or French.

You are not logged in.

#1 Code_Aster usage » Calculer la déformation dans le repère global » 2022-05-12 10:43:38

Groguiguy
Replies: 0

Bonjour,

Je souhaiterai comparer des résultats entre un modèle python et le modèle éléments finis de Code_Aster, pour cela j'ai besoin de la déformation dans le repère global mais je ne sais pas comment l'obtenir dans Code_Aster.

Si j'en crois la doc u4.81.04 (page 23) la déformation est calculée dans le repère de l'élément, comment la transposer dans le repère global ?

Soit en faisant un changement de repère, soit en calculant directement via Code_Aster les déformations dans le repère global.

Avez vous une idée de la procédure ?

La doc parle de direction principale, représentée par les vecteurs VECT_n_X etc, mais je ne comprends pas la correspondance avec le repère de l'élément, encore moins avec le repère global.

En vous remerciant,
Guillaume

#2 Re: Code_Aster usage » Force unique de réaction sur une face » 2022-04-01 15:38:32

C'est bon, j'ai fini par y arriver.
Effectivement il y avait bien une erreur dans mon .com pour ne pas sortir convenablement les valeurs, mais lorsque j'eus corrigé l'erreur j'avais alors le classique :

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   ! <S> Exception utilisateur levee mais pas interceptee.      !
   ! Les bases sont fermees.                                    !
   ! Type de l'exception : error                                !
   !                                                            !
   !  le calcul de l'option :  FORC_NODA                        !
   !  n'est possible pour aucun des types d'éléments du LIGREL. !
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Même en définissant convenablement CREA_GROUP_NO=_F( GROUP_MA='nom_du_groupe_de_maille_cible',),

Le problème est ailleurs :
Il faut dans CALC_CHAMP définir un groupe d'éléments volumiques, et non 2D. C'est trompeur sans avoir créé de groupe 3D dans l'outil de maillage Aster Study nous propose naturellement les groupes 2D qui servent à la définition des conditions limites.

#3 Re: Code_Aster usage » Force unique de réaction sur une face » 2022-03-22 17:28:25

Bonjour à tous,

Je n'ai toujours pas trouvé de solution pour calculer l'effort de réaction sur une face en déplacement imposée.
Dois je prévoir un programme Python prenant l'effort à chaque nœud et la pondération selon les lois de la MEF au sujet des efforts distribués ?
Il n'y a pas d'option de Code_Aster pour une réponse sur toute la face ?

En vous remerciant!

#4 Code_Aster usage » Force unique de réaction sur une face » 2022-03-04 17:42:17

Groguiguy
Replies: 4

Bonjour à tous,

Je souhaiterai obtenir l’effort de réaction orthogonal à une face contrainte en DDL, l'effort globale mesuré dans des mords d'un essai de traction simulé par exemple.
J'ai déjà lu ce sujet (code-aster.org/forum2/viewtopic.php?id=19987) et le livre de Jean-Pierre Aubry qui mentionne l'utilisation de :

#U4.81.21
sum_reac=POST_RELEVE_T (
ACTION=_F (
INTITULE=’sum reactions’ ,
GROUP_NO = ( ’groundS’ , ’groundN’ , ) ,
RESULTAT=stat ,
NOM_CHAM=’REAC_NODA’ ,
TOUT_ORDRE=’OUI’ ,
RESULTANTE = ( ’DX’ , ’DY’ , ’DZ’ , ) ,
MOMENT = ( ’DRX’ , ’DRY’ , ’DRZ’ , ) ,
POINT = ( 0 , 0 , 0 , ) ,
OPERATION=’EXTRACTION’ ,
               ) ,
) ;
#U4.91.03
IMPR_TABLE ( TABLE=sum_reac , FORMAT_R=’1PE12.3’ , )

Cependant cette méthode sort une valeur de force de réaction pour chaque noeud de la face sélectionnée.
Hors, la résultante n'est pas la somme de ces efforts, notamment parce que des efforts vont être négatifs et positifs dans un même élément (différence entre DLOAD et DSLOAD de mémoire).
Si j'en crois la doc u4/u4.81.04 il faut appliquer TOUT=OUI à la définition du groupe observé, mais dans ce cas je ne comprends pas la différence entre Reac_Noda et Force_Noda puisque le calcul est effectué partout.

J'ai forcément mal compris quelque chose!
En vous remerciant!

#5 Re: Code_Aster usage » Problème de convergence en Hyperelastic » 2019-01-17 17:33:19

Alors j'ai un truc qui cloche, car je viens de faire le test avec 2400, et il réécrit 196!!!

PS : je précise que j'utilise appli_V2017.0.2 mais je ne pense pas que ça change grand chose.

#6 Re: Code_Aster usage » Problème de convergence en Hyperelastic » 2019-01-17 17:16:05

Est ce que dans le fonctionnement de Code_Aster il fait appelle à ASTK en off, et ce serait ASTK qui réécrirait "tpsjob 196" dans le fichier export?
Si c'est ça, je ne vois pas comment contourner le problème, mais je ne vois pas pourquoi ça ne le faisait pas avant non plus.
Dans la doc u1.04.00 il est justement fait mention que l'emploi de as_run utilise ASTK... le problème vient peut être de là.

#7 Re: Code_Aster usage » Problème de convergence en Hyperelastic » 2019-01-17 15:31:11

Je n'utilise pas ASTK car j'intègre l'ensemble de mes calculs dans un futur logiciel piloté en python où j'écris directement dans les fichiers export et .comm. Je pense que mes calculs vont prendre pas loin de 80h, voir plus vu que je fais des tests de convergence.

Le truc bizarre ici, c'est qu'en lançant le calcul via shell 'as run fichier.export' mon fichier export est réécrit automatiquement et la valeur tpsjob est changée!!! Je comprends pas comment Code_aster peut avoir la main pour modifier le fichier export que j'ai moi même écrit et lancé.
Le comble étant que du coup, pas moyen d'augmenter le temps CPU vu qu'il remet tpsjob 196 qu'il me demande lui même d'augmenter!

La lecture du fichier .comm pourrait influencer le fichier export? Vraiment je pige pas. De plus je n'avais pas ce problème avant, donc c'est sans doute une erreur de ma part dans l'écriture...

#8 Re: Code_Aster usage » Problème de convergence en Hyperelastic » 2019-01-17 13:59:45

Je me permet de remonter le sujet car j'ai toujours un problème avec le temps CPU (Arrêt pour manque de temps CPU).

Dans mon fichier export, je le défini comme suis :

P time_limit 11700.0
P tpsjob 19600

Où il me semble bien que c'est tpsjob qui défini la limite. Seulement, dès que je lance le calcul par le terminal linux, le fichier export est réécrit avec tpsjob 196 ce qui n'est pas suffisant. Comment faire pour augmenter le temps CPU du coup?

J'ai aussi défini cela dans le fichier .comm lié :

RESERVE_CPU=_F(BORNE=90000,)

Mais je sais pas si c'est suffisant.

#9 Re: Code_Aster usage » Deformed shape with COQUE_3D » 2018-12-17 18:03:10

Ok, thank you very much, I will try this!

#10 Code_Aster usage » Coque + 3D elements » 2018-12-12 17:40:09

Groguiguy
Replies: 0

Hi everebody,

What's happend when I use COQUE and 3D elements together in same model?
At the connection between both elements types, I don't know how Code_Aster transfert energy of moment.

DOF Nodes in COQUE :
-DX
-DY
-DZ
-RX
-RY
-RZ

DOF Nodes in 3D elements :
-DX
-DY
-DZ

What's happend when COQUE in flexion is used with volume?
Doc R3.03.03 explains Connections 3D-Beam,  2D-Beam, but what about 3D-Coque?

Thank you.

PS : //www            .code-aster.org/V2/doc/v13/fr/man_v/v3/v3.05.100.pdf
That works, but I don't know how.

#11 Code_Aster usage » Deformed shape with COQUE_3D » 2018-12-10 14:14:53

Groguiguy
Replies: 2

Hello everyone,

I would like to see deformed shape (or Warp by Vector in Paravis) for my model, but I use MODELISATION = COQUE_3D and if I understand correctly, Paravis can't show Warp by Vector with elements TRIA7 who works with COQUE_3D.

How can I see the deformed shape? Other software?

Of course, I can realize python program to transform nodes coordonates plus displacements to new coordonates. But may be is it possible to avoid that.

Thank you!

#12 Re: Code_Aster usage » [SOLVE]Membrane et Solid en maillage conforme. » 2018-10-26 13:50:55

C'est bon j'ai trouvé!
Il me suffisait de n'utiliser qu'un seul maillage, et pas d'en créer un autre qui se retrouvait avec des mouvements de corps rigides.
En fichier joint on peut voir comment j'ai défini le problème, très simplement finalement.
Cela marche très bien, même en grandes déformations, c'est super!

#13 Re: Code_Aster usage » [SOLVE]Membrane et Solid en maillage conforme. » 2018-10-24 14:39:20

Oui oui bien sur, ce n'était qu'une erreur de normal à la surface ici, mais je parlais de l'écriture générale.
Depuis, par exemple, j'ai cru comprendre qu'il fallait mettre 2 COMPORTEMENTS pour la résolution d'un problème MEMBRANE + 3D.

J'ai l'impression d'avoir bien défini le problème, pourtant Code_Aster m'annonce ceci :

!---------------------------------------------------------------------------------------------!
   ! <EXCEPTION> <FACTOR_42>                                                                     !
   !                                                                                             !
   ! Solveur MUMPS :                                                                             !
   !   La Matrice est non factorisable. Elle est détectée comme étant singulière                 !
   !   (en structure ou numériquement).                                                          !
   !                                                                                             !
   ! Conseils :                                                                                  !
   !   Il peut s'agir d'une erreur de programmation ou d'un problème de mise en données (blocage !
   !   absent ou surabondant).                                                                   !
   !---------------------------------------------------------------------------------------------!

Et là, malgré la simplicité du cas, je ne comprends pourquoi ce message.
Je cherche juste à modéliser un Cube (GROUP_MA = Cube) où l'une des faces comprends une "peau" avec des éléments MEMBRANE (GROUP_MA = GROUP_10). Puis faire un essais de traction simple. Le but est juste de montrer qu'il est possible de faire ce genre de modélisation en maillage2D/3D conforme sans contact.

#14 Code_Aster usage » [SOLVE]Membrane et Solid en maillage conforme. » 2018-10-19 13:37:22

Groguiguy
Replies: 3

Bonjour à tous,

Avançant à grand pas sur nos modélisations on se rends compte que Code_Aster a peut être une méthode qui nous permettrait de nous passer de la modélisation avec des contacts entre éléments membrane et solide. Quand un maillage 3D est lu par Code_Aster, ce dernier distingue 3 types d'éléments :
- des éléments linéaires
- des éléments surfaciques,
- des éléments volumiques,

Logique, puisque le mailleur employé (NetGen 1D-2D-3D) crée bien ces 3 types d'éléments pour obtenir un maillage volumique.

Du coup, nous aimerions utiliser cette dissociation d'éléments pour affecter différentes propriétés au maillage, exemple : faire un cube en acier, avec une "peau" autour en aluminium. D'ordinaire nous utiliserions du contact, mais il me semble possible d'affecter une modélisation MEMBRANE au bord du maillage, et 3D dans le cube. Mais je n'y arrive pas (voir fichier joint).

Pensez vous cela possible?

PS : pour tester cela, un cube a été maillé en 3D et décomposé en 2 maillages ; un groupe contenant tous les noeuds (groupe CUBE) et un groupe contenant les noeuds à la surface du cube ('Group_7','Group_8','Group_9','Group_10','Group_11','Group_12').

#15 Re: Code_Aster usage » [SOLVE] MEMBRANE 'Erreur de programmation' » 2018-10-10 08:41:52

Effectivement, j'avais oublié ce détail, merci!

Comment émet-on une fiche sur bitbucket?

#16 Re: Code_Aster usage » [SOLVE] MEMBRANE 'Erreur de programmation' » 2018-10-09 08:45:54

Il faut émettre une fiche sur bitbucket pour qu'on pense à corriger ça.

Volontiers, mais comment faire? Un tuto/exemple?

En fait, de manière générale, la matrice élastique n'est pas programmée dans les éléments de membrane.
Donc si vous essayez MATRICE='ELASTIQUE', vous aurez le même problème

Même pas, j'utilise MATRICE='TANGENTE'.

#17 Re: Code_Aster usage » [SOLVE] MEMBRANE 'Erreur de programmation' » 2018-10-08 12:58:53

Trouvé!

Le problème "erreur de programmation" vient de cette ligne :

PAS_MINI_ELAS=0.1,

En l'effaçant le calcul passe sans soucis.
J'en conclu qu'avec la modélisation Membrane la fonction PAS_MINI_ELAS n'a pas dû être implémentée, du coup le message d'erreur non plus!

Je n'aurai pas trouvé sans la comparaison avec le fichier jeanpierreaubry! D'autant que j'utilise souvent les mêmes lignes de code pour définir mes Stat_nonlin dans différents cas. Merci encore!

#18 Re: Code_Aster usage » [SOLVE] MEMBRANE 'Erreur de programmation' » 2018-10-04 08:56:25

Parfait merci, je vais regarder ça attentivement!


EDIT : Erreur, je n'ai pas corrigé le problème contrairement à ce que disait un message suivant que j'ai effacé.

#19 Re: Code_Aster usage » [SOLVE] MEMBRANE 'Erreur de programmation' » 2018-10-03 16:41:37

Mais pourtant on rentre bien une épaisseur via EPAIS, donc pour moi cela conservé un sens. De toute façon, si je repasse avec un BC sous la forme :
BC3 = AFFE_CHAR_CINE(MECA_IMPO=(_F(DX=10.0, GROUP_MA=('Group_4', ))), MODELE=model)
ça ne change pas le message d'erreur.
Même en mettant un déplacement imposé.

Il y a forcément quelque chose que j'ai raté dans Code_Aster en général.
Il m'arrive assez souvent d'avoir cette erreur de mailles n'appartenant pas au model. Pourtant, dans la définition du model il me semble que TOUT='OUI' suffit à définir la physique sur l'ensemble.

Quoique, en l'expliquant ici je me rends compte que je ne peux définir l'élément MEMBRANE sur des SEG2... mais que faut il alors faire pour les conditions limites? Elles sont souvent définis sur des bords de maillage.

#20 Re: Code_Aster usage » [SOLVE] MEMBRANE 'Erreur de programmation' » 2018-10-03 16:34:20

Oui c'est ça, le Groupe_4 est une arête, puisque je souhaite appliquer une condition là dessus. Donc c'est bien un groupe de SEG2.

#21 Re: Code_Aster usage » [SOLVE] MEMBRANE 'Erreur de programmation' » 2018-10-03 15:45:59

Le voici.

Même si le problème vient de là, quelle est cette autre "erreur de programmation"? Je suppose que c'en est pas une.

Je cherche juste à étudier le comportement des membranes de Code_Aster en hyperélastic, notamment refaire ce modèle :
https:               //www.code-aster.org/spip.php?article997

#22 Re: Code_Aster usage » [SOLVE] MEMBRANE 'Erreur de programmation' » 2018-10-03 14:43:43

Où puis je trouver ces exemples?

Là même avec ELAS_MEMBRANE_SV ça ne fonctionne pas.

J'ai changé le chargement pour une pression mais j'obtiens alors :

BC4 = AFFE_CHAR_MECA(MODELE=model,
                       PRES_REP=_F(GROUP_MA=('Group_4', ),
                                   PRES=0.1,),
                       INFO=1,
                       VERI_NORM='OUI',)

  MAILLE : M15     
  MAILLE : M16     
  MAILLE : M17     
  MAILLE : M18     
  MAILLE : M19     
  MAILLE : M20     
  MAILLE : M21     
   
   !-------------------------------------------------------------------------!
   ! <EXCEPTION> <MODELISA6_96>                                              !
   !                                                                         !
   !  les 7 mailles imprimées ci-dessus n'appartiennent pas au modèle        !
   !  et pourtant elles ont été affectées dans le mot-clé facteur : PRES_REP !
   !-------------------------------------------------------------------------!

Et là non plus je ne comprends pas pourquoi cette erreur. Mon model est défini par TOUT='OUI' donc je ne vois pas comment ces éléments peuvent ne pas faire parti du modèle.

#23 Re: Code_Aster usage » [SOLVE] MEMBRANE 'Erreur de programmation' » 2018-10-03 13:51:07

Zut je me suis trompé de fichier en lien, désolé.
J'ai mis un autre fichier message en pièce jointe là, plus claire car plus court et où on voit le message d'erreur.
La plaque est décomposée en 5 groupes ; 4 arrêtes (groupe 1 à 4) et la plaque entière (Group_5).

#24 Code_Aster usage » [SOLVE] MEMBRANE 'Erreur de programmation' » 2018-10-03 13:24:32

Groguiguy
Replies: 20

Hi,

I don't know why my model with MEMBRANE doesn't work. It's just a tensile test plate.
And error message is :

!-------------------------------------------------------------------------------------------------------------------------------------!
   ! <EXCEPTION> <DVP_1>                                                                                                                 !
   !                                                                                                                                     !
   ! Erreur de programmation.                                                                                                            !
   !                                                                                                                                     !
   ! Condition non respectée:                                                                                                            !
   !     .false.                                                                                                                         !
   ! Fichier                                                                                                                             !
   ! /home/A21173/smeca/yamm_build/V2017_0_2_build/V2017_public/tools/src/Code_aster_stable-v134_smeca/bibfor/elements/mbgnlr.F90, ligne !
   ! 91                                                                                                                                  !
   !                                                                                                                                     !
   ! --------------------------------------------                                                                                        !
   ! Contexte du message :                                                                                                               !
   !    Option         : RIGI_MECA_ELAS                                                                                                  !
   !    Type d'élément : MEMBTR3                                                                                                         !
   !    Maillage       : mesh                                                                                                            !
   !    Maille         : M29                                                                                                             !
   !    Type de maille : TRIA3                                                                                                           !
   !    Cette maille appartient aux groupes de mailles suivants :                                                                        !
   !       Group_5                                                                                                                       !
   !    Position du centre de gravité de la maille :                                                                                     !
   !       x=-452.380952 y=-452.380952 z=0.000000                                                                                        !
   !                                                                                                                                     !
   !                                                                                                                                     !
   !                                                                                                                                     !
   ! Il y a probablement une erreur dans la programmation.                                                                               !
   ! Veuillez contacter votre assistance technique.                                                                                      !
   !-------------------------------------------------------------------------------------------------------------------------------------!

message_file is join.
My mesh is just a simple plate, with triangular elements (linear or quadratic). I don't understand what's happening!
I suppose that is matrix factorisation error.

Have you some example to test MEMBRANE modelisation?

#25 Re: Code_Aster usage » Element DKT + 3D » 2018-09-17 12:33:50

Merci, le lien vers le gonflement des éléments membranes est effectivement très intéressant.
Mais lorsque je veux tester ces éléments sur un cas simple de plaque en traction j'ai un message d'erreur me signalant :

Il y a probablement une erreur dans la programmation.
Veuillez contacter votre assistance technique.

Vous trouverez le ficher message en pièce jointe pour plus d'information. J'avoue ne pas comprendre le problème.