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

代码如下：set dyn=on large :o!bz>T  
ini xvel=0 yvel=0 zvel=0 IU f1N+-z  
ini xdis 0 ydis 0 zdis 0 ?]JTrv"zp  
free x y z range y -.1 .1 v%mAU3M  
def setup 1e Wl:S}  
omega = 4.0 * pi v37TDY3;  
alfa = 2.2 942(a  
beta = 3.75e-1 xwSi}.  
gamma = 8.0 _RzwE$+9  
end H;#3S<  
def wave SONv]));  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) + ?z=,')  
end ~;m3i3D  
setup Pfy2PpA  
apply xacc=1.0 hist=wave  range  y -.1 .1 D1Yc_  
apply ff fSe$
w#*I  
<"Ox)XG3]W  
hist gp xvel 140,0,0 H9d!-9I  
hist gp xvel 200,98,0 Z/y&;N4   
hist gp yvel 200,98,0 Q
ktj  
hist gp xdis 122,41,0 _;o)MTw|'  
hist gp xdis 122,41,40 80i-)a\n  
hist gp ydis 122,41,0 ;<aT|4  
hist gp ydis 122,41,40 ~Ibq,9i  
hist dytime )qgcz<p?W  
hist wave Gp0B^^H$  
set dyn damp local 0.15 ;(approx 5% damping) fjc8@S5x9j  
solve age 4 qs96($  
save dam3.sav IK:F~I  
=qQH,{]c6  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 i_ha^mq3  

声明下y是竖直方向