有限元课件12.pptVIP

热传递问题的有限元法 CONTENTS 场问题 用加权残量法形成有限元法 1D 热传递问题 2D 热传递问题 总结 实例研究 场问题 2-D线性稳态场问题系统方程的一般形式： 场问题 2D叶片中的热传递 场问题 2D长物体的热传递 场问题 1D叶片中的热传递 场问题 复合壁的热传递 场问题 杆的扭转变形 场问题 声学问题 用加权残量法形成有限元法 不需已知泛函而通过控制方程建立有限元方程。 用加权残量法形成有限元法 离散成小单元以保证得到良好的近似 在每个单元中， 使用 N 作为权函数 1D 热传递问题 1D叶片 1D叶片 1D叶片 1D叶片 1D叶片 1D叶片 1D叶片 1D叶片 1D叶片 1D叶片 直接组装方法 1D叶片 直接组装方法 1D叶片 直接组装方法 1D叶片 算例：1D叶片中的热传递 1D叶片 1D叶片 1D叶片 复合壁 复合壁 算例：穿过复合壁的热传递 复合壁 复合壁 复合壁 复合壁 算例：通过薄膜层的热传递 复合壁 复合壁 复合壁 2D 热传递问题 单元方程 单元方程 单元方程 单元方程 单元方程 三角形单元 三角形单元 三角形单元 矩形单元 矩形单元 矩形单元 边界条件和向量 b(e) 边界条件和向量 b(e) 边界条件和向量 b(e) 边界条件和向量 b(e) 边界条件和向量 b(e) 边界条件和向量 b(e) 边界条件和向量 b(e) 边界条件和向量 b(e) 边界条件和向量 b(e) 边界条件和向量 b(e) 边界条件和向量 b(e) 电热元或热汇 位于单元内部的点热源或热汇 总结 实例研究 实例研究 实例研究 单元1： 单元2： 组装上述两单元， (未知量，用于使方程得到平衡) 解得： 原材料 0.2 mm h =0.01 W/cm 2 / ° C f f = 150 ° C 0.2mm 加热器 k =0.1 W/cm/ ° C 2mm 玻璃 铁 铂 k =0.5 W/cm/ ° ° C k =0.4 W/cm/ ° ° C 基体 等离子区 · · · · 1 2 3 4 (1) (2) (3) 300?C · · · · 1 2 3 4 (1) (2) (3) 单元1： 单元2： 单元3： 由于 ?1 = 300°C, 解得： 单元方程 对于一个单元， 注意： W = N : 伽辽金方法 (在计算该残量的积分时需利用高斯散度定理。) 故， 高斯散度定理： (乘积的微分法则) 被积函数中的第2项： 故， 式中 定义 , (应变矩阵) ? 注意：常数矩阵 (或 Ni = Li) (面积坐标) 例如 故， 类似地， 注意： b(e) 将在后面讨论 注意：实际中通常利用高斯积分法计算该积分 域内 边界上 bI(e)为零 需要在边界上计算bB(e) 不需计算 需计算 bB(e) 位于自然边界 ?2 穿过边界的热流量 绝热边界： M = S = 0 ? 对流边界条件： 给定边界上的热流量： 正数 若热从边界流入 负数 若热从边界流出 0 绝热 对于其它M 和S ? 0的情况： 式中 , 对于矩形单元： (节点1和节点2平分) 平均分配原则对具有线性形函数的所有单元均成立 同样适用于三角形单元 对于矩形单元为 分配比率为 2/6, 1/6, 1/6, 2/6 对于三角形单元存在类似关系 在热源或热汇处最好设置节点 (德耳塔函数) 点热源/热汇 ? * * 第12章 (亥姆霍兹方程) 1D问题： 注： 热传导 热对流 热源 注： Dx = kx, Dy =ky, g = 0 and Q = q 热传导 热源 注: 热传导 热对流 热源 注: 热传导 热源 热对流 热对流 Note: Dx=1/G, Dy=1/G, g=0, Q=2q 理想无漩流体的流动 注: Dx = Dy = 1, g = Q = 0 (? - 应力函数) (? - 流线函数和 ? - 势函数) 注: , Dx = Dy = 1, Q = 0 P – 外界压力； w – 波动频率； c – 介质中的波速 近似解： 权函数 (强式) (积分形式) 其中 伽辽金方法 然后组装所有单元的残量并令其等于零。 1D 叶片 (给定的边界条