常用数值分析方法4有限差分法和有限单元法.pptVIP

4.1有限差分法 4.1.1概述 有限差分方法是数值计算中应用非常广泛的一种 方法,是求解微分方程的主要方法之一。其实质就是 以有限差分代替无限微分、以差分代数方程代替徼分 方程、以数值计算代替数学推导的过程,从而将连续 函数离散化。以有限的、离散的数值代替连续的函数 分布。 15:392018/11/23 以课又生院 X./ Lin  4.1.2有限差分法的主要步骤 构成差分格式 首先选择网格布局、差分形式和步长;其次,以有限差分 代替无限微分,即x2-x=A代替dx.以差当一=N代 ,-. 替微商(导数“,以差分方程代替微分力程皮齐件。 差分方程通常是一组数量较多的线性代数方程(即:线性方 程组)。其求解方法有下列两种:(1)精确法,又称直接法, 即消元法;(2)近似法,又称间接法,即迭代法 3、对所得到的数值解进行精度与收敛性分析和检验。 15:392018/11/23 以课又生院 X./ Lin  4.1.3差分方程的建立 建立差分方程是有限差分法的关键环节。导出差分 方程的途径可有两种 (1)从微分方程出发,以泰勒级数截断,从有限差分 的数学含义去建立有限差分和差分方程 (2)从由网格所划分的单元体的能量平衡分析出发、 由积分方程去建立差分方程,该方法又称单元体平衡法 两种方法各具特色,但无论采取何种差分方程的推 方法,在建立差分方程前,均需对所论区域进行离散 15:392018/11/23 以课又生院 X./ Lin  差分方程的建立过程(之一) 理选择网格布局及步 在实施有限差分法中.首先在如图41所示的求解区域内,将自 变量x,y分别沿x,y轴方向的连续变化,离散为xo,x1,x2,…, yn个不连续点.形成离散化网格;网格交点称 为结点(或节点),依次将结点编号,与区域自变量离散化相对应, 区域内函数也将同时被离散化 离散化后各相邻离散点之间的距离,或离散化单 元的长度称为步长,步长的大小可以是常量,也可以 是变量。网格的粗细与是否均匀,要根据求解区域物 理场的实际分布和对结果所要求的精确度而定。 般说来,对均质、形状简单且规则、物理量变 化不剧烈的物体.或求解精度要求不高时,可采用等 步长、大步长,即采用均匀网格;而对形状复杂、组 分不同、物理量变化剧烈的物体,或求解精度要求较 图4求解区城离散化厂高时,则采用小步长、变步长  合理选择网格布局及步长(续) 离散化网格的布局,要根据所要求解的问题的性质及求解 要求确定。一般说来,有两种方法 (1)物理划分法:这种方法是根据问题的物理特性划分,如建筑 物墙壁内外层面砖、普通砖和内灰泥层组成;若拟求各层界面壁 温,则离散化时应按不同材料组分划分区域。 另外,对一些较复杂的问题,在选择网格与步长前,往往要对所论区 域的物理场作出粗略估计,然后以较粗的网格、较大的步长计算出参考性 (2)几何划分法:以几何区域 形状为依据来划分,如对矩形区 域可采用矩形离散化网格,非矩 形区域可采用三角形、四角形或 其他形状的网格,以适应温度场 分布的要求。 图42扇形网格和三角形网格∠Mm