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

代码如下：set dyn=on large X:lPWz!7{  
ini xvel=0 yvel=0 zvel=0 L]d@D0.Z  
ini xdis 0 ydis 0 zdis 0 Z~c'h  
free x y z range y -.1 .1 /\rq$W_  
def setup -kbm$~P  
omega = 4.0 * pi wkD:i2E7  
alfa = 2.2 }4SSo)Uv/  
beta = 3.75e-1 (0W}e(D8  
gamma = 8.0 t1jlxK  
end 1>w^ q`P  
def wave ht)nx,e=  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) = O1;vc}AA  
end JtpY][}"~3  
setup Da[C'm=  
apply xacc=1.0 hist=wave  range  y -.1 .1
L\NZDkd  
apply ff N@6OQ:,[F  
A Vm{#^p[(  
hist gp xvel 140,0,0 Z=@)  
hist gp xvel 200,98,0 N?;o_^C  
hist gp yvel 200,98,0 u!o]Co>  
hist gp xdis 122,41,0 `mjx4Lb  
hist gp xdis 122,41,40 NRisr  
hist gp ydis 122,41,0 |xZcT4  
hist gp ydis 122,41,40 AwUc{h l<  
hist dytime | ",[C3Jg  
hist wave \oX8/-0f  
set dyn damp local 0.15 ;(approx 5% damping) OZD!#YI  
solve age 4 9T2A)a]0  
save dam3.sav R9h>I3F=c  
zpq
Gh  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 c^O#O  

声明下y是竖直方向