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

代码如下：set dyn=on large .A|@?p[  
ini xvel=0 yvel=0 zvel=0 (w zQ2Dk  
ini xdis 0 ydis 0 zdis 0 #rg6,.I)<  
free x y z range y -.1 .1 c1gQ cqF  
def setup DW3G
 
omega = 4.0 * pi JsS-n'gF'  
alfa = 2.2 f)<6  
beta = 3.75e-1 R^e'}+Z  
gamma = 8.0 H6gSO(U  
end &,)&%Sg[  
def wave [PbOfxxgA  
wave = sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime) &6k3*dq  
end o;RI*I  
setup .eC1qWZJpd  
apply xacc=1.0 hist=wave  range  y -.1 .1 UL9n-M=  
apply ff %SUQ9\SEs  
o,wUc"CE  
hist gp xvel 140,0,0 ;9'OOz|+1  
hist gp xvel 200,98,0 'E.w=7z&  
hist gp yvel 200,98,0 @KUWxFak  
hist gp xdis 122,41,0 =WJNWt>  
hist gp xdis 122,41,40 EBmt9S
 
hist gp ydis 122,41,0 nT)vNWT=  
hist gp ydis 122,41,40 bQ5\ ]5M  
hist dytime Ht&YC<X  
hist wave &>}5jC.I  
set dyn damp local 0.15 ;(approx 5% damping) wS3'?PRX  
solve age 4 a09<!0Rp  
save dam3.sav H%lVl8oQ  
W(/h Vt  
怎末会出现 unknown fatel exception 啊？困惑好长时间不能解决 希望高手帮忙 R/a*LSe@&  

声明下y是竖直方向