【讨论】求救！温度分布差别为何如此大！
这是我做的二维轴对称焊补温度场的命令流。焊补时，在中心上表面挖r10mm,深为h10mm的凹坑。对此区域进行补焊。焊接电压28v,电流200I,热效率0.75，有效加热半径rh=7.5mm。焊接加热时间：15s,采用高斯内热源加热。热生成率在厚度（y)方向一致。在径向(x)高斯分布。得到的温度场在上表面中心0~10mm相差太大。理想温度场(15s时)：中心2500度，r10mm处1500度左右。恳请各位大虾run 一下。帮忙分析！谢谢！
/prep7
et,1,plane55,
keyopt,1,1,1!根据固体表面温度计算换热系数
keyopt,1,3,1!轴对称
/units,si
mptemp,1,0,400,600,800,1000,2000
mpdata,kxx,1,1,18.9,28.35,31.08,32.76,31.5,25.2
mpdata,c,1,1,580,756,840,882,756,126
mpdata,enth,1,1,0,2.0174e9,3.2272e9,4.5272e9,5.7639e9,9.4574e9
mpdata,hf,1,1,9.66,59.43,84.32,109.2,180.6,537.6
mpdata,dens,1,1,7550,7550,7550,7550,7550,7550
!输入材料热物理参数
Qmax=7130000000 !最大热生成率
rectng,0,0.1,0,0.03
wpave,0.01
wprota,,,90
asbw,all
wpave,0.02
asbw,all
wpave,0.04
asbw,all
wpave,0.06
asbw,all
wpave,0,0.02
wprota,,90
asbw,all
wpave,0,0.01
asbw,all
WPSTYLE,,,,,,,,0
/pnum,area,1
aglue,all
!建模
! X向划分网格
lsel,s,,,9
lsel,a,,,24
lsel,a,,,16
lsel,a,,,8
lesize,all,0.001
allsel
lsel,s,,,12
lsel,a,,,28
lsel,a,,,35
lsel,a,,,11
lesize,all,0.0015
allsel
lsel,s,,,15
lsel,a,,,31
lsel,a,,,37
lsel,a,,,14
lesize,all,0.002
allsel
lsel,s,,,18
lsel,a,,,33
lsel,a,,,38
lsel,a,,,17
lesize,all,0.003
allsel
lsel,s,,,3
lsel,a,,,19
lsel,a,,,7
lsel,a,,,1
lesize,all,0.005
! Y向划分网格
lsel,s,,,25
lsel,a,,,26
lsel,a,,,29
lsel,a,,,21
lsel,a,,,20
lesize,all,0.001
lsel,s,,,39,44,
lesize,all,0.0012
lsel,s,,,10
lsel,a,,,13
lsel,a,,,34
lsel,a,,,36
lsel,a,,,4
lsel,a,,,2
lesize,all,0.0015
allsel
type,1
mat,1
mshkey,1
amesh,all
!划分网格
/solu
antype,trans
trnopt,full
nropt,auto, ,
solcontrol,on
autots,on
kbc,1 !step
outres,all,all,
tunif,20
tref,0
lsel,s,loc,y,0
lsel,a,loc,y,0.03
lsel,a,loc,x,0.1
nsll,s,1,1,
sf,all,conv,-1,20 !换热系数随材料温度变化
allsel
eplot
!设定初始温度求解
time,0.001
deltim,0.001,0.001,0.001
timint,off !time integration = off
lswrite,1
*dim,b,,11
*do,i,0,8
b(i+1)=Qmax*exp(-3*i*i*0.5625/56.25) !高斯热源
*enddo
*dim,hg,table,11,1,1,x,,
*do,i,1,11
hg(i,1,1)=b(i)
*enddo
*do,i,0,10,1
hg(i+1,0,1)=i*0.00075
*enddo
asel,s,,,9
esla,s,
nsle,s,1,1
BF,all,HGEN,%hg%
allsel
time,15
deltim,0.1,0.1,0.1
kbc,1 !step加载
timint,1 !瞬态分析
tintp,0.005, , ,1,0.5,0.2
outres,all,all
lnsrch,on
lswrite,2
bfdele,all,HGEN
time,30
deltim,0.1,0.1,0.1
autots,on
lswrite,3
time,100
deltim,1,0.5,5
autots,on
lswrite,4
lssolve,1,4,1
/post1
plnsol,temp
asel,s,,,9
esla,s
nsle,s,1,1
set,,,,,15
plnsol,temp
PATH,x-temp,2,100,100
PPATH,1,,0,0.03
PPATH,2,,0.01,0.03
PDEF,temp15,TEMP,
/grid,on
PLPATH,temp15
!X轴温度分布
页:
[1]