Example cyclical substructuring
===================================

.. code-block:: text


    PLAQUE ANNULAIRE ENCASTREE SUR A MOYEU - METHODE OF CRAIG - BAMPTON


    sector = LIRE_MAILLAGE ()

    model = AFFE_MODELE (MAILLAGE = sector,
AFFE =_F (TOUT =' OUI ',
PHENOMENE =' MECANIQUE ',
MODELISATION =' DKT '))

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

    chammat = AFFE_MATERIAU (MAILLAGE = sector,
AFFE =_F (TOUT =' OUI ', MATER = subdue))

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

    loaded = AFFE_CHAR_MECA (MODELE = model
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 FROM RIGIDITE AND FROM MASSE FROM SECTEUR FROM BASE

    #

    rigid = CALC_MATR_ELEM (MODELE = model, CHARGE = load,
CHAM_MATER = chammat, CARA_ELEM = chamcar,
OPTION = 'RIGI_MECA')

    mass = CALC_MATR_ELEM (MODELE = model, CHARGE = load,

    CHAM_MATER = chammat, CARA_ELEM = chamcar,

    OPTION =' MASS_MECA ')

    number = NUME_DDL (MATR_RIGI = rigid)

    matrigi = ASSE_MATRICE (MATR_ELEM = rigid, NUME_DDL = number)

    matmass = ASSE_MATRICE (MATR_ELEM = massed, NUME_DDL = number)

    #

    # CALCUL DES MODES DYNAMIQUES FROM SECTEUR FROM BASE

    #

    fashion = CALC_MODES (MATR_RIGI = matrix,

    MATR_MASS = matmass,

    CALC_FREQ = _F (NMAX_FREQ = 15)

    #

    # DEFINITION DES INTERFACES AND DES MODES STATIQUES ASSOCIES

    #

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

    #

    # CALCUL OF THE BASE OF PROJECTION = RECUPERATION DES MODES DYNAMIQUES

    # AND CALCUL DES MODES STATIQUES

    bamo = DEFI_BASE_MODALE (CLASSIQUE = _F (INTERF_DYNA = lint, IMPR = 2,
MODE_MECA = fashion,
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))