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

代码如下：set dyn=on large Fy-+? ~  
ini xvel=0 yvel=0 zvel=0 VNx|nP&  
ini xdis 0 ydis 0 zdis 0 jRXpEiM  
free x y z range y -.1 .1 jyF0asb  
def setup [07N<<  
omega = 4.0 * pi 84[T!cDk  
alfa = 2.2 ;u LD_1%  
beta = 3.75e-1 t5h]]TOz  
gamma = 8.0 ]saf<?fzr  
end |Ia3bVW  
def wave mLM$dk3  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) _%Ay\4H^\  
end (80#{4kl  
setup R4,j
  
apply xacc=1.0 hist=wave  range  y -.1 .1 C)> ])'S  
apply ff h'wOslyFa  
gBRhO^Sz  
hist gp xvel 140,0,0 =Haqr*PDx  
hist gp xvel 200,98,0 jnFCtCB  
hist gp yvel 200,98,0 3=xb%Upw  
hist gp xdis 122,41,0 B\&;eZY'G  
hist gp xdis 122,41,40 IC}?oXs5G  
hist gp ydis 122,41,0 't=\YFQ*v  
hist gp ydis 122,41,40 ;i`X&[y;  
hist dytime W(9fCDO;  
hist wave S5ka;g  
set dyn damp local 0.15 ;(approx 5% damping) nGA'\+zjL  
solve age 4 Xz5 aTJ&  
save dam3.sav c@:L7#8  
}9fH`C/m  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙  "m3:HS  

声明下y是竖直方向