温度分布为什么会这样?请高手指点!
这是我的一个例子,是一个壁和底都是四层的不同材料组成的圆柱,里层是恒定1600度,外面和空气对流传热,可是结果不对,跟实践不符合,不知道问题在哪里?我在各层体之间用vglue实现热量在公共边界的传递 ,但是在底壁外壳的温度竟然跟里面的温度一样,这是我的命令流和温度贴图,请高手指教,谢谢!/filename,ladle_mo_11
/title,ladle_mo_11
th1=170
th2=75
th3=32
h1=230
h2=130
h3=32
/prep7
!包壁尺寸
r=1900
h=3093
th4=36
h_1=1000
!包底尺寸
h4=75
!透气砖1尺寸(定位)
tou_posx=1200
tou_posy=0
tou_block=360
!透气砖2尺寸(定位)
tou_posx1=-1200
tou_posy1=0
!冲击块
chong_wide=800
chong_length=1600
distx=0
disty=800
!水口尺寸(定位)
water_posx=-318
water_posy=-1379
water_block=510
wpcsys,-1,0
!指定单元
et,1,plane55
et,2,solid70
!定义材料
!包底工作层 高铝质浇注料
mptemp,1,20,400,800,1200
mpdata,kxx,1,1,1.31,1.246,1.18,1.11 !!
mp,c,1,0.8e6 !!
mp,alpx,1,5.8e-6
mp,dens,1,2.8e-6 !!
mp,ex,1,5.7e3
mp,nuxy,1,0.21
!包底第一永久层! 高铝质砖
mpdata,kxx,2,1,1.31,1.246,1.18,1.11 !!
mp,c,2,0.8e6 !!
mp,alpx,2,5.8e-6 !!
mp,dens,2,2.8e-6 !!
mp,ex,2,5.7e3
mp,nuxy,2,0.21
!包底第二永久层!
mpdata,kxx,3,1,1.31,1.246,1.18,1.11 !!
mp,c,3,0.8e6 !!
mp,alpx,3,5.8e-6 !!
mp,dens,3,2.8e-6 !!
mp,ex,3,5.7e3 !!
mp,nuxy,3,0.21
!包壁工作层! 铝镁碳
mpdata,kxx,4,1,5.05,3.2,2.646,2.66 !!
mp,c,4,1.07e6 !!
mp,alpx,4,8.5e-6 !!
mp,dens,4,2.95e-6 !!
mp,ex,4,6.3e3
mp,nuxy,4,0.21
!包壁第一永久层! 高铝
mpdata,kxx,5,1,1.31,1.246,1.18,1.11 !!
mp,c,5,0.8e6 !!
mp,alpx,5,5.8e-6
mp,dens,5,2.8e-6 !!
mp,ex,5,5.7e3
mp,nuxy,5,0.21
!包壁第二永久层!
mpdata,kxx,6,1,1.31,1.246,1.18,1.11 !!
mp,c,6,0.8e6 !!
mp,alpx,6,5.8e-6
mp,dens,6,2.8e-6 !!
mp,ex,6,5.7e3
mp,nuxy,6,0.21
!透气砖!
mp,kxx,7,3.5 !!
mp,c,7,1.134e6 !!
mp,alpx,7,5e-6 !!
mp,dens,7,3e-6 !!
mp,nuxy,7,0.22 !!
mp,ex,7,6.2e3 !!
!水口砖!
mp,kxx,8,3.5 !!
mp,c,8,1.134e6 !!
mp,alpx,8,5e-6 !!
mp,dens,8,3e-6 !!
mp,ex,8,6.2e3 !!
mp,nuxy,8,0.22 !!
!冲击块! 刚玉
mpdata,kxx,9,1,1.8,2.46,3.08,3.6
mp,c,9,1.13e6 !!
mp,alpx,9,8.8e-6 !!
mp,dens,9,3.1e-6 !!
mp,ex,9,5.7e3 !!有30e3
mp,nuxy,9,0.21
!钢壳!
mp,kxx,10,50 !!
mp,c,10,0.470e6 !!
mp,alpx,10,13e-6 !!
mp,dens,10,7.8e-6 !!
mp,ex,10,175e3 !!
mp,nuxy,10,0.3 !!
!水口,透气周围的铝镁碳砖
mpdata,kxx,11,1,5.05,3.2,2.646,2.66 !!
mp,c,11,1.07e6 !!铝的含量65% 镁10%
mp,alpx,11,8.5e-6 !!
mp,dens,11,2.95e-6 !!
mp,ex,11,6.3e3
mp,nuxy,11,0.21
!上渣线层 镁碳砖
mpdata,kxx,12,1,5.2,4.3,3.3,2.4!!导热率
mp,c,12,1.08e6 !!比热
mp,alpx,12,15e-6 !!热膨胀系数
mp,dens,12,2.9e-6 !!密度
mp,ex,12,8e3 !!弹性模量 ,有120e3
mp,nuxy,12,0.21 !!泊松比
/nerr,0
!绘制包底
!第一层
pcirc,r
cm,di1,area !将选择的实体生成一个元件
cmsel,none
rectang,tou_posx-tou_block/2-65,tou_posx+tou_block/2+65,tou_posy-tou_block/2-65,tou_posy+tou_block/2+65
rectang,tou_posx1-tou_block/2-65,tou_posx1+tou_block/2+65,tou_posy1-tou_block/2-65,tou_posy1+tou_block/2+65
rectang,water_posx-water_block/2-65,water_posx+water_block/2+65,water_posy-water_block/2-65,water_posy+water_block/2+65
cm,hole1,area
cmsel,a,di1
asba,di1,hole1,,dele,dele
cm,di_sub_hole,area!减去透气,水口(包括周围的铝镁碳内衬砖)后的底
!冲击块
cmsel,none
cyl4,0,0,0,0,600,90
cm,temp_area1,area
cmsel,none
rectang,475,600,0,600
cm,temp_area2,area
cmsel,a,temp_area1
asba,temp_area1,temp_area2,,dele,dele
cm,chong_temp1,area
arsym,x,chong_temp1 !对面chong_temp1进行x轴镜像
aadd,all !面相加
cm,chong_temp2,area
arsym,y,chong_temp2
aadd,all
cm,chongji1_area,area
agen,,chongji1_area,,,distx,disty,,,,1 !平移chong_temp3
cmsel,a,di_sub_hole
asba,di_sub_hole,chongji1_area,,dele,keep
cmsel,u,chongji1_area
cm,di_sub_hole1,area
!水口砖和透气砖
cmsel,none
rectang,tou_posx-tou_block/2,tou_posx+tou_block/2,tou_posy-tou_block/2,tou_posy+tou_block/2
cm,touqi11_area,area
cmsel,none
rectang,tou_posx1-tou_block/2,tou_posx1+tou_block/2,tou_posy1-tou_block/2,tou_posy1+tou_block/2
cm,touqi12_area,area
cmsel,none
rectang, water_posx-water_block/2,water_posx+water_block/2,water_posy-water_block/2,water_posy+water_block/2
cm,water1_area,area
!水口砖和透气砖周围加了一圈铝镁碳材质
cmsel,none
rectang,tou_posx-tou_block/2-65,tou_posx+tou_block/2+65,tou_posy-tou_block/2-65,tou_posy+tou_block/2+65
cm,touqi11_wai,area
cmsel,a,touqi11_area
asba,touqi11_wai,touqi11_area,,delete,keep
cmsel,u,touqi11_area
cm,touqi11_wei,area
cmsel,none
rectang,tou_posx1-tou_block/2-65,tou_posx1+tou_block/2+65,tou_posy1-tou_block/2-65,tou_posy1+tou_block/2+65
cm,touqi12_wai,area
cmsel,a,touqi12_area
asba,touqi12_wai,touqi12_area,,delete,keep
cmsel,u,touqi12_area
cm,touqi12_wei,area
cmsel,none
rectang,water_posx-water_block/2-65,water_posx+water_block/2+65,water_posy-water_block/2-65,water_posy+water_block/2+65
cm,water1_wai,area
cmsel,a,water1_area
asba,water1_wai,water1_area,,delete,keep
cmsel,u,water1_area
cm,water1_wei,area
!复制到各层
!底
cmsel,,di_sub_hole1!减去透气,水口,冲击后的底
agen,2,di_sub_hole1,,,,,-h1
cmsel,u,di_sub_hole1
cm,di_sub_hole2,area
cmsel,,di_sub_hole1
agen,2,di_sub_hole1,,,,,-h1-h2
cmsel,u,di_sub_hole1
cm,di_sub_hole3,area
cmsel,,di_sub_hole1
agen,2,di_sub_hole1,,,,,-h1-h2-h3
cmsel,u,di_sub_hole1
cm,di_sub_hole4,area
!水口砖和透气砖,冲击块
cmsel,,touqi11_area
agen,2,touqi11_area,,,,,-h1
cmsel,u,touqi11_area
cm,touqi21_area,area
cmsel,,touqi11_area
agen,2,touqi11_area,,,,,-h1-h2
cmsel,u,touqi11_area
cm,touqi31_area,area
cmsel,,touqi11_area
agen,2,touqi11_area,,,,,-h1-h2-h3
cmsel,u,touqi11_area
cm,touqi41_area,area
cmsel,,touqi11_wei
agen,2,touqi11_wei,,,,,-h1
cmsel,u,touqi11_wei
cm,touqi21_wei,area
cmsel,,touqi11_wei
agen,2,touqi11_wei,,,,,-h1-h2
cmsel,u,touqi11_wei
cm,touqi31_wei,area
cmsel,,touqi11_wei
agen,2,touqi11_wei,,,,,-h1-h2-h3
cmsel,u,touqi11_wei
cm,touqi41_wei,area
cmsel,,touqi12_area
agen,2,touqi12_area,,,,,-h1
cmsel,u,touqi12_area
cm,touqi22_area,area
cmsel,,touqi12_area
agen,2,touqi12_area,,,,,-h1-h2
cmsel,u,touqi12_area
cm,touqi32_area,area
cmsel,,touqi12_area
agen,2,touqi12_area,,,,,-h1-h2-h3
cmsel,u,touqi12_area
cm,touqi42_area,area
cmsel,,touqi12_wei
agen,2,touqi12_wei,,,,,-h1
cmsel,u,touqi12_wei
cm,touqi22_wei,area
cmsel,,touqi12_wei
agen,2,touqi12_wei,,,,,-h1-h2
cmsel,u,touqi12_wei
cm,touqi32_wei,area
cmsel,,touqi12_wei
agen,2,touqi12_wei,,,,,-h1-h2-h3
cmsel,u,touqi12_wei
cm,touqi42_wei,area
cmsel,,water1_area
agen,2,water1_area,,,,,-h1
cmsel,u,water1_area
cm,water2_area,area
cmsel,,water1_area
agen,2,water1_area,,,,,-h1-h2
cmsel,u,water1_area
cm,water3_area,area
cmsel,,water1_area
agen,2,water1_area,,,,,-h1-h2-h3
cmsel,u,water1_area
cm,water4_area,area
cmsel,,water1_wei
agen,2,water1_wei,,,,,-h1
cmsel,u,water1_wei
cm,water2_wei,area
cmsel,,water1_wei
agen,2,water1_wei,,,,,-h1-h2
cmsel,u,water1_wei
cm,water3_wei,area
cmsel,,water1_wei
agen,2,water1_wei,,,,,-h1-h2-h3
cmsel,u,water1_wei
cm,water4_wei,area
cmsel,,chongji1_area
agen,2,chongji1_area,,,,,-h1
cmsel,u,chongji1_area
cm,chongji2_area,area
cmsel,,chongji1_area
agen,2,chongji1_area,,,,,-h1-h2
cmsel,u,chongji1_area
cm,chongji3_area,area
cmsel,,chongji1_area
agen,2,chongji1_area,,,,,-h1-h2-h3
cmsel,u,chongji1_area
cm,chongji4_area,area
allsel
/pnum,area,1!显示面!绘制成体
!绘制成体
!1
!底
/VIEW,1,1,1,1
/ANG,1
/REP,FAST
cmsel,,di_sub_hole1
vext,di_sub_hole1,,,,,-h1 !通过给定偏移量由面生成第一层底
csys,4 !将坐标系改为工作平面坐标系
wpstyle,,,,,,,,1 !显示工作平面坐标系
wpoff,0,0,-h1
cyl4,0,0,r,,r+th1,,h1
cyl4,0,0,r+th1,,r+th1+th2,,h1
cyl4,0,0,r+th1+th2,,r+th1+th2+th3,,h1
cyl4,0,0,r+th1+th2+th3,,r+th1+th2+th3+th4,,h1
!2
!底
cmsel,,di_sub_hole2
vext,di_sub_hole2,,,,,-h2 !通过给定偏移量由面生成第二层底
wpoff,0,0,-h2
cyl4,0,0,r,,r+th1,,h2
cyl4,0,0,r+th1,,r+th1+th2,,h2
cyl4,0,0,r+th1+th2,,r+th1+th2+th3,,h2
cyl4,0,0,r+th1+th2+th3,,r+th1+th2+th3+th4,,h2
!3
!底
cmsel,,di_sub_hole3
vext,di_sub_hole3,,,,,-h3 !通过给定偏移量由面生成第三层底
wpoff,0,0,-h3
cyl4,0,0,r,,r+th1,,h3
cyl4,0,0,r+th1,,r+th1+th2,,h3
cyl4,0,0,r+th1+th2,,r+th1+th2+th3,,h3
cyl4,0,0,r+th1+th2+th3,,r+th1+th2+th3+th4,,h3
!4
!底
cmsel,,di_sub_hole4
vext,di_sub_hole4,,,,,-h4 !通过给定偏移量由面生成第四层底
wpoff,0,0,-h4
cyl4,0,0,r,,r+th1,,h4
cyl4,0,0,r+th1,,r+th1+th2,,h4
cyl4,0,0,r+th1+th2,,r+th1+th2+th3,,h4
cyl4,0,0,r+th1+th2+th3,,r+th1+th2+th3+th4,,h4
!壁1
wpoff,0,0,h1+h2+h3+h4
cyl4,0,0,r,,r+th1,,h
cyl4,0,0,r+th1,,r+th1+th2,,h
cyl4,0,0,r+th1+th2,,r+th1+th2+th3,,h
cyl4,0,0,r+th1+th2+th3,,r+th1+th2+th3+th4,,h
!壁2
wpoff,0,0,h
cyl4,0,0,r,,r+th1,,h_1
cyl4,0,0,r+th1,,r+th1+th2,,h_1
cyl4,0,0,r+th1+th2,,r+th1+th2+th3,,h_1
cyl4,0,0,r+th1+th2+th3,,r+th1+th2+th3+th4,,h_1
!第一层透气,水口,冲击
cmsel,,touqi11_area
vext,touqi11_area,,,,,-h1
cmsel,,touqi12_area
vext,touqi12_area,,,,,-h1
cmsel,,water1_area
vext,water1_area,,,,,-h1
cmsel,,chongji1_area
vext,chongji1_area,,,,,-h1
!第二层透气,水口,冲击
cmsel,,touqi21_area
vext,touqi21_area,,,,,-h2
cmsel,,touqi22_area
vext,touqi22_area,,,,,-h2
cmsel,,water2_area
vext,water2_area,,,,,-h2
cmsel,,chongji2_area
vext,chongji2_area,,,,,-h2
!第三层透气,水口,冲击
cmsel,,touqi31_area
vext,touqi31_area,,,,,-h3
cmsel,,touqi32_area
vext,touqi32_area,,,,,-h3
cmsel,,water3_area
vext,water3_area,,,,,-h3
cmsel,,chongji3_area
vext,chongji3_area,,,,,-h3
!第四层透气,水口,冲击
cmsel,,touqi41_area
vext,touqi41_area,,,,,-h4
cmsel,,touqi42_area
vext,touqi42_area,,,,,-h4
cmsel,,water4_area
vext,water4_area,,,,,-h4
cmsel,,chongji4_area
vext,chongji4_area,,,,,-h4
cmsel,none
!第一层透气,水口,冲击周围
cmsel,,touqi11_wei
vext,touqi11_wei,,,,,-h1
cmsel,,touqi12_wei
vext,touqi12_wei,,,,,-h1
cmsel,,water1_wei
vext,water1_wei,,,,,-h1
!第二层透气,水口,冲击周围
cmsel,,touqi21_wei
vext,touqi21_wei,,,,,-h2
cmsel,,touqi22_wei
vext,touqi22_wei,,,,,-h2
cmsel,,water2_wei
vext,water2_wei,,,,,-h2
!第三层透气,水口,冲击周围
cmsel,,touqi31_wei
vext,touqi31_wei,,,,,-h3
cmsel,,touqi32_wei
vext,touqi32_wei,,,,,-h3
cmsel,,water3_wei
vext,water3_wei,,,,,-h3
!第四层透气,水口,冲击周围
cmsel,,touqi41_wei
vext,touqi41_wei,,,,,-h4
cmsel,,touqi42_wei
vext,touqi42_wei,,,,,-h4
cmsel,,water4_wei
vext,water4_wei,,,,,-h4
!mat,1
vsel,s,,,1,3,1,1
vatt,1,,2,0
vsel,s,,,5,7,1,1
vatt,1,,2,0
vsel,s,,,9,11,1,1
vatt,1,,2,0
vsel,s,,,13,15,1,1
vatt,1,,2,0
vsel,s,,,22,23,1,1
vatt,1,,2,0
vsel,s,,,26,27,1,1
vatt,1,,2,0
!mat,10
vsel,s,,,4
vatt,10,,2,0
vsel,s,,,8
vatt,10,,2,0
vsel,s,,,12
vatt,10,,2,0
vsel,s,,,16,20,1,1
vatt,10,,2,0
vsel,s,,,24
vatt,10,,2,0
vsel,s,,,28
vatt,10,,2,0
vsel,s,,,56
vatt,10,,2,0
!mat,4
vsel,s,,,21
vatt,4,,2,0
vsel,s,,,45,55,1,1
vatt,4,,2,0
!mat,12
vsel,s,,,25
vatt,12,,2,0
!mat,7
vsel,s,,,29,31,1,1
vatt,7,,2,0
vsel,s,,,33,35,1,1
vatt,7,,2,0
vsel,s,,,37,39,1,1
vatt,7,,2,0
vsel,s,,,41,43,1,1
vatt,7,,2,0
!mat,9
vsel,s,,,32
vatt,9,,2,0
vsel,s,,,36
vatt,9,,2,0
vsel,s,,,40
vatt,9,,2,0
vsel,s,,,44
vatt,9,,2,0
allsel
wpcsys,-1,0 !将工作平面与笛卡儿坐标对齐
vglue,all!粘接,相邻元素具有相同的边界条件,作用是传递温度
esize,240,0,
mshkey,0
mshape,1,3d
!vmesh,all!自由网格,为四面体网格
vsweep,all!六面体网格
numcmp,all
!求解
/solu
antype,static
csys,1
!指定边界条件
asel,none
asel,s,loc,x,r
asel,r,loc,z,0,h+h_1
da,all,temp,1600,1
asel,none
asel,s,loc,x,0,r
asel,r,loc,z,0
da,all,temp,1600,1
asel,none
asel,s,loc,z,-h1-h2-h3-h4,h+h_1
asel,r,loc,x,r+th1+th2+th3+th4
asel,a,loc,z,-h1-h2-h3-h4
sfa,all,,conv,10e-3,30
csys,0
allsel
eqslv,sparse,1e-5
TIME,5
AUTOTS,-1
NSUBST,1,1,1,1
KBC,0
solve
!save
!后处理
/post1
plnsol,temp!用等值线或云图显示计算结果(这里是温度)
页:
[1]