Code_Aster ®
Version
7.4
Titrate:
Operator CALC_CHAR_SEISME
Date:
08/02/05
Author (S):
Y. PONS, Fe WAECKEL, L. VIVAN
Key: U4.63.01-G Page:
1/6
Organization (S): EDF-R & D/AMA, SINETICS, CS IF
Handbook of Utilization
U4.6- booklet: Stamp/Vecteurs elementary and assembly
Document: U4.63.01
Operator CALC_CHAR_SEISME
1 Goal
To establish the seismic loading for a calculation of response moving relative compared to
supports. For a better theoretical comprehension of this loading, one will refer to the document
[R4.05.01].
The produced concept is directly usable during a direct transitory analysis with
DYNA_LINE_TRAN [U4.53.02] or by modal synthesis with DYNA_TRAN_MODAL [U4.53.21]. In
revenge for a nonlinear direct transitory analysis with DYNA_NON_LINE [U4.53.01], it is necessary
to transform this concept charges some starting from operator AFFE_CHAR_MECA [U4.44.01].
Product an assembled vector, concept of the cham_no_depl_R. type.
Handbook of Utilization
U4.6- booklet: Stamp/Vecteurs elementary and assembly
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titrate:
Operator CALC_CHAR_SEISME
Date:
08/02/05
Author (S):
Y. PONS, Fe WAECKEL, L. VIVAN
Key: U4.63.01-G Page:
2/6
2 Syntax
S [cham_no_depl_R] = CALC_CHAR_SEISME
(MATR_MASS = m
[matr_asse_DEPL_R]
DIRECTION = (d1, d2, d3, r1, r2, r3)
[l_R]
/MONO_APPUI = “YES”
/
MODE_STAT = mode [mode_stat]
/NOEUD
=
noeu
[l_noeud]
/
GROUP_NO
=
g_noeu [l_gr_noeud]
TITER = title [l_Kn]
);
Handbook of Utilization
U4.6- booklet: Stamp/Vecteurs elementary and assembly
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titrate:
Operator CALC_CHAR_SEISME
Date:
08/02/05
Author (S):
Y. PONS, Fe WAECKEL, L. VIVAN
Key: U4.63.01-G Page:
3/6
3 Operands
3.1 Operand
MATR_MASS
MATR_MASS = m
Stamp of mass of the system.
3.2 Operand
DIRECTION
DIRECTION = (d1, d2, d3, r1, r2, r3)
Components of a vector giving the direction of the seism in the total reference mark. It is a list of
three realities if the imposed accélérogrammes are only translations. If one imposes
also accelerations of rotations, one awaits a list of six realities (valid for
modelings with discrete elements).
3.3
Description of the movement of drive
3.3.1 Operand
MONO_APPUI
/MONO_APPUI = “YES”
The structure is uniformly excited in all the supports (movement of drive of body
solid).
3.3.2 Excitation multi supports
In this case, the accelerations undergone by the whole of the points of anchoring of the studied structure
are not inevitably identical and in phase.
3.3.2.1 Operand
MODE_STAT
/
MODE_STAT = mode
Static modes of the structure: concept of the mode_stat type produces by the operator
MODE_STATIQUE [U4.52.14] with option MODE_STAT. They correspond to the 6
nb_supports static modes where nb_supports is the number of accélérogrammes
different undergone by the structure.
Note:
If the structure is requested only by translations, there are then 3 nb_supports
static modes.
3.3.2.2 Operands
NODE/GROUP_NO
/NOEUD = noeu
/
GROUP_NO =
g_noeu
List nodes (noeu) or groups of nodes (g_noeu) of the structure subjected to
the seismic excitation: these nodes support the ddl supports of the structure
which are applied the imposed movements.
3.4 Operand
TITER
TITER
= title
Titrate attached to the concept produced by this operator [U4.03.01].
Handbook of Utilization
U4.6- booklet: Stamp/Vecteurs elementary and assembly
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titrate:
Operator CALC_CHAR_SEISME
Date:
08/02/05
Author (S):
Y. PONS, Fe WAECKEL, L. VIVAN
Key: U4.63.01-G Page:
4/6
4 Examples
The two examples which follow illustrate the employment of operator CALC_CHAR_SEISME in the cases
of a seismic excitation mono support and in the case multi supports (identical acceleration then
different on each support).
4.1
Calculation of a second member in mono support
One considers seismic modeling beam of the building engine 1300 MW of the nuclear thermal power station
CIVAUX (case test SDLL109B) as presented in the documentation of the operator
POST_ELEM [U4.81.22].
One wishes to determine efforts at the time of the transitory dynamic response of the structure to a seism
in direction X. the transitory calculation of answer is carried out here by modal recombination by
DYNA_TRAN_MODAL [U4.53.21].
One calculates the modes of vibrations of the structure of the model resting on only one elastic support
(comes out from ground):
# --- seek clean modes of vibration -------------------------
MODES = MODE_ITER_SIMULT (MATR_A = RIGIDITY, MATR_B = MASS,
CALC_FREQ =_F (OPTION = “PLUS_PETITE”,
NMAX_FREQ = 33));
One defines the accélérogramme seism:
# --- excitation -------------------------------------------------------
LBNSNL1 = DEFI_FONCTION (NOM_RESU = “ACCE”, NOM_PARA = “INST”,
PROL_GAUCHE = “EXCLUDED”, PROL_DROIT = “EXCLUDED”,
VALE = (
0.00000E+00 9.98700E-02 1.00000E-02 6.60700E-02
2.00000E-02 - 5.65000E-03 3.00000E-02 - 9.46800E-02
-----------------------------------------------------
1.19800E+01 1.68110E-01 1.19900E+01 8.80300E-02
1.20000E+01 0.00000E+00 9.98700E-02 0.00000E+00)
);
ACCELERO = CALC_FONCTION (COMB =_F (FONCTION= LBNSNL1, COEF= 1.47));
One calculates the second member (field with the nodes of the inertias of drive) and one defines
direction of the seism
DIRSEISM = CALC_CHAR_SEISME (MATR_MASS = MASS,
MONO_APPUI = “YES”,
DIRECTION = (- 1., 0., 0. ));
Handbook of Utilization
U4.6- booklet: Stamp/Vecteurs elementary and assembly
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titrate:
Operator CALC_CHAR_SEISME
Date:
08/02/05
Author (S):
Y. PONS, Fe WAECKEL, L. VIVAN
Key: U4.63.01-G Page:
5/6
One carries out the calculation of the transitory response in modal space
# --- projection of the matrices and vector assembled on the modes
MACRO_PROJ_BASE (BASE=MODE,
NB_VECT=33, MATR_ASSE_GENE= (
_F (
STAMP = CO (“MASSGENE”),
MATR_ASSE = MASS),
_F (
STAMP = CO (“RIGIGENE”),
MATR_ASSE = RIGIDITY)),
VECT_ASSE_GENE=_F (
VECTOR = CO (“VECTGENE”),
VECT_ASSE = SEISM)
);
# --- calculation by modal combination ------------------------------------
LISTAMOR= (0.055, 0.055, 0.070, 0.070, 0.071, 0.072, 0.157, 0.085, 0.086,
0.070, 0.076, 0.074, 0.071, 0.072, 0.115, 0.073, 0.076, 0.086,
0.081, 0.070, 0.072, 0.075, 0.074, 0.070, 0.152, 0.148, 0.074,
0.297, 0.074, 0.075, 0.089, 0.138, 0.118,)
TRANGENE=DYNA_TRAN_MODAL (MASS_GENE=MASSGENE, RIGI_GENE=RIGIGENE,
METHODE=METHODE,
AMOR_REDUIT=LISTAMOR,
INCREMENT=_F (INST_INIT = 0.,
INST_FIN = 4., PAS = 0.002,
VERI_PAS = “NOT”),
EXCIT=_F (VECT_GENE = VECTGENE,
FONC_MULT = ACCELERO),
))
One used the accélérogramme here imposed movement.
One passes by again in physical space:
# --- restitution in physical base -------------------------------------
LISTINST=DEFI_LIST_REEL (DEBUT=1.334,
INTERVALLE=_F (JUSQU_A = 3.154, A NUMBER = 1))
TRANPHYS=REST_BASE_PHYS (RESU_GENE=TRANGENE,
CRITERE=' RELATIF', PRECISION=1.E-06,
LIST_INST=LISTINST, TOUT_CHAM=' OUI')
4.2
Calculation of a second member in multi supports
The example which follows watch the calculation of the second members when the seisms are different on
supports.
Handbook of Utilization
U4.6- booklet: Stamp/Vecteurs elementary and assembly
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titrate:
Operator CALC_CHAR_SEISME
Date:
08/02/05
Author (S):
Y. PONS, Fe WAECKEL, L. VIVAN
Key: U4.63.01-G Page:
6/6
One considers modeling in elements of beam of the line of following piping on which one
seism is applied in direction X:
P1
Z
X
P2
P3
One calculates the static modes to express the vector of drive like a combination
linear of those.
MODST = MODE_STATIQUE (MATR_RIGI = RG,
MATR_MASS = MS,
MODE_STAT =_F (NODE = (“P1”, “P2”, “P3”),)
(AVEC_CMP = “DX”),));
One must calculate three second members different (field with the nodes from the inertias
of drive) for each node or node subjected to the same seismic signal groups.
SMP1 = CALC_CHAR_SEISME (MATR_MASS = MS,
DIRECTION = (1., 0., 0.),
MODE_STAT = MODST,
NOEUD = “P1”);
SMP2 = CALC_CHAR_SEISME (MATR_MASS = MS,
DIRECTION = (1., 0., 0.),
MODE_STAT = MODST,
NOEUD = “P2”);
SMP3 = CALC_CHAR_SEISME (MATR_MASS = MS,
DIRECTION = (1., 0., 0.),
MODE_STAT = MODST,
NOEUD = “P3”);
In accordance with what is explained in the reference document [R4.05.01] the operator
CALC_CHAR_SEISME calculates the following expression:
- M (
+
+
P1, X
P2, X
P,
3 X) S
with S (1., 0., 0.), direction of the seism.
Calculation of the transitory response (without damping) with three different seismic signals
ACCELP1, ACCELP2 and ACCELP3:
TRANGENE = DYNA_LINE_TRAN (MATR_MASS = MS,
MATR_RIGI = RG,
NEWMARK =_F (),
LIST_INST = LI,
EXCIT = (
_F (VECT_ASSE = SMP1, FONC_MULT = ACCELP1),
_F (VECT_ASSE = SMP2, FONC_MULT = ACCELP2),
_F (VECT_ASSE = SMP3, FONC_MULT = ACCELP3),
)
);
One will be able to find other examples of multimedia structures subjected to stresses
seismic by consulting tests SDLD103 and SDND102.
Handbook of Utilization
U4.6- booklet: Stamp/Vecteurs elementary and assembly
HT-66/05/004/A
Outline document