- 积分
- 0
- 注册时间
- 2006-4-15
- 仿真币
-
- 最后登录
- 1970-1-1
|
finish
/clear,nostart
KEYW,PR_CFD,1
/filname,THERMAL !指定文件名
/title,2D THERMAL SIMULATION !指定标题
/units,SI !国际单位制
!*
!*
!*定义截面几何参数
!定义入口尺寸
*set,INLET_HEIGHT,200E-3
*set,INLET_WIDTH,758E-3
*set,INLET_THICK,2E-3
*set,INLET_LENGTH,243E-3
!定义出口尺寸
*set,OUTLET_HEIGHT,200E-3
*set,OUTLET_WIDTH,758E-3
*set,OUTLET_THICK,2E-3
*set,OUTLET_LENGTH,243E-3
!定义热通道底尺寸
*set,BOT_WIDTH,INLET_WIDTH
*set,BOT_HEIGHT,10E-3
!定义热通道顶尺寸
*set,UP_WIDTH,OUTLET_WIDTH
*set,UP_HEIGHT,10E-3
!定义外片玻璃的尺寸
*set,WB_WIDTH,15E-3
*set,WB_HEIGHT,7111E-3
!定义内片玻璃的尺寸
*set,NB_WIDTH,24E-3
*set,NB_HEIGHT,7535E-3
!定义热通道空气层厚度
*set,A_WIDTH,500E-3
*SET,A_HEIGHT,WB_HEIGHT
!定义入口百叶中心间距及边距
*set,BY_DIST,40E-3
*set,BY_EDGE,20E-3
!定义通风口尺寸
*set,TFK_DIST,20E-3
*set,TFK_WIDTH,10E-3
*set,TFK_EDGE,50E-3
*set,TFK_HEIGHT,INLET_THICK
!*
!*设置材料及边界参数
*set,KXX_WB,0.76
*set,RDSF_WB,0.837
*set,HF_WBOUT,17.79
*set,KXX_NB,0.76
*set,RDSF_NB,0.837
*set,HF_NBIN,8
*set,KXX_AL,160
*set,RDSF_AL,0.9
!内热源
*SET,HGEN_WB,171.87/0.015
*SET,HGEN_NB,158.32/0.024
!设置入口温度和流速
*SET,V_INLET,0.22
*SET,TEMP_INLET,33.2
!设置出口的压强
*SET,PRES_OUTLET,0
!*进入前处理,建立几何模型
/PREP7
!定义单元为二维热流体单元
ET,1,FLUID141
!定义材料参数
MP,KXX,2,KXX_WB
MP,KXX,3,KXX_NB
MP,KXX,4,KXX_BY
!卡式坐标原点为模型左下角点
CSYS,4 !设定活动坐标系为工作平面
!建立几何模型
!(A1-A5)
WPOFF,INLET_LENGTH+WB_WIDTH/2,BOT_HEIGHT+BY_EDGE
k,1,-15E-3,-3E-3
K,2,-15E-3,3E-3
K,3,-10E-3,5E-3
K,4,10E-3,5E-3
K,5,15E-3,3E-3
K,6,15E-3,-3E-3
K,7,10E-3,-5E-3
K,8,-10E-3,-5E-3
FLST,2,8,3
FITEM,2,1
FITEM,2,2
FITEM,2,3
FITEM,2,4
FITEM,2,5
FITEM,2,6
FITEM,2,7
FITEM,2,8
A,P51X
AGEN,5,1,,,,BY_DIST
CSYS,0
WPAVE,0,0
CSYS,4
!(A6)
RECTNG,0,BOT_WIDTH,0,UP_HEIGHT
!(A7)
RECTNG,0,INLET_WIDTH,BOT_HEIGHT,BOT_HEIGHT+INLET_HEIGHT
!(A8)
RECTNG,0,INLET_WIDTH,BOT_HEIGHT+INLET_HEIGHT,BOT_HEIGHT+INLET_HEIGHT+INLET_THICK
!偏移工作面
WPOFF,INLET_LENGTH,BOT_HEIGHT+INLET_HEIGHT+INLET_THICK
!(A9)
RECTNG,0,WB_WIDTH,0,WB_HEIGHT
!(A10)
RECTNG,WB_WIDTH,WB_WIDTH+A_WIDTH,0,A_HEIGHT
CSYS,0
WPAVE,0,0
CSYS,4
!(A11)
RECTNG,BOT_WIDTH,BOT_WIDTH+NB_WIDTH,0,NB_HEIGHT
!(A12-A31)
WPOFF,INLET_LENGTH+WB_WIDTH,BOT_HEIGHT+INLET_HEIGHT+INLET_THICK
RECTNG,TFK_EDGE,TFK_EDGE+TFK_WIDTH,-INLET_THICK,0
AGEN,10,12,,,TFK_DIST
AGEN,2,12,21,1,10*TFK_DIST+TFK_WIDTH
AGEN,2,8,,,,INLET_THICK+WB_HEIGHT
AGEN,2,6,,,,BOT_HEIGHT+INLET_HEIGHT*2+INLET_THICK*2+WB_HEIGHT
AGEN,2,7,,,,INLET_HEIGHT+INLET_THICK*2+WB_HEIGHT
AGEN,2,1,5,1,,WB_HEIGHT+INLET_THICK*2+INLET_HEIGHT
AGEN,2,12,31,1,,WB_HEIGHT+INLET_THICK
ALLSEL,ALL
AOVLAP,ALL
/PNUM,AREA,1
ALLSEL,ALL
APLOT
ALLSEL,ALL
AGLUE,ALL
APLOT
!将流体面积加起来,赋予材料属性及单元类型
ALLSEL,ALL
ASEL,S,AREA,,104,106,1
ASEL,A,AREA,,12,31,1
ASEL,A,AREA,,40,59,1
AATT,1,,1
AADD,ALL
!
ALLSEL,ALL
ASEL,S,AREA,,9
AATT,2,,1
!
ALLSEL,ALL
ASEL,S,AREA,,103
AATT,3,,1
ALLSEL,ALL
APLOT,ALL
ALLSEL,ALL
ASEL,S,AREA,,7,9,2
ASEL,A,AREA,,103
ASEL,INVE
APLOT
AATT,4,,1
SMRTSIZE,6
ALLSEL,ALL
AMESH,ALL
!*设置边界条件
ALLSEL,ALL
ASEL,S,AREA,,9
APLOT
BFA,ALL,HGEN,HGEN_WB
!
ALLSEL,ALL
ASEL,S,AREA,,103
APLOT
BFA,ALL,HGEN,HGEN_NB
/PNUM,LINE,1
/PNUM,AREA,0
/NUMBER,0
LPLOT
!*
!*
!施加对流换热系数
ALLSEL,ALL
LSEL,S,LINE,,56
SFL,ALL,CONV,HF_WBOUT,,33.2
!******************************此处有问题********************
!施加面面辐射换热系数
ALLSEL,ALL
LSEL,S,LINE,,18
SFL,ALL,RDSF,RDSF_WB,,1
!施加面面辐射换热系数
ALLSEL,ALL
LSEL,S,LINE,,370
SFL,ALL,RDSF,RDSF_NB,,1
!******************************请大侠帮忙*********************
!施加对流换热系数
ALLSEL,ALL
LSEL,S,LINE,,62
SFL,ALL,CONV,HF_NBIN,,25
!设置入口的流速和温度
ALLSEL,ALL
LSEL,S,LINE,,372
DL,ALL,,VX,V_INLET
DL,ALL,,TEMP,TEMP_INLET
!设置出口的压强
ALLSEL,ALL
LSEL,S,LINE,,373
DL,ALL,,PRES,PRES_OUTLET
!*设置求解选项
ALLSEL,ALL
FLDATA1,SOLU,FLOW,T !设置求解选项包括流场
FLDATA1,SOLU,TEMP,T !设置求解选项包括温度场
FLDATA1,SOLU,TURB,T !设置求解选项为湍流
FLDATA2,ITER,EXEC,9000 !设置求解迭代次数
FLDATA13,VARY,DENS,T
FLDATA7,PROT,DENS,AIR-SI !设置流体密度
FLDATA7,PROT,VISC,AIR-SI !设置流体的粘度
FLDATA7,PROT,COND,AIR-SI !设置流体的导热系数
FLDATA7,PROT,SPHT,AIR-SI !设置流体的比热容
!*
ACEL,0,9.81,0 ! Set flow environment gravity ,y direction
FLDATA24,TURB,MODL,2 !设置湍流模型为零方程模拟
FLDATA24,TURB,RATI,2 !设置湍流比为2
FLDATA18,METH,TEMP,4 !设置温度方程的计算方法为PGMR
FINISH
/SOLU
CSYS,0
WPAVE,0,0
ALLSEL,ALL
APLOT
!设置辐射计算的参数
ALLSEL,ALL
STEF,5.67E-8 !设置玻尔兹曼常数
TOFFST,273.15 !摄式温度与开尔文温度的差值
IC,ALL,TEMP,25 !设置初始温度
ALLSEL,ALL
NSUBST,100
SOLVE !求解FLOTRAN CFD
FINISH
/POST1
!读入最后一步的计算结果
SET,LAST
!查看总体温度分布图
/title,TEMPERATURE DISTRIBUTION
PLNSOL,TEMP
!查看流体流速矢量图
/title,VELOCITY DISTRIBUTION
PLVECT,V,,,,VECT,ELEM,ON,0
!*查看表面换热系数
/EFACET,1
PLNSOL, HFLM,, 0 |
|