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

代码如下：set dyn=on large ]1MZ:]k  
ini xvel=0 yvel=0 zvel=0 ((M>To_l  
ini xdis 0 ydis 0 zdis 0 2s}G6'xE]P  
free x y z range y -.1 .1 <`NtTG  
def setup 4~2 9,  
omega = 4.0 * pi [kMWsiZ  
alfa = 2.2 @h(Z;  
beta = 3.75e-1 P1|3%#c  
gamma = 8.0 <@c@`K  
end i?,\>LTG  
def wave R0K{wY58  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) 0n\AUgVPF  
end +`!>lo{X  
setup .vd*~U"  
apply xacc=1.0 hist=wave  range  y -.1 .1 
`#l1  
apply ff 0qm CIcg  
AqHGBH0  
hist gp xvel 140,0,0 ,[3}t%
Da  
hist gp xvel 200,98,0 qRq4PQ@  
hist gp yvel 200,98,0 m ;wj|@cF  
hist gp xdis 122,41,0 <v?2p{U%  
hist gp xdis 122,41,40 ^2i$AM1t  
hist gp ydis 122,41,0 ux)Wh.5  
hist gp ydis 122,41,40 B;nIKZ  
hist dytime /v.<h*hxWy  
hist wave ]`|;ZQiD  
set dyn damp local 0.15 ;(approx 5% damping) =&~ K;=:  
solve age 4 hsJGly5H  
save dam3.sav 8Nx fYA  
Pfd%[C/vdm  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 ft~QVe!  

声明下y是竖直方向