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

代码如下：set dyn=on large z=U!D `]v  
ini xvel=0 yvel=0 zvel=0 Z-?9F`}  
ini xdis 0 ydis 0 zdis 0 m~K]|]iqQ  
free x y z range y -.1 .1 .,,73"  
def setup H6eGLg={  
omega = 4.0 * pi %Z? o]  
alfa = 2.2 #Grm-W9E  
beta = 3.75e-1 2P}RZvUd  
gamma = 8.0 y()( 8L  
end #wyS?FP-  
def wave 
uI[*uAR  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) A0ToX) |C  
end )em.KbsPPF  
setup !Z
ZAI_N  
apply xacc=1.0 hist=wave  range  y -.1 .1 '9gI=/29D  
apply ff m @%|Q;  
9lx
T5Wg  
hist gp xvel 140,0,0 wMoAvA_oS  
hist gp xvel 200,98,0 U{Z>y?V/  
hist gp yvel 200,98,0 @!da1jN  
hist gp xdis 122,41,0 ^J_hkw~gO  
hist gp xdis 122,41,40 dw|0K+-PH  
hist gp ydis 122,41,0 ik*_,51Zj  
hist gp ydis 122,41,40 ]lY9[~ v  
hist dytime >>r:L3<!  
hist wave loJ0PY'}=  
set dyn damp local 0.15 ;(approx 5% damping) *Y ZLQT  
solve age 4 K&h6#[^\d  
save dam3.sav P.:T zk6  
ihVQ,Cth  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 6>I.*Qt \l  

声明下y是竖直方向