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

代码如下：set dyn=on large -{L[Wt{1  
ini xvel=0 yvel=0 zvel=0 Wb*d`hzQ}  
ini xdis 0 ydis 0 zdis 0 p=7kFv  
free x y z range y -.1 .1 #Va@4<4r  
def setup -4hX-  
omega = 4.0 * pi :&TOQ<vM  
alfa = 2.2 &1
B)mj  
beta = 3.75e-1 k#&y
  
gamma = 8.0 b(~NqV!i  
end x%x[5.CT  
def wave 6Ajiz_~U  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime)
40q8,M  
end ?<F\S2W

 
setup %AW  
apply xacc=1.0 hist=wave  range  y -.1 .1 g<.VW0  
apply ff #j;&g1  
Q) FL|  
hist gp xvel 140,0,0 pu`|HaQaE  
hist gp xvel 200,98,0 g7d)YUc  
hist gp yvel 200,98,0 y f+/Kj< a  
hist gp xdis 122,41,0 ^QFjBQ-Hai  
hist gp xdis 122,41,40 6dg[  
hist gp ydis 122,41,0 t3bDi/m  
hist gp ydis 122,41,40 NrL%]dl3/  
hist dytime jRq>Sz{8  
hist wave a(BC(^1!  
set dyn damp local 0.15 ;(approx 5% damping) "=/XIM.  
solve age 4 QiwZk<rb  
save dam3.sav "$Rl9(}  
eKLxNw5  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 lWOB!l  

声明下y是竖直方向