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

代码如下：set dyn=on large n,Z B-"dW  
ini xvel=0 yvel=0 zvel=0 (:p&[HNuN  
ini xdis 0 ydis 0 zdis 0 '$cU\DTN6  
free x y z range y -.1 .1 /y\KLa  
def setup CDz-IQi  
omega = 4.0 * pi I&}Md73  
alfa = 2.2 Q\>Kd N{  
beta = 3.75e-1 |VC/(A  
gamma = 8.0 .-:R mYGR  
end mST/u>'  
def wave fYU-pdWPT  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) O*<,lq 0K  
end T]5Jsr
T  
setup ye9-%~sjX  
apply xacc=1.0 hist=wave  range  y -.1 .1 "ct_EPr`  
apply ff CENA!WWQ  
XOM@Pi#z  
hist gp xvel 140,0,0 TT(dCHft  
hist gp xvel 200,98,0 hE-u9i  
hist gp yvel 200,98,0 8^H <dR  
hist gp xdis 122,41,0 ;tR,w   
hist gp xdis 122,41,40 pGy]t  
hist gp ydis 122,41,0 FuYV}C  
hist gp ydis 122,41,40 XG5mfKMt+  
hist dytime ,UdTUw~F  
hist wave e/?>6'6 5  
set dyn damp local 0.15 ;(approx 5% damping) #nL&x3  
solve age 4 y`~[R7E  
save dam3.sav @Y#{[@Hp%  
F
afOd9>AO  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 s-Mzl?o  

声明下y是竖直方向