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

代码如下：set dyn=on large f4lC*nCN  
ini xvel=0 yvel=0 zvel=0 ky{-NrK  
ini xdis 0 ydis 0 zdis 0 k%h%mz  
free x y z range y -.1 .1 #RVN7-x  
def setup ?vocI  
omega = 4.0 * pi vF.Ml  
alfa = 2.2 )jm u*D5N  
beta = 3.75e-1 B$&&'i%  
gamma = 8.0 9p%8VDF=  
end Z)dE#A_X  
def wave )/4U]c{-  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) |s{[<;  
end "Rq)%o$Z  
setup =(]||1.  
apply xacc=1.0 hist=wave  range  y -.1 .1 {U7A&e0eW  
apply ff
'/
GZ,~q  
mqKr+  
hist gp xvel 140,0,0 O`2hTY\  
hist gp xvel 200,98,0 RtScv  
hist gp yvel 200,98,0 #_4JTGJ  
hist gp xdis 122,41,0 BV512+M  
hist gp xdis 122,41,40 =7 ${bp!  
hist gp ydis 122,41,0 N-<m/RS  
hist gp ydis 122,41,40 p'YNj3&u  
hist dytime 3PRK.vf  
hist wave tjt#2i8/  
set dyn damp local 0.15 ;(approx 5% damping) x L]Z3"p%  
solve age 4 {aYCrk1  
save dam3.sav )NCSO b  
/+{1;}AT  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 Qhsk09K_=4  

声明下y是竖直方向