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

代码如下：set dyn=on large F *FwRj  
ini xvel=0 yvel=0 zvel=0 )M6w5g  
ini xdis 0 ydis 0 zdis 0 <:SZAAoIV  
free x y z range y -.1 .1 &wD;SMr<  
def setup P:30L'.=[  
omega = 4.0 * pi S;$-''o?9  
alfa = 2.2 =|3fs7  
beta = 3.75e-1 Fa'k0/_j  
gamma = 8.0 h4q|lA6!k8  
end S S2FTb-m  
def wave &|f@$ff  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) q#<^^4U  
end j@YU|-\qh  
setup ]]Cb$$Td  
apply xacc=1.0 hist=wave  range  y -.1 .1 B!vmQR*1  
apply ff M$Zcn#A  
0TQ$C-%  
hist gp xvel 140,0,0 Z.92y  
hist gp xvel 200,98,0 su*'d:L  
hist gp yvel 200,98,0 S H!  
hist gp xdis 122,41,0 87D*-Gw  
hist gp xdis 122,41,40 N[s}qmPha  
hist gp ydis 122,41,0 ^$b Y,CE  
hist gp ydis 122,41,40 {q"OM*L(  
hist dytime !o:f$6EA~C  
hist wave `DV.+>O-1  
set dyn damp local 0.15 ;(approx 5% damping) 3AU;>D^5  
solve age 4 9I6
a"PGDb  
save dam3.sav :]\([Q+a  
:Zlwp6  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 i\,-oO  

声明下y是竖直方向