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

代码如下：set dyn=on large iMt3h8  
ini xvel=0 yvel=0 zvel=0 [zBi*%5O  
ini xdis 0 ydis 0 zdis 0 a_+?#m  
free x y z range y -.1 .1 ]+46r!r|  
def setup y+T[="W  
omega = 4.0 * pi ~uH_y
-  
alfa = 2.2 S :8  
beta = 3.75e-1 Pw| h`[h  
gamma = 8.0 =/_uk{  
end r8@]|`j  
def wave g9q}D-  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) y] c1x=x  
end |ZE^'e*k  
setup Db<#gH  
apply xacc=1.0 hist=wave  range  y -.1 .1 @J&korU  
apply ff WB?HY?[r  
:IU7dpwDl  
hist gp xvel 140,0,0 nX 9]dz
 
hist gp xvel 200,98,0 S\h5 D2G;  
hist gp yvel 200,98,0 HO['o{>BL  
hist gp xdis 122,41,0 JLn
v O  
hist gp xdis 122,41,40 ka!v(j{E  
hist gp ydis 122,41,0 A$r$g\
5+  
hist gp ydis 122,41,40 D/f4kkd  
hist dytime );':aXj  
hist wave ;<N:!$p  
set dyn damp local 0.15 ;(approx 5% damping) =$Mf:F@  
solve age 4 ?"r=08  
save dam3.sav Hc=QSP  
C<AW)|r_  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 C&>*~  

声明下y是竖直方向