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

代码如下：set dyn=on large 3> /K0N|$  
ini xvel=0 yvel=0 zvel=0 ktDC/8
  
ini xdis 0 ydis 0 zdis 0 d GP*O  
free x y z range y -.1 .1 Wu)>U  
def setup R *F l8
 
omega = 4.0 * pi dL|+d:v  
alfa = 2.2 0a"igq9t  
beta = 3.75e-1 !n^OM?.
4  
gamma = 8.0 u4Em%:Xj  
end <3,<\ub  
def wave ]({~,8s  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) ] }f9JNf$  
end >vo=]cw  
setup 4Vq%N  
apply xacc=1.0 hist=wave  range  y -.1 .1 ,^icPQSwc  
apply ff M
Qin"\  
{nU=%w"\  
hist gp xvel 140,0,0 {}:ToIp  
hist gp xvel 200,98,0 v9T_
&  
hist gp yvel 200,98,0 r H~" 4  
hist gp xdis 122,41,0 I@\OaUGr+  
hist gp xdis 122,41,40 }^B6yWUN  
hist gp ydis 122,41,0 Ytg
j|@jsp  
hist gp ydis 122,41,40 soCi[j$lH  
hist dytime F:mq'<Q  
hist wave BL&LeSa  
set dyn damp local 0.15 ;(approx 5% damping) (rg;IXAq%  
solve age 4 )?wJF<[_#  
save dam3.sav ;2Q~0a|  
A=Ss6-Je  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 %c[V  

声明下y是竖直方向