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

4.1 有限差分法 4.1.2 有限差分法的主要步骤 4.1.3 差分方程的建立 差分方程的建立过程（之一） 合理选择网格布局及步长（续） 差分方程的建立过程（之二） 差 分 差 商 4.2 有限单元法 4.2.2 有限元法的基本思想 实例1 实例2 4.2.3 有限元法的基本概念与一般步骤 直接刚度法（例） 1 前处理过程 （2）直接刚度法分析（结点分析法） （3）有限单元法分析（单元分析法） 2 求解阶段 3 后处理阶段 4.2.4 有限元程序的结构和特点 有限元程序的结构和特点 （续1） 有限元程序的结构和特点 （续2） 4.2.5 有限元软件简介 ANSYS 简介 ANSYS的主要分析功能 1．结构静力分析 用来求解稳态外载引起的系统或部件的位移、应变、应力和力。静力分析很适合求解惯性和阻尼对结构的影响并不显著的问题，如确定结构中应力集中现象。ANSYS程序中的静力分析不仅可以进行线性分析，而且可以进行非线性分析，如塑性、蠕变、膨胀、大 变形、大应变及接触分析。 2．结构动力学分析 结构动力学分析用来求解随时间变化的载荷对结构或部件的影响。与静力分析不同， 动力学分析要考虑随时间变化的力载荷以及对阻尼和惯性的影响。ANSYS可进行的结构动力学分析类型包括：瞬态动力学分析、模态分析、谐波响应分析及随机振动响应分析。 W Y * * */38 X.Z.Lin * * */45 X.Z.Lin 有限差分方法是数值计算中应用非常广泛的一种方法，是求解微分方程的主要方法之一。其实质就是以有限差分代替无限微分、以差分代数方程代替微分方程、以数值计算代替数学推导的过程，从而将连续函数离散化。以有限的、离散的数值代替连续的函数分布。 4.1.1 概述 1、构成差分格式 差分方程通常是一组数量较多的线性代数方程（即：线性方程组）。其求解方法有下列两种：（1）精确法，又称直接法，即消元法；（2）近似法，又称间接法，即迭代法。 2、求解差分方程 3、对所得到的数值解进行精度与收敛性分析和检验。 首先选择网格布局、差分形式和步长；其次，以有限差分代替无限微分，即以 代替dx．以差商 代替微商（导数） ，以差分方程代替微分方程及边界条件。 建立差分方程是有限差分法的关键环节。导出差分方程的途径可有两种： （1）从微分方程出发，以泰勒级数截断，从有限差分的数学含义去建立有限差分和差分方程。 （2）从由网格所划分的单元体的能量平衡分析出发、由积分方 程去建立差分方程，该方法又称单元体平衡法。 两种方法各具特色，但无论采取何种差分方程的推导方法，在建立差分方程前，均需对所论区域进行离散化 。 ——合理选择网格布局及步长 在实施有限差分法中．首先在如图4.1所示的求解区域内，将自变量x，y 分别沿x，y轴方向的连续变化，离散为x0，x1 ，x2，…，xn及y0，y1 ，y2，…，yn个不连续点．形成离散化网格；网格交点称为结点(或节点)，依次将结点编号，与区域自变量离散化相对应，区域内函数也将同时被离散化 。 离散化后各相邻离散点之间的距离，或离散化单元的长度称为步长，步长的大小可以是常量，也可以是变量。网格的粗细与是否均匀，要根据求解区域物理场的实际分布和对结果所要求的精确度而定。 一般说来，对均质、形状简单且规则、物理量变化不剧烈的物体．或求解精度要求不高时，可采用等步长、大步长，即采用均匀网格；而对形状复杂、组分不同、物理量变化剧烈的物体，或求解精度要求较高时，则采用小步长、变步长。 图4.1 求解区域离散化 另外，对一些较复杂的问题，在选择网格与步长前，往往要对所论区域的物理场作出粗略估计，然后以较粗的网格、较大的步长计算出参考性物理场，根据这一参考性物理场再选择合理的离散化网格。 离散化网格的布局，要根据所要求解的问题的性质及求解要求确定。一般说来，有两种方法： （1）物理划分法：这种方法是根据问题的物理特性划分，如建筑物墙壁内外层面砖、普通砖和内灰泥层组成；若拟求各层界面壁温，则离散化时应按不同材料组分划分区域。 图4.2 扇形网格和三角形网格 （2）几何划分法：以几何区域形状为依据来划分，如对矩形区域可采用矩形离散化网格，非矩形区域可采用三角形、四角形或其他形状的网格，以适应温度场分布的要求。 ——将微分方程转化为差分方程 微