有限元分析热分析

1、主要讲授三方面内容：主要讲授三方面内容： ANSYS热分析基础知识简介热分析基础知识简介 稳态热分析实例稳态热分析实例 瞬态热分析实例瞬态热分析实例 7.1.1 ANSYS热分析基础知识简介热分析基础知识简介 一、一、ANSYS热分析功能介绍热分析功能介绍 ANSYS热分析模块主要有：热分析模块主要有： ANSYS/Multiphysics ANSYS/Mechanical ANSYS/Thermal ANSYS/FLOTRAN ANSYS/ED 其中，其中，ANSYS/FLOTRAN不含相变热分析。不含相变热分析。 ANSYS热分析基于能量守恒原理的热平衡方热分析基于能量守恒原理的热平衡方

2、程，用有限元法计算物体内部各节点的温度，并程，用有限元法计算物体内部各节点的温度，并 导出其它热物理参数。导出其它热物理参数。 运用运用ANSYS软件可进行热传导、热对流、热软件可进行热传导、热对流、热 辐射、相变、热应力以及接触热阻等问题的分析辐射、相变、热应力以及接触热阻等问题的分析 求解求解 。 此外，此外，ANSYS不仅能解决纯粹的热分析问不仅能解决纯粹的热分析问 题，还能解决与热相关的其它问题，如热题，还能解决与热相关的其它问题，如热应应 力分析、热力分析、热电分析、热电分析、热磁分析等。一般称磁分析等。一般称 这类涉及两个或多个物理场相互作用的问题为这类涉及两个或多个物理场相互作用

3、的问题为 耦合场分析。耦合场分析。 ANSYS提供了两种分析耦合场的方法：直提供了两种分析耦合场的方法：直 接耦合法与间接耦合法。接耦合法与间接耦合法。 二、单位制问题：在二、单位制问题：在ANSYS热分析过程中，不一热分析过程中，不一 定都要采用国际单位制，但必须要使所有物理定都要采用国际单位制，但必须要使所有物理 量的单位统一起来。量的单位统一起来。 ANSYS中共有五种单位可供选择（命令流中共有五种单位可供选择（命令流 方式：方式：/UNITS；或；或Main menuPreprocessorMaterial PropsMaterial Library Select Units）：）：

4、SI（MKS）代表国际单位制，其基本单位）代表国际单位制，其基本单位 为为m，kg，s，K。 CGS代表厘米、克、秒单位制，其基本单代表厘米、克、秒单位制，其基本单 位为位为cm，g，s，。 BFT代表以英尺为主的英制单位制，其基代表以英尺为主的英制单位制，其基 本单位为本单位为ft，slug，s， 。 BIN代表以英寸为主的英制单位制，其基本代表以英寸为主的英制单位制，其基本 单位为单位为in，ibm，s， 。 USER代表用户自定义单位制，即用户可以代表用户自定义单位制，即用户可以 根据需要定义基本单位。根据需要定义基本单位。 三、热分析时的三类边界条件和初始条件：三、热分析时的三类边界条

5、件和初始条件： 第一类边界条件：物体边界上的温度函数已知；第一类边界条件：物体边界上的温度函数已知； 第二类边界条件：物体边界上的热流密度已知；第二类边界条件：物体边界上的热流密度已知； 第三类边界条件：与物体相接触的流体介质的温第三类边界条件：与物体相接触的流体介质的温 度和换热系数已知。度和换热系数已知。 初始条件：初始条件是指传热过程开始时，物体初始条件：初始条件是指传热过程开始时，物体 在整个区域中所具有的温度为已知值。在整个区域中所具有的温度为已知值。 四、热分析时的载荷：四、热分析时的载荷：ANSYS共提供了共提供了6种载荷，种载荷， 可以施加在实体模型或单元模型上。可以施加在实体

6、模型或单元模型上。 （1）温度：作为第一类边界条件，温度可以施加）温度：作为第一类边界条件，温度可以施加 在有限元模型的节点上，也可以施加在实体模在有限元模型的节点上，也可以施加在实体模 型的关键点、线段及面上。型的关键点、线段及面上。 （2）热流率：热流率（）热流率：热流率（Heat Flow）是一种）是一种节点节点 集中载荷集中载荷，只能施加在节点或关键点上，主要，只能施加在节点或关键点上，主要 用于线单元模型。用于线单元模型。 （3）对流：对流（）对流：对流（Convection）是一种）是一种面载荷面载荷， 用于计算流体与实体的热交换。它可以施加在用于计算流体与实体的热交换。它可以施加

7、在 有限元模型的节点及单元上，也可以施加在实有限元模型的节点及单元上，也可以施加在实 体模型的线段和面上。体模型的线段和面上。 （4）热流密度：又称热通量（）热流密度：又称热通量（Heat Flux）单位为）单位为 W/m2。热流密度是一种。热流密度是一种面载荷面载荷，表示通过单，表示通过单 位面积的热流率。当通过单位面积的热流率已位面积的热流率。当通过单位面积的热流率已 知时，可在模型相应的外表面施加热流密度。知时，可在模型相应的外表面施加热流密度。 若输入值为正，则表示热流流入单元；反之，若输入值为正，则表示热流流入单元；反之， 则表示热流流出单元。它可以施加在有限元模则表示热流流出单元。

8、它可以施加在有限元模 型的节点及单元上，也可以施加在实体模型的型的节点及单元上，也可以施加在实体模型的 线段和面上。线段和面上。 （5）生热率：生热率既可看成是材料的一种基本）生热率：生热率既可看成是材料的一种基本 属性，又可作为载荷施加在单元上。它可以施属性，又可作为载荷施加在单元上。它可以施 加在有限元模型的节点及单元上，也可以施加加在有限元模型的节点及单元上，也可以施加 在实体模型的关键点、线段、面及体上。在实体模型的关键点、线段、面及体上。 （6）热辐射率：热辐射率也是一种）热辐射率：热辐射率也是一种面载荷面载荷，通，通 常施加于实体的外表面。它可以施加在有限元常施加于实体的外表面。它

9、可以施加在有限元 模型的节点及单元上，也可以施加在实体模型模型的节点及单元上，也可以施加在实体模型 的线段和面上。的线段和面上。 五、热分析时的三种传热方式及材料基本属性五、热分析时的三种传热方式及材料基本属性 （1）热传导：当物体内部存在温差，即存在温度）热传导：当物体内部存在温差，即存在温度 梯度时，热量从物体的高温部分传递到低温部梯度时，热量从物体的高温部分传递到低温部 分；而且不同温度的物体相互接触时热量会从分；而且不同温度的物体相互接触时热量会从 高温物体传递到低温物体。这种热量传递的方高温物体传递到低温物体。这种热量传递的方 式称为热传导。式称为热传导。 Q/t=KA(Thot-T

10、cold)/d 式中：式中：Q为时间为时间t内的传热量或热流量；内的传热量或热流量；K为热为热 传导率或热传导系数；传导率或热传导系数； （2）对流：热对流是指固体的表面与它周围接触）对流：热对流是指固体的表面与它周围接触 的液体或气体（统称为流体）之间，由于温差的液体或气体（统称为流体）之间，由于温差 的存在而引起的热量交换。的存在而引起的热量交换。 高温物体表面（如暖气片）常常发生对流高温物体表面（如暖气片）常常发生对流 现象，这是因为高温物体表面附近的空气因受现象，这是因为高温物体表面附近的空气因受 热而膨胀，密度降低并向上流动。与此同时，热而膨胀，密度降低并向上流动。与此同时， 密度较

11、大的冷空气下降并代替原来的受热空气。密度较大的冷空气下降并代替原来的受热空气。 热对流可以分为两类：自然对流和强制对流。热对流可以分为两类：自然对流和强制对流。 q=h(TS-TB) 式中：式中： h为对流换热系数（或称膜传热系数、给热系数、为对流换热系数（或称膜传热系数、给热系数、 膜系数等）；膜系数等）； TS为固体表面的温度；为固体表面的温度； TB为周围流体的温度。为周围流体的温度。 （3）辐射：热辐射是指物体发射电磁能，并被其）辐射：热辐射是指物体发射电磁能，并被其 它物体吸收转变为热能的热量交换过程。物体它物体吸收转变为热能的热量交换过程。物体 温度越高，单位时间辐射的热量越多。热

12、传导温度越高，单位时间辐射的热量越多。热传导 和热对流都需要有传热介质，而热辐射无须任和热对流都需要有传热介质，而热辐射无须任 何介质。何介质。 实质上，在真空中的热辐射效率最高。在实质上，在真空中的热辐射效率最高。在 工程中通常考虑两个或两个以上物体之间的辐工程中通常考虑两个或两个以上物体之间的辐 射，系统中每个物体同时辐射并吸收热量，它射，系统中每个物体同时辐射并吸收热量，它 们之间的净热量传递可以用如下斯蒂芬们之间的净热量传递可以用如下斯蒂芬波尔波尔 兹曼方程来计算：兹曼方程来计算： q=A1F12(T41-T42) 式中：式中： q为热流率；为热流率； 为实际物体的辐射率，或称黑度，它

13、的数值为实际物体的辐射率，或称黑度，它的数值 处于处于01之间；之间； 为斯蒂芬为斯蒂芬波尔兹曼常数，约为波尔兹曼常数，约为 5.67108W/m2.K4； A1为辐射面为辐射面1的面积；的面积； F12为由辐射面为由辐射面1到辐射面到辐射面2的形状系数；的形状系数； T1为辐射面对为辐射面对1的绝对温度；的绝对温度； T2为辐射面为辐射面2的绝对温度。的绝对温度。 在在ANSYS热分析中，在确定分析选项，即热分析中，在确定分析选项，即 Main MenuSolutionAnalysis TypeAnalysis Options的对话框中有一个选项：的对话框中有一个选项：Temperature

14、 difference，该选项用于确定绝对零度，即需要，该选项用于确定绝对零度，即需要 将目前的温度值换算为绝对温度。如果在热分将目前的温度值换算为绝对温度。如果在热分 析过程中使用的温度单位是摄氏度，则该值应析过程中使用的温度单位是摄氏度，则该值应 设定为设定为273。 从上式可以看出，包含热辐射的热分析是从上式可以看出，包含热辐射的热分析是 高度非线性的。高度非线性的。 （4）比热容（）比热容（Specific Heat）：是指单位质量的）：是指单位质量的 物质每升高（或降低）物质每升高（或降低）1所吸收（或放出）的所吸收（或放出）的 热量，简称比热，其单位为热量，简称比热，其单位为J/(

15、Kg.)。其计算。其计算 公式为：公式为： C=Q/（m.T） 式中：式中：T= TE-TB，为，为TE为终止时刻温度；为终止时刻温度；TB 为开始时刻温度；为开始时刻温度；Q为该时间段内物体吸收或为该时间段内物体吸收或 放出的总热量；放出的总热量；m为质量。为质量。 （5）焓（）焓（Enthalpy）：）： 定义为定义为H=U+PV，式中，式中H为焓，为焓，U为内能，为内能，P、 V分别为压力和体积。分别为压力和体积。 对于常压情况，上式又可表示为：对于常压情况，上式又可表示为： H=U+P.V=Q 说明在常压条件下，焓的变化即为热量的说明在常压条件下，焓的变化即为热量的 变化。变化。 （6

16、）生热率（）生热率（Heat Generation Rate）：）： 生热率既可以材料属性的形式进行定义，生热率既可以材料属性的形式进行定义， 同时又可以体载荷的形式施加到单元上，用于同时又可以体载荷的形式施加到单元上，用于 模拟化学反应生热或电流生热，其单位是单位模拟化学反应生热或电流生热，其单位是单位 体积的热流率。体积的热流率。 六、稳态热分析与瞬态热分析：六、稳态热分析与瞬态热分析： （1）稳态传热：）稳态传热： 如果系统的净热流率为如果系统的净热流率为0，即流入系统的热，即流入系统的热 量加上系统自身产生的热量等于流出系统的热量加上系统自身产生的热量等于流出系统的热 量，则系统处于热

17、稳态。在稳态热分析中任一量，则系统处于热稳态。在稳态热分析中任一 节点的温度都不随时间变化。节点的温度都不随时间变化。 （2）瞬态传热：）瞬态传热： 瞬态传热过程是指一个系统的加热或冷却过瞬态传热过程是指一个系统的加热或冷却过 程。在这个过程中，系统的温度、热流率、热程。在这个过程中，系统的温度、热流率、热 边界条件以及系统内能随时间都有明显变化。边界条件以及系统内能随时间都有明显变化。 七、线性与非线性热分析七、线性与非线性热分析 ANSYS在热分析过程中，如果有下列情况中在热分析过程中，如果有下列情况中 的一种或几种出现，则该分析为非线性热分析：的一种或几种出现，则该分析为非线性热分析：

18、材料热性能随温度变化；材料热性能随温度变化； 边界条件随温度变化；边界条件随温度变化； 含有非线性单元；含有非线性单元； 考虑辐射传热。考虑辐射传热。 7.1.2 稳态热分析实例稳态热分析实例1长空心圆柱长空心圆柱 体的热传导过程体的热传导过程ANSYS分析分析 问题描述：问题描述： 有一空心钢圆柱体，内半径与外半径分别为有一空心钢圆柱体，内半径与外半径分别为 0.2 m、0.6 m，长度为，长度为10 m。钢的导热系数为。钢的导热系数为 70W/（m.)，现在柱体的外表面施加均匀温，现在柱体的外表面施加均匀温 度载荷度载荷80，假设柱体内表面温度为恒定值，假设柱体内表面温度为恒定值 20。试

19、求钢柱体内部的温度场分布。试求钢柱体内部的温度场分布。 简要分析：简要分析： 该问题属于稳态热力学问题。由于柱体的长度远大该问题属于稳态热力学问题。由于柱体的长度远大 于其直径，可忽略其终端效应，同时根据问题的对称性，于其直径，可忽略其终端效应，同时根据问题的对称性， 求解过程中取圆柱体横截面的求解过程中取圆柱体横截面的1/4建立几何模型。建立几何模型。 几何模型 1. 建立工作文件名和工作标题，并选择建立工作文件名和工作标题，并选择Preferences进行筛选：进行筛选： 2. 定义单元类型和材料属性定义单元类型和材料属性 3. 创建几何模型、划分网格创建几何模型、划分网格 3.1 几何建

20、模几何建模 3.2 划分网格：先对线进行标注，然后画线以便于操作。划分网格：先对线进行标注，然后画线以便于操作。 4 加载求解加载求解 4.1 选择分析类型：选择分析类型： 4.2 对线上各节点施加温度载荷：先对对线上各节点施加温度载荷：先对1线上的节点加温线上的节点加温 度载荷度载荷 4.3 定义载荷定义载荷 4.4 定义定义2线上的温度载荷：线上的温度载荷： 4.5 定义温度载荷值：定义温度载荷值： 4.6 载荷步及输出结果控制选项设置：载荷步及输出结果控制选项设置： 4.7 求解：求解： 5. 查看求解结果：查看求解结果： 7.1.3 瞬态热分析实例瞬态热分析实例1型材瞬态传型材瞬态传

21、热过程热过程ANSYS分析分析 问题描述：问题描述： 有一横截面为矩形的各向异性型材，其初有一横截面为矩形的各向异性型材，其初 始温度为始温度为500，现突然将其置于温度为，现突然将其置于温度为20的的 空气中，求空气中，求1分钟后该型材的温度场分布及其中分钟后该型材的温度场分布及其中 心温度随时间的变化规律。材料性能参数如下：心温度随时间的变化规律。材料性能参数如下： 密度为密度为2400 kg/m3，导热系数，导热系数KXX为为30 W/(m.)，KYY、KZZ为弹性模量为为弹性模量为10 W/(m.)，比热为，比热为352 J/(kg.)，对流系数为，对流系数为 110W/(m2.) 。

22、 简要分析：简要分析： 该问题属于瞬态热传导问题。由于材料沿长度方向该问题属于瞬态热传导问题。由于材料沿长度方向 的尺寸远大于其它两个方向的尺寸，将其简化为平面应的尺寸远大于其它两个方向的尺寸，将其简化为平面应 变问题。在分析过程中取型材横截面的变问题。在分析过程中取型材横截面的1/4建立模。建立模。 0.4 0.2 1. 建立工作文件名和工作标题，并选择建立工作文件名和工作标题，并选择Preferences进行筛选：进行筛选： 2. 定义单元类型和材料属性：定义单元类型和材料属性： 3. 创建几何模型、划分网格创建几何模型、划分网格 3.1 几何建模：几何建模： 对线进行标注并画线，结果如下

23、图所示对线进行标注并画线，结果如下图所示 ： 3.2 划分网格：先设置单元尺寸（对所选划分网格：先设置单元尺寸（对所选1号线进行定义）：号线进行定义）： 结果如下图所示：结果如下图所示： 对对3号线进行相似操作：号线进行相似操作： 对对2号线进行相似操作：号线进行相似操作： 对对4号线进行相似操作（号线进行相似操作（NDIV取取30，SPACE取取5）后得）后得 ： 3.3 用用Mesh Tool划分网格：划分网格： 结果如下图所示结果如下图所示 ： Change Title：ELEMENTS IN MODEL后点击后点击Plot Elements为为 ： 4. 加载求解加载求解 4.1 分析

24、类型选择：分析类型选择： 4.2 求解控制选择：求解控制选择：Basic选项（变化）和选项（变化）和Transient选项（没变）的选项（没变）的 设置设置 ： 4.3 施加温度载荷：先对全部施加均布温度载荷施加温度载荷：先对全部施加均布温度载荷500 再对线再对线2和线和线3上的各节点施加温度及对流载荷。首先点击上的各节点施加温度及对流载荷。首先点击Select Entities：如下图（选定所要加对流和温度的线：通过拾取获得）：如下图（选定所要加对流和温度的线：通过拾取获得）： 点击所要加对流及温度的两条线：点击所要加对流及温度的两条线： 然后选择从属于此线上的所有节点，按下图操作：然后选

25、择从属于此线上的所有节点，按下图操作： 对所选节点施加对流及温度载荷：按下图选择，然后点击对所选节点施加对流及温度载荷：按下图选择，然后点击Pick All： 对流及温度值的设置对流及温度值的设置 ： Select everything后得后得 ： 分析类型选项的设置分析类型选项的设置 ： 3.4 求解求解 4.查看结果查看结果 4.1 先点击先点击Read Results中的中的Last Set 4.2 然后然后Plot Results，按下图设置：，按下图设置： 最后，最后，1分钟后型材横截面上温度场等值线图分布如下分钟后型材横截面上温度场等值线图分布如下 ： 4.3 察看中心温度云纹图。坐标轴设置：察看中心温度云纹图。坐标轴设置： PlotCtrlsStyleGraphsModify Axes ： 曲线设置：曲线设置：PlotCtrlsStyleGraphsModify Curves ： 颜色设置：颜色设置：PlotCtrlsStyleColorsGraph Colors：在：在CURVE Graph curve number 1下拉列表中选择黄色：下拉列表中选择黄色： 4.4 定义时域后处理变量：定义时域后处理变量：Main