用Matlab求解一维非稳态周期性导热问题（有限单元法+隐式离散+高斯赛德尔迭代法）

原创

已于 2023-08-03 16:50:49 修改 · 7.8k 阅读

于 2022-05-28 12:48:45 首次发布

一、问题描述与分析

本次问题条件如下：

计算模拟如下一维常物性无内热源非稳态导热的温度场，以及内外壁面的热流密度，并进行温度场和热流的特点分析，相关参数如下。

室内温度恒定为20℃，无限大平壁的比热为1.0kJ/kg℃，室外气温周期性变化，周期为1天（即24小时）。

二、网格划分与节点方程式

由问题所给条件可知，两侧边界条件均为对流边界条件，室内侧温度恒定，室外侧温度周期性变化，无限大平壁内部没有热源，将该无限大平壁沿厚度方向均匀的划分为300个单元，左侧室内右侧室外，节点编号依次为1-2-3-...-300-301：

则该一维常物性无内热源非稳态导热问题的有限单元节点方程如下（注这里不是有限差分法，是有限单元法，离散格式为隐式）：

左端节点方程

右端节点方程

中间节点方程

其中，i表示空间维度节点，t表示时间维度节点。上述公司中，左端是t时刻的节点温度，是未知量需要求解的，右端是t-1时刻的节点温度，是已知的。这部分的公式要注意和后面高斯赛德尔迭代进行区分，下一节的迭代次数专门用了不同字母表示，就是为了防止非稳态传热的时间递进与高斯赛德尔本身的迭代递进弄混。

三、高斯赛德尔迭代法

上述节点方程式可以通过高斯赛德尔迭代法求解。节点方程式本质是一个n+1元的线性方程组（假定划分单元数为n）：

$\begin{cases} a_{1,1}x_1 + a_{1,2}x_2 + \cdots + a_{1,n}x_{n+1} = b_1 \\ a_{2,1}x_1 + a_{2,2}x_2 + \cdots + a_{2,n}x_{n+1} = b_2 \\ \cdots \\ a_{n+1,1}x_1 + a_{n+1,2}x_2 + \cdots + a_{n+1,n+1}x_{n+1} = b_{n+1} \end{cases}$

其矩阵形式如下：

$\begin{pmatrix} a_{1,1} & a_{1,2} & \cdots & a_{1,n+1} \\ a_{2,1} & a_{2,2} & \cdots & a_{2,n+1} \\ \vdots & \vdots & \ddots & \vdots \\ a_{n+1,1} & a_{n+1,2} & \cdots & a_{n+1,n+1} \end{pmatrix} \begin{pmatrix} x_1 \\ x_2 \\ \vdots \\ x_{n+1} \end{pmatrix} = \begin{pmatrix} b_1 \\ b_2 \\ \vdots \\ b_{n+1} \end{pmatrix}$