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

代码如下：set dyn=on large /Vpd*obMB  
ini xvel=0 yvel=0 zvel=0 Z_%}pe39B  
ini xdis 0 ydis 0 zdis 0 DSwF }  
free x y z range y -.1 .1 ?FRR";  
def setup tVx.J'"Y  
omega = 4.0 * pi >K`.!!av,Y  
alfa = 2.2 '-jKv=D+  
beta = 3.75e-1 %iPu51+=  
gamma = 8.0 Q9K+k*?{N  
end Isq3YY  
def wave _/[n/"gn  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) 'e6W$?z  
end y)3(  
setup `92 D]^g  
apply xacc=1.0 hist=wave  range  y -.1 .1 c$f|a$$b   
apply ff `R@24 )  
lY}mrb  
hist gp xvel 140,0,0 <_*8a(j3  
hist gp xvel 200,98,0 ;WIL?[;w  
hist gp yvel 200,98,0 =q|fe%#  
hist gp xdis 122,41,0 *$(=I6b  
hist gp xdis 122,41,40 p71%-nV  
hist gp ydis 122,41,0 <$liWAGX\  
hist gp ydis 122,41,40 >J=x";,D|~  
hist dytime (PYUfiOf  
hist wave m[^;H
wJ  
set dyn damp local 0.15 ;(approx 5% damping) X.0/F6U  
solve age 4 ,8(%J3J  
save dam3.sav _ED1".&#f  
:,F^{  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 Vvx(7p-GQ  

声明下y是竖直方向