热分析

分析日照等外界环境对某一混凝土箱梁的影响，所遇到的问题是如何施加和确立初始温度场。一 些参考文献中指出，凌晨6点左右箱梁的温度场为稳定温度场，因此，在分析前施加命令IC（TUNIT）命令，给全部单元施加温度20°，问题1，是不是用了这个命令之后，就把初始的温度场给确定了？下面是具体命令流，有做这个的虫友可以参考下，同时也希望大神给予指点和改正，谢谢！命令中有不明白之后，可以留言，尽力解答！！
FINISH
/CLEAR
KEYW,PR_THERM,1
/PREP7
ET,1,PLANE35
MP,EX,1,3.45E11
MP,KXX,1,1.54
MP,C,1,920
MP,DENS,1,2500
K,1,0,0
K,2,4.25
K,3,4.25,5.27
K,4,6.785,5.67
K,5,6.785,5.97
K,6,-3.625,5.97
K,7,-3.625,5.67
K,8,0,5.27
L,1,2
L,2,3
L,3,4
L,4,5
L,5,6
L,6,7
L,7,8
L,8,1
K,9,1.4,0.68
K,10,2.85,0.68
K,11,3.75,1.68
K,12,3.75,4.34
K,13,2.85,5.14
K,14,1.4,5.14
K,15,0.5,4.34
K,16,0.5,1.68
L,9,10
L,10,11
L,11,12
L,12,13
L,13,14
L,14,15
L,15,16
L,16,9
/PNUM,LINE,1
LPLOT
AL,9,10,11,12,13,14,15,16
AL,1,2,3,4,5,6,7,8
ASBA,2,1
NUMCMP,ALL
NUMMRG,ALL
LESIZE,ALL,0.2
MSHAPE,2D
MSHKEY,0
MAT,1
AMESH,1
FINISH
/SOLU
ANTYPE,TRANS
TRNODT,FULL
LUMPM,0
HEMIOPT,30
IC,ALL,TEMP,20(TUNIF,20）！！！！这个用哪个比较好？有什么区别？
LSEL,S,,,9,16,1
SFL,ALL,CONV,5.8,,27.2
LSEL,S,,,5
SFL,ALL,CONV,14.8,,39.0
LSEL,S,,,6,8,1
SFL,ALL,CONV,14.8,,38.6
LSEL,S,,,2,4,1
SFL,ALL,CONV,14.8,,37.7
LSEL,S,,,1
SFL,ALL,CONV,14.8,,31.3
ALLSEL
SFTRAN
AUTOTS,ON
TIMINT,ON
KBC,1
TIME,3600
NSUBST,20,40,20
OUTRES,ALL,ALL
SOLVE
FINISH
/POST1
SET,LAST
/SHOW,WIN32C
/CONTOUR,1,15
PLNSOL,TEMP

TB,MISO,1,5,,    !Tb, lab, mat, ntemp,npts,tbopt,eosopt 定义非线性材料特性表
!TBTEMP,20        !Mat: 材料号   Ntemp: 数据的温度数   Npts: 对某一给定温度数据的点数
!TBPT,defi,0.0017627,352.54   !! Miso: 多段线性等向强化（最多100个点）

1  先检查一下材料数据是否正确，比如说：
    你的导热系数定义的温度最i高达到 1200 °，其它材料性能只定义到 750 °，如果温度超过 750 °，结果可能会有问题。
    看你的命令流，你使用的长度单位好像是 mm，如是，则你的密度值不对。因为对这种情况，密度单位应该是 T/mm^3，其中 T - 吨；因此你的密度值应该是 7.81e-9，不是 7.81e-6 （Kg/mm^3）；
    此外，看你的密度和弹性模量值，你的材料应该是钢，那么你定义的导热系数值也就不对了，小了好几个数量级。此外，你的对流换热系数也需要检查一下。
    材料性能出问题，多半是单位制不协调造成的，所以你需要仔细检查一下各材料属性的单位，看是否与：
       mm - T - s   单位制一致。

查看温度结果是否合理
    根据你定义的材料性能，估计只在对流换热的线段附近温度有明显变化，远一点的地方温度变化不会很大。
    如确实这样，需要好好检查和修改材料性能数据。

我定义的是45钢的，材料参数能查到的都用上了，我还想在问几个问题：
1  导热系数和比热部分我认为是在温度场用，不影响结构场，所以就定义的多了点，那是不是最好和结构场的应力应变温度对应？
2  我的单位是mm，密度的确对应的是kg/mm^3，我见参考书是这个量级，我就这样定义了，我可以改一下试试，导热系数单位是（W/mm.℃），这个应该没错，换算成一般的（W/m.K）就是乘以10*3了
3  对流换热系数是我自己仿照做实验的那个情况、尺寸算的（文献上没查到:sweat:），比一般的例子所用的钢材小了一个数量级，我理解是因为我的尺寸小，换热面小，可以吗？我再重新算一下
4  温度结果如图，到505°后就基本均匀分布了，你一说我也不敢肯定说正确不正确了
现在问题是我把各温度下本构关系直接在结构场拉伸，结果正确，那就是本构定义没问题了吧，就是一耦合温度场，读取了温度结果就不对，我改了载荷步大小，能计算，也有结果，就是附图，但是结果不正确，从开始颈缩处应力竟然一直还增大，而且比我定义的强度极限小，更改网格划分竟然出现中间缩回去的结果，我还应该再从哪里修改模拟了？我觉得还是哪里有问题，就是不知道再从哪里下手，稳态的我也算过了，基本上结果分布一样，谢谢，收敛图上传错了，忽略它

原因还是你的材料性能有问题：导热系数太小，热不能传到结构内部，仅对表面有影响；比热太大，结构温度难以变化。换热系数也偏小，热交换很少，也是的结构温度上升很小。
给你一个 45 号钢的材料热性能数据看看吧，未必准确，可以参考一下：
表1 45号钢的比热容. 热导率与温度关系 (长度单位 m； 质量单位 kg)
    T/℃        100     200     300    400     500     600     700    755     800     900     1000
Cp/(J/kg·℃)   480     498     524    560     615     700     854    1 064   806     637     602
λ/(W/m·℃)   43.53  40.44  38.13  36.02  34.16  31.98  28.66  2 5.14   26.49  25.92  24.02
  改变比热和导热系数的单位 (m -> mm;  Kg -> t : 吨)：
      并注意到：  J/kg/℃ = m^2/ s^2/℃ = 1e6 · mm^2/ s^2/℃
                          W/m/℃  = kg·m/s^3/℃ = t·mm/s^3/℃
    T/℃              100      200     300     400     500     600     700    755     800     900     1000
Cp/(J/kg·℃)       480      498     524     560     615     700     854    1 064   806     637      602  - 全部都要乘以 1e6
λ/(t·mm/s^3/℃) 43.53  40.44 3 8.13  36.02  34.16  31.98   28.66  25.14   26.49  25.92   24.02  - 全部不变
表2 45号钢换热系数与温度关系   (长度单位 m； 质量单位 kg)
T/℃                 50    100     200     300      350      400      500    600     700    800
h/(W/m^2·℃) 2000  3800   6000   13500  14000   12500   7000   4200   1500   500
改变换热系数单位：   (m -> mm;  Kg -> t : 吨)：
T/℃                     50     100      200       300        350      400       500      600      700     800
h/(W/mm^2·℃)   2e-3   3.8e-3   6e-3    13.5e-3   14e-3   12.5e-3   7e-3   4.2e-3   1.5e-3   5e-3

是不是把应力图也看成温度图了？按上面的数据不能计算:sweat:

你只给了一张温度图，不至于看错吧？！
给你的数据只是给你参考的，说明你的材料数据有问题。
不能计算又是出现什么情况？

额:sweat:忘了到底是怎么不能计算了:sweat::sweat:来来回回改参数改模型改的糊涂了，我改了个三维的差不多做出来了，反正出问题最多的提示是说某些单元过度扭曲，可能我拉的厉害而且没有定义蠕变和破坏准则吧。不过有个问题挺郁闷，我输入多条应力应变曲线计算几步后面就不收敛，要是一条的话，连刚开始颈缩的阶段都可以模拟，abaqus里面反正不能输入应力应变下降段，我查了半天也没查出来到底ansys能不能定义，反正是那样了

ANSYS 中好像也不能够输入应力 - 应变下降段；因为应力 - 应变 下降段意味着材料已经破坏了。