热分析C的问题

baochaoxiong

本人新手，最近做热分析，加了对流和外界温度，但是如果加上了比热容C材料属性，节点温度就不变，不加就会变化。命令如下：
/CONFIG,NRES,99999
/prep7                           !进入前处理

ET,1,PLANE55                     !选用2维热实体平面单元
MP,KXX,1,3.3                     !定义导热系数
MP,DENS,1,2500                   !密度

!*定义外部左边温度随时间变化函数
*DIM,wz,TABLE,4,1,1,time,temp,
*SET,wz(1,0,1) , 0  
*SET,wz(1,1,1) , 27   
*SET,wz(2,0,1) , 1.67
*SET,wz(2,1,1) , 33   
*SET,wz(3,0,1) , 8.67  
*SET,wz(3,1,1) , 33
*SET,wz(4,0,1) , 10.34
*SET,wz(4,1,1) , 34
!*结果显示
*vplot,wz(1,0),wz(1,1)
!*定义外部右边温度随时间变化函数
*DIM,wy,TABLE,4,1,1,time,temp,
*SET,wy(1,0,1) , 0  
*SET,wy(1,1,1) , 27   
*SET,wy(2,0,1) , 1.67
*SET,wy(2,1,1) , 29   
*SET,wy(3,0,1) , 8.67  
*SET,wy(3,1,1) , 35
*SET,wy(4,0,1) , 10.34
*SET,wy(4,1,1) , 40
!*结果显示
*vplot,wy(1,0),wy(1,1)
!*定义外部上下边温度随时间变化函数
*DIM,wsx,TABLE,4,1,1,time,temp,
*SET,wsx(1,0,1) , 0  
*SET,wsx(1,1,1) , 27   
*SET,wsx(2,0,1) , 1.67
*SET,wsx(2,1,1) , 31   
*SET,wsx(3,0,1) , 8.67  
*SET,wsx(3,1,1) , 34
*SET,wsx(4,0,1) , 10.34
*SET,wsx(4,1,1) , 37
!*结果显示
*vplot,wsx(1,0),wsx(1,1)
!*定义外部内部温度随时间变化函数
*DIM,n,TABLE,4,1,1,time,temp,
*SET,n(1,0,1) , 0  
*SET,n(1,1,1) , 28   
*SET,n(2,0,1) , 1.67
*SET,n(2,1,1) , 30   
*SET,n(3,0,1) , 8.67  
*SET,n(3,1,1) , 32
*SET,n(4,0,1) , 10.34
*SET,n(4,1,1) , 32
!*结果显示
*vplot,n(1,0),n(1,1)


k,1,0,0,0                        !建立节点
k,2,7.6,0,0
k,3,7.6,-0.2,0
k,4,7.1,-0.25
k,5,5.9,-0.57,0
k,6,5.9,-6.032,0
k,7,1.7,-6.032,0
k,8,1.7,-0.57,0
k,9,0.5,-0.25,0
k,10,0,-0.2,0
k,11,3.1,-0.32,0
k,12,4.5,-0.32,0
k,13,5.25,-0.57,0
k,14,5.25,-4.852,0
k,15,5.05,-5.252,0
k,16,2.55,-5.252,0
k,17,2.35,-4.852,0
k,18,2.35,-0.57,0
l,1,2                            !由点生成线
l,2,3
l,3,4
l,4,5
l,5,6
l,6,7
l,7,8
l,8,9
l,9,10
l,10,1
l,11,12
l,12,13
l,13,14
l,14,15
l,15,16
l,16,17
l,17,18
l,18,11
lplot
a,1,2,3,4,5,6,7,8,9,10           !由线连成面
a,11,12,13,14,15,16,17,18
aplot
asba,1,2                         !差集
/replot
CM,_Y,AREA  
ASEL,,,,3
CM,_Y1,AREA
CMSEL,S,_Y  
CMSEL,S,_Y1
AATT,1,,1,0,  
CMSEL,S,_Y  
CMDELE,_Y   
CMDELE,_Y1
FLST,2,1,5,ORDE,1   
FITEM,2,3   
AESIZE,P51X,0.1,
MSHAPE,0,2D
MSHKEY,0
CM,_Y,AREA  
ASEL,,,,3
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
AMESH,_Y1   
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  

/SOL
ANTYPE,4                         !选择类型：瞬态分析
TIMINT,off                       !首先进行稳态分析，确定初始条件
TINTP,0.005
TIME,0.01                        !时间  
AUTOTS,-1                        !自动时间步长
DELTIM,0.01                      !时间步长设为0.01   

lsel,s,,,7,10
nsll,s,
sf,all,conv,13,%WZ%,             !赋予左边界对流系数
lsel,s,,,2,5
nsll,s,
sf,all,conv,13,%WY%,             !赋予右边界对流系数  
lsel,s,,,11,18
nsll,s,
sf,all,conv,13,%N%,              !赋予内边界对流系数
lsel,s,,,1
lsel,a,,,6
nsll,s,
sf,all,conv,13,%WSX%,            !赋予上下边界对流系数


FLST,2,77,1,ORDE,4               !wz边界条件赋予外面左边边界   
FITEM,2,1   
FITEM,2,154
FITEM,2,196
FITEM,2,-270
D,P51X,,%WZ%,,,,TEMP,,,,,         
FLST,2,77,1,ORDE,3               !wy边界条件赋予外面右边边界  
FITEM,2,2   
FITEM,2,78  
FITEM,2,-153
D,P51X,,%WY%,,,,TEMP,,,,,         
FLST,2,114,1,ORDE,4              !wsx边界条件赋予外面上下边边界
FITEM,2,4   
FITEM,2,-77
FITEM,2,155
FITEM,2,-194
D,P51X,,%WSX%,,,,TEMP,,,,,      
FLST,2,151,1,ORDE,2              !n边界条件赋予内部边界
FITEM,2,271
FITEM,2,-421
D,P51X,,%N%,,,,TEMP,,,,,

/pbc,temp,,1
/psf,conv,,2
allsel
solve                            !得到初始温度分布

time,12                          !进行瞬态分析
timint,on                        !打开时间积分
deltim,0.1,0.01,0.25             !设置时间步长，最大及最小时间步长
autots,on                        !打开自动时间步长
ddelet,all,temp                  !删除稳态分析中定义的接点温度
outres,all,all                   !将每个子步的值写入数据库文件
solve
finish                           !得到rth文件
save

/post26                          !进入post26后处理
solu,2,dtime,dtime               !2~每一子步采用的时间步长
nsol,3,430,temp,,T_one           !3~一号测点
nsol,4,1218,temp,,T_two          !4~二号测点
nsol,5,559,temp,,T_three         !5~三号测点
nsol,6,422,temp,,T_four          !6~四号测点
nsol,7,718,temp,,T_five          !7~五号测点
nsol,8,425,temp,,T_six           !8~六号测点
nsol,9,938,temp,,T_ten           !9~十号测点
Plvar,2
plvar,3,4,5,6,7,8,9
finish

!进入POST1后处理
/post1     
set,last                         !设置为最后一个子步
plnsol,temp                      !温度等值线图
finish

/prep7
etchg,tts                        !将热单元转化为PLANE42结构单元
MP,EX,1,3.5e10                   !弹性模量
MP,PRXY,1,0.167                  !泊松比
MP,ALPX,1,1e-5                   !热膨胀系数
MP,C,1,925                       !比热
finish

/solu                            !进入结构分析求解
ldread,temp,,,,,,rth             !读入节点温度文件rth
/pbc,all,1
/psf,pres,,1
/pbf,temp,,1
eplot
solve
finish

!进入POST1后处理
/post1                           !设置为最后一个子步
set,last                     
esel,s,mat,,1
nsle,s
plnsol,s,EQV                     !等效应力等值线图
finish

baochaoxiong

如果材料加上C=925，节点温度就没有变化基本上