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

代码如下：set dyn=on large '3,JL!  
ini xvel=0 yvel=0 zvel=0 [+2^n7R  
ini xdis 0 ydis 0 zdis 0 u(Q(UuI  
free x y z range y -.1 .1 T%b^|="@  
def setup Wa<NId  
omega = 4.0 * pi (?=(eo<N  
alfa = 2.2 Y+PxV*"a  
beta = 3.75e-1 .Sth  
gamma = 8.0 p~y 4q4  
end %JU23c*  
def wave yOm6HA``hT  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) WxI]Fcb<  
end /IR5[67  
setup IQ`aDo-V  
apply xacc=1.0 hist=wave  range  y -.1 .1 ~wV98u-N  
apply ff Lgw@y!Llij  
X>YOo~yS5  
hist gp xvel 140,0,0 m=b+V#4i(  
hist gp xvel 200,98,0 wH5O>4LO  
hist gp yvel 200,98,0 8IcQpn#  
hist gp xdis 122,41,0 JQv ZTwSI  
hist gp xdis 122,41,40 5_T>HHR6  
hist gp ydis 122,41,0 '<jp.sZQ  
hist gp ydis 122,41,40 &?6~v  
hist dytime Gf$>!zXr  
hist wave A.*nDl`H  
set dyn damp local 0.15 ;(approx 5% damping) {twf7.eY  
solve age 4 Hq
y>!1!  
save dam3.sav {+59YO  
o?x|y   
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 Nr7.BDA  

声明下y是竖直方向