Exemple sous-structuration cyclique
===================================

.. code-block:: text


    PLAQUE ANNULAIRE ENCASTREE SUR UN MOYEU - METHODE DE CRAIG-BAMPTON


    secteur = LIRE_MAILLAGE        (    )

    modele = AFFE_MODELE             (    MAILLAGE= secteur,
AFFE =_F(    TOUT    ='OUI',
PHENOMENE ='MECANIQUE',
MODELISATION='DKT')        )

    mater = DEFI_MATERIAU        (ELAS =_F(E=2.E11, NU=0.3, RHO=7800.0)        )

    chammat = AFFE_MATERIAU        (MAILLAGE= secteur,
AFFE =_F(TOUT ='OUI', MATER= mater)        )

    chamcar = AFFE_CARA_ELEM    (MODELE = modele,
COQUE =(TOUT ='OUI', EPAIS= 0.001)    )

    charge = AFFE_CHAR_MECA    (MODELE = modele
DDL_IMPO=(TOUT='OUI',DX=0.,DY=0.,DRZ=0.),
DDL_IMPO=(GROUP_NO=’AXE’,DZ=0.,DRX=0.,DRY=0.),
DDL_IMPO=(GROUP_NO=’DROIT’,DZ=0.,DRX=0.,DRY=0.),
DDL_IMPO=(GROUP_NO=’GAUCH’,DZ=0.,DRX=0.,DRY=0.))

    #

    #         CONSTRUCTION DES MATRICES DE RIGIDITE ET DE MASSE DU SECTEUR DE BASE

    #

    rigiele = CALC_MATR_ELEM     (MODELE = modèle,     CHARGE = charge,
CHAM_MATER= chammat,    CARA_ELEM = chamcar,
OPTION = 'RIGI_MECA'        )

    massele = CALC_MATR_ELEM     (MODELE = modele,    CHARGE = charge,

    CHAM_MATER= chammat,    CARA_ELEM = chamcar,

    OPTION ='MASS_MECA'        )

    numerot = NUME_DDL            (MATR_RIGI = rigiele            )

    matrigi = ASSE_MATRICE        (MATR_ELEM = rigiele, NUME_DDL = numerot        )

    matmass = ASSE_MATRICE        (MATR_ELEM = massele, NUME_DDL = numerot        )

    #

    #                     CALCUL DES MODES DYNAMIQUES DU SECTEUR DE BASE

    #

    modes = CALC_MODES (MATR_RIGI = matrigi,

    MATR_MASS = matmass,

    CALC_FREQ= _F(NMAX_FREQ= 15)    )

    #

    #         DEFINITION DES INTERFACES ET DES MODES STATIQUES ASSOCIES

    #

    lint = DEFI_INTERF_DYNA    (NUME_DDL = numerot, IMPR= 2,
INTERFACE= _F(NOM='DROITE', TYPE='CRAIGB',
GROUP_NO= ’DROIT’,
MASQUE= ('DX', 'DY', 'DRZ'), ),
INTERFACE= _F(NOM='GAUCHE', TYPE='CRAIGB',
GROUP_NO= ’GAUCH’,
MASQUE= ('DX', 'DY', 'DRZ') ) )

    #

    #         CALCUL DE LA BASE DE PROJECTION = RECUPERATION DES MODES DYNAMIQUES

    # ET CALCUL DES MODES STATIQUES

    bamo = DEFI_BASE_MODALE    (CLASSIQUE= _F(INTERF_DYNA= lint, IMPR= 2,
MODE_MECA = modes,
NMAX_MODE= 15 ) )

    #

    #                                 CALCUL DES MODES CYCLIQUES

    #

    modcyc = MODE_ITER_CYCL (BASE_MODALE= bamo, NB_MODE=15, NB_SECTEUR=18,

    LIAISON=_F(DROITE= 'DROITE',    GAUCHE= 'GAUCHE'),

    CALCUL =_F(NB_DIAM=(0, 1, 2, 3), NMAX_FREQ=2 ))