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

代码如下：set dyn=on large 3#>%_@<  
ini xvel=0 yvel=0 zvel=0 Qc PU{#6  
ini xdis 0 ydis 0 zdis 0 h%4UeL &F  
free x y z range y -.1 .1 &%qDi_UD  
def setup R(cg`8  
omega = 4.0 * pi Tm7LaM  
alfa = 2.2 .c__T{<)[  
beta = 3.75e-1 +*Uv+oC|  
gamma = 8.0 5ddfdIp  
end KU+\fwYpnk  
def wave Ld/6{w4ir  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) 7-~)/7L  
end Y"yrc0'&T  
setup *">CEQ[MT  
apply xacc=1.0 hist=wave  range  y -.1 .1 IA]wO%c  
apply ff 9d(#/n  
Hr!$mf)h  
hist gp xvel 140,0,0 D<<q5gG  
hist gp xvel 200,98,0 R osU~OK  
hist gp yvel 200,98,0 M_2>b:#A*  
hist gp xdis 122,41,0 83{P7PBQ;]  
hist gp xdis 122,41,40 `Lavjmfr2V  
hist gp ydis 122,41,0 -!li,&,A1  
hist gp ydis 122,41,40 LEOa=(mN\  
hist dytime {d0-.  
hist wave l+hOD{F4pS  
set dyn damp local 0.15 ;(approx 5% damping) 7y)Ar 8!D  
solve age 4 \+l_H4\`K  
save dam3.sav _VmXs&4  
iDhC_F|  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 bQwG"N  

声明下y是竖直方向