命令流如下:
finish
/clear,nostart
KEYW,PR_CFD,1
/filname,THERMAL ANALYSIS !指定文件名
/title,2D THERMAL SIMULATION OF DOUBLE Curtain WALL !指定标题
/units,SI !国际单位制
!*
!*
!*定义模型的几何参数
*set,WB_WIDTH,24E-3
*set,WB_HEIGHT,4
*set,AIR_OUT_WIDTH,109E-3
*set,AIR_OUT_HEIGHT,4
*set,BY_WIDTH,1E-3
*set,BY_HEIGHT,4
*set,AIR_IN_WIDTH,109E-3
*set,AIR_IN_HEIGHT,4
*set,NB_WIDTH,8E-3
*set,NB_HEIGHT,4
*set,KXX_AIR,0.024
*set,DENS_AIR,1.11
*set,VISC_AIR,1.7e-5
*set,KXX_WB,0.0564
!*
*set,KXX_NB,0.9987
*set,KXX_BY,160
!*
*SET,V_INLET,0.04
*SET,TEMP_INLET,24
!*
*SET,PRES_OUTLET,0
/PREP7
!定义单元为二维热流体单元
ET,1,FLUID141
KEYOPT,1,4,1
!定义材料参数
MP,KXX,1,KXX_AIR
MP,KXX,2,KXX_WB
MP,C,2,840
MP,KXX,3,KXX_NB
MP,C,3,840
MP,KXX,4,KXX_BY
MP,C,4,2660
!卡式坐标原点为模型左下角点
CSYS,4 !设定活动坐标系为工作平面
!建立几何模型
RECTNG,0,WB_WIDTH,0,WB_HEIGHT
WPOFF,WB_WIDTH
RECTNG,0,AIR_OUT_WIDTH,0,AIR_OUT_HEIGHT
WPOFF,AIR_OUT_WIDTH
RECTNG,0,BY_WIDTH,0,BY_HEIGHT
WPOFF,BY_WIDTH
RECTNG,0,AIR_IN_WIDTH,0,AIR_IN_HEIGHT
WPOFF,AIR_IN_WIDTH
RECTNG,0,NB_WIDTH,0,NB_HEIGHT
WPOFF,-(WB_WIDTH+AIR_OUT_WIDTH+BY_WIDTH+AIR_IN_WIDTH)
!*
ALLSEL,ALL
/PNUM,AREA,1
APLOT
!*
ALLSEL,ALL
ASEL,S,AREA,,2,4,2
AATT,1,,1
!*
ALLSEL,ALL
ASEL,S,AREA,,1
AATT,2,,1
ALLSEL,ALL
ASEL,S,AREA,,5
AATT,3,,1
ALLSEL,ALL
ASEL,S,AREA,,3
AATT,4,,1
!*
ALLSEL,ALL
/PNUM,AREA,0
/PNUM,LINE,1
WPSTYLE,,,,,,,,0
/TRIAD,OFF
LPLOT
!*
ALLSEL,ALL
LSEL,S,LINE,,2,20,2
LESIZE,ALL,,,300
!*
ALLSEL,ALL
LSEL,S,LINE,,1,3,2
LESIZE,ALL,,,10
!*
ALLSEL,ALL
LSEL,S,LINE,,5,7,2
LSEL,A,LINE,,13,15,2
LESIZE,ALL,,,50
!*
ALLSEL,ALL
LSEL,S,LINE,,9,11,2
LESIZE,ALL,,,2
!*
ALLSEL,ALL
LSEL,S,LINE,,17,19,2
LESIZE,ALL,,,5
!*
ALLSEL,ALL
MSHKEY,1
MSHAPE,0,2D
AMESH,ALL
!*
ALLSEL,ALL
LPLOT
!*
!*********标记流固耦合面************
!*
ALLSEL,ALL
LSEL,S,LINE,,2,8,6
NSLL,S,1
SF,ALL,FSIN,1
!*
ALLSEL,ALL
LSEL,S,LINE,,6,12,6
NSLL,S,1
SF,ALL,FSIN,2
!*
ALLSEL,ALL
LSEL,S,LINE,,10,16,6
NSLL,S,1
SF,ALL,FSIN,3
!*
ALLSEL,ALL
LSEL,S,LINE,,14,20,6
NSLL,S,1
SF,ALL,FSIN,4
!*
ALLSEL,ALL
LSEL,S,LINE,,4
NSLL,S,0
SF,ALL,CONV,17.6,,32
ALLSEL,ALL
LSEL,S,LINE,,18
NSLL,S,0
SF,ALL,CONV,8,,24
!*
!设置入口的流速和温度
ALLSEL,ALL
LSEL,S,LINE,,5,13,8
NSLL,S,0
D,ALL,VY,V_INLET
D,ALL,TEMP,TEMP_INLET
!*
!设置出口的压强
ALLSEL,ALL
LSEL,S,LINE,,7,15,8
NSLL,S,0
D,ALL,PRES,PRES_OUTLET
!*
!*设置求解选项
ALLSEL,ALL
FLDATA1,SOLU,FLOW,T !设置求解选项包括流场
FLDATA1,SOLU,TEMP,T !设置求解选项包括温度场
FLDATA1,SOLU,TURB,T !激活湍流模型
FLDATA1,SOLU,ALE,T
!*
FLDATA2,ITER,EXEC,90000 !设置求解迭代次数
FLDATA3,TERM,PRES,1E-5 !指定压强收敛精度
FLDATA3,TERM,TEMP,1E-5 !指定温度收敛精度
!*
FLDATA7,PROT,DENS,AIR-SI !设置流体密度
!FLDATA13,VARY,DENS,T
FLDATA7,PROT,VISC,AIR-SI !设置流体的粘度
FLDATA7,PROT,COND,AIR-SI !设置流体的导热系数
FLDATA7,PROT,SPHT,AIR-SI !设置流体的比热容
!*
FSAN,ON
FSCO,ALL,0.1
FSRE,ALL,0.5
FSOR,FLUID
FSTR,FLUID,TRAN
FSTR,SOLID,TRAN
!*
FSTI,7200
FSDT,1
FSIT,100
FSIN,CONS
!*
!*
FLDATA5,OUTP,HFLM,T
!*
FINISH
/SOLU
!*
ANTYPE,TRANS
!*
ALLSEL,ALL
SOLVE !求解FLOTRAN CFD
!*
FINISH
/POST1
!*
!读入最后一步的计算结果
SET,LAST
!*
!查看总体温度分布图
ALLSEL,ALL
/title,TEMP DISTRIBUTION
PLNSOL,TEMP
!*
出错提示信息:
*** ERROR *** CP = 10.922 TIME= 11:39:04
Number of FLUID/SOLID interfaces do not match for FSIN boundary
condition. |