【求助】大家看看我的流场程序问题出在哪里

pengkun_laser20

我做的焊接模拟，材料是铝合金，电子束焊接，加的是双椭球热源"tuoqiu1"和"tuoqiu2"
我想模拟出大于700度左右区域的熔池流动情况
为什么流场算出来的速度非常乱呢，而且速度很小，感觉粘度差不多啊
下面是我的程序，刚开始算可能比较慢，要稍等会就好了
/PREP7  
/com,***************input material property*************************
/UNITS,SI   
*set,rb,0.001   
et,1,142
et,2,141
fldata7,prot,dens,constant  
fldata7,prot,visc,table
fldata7,prot,spht,table
fldata7,prot,cond,table
fldata8,nomi,dens,2690  
fldata13,vary,cond,t
fldata13,vary,spht,t
fldata13,vary,visc,t
!define material property   
mptemp  
mptemp,1,298,700,800,900,1200,1500  
mpdata,dens,1,1,2690,2690,2690  
mpdata,dens,1,4,2690,2690,2690  
MPTEMP  
MPTEMP,1,298,700,800,900,1200,1500  
MPDATA,C,1,1,963,1063,1163  
mpdata,c,1,4,1263,1363,1463
MPTEMP  
MPTEMP,1,298,700,800,900
mptemp,5,1200,1500  
MPDATA,VISC,1,1,0.6E+10,1.3,1.1,1.0
mpdata,visc,1,5,0.98,0.92  
MPTEMP  
MPTEMP,1,298,700,800,900,1200,1500  
MPDATA,kxx,1,1,167,157,147  
mpdata,kxx,1,4,137,127,117  
MPDATA,kyy,1,1,167,157,147  
mpdata,kyy,1,4,137,127,117  
MPDATA,kzz,1,1,167,157,147  
mpdata,kzz,1,4,137,127,117  
block,0,0.1,0,0.05,-0.008,0
TYPE,2  
MAT,1   
REAL,   
ESYS,0  
SECNUM,
!*  
FLST,5,2,4,ORDE,2   
FITEM,5,1   
FITEM,5,3   
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1,0.005, , , , , , ,1  
!*  
FLST,5,2,4,ORDE,2   
FITEM,5,2   
FITEM,5,4   
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1,0.003, , , , , , ,1  
!*  
MSHAPE,0,2D
MSHKEY,1
!*  
CM,_Y,AREA  
ASEL, , , ,       1
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
!*  
AMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
FLST,5,68,2,ORDE,2  
FITEM,5,1   
FITEM,5,-68
CM,_Y,ELEM  
ESEL, , , ,P51X
CM,_Y1,ELEM
CMSEL,S,_Y  
CMDELE,_Y   
!*  
!*  
EREF,_Y1, , ,1,0,0,1
CMDELE,_Y1  
!*  
vplot   
FLST,5,4,4,ORDE,2   
FITEM,5,9   
FITEM,5,-12
CM,_Y,LINE  
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y   
!*  
LESIZE,_Y1,0.001, , , , , , ,1  
!*  
CM,_Y,VOLU  
VSEL, , , ,       1
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VSWEEP,_Y1  
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
!*  
SAVE
ACLEAR,       1
vplot   
!*  
ETDEL,2
!*  
SAVE
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNCCSYS   
*SET,_FNCNAME,'tuoqiu1'
*SET,_FNCCSYS,0
! /INPUT,tuoqiu1.func,,,1   
*DIM,%_FNCNAME%,TABLE,6,53,1,,,,%_FNCCSYS%  
!   
! Begin of equation: (6*sqrt(3)*0.6*6000*exp(-3*(({X}-0.02*{TIME})**2)/(0.0009*
! *2)-3*({Y}**2)/(0.002**2)-3*({Z}**2)/(0.035**2)))/(0.0009*0.002*0.035*PI*sqrt
! (PI))
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 16, 3, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 6, 0, 0, -1
*SET,%_FNCNAME%(0,3,1),   0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 0.6, 0, 0, -3   
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 6000, 0, 0, -2  
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 0, 0, 0, 0  
*SET,%_FNCNAME%(0,9,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,10,1), 0.0, -4, 0, 1, -1, 2, -2   
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,12,1), 0.0, -2, 0, 1, -4, 3, -1   
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 0.02, 0, 0, 1  
*SET,%_FNCNAME%(0,14,1), 0.0, -4, 0, 1, -1, 3, 1
*SET,%_FNCNAME%(0,15,1), 0.0, -1, 0, 1, 2, 2, -4
*SET,%_FNCNAME%(0,16,1), 0.0, -4, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,17,1), 0.0, -5, 0, 1, -1, 17, -4  
*SET,%_FNCNAME%(0,18,1), 0.0, -1, 0, 1, -2, 3, -5   
*SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 0.0009, 0, 0, 0
*SET,%_FNCNAME%(0,20,1), 0.0, -4, 0, 2, 0, 0, -2
*SET,%_FNCNAME%(0,21,1), 0.0, -5, 0, 1, -2, 17, -4  
*SET,%_FNCNAME%(0,22,1), 0.0, -2, 0, 1, -1, 4, -5   
*SET,%_FNCNAME%(0,23,1), 0.0, -1, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,24,1), 0.0, -4, 0, 1, 3, 17, -1   
*SET,%_FNCNAME%(0,25,1), 0.0, -1, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,26,1), 0.0, -5, 0, 1, -1, 3, -4   
*SET,%_FNCNAME%(0,27,1), 0.0, -1, 0, 0.002, 0, 0, 0
*SET,%_FNCNAME%(0,28,1), 0.0, -4, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,29,1), 0.0, -6, 0, 1, -1, 17, -4  
*SET,%_FNCNAME%(0,30,1), 0.0, -1, 0, 1, -5, 4, -6   
*SET,%_FNCNAME%(0,31,1), 0.0, -4, 0, 1, -2, 2, -1   
*SET,%_FNCNAME%(0,32,1), 0.0, -1, 0, 2, 0, 0, 4
*SET,%_FNCNAME%(0,33,1), 0.0, -2, 0, 1, 4, 17, -1   
*SET,%_FNCNAME%(0,34,1), 0.0, -1, 0, 3, 0, 0, -2
*SET,%_FNCNAME%(0,35,1), 0.0, -5, 0, 1, -1, 3, -2   
*SET,%_FNCNAME%(0,36,1), 0.0, -1, 0, 0.035, 0, 0, 0
*SET,%_FNCNAME%(0,37,1), 0.0, -2, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,38,1), 0.0, -6, 0, 1, -1, 17, -2  
*SET,%_FNCNAME%(0,39,1), 0.0, -1, 0, 1, -5, 4, -6   
*SET,%_FNCNAME%(0,40,1), 0.0, -2, 0, 1, -4, 2, -1   
*SET,%_FNCNAME%(0,41,1), 0.0, -1, 7, 1, -2, 0, 0
*SET,%_FNCNAME%(0,42,1), 0.0, -2, 0, 1, -3, 3, -1   
*SET,%_FNCNAME%(0,43,1), 0.0, -1, 0, 0.0009, 0, 0, 0
*SET,%_FNCNAME%(0,44,1), 0.0, -3, 0, 0.002, 0, 0, -1
*SET,%_FNCNAME%(0,45,1), 0.0, -4, 0, 1, -1, 3, -3   
*SET,%_FNCNAME%(0,46,1), 0.0, -1, 0, 0.035, 0, 0, -4
*SET,%_FNCNAME%(0,47,1), 0.0, -3, 0, 1, -4, 3, -1   
*SET,%_FNCNAME%(0,48,1), 0.0, -1, 0, 3.14159265358979310, 0, 0, -3  
*SET,%_FNCNAME%(0,49,1), 0.0, -4, 0, 1, -3, 3, -1   
*SET,%_FNCNAME%(0,50,1), 0.0, -1, 16, 3.14159265358979310, 0, 0, 0  
*SET,%_FNCNAME%(0,51,1), 0.0, -3, 0, 1, -4, 3, -1   
*SET,%_FNCNAME%(0,52,1), 0.0, -1, 0, 1, -2, 4, -3   
*SET,%_FNCNAME%(0,53,1), 0.0, 99, 0, 1, -1, 0, 0
! End of equation: (6*sqrt(3)*0.6*6000*exp(-3*(({X}-0.02*   
! {TIME})**2)/(0.0009**2)-3*({Y}**2)/(0.002**2)-3*({Z}**2)/(0.035**2)))/(0.0009
! *0.002*0.035*PI*sqrt(PI))
!-->
*DEL,_FNCNAME   
*DEL,_FNCMTID   
*DEL,_FNCCSYS   
*SET,_FNCNAME,'tuoqiu2'
*SET,_FNCCSYS,0
! /INPUT,tuoqiu2.func,,,1   
*DIM,%_FNCNAME%,TABLE,6,53,1,,,,%_FNCCSYS%  
!   
! Begin of equation: (6*sqrt(3)*1.4*6000*exp(-3*(({X}-0.02*{TIME})**2)/(0.0027*
! *2)-3*({Y}**2)/(0.002**2)-3*({Z}**2)/(0.035**2)))/(0.0027*0.002*0.035*PI*sqrt
! (PI))
*SET,%_FNCNAME%(0,0,1), 0.0, -999   
*SET,%_FNCNAME%(2,0,1), 0.0
*SET,%_FNCNAME%(3,0,1), 0.0
*SET,%_FNCNAME%(4,0,1), 0.0
*SET,%_FNCNAME%(5,0,1), 0.0
*SET,%_FNCNAME%(6,0,1), 0.0
*SET,%_FNCNAME%(0,1,1), 1.0, -1, 16, 3, 0, 0, 0
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 6, 0, 0, -1
*SET,%_FNCNAME%(0,3,1),   0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,4,1), 0.0, -1, 0, 1.4, 0, 0, -3   
*SET,%_FNCNAME%(0,5,1), 0.0, -2, 0, 1, -3, 3, -1
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 6000, 0, 0, -2  
*SET,%_FNCNAME%(0,7,1), 0.0, -3, 0, 1, -2, 3, -1
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 0, 0, 0, 0  
*SET,%_FNCNAME%(0,9,1), 0.0, -2, 0, 1, 0, 0, -1
*SET,%_FNCNAME%(0,10,1), 0.0, -4, 0, 1, -1, 2, -2   
*SET,%_FNCNAME%(0,11,1), 0.0, -1, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,12,1), 0.0, -2, 0, 1, -4, 3, -1   
*SET,%_FNCNAME%(0,13,1), 0.0, -1, 0, 0.02, 0, 0, 1  
*SET,%_FNCNAME%(0,14,1), 0.0, -4, 0, 1, -1, 3, 1
*SET,%_FNCNAME%(0,15,1), 0.0, -1, 0, 1, 2, 2, -4
*SET,%_FNCNAME%(0,16,1), 0.0, -4, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,17,1), 0.0, -5, 0, 1, -1, 17, -4  
*SET,%_FNCNAME%(0,18,1), 0.0, -1, 0, 1, -2, 3, -5   
*SET,%_FNCNAME%(0,19,1), 0.0, -2, 0, 0.0027, 0, 0, 0
*SET,%_FNCNAME%(0,20,1), 0.0, -4, 0, 2, 0, 0, -2
*SET,%_FNCNAME%(0,21,1), 0.0, -5, 0, 1, -2, 17, -4  
*SET,%_FNCNAME%(0,22,1), 0.0, -2, 0, 1, -1, 4, -5   
*SET,%_FNCNAME%(0,23,1), 0.0, -1, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,24,1), 0.0, -4, 0, 1, 3, 17, -1   
*SET,%_FNCNAME%(0,25,1), 0.0, -1, 0, 3, 0, 0, -4
*SET,%_FNCNAME%(0,26,1), 0.0, -5, 0, 1, -1, 3, -4   
*SET,%_FNCNAME%(0,27,1), 0.0, -1, 0, 0.002, 0, 0, 0
*SET,%_FNCNAME%(0,28,1), 0.0, -4, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,29,1), 0.0, -6, 0, 1, -1, 17, -4  
*SET,%_FNCNAME%(0,30,1), 0.0, -1, 0, 1, -5, 4, -6   
*SET,%_FNCNAME%(0,31,1), 0.0, -4, 0, 1, -2, 2, -1   
*SET,%_FNCNAME%(0,32,1), 0.0, -1, 0, 2, 0, 0, 4
*SET,%_FNCNAME%(0,33,1), 0.0, -2, 0, 1, 4, 17, -1   
*SET,%_FNCNAME%(0,34,1), 0.0, -1, 0, 3, 0, 0, -2
*SET,%_FNCNAME%(0,35,1), 0.0, -5, 0, 1, -1, 3, -2   
*SET,%_FNCNAME%(0,36,1), 0.0, -1, 0, 0.035, 0, 0, 0
*SET,%_FNCNAME%(0,37,1), 0.0, -2, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,38,1), 0.0, -6, 0, 1, -1, 17, -2  
*SET,%_FNCNAME%(0,39,1), 0.0, -1, 0, 1, -5, 4, -6   
*SET,%_FNCNAME%(0,40,1), 0.0, -2, 0, 1, -4, 2, -1   
*SET,%_FNCNAME%(0,41,1), 0.0, -1, 7, 1, -2, 0, 0
*SET,%_FNCNAME%(0,42,1), 0.0, -2, 0, 1, -3, 3, -1   
*SET,%_FNCNAME%(0,43,1), 0.0, -1, 0, 0.0027, 0, 0, 0
*SET,%_FNCNAME%(0,44,1), 0.0, -3, 0, 0.002, 0, 0, -1
*SET,%_FNCNAME%(0,45,1), 0.0, -4, 0, 1, -1, 3, -3   
*SET,%_FNCNAME%(0,46,1), 0.0, -1, 0, 0.035, 0, 0, -4
*SET,%_FNCNAME%(0,47,1), 0.0, -3, 0, 1, -4, 3, -1   
*SET,%_FNCNAME%(0,48,1), 0.0, -1, 0, 3.14159265358979310, 0, 0, -3  
*SET,%_FNCNAME%(0,49,1), 0.0, -4, 0, 1, -3, 3, -1   
*SET,%_FNCNAME%(0,50,1), 0.0, -1, 16, 3.14159265358979310, 0, 0, 0  
*SET,%_FNCNAME%(0,51,1), 0.0, -3, 0, 1, -4, 3, -1   
*SET,%_FNCNAME%(0,52,1), 0.0, -1, 0, 1, -2, 4, -3   
*SET,%_FNCNAME%(0,53,1), 0.0, 99, 0, 1, -1, 0, 0
! End of equation: (6*sqrt(3)*1.4*6000*exp(-3*(({X}-0.02*   
! {TIME})**2)/(0.0027**2)-3*({Y}**2)/(0.002**2)-3*({Z}**2)/(0.035**2)))/(0.0027
! *0.002*0.035*PI*sqrt(PI))
!-->
FINISH  
/SOL
allsel  
*do,t,0,30,0.005
*SET,xcenter,0.02*t
nsel,s,loc,x,xcenter,0.01   
nsel,a,loc,y,0,0.05
nsel,a,loc,z,-0.008,0   
esln,s,0,ACTIVE
bfe,all,hgen,,%tuoqiu1%
allsel  
nsel,s,loc,x,0,xcenter  
nsel,a,loc,y,0,0.05
nsel,a,loc,z,-0.008,0   
esln,s,0,ACTIVE
bfe,all,hgen,,%tuoqiu2%
*enddo  
allsel  
FLDATA1,SOLU,TRAN,1    !
FLDATA1,SOLU,FLOW,1
FLDATA1,SOLU,TEMP,1    !
FLDATA1,SOLU,TURB,0
FLDATA1,SOLU,COMP,0
FLDATA1,SOLU,VOF,0  
FLDATA1,SOLU,SFTS,0
FLDATA1,SOLU,IVSH,0
FLDATA1,SOLU,SWRL,0
FLDATA1,SOLU,SPEC,0
FLDATA1,SOLU,ALE,0  
!*  
FLDATA12,PROP,DENS,0   !
FLDATA13,VARY,DENS,0   !
FLDATA12,PROP,VISC,16  !粘度迭代次数为16
FLDATA13,VARY,VISC,1   !
FLDATA12,PROP,COND,16  !导热系数迭代次数16  
FLDATA13,VARY,COND,1   !
FLDATA12,PROP,SPHT,16  !比热系数迭代次数16  
FLDATA13,VARY,SPHT,1   !
!*  
FLDATA7,PROT,DENS,TABLE
FLDATA8,NOMI,DENS,-1,   
FLDATA9,COF1,DENS,0
FLDATA10,COF2,DENS,0
FLDATA11,COF3,DENS,0
FLDATA7,PROT,VISC,TABLE
FLDATA8,NOMI,VISC,-1
FLDATA9,COF1,VISC,0
FLDATA10,COF2,VISC,0
FLDATA11,COF3,VISC,0
FLDATA12,PROP,IVIS  
FLDATA7,PROT,COND,TABLE
FLDATA8,NOMI,COND,-1
FLDATA9,COF1,COND,0
FLDATA10,COF2,COND,0
FLDATA11,COF3,COND,0
FLDATA7,PROT,SPHT,TABLE
FLDATA8,NOMI,SPHT,-1
FLDATA9,COF1,SPHT,0
FLDATA10,COF2,SPHT,0
FLDATA11,COF3,SPHT,0
!*  
ACEL,0,0,-9.8,  
asel,s,,,1  
asel,a,,,2  
asel,a,,,3  
asel,a,,,4  
asel,a,,,5  
asel,a,,,6  
DA,all,PRES,1013500,1   
allsel  
FLDATA15,PRES,REFE,101350000000,
FLDATA16,BULK,BETA,1e+015,  
FLDATA17,GAMM,COMP,1.4,
FLDATA14,TEMP,NOMI,293,
FLDATA14,TEMP,TTOT,293,
FLDATA14,TEMP,BULK,293,
TOFFSET,273,
fldata31,capp,temp,T           
fldata31,capp,tmax,3500
FLDATA4,TIME,STEP,0.005
FLDATA4,TIME,ISTEP,0,   
FLDATA4,TIME,NUMB,20,   
FLDATA4,TIME,TEND,0.5   
FLDATA4,TIME,GLOB,10,  
FLDATA4,TIME,PRES,1e-006,   
FLDATA4,TIME,TEMP,1e-006,   
FLDATA4,TIME,APPE,0.002   
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
allsel  
solve

pengkun_laser20

大牛赶紧出现啊！！！啊啊啊啊

crruose

非真空
FLDATA15,PRES,REFE,101350000000,   为什么这么大   麻烦LZ解释下
关键在于没有加载表面张力

ansys平安

太有难度了。

pengkun_laser20

电子束是真空环境啊
我知道FLDATA15的缺省值是101350
但程序是从师兄那里的改过来的 这个值也就没变了  
我变过 好像作用不大
流场一定要加表面张力么
我看有些流场不加力也可以算的啊
表面张力要加的话怎么加呢
万分感谢！
3# crruose