03 材料科学研究中常用的数值分析方法.ppt

工程材料学 第三章 材料科学研究中常用的数值分析方法 3.1 概述 许多力学问题和物理问题已经得到了它们应遵循的基本规律（微分方程）和相应的定解条件。但是只有少数性质比较简单、边界比较规整的问题能够通过精确的数学计算得出其解析解。大多数问题很难得到解析解。 面临的问题是如何对我们所建立的方程进行求解 解决这类问题通常有两种途径：（1）对方程和边界条件进行简化，从而得到问题在简化条件下的解答；（2）采用数值解法。 第一种方法只在少数情况下有效，因为过多的简化会引起较大的误差，甚至得到错误的结论。 目前，常用的数值解法大致可以分为两类：有限差分法和有限元法。 数值模拟通常由前处理、数值计算、后处理三部分组成 前处理 实体造型、物性赋值、定义单元类型、网格划分 数值计算 施加载荷、设定时间步、确定计算控制条件、求解计算 后处理 显示和分析计算结果、分析计算误差 3.2有限差分法 有限差分法是数值计算中应用非常广泛的一种方法。其实质是以有限差分代替无限微分、以差分代数方程代替微分方程、以数值计算代替数学推导的过程，从而将连续函数离散化，以有限的、离散的数值代替连续的函数分布。 差分方程的建立：首先选择网格布局、差分形式和布局；其次，以有限差分代替无限微分，即以代替，以差商代替微商，并以差分方程代替微分方程及其边界条件。 1.差分方程的建立 合理选择网格布局及步长 将离散后各相邻离散点之间的距离，或者离散化单元的长度称为步长。 将微分方程转化为差分方程 向前差分 向后差分 中心差分 差分格式的物理意义 差分格式的误差分析 2.差分方程的求解方法 直接法－Gauss列主元素消元法 其解为： 间接法－迭代法 对于线性方程组，构造一个值，将代入上式，得出新的值，再将结果代入得到更新的，依次迭代下去，即可使其迭代值收敛于该方程组的精确解。根据选择的方法不同，又可以分为简单迭代法（同步迭代法）和Guass-Seidel迭代法。 对于线性方程组，当，则可表示为下式： 欲求解方程组，首先假设一个解，代入式子的右端，计算出解的一次迭代值，即 4. 有限差分法解题示例 利用差分法解Laplace方程第一边值问题。 3.3 有限单元法 有限元法（FEMA）也称为有限单元法或有限元素法，基本思想是将求解区域离散为一组有限个且按一定方式相互连接在一起的单元的组合体。它是随着电子计算机的发展而迅速发展起来的一种现代计算方法。 把物理结构分割成不同大小、不同类型的区域，这些区域称为单元。 根据不同分析科学，推导出每一个单元的作用力方程，组集成整个结构的系统方程，最后求解该系统方程，就是有限元法。 一.有限元法的基本概念——直接刚度法 例：考虑一个变截面杆，如图所示。杆的一端固定，另一端承受P=1000N的载荷，杆的顶部宽w1=2cm，杆的底部宽w2=1cm，杆的厚度t=0.125cm，长度L=10cm，杆的弹性模量E=10.4×106MPa。试分析该杆沿长度方向不同位置的变形情况，假设杆的质量可以忽略不计。 前处理阶段 将求解区域离散化 先将求解问题分解为结点和单元，如图所示。 杆的横截面积沿y轴方向变化，作为近似，将杆模型化为不同截面的等截面杆的串联，这样，杆可以由一个四个弹簧串联组成的模型来表示，每个单元模型的弹性行为可以用等价线性弹簧来表述： 其中，等价单元刚度为 ，Ai和Ai+1分别是结点i和i+1处的横截面积，l是单元的长度。 应用上述模型，静态平衡要求作用在每个结点上的力的总和为零，因此可得到如下方程： 结点1： 结点2： 结点3： 结点4： 结点5： 将这些方程进行变化，得到： 将这些方程写成矩阵的形式，有： = 求解阶段 杆的截面沿轴的变化可以描述为： 利用上式可以计算出各结点处截面的大小为：A1=0.25， A2=0.21875 ， A3=0.1875 ， A4=0.15625 ， A5=0.125 后处理阶段 分析结果，得到其它信息，如平均正应力。 二.有限元法的基本理论 加权余量法 工程或物理学中许多问题，通常是以未知函数应满足的微分方程和边界条件的形式提出来的，可以一般的表示为未知函数u应满足微分方程组 变分原理和能量法 若设自变量为x，因变量为y，则因变量y是自变量x的函数，可表示为y=y(x) 泛函J的自变量是函数y=y(x)，即泛函是函数的函数，可以表示为J=J(y)=J[y(x)