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

代码如下：set dyn=on large F=PBEaX  
ini xvel=0 yvel=0 zvel=0 Ft}nG&D  
ini xdis 0 ydis 0 zdis 0 TH55@1W,[  
free x y z range y -.1 .1 )TBm?VMe  
def setup 79D;0  
omega = 4.0 * pi ~Q]/=HK  
alfa = 2.2 ?T,a(m<i{  
beta = 3.75e-1 xg
tx5tg  
gamma = 8.0 N7GZ'-t^Er  
end 5g;mc.Cvt  
def wave sxdDI?W4  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) N/1xc1$SB  
end ys:F  
setup I|2dV9y  
apply xacc=1.0 hist=wave  range  y -.1 .1 :G=ol2Q  
apply ff E 0OHl  
m:k;?p:x  
hist gp xvel 140,0,0 hpq\  
hist gp xvel 200,98,0 32pPeYxB!-  
hist gp yvel 200,98,0 IR&b2FTcU  
hist gp xdis 122,41,0 1c*:" k  
hist gp xdis 122,41,40 5A%Uv*  
hist gp ydis 122,41,0 wzcv[C-x  
hist gp ydis 122,41,40 6KE64: \;  
hist dytime 2_Zn?#G8dl  
hist wave j'Gezx^.<e  
set dyn damp local 0.15 ;(approx 5% damping) \5a;_N[Ed  
solve age 4 ^* CKx  
save dam3.sav 0d89>UB-8q  
.3SP#mI  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 HIvSh6|0p  

声明下y是竖直方向