P06_Chap_01_一般热传递分析.doc

midas FEA第六部分 midas FEA 第六部分 热传递分析 midas FEA第一章 midas FEA 第一章 | 一般热传递分析 We Analyze and Design the Future PAGE We Analyze and Design the Future PAGE 384 PAGE 385 PAGE 385 We Analyze and Design the Future 第一章 一般热传递分析1-1 第一章 一般热传递分析 在高温下，结构随时间的变化会发生较大的温度变化和梯度变化；热传递分析就是通过计算结构温度场随时间的变化来分析结构的热传递特性以及温度变化对结构的影响(变形、应力)；在短时间内浇注大量混凝土时由于产生的水化热或因为火灾等原因导致结构局部加热就是典型需要进行热传递分析的例子。 热传递主要分为热传导、对流以及辐射，热传导是指大量的原子、分子的相互碰撞而使热能由温度较高的部分传向温度较低部分的过程，液体或者气体内的传递叫做对流；辐射是指受热物体以电磁辐射的形式向外界发送并传送能量的过程。 midas FEA是以固体里常用的传导理论为基础，用有限元的方法进行热传递分析，对于对流用边界条件和荷载的形式来考虑。 1-2 热传递方程式 传导里的热传递方程式遵循能量守恒定律，比如一定的体积和其周边的面积为时能量守恒如下： 热增加率 () = 热流量 () + 发热率() 体积里的热增加率根据比热和密度按照如下公式计算： (1.2.1) 其中： : 比热 () : 密度 () 对于热传导传递的的热流量可根据傅立叶原理利用温度梯度和热传导率来计算： (1.2.2 其中： : 热传导率 () : 温度梯度 () 假设单位体积的发热率为，热传导率在空间里一定并且各向相同，那么通过式(1.2.1)和(1.2.2)就可以得出如下的热传递方程式。 (1.2.3 上式根据有限元计算方法进行空间离散之后会变成时间的矩阵微分方程式，其分析结果就是各时刻的节点温度。 (1.2.4 其中： : 热容量(capacitance)矩阵 : 热传导(conduction)矩阵 : 热荷载向量 热容量矩阵是根据密度和比热计算的各个单元的矩阵，midas FEA里使用的是分布矩阵 (consistent matrix)；热传导矩阵是由热传导率计算的各个单元的矩阵；热荷载向量是根据热源、热通量、热流量荷载计算的，主要反映对流边界条件的影响。 midas FEA里可以计算稳态(steady state)温度响应和瞬态(transient)温度响应。稳态分析可以通过忽略(1.2.4)公式中热容量矩阵的效果简单的计算得出： (1.2.5 瞬态的解是通过时间的数值积分来得出的；假设时刻的温度状态()和时刻的温度状态()分别满足式(1.2.4)就可表示如下： (1.2.6 式(1.2.6)若用和加权平均则得 (1.2.7 将式(1.2.7)应用如下关系，可得和的方程式如下： (1.2.8 (1.2.9 若用如下所示的等价矩阵替换则可简化为： (1.2.10 式(1.2.10)根据使用的积分参数()的不同， = 0.0 : 向前差分(随的大小收敛) = 0.5 : Crank-Nicolson (无条件收敛) = 2/3 : Galerkin法(无条件收敛) = 1.0 : 向后差分(无条件收敛) midas FEA里除了以外可以使用所有的值;基于的一次微分方程式的非稳态分析结果的特征可以通过与如下所示的简化的一维扩散方程式的结果比较来确认： 一维扩散方程式 :  边界条件 : (), (), 初始条件 : 下面的图表是三种具有无条件收敛特性的原理在不同时间步长下的分析结果的比较，可简单整理如下： - Crank-Niconsol方法 (= 0.5)时间步长大的时候虽然振动会逐渐消失，但是刚开始计算的时候会有比较大的振动现象，减小时间步长时会提高结果的精确度。 - Galerkin法 (= 2/3)时间步长大的时候虽然计算初期也会有振动现象，但是相对Crank-Nicolson法振动会弱一些，计算初期的结果与时间步长关联不大而且非常精确。 - 向后差分法 (= 1.0)虽然没有振