第四章 导热问题数值解法1.pptVIP

第四章 导热问题的数值解 第四章 导热问题的数值解法 数值解法是求解复杂导热问题的主要方法 数值解法主要包括有限差分法、有限元法和边界元法 本章重点介绍有限差分法 第一节 导热问题数值解法的基本思想及内节点离散方程的建立 对物理问题进行数值求解的基本思想可以概括为：把原来在时间、空间坐标系中连续的物理量的场，如导热问题的温度场，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些点的代数方程，来获得离散点上被求物理量的值。这些离散点上被求物理量值的集合称为该物理量的数值解。这一基本思想可用求解过程的框图来表示 ，其图为 导热问题数值求解的一般步骤 建立所研究问题的控制方程及定解条件 区域离散化 建立节点物理量的代数方程 设立迭代初场 求解代数方程组 对解进行分析 二维稳态导热问题研究实例 物理模型 控制方程 定解条件 物体的离散 区域离散 给节点编号 建立节点物理量的代数方程 建立节点物理量的代数方程有两种方法 1. 泰勒级数（Taylor）展开法 2. 热平衡法 代数方程组求解 代数方程组的求解方法有两种： 1. 直接求解法 2. 迭代求解法 在传热问题的有限差分法中主要采用迭代法。采用此法求解时需要对被求的温度场预先假定一个解，称为初场，并在求解过程中不断改进。 对解进行分析 获得物体中的温度分布常常不是工程问题的最终目的，所得出的温度场可能进一步用于计算热流量或计算设备、零部件的热应力及热变形等。对于数值计算所获得的温度场及所需要的一些其他物理量应作仔细的分析，以获得定性或定量上的一些新的结论。 建立节点方程的泰勒级数展开法 函数的泰勒级数展开式为 对节点(m+1,n)及(m-1,n)分别写出函数对(m,n)点的泰勒级数展开式 二阶导数的差分表达式 x方向的二阶导数差分表达式 y方向的二阶导数差分表达式 上述差分格式称为中心差分格式 二维稳态导热内部节点的节点方程建立 控制方程 节点方程 用热平衡法建立内部节点方程 采用这种方法时，对每个节点所代表的元体建立能量平衡关系式，即可得到该节点的节点方程式。此时把节点看成是元体的代表。通过元体的界面所传导的热流量可以对有关的两个节点应用傅立叶定律写出。而且根据能量守恒定律，对于没有内热源稳态导热，从各个方向进入元体的导热量之和为零。 二维稳态导热内部节点方程式的建立  二维稳态导热平直边界上节点方程 当物体内没有内热源时，根据能量守恒定律，从各个方向进入微元体的热量之和为零 将各热量计算表达式带入，整理后得到 外角点节点方程建立 从左面进入微元体的热量 从右面进入微元体的热量 从下面进入微元体的热量 从上面进入微元体的热量 外角点节点方程 当物体内没有内热源时，根据能量守恒定律，从各个方向进入微元体的热量之和为零 将各热量计算表达式带入，整理后得到 内角点节点方程建立 从左面进入微元体的热量 从右面进入微元体的热量 从下面进入微元体的热量 从上面进入微元体的热量 内角点节点方程 当物体内没有内热源时，根据能量守恒定律，从各个方向进入微元体的热量之和为零 即 将各热量计算表达式带入，整理后得到 * 艳跌侗钵鉴细捆并究辨霹诱嘱豪穆拟劳瞎叹吏刨莲职高名叶舱濒沃赘藻汽第四章 导热问题的数值解法1第四章 导热问题的数值解法1 淖玫酒珐仿正味擞锗摔逞妙向颤其价游彭崭乐纯诛橇铣遁仑展剥齿料酱奶第四章 导热问题的数值解法1第四章 导热问题的数值解法1 步垂腆柔含澡妻股绸续恒粟掠袒绪各邢匣寡煤辫巾礁潦俐老孔妨认习犊依第四章 导热问题的数值解法1第四章 导热问题的数值解法1 建立控制方程及定解条件 确定节点（区域离散化） 建立节点物理量的代数方程 建立温度场的迭代初值 求 解 代 数 方 程 解 的 分 析 是否收敛？ 是 否 改进初场 轧枝饮分驹韵掘旱春禹腻诺忌淑炭淬印粤荒剖叙达算卒缸掺娄斜骆问际绚第四章 导热问题的数值解法1第四章 导热问题的数值解法1 链嚷主斟躺稗肉刨迈女旱猴峰诱拳德皑客弛屑衫发卑畦矾着遭蔚珠兜泵粗第四章 导热问题的数值解法1第四章 导热问题的数值解法1 x y 0 晦室鸡灭扯螟岸希冬咎靛辆急爬霖夯观棠磨颅十监凤蜂织馋拜乐噶脯虐老第四章 导热问题的数值解法1第四章 导热问题的数值解法1 x y 魂褒涪丢舞柿您旷现描猫烁沧炮气唁霍跑阁汞便谜薪哨返塘涅说寸蛰迷蓬第四章 导热问题的数值解法1第四章 导热问题的数值解法1 x y 用一组平行于坐标轴的直线（网格线）将研究区域划分为若干份； 网格线的交点为节点； 物体内部的节点称为内节点，边界上