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

代码如下：set dyn=on large c0!Te'?  
ini xvel=0 yvel=0 zvel=0 4'+d"Ok  
ini xdis 0 ydis 0 zdis 0 &q," !:L]  
free x y z range y -.1 .1 QE7+rBa  
def setup 9O),/SH;:  
omega = 4.0 * pi 0=N4O!X9
 
alfa = 2.2 g>6:CG"  
beta = 3.75e-1 p&k0Rx0Q3  
gamma = 8.0 iyA*JCD  
end 6obQ9L c  
def wave 4/*]`  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) 7j@^+rkr3f  
end Ep^B
,;~  
setup w"fCI13  
apply xacc=1.0 hist=wave  range  y -.1 .1 M,b<B_$  
apply ff lSZ"y Q+  
9>A-$a4R>  
hist gp xvel 140,0,0 + $k07mb\  
hist gp xvel 200,98,0 $x;tSJ)m~  
hist gp yvel 200,98,0  O]e6i%?  
hist gp xdis 122,41,0 Nf=C?`L  
hist gp xdis 122,41,40 t!6\7Vm/  
hist gp ydis 122,41,0  )x$!K[=  
hist gp ydis 122,41,40 gzl%5`DBw  
hist dytime tpv?`(DDU  
hist wave ^z[_U}N\}  
set dyn damp local 0.15 ;(approx 5% damping) oS[W*\7'!  
solve age 4 q1N4X7<_  
save dam3.sav [TRGIGtq  
P\D[n-&  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 Bv;I0i:_  

声明下y是竖直方向