nichtlineare Materialkennlinie wird ignoriert

Hallo,

habe folgenden Code. Dabei wird die Materialkennlinie ignoriert, warum nur?

[code:1]
DEBUT();
MAILIMPO=LIRE_MAILLAGE(UNITE=20,
FORMAT='MED',
NOM_MED='H20Ausgl',
INFO=1,);
MAIL=CREA_MAILLAGE(MAILLAGE=MAILIMPO,
MODI_MAILLE=_F(GROUP_MA='alleElem',
OPTION='TRIA6_7',
PREF_NOEUD='NSS',),);
MAIL=MODI_MAILLAGE(reuse =MAIL,
MAILLAGE=MAIL,
ORIE_LIGNE=_F(GROUP_MA='alleElem',),);
T10012=DEFI_FONCTION(NOM_PARA='EPSI',VALE=(0.07,14.0,
0.2,14.2,
0.8,14.25,
2.0,14.27,
50,14, ),INTERPOL='LIN',PROL_DROITE='LINEAIRE',PROL_GAUCHE='EXCLU',);
12T100=DEFI_MATERIAU(ELAS=_F(E=200.0,
NU=0.4,
RHO=1000000000.0,),
TRACTION=_F(SIGM=T10012,),);
MatMod=AFFE_MATERIAU(MAILLAGE=MAIL,
AFFE=_F(TOUT='OUI',
MATER=12T100,),);
EleMod=AFFE_MODELE(MAILLAGE=MAIL,
INFO=2,
AFFE=_F(GROUP_MA='alleElem',
PHENOMENE='MECANIQUE',
MODELISATION='COQUE_3D',),);
EleProp=AFFE_CARA_ELEM(MODELE=EleMod,
COQUE=_F(GROUP_MA='alleElem',
EPAIS=1,
ANGL_REP=(0.,0.,),
COQUE_NCOU=5,),);
Ground=AFFE_CHAR_MECA(VERI_NORM='OUI',
MODELE=EleMod,
DDL_IMPO=_F(GROUP_MA='LagFace',
DX=0.0,
DY=0.0,
DZ=0.0,
DRX=0.0,
DRY=0.0,
DRZ=0.0,),);
Pilot=AFFE_CHAR_MECA(VERI_NORM='OUI',
MODELE=EleMod,
PRES_REP=_F(GROUP_MA='alleElem',
PRES=-1.0,),);
dt = 0.1;
t0 = 0.0;
t1 = 1.0;
tsteps = 10;
ramp=DEFI_FONCTION(NOM_PARA='INST',VALE=(t0,0.00,
t1,1.00,
),INFO=2,TITRE='ramp',);
schritte=DEFI_LIST_REEL(DEBUT=0.0,
INTERVALLE=_F(JUSQU_A=t1,
NOMBRE=tsteps,),
INFO=2,
TITRE='time',);
Solution=STAT_NON_LINE(MODELE=EleMod,
CHAM_MATER=MatMod,
CARA_ELEM=EleProp,
EXCIT=(_F(CHARGE=Ground,),
_F(CHARGE=Pilot,
FONC_MULT=ramp,),),
COMP_INCR=_F(RELATION='ELAS',
DEFORMATION='GREEN',),
INCREMENT=_F(LIST_INST=schritte,),
CONVERGENCE=_F(RESI_GLOB_MAXI=1e-8,),
NEWTON=_F(REAC_ITER=1,),);
Solution=CALC_ELEM(reuse =Solution,
MODELE=EleMod,
RESULTAT=Solution,
REPE_COQUE=_F(NIVE_COUCHE='MOY',),
TYPE_OPTION='SIGM_STRUCT',
OPTION=('SIGM_ELNO_COQU','SIEF_ELNO_ELGA',),);
Solution=CALC_NO(reuse =Solution,
RESULTAT=Solution,
OPTION=('FORC_NODA','REAC_NODA','SIGM_NOEU_COQU',),);
IMPR_RESU(MODELE=EleMod,
FORMAT='MED',
UNITE=80,
RESU=_F(MAILLAGE=MAIL,
RESULTAT=Solution,
INFO_MAILLAGE='OUI',),
INFO=1,);
FIN(FORMAT_HDF='OUI',);
[/code:1]

Über Tipps und Infos würde ich mich freuen.

Gruß Markus<br /><br />Post edited by: MGolbs, at: 2011/03/18 14:41

Hallo Claus,

danke für den guten Tipp. Jetzt scheint es zu laufen. Nur ist diese Kombination sinnvoll?

[code:1]
DEBUT();
MAILIMPO=LIRE_MAILLAGE(UNITE=20,
FORMAT='MED',
NOM_MED='H20Ausgl',
INFO=1,);
MAIL=CREA_MAILLAGE(MAILLAGE=MAILIMPO,
MODI_MAILLE=_F(GROUP_MA='alleElem',
OPTION='TRIA6_7',
PREF_NOEUD='NSS',),);
MAIL=MODI_MAILLAGE(reuse =MAIL,
MAILLAGE=MAIL,
ORIE_LIGNE=_F(GROUP_MA='alleElem',),);
T10012=DEFI_FONCTION(NOM_PARA='EPSI',VALE=(0.07,14.0,
0.2,14.2,
0.8,14.25,
2.0,14.27,
50.0,14.0,
),INTERPOL='LIN',PROL_DROITE='LINEAIRE',PROL_GAUCHE='EXCLU',);
12T100=DEFI_MATERIAU(ELAS=_F(E=200.0,
NU=0.4,
RHO=1000000000.0,),
TRACTION=_F(SIGM=T10012,),);
MatMod=AFFE_MATERIAU(MAILLAGE=MAIL,
AFFE=_F(TOUT='OUI',
MATER=12T100,),);
EleMod=AFFE_MODELE(MAILLAGE=MAIL,
INFO=2,
AFFE=_F(GROUP_MA='alleElem',
PHENOMENE='MECANIQUE',
MODELISATION='COQUE_3D',),);
EleProp=AFFE_CARA_ELEM(MODELE=EleMod,
COQUE=_F(GROUP_MA='alleElem',
EPAIS=5.0,
VECTEUR=(12.0,19.0,27.0,),
COQUE_NCOU=5,),);
Ground=AFFE_CHAR_MECA(VERI_NORM='OUI',
MODELE=EleMod,
DDL_IMPO=_F(GROUP_MA='LagFace',
DX=0.0,
DY=0.0,
DZ=0.0,
DRX=0.0,
DRY=0.0,
DRZ=0.0,),);
Pilot=AFFE_CHAR_MECA(VERI_NORM='OUI',
MODELE=EleMod,
PRES_REP=_F(GROUP_MA='alleElem',
PRES=-1,),);
dt = 0.025;
t0 = 0.0;
t1 = 1.0;
tsteps = 40;
ramp=DEFI_FONCTION(NOM_PARA='INST',VALE=(t0,0.00,
t1,1.00,
),INFO=2,TITRE='ramp',);
schritte=DEFI_LIST_REEL(DEBUT=0.0,
INTERVALLE=_F(JUSQU_A=t1,
NOMBRE=tsteps,),
INFO=2,
TITRE='time',);
Solution=STAT_NON_LINE(MODELE=EleMod,
CHAM_MATER=MatMod,
CARA_ELEM=EleProp,
EXCIT=(_F(CHARGE=Ground,),
_F(CHARGE=Pilot,
FONC_MULT=ramp,),),
COMP_INCR=_F(RELATION='VMIS_ISOT_TRAC',
DEFORMATION='GREEN',
TOUT='OUI',),
INCREMENT=_F(LIST_INST=schritte,),
CONVERGENCE=_F(RESI_GLOB_MAXI=1e-8,),
NEWTON=_F(REAC_ITER=1,),);
Solution=CALC_ELEM(reuse =Solution,
MODELE=EleMod,
RESULTAT=Solution,
REPE_COQUE=_F(NIVE_COUCHE='MOY',),
TYPE_OPTION='SIGM_STRUCT',
OPTION=('SIGM_ELNO_COQU','SIEF_ELNO_ELGA',),);
Solution=CALC_NO(reuse =Solution,
RESULTAT=Solution,
OPTION=('FORC_NODA','REAC_NODA','SIGM_NOEU_COQU',),);
IMPR_RESU(MODELE=EleMod,
FORMAT='MED',
UNITE=80,
RESU=_F(MAILLAGE=MAIL,
RESULTAT=Solution,
INFO_MAILLAGE='OUI',),
INFO=1,);
FIN(FORMAT_HDF='OUI',);
[/code:1]

Gruß und Dank Markus

Wenn Sie haben ein Problem mit der Konvergenz, vielleicht die Traktion Kurve muss immer erhöhen. Löschen Sie die Paar (50.0,14.0) oder versuchen (50.0,14.28) wenn es hilft.

Mit freundlichen Grüßen,

Tim Donley<br /><br />Post edited by: T. Donley, at: 2011/03/19 02:03

Hallo,

danke für den Tipp. Habe nun eine steigende Materialkurve drin (10.0, 14.3). Das Modell bricht aber um den Faktor 7 zu früh ab. Dieser Faktor kommt aus dem Abgleich Rechnung Versuch. Gibt es das Coque_3D Element auch als lineares Element? Besonders die Biegeschubspannungen sind in dem Modell wichtig. Deshalb habe ich mich für Coque_3D entscheiden. Kann ich mit anderen Elementen z.B. DKT auch in Schichten rechnen? Bekomme ich mit anderen 2D Elementen diese Biegeschubspannungen auch gut nichtlinear abgebildet?

Gruß Markus