地震波動(dòng)數(shù)值模擬及穩(wěn)定性研究.pdf

1、本文針對(duì)顯式有限元-有限差分法進(jìn)行了系統(tǒng)的研究,包括理論分析,算法改進(jìn)以及實(shí)際應(yīng)用。通過理論推導(dǎo)定量地研究了采用這種數(shù)值模擬方法對(duì)波動(dòng)傳播規(guī)律的影響,以及其在結(jié)合局部透射人工邊界進(jìn)行計(jì)算時(shí)的穩(wěn)定條件。分別就采用該方法結(jié)合局部透射人工邊界求解波源問題和散射問題的情況分別提出了改進(jìn)計(jì)算的措施,并通過數(shù)值試驗(yàn)驗(yàn)證其有效性，提出了一種有限元和顯式有限元-有限差分法耦合方法,并采用其對(duì)可控震源激震過程進(jìn)行了數(shù)值模擬。利用在汶川地震中獲得的數(shù)據(jù),進(jìn)

2、行場地效應(yīng)分析。主要工作如下：
　　 ⑴根據(jù)由李小軍等人提出的一種求解動(dòng)力方程的顯式積分格式應(yīng)用于顯式有限元-有限差分法進(jìn)行波動(dòng)數(shù)值模擬時(shí)的遞推形式,給出了波動(dòng)傳遞函數(shù)的形式；以一維均勻離散體系為模型,推導(dǎo)了該情況下數(shù)值波動(dòng)的傳遞函數(shù),并進(jìn)行數(shù)值試驗(yàn),驗(yàn)證了傳遞函數(shù)求解的正確性；并通過分析求解所得的傳遞函數(shù),探討了這種顯式積分格式對(duì)離散網(wǎng)格中波動(dòng)傳播規(guī)律包括截止頻率、頻散現(xiàn)象和能量耗散等方面的影響。
　　 ⑵針對(duì)顯式有限

3、元-有限差分法結(jié)合局部透射邊界進(jìn)行波動(dòng)數(shù)值模擬的情況,推導(dǎo)并計(jì)算了循環(huán)系數(shù),給出了計(jì)算過程的穩(wěn)定條件,并通過數(shù)值試驗(yàn)對(duì)該判定準(zhǔn)則進(jìn)行驗(yàn)證；在給出的循環(huán)系數(shù)的基礎(chǔ)上,定量地分析了積分格式的能耗特性對(duì)由邊界條件引入的高頻失穩(wěn)的抑制和消除作用；并探討了如何更有效地利用能耗特性保證計(jì)算穩(wěn)定性,即實(shí)際計(jì)算模型中時(shí)間離散步距以及空間模型大小等參數(shù)如何取值既能確保計(jì)算穩(wěn)定又能最大程度減小計(jì)算量。
　　 ⑶分別針對(duì)采用顯式有限元-有限差分法結(jié)合

4、局部透射人工邊界求解波源問題和散射問題的情況分別提出了改進(jìn)措施。對(duì)于波源問題,建議針對(duì)每個(gè)邊界節(jié)點(diǎn)可以采用不同的人工波速,該人工波速可以通過試算一定時(shí)間長度,然后通過相關(guān)處理獲得。通過數(shù)值實(shí)驗(yàn)證明了,在波源問題的求解中采用這種方法選擇人工波速可以在不增加計(jì)算失穩(wěn)風(fēng)險(xiǎn)的同時(shí)提高數(shù)值計(jì)算精度。對(duì)于散射問題,則給出了一種關(guān)于波動(dòng)輸入方法的改進(jìn)措施,即在分區(qū)分離波場中入射波不再采用解析解而是采用離散數(shù)值解代替,數(shù)值實(shí)驗(yàn)證明了采用該方法可以有效地

5、抑制漂移失穩(wěn)。
　　 ⑷給出了一種基于有限元和顯式有限元.有限差分法相耦合的區(qū)域分裂法。將計(jì)算模型分為兩個(gè)區(qū)域,對(duì)于廣義結(jié)構(gòu)及其附近區(qū)域采用有限元法模擬,而土體其它部分采用顯式有限元-有限差分法,兩個(gè)區(qū)域的模擬通過搭接區(qū)域?qū)崿F(xiàn)數(shù)據(jù)的交換。分別代表波源問題和散射問題的數(shù)值試驗(yàn)都驗(yàn)證了該方法的有效性。采用提出的耦合方法針對(duì)可控震源激振過程建立的數(shù)值計(jì)算模型,并對(duì)不同的工況進(jìn)行計(jì)算,研究分析了各因素包括接觸土體密度、接觸土體剛度、接觸

6、土層厚度、軟、硬夾層以及底平板半徑對(duì)可控震源能量輸入大地體系的影響,為如何能高效地將激震能量輸入土體提供了參考。
　　 ⑸利用自貢地形影響臺(tái)陣在2008年汶川地震中獲得的主震強(qiáng)震記錄,采用水平/垂直譜比法(HVSR)對(duì)該山脊地形場地對(duì)地震動(dòng)的放大效應(yīng)進(jìn)行了分析,并探討了所獲記錄結(jié)果的可靠性；采用顯式有限元-有限差分法結(jié)合局部透射人工邊界的二維場地影響分析模型與方法,模擬了該場地的地震動(dòng)響應(yīng)。對(duì)模擬結(jié)果與HVSR法分析結(jié)果進(jìn)行對(duì)比