本站小編為你精心準備了非平穩地震動模擬和設計譜擬合參考范文，愿這些范文能點燃您思維的火花，激發您的寫作靈感。歡迎深入閱讀并收藏。

《大連理工大學學報》2016年第三期

摘要：

提出了一種基于小波包分解和重構的模擬非平穩地震動并擬合設計譜的方法．該方法模擬出工程用地震動并調整頻率分量，使其以一定的精度擬合設計譜．首先，利用5個參數生成了構成地震動加速度時程的小波包系數矩陣，并通過小波包重構獲得加速度時程，該系數矩陣模擬了實際地震記錄中頻率隨持時逐漸減小的特點．其次，將模擬的地震動利用小波包分解方法分解為具有高分辨率非重疊的小波包系數矩陣，然后根據設計譜調整頻率分量．數值算例表明，該方法模擬的不同持時的加速度時程均能吻合同一設計譜，且迭代后仍然能保留地震動的非平穩特性．另外，該方法具有較穩定和較快的收斂過程，能在有限次調整迭代中實現較高的擬合精度．

關鍵詞：

地震動模擬；譜兼容加速度時程；小波包；非平穩時程

在工程抗震設計中，反應譜技術自提出以來，一直是抗震設計規范中的重要內容之一［1－3］．而隨著設計方法和計算機技術的發展，人們已經不再滿足于抗震設計的反應譜靜力分析方法，而開始考慮全部地震過程，尤其對于重要工程結構，如核電站、大壩、大型儲物裝置、大跨橋梁、超高層建筑等．輸入結構的地震動時程的反應譜必須在一定精度下與設計加速度反應譜（以下簡稱設計譜）相擬合，也就是所謂的譜兼容時程，通常情況下，擬合方法是對實際強震記錄做修正，讓其滿足設計譜．然而，由于強震記錄的數量有限且分布不均勻，不能滿足多方面要求，地震動的模擬越來越受到重視．一條加速度時程的反應譜是唯一的且是很容易計算的，但是由給定的設計譜來尋求加速度時程這一反問題求解卻是非常困難的．在過去的三十多年里，許多學者對此問題給出了不同的解法．Scanlan等在1974年提出了著名的三角級數模型來模擬具有給定功率譜的地震動加速度時程［4］，這一方法也是工程上應用最廣泛的方法之一．隨后，就開啟了利用傅里葉譜幅值和相位迭代的研究思路．

傅里葉變換是非常直接的反應譜擬合方法，多數學者把重點放在了時域幅值非平穩的模擬上［5－7］．然而該方法存在著局限性：一是沒有很好的收斂特性；二是如果傅里葉譜修正得過多，則加速度時程的非平穩特性將發生較大的改變．之后，演化功率譜方法［8－9］被提出來反映模擬地震動的頻率特性，該方法假定相位和相位差遵循一定的概率分布．從此，頻率非平穩性的模擬受到了學者的重視．Mukherjee等［10］基于小波變換提出了一種直接對實際記錄進行頻率修正來擬合設計譜的方法，從而可以保留實際記錄的大部分頻率特性．但該方法精度不高且沒有給出擬合精度的分析結果．其他方法，如基于諧小波的方法［11］和希爾伯特黃變換法［12］也在保留地震動的瞬時頻率方面取得了較為理想的效果．國內對反應譜擬合技術的研究多采用基于相位譜和相位差譜的方法［13－16］．隨后，趙鳳新等提出了多阻尼反應譜擬合的時域疊加法［17］并引入到小波變換中［18］，其在精度上取得了較好的效果．本文提出一種基于小波包變換的反應譜擬合方法，首先模擬具有時間和頻率非平穩性的地震動；然后采用小波包分解的方法將地震動加速度時程分解為有限頻率段；最后根據設計譜，對各個頻率段進行調整并迭代，從而達到較高的擬合精度．通過對一個核電站的設計譜進行擬合，分別得到兩條不同的地震動加速度時程，并且討論該方法的精度．

1基本原理

小波方法同時在時間域和頻域對信號做分解，從而將頻譜在時間軸上展開．小波包是小波方法的延伸，它不僅對低頻部分做分解，同時也對高頻部分繼續分解，提高了頻率的分辨率．

1．1小波包分解信號的時間和頻率關系叫時頻圖，能夠表征信號這一關系的方法有很多，如短時傅里葉變換、離散與連續小波變換、小波包變換和偽韋格納分布．然而地震動模擬不僅需要對信號進行分解，同時需要能夠方便、完整地逆變換回時域，小波包分解與重構即符合這一要求．

1．2地震動模擬從式（2）可以看出，如果得到了小波包系數矩陣，就可以模擬出一條地震動時程．而如果該小波包系數矩陣擁有與實際記錄相似的特征，則通過小波包重構得到的模擬時程同樣和實際記錄相似．本文首先通過對一條高斯白噪聲進行窗口化，來模擬地震動的時域非平穩性．所用的Saragoni窗函數［19］表達如下。然后，對窗口化的高斯白噪聲做小波包分解．這里采用的模擬時間間隔為0．005s，分解層數j為9，所以分解后的小波包系數頻率分辨率為0．1953Hz，時間分辨率為2．5600s．此時小波包系數在時頻圖的頻率帶上是均勻分布的，需要對每一列系數做修正．圖3（a）展示了美國1992年BigBear地震中SealBeach－OfficeBldg臺站記錄到的一條加速度時程的時頻圖，可以看出其頻率在時域上是不斷變化的．在地震波開始階段，主要包含了高頻分量，隨著時間的進行和地震動幅值的減弱，頻率逐漸衰減．從圖1也可以看出該特性．通過對大量強震記錄的分析，發現大多數地震動均存在該現象．因此，在對每一列小波包系數進行調整時，本文采用了不同的歸一化函數．這里采用對數正態分布密度函數作為基準函數。

1．3設計譜擬合當得到模擬時程x（t）后，重新將其進行小波包分解．該步中，分解層數j增加到12層，這樣頻率分辨率將增大到0．0244Hz，從而保證了反應譜在2～4s周期上的分辨率．分解后得到的小波包系數矩陣c具有2j個行向量，每個行向量cij代表了對應頻率下的能量。

2擬合實例

在實際的設計地震動模擬和擬合中，通常要求地震動的反應譜包絡設計譜．以核電站設計為例［1－2］，在所有控制周期（T）點中，模擬地震動反應譜低于設計譜的點數不超過5個，并且模擬反應譜與設計譜之間的相對誤差應控制在10％以內．根據上述要求，本章在周期0．03～4．00s范圍內選取了76個周期控制點，且設計譜的峰值設置為1g．在該例中，針對設計譜模擬了兩種不同持時的地震動加速度時程，即x1（t）和x2（t）．圖4顯示了迭代后的加速度時程x1（t）和x2（t）以及相應的速度和位移時程．可以看出，加速度時程包含了時間和頻率的非平穩性．x2（t）在整個持時前半程具有相對較高的頻率，而在后半程頻率逐漸降低．由于兩個時程均吻合設計譜，其速度與位移時程也具有相似性．圖5顯示了模擬地震動時程的反應譜（PSA）對設計譜的擬合情況．可以看出，模擬時程的反應譜較好地包絡了設計譜，與設計譜非常接近．由于設定開始迭代的條件為初始模擬出的地震動的反應譜在各個周期上與設計譜誤差在一定范圍內，圖中兩個時程迭代前，其反應譜已經接近設計譜，這是通過設定合理的μwm和σwm來達到的．另外，通過將設計譜譜值增大一定比例的方法，就能達到規范中模擬反應譜低于設計譜點數不大于5的要求．

3結語

本文基于小波包方法提出了一種新的模擬非平穩地震動且擬合設計譜的方法．該方法首先通過模擬小波包系數矩陣來生成地震動，然后根據設計譜來調整小波包系數，用迭代的方法達到模擬地震動的反應譜吻合設計譜的目的．本文提出的方法具有以下特點：（1）本文所模擬的地震動具有實際地震記錄的特點，即具有時間和頻率非平穩性．（2）模擬地震動的過程簡單實用，通過指定5個參數———Td、θ、τ、μwm和σwm，可以模擬出工程需要的地震動．（3）該方法同時包含了隨機性和確定性兩方面．隨機性表現在每次迭代失敗時將重新隨機產生高斯白噪聲；確定性表現在5個初始設定的參數，將限制模擬出的地震動在有限的范圍內變動．（4）該方法收斂過程較快，且比較穩定，擬合設計譜相對誤差較小．

作者：李亞楠 王國新 單位：西北工業大學航海學院