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

代码如下：set dyn=on large QWEK;kUa@  
ini xvel=0 yvel=0 zvel=0 |zE7W  
ini xdis 0 ydis 0 zdis 0 fkM4u<R^  
free x y z range y -.1 .1 |QdS;  
def setup
Pc<ZfO #  
omega = 4.0 * pi WRCi!  
alfa = 2.2 P+a&R<Dj4  
beta = 3.75e-1 ^M"g5+q  
gamma = 8.0 }$ der  
end RP$A"<goP  
def wave 7=9jXNk Y  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) Q@R8qc=*  
end |HLh?AcX  
setup (%1*<6ka  
apply xacc=1.0 hist=wave  range  y -.1 .1 C{-pVuhK+  
apply ff DS yE   
3@PVUJ0B|  
hist gp xvel 140,0,0 \b->AXe8  
hist gp xvel 200,98,0 BlCKJp{m$  
hist gp yvel 200,98,0 Y/gCtSF  
hist gp xdis 122,41,0 QPnc "!  
hist gp xdis 122,41,40 04:Dbt~=?p  
hist gp ydis 122,41,0 A6VkVJZx  
hist gp ydis 122,41,40 4Ki'r&L\  
hist dytime >e%Po,Fg$  
hist wave rxA<\h,A  
set dyn damp local 0.15 ;(approx 5% damping) <V{BRRx  
solve age 4 P^UcpU,  
save dam3.sav X+iULr.^`~  
7w|s8B  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 t<tBOesQ  

声明下y是竖直方向