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

代码如下：set dyn=on large a<]vHC7  
ini xvel=0 yvel=0 zvel=0 eYn/F~5-  
ini xdis 0 ydis 0 zdis 0 5E 9
R+N  
free x y z range y -.1 .1 6OJhF7\0&  
def setup xX0wn?,~  
omega = 4.0 * pi #s#BYbF  
alfa = 2.2 {iCX?Sb  
beta = 3.75e-1 *5\'$;Rg  
gamma = 8.0 E$A=*-u  
end hQz1zG`z7  
def wave GuaF B[4  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) =s*4y$%
I  
end ({$rb-  
setup DGw*
BN%`  
apply xacc=1.0 hist=wave  range  y -.1 .1 Lo Y*,Aa&  
apply ff }IdkXAB.  
(=Oo=8\  
hist gp xvel 140,0,0 +.rE|)BPy  
hist gp xvel 200,98,0 _;/+8=  
hist gp yvel 200,98,0 -G#m'W&  
hist gp xdis 122,41,0 (]VY==t~  
hist gp xdis 122,41,40 Y4lNxvY  
hist gp ydis 122,41,0 0bD\`Jiv,  
hist gp ydis 122,41,40 <]_[o:nOP  
hist dytime `,Y3(=3Xe?  
hist wave ^rO!-  
set dyn damp local 0.15 ;(approx 5% damping) rmFcSolt,f  
solve age 4 snNB;hkj  
save dam3.sav 0-uVmlk=/  
;TK$?hrv*1  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 iSfRo31  

声明下y是竖直方向