生成地震波的过程和结果.doc

地震反应谱是用实际地震波时域信号计算得到的，而该运算是不可逆的。在生成地震波前，先用给定的期望地震反应谱来计算功率谱，再根据功率谱生成人工地震波。 ??? 对于给定自振频率和阻尼的单自由度振动系统，建立其标准差 与反应谱 关系式 ? (1) ? ??? 对于系数项 到目前为止还没有准确的方法来计算，所以我们用近似方发来逼近 。 ??? 首先解释标准差 ，对于窄带白噪声（均值为零）的方差表示为 ? (2) ??? 对在平稳随机信号激励下，线性系统的频响函数 为 (3) ??? 其中 是反应信号的功率谱， 是激励信号的功率谱, 这样反应信号的方差表示为 ? (4) ??? 由标准差 ，反应信号的加速度伪标准差 可用下式计算（ 为自振角频率） ? ??????????????????????????????????????? (5) ??? 当输入为加速度信号，输出为相对位移时，单自由度振动系统的频响函数（增益因子）按下式计算 ? (6) ??? 其中：为自振角频率 ????????? ω为强迫振动角频率 ????????? ξ为阻尼比 对于小阻尼振动系统，在伯德（Bode）图上，频响函数 在 处显示出很尖的峰值。以此为先决条件，反应加速度伪标准差由下式计算（把(6)和(3)式代入(5)整理得到 ? ? ????????????????????? (7) ?? 上式是以平稳随机过程条件的。对比较短的地震过程就不能满足此条件。 ??? Corotic和Vanmarcke用与时间相关的谱密度函数Ｇ(ω,t)来近似计算标准方差 ? (8) ???? 函数 依据于输入信号的谱密度函数 和系统的特性。对线性系统，单位脉冲的频响函数定义为 ? (9) ??? 当 时，函数 收敛于频响函数 。则与时间相关的谱密度函数表示为 ? (10) ??? 对 有如下的近似计算公式 ? (11) ??? 当ｔ远大于 时，上式的积分结果接近于 。把这个值与(11)式比较，假设与时间相关的阻尼定义为 (12) ??? 将 代入(6)式，即得出与时间相关的频响函数的近似计算公式 ? (13) ??? 对于标准差 的系数项 也要进行近似计算。Vanmarcke用下式来计算 ? ? (14) ??? 式中： 为在单自度周期内的振动量 ? 为反应谱谱线 i的角频率 ? 为地震振动周期，取决于包络线的形式和周期 ?? 概率系数(0-1.0) ? ??? 用概率系数可以控制计算过程中得到的反应谱的平均幅值是低于或高于期望反应谱的幅值。当时，相对于期望反应谱的计算反应谱幅值上升和下降的概率是相等的。当概率系数时，实际计算反应谱的幅值低于期望反应谱。当概率系数时，项中的 。 ??? 根据以上公式的推导过程，把(1)和(7)式整理得出地震波功率谱的计算公式 ? (15) ??? 由期望反应谱计算出地震波功率谱以后，用随机信号发生器来生成随机相位，再经富氏逆变换来生成地震波。