地震波模擬的案例
地震波數值模擬技術
地震數值模擬是地震勘探和地震學的重要基礎。地震數值模擬就是在假定地下介質結構模型和相應物理參數已知的情況下，模擬研究地震波在地下各種介質中的傳播 規律，并計算在地面或地下各觀測點所觀測到的數值地震記錄的一種地震模擬方法。這種地震數值模擬方法已在地震勘探和天然地震領域中得到廣泛的應用，它不但 在石油、天然氣、重金屬和非金屬等礦產資源及工程和環境地球物理中得到普遍的應用，而且在地震災害預測、地震區帶劃分以及地殼構造和地球內部結構研究中， 也得到相當廣泛的應用。
地震數值模擬在地震勘探和地震學各工作階段中都有重要的作用。在地震數據采集設計中，地震數值模擬可用于野外觀測系統的設計和評估，并進行地震觀測系統的 優化。在地震數據處理中，地震數值模擬可以檢驗各種反演方法的正確性。在地震數據處理結果的解釋中，地震數值模擬又可以對地震解釋結果的正確性進行檢驗。
由于實際工作中所模擬的介質不同，所用的模擬方程也不一樣。根據模擬方程的不同，波動方程數值模擬主要有：聲波模擬、彈性波模擬、粘彈性波模擬以及裂隙和孔隙彈性模擬等。由于可以用射線理論、積分方程、微分方程來描述地震波的傳播，模擬方法也相應地有射線追蹤法、積分方程數值求解方法以及微分方程數值求解方法。
射線追蹤方法通過求解程函方程計算地震波旅行時，通過求解傳播方程計算地震波振幅。該方法以高頻近似為前提，適合于物性緩變模型中地震波傳播模擬。模型簡 單時該方法具有計算速度快的突出優點，正因為如此，它在地震成像、旅行時層析等方面得到廣泛應用。也正是高頻近似，該方法不適合物性參數變化較大模型中地 震波的傳播模擬。
積分方程數值求解地震波數值模擬方法是基于惠更斯原理而得到的一種波場計算方法，它又可以分為體積分方法和邊界積分方法。
展開 PEER下載地震波方法介紹
在我們做抗震研究時，總是離不開地震波的數據。在進行有限元模擬時，我們往往需要輸入地震波來模擬構件在地震作用下的破壞。PEER（太平洋地震工程研究中心）是一個非常友好的網站，上面提供了很多地區地震時的記錄。那么怎么在peer上下載地震波呢？
（1）賬號注冊及授權登錄
直接輸入peer網址
https://peer.berkeley.edu/。網頁界面如下圖所示：
首次進入時自己需要使用郵箱注冊賬號，點擊注冊入口。填寫相關信息（注冊信息要記錄下來）
注冊好以后，還不能登陸使用，網站會向剛才注冊的郵箱發送授權確認的請求，收到郵件以后確認授權即可，再次進入網站進行登陸。
賬號密碼登錄界面
點擊NGA West2進入
（2）地震波下載
進去以后網站提示要選擇獲取目標譜的方法，這里采用自定義譜的方法。
選擇自定義譜方法
需要上傳反應譜曲線數據
注：反應譜數據可以使用相應軟件獲取。比如崔繼東博士的小軟件。這個在這里不再講解，大家自行學習。
我們選擇相應的反應譜曲線
確認選擇好文件以后，點擊上傳。國外網站有時候進入較慢，等待片刻，顯示上傳成功，點擊下面的提交選項。提交文件，等待網站識別上傳的文件。
在等待一些時間之后，就可以看到相應的反應譜曲線。
但是這里采用的是對數坐標，選擇Loglog下面的下拉箭頭，將對數坐標系改成線性Liner坐標系，這樣將顯示與小軟件中相同的反應譜曲線。
點擊下面的搜索地震記錄，查找與規范反應譜相吻合的地震名稱。
頁面跳轉，繼續點擊搜索，之后會顯示查找到100條的地震記錄。或者可以在選項中填入具體的參數進行細致的篩選。
展開 結構地震波作用下的頻率響應模擬 ￥400
本案例仿真了一結構在一側受到低頻振動作用下的頻率響應結果，如圖1所示。
圖1 仿真結果
一文讀懂地震波 附Ei-center地震波下載
因為面波隨著距離衰減的速率比P波或S波慢，在距震源距離大時感知的或長時間記錄下來的主要是面波。
最后 在面波波列之后構成地震記錄的重要部分，稱之為地震尾波。地震波的尾部事實上包含著沿散射的路徑穿過復雜巖石構造的P波、S波、勒夫波和瑞利波的混合波。尾波中繼續的波動旋回對于建筑物的破壞可能起到落井下石的作用，促使已被早期到達的較強S波和面波削弱的建筑物進一步破壞以至于倒塌
下載地址：Ei-center地震波
展開 
【JY】建筑結構施加地震波的方法與理論機理
最近有人問在軟件中施加地震波應該怎么加？
其實，在對結構模型進行地震時程分析時，最常見的方式有三種方式。
地震波的輸入有三種方式：
第一種是是質量加速度施加法，通過達朗貝爾原理，將地震作用轉化為施加在質點上慣性力；
第二種是底部位移法，在結構底部直接輸入位移地震波，模擬地面震動反應，計算結構響應；
第三種是底部加速度法
，將地震加速度從結構模型底部輸入。
在大部分軟件中，采用的是第一、二種方式進行加載，而對于少部分軟件也可采用第三種方式進行施加，具體如圖所示。
二、
案例詳述
對于施加加速度或者位移，通常我們可以收集得到相關的地震動加速度，首先將地震動加速度進行轉換得到位移，可采用
自編的SignalData軟件進行轉換。
詳情的安裝和使用方法請點擊以下文章
【JY】SignalData軟件開發應用分享
1、 在Sap2000中：
SAP2000軟件采用兩種施加方式，分別為：
質量加速度施加法(在結構質量上施加加速度荷載)；
底部位移法(以地面位移方式施加與結構底部)。
展開 六大常用地震數據庫，地震波庫特點及選波建議
近年的中強震包括 2003 年伊朗 Bam地震（6.6）、2004 年加利福尼亞 Parkfield 地震（M7.2）、2010 年新西蘭 Darfield 地震（M7.0）和 2011 年新西蘭 Christchurch 地震（ M6.2）。斷層距的分布從 0.05 km 至 1533 km不等，常用數據基本分布于 400 km 以內，共涉及 4230個臺站。其中有大量的常用地震波，如 1940 年帝谷地震獲得的 El Centro 波。從 NGA-West2可以得到其相關信息，包括該地震的發生時間、地點、震級、地震記錄的距離信息、反應譜、場地條件等。
NGA-West2 數據庫的數據具有極高的格式一致性，所有的數據都具有相同的格式。加速度時程記錄的文件后綴為“.AT2”。其數據格式均為 4 行表頭和加速度時程記錄的組合。其數據格式的一致性對于大批量數據的處理非常便利。4 行表頭包含地震名稱、地震發生時間、臺站名稱、記錄分量的角度、時程記錄的單位、總記錄數據的數量和采樣時間間隔。PEER 數據庫中的數據均為未調幅（un- scaled）和真實的（as-record）未經旋轉（unrotated）的數據記錄。但是數據已經經過了校正，包括基線初始化、濾波等，可以直接使用 PEER提供的加速度數據得到的速度、位移數據。但 PEER 記錄的事前和事后時間并不固定， 有時甚至沒有事前和事后。此數據目前能夠在 PEER 的網站上的 NGA-West2 數據庫獲得，并且可以通過自行指定地震事件、記錄編號、各地震動參數區間等篩選條件挑選合適的記錄進行下載。
展開 EI 地震波數據 0.02
EI地震波.xls
人工地震波的生成
請問大家simqke-gr這個軟件怎么安裝啊。安裝包下載下來卻安裝不上。
地震波加載問題
得到的地震波有 up ew ns 三種，這三種對應的面波和體波的對應關系是怎樣的？還有在 ABAQUS 中地層中這三個波都要加載還是只加載一種
如何利用反應譜生成地震波
NCYCLE：迭代的次數+1，如果 NCYCLE=1 則不進行迭代，使得在計算過
程中得到的反應譜曲線更加平滑
AGMAX：峰值加速度，由程序自動計算(這里的單位是g，跟我們輸入的是一樣的)
NPA：要生成的時程波條數，這里選擇 3 條
IIX：任意奇數，缺省值為 1235
AMOR：反應譜的阻尼比，這里選擇 0.05
===============================================================
然后點擊Simulated Earthquake,生成地震波如下圖所示：
查看相應的數據：
導出的數據如下所示：
這期的利用反應譜生成地震波的教程就結束了，這個軟件網上可以直接下載，也可以關注公眾號“生活中的力學仿真”，后臺回復“地震波生成軟件”獲取。
展開 如何確定所施加地震波的方向
在ANSYS中導入地震波（txt文件，第一列時間，第二列加速度），怎么看地震波的施加方向啊？（水平面為x和y方向）
求汶川地震波數據
哪位大哥有地震波數據，小弟需要汶川波數據，有償
PFC導入地震波加載
邊坡的反應并不是這篇文章的重點，提出的這種加載地震波的方式是可以參考一下的，感興趣的話積分方法可以采用數值分析中介紹的方法去優化。
ansys之——地震波的輸入和求解
對于地震波的輸入，可以把荷載記錄做成文件，利用apdl的讀取功能讀入倒數據庫中。下面的例子是自己編的一個小文件。修改一下可以更簡潔。有用到的朋友自己作一下把。
fini
/config,nres,1000
*dim,aceX,TABLE,3000,1
*dim,aceY,TABLE,3000,1
*dim,aceZ,TABLE,3000,1
*creat,ff
*vread,aceX(1,1),acex,txt,,1
(e16.6)
*vread,aceX(1,0),ACETT,,,1
(e17.6)
ACEX(0,1)=1
*end
/input,ff
*creat,ff
*vread,aceY(1,1),acey,txt,,1
(e16.6)
*vread,aceY(1,0),ACETT,,,1
(e17.6)
ACEY(0,1)=1
*end
/input,ff
*creat,ff
*vread,aceZ(1,1),acez,txt,,1
(e16.6)
*vread,aceZ(1,0),ACETT,,,1
(e17.6)
ACEZ(0,1)=1
*end
/input,ff
！地震波時程記錄分成了3個文件，每個文件是一列。分別記錄x,y,z方向的加速度。acett是時間記錄。
這樣就可以把加速度記錄讀取倒ansys數據庫中作為數組。
也可以把加速度記錄做成一個文件，這樣程序就簡單多了。大家可以試看看修改一下。
下面是計算部分語句：
/SOLU
ANTYPE,trans
！
展開 地震波的90%的python腳本持時計算 ￥2
話不多說，先展示計算后的結果示意圖，包括下面給的代碼，運行計算后都會給出相應的持時圖！
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
# 設置中文字體和負號顯示
plt.rcParams['font.sans-serif'] = ['SimHei'] # 顯示中文
plt.rcParams['axes.unicode_minus'] = False # 正常顯示負號
# 讀取 Excel 文件
file_path = 'x.xlsx' # 替換為你的文件路徑
df = pd.read_excel(file_path, engine='openpyxl')
# 提取時間和加速度列
time = df['Time'].values
acc = df['Acceleration'].values
# 檢查數據采樣是否均勻
dt = np.diff(time)
if not np.allclose(dt, dt[0], atol=1e-5):
raise ValueError("時間列采樣不均勻，需預處理為均勻采樣數據")
dt = dt[0] # 時間步長
展開 