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

代码如下：set dyn=on large W~d^ *LZt  
ini xvel=0 yvel=0 zvel=0 U}^`R,C  
ini xdis 0 ydis 0 zdis 0 Pm==m9  
free x y z range y -.1 .1 )UN_,'H/V  
def setup zp:EssO=Q  
omega = 4.0 * pi R-OQ(]<*  
alfa = 2.2 "eZ~]m}L0  
beta = 3.75e-1 0+6=ag%  
gamma = 8.0 UB3hC`N\  
end @\|Fd)  
def wave cs7K^D;.V  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) Wz)@k2  
end G}#p4\/  
setup HA]5:ck  
apply xacc=1.0 hist=wave  range  y -.1 .1 :[!b";pR  
apply ff T/iZ"\(~w  

/5"RedP<  
hist gp xvel 140,0,0 )kvrQ6  
hist gp xvel 200,98,0 NXSjN~aG2  
hist gp yvel 200,98,0 `tH:oP0=  
hist gp xdis 122,41,0 (=t41-l  
hist gp xdis 122,41,40 A!IZIT5)m  
hist gp ydis 122,41,0 ;>r E+k%_  
hist gp ydis 122,41,40 y[d>7fc
f  
hist dytime p}(pIoyUF  
hist wave KkyZd9  
set dyn damp local 0.15 ;(approx 5% damping) ZfnJ&H'  
solve age 4 'QQa :3<x  
save dam3.sav JlhI3`X;/  
W
WN2  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 uh&Qdy!I  

声明下y是竖直方向