andylau436 发表于 2006-7-7 21:52:00

温度分布为什么会这样?请高手指点!

这是我的一个例子,是一个壁和底都是四层的不同材料组成的圆柱,里层是恒定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]
查看完整版本: 温度分布为什么会这样?请高手指点!