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

代码如下：set dyn=on large -'*B
%yy  
ini xvel=0 yvel=0 zvel=0 }c:s+P+/  
ini xdis 0 ydis 0 zdis 0 [yW0U:m  
free x y z range y -.1 .1 Dc}-wnga  
def setup a>ZV'~zTf  
omega = 4.0 * pi LW{7|g  
alfa = 2.2 18zv]v %  
beta = 3.75e-1 z^ rf;  
gamma = 8.0 ~]O~a}]g(  
end 69p>?zn  
def wave 5,Mc`IIK1  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) gC> A*~J;  
end
j[F
\f>  
setup 6H53FMqr  
apply xacc=1.0 hist=wave  range  y -.1 .1 J)"2^?!&B  
apply ff T9}G:6  
9NBFG~)|l[  
hist gp xvel 140,0,0 f)Z$,&  
hist gp xvel 200,98,0 qm{(.b^  
hist gp yvel 200,98,0 B#OnooJI  
hist gp xdis 122,41,0 to|O]h2*U2  
hist gp xdis 122,41,40 >h(n8wTP  
hist gp ydis 122,41,0 [?Cv^t${+  
hist gp ydis 122,41,40 :uSo2d  
hist dytime [h {zT)[  
hist wave H"=%|/1M0  
set dyn damp local 0.15 ;(approx 5% damping) :54ik,l  
solve age 4 U_K"JOZ  
save dam3.sav  zN:VT&  
(1^AzE%U+Z  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 5N@k9x  

声明下y是竖直方向