同样的热力耦合问题用PLANE13和PLANE223分别计算,出现了不同的结果,并且用PLANE223单元分别计算了强耦合和弱耦合两种情况,发现结果基本一样,不知道应该以哪个为准?从趋势上看PLANE13计算的结果是正确的,但不是应该八节点的PLANE223单元应该更精确些吗?我应该相信哪个结果呢?在面上沿径向分别取了5个点画出了随时间的变化情况(见附件),请大家帮忙分析一下。
命令流如下:
/filename,oxidescale
/title,temperature change along the radius
/PREP7
!选择单元类型
ET,1,PLANE223,11
KEYOPT,1,2,1
!强耦合(KEYOPT =0)弱耦合(KEYOPT =1)
KEYOPT,1,3,2
!平面应变
!设置材料属性
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPDATA,EX,1,,2.18e11
MPDATA,EX,1,,2.13e11
MPDATA,EX,1,,2.06e11
MPDATA,EX,1,,1.99e11
MPDATA,EX,1,,1.91e11
MPDATA,EX,1,,1.81e11
MPDATA,EX,1,,1.70e11
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPDATA,PRXY,1,,0.3
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,7840
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
UIMP,1,REFT,,,517
MPDATA,ALPX,1,,1.1e-5
MPDATA,ALPX,1,,1.2e-5
MPDATA,ALPX,1,,1.3e-5
MPDATA,ALPX,1,,1.3e-5
MPDATA,ALPX,1,,1.4e-5
MPDATA,ALPX,1,,1.4e-5
MPDATA,ALPX,1,,1.4e-5
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPDATA,C,1,,460
MPDATA,C,1,,500
MPDATA,C,1,,500
MPDATA,C,1,,540
MPDATA,C,1,,630
MPDATA,C,1,,710
MPDATA,C,1,,800
MPTEMP,,,,,,,,
MPTEMP,1,20
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,600
MPDATA,KXX,1,,34.8
MPDATA,KXX,1,,37.3
MPDATA,KXX,1,,38.5
MPDATA,KXX,1,,38.5
MPDATA,KXX,1,,37.3
MPDATA,KXX,1,,34.8
MPDATA,KXX,1,,32.7
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,2.1E+011
MPDATA,PRXY,2,,0.29
MPTEMP,,,,,,,,
MPTEMP,1,50
MPTEMP,2,100
MPTEMP,3,200
MPTEMP,4,300
MPTEMP,5,400
MPTEMP,6,500
MPTEMP,7,530
MPTEMP,8,600
UIMP,2,REFT,,,517
MPDATA,ALPX,2,,11e-6
MPDATA,ALPX,2,,12e-6
MPDATA,ALPX,2,,13e-6
MPDATA,ALPX,2,,14e-6
MPDATA,ALPX,2,,16e-6
MPDATA,ALPX,2,,24e-6
MPDATA,ALPX,2,,25e-6
MPDATA,ALPX,2,,16e-6
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,3598
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,2,,6
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,2,,630
!创建几何模型
PCIRC,0.0175,0.0174,0,90,
PCIRC,0.0255,0.0175,0,90,
AGLUE,ALL
NUMCMP,AREA
NUMCMP,LINE
!划分网格
/PNUM,LINE,1
/PNUM,AREA,1
/NUMBER,0
/REPLOT
LSEL,S,LINE,,2,4,2
LESIZE,ALL,,,2
LSEL,S,LINE,,1,3,2
LESIZE,ALL,,,50
MAT,2
AMESH,1
APLOT
LSEL,S,LINE,,6,7
LESIZE,ALL,,,15,5
LSEL,S,LINE,,5
LESIZE,ALL,,,50
MAT,1
AMESH,2
CSWPLA,11,1,1,1, !创建圆柱局部坐标
CSYS,1
ALLSEL
NROTAT,ALL
LSEL,S,LINE,,3
!耦合内壁径向位移
NSLL,S,1
CP,8,UX,ALL
LSEL,S,LINE,,5
!耦合外壁径向位移
NSLL,S,1
CP,9,UX,ALL
ALLSEL
LPLOT
FINISH
!加载求解
/SOL
ANTYPE,TRANS
!指定分析类型为瞬态分析
TIMINT,OFF
!先做稳态分析,确定初始条件
TIME,0.01
!指定计算终止时间
DELTIM,0.01,0,0
!指定时间步长
!TREF,517,
!设置参考温度
ALLSEL
!设置位移自由度约束
D,ALL,UY,0
LSEL,S,LINE,,3
!施加压力载荷
SFL,ALL,PRES,1. 7E7,
!施加内外壁对流载荷
SFL,ALL,CONV,4359, ,500,
LSEL,S,LINE,,5
SFL,ALL,CONV,125, ,932,
ALLSEL
SOLVE
!求解,得到初始温度场和应力场分布
TIMINT,ON
!打开时间积分选项,进行瞬态分析
TIME,120
!指定计算终止时间
KBC,1
!指定加载方式为阶跃载荷
NSUBST,60,240,30
!设置时间子步数
AUTOTS,ON
OUTRES,ALL,ALL
LSEL,S,LINE,,3
SFL,ALL,CONV,4359, ,400, !修改对流边界条件,水蒸汽温度降低100
ALLSEL
SOLVE
!求解
SAVE,'oxidescale', 'db', |