5. Exemple

On traite ici le cas d’une structure 3D soumise à un chargement dynamique transitoire (cf. cas test sdls114b [V2.03.114]). Après construction des matrices de masses et de rigidité, on peut calculer les modes propres de la structure:

MODE=CALC_MODES( MATR_RIGI=RIG_ASS,

MATR_MASS=MA_ASS,

CALC_FREQ=_F(NMAX_FREQ=60,),);

On peut alors calculer le déplacement de la structure soumise à un chargement dynamique:

RES_DYNA=DYNA_VIBRA(TYPE_CALCUL=”TRAN”,BASE_CALCUL=”GENE”,

MATR_MASS=MASS_GE,

MATR_RIGI=RIGI_GE,

MATR_AMOR=AMOR_GE,

SCHEMA_TEMPS=_F(SCHEMA=”EULER”,),

INCREMENT=_F(INST_INIT=0.,

INST_FIN=tfin,

VERI_PAS = “OUI”,

PAS=pas,),

EXCIT=_F(VECT_ASSE=CHA_ASS,

FONC_MULT=RAMPE,),

ARCHIVAGE=_F(PAS_ARCH = nbpas,),

IMPRESSION = _F(TOUT=”OUI”,),);

Le calcul en mécanique de la rupture commence par la définition du fond de fissure:

FF2=DEFI_FISS_XFEM(MODELE=MO,

DEFI_FISS=_F(GROUP_MA_FISS=”LEV_SUP”,

GROUP_MA_FOND=”FN_FS”, ),

GROUP_MA_ENRI=”VVOLTOT”,);

Le calcul des fonctions \({K}_{I}(t)\) , \({K}_{\mathit{II}}(t)\) et \({K}_{\mathit{III}}(t)\) se fait après le calcul préalable des facteurs d’intensité des contraintes modaux par CALC_G:

KT2 = POST_K_TRANS( RESU_TRANS = RES_DYNA,

K_MODAL = _F(TABL_K_MODA =GLMO1,

FISSURE=FF2, ),);