Exemple
=======

On traite ici le cas d’une structure 3D soumise à un chargement dynamique transitoire (cf. cas test sdls114b [:external:ref:`V2.03.114 <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:


.. code-block:: text
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 :math:`{K}_{I}(t)` , :math:`{K}_{\mathit{II}}(t)` et :math:`{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, ),);