Atom topic feed | site map | contact | login | Protection des données personnelles | Powered by FluxBB | réalisation artaban
You are not logged in.
Hello,
I have a big model, with different AFFE_CARA_ELEM and different AFFE_MATERIAU. Also I have test element quality. No I get following error.
! <EXCEPTION> <COMPOR5_1> !
! !
! On ne trouve pas la courbe de traction (mot-clef TRACTION) dans le mat?riau fourni. !
! !
! - !
! Contexte du message : !
! Option : RIGI_MECA_TANG !
! Type d'?l?ment : MEDKTR3 !
! Maillage : NETZ !
! Maille : M30474 !
! Type de maille : TRIA3 !
! Cette maille appartient aux groupes de mailles suivants : !
! EL2DPR01 ALLE2DEL EL2DTRIA !
! Position du centre de gravit? de la maille : !
! x=1312.662005 y=665.706398 z=1507.902874 !
!-!
!-!
! <EXCEPTION> <COMPOR5_1> !
! !
! On ne trouve pas la courbe de traction (mot-clef TRACTION) dans le mat?riau fourni. !
! !
! - !
! Contexte du message : !
! Option : RIGI_MECA_TANG !
! Type d'?l?ment : MEDKTR3 !
! Maillage : NETZ !
! Maille : M3004 !
! Type de maille : TRIA3 !
! Cette maille appartient aux groupes de mailles suivants : !
! EL2DPR01 ALLE2DEL EL2DTRIA ... !
! Position du centre de gravit? de la maille : !
! x=1019.073286 y=708.970468 z=37.550732 !
!-!
The listed Maille are over all good mesh quality. Also only "0,1%" this elements in GROUP_MA="..." brings <EXCEPTION> <COMPOR5_1> the rest off 99,9% are correct.
Why brings only less elements in the GROUP_MA="..." this "..traction (mot-clef TRACTION) dans le mat?riau fourni."? About informations I'm very happy.
Thanks Markus
The superfluous chase, is to miss the essential.
Jules Saliège
Offline
Hallo.
I'm not sure about myself but:
Why do you think this is a mesh issue ?
Reading u4.43.01 §4.1.1, isn't the error message related to the TRACTION keyword of DEFI_MATERIAU ?
I read:
TRACTION = _F(SIGM =sigm_f),
Do you have defined any sigm_f fonction ?
I can translate:
! On ne trouve pas la courbe de traction (mot-clef TRACTION) dans le mat?riau fourni. !
to:
'We cannot find the traction curve (keyword TRACTION) of the given material."
if it can help.
Last edited by Voulet (2019-05-06 15:56:56)
Offline
Hello,
very thanks. My functions for material looks like:
S355Funk=DEFI_FONCTION(NOM_PARA='EPSI',
VALE=(0.00169 ,355.0 ,
0.08 ,500.0 ,
5.0 ,510.0 ,),
INTERPOL='LIN',
PROL_DROITE='LINEAIRE',
PROL_GAUCHE='EXCLU',);
S355=DEFI_MATERIAU(ELAS=_F(E=210000.0,
NU=0.3,
RHO=7.85e-09,),
TRACTION=_F(SIGM=S355Funk,),);
Is this a mistake? But why by so less elements?
Thanks Markus
The superfluous chase, is to miss the essential.
Jules Saliège
Offline
S355Funk=DEFI_FONCTION(NOM_PARA='EPSI',
VALE=(0.00169 ,355.0 ,
0.08 ,500.0 ,
5.0 ,510.0 ,),
INTERPOL='LIN',
PROL_DROITE='LINEAIRE',
PROL_GAUCHE='EXCLU',);
Maybe the curve is not defined for epsi < 0.00169 because of PROL_GAUCHE = 'EXCLU' ?
Maybe these elements have a lower epsi ?
Offline
Hello,
very thanks. I dont understand the problem, because the rest of 99,9% in GROUP_MA('EL2DPR01','EL2DPR02','EL2DPR04',) are correct.
Maybe the curve is not defined for epsi < 0.00169 because of PROL_GAUCHE = 'EXCLU' ?
I think epsi < 0.00169 then ELAS definition works?
Thanks Markus
Last edited by m_golbs (2019-05-06 16:54:51)
The superfluous chase, is to miss the essential.
Jules Saliège
Offline
hello
i would say that [as usual] the problem is probably not lying where one thinks it is at first
the .comm and mesh file would help to understand
jean pierre aubry
consider reading my book
freely available here https://framabook.org/beginning-with-code_aster/
Offline
Hallo
is it possible to use test alternativ? How works the ELAS sector on this end?
S355Funk=DEFI_FONCTION(NOM_PARA='EPSI',
VALE=(0.00169 ,355.0 ,
0.08 ,500.0 ,
5.0 ,510.0 ,),
INTERPOL='LIN',
PROL_GAUCHE='LINEAIRE',
PROL_DROITE='LINEAIRE',);
# PROL_GAUCHE='EXCLU',);
S355=DEFI_MATERIAU(ELAS=_F(E=210000.0,
NU=0.3,
RHO=7.85e-09,),
TRACTION=_F(SIGM=S355Funk,),);
Can I use DEFI_MATERIAU only with TRACTION function, without ELAS?
It is a big model with 20 AFFE_CARA_ELEM and 5 DEFI_MATERIAU 2D DKT, STAT_NON_LINE, PETIT_REAC, MPI .... I will try different employments.
[as usual] the problem is probably not lying where one thinks it is at first
the .comm and mesh file would help to understand..
Yes this is the correct way. But I think the *.comm is not the primary starting point. I have test the *comm with a big Modell only one DKT group, equable CAD Modell >> equable mesh. It works all right.
How can I find problematic element quality with *.comm file implement?
Thanks Markus
Last edited by m_golbs (2019-05-06 17:46:14)
The superfluous chase, is to miss the essential.
Jules Saliège
Offline
How can I find problematic element quality with *.comm file implement?
looking at the message i am ready to bet this is not a quality problem
but more probably a problem of group assignement
but as i said only .comm and mesh file would help
It is a big model with 20 AFFE_CARA_ELEM and 5 DEFI_MATERIAU 2D DKT, STAT_NON_LINE, PETIT_REAC, MPI
i think i have already seen or made problem of that size or complexity
i beleive i could understand it!
Can I use DEFI_MATERIAU only with TRACTION function, without ELAS?
i do ot think so
consider reading my book
freely available here https://framabook.org/beginning-with-code_aster/
Offline
Hallo,
many Thanks! The *.unv mesh file is nearly 100MB. It is a development project. I Cant post it complete. I have modify the last points of DEFI_FONCTION, upper "E-modul". Then I dont have the error "On ne trouve pas la courbe de traction (mot-clef TRACTION) dans le mat?riau fourni". It is not the absolut solution, but works now without this error.
i think i have already seen or made problem of that size or complexity..
This is also my problem. Little models works, but big modells brings question. The project works fine with MECA_STATIQUE little depl, little stress/stains, also with STAT_NON_LINE without COMPORTEMENT= .. but with STAT_NON_LINE and COMPORTEMENT= I have a convergence offset.
» STAT_NON_LINE basic convergence question
Thanks and greeting Markus
Last edited by m_golbs (2019-05-09 14:20:23)
The superfluous chase, is to miss the essential.
Jules Saliège
Offline
STAT_NON_LINE without COMPORTEMENT= ..
what is the use of that?
i do not think this is allowed or it uses default values
consider reading my book
freely available here https://framabook.org/beginning-with-code_aster/
Offline
Hello,
many thanks! I wonder, that in u4.51.03 no information to COMPORTEMENT. Therefore the question to default in STAT_NON_LINE to the feature COMPORTEMENT. STAT_NON_LINE without COMPORTEMENT runs basically.
Total 177785.81 16693.26 194479.07 5077.32
---------------------------------------------------------------------------------
(*) cpu and system times may be not correctly counted using mpirun.
as_run 2018.0
------------------------------------------------------------
--- DIAGNOSTIC JOB : <A>_ALARM
------------------------------------------------------------
EXIT_CODE=0
But under with conditions to come the results?
We have real test results by this project. STAT_NON_LINE with
...
...
KD2EVC=AFFE_MODELE(MAILLAGE=NETZ,
AFFE=(_F(TOUT='OUI',
PHENOMENE='MECANIQUE',
MODELISATION='3D',),
_F(GROUP_MA=('EL2DPR01','EL2DPR02','EL2DPR03','EL2DPR04','EL2DPR05','EL2DPR06','EL2DPR07','EL2DPR08',
'EL2DPR09','EL2DPR10','EL2DPR11','EL2DPR12','EL2DPR13','EL2DPR14','EL2DPR15',
'EL2DPR16','EL2DPR18','EL2DPR19','EL2DPR30','EL2DPR50',),
PHENOMENE='MECANIQUE',
MODELISATION='DKTG',),),);
..
..
..
S355Funk=DEFI_FONCTION(NOM_PARA='EPSI',
VALE=(0.00169 ,355.0 ,
0.08 ,500.0 ,
0.5 ,510.0 ,),
INTERPOL='LIN',
PROL_GAUCHE='LINEAIRE',
PROL_DROITE='LINEAIRE',);
# PROL_GAUCHE='EXCLU',);
S355=DEFI_MATERIAU(ELAS=_F(E=210000.0,
NU=0.3,
RHO=7.85e-09,),
TRACTION=_F(SIGM=S355Funk,),);
....
S420Funk...
S550Funk...
DC04Funk...
SG03Funk...
....
....
ERGEBNIS=STAT_NON_LINE(MODELE=KD2EVC,
CHAM_MATER=MATERIAL,
CONTACT = KONTAKT,
CARA_ELEM=CARA,
EXCIT=(_F(CHARGE=LAGER01,
TYPE_CHARGE='FIXE_CSTE',),
_F(CHARGE=Lastseit,
FONC_MULT=Lastfunk,
TYPE_CHARGE='FIXE_CSTE',),),
COMPORTEMENT=(_F(DEFORMATION='GROT_GDEP',
RELATION='ELAS_VMIS_TRAC',
TOUT='OUI',),
_F(DEFORMATION='GROT_GDEP',
GROUP_MA=('ALLE2DEL', ),
RELATION='ELAS_VMIS_TRAC',),),
INCREMENT=_F(LIST_INST=SolSchrX,
PRECISION=1e-06,),
METHODE='NEWTON',
# NEWTON=_F(MATRICE='ELASTIQUE',
# MATRICE='TANGENTE',
# PREDICTION='EXTRAPOLE',),
# PREDICTION='TANGENTE',
CONVERGENCE=_F(RESI_GLOB_MAXI=10,
# RESI_GLOB_RELA=1e-06,
ARRET='NON',
ITER_GLOB_MAXI=25,),
SOLVEUR=_F(METHODE='MUMPS',
RENUM='METIS',
NPREC=12,
ELIM_LAGR='NON',
STOP_SINGULIER='NON',),
ARCHIVAGE=_F(LIST_INST=ArcSchri,
CRITERE='RELATIF',
PRECISION=1e-06,),);
bring only 5% of real test fail force. In my the code is also an error, but where? I think my mistake located in COMPORTEMENT, possibly in DEFI_MATERIAU...
Thanks and greeting Markus
Last edited by m_golbs (2019-05-09 17:30:39)
The superfluous chase, is to miss the essential.
Jules Saliège
Offline
Hi Markus,
Did your Problem get solved? I have exactly the same Problem...
I tried to remesh and tried other Things. So far no success.
Any help would be great.
Cheers,
FEUser
Offline
Hello,
unfortunately no. Next I will Mr. jeanpierreaubry hint test "..but more probably a problem of group assignement.."
TL61GRA3=AFFE_MODELE(MAILLAGE=NETZ,
INFO=1,
VERI_JACOBIEN='OUI',
DISTRIBUTION=_F(METHODE='CENTRALISE',),
AFFE=(_F(TOUT='OUI',
PHENOMENE='MECANIQUE',
MODELISATION='3D',),
_F(GROUP_MA=('EL2DPR01','EL2DPR02','EL2DPR03','EL2DPR04','EL2DPR05','EL2DPR06','EL2DPR07','EL2DPR08',
'EL2DPR09','EL2DPR10','EL2DPR11','EL2DPR12','EL2DPR13','EL2DPR14','EL2DPR15',
'EL2DPR16','EL2DPR18','EL2DPR19','EL2DPR30','EL2DPR50',),
PHENOMENE='MECANIQUE',
MODELISATION='COQUE_3D',),),);
name definition of groups
EL: Elements (NO for nodes,..)
2D: Typ of dimension
PR: properties ( thickness >> material...)
01: consecutive number
Means "..but more probably a problem of group assignement.." that are one element in two groups and so on...?
Thanks and greeting Markus
The superfluous chase, is to miss the essential.
Jules Saliège
Offline
Hi all,
Did anyone got problem resolved? I got the same error when I'm trying to simulate the contact between the cylinder and a concrete beam under the direct shear test. I assigned all the metal cylinder by elasto-plastic material (COMPORTEMENT = 'VMIS_ISOT_TRAC') and I got this error. The model with elastical material ran very well.
I join here the .comm file and the meshes.
Thanks.
Offline
hello,
you didn't affect traction curve on the concrete
Code_Asterの開発者
Offline
Hello,
Thank for your reply, but I want that there is only cylinder (steel) is affected by the elasto-plastical behavior, the beam in concrete still work in elastic. So how can I do that?
Thanks.
Offline
Hello,
you must add the 2 behaviors in STAT_NON_LINE:
COMPORTEMENT=(_F(GROUP_MA=('steel'),
RELATION='VMIS_ISOT_TRAC'),
_F(GROUP_MA=('concr'),
RELATION='ELAS')),
Konyaro
Last edited by konyaro (2020-07-09 05:43:39)
失敗は成功のもと (L'échec est la base de la réussite)
Offline
Oh well, I didn't paid attention on that, now it works. That great!
Thank you a lot.
Offline