ansys焊接热分析直接提示不收敛。大佬看下哪里有错。到底哪里有错啊?
/prep7
/units,si
mptemp,1,20,250,500,750,1000,1500
mptemp,7,1700,2500
mpdata,kxx,1,1,50,47,40,27,30,35
mpdata,kxx,1,7,140,142
mpdata,ex,1,1,2.05e+11,1.87e+11,1.50e+11,7.00e+10,2.00e+10,1.90e+10
mpdata,ex,1,7,1.80e+10,1.20e+10
mpdata,dens,1,1,7820,7700,7610,7550,7490,7350
mpdata,dens,1,7,7300,7090
mpdata,c,1,1,460,480,530,675,670,660
mpdata,c,1,7,780,820
mpdata,prxy,1,1,0.28,0.29,0.31,0.35,0.4,0.45
mpdata,prxy,1,7,0.48,0.5
mpdata,alpx,1,1,1.1e-05,1.22e-05,1.39e-05,1.48e-05,1.34e-05,1.33e-05,1.32e-05,1.31e-05
mpdata,alpx,1,7,1.32e-05,1.31e-05
mptemp,,,,,,,,
mptemp,1,20,1533,1590,1670
mpdata,enth,1,1,0,7.5e+009,9.6e+009,1.1e+010
rb=0.002 !单元尺寸
length=0.12 !长方体几何参数
width=0.06
height=0.006
et,1,solid70
block,0,length,0,width,0,height
Block,0,length,0,3*Rb,0,height
Block,0,length,0,6*Rb,0,height
Block,0,length,0,12*Rb,0,height
VOVLAP,ALL
Lesize,4,rb/2 !mapped划分单元尺寸
lesize,5,rb/2
lesize,19,rb/2
lesize,14,rb/2
lesize,50,rb/2
lesize,10,rb/2
lesize,52,rb/2
lesize,23,rb/2
lesize,51,rb/2
lesize,9,rb/2
lesize,49,rb/2
lesize,24,rb/2
lesize,31,rb
lesize,26,rb
lesize,43,4*rb
lesize,38,4*rb
lesize,2,4*rb
lesize,7,4*rb
lesize,53,4*rb
lesize,54,4*rb
lesize,55,4*rb
lesize,56,4*rb
lesize,11,0.003
lesize,12,0.003
lesize,47,0.003
lesize,48,0.003
lesize,35,0.003
lesize,36,0.003
lesize,64,0.004
lesize,61,0.004
lesize,62,0.004
lesize,63,0.004
lesize,60,0.002
lesize,57,0.002
lesize,58,0.002
lesize,59,0.002
type,1
mshkey,1
mshape,0
vmesh,5
vmesh,6
type,1
mshkey,0
mshape,1
vmesh,7
vmesh,8
/solu !求解设置
antype,trans
timint,on,therm
tref,25
nropt,full,
kbc,1
autots,on
*do,i,1,60,1 !高斯移动热源加载
b=0.2*i !时间步
time,b
deltim,0.2,0.01,0.2
xcenter=i*0.002 !加载中心点
*DEL,_FNCNAME
*DEL,_FNCMTID
*DEL,_FNCCSYS
*SET,_FNCNAME,'flux2'
*SET,_FNCCSYS,0
! /INPUT,C:\Users\Administrator\Desktop\gaosi.func,,,1
*DIM,%_FNCNAME%,TABLE,6,19,1,,,,%_FNCCSYS%
!
! Begin of equation: 89524655/exp(222222*({x}-0.01*{time})^2+({y})^2+
! ({z}-0.006)^2)
*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, 0, 0.01, 0, 0, 1
*SET,%_FNCNAME%(0,2,1), 0.0, -2, 0, 1, -1, 3, 1
*SET,%_FNCNAME%(0,3,1), 0, -1, 0, 1, 2, 2, -2
*SET,%_FNCNAME%(0,4,1), 0.0, -2, 0, 2, 0, 0, -1
*SET,%_FNCNAME%(0,5,1), 0.0, -3, 0, 1, -1, 17, -2
*SET,%_FNCNAME%(0,6,1), 0.0, -1, 0, 222222, 0, 0, -3
*SET,%_FNCNAME%(0,7,1), 0.0, -2, 0, 1, -1, 3, -3
*SET,%_FNCNAME%(0,8,1), 0.0, -1, 0, 2, 0, 0, 3
*SET,%_FNCNAME%(0,9,1), 0.0, -3, 0, 1, 3, 17, -1
*SET,%_FNCNAME%(0,10,1), 0.0, -1, 0, 1, -2, 1, -3
*SET,%_FNCNAME%(0,11,1), 0.0, -2, 0, 0.006, 0, 0, 4
*SET,%_FNCNAME%(0,12,1), 0.0, -3, 0, 1, 4, 2, -2
*SET,%_FNCNAME%(0,13,1), 0.0, -2, 0, 2, 0, 0, -3
*SET,%_FNCNAME%(0,14,1), 0.0, -4, 0, 1, -3, 17, -2
*SET,%_FNCNAME%(0,15,1), 0.0, -2, 0, 1, -1, 1, -4
*SET,%_FNCNAME%(0,16,1), 0.0, -1, 7, 1, -2, 0, 0
*SET,%_FNCNAME%(0,17,1), 0.0, -2, 0, 89524655, 0, 0, -1
*SET,%_FNCNAME%(0,18,1), 0.0, -3, 0, 1, -2, 4, -1
*SET,%_FNCNAME%(0,19,1), 0.0, 99, 0, 1, -3, 0, 0
! End of equation: 89524655/exp(222222*({x}-0.01*{time})^2+({y})^2+({z}-0.006)^
! 2)
allsel,all !选取加载中心点半径6mm内的点施加热流密度
nsel,s,loc,z,0.006
nsel,r,loc,x,xcenter-0.006,xcenter+0.006
nsel,r,loc,y,0,0.006
local,11,1,xcenter,0,0.006
csys,11
nsel,r,loc,x,0,0.006
csys,0
sf,all,hflux,%flux2%allsel,all !高度方向4个面施加对流载荷
nsel,s,loc,x,0
sf,all,conv,10,20
allsel,all
nsel,s,loc,x,0.12
sf,all,conv,10,20
allsel,all
nsel,s,loc,y,0.06
sf,all,conv,10,20
allsel,all
nsel,s,loc,y,0
sf,all,conv,10,20
allsel,all
solve
*enddo
你的命令:
sf,all,hflux,%flux2%allsel,all !高度方向4个面施加对流载荷
是在节点上施加热流载荷,不是在面上施加对流载荷,不知是否你需要的。检查一下载荷值是否太大了。否则,温度结果太大,难以收敛。
而且该命令最后的 allsel,all 也应该另起一行,放在这里不知会造成什么后果。
另外,你的求解设置过于简单。可以在求解前增加以下命令试试:
nsubst, 10, 100,2
autots,on
kbc,0
outres,all,all
/gst,1
ncnv,2