動(dòng)力分析入門(mén)：一個(gè)簡(jiǎn)單的動(dòng)力分析模型

數(shù)據(jù)文件：
conf dyn ;ext 5
grid 1 1 50
gen zone copy 15 0 0 ;建立兩個(gè)完全相同的模型
m e
prop shear 1e7 bulk 2e7 dens 1000 ;為了得到剪切波速100，所以shear＝1e7,dens=1000（因?yàn)椋篊s=sqrt(shear/dens)
def wave
if dytime > 1.0/freq
wave = 0.0
else
wave = 0.5 * (1.0 - cos(2.0*pi*freq*dytime))
endif
end
def dyn_time
array vec(3)
vec(1) = 18.0
vec(2) = 0.0
vec(3) = 5.0
oo = set_fontsize(1.5) ;設(shè)置字體大小
oo = draw_string(vec,'Time') ;添加‘Time’
vec(3) = 0.0
oo = draw_string(vec,string(dytime))
vec(3) = 45.0
oo = set_fontsize(1.0)
oo = draw_string(vec,'stress')
vec(1) = 3.0

oo = draw_string(vec,'vel.')
end ;給圖形添加標(biāo)示
set freq=4.0 dyn dt=2e-4 ;設(shè)置dt＝2e-4，表示每cyc 1 個(gè)時(shí)步，實(shí)際時(shí)間為2e-4秒
fix z y range z .5 100
apply xvel 1.0 hist wave yvel 0 zvel 0 range z -.1 .1 ;給模型底部施加剪力波
his gp xvel 0 0 0
his gp xvel 0 0 25
his gp xvel 0 0 50 ;記錄模型底、中、頂部的波速
plo crea qqqj7ZO+X
plo add surf red
plo add vel max 0.5 range x -1 2 ;左邊的模型顯示剪力波
plo add bcon sxz int 1e4 max 1e5 range x 5 18 ; 右邊的模型同步顯示剪力
plo add fish dyn_time black ;添加fish函數(shù)dyn_time設(shè)置的標(biāo)示
plo set per off
plo sho
solve age 2.0 ;算2秒