ACCELERATION GENERATING
CORRESPONDING TO THE PROPOSED RESPONSE SPECTRA
Jan Dupal
1a1University of West Bohemia in Pilsen, Faculty of Applied Sciences, Dept. of Mechanics e-mail: adupal@kme.zcu.cz
Abstract
Nuclear devices are checked also from the seismic safety point of view. The response spectra method is one of the most important approaches to obtain an upper estimation of some output variables like displacements, velocities, accelerations, stresses etc. The proposed response spectra of the base especially accelerations serves as an input for extreme estimations of mentioned output quantities depending on frequency [1]. The alternative approach for functionality assessment can be used the base acceleration as an input for time integration of the mathematical model and the output quantities can be obtained in the direct integration way of the model. The main require- ment on the input acceleration is posed on the desired response spectrum satisfaction. The paper deals with an approach to generation of the base acceleration in such a way to be satisfied the properties of the proposed pre- scribed response spectrum of the acceleration.
GENEROVÁNÍ AKCELEROGRAMU NA ZÁKLADĚ NÁVRHOVÉHO SPEKTRA ODEZVY
Summary
Bezpečnost jaderných zařízení je kontrolována mimo jiné také ze seizmického hlediska. Metoda spekter odezvy patří mezi nejdůležitější metody pro určení horního odhadu výstupních veličin, za které považujeme posuvy, rychlosti, zrychlení, napětí atd. Presentovaný přístup pro posouzení funkčnosti využívá zrychlení základu odpovídající návrhovému spektru odezvy jako kinematické buzení matematického modelu a výstupní veličiny lze získat přímou integrací.
1. INTRODUCTION
What is it the response spectrum? It can be defined as a maximum of absolute value of one DOF vibrating system response to a prescribed excitation. The corresponding mathematical model has form
( ) t b q ( ) t kq ( ) t f ( ) t , q
m & & + & + =
(1)where
f ( ) t
is the mentioned excitation. The response spectrum of such system is expressed as( ) max 1 ( )
( )sin ( ) ,
max
max 0
,
0
Ω ∫ Ω −
=
Ω
− Ω −∈
T
D t
D
T D
q t
f e t d
S m τ
ττ τ
(2)where
Ω = k / m
is natural frequency of the un-damped system,Ω
D= Ω 1 − D
2 andkm D b
2
=
is damping ratio which is supposed to be input constant parameter. When the system is excited by the base movement, the excitation can be expressed asf ( ) t = − m u & & ( ) t
and the response spectrum of displacement has form( ) max 1 ( )
( )sin ( ) .
max
max 0
,
0
− Ω ∫ Ω −
=
Ω
− Ω −∈
T
D t
D
T D
q t
u e t d
S & & τ
ττ τ
(3)Between response spectrum of displacement and response spectrum of acceleration holds simple relation
( ) Ω = Ω
2 q( ) Ω .
q
S
S
&& (4)2. SOLUTION OF THE RESPONSE
The relation (3) expresses a maximum of the absolute value of the system described by eq. of motion
( ) t 2 D q ( ) t
2q ( ) t u ( ) t .
q & & + Ω & + Ω = − & &
(5)As a start point of the presented method is replacement of the kinematic excitation by periodic excitation with fre- quency step
∆ ω ,
it means by Fourier series with the basic frequency∆ ω
. In case the frequency step is variable the substitutive excitation will be poly-harmonic. The example of the proposal response spectrum is depicted in Fig. 1Fig. 1 Example of the proposal response spectrum
Having constant frequency step
∆ ω
we can the frequency values chose in form, ,
1 ..., , 1 , 0
,
minmin
ω ω ω
ω
ω
k= + k ∆ k = N − = ∆
(6)it means N values. It means that
ω
max= ω
min+ ( N − 1 ) ∆ ω .
In case the step∆ ω
is variable the frequency values of poly-harmonic excitationω
k are chosen from the set{
1,
2, ... ,
N} ,
k
ω ω ω
ω ∈
(7)
where
ω
1= ω
min, ω
max= ω
N.
In both cases the eq. (5) takes a form( ) 2 ( ) ( ) ( cos sin ) .
1
2
∑
=
+
−
= Ω + Ω +
N
k
k k k
k
b
a t
q t
q D t
q & & ω τ ω τ
&
(8)
The goal is to find such values
a
k, b
k, k = 1 , 2 , ..., N ,
to achieve the result acceleration response spectrum of the substitutive excitation be identical with the proposal response spectrum in the defined frequency pointsω
k. As theexperience shown the calculation of convolution integrals in (3) and (4) is very time consuming and for this reason the analytical solution by means of Laplace transform of the eq. (8) was taken which respects, despite of Fourier trans- form, the zeros initial conditions. The Laplace transform of (8) can be expressed as
( ) ( ) ( ) ∑
=
+ +
− +
= Ω
+ Ω
+
N
k k
k k k
k
b p
p a p p
Q p
pQ D p Q p
1
2 2 2
2 2
2
2
ω ω
ω
, (9)with solution
( )
1( )( 2 ) ( )( 2 ) ,
2 2
2 2 2
2 2
∑
2=
Ω + Ω + + +
Ω + Ω +
− +
=
N
k k
k k
k
k
b p p D p
p D p p
a p p
Q ω
ω
ω
(10)where
L { q ( ) t } = Q ( ) p .
To obtain an inverse Laplace transform of (10) it is necessary to decompose the relation (10) into so called partial fractions according to( )
( ) ∑
=
−
=
n
j j
j
p p
p c Z
p P
1
1
, (n is number of roots of the Eq.Z ( ) p = 0
) (11)where
( ) ( )
jj j
p dp Z
d p c = P
. (12)
It can be written for both fractions in (10)
( ) p ( p
k)( p D p ) p i k p D i
D
Z =
2+ ω
2 2+ 2 Ω + Ω
2= 0 ⇒
1,2= ± ω ,
3,4= − Ω ± Ω
(13)
( ) p = p ( p + D Ω p + Ω ) ( + p + ) ( p + D Ω )
dp Z d
k
2 2
2
2
2 2 2ω
2. (14)
Now we can come to the coefficients of the first fraction in (10)
( 2 ) ,
2
1
2 , 2
1
= − + Ω + Ω
k k
k
i D
c ω ω 2 ( 2 ) ,
1
, 2 1 , 2
2 k
k k
k
c
D
c i =
Ω + Ω
−
= −
ω
ω
(15)( 2 ) ,
2
2 2 2 2, 3
k D D
D
D
k
i D i D
i c D
ω + Ω
− Ω Ω
− Ω Ω
Ω + Ω
= −
(
D D k)
kD
D
k
c
D i D i
i
c
4, 2 2D
2 2 3,2
2 =
+ Ω
− Ω Ω + Ω Ω
−
Ω
− Ω
= −
ω
(16) and for the second fraction in (10)
( 2 ) ,
2
1
, 2 1
2
* ,
1 k
k k
k
ic
D i
c i = −
Ω + Ω +
= −
ω
ω 2 ( 2 ) ,
1
*, 2 1 2
* ,
2 k
k k
k
c
D i
c i =
Ω + Ω
−
−
= −
ω
ω
(17)( 2 ) ,
2
2 2 2 2* , 3
k D D
D
k
k
i D i D
c ω
ω
+ Ω
− Ω Ω
− Ω
= Ω
( 2 ) .
2
* , 2 3 2 2
2
* ,
4 k
k D D
D
k
k
c
D i D
c i =
+ Ω
− Ω Ω + Ω Ω
= −
ω ω
(18)
t p
j
e
jp
p =
−
−1
1
L
(19)we can obtain the total solution of Eq. (8) (inverse Laplace transform of (10)) in form
( ) { [ ( ) ( ) ]
( ) ( )
[ ]} .
,
* , 3
* , 3
* , 1
* , 1 1
, 3 ,
3 ,
1 ,
1
t i D k t i D k t i k t i k k N
k
t i D k t i D k t i k t i k k
D D
k k
D D
k k
e c e
c e c e c b
e c e
c e c e c a t
q
Ω
− Ω
− Ω
+ Ω
−
−
=
Ω
− Ω
− Ω
+ Ω
−
−
+ +
+ +
+ +
+ +
−
=
Ω ∑
ω ω
ω ω
(20)
From the analytical solution (20) it can be simply determined response spectrum of displacement and acceleration
( ) ( )
( ) ( )
[ ]
{
( ) ( )
[ ]}
, max
, max ,
* , 3
* , 3
* , 1
* , 1 1
, 3 ,
3 ,
1 ,
1 ,
0 ,
0 max max
t i D k t i D k t i k t i k k N
k
t i D k t i D k t i k t i k k T
t T
q t
D D
k k
D D
k k
e c e
c e
c e c b
e c e
c e
c e c a t
q S
Ω
− Ω
− Ω
+ Ω
− −
=
Ω
− Ω
− Ω
+ Ω
−
−
∈
∈
+ +
+ +
+ +
+ +
=
=
Ω ∑
ω ω
ω ω
p p
(21)
( ) ( )
( ) ( )
[ ]
{
( ) ( )
[ ] }
. max
, ,
* , 3
* , 3
* , 1
* , 1 1
, 3 ,
3 ,
1 ,
1 ,
0 2 2
max D i t
k t i D k t i k t i k k N
k
t i D k t i D k t i k t i k k T
q t q
D D
k k
D D
k k
e c e
c e
c e c b
e c e
c e
c e c a S
S
Ω
− Ω
− Ω
+ Ω
−
−
=
Ω
− Ω
− Ω
+ Ω
−
−
∈
+ +
+ +
+ +
+ +
Ω
= Ω Ω
=
Ω ∑
ω ω
ω ω
p
&
p
&
(22) The terms in (20) are complex conjugated pairs and for this reason the response (20) is of course real. In relations (21) and (22) was used a vector symbol p which means
, ,
,
21
2 1
2 N
N N
N N
b b b
a a a
R b
R a
b R
p a ∈
=
∈
=
∈
=
M
M
(23)finding of which will be goal of optimization.
3. FORMULATION OF THE OPTIMIZATION PROBLEM [2]
The goal of optimization process will be the minimization of square difference sum proposal response spectrum
( ) Ω
*
S
q&& and current spectrumS
&q&( Ω , p )
. Both spectra will be expressed in pointsω
k∈ { ω
1, ω
2, ... , ω
N} .
Firstly we assemble the discrete values of both spectra into the vectors
( )
( )
( )
( )
( )
( ) ( ) ( )
. ,
, , ,
,
2 ,11
1 ,
* 2
* 1
*
* 1 2 ,
1
N
N N
N q q q
q N
N q
q q
q
S S S
S S S
Ω R Ω R
S R p p p Ω p
S ∈
=
∈
=
∈
=
ω ω ω
ω ω ω
ω ω ω
M M M
&
&
&
&
&
&
&
&
&
&
&
&
&
&
&
&
(24) The cost function can be written down in form
( ) p = [ S
&q&( Ω , p ) − S
*&q&( ) Ω ]TW [ S
&q&( Ω , p ) − S
*q&&( ) Ω ] / ( S
*q&&T( ) Ω WS
q*&&( ) Ω ) ,
ψ
(25)where
W ∈ R
N ,N is diagonal matrix of weighting parameters which will be chosen as identity matrix. In the opposite case this matrix enables the preference of representative frequencies for better coincidence of both spectra in these frequencies. Let us introduce the functional vectors( )
( ) ( )
( ) ( )
( ) ( )
, ,
3 3
1 1
32 32
12 12
31 31
11 11
+ +
+
+ +
+
+ +
+
= Ω
Ω
− Ω
− Ω
+ Ω
−
−
Ω
− Ω
− Ω
+ Ω
−
−
Ω
− Ω
− Ω
+ Ω
−
−
t i D N t i D N t i N t i N
t i D t
i D t
i t
i
t i D t
i D t
i t
i
D D
k k
D D
k k
D D
k k
e c e
c e
c e c
e c e
c e
c e c
e c e
c e
c e c t
ω ω
ω ω
ω ω
a
M f
(26)
( )
( ) ( )
( ) ( )
( ) ( )
+ +
+
+ +
+
+ +
+
= Ω
Ω
− Ω
− Ω
+ Ω
− −
Ω
− Ω
− Ω
+ Ω
−
−
Ω
− Ω
− Ω
+ Ω
−
−
t i D N t i D N t i N t i N
t i D t
i D t
i t
i
t i D t
i D t
i t
i
D D
k k
D D
k k
D D
k k
e c e
c e
c e c
e c e
c e
c e c
e c e
c e
c e c t
* , 3
* , 3
* , 1
* , 1
* 2 , 3
* 2 , 3
* 2 , 1
* 2 , 1
* 1 , 3
* 1 , 3
* 1 , 1
* 1 , 1
,
ω ω
ω ω
ω ω
b
M f
(27)
and the total functional vector
( ) ( )
( , ) .
, ,
Ω
= Ω
Ω t
t t
b a
f f f
(28)
Then we can rewrite the relation (22) into form
( , ) max ( , ) ( , ) max ( , ) .
max
max 0,
2 ,
0
2
− Ω − Ω = Ω − Ω
Ω
= Ω
∈
∈
t t t
S
TT t T
T T
q&&
p
ta f
ab f
bp f
(29)
The proposal response spectrum of acceleration of the nuclear base for
D = 0 . 02
is depicted in Fig. 2 (blue color).Fig. 2 Response spectra of acceleration
The generated acceleration corresponding to the achieved response spectrum (Fig. 2-green line) is depicted in Fig. 3
Fig. 3 Generated acceleration
4. CONCLUSION
The presented approach will be used for safety assess- ment of the nuclear power station Temelin especially its reactor and pipeline system. The requirement of three acceleration functions for kinematic excitation of base in three directions (horizontal x, z and vertical y) was posed by Nuclear Research Institute in Rez. Only one of these three accelerations was presented in this paper for the reason of space lack.
This contribution was supported by the Europe- an Regional Development Fund (ERDF), project
“NTIS - New Technologies for Information Society”, European Centre of Excellence, CZ.1.05/1.1.00/02.0090.
REFERENCES
1. Dimarogonas A.: Vibration for engineers. Prentice Hall inc., 1996.
2. Dupal J.: Generating of the synthetical accelerogram from the proposal acceleration response spectrum, report 52 220-01-12, Univerzity of West Bohemia in Pilsen, 2012. (In Czech)