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

代码如下：set dyn=on large !do?~$Og  
ini xvel=0 yvel=0 zvel=0 pH@]Y+W  
ini xdis 0 ydis 0 zdis 0 !w39FfU{  
free x y z range y -.1 .1 X4|4QgY  
def setup \ZI'|Ad  
omega = 4.0 * pi x=q;O+7]  
alfa = 2.2 ;# uZhd  
beta = 3.75e-1 -0C@hM,wm  
gamma = 8.0 BU|#e5  
end @-&MA
)SN  
def wave HKDID[d0  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) lpB3&H8&  
end Yg9
joNBh  
setup 5jB*fIz  
apply xacc=1.0 hist=wave  range  y -.1 .1 L|hx arJ  
apply ff UUc8*yU)  
BlA[T%  
hist gp xvel 140,0,0 Cbm\h/PXl  
hist gp xvel 200,98,0 |7k_N|E  
hist gp yvel 200,98,0 `aC){&
AP(  
hist gp xdis 122,41,0 ;30nd=  
hist gp xdis 122,41,40 5PT5#[  
hist gp ydis 122,41,0 XH}'w9VynR  
hist gp ydis 122,41,40 MGJ.,tK1  
hist dytime ntVS:F
 
hist wave k8AW6oO/i  
set dyn damp local 0.15 ;(approx 5% damping) vBcq_sbo  
solve age 4 r^Zg-|gr  
save dam3.sav he(A3{'  
Ztr Cv?  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 `=lc<T^  

声明下y是竖直方向