Code_Aster ®
Version
7.4
Titre :
Opérateur MODE_ITER_CYCL
Date :
08/03/05
Auteur(s) :
O. NICOLAS, E. BOYERE Clé
:
U4.52.05-G Page
: 1/6
Organisme(s) : EDF-R&D/AMA
Manuel d'Utilisation
Fascicule U4.5- : Méthodes de résolution
Document : U4.52.05
Opérateur MODE_ITER_CYCL
1 But
Calculer les modes propres d'une structure à symétrie cyclique.
On calcule les composantes généralisées des modes propres de la structure entière, par une méthode
de sous-structuration cyclique, à partir de la base modale d'un secteur de référence. L'axe de symétrie
est l'axe OZ. La base modale du secteur doit être de type CLASSIQUE. Les interfaces DROITE,
GAUCHE et éventuellement AXE doivent être de même type. Les côtés droit et gauche sont définis par
le sens trigonométrique dans le plan OXY.
Produit une structure de données de type mode_cycl.
Manuel d'Utilisation
Fascicule U4.5- : Méthodes de résolution
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur MODE_ITER_CYCL
Date :
08/03/05
Auteur(s) :
O. NICOLAS, E. BOYERE Clé
:
U4.52.05-G Page
: 2/6
2 Syntaxe
mocy [mode_cycl] = MODE_ITER_CYCL
(
BASE_MODALE
=
bamo,
[base_modale]
NB_MODE =
/
nbmo,
[I]
/ 999
,
[DEFAUT]
NB_SECTEUR
=
nbsec,
[I]
LIAISON = _F(
DROITE
=
'nom_int',
[Kn]
GAUCHE
=
'nom_int',
[Kn]
AXE
=
'nom_int', [Kn]
),
CALCUL = _F(
TOUT_DIAM
= / 'OUI',
NB_DIAM = li, [l_I]
OPTION = / 'PLUS_PETITE',
[DEFAUT]
/
'CENTRE',
/
'BANDE',
FREQ = lifreq,
[l_R]
NMAX_FREQ = / nbfreq,
[I]
/
10,
[DEFAUT]
PREC_SEPARE =
/
pre_sep,
[R]
/
1.E+2,
[DEFAUT]
PREC_AJUSTE =
/
pre_ajus,
[R]
/
1.E-6,
[DEFAUT]
NMAX_ITER
=
/
niter,
[I]
/
50,
[DEFAUT]
),
VERI_CYCL = _F (
PRECISION
= / prec, [R]
/
1.D-3,
[DEFAUT]
CRITERE =
/ 'RELATIF',[DEFAUT]
DIST_REFE =
dist_ref, [R]
),
INFO =
/
1,
[DEFAUT]
/ 2,
)
Manuel d'Utilisation
Fascicule U4.5- : Méthodes de résolution
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur MODE_ITER_CYCL
Date :
08/03/05
Auteur(s) :
O. NICOLAS, E. BOYERE Clé
:
U4.52.05-G Page
: 3/6
3 Opérandes
3.1 Opérande
BASE_MODALE
BASE_MODALE = bamo
Nom de la base modale du secteur construite par DEFI_BASE_MODALE [U4.64.02].
3.2 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.
3.3 Opérande
NB_SECTEUR
NB_SECTEUR = nbsec
Nombre de secteurs de base nécessaires à la construction de la structure globale.
3.4 Mot
clé
LIAISON
LIAISON
Mot clé facteur pour la définition des liaisons entre les secteurs.
3.4.1 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.
3.5 Mot
clé
CALCUL
CALCUL
Mot clé facteur pour la définir le mode de recherche des modes propres.
3.5.1 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.
Manuel d'Utilisation
Fascicule U4.5- : Méthodes de résolution
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur MODE_ITER_CYCL
Date :
08/03/05
Auteur(s) :
O. NICOLAS, E. BOYERE Clé
:
U4.52.05-G Page
: 4/6
3.5.2 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.
3.5.3 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 ( f f
1
2) 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).
3.5.4 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).
Manuel d'Utilisation
Fascicule U4.5- : Méthodes de résolution
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur MODE_ITER_CYCL
Date :
08/03/05
Auteur(s) :
O. NICOLAS, E. BOYERE Clé
:
U4.52.05-G Page
: 5/6
3.6 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.
Interface gauche
Y
Interface axe
Z
X
Interface droite
Figure 3.6-a
3.6.1 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.
3.7 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.
Manuel d'Utilisation
Fascicule U4.5- : Méthodes de résolution
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur MODE_ITER_CYCL
Date :
08/03/05
Auteur(s) :
O. NICOLAS, E. BOYERE Clé
:
U4.52.05-G Page
: 6/6
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 = MODE_ITER_SIMULT
(MATR_A = matrigi, MATR_B = 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 ))
Manuel d'Utilisation
Fascicule U4.5- : Méthodes de résolution
HT-66/05/004/A
Document Outline