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

代码如下：set dyn=on large Z@umbyM  
ini xvel=0 yvel=0 zvel=0 :S2MS{>Mo  
ini xdis 0 ydis 0 zdis 0 D(GHkS*0q  
free x y z range y -.1 .1 4<HJD&@V  
def setup >FhBl\oIi  
omega = 4.0 * pi $ {"St&(  
alfa = 2.2 c9={~  
beta = 3.75e-1 7dW&|U  
gamma = 8.0 Q&;qFv5-l  
end ,~w)@.  
def wave Q:=/d$*xd  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) 06O  
end t)+dW~g  
setup 0\;a:E.c  
apply xacc=1.0 hist=wave  range  y -.1 .1 &(7
Io?  
apply ff rrq-so1u}  
f>PU# D@B  
hist gp xvel 140,0,0 'D{abm0  
hist gp xvel 200,98,0 7 {<lH%Tn  
hist gp yvel 200,98,0 *mtv[  
hist gp xdis 122,41,0 ]d(}b>gR~(  
hist gp xdis 122,41,40 r4zS,J;,  
hist gp ydis 122,41,0 4h(Hy&1C  
hist gp ydis 122,41,40 GT0'bge  
hist dytime hQeZI+  
hist wave (q7mzZY  
set dyn damp local 0.15 ;(approx 5% damping) ?uv%E*TU  
solve age 4 9)X<}*(qo  
save dam3.sav zvg&o)/[  
4\RuJx  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 {S~$\4vC!  

声明下y是竖直方向