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

第二章材料科学研究中常用的数值分析方法;; 材料模拟技术中对应的时间－空间尺度;在科学技术及工程领域，许多力学及物理问题已经得到了反映其规律的基本方程（微分方程）和相应的定解条件。但只有少数性质较简单、边界较规整的问题能通过精确的数值计算得到解析解。; 解决此类问题的通常途径： ①对方程及边界条件加以简化，得到问题在简化条件下的解答；（只在少数情况下有效，因为过多的简化会引起模型失真，导致较大的误差，甚至得到错误的结论） ②采用数值解法，得到模型的数值解。随着计算机技术的发展和应用，数值分析方法已经成为求解科学技术问题的主要工具。 目前，常用的数值分析方法大致分为两大类：有限差分法和有限元法。 ;直接法和间接法。 直接法：精度高，重复工作量小，但编制计算程序复杂，对计算机资源占用较多。 间接法??即迭代法。计算程序简单，占用内存小，但重复工作量大，计算精度取决于迭代次数。;一、直接法：;;;x1=1,x2=2,x3=3.;课堂作业：;2、高斯列主元消去法 ;;;课 后 作 业（高斯列主元消去法）;3、追赶法;方程组的解可用递推公式表示为：;;课堂作业（追赶法）;二、间接法（迭代法）;;1、雅可比(Jacobi)迭代法（简单）;课堂练习;课 后 作 业;2、高斯—塞德尔（Gauss-Seidel)迭代法;;;3、超松弛迭代法;;课 后 作 业;*/24;有限差分方法是计算机数值模拟最早采用的方法，至今仍被广泛运用。 该方法将求解区域划分为差分网格，用有限个网格节点代替连续的求解域。 有限差分法以Taylor级数展开等方法，把控制方程中的导数用网格节点上的函数值的差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。 该方法是一种直接将微分问题变为代数问题的近似数值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。 ;差分的概念：某物理量的增量。 1、向前差分： (1)一阶差分: (2)二阶差分： 一般的m阶差分用m-1阶差分定义： 称为向前差分算子。 ;;3、中心差分;概述 实质是以有限差分代替无限微分、以差分代数方程代替微分方程、以数值计算代替数学推导的过程，将连续函数离散化，以有限的、离散的数值代替连续的函数分布。 ;有限差分法的主要步骤：;差分方程的建立 在建立差分方程前，均需对所求解区域进行离散化。 1、合理选择网格布局及步长 将自变量x,y分别沿x,y轴方向离散，形成离散化网格。网格交点称为节点（或结点），依次将节点编号。;;2、将微分方程转化为差分方程;有限差分法计算的误差分析;;例：利用差分法求解方程;;对此高阶方程组求解，即可得各结点的数值解。;第三节 有限元法;一、有限元法的基本概念－直接刚度法;将求解区域离散化，划分为5个结点和4个单元，每个单元的面积为组成该单元的结点处面积的平均值。 长度为L，具有均一截面A的固体单元在受到外力F时，;广义胡克定律： ;静态平衡要求作用在每个结点上的力的总和为零，得到以下方程：;将以上方程写成矩阵形式，且R1=P。;二、有限元程序的结构和特点;有限元方法是求解偏微分方程的一种数值方法　　偏微分方程是描述客观世界数量关系的一种重要的数学方法。大量的工程、科学、技术和生产问题常常归结为微分方程的求解。以前只能采用解析法求解微分方程，远远不能满足实际的需要。有??元方法原则上可以求解任何复杂的偏微分方程和任何复杂的求解区域问题。 有限元软件是高科技大厦的重要支柱? 有限元软件适用面广，通用性强。随着科学技术和生产的迅速发展，有限元方法的应用越来越广泛，今天各企业、设计、科研部门已普遍采用有限元方法进行生产过程的数值模拟、科研的数值试验和产品的优化设计。;读输入数据、空间分配、数据检查;Two way flow of statistical information;应用领域1—电子行业;应用领域2—水工大坝;应用领域3—桥梁工程;应用领域4—隧道边坡;应用领域5—高层建筑;应用领域6—加工成形;应用领域7—核工业;应用领域8—汽车工业;应用领域9—航空航天;应用领域10—军事工业;应用领域11—石油化工;应用领域12—民用医疗;有限元（FEM）能够准确的拟合构件的边界形状，算法复杂。 有限差分（FDM）不能进行应力分析