u4.52.05 Opérateur MODE_ITER_CYCL#

Syntaxe#

Détail de la syntaxe
mode_cycl = MODE_ITER_CYCL(
    ◆ BASE_MODALE = mode_meca,
    ◇ NB_MODE = int,
    ◆ NB_SECTEUR = int,
    ◆ LIAISON = _F(
         ◆ DROITE = text,
         ◆ GAUCHE = text,
         ◇ AXE = text,
      ),
    ◇ VERI_CYCL = _F(
         ◇ PRECISION = float (défaut: 0.001),
         ◇ CRITERE = "RELATIF",
         ◇ DIST_REFE = float,
      ),
    ◆ CALCUL = _F(
         ◆ / TOUT_DIAM = "OUI" (ou non renseigné),
           / NB_DIAM = int,
         ◇ OPTION = / "PLUS_PETITE" (par défaut),
                    / "CENTRE",
                    / "BANDE",
         # Si: equal_to("OPTION", 'CENTRE')
             ◆ FREQ = float,
         # Si: equal_to("OPTION", 'BANDE')
             ◆ FREQ = float,
         ◇ NMAX_FREQ = int (défaut: 10),
         ◇ PREC_SEPARE = float (défaut: 100.0),
         ◇ PREC_AJUSTE = float (défaut: 1e-06),
         ◇ NMAX_ITER = int (défaut: 50),
      ),
    ◇ INFO = / 1 (par défaut),
             / 2,
)


◆ : obligatoire
◇ : optionnel
⟐ : présent par défaut
& : ensemble
/ : un seul parmi
| : plusieurs choix possibles

Opérandes#

Opérande BASE_MODALE#

♦    BASE_MODALE= bamo

Nom de la base modale du secteur construite par DEFI_BASE_MODALE [U4.64.02].

Opérande NB_MODE#

◊    NB_MODE= nbmo

Nombre de modes propres du secteur à utiliser pour le calcul cyclique. Par défaut, si le mot clé n’apparaît pas, tous les modes propres de la base modale sont utilisés.

Opérande NB_SECTEUR#

♦    NB_SECTEUR= nbsec

Nombre de secteurs de base nécessaires à la construction de la structure globale.

Mot clé LIAISON#

♦    LIAISON

Mot clé facteur pour la définition des liaisons entre les secteurs.

Opérandes DROITE / GAUCHE / AXE#

Voir [Figure 3.6-a].


♦ DROITE = ’nom_int’

Nom de l’interface droite du secteur.


♦ GAUCHE = ’nom_int’

Nom de l’interface gauche du secteur.


◊ AXE = ’nom_int’

Nom de l’interface de l’axe du secteur.

Ce sont des points communs à tous les secteurs.

Mot clé CALCUL#

♦    CALCUL

Mot clé facteur pour définir le mode de recherche des modes propres.

Opérandes TOUT_DIAM / NB_DIAM#


◊ TOUT_DIAM = “OUI”

Les modes associés à tous les nombres de diamètres nodaux seront calculés.


◊ NB_DIAM = li

Liste des nombres de diamètres nodaux à calculer. Par défaut, tous les nombres de diamètres nodaux possibles sont étudiés.

Opérande OPTION#


◊ OPTION=

“PLUS_PETITE”:

calculer par une méthode d’itération inverse les modes propres correspondant aux plus petites fréquences pour chaque nombre de diamètres demandés.

“CENTRE”:

calculer les modes propres centrés autour d’une fréquence demandée par le mot clé LIST_FREQ.

“BANDE”:

calculer les modes propres entre deux fréquences données par l’utilisateur par le mot clé LIST_FREQ. Les fréquences propres sont séparées par dichotomie puis les modes propres calculés par itérations inverses centrées sur les fréquences issues de l’étape de séparation.

Opérandes FREQ / NMAX_FREQ#


◊ FREQ = lifreq

Liste des fréquences dont l’utilisation dépend de l’option choisie :


OPTION = “BANDE”

On attend 2 valeurs

../../../../_images/Object_1107.svg

qui définissent la bande.


OPTION = “CENTRE”

On attend 1 valeur qui est la fréquence centrale de l’intervalle.


OPTION = “PLUS_PETITE”

On calcule les plus petites fréquences propres de la structure. Par défaut, on calcule les 10 premières. Le mot clé FREQ n’a alors pas de sens dans ce cas, il n’a pas à être renseigné.


◊ NMAX_FREQ = nbfreq

Nombre de fréquences à calculer pour chaque nombre de diamètres nodaux demandé. Si ce mot clé n’apparaît pas, on calcule autant de fréquences, pour chaque diamètre nodal, qu’il y a de modes propres utilisés dans la base modale (mot clé NB_MODE).

Opérandes PREC_SEPARE / PREC_AJUSTE / NMAX_ITER#


◊ PREC_SEPARE = pre_sep

Précision de séparation des fréquences pour option “BANDE”.


◊ PREC_AJUSTE = pre_ajus

Précision utilisée pour le calcul des modes (toutes OPTIONS).


◊ NMAX_ITER = niter

Nombre maximum d’itérations inverses (toutes OPTIONS).

Mot clé VERI_CYCL#

♦    VERI_CYCL

Mot clé pour vérification de la cohérence des interfaces données en terme de répétitivité cyclique.

../../../../_images/1000096A00002C2200000F9B40247FADF9EC8BF2.svg

Figure 3.6-a

Opérandes PRECISION / DIST_REFE#

◊    PRECISION    =    prec

◊    DIST_REFE    =    dist_ref

Le test de cohérence entre 2 secteurs contigus sera déterminé par le produit prec*dist_ref. Si DIST_REFE n’est pas renseigné, il sera automatiquement calculé proportionnellement à prec et à une valeur maximale de coordonnée d’un secteur.

Opérande INFO#

◊    INFO    =

Niveau d’impression

1

pas d’impression,

2

écriture des fréquences et paramètres généralisés obtenus et des participations relatives des différents modes de la base.

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