【推荐】TIG焊接温度场和流场的模拟

转自仿真  
  
https://www.simwe.com/cgi-bin/ut/topic_show.cgi?id=27450&h=1#148159  
----------------------------  
这是我做的模拟的命令流，问题一定不少，希望高手们不吝赐教！  
/FILNAME,flotran,0  
/TITLE,flotran transient analysis  
/PREP7  
smrt,off  
ET,1,fluid142  
KEYOPT,1,1,2  
BLOCK,-0.025,0.025,-0.01,0,-0.25,0.25,  
block,-0.004,0.004,-0.01,0,-0.25,0.25,  
vsel,s,volu,,1,2  
vovlap,all  
vglue,all  
NUMCMP,LINE  
NUMCMP,AREA  
NUMCMP,VOLU  
/REPLOT  
AESIZE,ALL,0.002,  
LESIZE,ALL,0.002, , , ,1, , ,1,  
ESIZE,0.002,0,  
MSHKEY,1  
MSHAPE,0,3d  
vmesh,2,3  
AESIZE,ALL,0.002,  
LESIZE,ALL,0.002, , , ,1, , ,1,  
ESIZE,0.002,0,  
MSHKEY,1  
MSHAPE,0,3d  
vmesh,1  
/VIEW, 1 ,1,1,1  
/ANG, 1  
/REP,FAST  
/psf,hflu,,2,0  
/PNUM,LINE,1  
asel,s,,,12,16,4!最上层节点的速度vy设为0，其他的速度为0.其他表面节点仍设为三个方向的速度为0  
asel,a,,,6  
nsla,s,1  
d,all,vy,0  
allsel,all  
asel,s,ext,,all  
asel,u,,,12,16,4  
asel,u,,,6  
nsla,s,1  
D,all,,,,,,VX,vy,vz  
asel,s,,,12,16,4  
asel,a,,,1,2,1  
nsla,s,1  
sf,all,conv,50e-6,293  
asel,s,,,5  
asel,a,,,11,15,4  
nsla,s,,,1  
sf,all,conv,25e-6,293  
allsel,all  
!*  
!flow environment  
FLDATA15,PRES,REFE,101350,  
FLDATA16,BULK,BETA,1e+015,  
FLDATA17,GAMM,COMP,1.4,  
FLDATA14,TEMP,NOMI,293,  
FLDATA14,TEMP,TTOT,293,  
FLDATA14,TEMP,BULK,293,  
TOFFSET,0,  
nsel,s,,,all  
ic,all,sp01,1.0  
ic,all,sp02,0.0  
finish  
!*  
/solu  
allsel,all  
!*solution options  
FLDATA1,SOLU,TRAN,1 !瞬态分析  
FLDATA1,SOLU,FLOW,1  
fldata1,solu,temp,1  
/COM,,Transient Analysis,1  
FLDATA4,TIME,STEP,2!自定义时间步时间间隔2秒  
FLDATA4,TIME,ISTEP,0,  
FLDATA4,TIME,NUMB,10,  
FLDATA4,TIME,TEND,8!设置结束时间8秒  
FLDATA4,TIME,GLOB,10, !设置每个时间步多少次运算  
FLDATA4,TIME,PRES,1e-006,  
FLDATA4,TIME,TEMP,1e-006,  
FLDATA4,TIME,APPE,2 !设置记录时间间隔  
FLDATA4A,STEP,SUMF,10,  
FLDATA4,TIME,SUMF,1.0e6  
FLDATA4,TIME,BC,0  
fldata2,iter,exec,30  
fldata2,iter,over,2  
fldata5,outp,sumf,30  
!*fluid properties  
fldata13,vary,dens,1  
FLDATA7,PROT,DENS,LIQUID  
FLDATA8,NOMI,DENS,7800,  
FLDATA9,COF1,DENS,300,  
FLDATA10,COF2,DENS,0.0001,  
FLDATA11,COF3,DENS,0,  
FLDATA7,PROT,visc,1  
FLDATA7,PROT,VISC,LIQUID  
FLDATA8,NOMI,VISC,1,  
FLDATA9,COF1,VISC,300,  
FLDATA10,COF2,VISC,1644,  
FLDATA11,COF3,VISC,0,  
FLDATA7,PROT,cond,constant  
flda,nomi,cond,70  
FLDATA7,PROT,spht,constant  
fldata8,nomi,spht,470  
!*flow environment  
acel,0,-9.81,0  
fldata18,meth,pres,1  
fldata19,tdma,pres,100  
FLDA,OUTP,SP01,T  
FLDA,OUTP,SP02,T  
csys,0  
nsel,s,loc,y,0  
csys,2  
nsel,r,loc,x,0,0.006  
csys,0  
sf,all,hflux,50000000!%ga%!热源是固定的，此处热流应该加高斯分布的，不过我加了以后得出的结果和实际温度场差别很大，就用这个代替了，如果是移动热源，ga为200*24*0.6*100/{PI}/{PI}/  
(0.003^2)*exp(-1e6*100*({X}^2+({Z}-8.3e-3* {TIME})^2))，不知哪个参数取得不对？  
allsel,all  
solve  
!下面是一位师姐教我的用一个循环加载表面张力，对此还是一窍不通，还请高手指点！  
*do,i,1,8  
time=i*0.01  
/post1  
set,last  
!此处应该获取节点温度，不知该如何做？  
/solu  
csys,0  
nsel,s,loc,y,0  
csys,2  
nsel,r,loc,x,0.001,0.006!加表面张力我用了一个近似  
csys,0  
F,all,FX,%tx%!tx为0.00001*{X}/sqrt({X}^2+{Z}^2)  
F,all,FZ,%tz%!tz为0.00001*{Z}/sqrt({X}^2+{Z}^2)  
allsel,all  
solve  
*enddo