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

代码如下：set dyn=on large ry~3YYEMI0  
ini xvel=0 yvel=0 zvel=0 ( fdDFb#1  
ini xdis 0 ydis 0 zdis 0 ;lYO)Z`3\  
free x y z range y -.1 .1 Mh~T.;f.qq  
def setup }[LK/@h  
omega = 4.0 * pi 3qpkMu3  
alfa = 2.2 Ql
&P1|&  
beta = 3.75e-1 <>j,Q  
gamma = 8.0 x-tA{
_:  
end mG?a)P  
def wave }Q\yem  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime)
-)y"EJ(N  
end ZmmuP/~2K  
setup CvbY2_>Nh  
apply xacc=1.0 hist=wave  range  y -.1 .1 X@~sIUXx9  
apply ff AQQj]7Y  
u52;)"&=)  
hist gp xvel 140,0,0 t&eD;lg :  
hist gp xvel 200,98,0 Z
NCq/  
hist gp yvel 200,98,0 *2
:)Rf  
hist gp xdis 122,41,0 Ae ue:u>  
hist gp xdis 122,41,40 (a^F`#]  
hist gp ydis 122,41,0 Nz!AR$  
hist gp ydis 122,41,40 _F@FcFG1Z*  
hist dytime HowlJ[km%  
hist wave tCc}}2bC&  
set dyn damp local 0.15 ;(approx 5% damping) a#uJzYB0  
solve age 4 QRdb~f;<hj  
save dam3.sav &v 5yo}s  
^
f<f&V  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 5.GBd_;  

声明下y是竖直方向