[热点讨论]free-field边界问题 [复制链接]

代码如下：set dyn=on large y$7@~NH,d  
ini xvel=0 yvel=0 zvel=0 2\1\Jn#q  
ini xdis 0 ydis 0 zdis 0 M"$TXXe  
free x y z range y -.1 .1 w=H4#a?fc  
def setup .`Ts'0vVy  
omega = 4.0 * pi )QiHe}  
alfa = 2.2 V7 dAB,:  
beta = 3.75e-1 t?j2Rw3f`I  
gamma = 8.0 Cb/?hT  
end ALY% h!L  
def wave 7f`x-iH!]7  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) r]vD]  
end Q|W!m0XO  
setup *cn#
W]AE  
apply xacc=1.0 hist=wave  range  y -.1 .1 sw@*N  
apply ff Bt^];DjH  
R(sa.Q\D4  
hist gp xvel 140,0,0 _*w}"\4_  
hist gp xvel 200,98,0 %1p4K)  
hist gp yvel 200,98,0 Fttny]  
hist gp xdis 122,41,0 fMFlY%@t  
hist gp xdis 122,41,40 P!)k4n  
hist gp ydis 122,41,0 f3]u-e'b  
hist gp ydis 122,41,40 ,.+"10=N.  
hist dytime k^PqB+P!  
hist wave aE2 3[So  
set dyn damp local 0.15 ;(approx 5% damping) XT5Vo  
solve age 4 7J%v""\1!  
save dam3.sav W<uL{k.Kpd  
WsCzC_'j.  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 tKUy&]T  

声明下y是竖直方向