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

代码如下：set dyn=on large ?kG#qt]Q5  
ini xvel=0 yvel=0 zvel=0 tOp:e KN  
ini xdis 0 ydis 0 zdis 0 ZKiL-^dob  
free x y z range y -.1 .1 '2i)#~YO<  
def setup !rN#PF>  
omega = 4.0 * pi `t/@ L:  
alfa = 2.2 g6k&c"%IQ(  
beta = 3.75e-1 ]y$V/Ij=qK  
gamma = 8.0 C
>\h?<s  
end >"m@qkh  
def wave pfT`WT  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) bHlDm~5  
end -O5(%  
setup 9, sCJ5bb"  
apply xacc=1.0 hist=wave  range  y -.1 .1 V8|q"UX  
apply ff %s&E-*X  
vi()1LS/!  
hist gp xvel 140,0,0 >V ]*mS%K  
hist gp xvel 200,98,0 }(O D<  
hist gp yvel 200,98,0 HCn]#  
hist gp xdis 122,41,0 NC[GtAPD3  
hist gp xdis 122,41,40 SFXfo1dqH  
hist gp ydis 122,41,0 A(_^_p.|  
hist gp ydis 122,41,40 av|6r#  
hist dytime 1'@lg*^9  
hist wave o 0cc+  
set dyn damp local 0.15 ;(approx 5% damping) H8Z Z@@ qm  
solve age 4 !EyGJa[i  
save dam3.sav yScov)dp(  
F"HI>t)>  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 i32_ZBZ?y  

声明下y是竖直方向