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

代码如下：set dyn=on large ; k.
@=  
ini xvel=0 yvel=0 zvel=0 $Q#n'#c  
ini xdis 0 ydis 0 zdis 0 vS~tr sI  
free x y z range y -.1 .1 -j<E_!t  
def setup =dNE1rdzNa  
omega = 4.0 * pi s) s9Z,HY  
alfa = 2.2 jS
vq1$U  
beta = 3.75e-1 uVD^X*  
gamma = 8.0 f:
\)! &W  
end `)8~/G%  
def wave In?+  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) BV@q@C  
end HRE?uBkjf  
setup ~;$QSO\2h  
apply xacc=1.0 hist=wave  range  y -.1 .1 dh6kj-^;Cf  
apply ff L3oL>r'|  
X);'[/]E*  
hist gp xvel 140,0,0 3xz~##  
hist gp xvel 200,98,0 V*rLGY#  
hist gp yvel 200,98,0 W"@'}y  
hist gp xdis 122,41,0 {,Vvm*L/  
hist gp xdis 122,41,40 jq]5Y^e  
hist gp ydis 122,41,0 Co
nik`  
hist gp ydis 122,41,40 j"~"-E(79  
hist dytime =\2gnk~  
hist wave ~{{S<S v  
set dyn damp local 0.15 ;(approx 5% damping) ^wZx=kas  
solve age 4 B<BS^waU  
save dam3.sav LbJf5xdi  
0/DO"pnL@  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 2Cy,#X%j>  

声明下y是竖直方向