4. Exemple

4.1. OPTION=”MECA”

Voir les tests COMP001 [V6,07,101] , COMP002 [V6,07,102]

#unités en Pa

ACIER[0]=DEFI_MATERIAU(ELAS=_F(E=YOUNG_Pa,NU=POISSON,ALPHA=11.8e-6),ECRO_LINE=_F(D_SIGM_EPSI=pente_Pa,SY=SY_Pa,),);#unités en MPaACIER[1]=DEFI_MATERIAU(ELAS=_F(E=YOUNG,NU=POISSON,ALPHA=11.8e-6),ECRO_LINE=_F(D_SIGM_EPSI=pente,SY=SY,),)compor='VMIS_ISOT_LINE'tabresu=TEST_COMPOR(OPTION='MECA',COMPORTEMENT=_F(RELATION=compor,),NEWTON=_F(REAC_ITER=1),LIST_MATER=ACIER,VARI_TEST=('V1','VMIS','TRACE'),YOUNG=YOUNG,POISSON=POISSON,)

4.2. OPTION=”THER”

Voir les tests COMP008*

TREF = 0.

Tmax = 500.

YOUN=DEFI_FONCTION(NOM_PARA='TEMP',VALE=(TREF, 200000.,

Tmax, 100000.,

),);

ALPH=DEFI_FONCTION(NOM_PARA='TEMP',VALE=(TREF, 1.E-5,

Tmax, 2.E-5,

),);

SIGY=DEFI_FONCTION(NOM_PARA='TEMP',VALE=(TREF, 100.,

Tmax, 50.,

),);

DSDE=DEFI_FONCTION(NOM_PARA='TEMP',VALE=(TREF, 10000.,

Tmax, 5000.,

),);

MATERI=DEFI_MATERIAU(ELAS_FO=_F(E=YOUN, NU=ZERO,

TEMP_DEF_ALPHA=TREF,

ALPHA=ALPH,),

ECRO_LINE_FO=_F(D_SIGM_EPSI=DSDE,

SY=SIGY,),

);

LMAT2 = [None]*(NCAL)

time=0.

for i in range(NCAL):

timem = time

time = timem + tfin/NCAL

Ti = T0 + time /tfin * (Tmax - T0)

LMAT2[i]=DEFI_MATERIAU(ELAS =_F(E=YOUN(Ti),

NU=0.,

ALPHA=0.,),

ECRO_LINE=_F(D_SIGM_EPSI=DSDE(Ti),

SY=SIGY(Ti),),)

compor='VMIS_ISOT_LINE'


tabresu=TEST_COMPOR(

MATER=MATERI,COMPORTEMENT=_F(RELATION=compor),

LIST_MATER=LMAT2,ALPHA=ALPH,YOUNG=YOUN,TEMP_INIT=TREF,TEMP_FIN=Tmax,

NEWTON=_F(REAC_ITER=1),

NB_VARI=2, VARI_TEST=('V1','V2'),

)