2D Tria3 DKT-Modell mit nichtlinearem Materialverhalten - Zwischenschrittresultate möglich?

Hallo,

habe einen Hohlkörper (Tria3 DKT-Elemente) mit Innendruck nichtlinear (Materialnichtlinearität) am laufen. Kann ich in Code Aster eine Einstellung tätigen, wo ich Mises Zwischenergebnisse herausschreiben kann, auch wenn mein Innendruck zu hoch ist und keine Konvegenz des Zieldrucks erreicht wird? Sind die DKT für große Verformungen überhaupt geeignet? Coque_3D sind es ja nicht, oder?

Ist der folgende Code überhaupt sinnvoll?

[code:1]
# leider von Anfänger und deshalb fehlerhaft !!!!!
DEBUT();

Traction=DEFI_FONCTION(NOM_PARA='EPSI',VALE=(0.003143,220.0,
0.005199,290.0,
0.010668,330.0,
0.022036,365.0,
0.041728,405.0,
0.072294,440.0,
0.116391,470.0,
0.176765,510.0,
0.256243,545.0,
0.357725,585.0,
0.484174,620.0,
0.638615,650.0,
0.824129,660.0,
),INTERPOL='LIN',PROL_DROITE='CONSTANT',PROL_GAUCHE='EXCLU',);

Stahl=DEFI_MATERIAU(ELAS=_F(E=210000.0,
NU=0.3,),
TRACTION=_F(SIGM=Traction,),);

MAIL=LIRE_MAILLAGE(UNITE=20,
FORMAT='MED',
NOM_MED='ausgleic',
INFO=1,);

MAIL=MODI_MAILLAGE(reuse =MAIL,
MAILLAGE=MAIL,
ORIE_LIGNE=_F(GROUP_MA='EleTria3',),);

MatMod=AFFE_MATERIAU(MAILLAGE=MAIL,
AFFE=_F(TOUT='OUI',
MATER=Stahl,),);

EleMod=AFFE_MODELE(MAILLAGE=MAIL,
INFO=2,
AFFE=_F(GROUP_MA='EleTria3',
PHENOMENE='MECANIQUE',
MODELISATION='DKT',),);

EleProp=AFFE_CARA_ELEM(MODELE=EleMod,
COQUE=_F(GROUP_MA='EleTria3',
EPAIS=5.0,
VECTEUR=(12.0,19.0,27.0,),),);

Ground=AFFE_CHAR_MECA(VERI_NORM='OUI',
MODELE=EleMod,
DDL_IMPO=_F(GROUP_MA='LagerEle',
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='EleTria3',
PRES=-1,),);

Time=DEFI_LIST_REEL(DEBUT=0.0,
INTERVALLE=_F(JUSQU_A=1.0,
NOMBRE=2,),);

ramp=DEFI_FONCTION(NOM_PARA='INST',VALE=(0.0,0.0,
1.0,1.0,
),);

autostep=DEFI_LIST_INST(DEFI_LIST=_F(METHODE='AUTO',
LIST_INST=Time,
PAS_MINI=1.E-03,),
ECHEC=_F(EVENEMENT='DIVERGENCE',
SUBD_METHODE='UNIFORME',
SUBD_COEF_PAS_1=0.5,
SUBD_PAS=8,),
ADAPTATION=_F(EVENEMENT='SEUIL',
NB_INCR_SEUIL=2,
NOM_PARA='NB_ITER_NEWTON',
CRIT_COMP='LE',
MODE_CALCUL_TPLUS='FIXE',
PCENT_AUGM=100.0,),);

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',
RESI_DEBO_RELA=1e-03,
ITER_INTE_PAS=0,
GROUP_MA='EleTria3',),
INCREMENT=_F(LIST_INST=autostep,
INST_FIN=1.0,),
NEWTON=_F(REAC_INCR=1,
MATRICE='TANGENTE',
REAC_ITER=1,),);

Solution=CALC_ELEM(reuse =Solution,
MODELE=EleMod,
RESULTAT=Solution,
OPTION=('SIEF_ELNO_ELGA','EPSG_ELNO_DEPL','EPSI_ELNO_DEPL',),);

IMPR_RESU(MODELE=EleMod,
FORMAT='MED',
RESU=_F(MAILLAGE=MAIL,
RESULTAT=Solution,
INFO_MAILLAGE='OUI',
PARTIE='REEL',
TOUT_CHAM='OUI',
TOUT_CMP='OUI',
TOUT='OUI',),
INFO=2,);
FIN(FORMAT_HDF='OUI',);

[/code:1]

Gruß und Dank Markus<br /><br />Post edited by: MGolbs, at: 2011/01/13 17:29