[SOLVE] MEMBRANE 'Erreur de programmation'

Groguiguy · 2018-10-03 13:24:32

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?

Last edited by Groguiguy (2018-10-08 12:59:19)

Attachments:
Plaque_hyp.comm, Size: 2.03 KiB, Downloads: 418

jeanpierreaubry · 2018-10-03 13:36:46

bonjour

je vois que le message d'erreur fait référence à Group_5
je ne vois pas de définition de propriétés d'éléments pour ce groupe mais uniquement pour Group_10

jean pierre aubry

Groguiguy · 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).

Attachments:
message_hyp, Size: 37.06 KiB, Downloads: 414

jeanpierreaubry · 2018-10-03 14:02:26

j'ai fais des exemples fonctionnels, tous avec
des pressions sur la zone et
COMPORTEMENT=_F(RELATION='ELAS_MEMBRANE_SV',DEFORMATION='GROT_GDEP',),

Groguiguy · 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.

jeanpierreaubry · 2018-10-03 15:29:43

ça doit être tout simple
mais
sans le maillage il est impossible de dire

Groguiguy · 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

jeanpierreaubry · 2018-10-03 16:12:39

il arrive que code_aster génère ce genre de message "erreur de programmation"
lorsqu'il rencontre un cas non prévu dans le traitement d'erreur
cela m'est souvent arrivé avec des bourdes "énormes" dans la définition du maillage modèle
tellement énormes que l'on peut passer des heures sans les voir
je redis sans maillage ......

jeanpierreaubry · 2018-10-03 16:21:58

BC4 applique une PRES_REP sur le Group_4 qui me semble être un groupe de ligne SEG2

Groguiguy · 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.

Attachments:
Maillage_plaque.med, Size: 16.52 KiB, Downloads: 413

jeanpierreaubry · 2018-10-03 16:38:46

dans ce cas une PRES_REP qui est une pression (autrement dit en N/m^2) n'a aucun sens d'ou "erreur de programmation"

Groguiguy · 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.

Last edited by Groguiguy (2018-10-03 16:45:03)

jeanpierreaubry · 2018-10-04 08:15:45

ci joint les fichiers d'un problème similaire qui fonctionne
il compare le comportement membrane et le comportement plaque
pour une plaque carrée

la fixation sur deux cotés seulement se révèle problématique

noter ANGL_REP=(0,0) indispensable dans AFFE_CARA_ELEM très mal signalé comme erreur

Attachments:
membranegg.zip, Size: 4.62 KiB, Downloads: 424

Groguiguy · 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é.

Last edited by Groguiguy (2018-10-04 14:48:45)

Groguiguy · 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!

jeanpierreaubry · 2018-10-08 16:55:16

très bien

d'une manière générale je spécifie les paramètes a minima dans le .comm et je laisse code_aster choisir ses valeurs par défaut
tout au moins au premier run

Last edited by jeanpierreaubry (2018-10-08 16:59:48)

AsterO'dactyle · 2018-10-08 18:48:17

Bonjour,

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

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

Groguiguy · 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'.

AsterO'dactyle · 2018-10-09 20:00:12

Bonjour,

En fait je voulais dire que PAS_MINI_ELAS indique que si l'incrément de temps est inférieur à la valeur donnée, on passe de MATRICE='TANGENTE' à MATRICE='ELASTIQUE'.
Comme le pas de temps initial est inférieur à la valeur donnée par PAS_MINI_ELAS, dès le début aster bascule en MATRICE='ELASTIQUE'

Pour reproduire le bug sans PAS_MINI_ELAS, il suffit simplement de mettre MATRICE='ELASTIQUE'.

Groguiguy · 2018-10-10 08:41:52

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

Comment émet-on une fiche sur bitbucket?

Last edited by Groguiguy (2018-10-10 08:42:36)

Thomas DE SOZA · 2018-11-03 13:13:09

https://bitbucket.org/code_aster/codeaster-src/issues

#1 2018-10-03 13:24:32

[SOLVE] MEMBRANE 'Erreur de programmation'

#2 2018-10-03 13:36:46

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

#3 2018-10-03 13:51:07

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

#4 2018-10-03 14:02:26

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

#5 2018-10-03 14:43:43

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

#6 2018-10-03 15:29:43

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

#7 2018-10-03 15:45:59

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

#8 2018-10-03 16:12:39

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

#9 2018-10-03 16:21:58

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

#10 2018-10-03 16:34:20

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

#11 2018-10-03 16:38:46

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

#12 2018-10-03 16:41:37

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

#13 2018-10-04 08:15:45

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

#14 2018-10-04 08:56:25

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

#15 2018-10-08 12:58:53

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

#16 2018-10-08 16:55:16

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

#17 2018-10-08 18:48:17

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

#18 2018-10-09 08:45:54

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

#19 2018-10-09 20:00:12

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

#20 2018-10-10 08:41:52

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

#21 2018-11-03 13:13:09

Re: [SOLVE] MEMBRANE 'Erreur de programmation'

Board footer