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

代码如下：set dyn=on large Ymu=G3-  
ini xvel=0 yvel=0 zvel=0 v
(O.GhJ@  
ini xdis 0 ydis 0 zdis 0 QFYy$T+W  
free x y z range y -.1 .1 G;;~xfE'  
def setup 68*{Lo?U  
omega = 4.0 * pi oX~CTunP  
alfa = 2.2 :Eyv==  
beta = 3.75e-1 86g+c  
gamma = 8.0 qu%s 7+  
end h(-&.Sm")H  
def wave di5_5_$`o  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) 7fqYSMHR  
end a5M>1&j/eC  
setup Ul Iw&U  
apply xacc=1.0 hist=wave  range  y -.1 .1 ~}*;Ko\  
apply ff _1'Pb/1  
as4NvZ@+r  
hist gp xvel 140,0,0 }t'^Au`X  
hist gp xvel 200,98,0 ~aG-^BAS  
hist gp yvel 200,98,0 U:9vjY  
hist gp xdis 122,41,0 UJ[a&b  
hist gp xdis 122,41,40 ?s]`G'=>V`  
hist gp ydis 122,41,0 |_*$+  
hist gp ydis 122,41,40 $5r,Q{;$  
hist dytime A2d2V**Z  
hist wave C.-a:oQ[  
set dyn damp local 0.15 ;(approx 5% damping) s_LSsyqo  
solve age 4 @\?f77Of6  
save dam3.sav B(LV22#  
@
UJmbD{  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 ;>9pJ72r  

声明下y是竖直方向