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

代码如下：set dyn=on large h^kNM8  
ini xvel=0 yvel=0 zvel=0 e
.|RC  
ini xdis 0 ydis 0 zdis 0 Pl&`&N;  
free x y z range y -.1 .1 iCRw}[[  
def setup =v$s+`cP  
omega = 4.0 * pi IC.<)
I  
alfa = 2.2 zy6(S_j  
beta = 3.75e-1 ESjJHZoD(  
gamma = 8.0 8:?Q(M7  
end cqL7dlhIl  
def wave sJK:xk.6!  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) ."Ix#\|x  
end (Zg'pSs)  
setup IPi<sE  
apply xacc=1.0 hist=wave  range  y -.1 .1 Fi%W\Y'  
apply ff p]z54 ~  
~Z6p3# !o  
hist gp xvel 140,0,0 /3Ix,7  
hist gp xvel 200,98,0 XqS*;Zj0  
hist gp yvel 200,98,0 DPQGh`J  
hist gp xdis 122,41,0 Ty0T7D  
hist gp xdis 122,41,40 np\2sa`  
hist gp ydis 122,41,0 XW_xNkpL5c  
hist gp ydis 122,41,40 *M<BPxh0w]  
hist dytime 8t:&#h  
hist wave wx%nTf/Oa  
set dyn damp local 0.15 ;(approx 5% damping) 0$Y 9>)O  
solve age 4 ^@lg5d3F  
save dam3.sav a& aPBv1  
'oZn<c`  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 >"g<-!p@  

声明下y是竖直方向