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

代码如下：set dyn=on large "hsT^sy  
ini xvel=0 yvel=0 zvel=0 .e~"+Pe6b  
ini xdis 0 ydis 0 zdis 0 }UhYwJf89  
free x y z range y -.1 .1 5RP kAC  
def setup [8iY0m_Qe  
omega = 4.0 * pi #CC5+  
alfa = 2.2 k;l3^kTy  
beta = 3.75e-1 %j7b0pb  
gamma = 8.0 ]q]xU,  
end n=.P46|  
def wave "n%s>@$  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) sK}AS;:  
end ^#Ii=K-[^  
setup I^
y<W%Et  
apply xacc=1.0 hist=wave  range  y -.1 .1 oW3"J6,S  
apply ff |[*Bn3E:  
Y sM*d  
hist gp xvel 140,0,0 6cH8Jr _  
hist gp xvel 200,98,0 T`&zQQ6F'  
hist gp yvel 200,98,0 /WuYg OI  
hist gp xdis 122,41,0 xlI=)ak{  
hist gp xdis 122,41,40 <Riz!(G  
hist gp ydis 122,41,0 j6m;0
3<|  
hist gp ydis 122,41,40 BQ-x#[%s  
hist dytime 6_j |@  
hist wave &$MC!iMh  
set dyn damp local 0.15 ;(approx 5% damping) yb`PMjj15  
solve age 4 bv+u7B6,  
save dam3.sav (\AN0_  
z#*.9/y\^R  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 :"%/u9<A  

声明下y是竖直方向