Code - Aster: FORMA11 - Practical Works of The Formation "Analyzes Dynamic": Modal Analysis
Code - Aster: FORMA11 - Practical Works of The Formation "Analyzes Dynamic": Modal Analysis
Code_Aster default
Titre : FORMA11 - Travaux pratiques de la formation « Anal[...] Date : 30/07/2015 Page : 1/16
Responsable : BOYERE Emmanuel Clé : V2.08.011 Révision :
4766d086a198
Summary:
These tests correspond to practical works of the dynamic formation of Code_Aster. They are four geometries,
very made up of a linear elastic material for which one carries out a modal analysis: research of the Eigen
frequencies and the associated modes.
The statements of practical works of the dynamic formation are presented in this document:
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : FORMA11 - Travaux pratiques de la formation « Anal[...] Date : 30/07/2015 Page : 2/16
Responsable : BOYERE Emmanuel Clé : V2.08.011 Révision :
4766d086a198
1 Modeling A
1.1 Description of the problem
1.1.1 Objective
The objective of this modeling is to determine the first 11 Eigen frequencies of the fixed beam, in the
presence of multiple modes, by using the four following methods:
• Method of SORENSEN,
• Method of LANCZOS,
• Method of BATHE and WILSON
• Method of opposite iteration
1.1.2 Geometry
x
L = 0,5 m R(x)
X j
A B z
i
a = 0,05 m
y
a
The telegraphic grid could be built interactivement using Salomé. It is enough to define the points A ,
B then line AB . The size of the elements is identical. In the module of grid of Salomé, one will be
able to declare, while basing oneself on the geometry, like groups node the point A , then like meshs
the line groups AB . Then the grid will be saved with the format MED.
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : FORMA11 - Travaux pratiques de la formation « Anal[...] Date : 30/07/2015 Page : 3/16
Responsable : BOYERE Emmanuel Clé : V2.08.011 Révision :
4766d086a198
Notice : to go more quickly one can use the macro one ASSEMBLY to build the matrices!
N°1 question:
• To calculate the 11 smaller Eigen frequencies and the first associated modes (CALC_MODES
[U4.52.02]).
• To print the clean modes (IMPR_RESU) with the format MED for a visualization in Salomé.
N°2 question:
• To calculate the Eigen frequencies and the first associated modes present in the waveband
0 Hz and 6000 Hz , with the four methods SORENSEN, LANCZOS, BATHE, QZ
(CALC_MODES [U4.52.02], keyword factor SOLVEUR_MODAL=_F (METHODE=…)) then by
iterations opposite (simple keyword OPTION=' AJUSTE').
• To print the clean modes (IMPR_RESU) with the format MED for a visualization in Salomé.
right button of the mouse and to select sweep to animate the deformation.
1.4 Remarks
The good practice consists in using the option ‘BAND’ rather than ‘PLUS_PETITE’ or ‘CENTER’.
The number of frequencies is automatically given by the test of Sturm and thus all the multiplicities
are taken into account.
Method of ‘SORENSEN‘is fastest (in spite of its iterative character which allows D’ to obtain a better
precision)
One projects on a space moreover small (compared to ‘TRI_DIAG’ and with ‘JACOBI‘) and one does
not make qu’ a factorization (compared to the algorithms of the powers opposite implemented when
one uses ‘SEPARATE’ OPTION=' AJUSTE'//‘NEAR’).
In the presence of multiple modes, least modification of the parameter setting of the modal method (or
associated tripaille NUME_DDL…) product of the different clean vectors.
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : FORMA11 - Travaux pratiques de la formation « Anal[...] Date : 30/07/2015 Page : 5/16
Responsable : BOYERE Emmanuel Clé : V2.08.011 Révision :
4766d086a198
To be convinced some one can visualize the modal deformations or summon the effective modal
masses.
TRI_DIAG SORENSEN
MASS_EFFE_DY MASS_EFFE_DZ MASS_EFFE_DY MASS_EFFE_DZ
Mode 1 3.16920 2.80830 1.13915 4.83834
Mode 2 2.80672 3.17077 4.83834 1.13915
TOTAL 5.97592 5.97907 5.97749 5.97749
Frequency in Hz
Mode
(Method SORENSEN)
1 1.67638E+02
2 1.67638E+02
3 1.05060E+03
4 1.05060E+03
5 1.48054E+03
6 2.59704E+03
7 2.94237E+03
8 2.94237E+03
9 4.47822E+03
10 5.76991E+03
11 5.76991E+03
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : FORMA11 - Travaux pratiques de la formation « Anal[...] Date : 30/07/2015 Page : 6/16
Responsable : BOYERE Emmanuel Clé : V2.08.011 Révision :
4766d086a198
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : FORMA11 - Travaux pratiques de la formation « Anal[...] Date : 30/07/2015 Page : 7/16
Responsable : BOYERE Emmanuel Clé : V2.08.011 Révision :
4766d086a198
2 Modeling B
2.1 Description of the problem
2.1.1 Objective
The objective of this modeling is to determine for a structure of type sphere into free free, (presence of
multiple modes of rigid body):
• Eigen frequencies located in a waveband with the method of SORENSEN,
• Eigen frequencies located in a waveband with the method of LANCZOS with or without rigid
mode (OPTION=MODE_RIGIDE),
• 16 smaller Eigen frequencies with the method of SORENSEN
2.1.2 Geometry
Y
R = 0.01
Z X
Notice : to go more quickly one can use the macro one ASSEMBLY to build the matrices!
N°1 question:
• To calculate with the method of SORENSEN (under the keyword factor SOLVEUR_MODAL),
frequencies located in the waveband 0. Hz , 2880 Hz as well as the associated modes
(CALC_MODES).
• If calculation fails one will be able to extend the waveband on a slightly negative margin.
• To print the clean modes (IMPR_RESU) with the format MED for a visualization in Salomé.
N°2 question:
• To calculate with the method of LANCZOS with or without option MODE_RIGIDE frequencies
located in the waveband 0. Hz , 2880 Hz as well as the associated modes
(CALC_MODES).
N°3 question:
• To calculate with the method of SORENSEN 16 smaller frequencies as well as the associated
modes (CALC_MODES). One will be able to use the parameter PREC_SHIFT (under the
keyword factor CALC_FREQ) to circumvent the problem of the worthless frequencies.
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : FORMA11 - Travaux pratiques de la formation « Anal[...] Date : 30/07/2015 Page : 9/16
Responsable : BOYERE Emmanuel Clé : V2.08.011 Révision :
4766d086a198
2.4 Remarks
The factorization of the shiftées matrices K − −1= LDL T control by the parameter setting:
NMAX_ITER_SHIFT, SOLVEUR/NPREC[U4.50.01] and PREC_SHIFT, can be carried out only if those
are regular.
That often poses problem when magnitude of the terms of K is large in front of that of the terms of
M and that is one good approximate eigenvalue. The policy of the operator is to emit one ALARM
when it is about a test of Sturm and one ERREUR_FATALE when it is the matrix of work of the
operator. In the event of problem, one can always change option or of shift.
The detection of modes of rigid body is often problematic for the classical modal solveurs! The good
practice consists in using the method of SORENSEN by using a waveband whose lower limit is
worthless, even slightly negative. The solvor normally the capture without any problem (they are just
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : FORMA11 - Travaux pratiques de la formation « Anal[...] Date : 30/07/2015 Page : 10/16
Responsable : BOYERE Emmanuel Clé : V2.08.011 Révision :
4766d086a198
a little particular multiple modes). Failing this, with TRI_DIAG, the option MODE_RIGIDE is
recommended.
CALC_MODES
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : FORMA11 - Travaux pratiques de la formation « Anal[...] Date : 30/07/2015 Page : 11/16
Responsable : BOYERE Emmanuel Clé : V2.08.011 Révision :
4766d086a198
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : FORMA11 - Travaux pratiques de la formation « Anal[...] Date : 30/07/2015 Page : 12/16
Responsable : BOYERE Emmanuel Clé : V2.08.011 Révision :
4766d086a198
3 Modeling C
3.1 Description of the problem
The objective of this modeling is the use of the macro-order CALC_MODES [U4.52.02] with the option
‘BAND’ cut out in several sub-bands. This order allows to launch a succession of calculations of real
clean modes.
The modes are calculated by the order CALC_MODES [U4.52.02] with the option ‘BAND’ and
normalized by the order NORM_MODE [U4.52.11]. The calculated modes are filtered and concaténés by
means of the order EXTR_MODE [U4.52.12].
3.1.1 Objective
The objective of this modeling is:
• “To weigh” the model finite elements and to evaluate the number of modes,
• To use the macro-order CALC_MODES with the option ‘BAND’ cut out in several sub-bands to
calculate the Eigen frequencies.
3.1.2 Geometry
One will be able to recover the grid in the repertoire of the cases tests.
The thickness of the hull is of e=0,3045 m
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : FORMA11 - Travaux pratiques de la formation « Anal[...] Date : 30/07/2015 Page : 13/16
Responsable : BOYERE Emmanuel Clé : V2.08.011 Révision :
4766d086a198
Notice : to go more quickly one can use the macro one ASSEMBLY to build the matrices!
N°1 question:
• To weigh the model (POST_ELEM) and to evaluate the number of mode whose frequency is
lower than 4 Hz (INFO_MODE).
N°2 question:
• To calculate the Eigen frequencies and the first associated modes present in the waveband
0. Hz with 4 Hz (CALC_MODES)
• To standardize with the infinite standard, on all the components of the physical nodes while
requiring the calculation of the unit effective masses (NORM_MODE, EXTR_MODE).
• To print the clean modes (IMPR_RESU) with the format MED for a visualization in Salomé.
N°3 question:
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : FORMA11 - Travaux pratiques de la formation « Anal[...] Date : 30/07/2015 Page : 14/16
Responsable : BOYERE Emmanuel Clé : V2.08.011 Révision :
4766d086a198
N°4 question:
N°5 question:
• One carries out the same thing as that requested in question 4 but with the method of
LANCZOS.
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Version
Code_Aster default
Titre : FORMA11 - Travaux pratiques de la formation « Anal[...] Date : 30/07/2015 Page : 15/16
Responsable : BOYERE Emmanuel Clé : V2.08.011 Révision :
4766d086a198
3.4 Remarks
Good practices: it is necessary as a preliminary:
The order CALC_MODES, with the option ‘BAND’ cut out in several sub-bands, is more economic to
search a large spectrum.
2 2.80065 Hz
Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2017 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)