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

代码如下：set dyn=on large ljPq2v ]  
ini xvel=0 yvel=0 zvel=0 Q,.By&  
ini xdis 0 ydis 0 zdis 0 HG2GZ}~^1  
free x y z range y -.1 .1 *hVb5CS  
def setup 5YiZ-CQ>  
omega = 4.0 * pi x;]x_fz  
alfa = 2.2 [pii  
beta = 3.75e-1 &%^K,Q"  
gamma = 8.0 3Y z]8`C  
end XOrcygb2  
def wave 5W+{U8\  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) akT|Y4KxD  
end XGfzEld2"  
setup s^w
\zzYb  
apply xacc=1.0 hist=wave  range  y -.1 .1 D_d|=
i  
apply ff _(C^[:s  
Ic'Q5kfM  
hist gp xvel 140,0,0 ITyzs4"VV  
hist gp xvel 200,98,0 -(*nSD9  
hist gp yvel 200,98,0 sVu k  
hist gp xdis 122,41,0 vwKw?Z0%J  
hist gp xdis 122,41,40 .H8mRvd?  
hist gp ydis 122,41,0 v f`9*xF  
hist gp ydis 122,41,40 L=,OZ9aA  
hist dytime P##Z[$IJ3  
hist wave }YQ:6I  
set dyn damp local 0.15 ;(approx 5% damping) .APVjqG  
solve age 4 0 y<k][  
save dam3.sav K@HQrv<  
Wo{K}  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 +W9]ED  

声明下y是竖直方向