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

代码如下：set dyn=on large y;0Zk~R$  
ini xvel=0 yvel=0 zvel=0 HJr*\%D}1  
ini xdis 0 ydis 0 zdis 0 =5:vKL j  
free x y z range y -.1 .1 ;[zZI~wh  
def setup XW'7  
omega = 4.0 * pi :A#'8xE/  
alfa = 2.2 5m;pHgkb  
beta = 3.75e-1 )p!") :'fv  
gamma = 8.0 ;J?fK69%  
end aorL,l  
def wave q89yW)XG  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) 7F~Jz*,B*W  
end qq^[(n  
setup 3F6'3NvVc2  
apply xacc=1.0 hist=wave  range  y -.1 .1 Q0PqyobD  
apply ff r
I)&.5^  
j%vxCs>  
hist gp xvel 140,0,0 ,0L< w
a  
hist gp xvel 200,98,0 z:n JN%Qb  
hist gp yvel 200,98,0 EZQ!~  
hist gp xdis 122,41,0 ^')4R
U  
hist gp xdis 122,41,40 Pi&\GMzd  
hist gp ydis 122,41,0 54JI/!a  
hist gp ydis 122,41,40 PHsM)V+  
hist dytime <!u(_Bxw/  
hist wave ^[HX#JJ~  
set dyn damp local 0.15 ;(approx 5% damping) r|jBKq~  
solve age 4 { F0"U=  
save dam3.sav L8;`*H  
]~WIGl"g  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 LQ$dT#z2A  

声明下y是竖直方向