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

代码如下：set dyn=on large L{-LX=G^  
ini xvel=0 yvel=0 zvel=0 u\yVR$pQ  
ini xdis 0 ydis 0 zdis 0 #%0Bx3uM  
free x y z range y -.1 .1 Fh$slow4!  
def setup GLMm(  
omega = 4.0 * pi \3f&7wU  
alfa = 2.2 .B2]xfo"`  
beta = 3.75e-1 ]`g@UtD9`  
gamma = 8.0
=FP0\cQ.  
end mHc>"^R  
def wave 4GdX/6C.  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) FS6`6M.K  
end ').}Nz  
setup hVUh0XeO  
apply xacc=1.0 hist=wave  range  y -.1 .1 tBbOY}.VD  
apply ff 4<dcB@v  
Qn@Pd*DR  
hist gp xvel 140,0,0 (t+;O;  
hist gp xvel 200,98,0 H,unpZ(  
hist gp yvel 200,98,0 ZBT1Y.qA  
hist gp xdis 122,41,0 I#F!N6;  
hist gp xdis 122,41,40 K<`osdp=&  
hist gp ydis 122,41,0 ED>prE0  
hist gp ydis 122,41,40 `F YjQe"p  
hist dytime tJViA`@x  
hist wave W !w,f;  
set dyn damp local 0.15 ;(approx 5% damping) i:]*P  
solve age 4 XRx+Dddt;  
save dam3.sav b;sVls  
T;TA7{B  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 :KJ pk:<  

声明下y是竖直方向