（1）前言

 地震及其破坏性是众所周知的，在工程设计中，建筑物是否能抵御地震的破坏是不可忽视的设计技术指标，有些建筑工程在设计工作完成以后往往还需进行原型或模型的抗震试验，以检验地震对建筑物的破坏程度。

http://s15/mw690/002XKryizy7kh39spFcee&690

   
    目前我国已制造及引进了各种振动试验设备，特别是大型模拟地震振动台的建造和引进，使在实验室内模拟实际地震成为现实。为了模拟地震，一般可采用由地震现场记录到的实际地震波来做为振动试验设备的输入信号。但在实际工程设计、试验及计算中仍希望根据预定条件来生成人工模拟地震波。
 
    我们知道地震振动可描述为一个时间函数或反应谱（富氏谱或功率谱）。反应谱可以这样计算：几个具有相同阻尼比和不同自振频率的单自由度振动系统（弹簧--阻尼系统），在相同地震激励信号的作用下，记录各振动系统的最大反应（例如：加速度），并将其绘制在自振频率图上，这就是关于特定阻尼下振动系统的反应谱。
 
    本程序是用期望地震反应谱做为生成地震波的基础，期望地震反应谱可以根据建筑工程结构、地质条件、设计规范及综合各种地震反应谱等设计。
 
    生成人工模拟地震波程序是在我们编制的Vib'SYS数字信号分析、处理程序支持下的应用程序，程序由两部分组成：

    输入地震波期望反应谱程序
    生成人工模拟地震波程序

 http://www.earthquake.com.cn/software/sqk/resp.gif
（2） 生成人工模拟地震波原理

    地震反应谱是用实际地震波时域信号计算得到的，而该运算是不可逆的。在生成地震波前，先用给定的期望地震反应谱来计算功率谱，再根据功率谱生成人工地震波。

    对于给定自振频率和阻尼的单自由度振动系统，建立其标准差 与反应谱 关系式
 

http://www.earthquake.com.cn/software/sqk/imageUFK.JPG

(1)

 
    对于系数项到目前为止还没有准确的方法来计算，所以我们用近似方发来逼近 。
    首先解释标准差，对于窄带白噪声（均值为零）的方差表示为
 

http://www.earthquake.com.cn/software/sqk/imageJE9.JPG

(2)

    对在平稳随机信号激励下，线性系统的频响函数 为

http://www.earthquake.com.cn/software/sqk/imageTAF.JPG

(3)

    其中是反应信号的功率谱， 是激励信号的功率谱, 这样反应信号的方差表示为
 

http://www.earthquake.com.cn/software/sqk/image6LF.JPG

(4)

    由标准差，反应信号的加速度伪标准差 可用下式计算（ 为自振角频率）
 

http://www.earthquake.com.cn/software/sqk/imageNHF.JPG                                        (5)

    当输入为加速度信号，输出为相对位移时，单自由度振动系统的频响函数（增益因子）按下式计算
 

http://www.earthquake.com.cn/software/sqk/imageMBP.JPG

(6)

    其中：http://www.earthquake.com.cn/software/sqk/imageQF2.JPG为自振角频率
          ω为强迫振动角频率
          ξ为阻尼比

    对于小阻尼振动系统，在伯德（Bode）图上，频响函数 在处显示出很尖的峰值。以此为先决条件，反应加速度伪标准差由下式计算（把(6)和(3)式代入(5)整理得到
 

http://www.earthquake.com.cn/software/sqk/imageN6L.JPG
http://www.earthquake.com.cn/software/sqk/imageTAV.JPG	(7)

   上式是以平稳随机过程条件的。对比较短的地震过程就不能满足此条件。
    Corotic和Vanmarcke用与时间相关的谱密度函数Ｇ(ω,t)来近似计算标准方差
 

http://www.earthquake.com.cn/software/sqk/imageP7J.JPG

(8)

     函数 依据于输入信号的谱密度函数 和系统的特性。对线性系统，单位脉冲的频响函数定义为
 

http://www.earthquake.com.cn/software/sqk/image9S6.JPG

(9)

    当 时，函数收敛于频响函数 。则与时间相关的谱密度函数表示为
 

http://www.earthquake.com.cn/software/sqk/image8OR.JPG

(10)

    对有如下的近似计算公式
 

http://www.earthquake.com.cn/software/sqk/image6DJ.JPG

(11)

    当ｔ远大于http://www.earthquake.com.cn/software/sqk/imageMGH.JPG 。把这个值与(11)式比较，假设与时间相关的阻尼定义为

http://www.earthquake.com.cn/software/sqk/imageIKG.JPG

(12)

    将http://www.earthquake.com.cn/software/sqk/imageH8D.JPG 代入(6)式，即得出与时间相关的频响函数的近似计算公式
 

http://www.earthquake.com.cn/software/sqk/imageG0K.JPG

(13)

    对于标准差的系数项 也要进行近似计算。Vanmarcke用下式来计算
 

http://www.earthquake.com.cn/software/sqk/image9G7.JPG
http://www.earthquake.com.cn/software/sqk/imageMBJ.JPG	(14)

   

式中：	http://www.earthquake.com.cn/software/sqk/imageMV3.JPG 为在单自度周期内的振动量
	http://www.earthquake.com.cn/software/sqk/imageO1E.JPG 为反应谱谱线 i的角频率  http://www.earthquake.com.cn/software/sqk/imageN3P.JPG 为地震振动周期，取决于包络线的形式和周期  http://www.earthquake.com.cn/software/sqk/imageE9M.JPG  概率系数(0->1.0)

 
    用概率系数http://www.earthquake.com.cn/software/sqk/image86N.JPG 。
  
    根据以上公式的推导过程，把(1)和(7)式整理得出地震波功率谱的计算公式
 

http://www.earthquake.com.cn/software/sqk/imageOSC.JPG

(15)

    由期望反应谱计算出地震波功率谱以后，用随机信号发生器来生成随机相位，再经富氏逆变换来生成地震波。

 

 

（3） 包络

          为了使生成的人工地震波的形式更接近实际地震波，在程序中我们设计了四种包络线形式。根据需要可以用任何一种包络线来乘生成的地震波。

	无包络线(No Envelope)
	这种情况生成的地震波信号的包络线取省略梯型包络线形式，其包络线参数AT=Tmax/6、UT=4×Tmax/6、DT=Tmax/6，其中Tmax为地震波的最大时程。
	梯型包络线(Trapezoid)〔梯型包络线见附图一所示〕 参数上升时间AT，持续时间UT和下降时间DT是可以根据需要自由选择的。AT+UT+DT等于地震波最大持续时间。
	指数函数包络线(Exponential)〔指数函数包络线见附图二所示〕     http://www.earthquake.com.cn/software/sqk/imageQDG.JPG      参数AM,F1,F2和信号总时间长度Ｔ是可以根据需要自由选择的。
	复和包络线(Combined Envelope)〔复和包络线见附图三所示〕

 

（4）地震波基线和幅值校正

      为了保证所生成的地震波加速度信号在其结束处的位移和速度为零，我们用三次多项式对已生成的地震波进行校正。 

http://www.earthquake.com.cn/software/sqk/imageT8M.JPG

    一般来说，生成的地震波加速度信号的幅值应达到期望加速度幅值，但实际是做不到的。当已生成的地震波的幅值超过期望幅值时，我们对整个地震波用下面系数来乘
 

http://www.earthquake.com.cn/software/sqk/image5AE.JPG

   

其中：	http://www.earthquake.com.cn/software/sqk/imageELP.JPG 是期望地震波加速度幅值
	http://www.earthquake.com.cn/software/sqk/imageSN1.JPG  是实际地震波加速度幅值

    如果已生成的地震波的幅值低于期望幅值时，我们把http://www.earthquake.com.cn/software/sqk/imageCNU.JPG，这似乎是不可理解的，但在以后生成地震波的迭代过程中会逐步修正的。

 

（5）  迭代─地震波功率谱校正

 

      根据期望地震反应谱，我们近似计算出地震波的功率谱，经氏富逆变换生成出地震波，再用包络线乘地震波、基线校正和幅值校正以后，我们就来计算已生成的地震波的绝对加速度反应谱，然后用期望反应谱Ｓdes(ω)和计算反应谱Ｓach(ω)按下式来修正地震波的功率谱
 

http://www.earthquake.com.cn/software/sqk/image2M9.JPG

i=1,2,...,N

    这样的一次修正过程称为迭代一次，整个迭代过程根据需要进行Ｎ次，当然计算得到的地震反应谱是永远达不到期望地震反应谱的，但误差随迭代的次数的增加会逐渐减小的。

 

（六）计算地震波绝对加速度反应谱

       地震反应谱的计算是采用了Duhamel积分法， 对于单自由度振动体系我们有下面振动方程：
 

http://www.earthquake.com.cn/software/sqk/imageVRO.JPG

(1)

    式中，x 为相对位移；   ω   为自振角频率；
          ξ 为临界阻尼比；a(t) 为地震波加速度记录。

    a(t)可以看成由一系列等步长(t=1/Fsample)的彼此相连的直线段构成。 在ti≤t≤ti+1，有
 

http://www.earthquake.com.cn/software/sqk/image39F.JPG
(i=1,2,...,N)	(2)

于是，式(1)的Duhamel积分可表示为：
 

http://www.earthquake.com.cn/software/sqk/image64O.JPG

(3)

式中，Ｍ和Ｌ分别为二阶矩阵，每个矩阵元素都是ξ，ω和ｔ的函数。对于给定的ξ和ω，且ｔ为常数时，Ｍ和Ｌ的每一个元素都是常数。

可以证明，这时矩阵Ｍ的诸元素为：
 

http://www.earthquake.com.cn/software/sqk/imageSHR.JPG   http://www.earthquake.com.cn/software/sqk/image0V3.JPG   http://www.earthquake.com.cn/software/sqk/image613.JPG

Ｌ矩阵的诸元素分别为
 

http://www.earthquake.com.cn/software/sqk/image5VH.JPG   http://www.earthquake.com.cn/software/sqk/image2KR.JPG   http://www.earthquake.com.cn/software/sqk/imageG7G.JPG    http://www.earthquake.com.cn/software/sqk/imageOLG.JPG   http://www.earthquake.com.cn/software/sqk/imageA8I.JPG    http://www.earthquake.com.cn/software/sqk/imageGKM.JPG

式中, http://www.earthquake.com.cn/software/sqk/imageMU3.JPG

    给定频率 f，阻尼比ξ和步长t，按照以上Ｍ和Ｌ矩阵中诸元素的计算公式则可求出常数矩阵Ｍ和Ｌ。然后，按照式(2)，并假定初始速度和初
始位移为0，就可逐步计算出相对位移X(ti)和相对速度X(ti)(i=1,2,...,N)

    由式(1) 既可得到绝对加速度反应：
 

http://www.earthquake.com.cn/software/sqk/imageMGT.JPG

(4)

从Ｎ个反应中，挑出绝对值最大者，就得到了反应谱上的一个点。

    生成地震波的期望反应谱是绝对加速度反应谱，绝对加速度反应谱与工程单位动力放大倍数β之间的关系：
 

http://www.earthquake.com.cn/software/sqk/imageFTN.JPG

其中：	http://www.earthquake.com.cn/software/sqk/imageQOO.JPG  绝对加速度反应谱 [g]
	http://www.earthquake.com.cn/software/sqk/image4E4.JPG    为地震波最大加速度幅值 [g]

 

 

本文来源：http://www.earthquake.com.cn/software/sqk/sqk.html