热分析

类似于冷却问题，先用稳态求解加载温度载荷，然后瞬态求解，先删掉之前的温度载荷，施加对流换热，把之前算的结果作为此时瞬态的开始数据，观察温度变化，其中1图示稳态的结果，二图是增加对流换热的冷却结果、
问题：我设置的时间很短的，不应该全是室温呀，应该有一部分高于室温的，怎么全是红的？
以下是我瞬态程序：
/SOLU
ALLSEL,ALL
DDELET,ALL,TEMP
TIME,0.000001
DELTIM,0.0000001,0.0000001,0.0000001
ALLSEL,ALL
OUTRES,ERASE
OUTRES,ALL,1
ESEL,S,TYPE,,2
NSLE
SF,ALL,CONV,20,20
ALLSEL,ALL
/SOLU
ANTYPE,4
TRNOPT,FULL
LUMPM,0
LDREAD,TEMP,,,O.O1,2,'TEMPCHANGE','RTH',''
SOLVE

你提供的命令流中并没有稳态热分析的内容，由于稳态分析后离开了，再进入  SOLUTION 时，后续计算被当作新的分析，因而没有从稳态分析结果的基础上继续。

正确的做法是：
    在 SOLUTION 中，做完稳态分析后，不要离开 SOLUTION，删除原来的热载荷，增加对流载荷，然后立即做瞬态分析 (不要再读取稳态分析的结果，因为现在结构中的温度就是刚才计算得到的稳态温度分布)。

不离开solution，删除载荷是什么意思呢？如果命令流不复杂，你可以稍微用命令流做一下吗？

大概意思如下：
/solu
!  第一个载荷步
定义第一个载荷步的求解设置，施加载荷、约束，等
求解第一个载荷步
!*
!  第二个载荷步
删除第一个载荷步的载荷
定义第二个载荷步的求解设置，施加载荷 (约束不变)，等
求解第二个载荷步
......
求解完所有载荷步后在进行后处理

/PREP7
TOFFST,273
ET,1,PLANE55
ET,2,SURF151
KEYOPT,2,1,0
KEYOPT,2,2,0
KEYOPT,2,3,0
KEYOPT,2,4,0
KEYOPT,2,5,0
KEYOPT,2,6,0
KEYOPT,2,7,0
KEYOPT,2,8,2
KEYOPT,2,9,0
MP,DENS,1,4510
MPTEMP,1,20,100,200,300,400,500
MPTEMP,7,1000,1500,2500,,,
MPDATA,KXX,1,1,6.8,7.4,8.7,9.8,10.3,11.8
MPDATA,KXX,1,7,15.5,22.1,22.2,,,
MPTEMP,1,20,100,200,300,400,500
MPTEMP,7,1000,1500,2500,,,
MPDATA,C,1,1,611,624,653,674,691,703
MPDATA,C,1,7,1030,1850,1852,,,
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,2,,1E-9
MPDATA,KYY,2,,1E+9
MPDATA,KZZ,2,,1E-9
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,2,,1/4510
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,4510
MP,DENS,3,4110
MPTEMP,1,20,100,200,300,400,500
MPTEMP,7,1000,1500,2500,,,
MPDATA,KXX,3,1,6.8,7.4,8.7,9.8,10.3,11.8
MPDATA,KXX,3,7,15.5,22.1,22.2,,,
MPTEMP,1,20,100,200,300,400,500
MPTEMP,7,1000,1500,2500,,,
MPDATA,C,3,1,611,624,653,674,691,703
MPDATA,C,3,7,1030,1850,1852,,,
MPTEMP,1,20,100,200,300,400,500
MPTEMP,7,1000,1500,2500,,,
MPDATA,ENTH,3,1,7.8886E-31,2.216448E+08,5.218776E+08,8.379168E+08,1.1658552E+09,1.4982336E+09
MPDATA,ENTH,3,7,4.481736E+09,1.215672E+10,2.03883468E+10,,,
/REPLOT
RECTNG,0,0.002,0,0.002
RECTNG,0.0008,0.0012,0.002-5*1E-6,0.002
RECTNG,0.0008,0.0012,0,0.002-5*1E-6
APTN,ALL
ALLSEA,ALL
ASEL,S,AREA,,5
AATT,1,1,1
ALLSEA,ALL
ASEL,S,AREA,,6
AATT,1,1,1
ALLSEA,ALL
ASEL,S,AREA,,7
AATT,1,1,1
ALLSEA,ALL
ASEL,S,AREA,,4
AATT,3,1,1
ALLSEL,ALL
ESIZE,10E-06
MSHAPE,0,2D
MSHKEY,2
AMESH,ALL
ALLSEL,ALL
LSEL,S,LINE,,2
NSLL,S,
TYPE,2
ESURF,ALL
LSEL,S,LINE,,4
NSLL,S,
TYPE,2
ESURF,ALL
LSEL,S,LINE,,7
NSLL,S,
TYPE,2
ESURF,ALL
LSEL,S,LINE,,9
NSLL,S,
TYPE,2
ESURF,ALL
LSEL,S,LINE,,15
NSLL,S,
TYPE,2
ESURF,ALL
LSEL,S,LINE,,16
NSLL,S,
TYPE,2
ESURF,ALL
LSEL,S,LINE,,17
NSLL,S,
TYPE,2
ESURF,ALL
LSEL,S,LINE,,18
NSLL,S,
TYPE,2
ESURF,ALL
ASEL,S,AREA,,5
NSLA,,1
D,ALL,TEMP,20
ASEL,S,AREA,,6
NSLA,,1
D,ALL,TEMP,20
ASEL,S,AREA,,7
NSLA,,1
D,ALL,TEMP,20
ASEL,S,AREA,,4
NSLA,,1
D,ALL,TEMP,2000
TREF,20
ALLSEL,ALL
/SOLU
ANTYPE,TRANS
TIMINT,OFF
TIME,0.01
DELTIM,0.01
SOLVE
/POST1
SET,LAST
PLNSOL,TEMP          ！上面没有问题
/SOLU
ALLSEL,ALL
DDELET,ALL,TEMP
TIME,0.000001
DELTIM,0.0000001,0.0000001,0.0000001
ALLSEL,ALL
OUTRES,ERASE
OUTRES,ALL,1
ESEL,S,TYPE,,2
NSLE
SF,ALL,CONV,20,20
ALLSEL,ALL
/SOLU
ANTYPE,4
TRNOPT,FULL
LUMPM,0
LDREAD,TEMP,,,O.O1,2,'TEMPCHANGE','RTH',''
SOLVE
/POST1
SET,LAST
PLNSOL,TEMP
这是我的程序，你看看吧，我觉得想法应该没问题，就是不知道哪里有问题。谢谢你呀，我有点急。

因为我要做一百组数据，所以用apdl比较方便。麻烦你了。

/POST1
SET,LAST
PLNSOL,TEMP          ！上面没有问题
这几行不要
下一步前加一行
lsclear,all

我用的DDELET,ALL,TEMP,这两个有什么区别吗？