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

代码如下：set dyn=on large F+m[&MKL  
ini xvel=0 yvel=0 zvel=0 8!&nKy<Y  
ini xdis 0 ydis 0 zdis 0 -ey)J +?t  
free x y z range y -.1 .1 ?OE#q$g  
def setup v9=}S\=Cd  
omega = 4.0 * pi um7o!yg,  
alfa = 2.2 s.VA!@F5  
beta = 3.75e-1 {Bh("wg$Lk  
gamma = 8.0 /Y,r@D  
end Ea-bC:>  
def wave F|Q H
 
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) (ew} gJ  
end 3V?817&6z  
setup  A^ViDP  
apply xacc=1.0 hist=wave  range  y -.1 .1 w]nX?S8  
apply ff !siWEzw  
Z&Ue|Z4Qt  
hist gp xvel 140,0,0 Z0$] tS  
hist gp xvel 200,98,0 +c--&tBo  
hist gp yvel 200,98,0 Z0-ytODII  
hist gp xdis 122,41,0 i\)3l%AK]T  
hist gp xdis 122,41,40 \@K
~L4>  
hist gp ydis 122,41,0 Ql8bt77eI-  
hist gp ydis 122,41,40 gw^'{b  
hist dytime b._m8z ~  
hist wave V>Fesm"aq  
set dyn damp local 0.15 ;(approx 5% damping) fxd0e;NAAh  
solve age 4 WWNu:,  
save dam3.sav B8H75sz  
kx:j
I^  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 YGp)Oy}:  

声明下y是竖直方向