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

代码如下：set dyn=on large a(>oQG8F  
ini xvel=0 yvel=0 zvel=0 k@QU<cvI  
ini xdis 0 ydis 0 zdis 0 
Nm;(M=  
free x y z range y -.1 .1 HPM ggRs  
def setup 6
6S I  
omega = 4.0 * pi )+ }\NCFh  
alfa = 2.2 @&!`.Y oy  
beta = 3.75e-1 2xy &mNx  
gamma = 8.0  -58  
end -P"9KnsO  
def wave \va'>?#o1  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) .Y!]{c  
end sp%EA=: E  
setup C4vmgl&  
apply xacc=1.0 hist=wave  range  y -.1 .1 v5J% p4  
apply ff C>\0 "}iD  
d&mSoPf  
hist gp xvel 140,0,0 yU"pU>fV@  
hist gp xvel 200,98,0 r)S tp`p  
hist gp yvel 200,98,0 ktQMkEj#  
hist gp xdis 122,41,0 8^bc4(H  
hist gp xdis 122,41,40 09FHE/L  
hist gp ydis 122,41,0 K/)*P4C-  
hist gp ydis 122,41,40 "[W${q+0x  
hist dytime 'F3cvpc`  
hist wave |l5ol@2*  
set dyn damp local 0.15 ;(approx 5% damping) ,l}mCY  
solve age 4 AUCk]  
save dam3.sav lP$bxUNt  
ID+o6/V8  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 F$[1KjS  

声明下y是竖直方向