Code_Aster
Titre :
Mise en œuvre d'un calcul de modes propres d'une s[...]
Responsable :
Nicolas BRIE
Version default
Date :
14/05/2013
Page :
13/18
Clé :
U2.06.01
Révision :
11026
•Dans le fichier de commande Code_Aster :
• pour utiliser le mode basique avec INFO_MODE ou MACRO_MODE_MECA (parallélisation des sous-bandes uniquement), il n'y a rien à faire : les mot-clés par défaut activent la parallélisation des sous-bandes ;
• pour utiliser le mode avancé (parallélisation des sous-bandes pour INFO_MODE et
MACRO_MODE_MECA, et du solveur linéaire pour tous les opérateurs modaux) : utiliser le solveur linéaire MUMPS (mot-clé facteur SOLVEUR, opérande METHODE='MUMPS' ; on recommande également de paramétrer les opérandes RENUM='QAMD' et
GESTION_MEMOIRE='IN_CORE').
Exemple :
identique au paragraphe 3.2.4.1 en parallélisant à la fois les calculs sur les sous-bandes et le
solveur linéaire : modes = MACRO_MODE_MECA( MATR_RIGI = matr_k,
MATR_MASS = matr_m,
CALC_FREQ =_F( FREQ = (20., 100., 200., 300.) ),
NIVEAU_PARALLELISME = 'COMPLET',
SOLVEUR =_F( METHODE = 'MUMPS',
RENUM = 'QAMD',
GESTION_MEMOIRE = 'IN_CORE' ),
);
Remarque :
Pour des performances optimales, il est conseillé d'avoir des sous-bandes les plus équilibrées
possibles (i.e. : avec un nombre de modes recherchés par sous-bande relativement uniforme).
La mise en œuvre du parallélisme est présentée de manière plus détaillée dans la documentation générique [U2.08.06] et les documentations d'utilisation de INFO_MODE [U4.52.01] et de
MACRO_MODE_MECA [U4.52.02].
3.2.4.3 Réduction de modèle : calcul par sous-structuration
Lorsque le modèle numérique comporte un nombre élevé de degrés de liberté ou que la structure
étudiée est un assemblage de composants maillés séparément, on peut utiliser des méthodes de réduction de modèle par sous-structuration, qui reposent sur un partitionnement géométrique de la structure globale. Sur des grands modèles, ces méthodes présentent de meilleures performances CPU qu'un calcul direct.
•Sous-structuration dynamique
Cette méthode a un champ d'application très général. La documentation [U2.07.05] détaille sa mise en
œuvre.
•Sous-structuration cyclique
Cette méthode a un champ d'application beaucoup plus restrictif que la précédente : elle permet de traiter uniquement des structures à répétitivité cyclique (par exemple : roue aubagée, …). Le cas-test
SDLV301 donne un exemple de mise en œuvre.
4 Paramètres contenus dans un résultat de calcul modal
L'exécution de l'un des opérateurs de calcul modal s'accompagne de l'impression automatique de certains paramètres dans le fichier RESULTAT :
------------------------------------------------------------------------
LE NOMBRE DE DDL
Manuel d'utilisation
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.06 : Dynamique
Code_Aster
Titre :
Mise en œuvre d'un calcul de modes propres d'une s[...]
Responsable :
Nicolas BRIE
Version default
Date :
14/05/2013
Page :
14/18
Clé :
U2.06.01
Révision :
11026
TOTAL EST: 234
DE LAGRANGE EST: 120
LE NOMBRE DE DDL ACTIFS EST: 54
------------------------------------------------------------------------
L'OPTION CHOISIE EST: CENTRE
LA VALEUR DE DECALAGE EN FREQUENCE EST : 5.00000E+01
------------------------------------------------------------------------
INFORMATIONS SUR LE CALCUL DEMANDE:
NOMBRE DE MODES RECHERCHES : 1
LA DIMENSION DE L'ESPACE REDUIT EST : 0
ELLE EST INFERIEURE AU NOMBRE DE MODES, ON LA PREND EGALE A
4
=============================================
= METHODE DE SORENSEN (CODE ARPACK) =
= VERSION : 2.4 =
= DATE : 07/31/96 =
=============================================
NOMBRE DE REDEMARRAGES = 2
NOMBRE DE PRODUITS OP*X = 7
NOMBRE DE PRODUITS B*X = 20
NOMBRE DE REORTHOGONALISATIONS (ETAPE 1) = 6
NOMBRE DE REORTHOGONALISATIONS (ETAPE 2) = 0
NOMBRE DE REDEMARRAGES DU A UN V0 NUL = 0
------------------------------------------------------------------------
LES FREQUENCES CALCULEES INF. ET SUP. SONT:
FREQ_INF : 5.22037E+01
FREQ_SUP : 5.22037E+01
------------------------------------------------------------------------
CALCUL MODAL: METHODE D'ITERATION SIMULTANEE
METHODE DE SORENSEN
fréquences propres
NUMERO FREQUENCE (HZ) NORME D'ERREUR
2 5.22037E+01 7.02498E-10
3 6.74211E+01 9.12843E-10
NORME D'ERREUR MOYENNE: 0.80767E-09
position du mode dans le spectre global
------------------------------------------------------------------------
VERIFICATION A POSTERIORI DES MODES
DANS L'INTERVALLE ( 5.20730E+01, 5.23340E+01)
IL Y A BIEN 1 FREQUENCE(S)
------------------------------------------------------------------------
Si les modes calculés sont complexes, il y a en plus une colonne donnant les amortissements modaux :
Manuel d'utilisation Fascicule u2.06 : Dynamique
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Titre :
Mise en œuvre d'un calcul de modes propres d'une s[...]
Responsable :
Nicolas BRIE
Version default
Date :
14/05/2013
Page :
15/18
Clé :
U2.06.01
Révision :
11026
------------------------------------------------------------------------
LE NOMBRE DE DDL
TOTAL EST: 74
DE LAGRANGE EST: 44
LE NOMBRE DE DDL ACTIFS EST: 8
------------------------------------------------------------------------
INFORMATIONS SUR LE CALCUL DEMANDE:
NOMBRE DE MODES RECHERCHES : 5
le problème traité étant quadratique, on double l'espace de recherche
Méthode QZ dans MODE_ITER_SIMULT: On trouve un nombre de valeurs propres
17 différent du nombre de ddls physiques actifs 8 !
votre problème est fortement amorti.
valeur(s) propre(s) réelle(s) : 14
valeur(s) propre(s) complexe(s) avec conjuguée : 10
valeur(s) propre(s) complexe(s) sans conjuguée : 0
CALCUL MODAL: METHODE GLOBALE DE TYPE QR
fréquences propres (amorties)
NUMERO FREQUENCE (HZ) AMORTISSEMENT NORME D'ERREUR
1 5.52718E+00 8.68241E-03 4.00918E-13
2 1.08852E+01 1.71010E-02 7.31808E-14
3 1.59105E+01 2.50000E-02 5.40182E-14
4 2.04500E+01 3.21394E-02 4.03817E-14
5 2.43661E+01 3.83022E-02 3.48265E-14
NORME D'ERREUR MOYENNE: 0.12067E-12
amortissements modaux
Attention : pour l'instant, il n'y a pas de vérification de type STURM
(comptage du bon nombre des valeurs propres calculées)
lorsqu'on est dans le plan complexe :
problème modal généralisé avec MATR_RIGI complexe,
ou problème modal généralisé avec matrice(s) non symétrique(s),
ou problème modal quadratique (préence du mot-clé MATR_AMOR).
------------------------------------------------------------------------
VERIFICATION A POSTERIORI DES MODES
------------------------------------------------------------------------
En outre, la structure de données informatique produite lors d'un calcul modal peut contenir les paramètres suivants :
Manuel d'utilisation
Copyright 2015 EDF R&D - Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Fascicule u2.06 : Dynamique
