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

代码如下：set dyn=on large gQJLqs"F  
ini xvel=0 yvel=0 zvel=0 E&5S[n9{3  
ini xdis 0 ydis 0 zdis 0 oK$Krrs0&  
free x y z range y -.1 .1 VL'wrgk  
def setup eZ-fy,E  
omega = 4.0 * pi PsD)]V9%:  
alfa = 2.2 e,lLHg  
beta = 3.75e-1 w.#z>4#3-  
gamma = 8.0 #9~,d<H  
end TQ0ZBhd  
def wave Mc!LC .8  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) 7AWq3i{  
end c27(en(  
setup $Sa7N%D  
apply xacc=1.0 hist=wave  range  y -.1 .1 P ?nk>  
apply ff pS}IU{#;  
UB
k:B  
hist gp xvel 140,0,0 .w'b%M  
hist gp xvel 200,98,0 "fUNrhCx  
hist gp yvel 200,98,0 _*b`;{3  
hist gp xdis 122,41,0 axq~56"7E  
hist gp xdis 122,41,40 )VqPaKZl  
hist gp ydis 122,41,0 \fuz`fK:  
hist gp ydis 122,41,40 ]yL+lv  
hist dytime `"<hO 'WU  
hist wave He}?\C Bo  
set dyn damp local 0.15 ;(approx 5% damping) ~*<`PDO?  
solve age 4 `l/nAKg?W  
save dam3.sav 9)o@d`*  
sLXM$SMBh  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 B692Mn  

声明下y是竖直方向