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

代码如下：set dyn=on large k<Gmb~Tg1  
ini xvel=0 yvel=0 zvel=0 j3;W-c`5  
ini xdis 0 ydis 0 zdis 0 MnqT?Cc4$j  
free x y z range y -.1 .1 BKtb@o~(  
def setup _q#pEv  
omega = 4.0 * pi aT IzfqCM  
alfa = 2.2 GvG8s6IZ  
beta = 3.75e-1 No6-i{HZ  
gamma = 8.0 L~{(9J'(  
end (wq8[1Wzup  
def wave w^ z ftm  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) "5<!   
end SyI~iW#Y1  
setup ><D2of|  
apply xacc=1.0 hist=wave  range  y -.1 .1 Qt{){u
E  
apply ff  B>:U  
bAH<h   
hist gp xvel 140,0,0 -K?lhu  
hist gp xvel 200,98,0 YcX"Z~O6j=  
hist gp yvel 200,98,0 ^*`#+*C  
hist gp xdis 122,41,0 au+Jz_$)  
hist gp xdis 122,41,40 5NUaXQ  
hist gp ydis 122,41,0 ,'HjL:r  
hist gp ydis 122,41,40 |yO%w#  
hist dytime RHn3\N  
hist wave /eH37H  
set dyn damp local 0.15 ;(approx 5% damping) \P+^BG!  
solve age 4 ]tT=jN&(  
save dam3.sav ]  &"`  
y[85eM  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 ,
(Zxd4?y  

声明下y是竖直方向