4. Exemple sous-structuration cyclique

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 ))