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

代码如下：set dyn=on large Pi6C/$ K  
ini xvel=0 yvel=0 zvel=0 GlYly5F  
ini xdis 0 ydis 0 zdis 0 3Kq`<B~%  
free x y z range y -.1 .1 \{|ImCH  
def setup r#876.JK  
omega = 4.0 * pi w<wV]F*  
alfa = 2.2 Q4'C;<\@(Q  
beta = 3.75e-1 dDcZ!rRaL@  
gamma = 8.0 =yiOJyx  
end 7qIB7_K5  
def wave ]?l{j  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) 0%C^8%(x  
end A*]$
v  
setup HOW7cV'X  
apply xacc=1.0 hist=wave  range  y -.1 .1 .J.vC1 4gi  
apply ff >l1r,/\\
 
NlPS#  
hist gp xvel 140,0,0 *:"^[Ckc  
hist gp xvel 200,98,0 }p3b#fAr  
hist gp yvel 200,98,0 j BS4vvX?  
hist gp xdis 122,41,0 %e%7oqR?  
hist gp xdis 122,41,40 *> 3Qd7  
hist gp ydis 122,41,0 'ig, ATY  
hist gp ydis 122,41,40 v"lf-c  
hist dytime D,;\F,p  
hist wave ]KK`5Dv|,e  
set dyn damp local 0.15 ;(approx 5% damping) +&v\ /  
solve age 4 f?UzD#50D  
save dam3.sav L10IF  
%_)zWlN  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 [s6C ZcL  

声明下y是竖直方向