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

代码如下：set dyn=on large UE8j
8U'L  
ini xvel=0 yvel=0 zvel=0 Znr@-=xZO*  
ini xdis 0 ydis 0 zdis 0 YLJ^R$pi  
free x y z range y -.1 .1 DK)T2{:  
def setup :aQ.:b(n  
omega = 4.0 * pi mC2K &'[  
alfa = 2.2 G?kK:eV  
beta = 3.75e-1 MH=7(15R  
gamma = 8.0 ;NU-\<Q{  
end VWq]w5oQO  
def wave vMd3#@  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) 4>A|2+K\  
end !]5}N^X  
setup !7Eodq-0  
apply xacc=1.0 hist=wave  range  y -.1 .1 &^2SdF  
apply ff ZtyDip'x  
J0V`sK  
hist gp xvel 140,0,0 0civXZgj  
hist gp xvel 200,98,0 0vETg'r  
hist gp yvel 200,98,0 {ET
M >  
hist gp xdis 122,41,0 Lq|>nY  
hist gp xdis 122,41,40 J3`0i@  
hist gp ydis 122,41,0 ZDg(D"  
hist gp ydis 122,41,40 KpA1Ac)T  
hist dytime n{u\t+f  
hist wave B*Q9g r  
set dyn damp local 0.15 ;(approx 5% damping) G %'xEr0n  
solve age 4 %UAF~2]g  
save dam3.sav FA%_jM  
27k(`{K  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 _Jy7` 4B.  

声明下y是竖直方向