微分方程的数值解法.PPTVIP

加权残值法思想（续） （3） 将（3）式代入（1）和（2）式之后，一般不会满足，于是分别出现了内部和边界残差： 为了消除残差，通常引进内部权函数 和边界权函数 ，将它们分别与 和 相乘，列出消除内部残值方程式及消除边界方程式分别如下： 无网格方法思想 无网格方法（Mesh-less method）是在数值计算中不需要生成网格，而是按照一些任意分布的坐标点构造插值函数离散控制方程，问题域由一系列任意分布的节点来代替,不需要用单元或网格来进行场变量插值,也无须描述节点之间的关系。节点的生成可完全由计算机自动完成,这大大节省了分析人员的时间,也相对较容易在分析过程中对节点进行重新划分。 有限元方法（优点） 有限元可以运用于任何场问题： 没有几何形状的限制 边界条件和载荷没有限制 材料性质并不限于各向同性 具有不同行为和不同数学描述的分量可以结合。 有限元结构和被分析的物体或区域很类似 通过网格细分可以很容易地改善解的逼近度。 有限元方法 定义 应用领域 发展历史 一般步骤 计算机程序 预期目标 有限元方法（定义） 有限单元法（或有限元分析）是以剖分插值和能量原理为基础、以计算机为工具的结构分析数值方法。 基本思想：把一个大的结构划分为有限个称为单元的小区域，在每一个小区域里，假定结构的变形和应力都是简单的，小区域内的变形和应力都容易通过计算机求解出来，进而可以获得整个结构的变形和应力。 有限元方法（应用领域） 机械/航空航天/土木工程/自动化工程 结构分析（静力/动力分析，线形/非线性分析） 热分析/流体力学分析 电磁场分析 地质力学分析 生物医学分析 有限元方法（发展史） “有限单元”的名字是Clough在1960年提出的，其后不久，许多用来作应力分析的新型单元都开发出来，在1963年，有限元被认为是经典近似技术——瑞雷-里兹法的一种形式。1965年，出现了利用有限元方法分析热传导和渗流的论文。20世纪60年代末、70年代初，出现了用途广泛的有限分析计算机程序。自20世纪70年代末以来，不断增强功能的计算机图形学被用于有限元程序中，使有限元分析在实际设计和工程分析中更具吸引力。 有限元方法（分析步骤） 第一步：问题及求解域定义 第二步：求解域离散化 第三步：确定状态变量及控制方法 第四步：单元推导 第五步：总装求解。 第六步：联立方程组求解和结果解释 　　简言之，有限元分析可分成三个阶段：前处理、处理和后处理。前处理是建立有限元模型，完成单元网格划分;后处理则是采集处理分析结果，使用户能简便提取信息，了解计算结果。 有限元方法（计算机程序） 专业化小程序有 Truss Frame Plane Stress Heat Transfer 大型通用商业化程序 ANSYS ADINA ABAQUS DEFORM PANTRAN SAP 有限元方法（学习目标） 理解有限元方法的基本思想 认识不同类型单元的行为和应用范围 根据实际问题，能够给出合适的有限元模型 掌握一两中有限元分析软件的使用方法 能够解释并正确评估结果的合理性 * 假设： 注：bc_t表示bc边与厚度t组成的面，其他记号类似。 如果不区别τxy 和τyx ，则方程1-10可以写成下列形式： 如果不区别τxy 和τyx ，则方程1-10可以写成下列形式： 式1-19也可以写成下列形式： （1-20） （1-25） （1-25） 微分方程解题思路及主要解法 一阶方程 高阶方程 分离变量法 全微分方程 常数变易法 特征方程法 待定系数法 非全微分方程 非变量可分离 幂级数解法 降阶 作变换 作变换 积分因子 数值解法 微分方程的幂级数解法（概述） 例如求一阶微分方程： 满足初始条件 的特解，其中函数 这时我们可以将特解可展开为 是 和 的多项式： 的幂级数 (1) (2) 其中 是待定的系数 把（2）代入（1）中，便得一恒等式，比较这恒等式两端 的 同次幂的系数，就可定出常数 , 以这些常数为系数的级数（2）在其收敛区间内就是方程（1）满足初始条件 的特解。 微分方程的数值解法 有限差分法 有限元方法 边界元方法 加权残值方法 无网格法 在微分方程的求解中，除了采用级数和逐步逼近等方法得到解的近似表达式外，通常还有一类近似方法称为数值方法，它可以给出解在一些离散点上的近似值，这类方法通