Version

Code_Aster default
Titre : SDLD33 - Réponse sismique spectrale d'un système 2[...] Date : 09/09/2020 Page : 1/11
Responsable : Clé : V2.01.033 Révision :
4ea741f2e6e9

SDLD33 - Spectral seismic answer of a system 2

masses - 4 springs multimedia by 2 décorrélés
groups of supports

Summary:

The problem consists in calculating the spectral response of a system 2 masses – 4 springs multimedia,
subjected to a multiple seismic excitation, by regarding the 3 supports as 2 groups of décorrélés supports.

One tests the discrete element in traction, the calculation of the clean modes, the static modes and the spectral
response by modal superposition via the operator COMB_SISM_MODAL. This CAS-test makes it possible to
validate the order of the combinations to be considered in this case, namely the office plurality intra-group then
the office plurality joint committee.

The got results are in very good agreement with the analytical results of reference.

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 2020 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : SDLD33 - Réponse sismique spectrale d'un système 2[...] Date : 09/09/2020 Page : 2/11
Responsable : Clé : V2.01.033 Révision :
4ea741f2e6e9

1 Problem of reference
1.1 Geometry
The system is composed of a set of 4 springs, 2 specific masses, is supported by 3 supports.

  m 1 m 2
k1 k2 k3 k4

x 
N O 1  N O 2  N O 3  N O 4  N O 5 

1.2 Material properties

k =k 1 =k 2  =103  N / m
Stiffness of connection: k 3  =k 4 =2 k=2.10 3  N / m ;
specific mass: m=m 1 =m 2  =10  kg .

1.3 Boundary conditions and loadings

Boundary conditions :
Only authorized displacements are the translations according to the axis x .
Points NO1 , NO3 and NO5 are embedded: dx=dy=dz=drx=dry=drz=0 .
The other points are free in translation according to the direction x : dy =dz =drx=dry=drz=0 .

Loading :
The structure is subjected to a multiple spectral seismic excitation and differential displacements.
The spectra of answers of oscillator in pseudo-acceleration are simplified. Only the values
corresponding to the 2 Eigen frequencies of the system are mentioned. They do not depend on
damping:

with the node NO1 :

SRO NO1  f 1 =A11=7 m/ s2
SRO NO1  f 2 = A21 =5 m/ s 2
DDS NO1 =D 1=−0.04 m

with the node NO3 :

SRO NO3  f 1 =A11=7.7 m/ s2
SRO NO3  f 2 = A12 =5.5 m/ s 2
DDS NO3 =D 2=−0.044 m

with the node NO5 :

SRO NO5  f 1 = A21=12 m/ s 2
SRO NO5  f 2= A22 =6 m/s 2
DDS NO5 = D 3=0.06 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 2020 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : SDLD33 - Réponse sismique spectrale d'un système 2[...] Date : 09/09/2020 Page : 3/11
Responsable : Clé : V2.01.033 Révision :
4ea741f2e6e9

Excitations with the nodes NO1 and NO3 are correlated. One sets up 2 groups of décorrélés
supports: group 1 is composed of the nodes NO1 and NO3 ; group 2 is made up by the only node
NO5 .

1.4 Initial conditions

The system is at rest.

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 2020 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : SDLD33 - Réponse sismique spectrale d'un système 2[...] Date : 09/09/2020 Page : 4/11
Responsable : Clé : V2.01.033 Révision :
4ea741f2e6e9

2 Reference solution
2.1 Method of calculating used for the reference solution
One calculates the spectral response by modal superposition of a system 2 masses – 4 springs
subjected to three distinct excitations. One determines the displacement of the masses to the nodes
NO2 and NO4 along the axis x .
One calculates analytically:
• Eigen frequencies fi,
• associated clean vectors  N i standardized compared to the modal mass,
• static modes of supports  i system,
• factors of modal participation P ij relating to the supports,
• Rm ij the maximum of answer of each mode starting from the spectra of excitation,
• Re ij the contribution of the movement of training of each support starting from
differential displacements,

2.2 Results of reference

2.2.1 characteristic Matrices and vectors

• matrix of rigidity K

[ ]
k1 −k 1 0 0 0
−k 1 k 1k 2 −k 2 0 0
K= 0 −k 2 k 2k 3 −k 3 0 matrix relating to the degrees of freedom 1,2,3,4,5
0 0 −k 3 k 3k 4 −k 4
0 0 0 −k 4 k4

[ ]
k 1k 2 0 −k 1 −k 2 0
0 k 3k 4 0 −k 3 −k 4
p
K = −k 1 0 k1 0 0
−k 2 −k 3 0 k 2k 3 0
0 −k 4 0 0 k4

partitionnée matrix degrees of freedom of structure 2.4,

degrees of freedom of support 1,3,5

K p=
[ k k xs
k sx k ss ] [
k=
k 1k 2
0
0
k 3k 4 ] [ k xs =
−k 1 0 −k 2 0
0 −k 3 0 −k 4 ]
• matrix of mass 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 2020 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : SDLD33 - Réponse sismique spectrale d'un système 2[...] Date : 09/09/2020 Page : 5/11
Responsable : Clé : V2.01.033 Révision :
4ea741f2e6e9

[ ]
0 0 0 0 0
0 m1 0 0 0
M= 0 0 0 0 0
0 0 0 m2 0
0 0 0 0 0
matrix relating to the degrees of freedom 1,2,3,4,5

[ ]
m1 0 0 0 0
0 m2 0 0 0
p
M = 0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
partitionnée matrix degrees of freedom of structure 2.4, degrees of freedom of support 1,3,5

• modal calculation in embedded base

 K −i M  i=0 with i= 2i

•
dét  K − i M =0⇔ 2i −
 k 1k 2 k 3k 4
m1

m2   k k  k k 
i 1 2 3 4 =0
m1 m 2

2k 4k
 1=  2=
m m

• Eigen frequencies:

1 2
⇒ f 1= f 2=
2 2

• not normalized clean modes:

[] []
0 0
1 0
1= 0  2= 0
0 1
0 0

• generalized modal masses  i=Ti M i :

•  1= 2=m

• own standards modes with the unit generalized modal mass  Nk :

1 2
• ⇒  N1 =  N2 =
 1  2
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 2020 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : SDLD33 - Réponse sismique spectrale d'un système 2[...] Date : 09/09/2020 Page : 6/11
Responsable : Clé : V2.01.033 Révision :
4ea741f2e6e9

• static modes of supports  Sj

Matrix of the static modes reduced to the ddls of structure  S =−k

−1
k xs

[ ]
1 1
0
S= [ ][
−1 2 0 −k −k 0 =
4 k 0 1 0 −2 k −2 k ] 0
2
1
2
1
2 2

• static solution with a unit displacement of the node NO1 :

[]
1
1
displacements:  S1 = 2
0
0
0

• static solution with a unit displacement of the node NO3 :

[]
0
1
2
displacements:  S2 = 1
1
2
0
• static solution with a unit displacement of the node NO5 :

[]
0
0
0
displacements:  S3 =
1
2
1

• rigid mode of body  R1

Matrix of the rigid modes reduced to the ddls of structure:  R= S S R

[]
1
1
Rigid mode of body  R1= 1
1
1
It is checked well that:  R1= S1 S2 S3

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 2020 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : SDLD33 - Réponse sismique spectrale d'un système 2[...] Date : 09/09/2020 Page : 7/11
Responsable : Clé : V2.01.033 Révision :
4ea741f2e6e9

2.2.2 Loading 1 multi-supports

• factors of modal participation P kj = TNk M  j :

• contribution of the dynamic mode 1 to the movement imposed on the node NO1 :
⇒ P 11= TN1 M  1=
m
2
• contribution of the dynamic mode 1 to the movement imposed on the node NO3 :
⇒ P 12= TN1 M  2 =
m
2
• contribution of the dynamic mode 1 to the movement imposed on the node NO5 :
T
⇒ P 13= M  3=0N1

• contribution of the dynamic mode 2 to the movement imposed on the node NO1 :
T
⇒ P 21= M  1 =0
N2
• contribution of the dynamic mode 2 to the movement imposed on the node NO3 :
⇒ P 22=TN2 M  2=
m
2
• contribution of the dynamic mode 2 to the movement imposed on the node NO5 :
⇒ P 23=TN2 M  3 =
m
2

• answer of the mode i with the movement of the support j

Akj
Rm kj = Nk P kj
i2

Combined answers of the modal oscillators

A11 1 A11

Response of mode 1 to the movement of support 1:
Rm11= N1 P11 =
1 2 12 0 2

A A

Response of mode 1 to the movement of support 2:
Rm12= N1 P12 122 = 122 1
1 2 1 0
A

Response of mode 1 to the movement of support 3:
Rm13= N1 P13 132 = 0
1 0

A21

Response of mode 2 to the movement of support 1:
= 0
Rm 21 = N2 P 21 2
 0 2
A A

Response of mode 2 to the movement of support 2:
Rm 22 =N2 P 22 222 = 222 0
2 2  2 1

A23 A23 0

Response of mode 2 to the movement of support 3:
Rm 23 = N2 P 23 2
=
 2 2 22 1

• office plurality intra-group (algebraic sum)

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 2020 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : SDLD33 - Réponse sismique spectrale d'un système 2[...] Date : 09/09/2020 Page : 8/11
Responsable : Clé : V2.01.033 Révision :
4ea741f2e6e9

A11A12 1

mode 1:
Rm1groupe1 =Rm 11 Rm12=
2 21 0
A

mode 2:
Rm 2groupe1=Rm 21Rm 22= 222 0
2 2 1

• contribution of the support j with the movement of training

Re j = j D j
contributions of groups 1 and 2 to the movement of training

Re groupe1= Re 1 Re 2=
2 
1 D1  D2
D2  Re groupe2 =Re 3=
 
1 0
2 D3

• office plurality on the modes (quadratic)

 
A11 A12
2
2 1
Rm groupe1 =  Rm11 Rm12   Rm 21Rm 22 =
2 2

A22
2  22

A23 0
Rm groupe2 = Rm213 Rm 223 = 
2 22 1

• answer of the groups of supports 1 and 2

[ ]
R 1=  Rm2groupe1 Re 2groupe1  A11 A12 
2
2
2
 D1 D2 
1 
R21 = 1
4 A222
2
D 22
 2

[ 
R 2 = Rm2groupe2  Re 2groupe2 0
21 2
R=2 A
4 232 D32
2

• office plurality joint committee (quadratic)

[ ]
R=  R12R22  A11A122 2
2
 D1D 2
2 1  1
R=
4 A222 A223
2
D 22D23
 2

2.3 Uncertainty on the solution

No (analytical solution)

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 2020 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : SDLD33 - Réponse sismique spectrale d'un système 2[...] Date : 09/09/2020 Page : 9/11
Responsable : Clé : V2.01.033 Révision :
4ea741f2e6e9

3 Modeling A
3.1 Characteristics of the grid
The grid consists of 5 nodes, 4 meshs SEG2 and 2 meshs POI1.

3.2 Characteristics of modeling

The system is modelled by:

• 4 discrete elements DIS_T characteristics K_T_D_L,

• 2 discrete elements DIS_T characteristics M_T_D_N.

3.3 Parameters of modeling

Calculation of the static modes for the movements of the points NO1 , NO3 and NO5 .
Answer modelled spectral multi-supports décorrélés on the first 2 modes without static correction
(combination of modal answers SRSS), without damping.

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 2020 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : SDLD33 - Réponse sismique spectrale d'un système 2[...] Date : 09/09/2020 Page : 10/11
Responsable : Clé : V2.01.033 Révision :
4ea741f2e6e9

4 Results of modeling A
4.1 Eigen frequencies
MODE Reference Code_Aster Relative error (%)
1 2.2507900E+00 2.25079079E+00 3.5E-05
2 3.1830000E+00 3.18309886E+00 0,003

4.2 Total answer on complete modal basis

Modes 1 and 2 are taken into account. The components inertial (primary education) and statics (secondary) of
the answer are directly cumulated on the level of the supports.

COMB_MODE=' SRSS'

• answer of the support j=1 (node NO1 ): R 1=  Rm21  Re12 with

Rm1 = Rm211Rm 221

• answer of the support j=2 (node NO4 ): R 2 = Rm22 Re 22 with

Rm 2= Rm12
2
Rm 222
R=  R 1 R 2 .
2 2
• total answer:

• absolute displacements: DEPL

NODE Reference Code_Aster Relative error (%)

NO1 4.00000E-02 4.00000E-02 0,000
NO2 5.58000E-02 5.58083E-02 0,015
NO3 4.40000E-02 4.40000E-02 0,000
NO4 3.85600E-02 3.85683E-02 0,022
NO5 6.00000E-02 6.00000E-02 0,000

• nodal reactions: REAC_NODA

NODE Reference Code_Aster Relative error (%)

NO1 3.68000E+01 3.68044E+01 0,012
NO2 0.00000E+00 2.23756E-14 2.2E-14
NO3 8.64900E+01 8.64906E+01 7.0E-04
NO4 0.00000E+00 3.071843E-14 3.1E-14
NO5 7.71400E+01 7.71367E+01 -0,004

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 2020 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : SDLD33 - Réponse sismique spectrale d'un système 2[...] Date : 09/09/2020 Page : 11/11
Responsable : Clé : V2.01.033 Révision :
4ea741f2e6e9

5 Summary of the results

Perfect agreement of the results Code_Aster with the analytical values of reference.

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 2020 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)