地震波ansys模擬的案例
地震波數(shù)值模擬技術(shù)
地震數(shù)值模擬是地震勘探和地震學(xué)的重要基礎(chǔ)。地震數(shù)值模擬就是在假定地下介質(zhì)結(jié)構(gòu)模型和相應(yīng)物理參數(shù)已知的情況下，模擬研究地震波在地下各種介質(zhì)中的傳播 規(guī)律，并計(jì)算在地面或地下各觀測點(diǎn)所觀測到的數(shù)值地震記錄的一種地震模擬方法。這種地震數(shù)值模擬方法已在地震勘探和天然地震領(lǐng)域中得到廣泛的應(yīng)用，它不但 在石油、天然氣、重金屬和非金屬等礦產(chǎn)資源及工程和環(huán)境地球物理中得到普遍的應(yīng)用，而且在地震災(zāi)害預(yù)測、地震區(qū)帶劃分以及地殼構(gòu)造和地球內(nèi)部結(jié)構(gòu)研究中， 也得到相當(dāng)廣泛的應(yīng)用。
地震數(shù)值模擬在地震勘探和地震學(xué)各工作階段中都有重要的作用。在地震數(shù)據(jù)采集設(shè)計(jì)中，地震數(shù)值模擬可用于野外觀測系統(tǒng)的設(shè)計(jì)和評(píng)估，并進(jìn)行地震觀測系統(tǒng)的 優(yōu)化。在地震數(shù)據(jù)處理中，地震數(shù)值模擬可以檢驗(yàn)各種反演方法的正確性。在地震數(shù)據(jù)處理結(jié)果的解釋中，地震數(shù)值模擬又可以對地震解釋結(jié)果的正確性進(jìn)行檢驗(yàn)。
由于實(shí)際工作中所模擬的介質(zhì)不同，所用的模擬方程也不一樣。根據(jù)模擬方程的不同，波動(dòng)方程數(shù)值模擬主要有：聲波模擬、彈性波模擬、粘彈性波模擬以及裂隙和孔隙彈性模擬等。由于可以用射線理論、積分方程、微分方程來描述地震波的傳播，模擬方法也相應(yīng)地有射線追蹤法、積分方程數(shù)值求解方法以及微分方程數(shù)值求解方法。
射線追蹤方法通過求解程函方程計(jì)算地震波旅行時(shí)，通過求解傳播方程計(jì)算地震波振幅。該方法以高頻近似為前提，適合于物性緩變模型中地震波傳播模擬。模型簡 單時(shí)該方法具有計(jì)算速度快的突出優(yōu)點(diǎn)，正因?yàn)槿绱耍?em>地震成像、旅行時(shí)層析等方面得到廣泛應(yīng)用。也正是高頻近似，該方法不適合物性參數(shù)變化較大模型中地 震波的傳播模擬。
積分方程數(shù)值求解地震波數(shù)值模擬方法是基于惠更斯原理而得到的一種波場計(jì)算方法，它又可以分為體積分方法和邊界積分方法。
展開 結(jié)構(gòu)地震波作用下的頻率響應(yīng)模擬 ￥400
本案例仿真了一結(jié)構(gòu)在一側(cè)受到低頻振動(dòng)作用下的頻率響應(yīng)結(jié)果，如圖1所示。
圖1 仿真結(jié)果
ansys之——地震波的輸入和求解
對于地震波的輸入，可以把荷載記錄做成文件，利用apdl的讀取功能讀入倒數(shù)據(jù)庫中。下面的例子是自己編的一個(gè)小文件。修改一下可以更簡潔。有用到的朋友自己作一下把。
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
！地震波時(shí)程記錄分成了3個(gè)文件，每個(gè)文件是一列。分別記錄x,y,z方向的加速度。acett是時(shí)間記錄。
這樣就可以把加速度記錄讀取倒ansys數(shù)據(jù)庫中作為數(shù)組。
也可以把加速度記錄做成一個(gè)文件，這樣程序就簡單多了。大家可以試看看修改一下。
下面是計(jì)算部分語句：
/SOLU
ANTYPE,trans
！
展開 