有限元热分析一.ppt

第一部分:热分析的基本概念 符号 ANSYS中标准单位 ( 英制 ) 温度 热流量 热传导率 密度 比热 对流换热系数 热流 温度梯度 内部热生成 ANSYS中标准单位 (SI) 温度 热流量 热传导率 密度 比热 对流换热系数 热流 温度梯度 内部热生成 热传递的类型 热传递有三种基本类型: 传导 : 两个良好接触的物体之间的能量交换或一个物体内由于温度梯度引起的内部能量交换。 对流 : 在物体和周围介质之间发生的热交换。 辐射 : 一个物体或两个物体之间通过电磁波进行的能量交换。 在绝大多数情况下，我们分析的热传导问题都带有对流和/或辐射边界条件。 传导 传导的热流由传导的傅立叶定律决定： 负号表示热沿梯度的反向流动(i.e. ,热从热的部分流向冷的)。 对流 对流的热流由冷却的牛顿准则得出: 对流一般作为面边界条件施加 辐射 从平面 i 到平面 j 的辐射热流由施蒂芬-玻斯曼定律得出: 在ANSYS中将辐射按平面现象处理(i.e., 体都假设为不透明的)。 有限元方法 将区域分解或划分为简单的形状: 2-D模型中的四边形和/或三角形，3-D模型中的四面体，金字塔形或六面体。 有限元方法 ( 续 ) 假设单元内温度变化可以用多项式表示。一般情况下，根据单元类型的不同，应当包含不同的一次项, 平方和混合的立方项。多项式假设保证了温度在单元内部和单元边界上都是连续的。 写出以单元节点温度为未知数的多项式: 有限元方法 ( 续 ) 由单元节点温度得出每个单元的温度梯度和热流。 有限元方法 ( 续 ) 将假设的温度变化代入积分方程，注意到每项都乘上了实际的温度数值，将两边约去得到: 有限元方法 ( 续 ) 将方程可以重新写为简化形式: 有限元方法 ( 续 ) 其中: 有限元方法 ( 续 ) 系统方程是将单元的贡献组装而成: 例子: 3节点三角形单元 热传递的有限元方法可以用简单的3节点三角形实体单元来说明。使用4节点实体单元更好一些，但在本题中，使用形函数更加简单的线性三角形单元。 例子: 3节点三角形单元 (续) 有限元单元模型: 2 三角形单元 4 节点 推导单元 1 矩阵: 例子: 3节点三角形单元 (续) 推出梯度-温度矩阵 定义各向同性材料 特性矩阵 单元传导矩阵 例子: 3节点三角形单元 (续) 对流对传导矩阵的贡献 对流节点热流向量 例子: 3节点三角形单元 (续) 同样得到单元 2的矩阵并组合成为总体矩阵 矩阵可以分块如图，因为T3 = T4 = 0 同时求解得到未知的温度 求解单元1节点3的响应热流 例子: 3节点三角形单元 (续) 计算单元1的温度梯度向量* 计算单元1的热流向量* 有限元热分析中的基本特性 求解连续性 温度在一个单元中和单元内部边界上是连续的(i.e. ,单值的) 温度梯度和热流在一个单元中是连续的，在单元内部边界上是不连续的 能量平衡在每个节点上都能够满足，因为基本方程表示了节点能量平衡。 热传导的傅立叶定律满足，因为它用于推导基本方程并用于从单元温度梯度中求解单元热流。 有限元热分析中的基本特性(续) 一般来说，稳态分析中网格上节点温度比实际温度要低。也就是说，如果加密网格，温度将增加，但加密到一定程度，结果将不显著增加(i.e. , 结果收敛)。 引起奇异性的原因 整体求解的奇异性 在稳态分析中当有热量输入(e.g. , 施加节点热流，热流，内部热源)而无热流流出(指定的节点温度，对流载荷等)，稳态的温度将是无限大的。 等同于结构分析中的刚体位移。 温度梯度/热流奇异性 如果对点热源处的网格细分下去的话，梯度/热流将无限增加。 凹角和网格中的“裂缝”。 形状不好的单元。 网格划分误差 实际上任何产生不连续热流区域的有限元模型都是有误差的。在单元内部边界上热流不连续的大小将作为ANSYS进行误差估计的基础。 网格划分误差估计一般用于实体和壳单元，而且单元所在区域的单元类型是均一的(e.g., 具有共同的特性) ，热流在该区域中也就是连续的。 误差计算的细节在ANSYS 理论手册 ANSYS 误差度量(续) ANSYS 计算了几个数值，可以用来评估网格划分误差。误差计算可以用于线性和非线性的稳态分析，在通用后处理器- POST1中进行。 ANSYS中的网格划分误差度量功能: TERR - 估计选定单元中的热耗散能。在POST1中可以使用ETABLE命令存储，排序和列表 。TERR的云图可以使用Contour Plot Element Solution来完成。 TDSG - 单元中最大的热流偏差。计算单元中每个节点在各方向上平均热流和非平均热流之间最大差值。存储，排序，列表和绘图方法与TERR类