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

代码如下：set dyn=on large #isBE}sT{  
ini xvel=0 yvel=0 zvel=0 he/UvMu  
ini xdis 0 ydis 0 zdis 0 #GYCU!  
free x y z range y -.1 .1 PT|W{RlNl  
def setup a"{tqNc  
omega = 4.0 * pi $ #C$V>  
alfa = 2.2 m#'2 3  
beta = 3.75e-1 EV@xUq!x.  
gamma = 8.0 V$wf;v0d(  
end I hvL2zB  
def wave J0}OmNTzD  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) rBP!RSl1  
end *fq=["O  
setup Ywf.,V  
apply xacc=1.0 hist=wave  range  y -.1 .1 $b#"Rv  
apply ff l|fOi A*K  
FxC@KZG  
hist gp xvel 140,0,0 qPpC)6-Q  
hist gp xvel 200,98,0 5vL]Y)l  
hist gp yvel 200,98,0 6|05-x|  
hist gp xdis 122,41,0 ti'a^(  
hist gp xdis 122,41,40 "YGs<)S  
hist gp ydis 122,41,0 Xv'M\T}6C+  
hist gp ydis 122,41,40 ztG_::QtG]  
hist dytime ?EeHeN_  
hist wave M[0@3"}}  
set dyn damp local 0.15 ;(approx 5% damping) EM*YN=So  
solve age 4 )#r]x1[Kn  
save dam3.sav m?_S&/+*  
PCs+` WP!M  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 k[N46=u  

声明下y是竖直方向