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

代码如下：set dyn=on large !e(ZEV g  
ini xvel=0 yvel=0 zvel=0 wA>bLPTw  
ini xdis 0 ydis 0 zdis 0 bcy( ?(  
free x y z range y -.1 .1 "K
$ y(}C  
def setup Co^a$K  
omega = 4.0 * pi m8ydX6~max  
alfa = 2.2 hR3Pa'/i  
beta = 3.75e-1 "K*+8IO2  
gamma = 8.0 C
%+>uzVIw  
end wfc[B;K\  
def wave > [|SF%  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) ky5gU[  
end !kg)84C[  
setup !B\[Q$  
apply xacc=1.0 hist=wave  range  y -.1 .1 ~*Ir\w
E  
apply ff !WReThq  
$/uNV1]o  
hist gp xvel 140,0,0 q{t"=@lX01  
hist gp xvel 200,98,0 #pz{,  
hist gp yvel 200,98,0
bUi@4S  
hist gp xdis 122,41,0 |~3$L\X  
hist gp xdis 122,41,40 rO`nS<G  
hist gp ydis 122,41,0 sw@*N  
hist gp ydis 122,41,40 tHo0q<.oX  
hist dytime )iG+pP@.@  
hist wave 4D\+_Ic3  
set dyn damp local 0.15 ;(approx 5% damping)  p3r1lUw  
solve age 4 lZu
pn?  
save dam3.sav sP NAG  
k^PqB+P!  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 Ev\kq>2O  

声明下y是竖直方向