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

代码如下：set dyn=on large Wzqb
>.   
ini xvel=0 yvel=0 zvel=0 P-+M,>vNy[  
ini xdis 0 ydis 0 zdis 0 F@$RV_M  
free x y z range y -.1 .1 ,MM>cOQ  
def setup 8pPC 9ew\=  
omega = 4.0 * pi Aw4?y[{H  
alfa = 2.2 ^.#X<8hr  
beta = 3.75e-1 gr>o E#7  
gamma = 8.0 >&;>PZBPCO  
end z
]c,}Q  
def wave TH>?Gi)"  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) Q)Iv_N/  
end o8'Mks  
setup 2]D$|M?$~  
apply xacc=1.0 hist=wave  range  y -.1 .1 hDljY!P>p  
apply ff qB F!b0lr  
9$+^"ilk  
hist gp xvel 140,0,0 R6!cK[e]4  
hist gp xvel 200,98,0 bEBBwv  
hist gp yvel 200,98,0 \- =^]]b=  
hist gp xdis 122,41,0 yQZ/,KX  
hist gp xdis 122,41,40 sm;E2BR$ `  
hist gp ydis 122,41,0 S\9t4Ki_'  
hist gp ydis 122,41,40 u`Z0{d  
hist dytime @0z0m;8  
hist wave zr.+'  
set dyn damp local 0.15 ;(approx 5% damping) y~ _za(k  
solve age 4 "J pTE \/  
save dam3.sav q#99iiG1  
{?*<B=c  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 qByNHo7Tb  

声明下y是竖直方向