第二章材料科学研究中常用的数值分析方法.pptVIP

第二章材料科学研究中常用的数值分析方法 解决此类问题的通常途径： ①对方程及边界条件加以简化，得到问题在简化条件下的解答；（只在少数情况下有效，因为过多的简化会引起模型失真，导致较大的误差，甚至得到错误的结论） ②采用数值解法，得到模型的数值解。随着计算机技术的发展和应用，数值分析方法已经成为求解科学技术问题的主要工具。 目前，常用的数值分析方法大致分为两大类：有限差分法和有限元法。 第一节 线性方程组的数值解法 直接法和间接法。 直接法：精度高，重复工作量小，但编制计算程序复杂，对计算机资源占用较多。 间接法：即迭代法。计算程序简单，占用内存小，但重复工作量大，计算精度取决于迭代次数。 一、直接法： 课 后 作 业（高斯列主元消去法） 在四位十进制的限制下，分别用顺序高斯消去法和列主元消去法求解下列线性方程组。 3、追赶法 课堂作业（追赶法） 二、间接法（迭代法） 直接法比较适用于中小型方程组。对高阶方程组，既使系数矩阵是稀疏的，但在运算中很难保持稀疏性，因而有存储量大，程序复杂等不足。 迭代法则能保持矩阵的稀疏性，具有计算简单，编制程序容易的优点，并在许多情况下收敛较快。故能有效地解一些高阶方程组。 1、雅可比(Jacobi)迭代法（简单） 课 后 作 业 用Jacobi迭代法求解方程组，要求 2、高斯—塞德尔（Gauss-Seidel)迭代法  用Gause-Seidel迭代法求解方程组，要求 3、超松弛迭代法 课 后 作 业 试用ω＝1.25的超松弛迭代法求解方程组，要求 第二节 有限差分法（FDM） 有限差分方法是计算机数值模拟最早采用的方法，至今仍被广泛运用。 该方法将求解区域划分为差分网格，用有限个网格节点代替连续的求解域。 有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。 该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 差分的概念：某物理量的增量。 1、向前差分： (1)一阶差分: (2)二阶差分： 一般的m阶差分用m-1阶差分定义： 称为向前差分算子。 3、中心差分 概述 实质是以有限差分代替无限微分、以差分代数方程代替微分方程、以数值计算代替数学推导的过程，将连续函数离散化，以有限的、离散的数值代替连续的函数分布。 有限差分法的主要步骤： 差分方程的建立 在建立差分方程前，均需对所求解区域进行离散化。 1、合理选择网格布局及步长 将自变量x,y分别沿x,y轴方向离散，形成离散化网格。网格交点称为节点（或结点），依次将节点编号。 2、将微分方程转化为差分方程 实际上是以差分代替微分，以差商代替微商，是以有限小量代替无限微量的近似化过程。 有限差分法计算的误差分析 ui＝u(x) 将离散化的i+1,i-1结点的函数值ui+1,ui-1分别按Talor级数展开： 例：利用差分法求解方程 第三节 有限元法 finite?element?method 有限元分析技术是最重要的工程分析技术之一。 用有限个单元将连续体离散化，是通过对有限个单元作分片插值求解各种力学、物理问题的一种数值方法，单元之间通过有限个节点相互连接，然后根据能量方程或加权残量方程可建立有限个待定参量的代数方程组，求解此离散方程组就能得到有限元法的数值解。 在有限小的单元内设定试探函数。 基本理论要用到数学、力学方面的各种知识。 有限元法的发展离不开计算机。 一、有限元法的基本概念－直接刚度法 只能处理一些比较简单的实际问题，但对理解和明确有限元法的一些物理概念很有帮助。 例：杆的一端固定，另一端承受P＝1000N载荷，顶部宽w1=2cm，底部宽w2=1cm，厚度为t=0.125cm，长度L=10cm，弹性模量E=10.4×106MPa。分析杆沿长度分析不同位置的变形情况，杆的质量忽略不计。 将求解区域离散化，划分为5个结点和4个单元，每个单元的面积为组成该单元的结点处面积的平均值。 长度为L，具有均一截面A的固体单元在受到外力F时， 未来发展趋势——多尺度模拟 有限元（FEM）能够准确的拟合构件的边界形状，算法复杂。 有限差分（FDM）不能进行应力分析 应用领域1—电子行业 按键舒适性分析 电视机线性动力学分析 热传导和热应力问题分析 共晶焊点热疲劳失效分析 复杂几何体的冲击损坏 断裂和材料失效分析 应用领域2—水工大坝 大坝施工过程模拟 蓄水模拟 快速降水模拟 ?抗震模拟 断裂和材料失效分析 应用领域3—桥梁工程 桥梁