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

代码如下：set dyn=on large G lz0`z  
ini xvel=0 yvel=0 zvel=0 =~OH.=9\  
ini xdis 0 ydis 0 zdis 0 f{b$Y3  
free x y z range y -.1 .1 (BB&ZUdyv  
def setup QbF!V%+a's  
omega = 4.0 * pi SMMV$;O{9  
alfa = 2.2 'u
\my
 
beta = 3.75e-1 *1}UK9X;  
gamma = 8.0 Rmq8lU  
end X&
B2&e;  
def wave ,?OV39h  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) k/"^W.B aj  
end Ft7{P.g  
setup &8l"Dl  
apply xacc=1.0 hist=wave  range  y -.1 .1 j^t#>tZS  
apply ff $\W|{u`  
?,_$;g  
hist gp xvel 140,0,0 VSK!Pc.G}
 
hist gp xvel 200,98,0 v<*ga7'S  
hist gp yvel 200,98,0 WBo|0(#  
hist gp xdis 122,41,0 )FNvtLZ  
hist gp xdis 122,41,40 $.a4Og2  
hist gp ydis 122,41,0 W[5a'}OV  
hist gp ydis 122,41,40 %n^jho5  
hist dytime h";0i:  
hist wave i2a""zac  
set dyn damp local 0.15 ;(approx 5% damping) %R18  
solve age 4 E3pnu.;U:_  
save dam3.sav 9@nDXZPY&  
GI&XL'K&  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 sH//*y
  

声明下y是竖直方向