(计算机在材料科学与工程中的应用)第三章材料科学与工过程中的数学模型和数值分析.ppt

3.3 常用的数值计算方法—3.3.2 有限差分法求解 3.3 常用的数值计算方法—3.3.2 有限差分法求解 3.3 常用的数值计算方法—3.3.2 有限差分法求解 3.3 常用的数值计算方法—3.3.2 有限差分法求解 3.3 常用的数值计算方法—3.3.2 有限差分法求解 间接法 对于线性方程组Ax=b，构造一个 值，将 代入， 得出新的值 ，再将结果代入得到更新的 ，依次迭代下去，即可使其迭代值收敛于该方程组的精确解 。根据选择 的方法不同，又可以分为简单迭代法（同步迭代法）和Guass-Seidel迭代法。 对线性方程组Ax=b，当 时，可表示为 （3.89） 可改写为 （3.90） 3.3 常用的数值计算方法—3.3.2 有限差分法求解 欲求解方程组，首先假设一个解 ，代入式（3.90）的右端，计算出解的一次迭代值，即 再将 代入式子的右端，得到第二次迭代值，以此类推，得到第k次迭代值 即可认为迭代已经满足精度要求。其中c为某适当小的量，其具体大小取决于精度要求。 迭代次数无限增多时， 将收敛于方程组的精确解 。一般满足 3.3 常用的数值计算方法—3.3.2 有限差分法求解 Eg. 边界点温度已知，且区域内无内热源。利用有限差分方法来计算节点1、2、3的温度。 3.3 常用的数值计算方法—3.3.2 有限差分法求解 Eg. 边界点温度已知，且区域内无内热源。利用有限差分方法来计算节点1、2、3的温度。 解方程得 3.3 常用的数值计算方法—3.3.2 有限差分法求解 有限元法（Finite Element Method，FEM），也称为有限单元法或有限元素法，基本思想是将求解区域离散为一组有限个、且按一定方式相互连接在一起的单元的组合体。它是随着电子计算机的发展而迅速发展起来的一种现代计算方法。 把物理结构分割成不同大小、不同类型的区域，这些区域就称为单元。根据不同分析科学，推导出每一个单元的作用力方程，组集成整个结构的系统方程，最后求解该系统方程，就是有限元法。简单地说，有限元法是一种离散化的数值方法。离散后的单元与单元间只通过节点相联系，所有力和位移都通过节点进行计算。对每个单元，选取适当的插值函数，使得该函数在子域内部、子域分界面上（内部边界）以及子域与外界分界面（外部边界）上都满足一定的条件。然后把所有单元的方程组合起来，就得到了整个结构的方程。求解该方程，就可以得到结构的近似解。 离散化是有限元方法的基础。必须依据结构的实际情况，决定单元的类型、数目、形状、大小以及排列方式。这样做的目的是：将结构分割成足够小的单元，使得简单位移模型能够足够近似地表示精确解。同时，又不能太小，否则计算量很大。 有限元法 3.3 常用的数值计算方法—3.3.3 有限元法 有限差分法：直观，理论成熟，精度可选。但是不规则区域处理繁琐，虽然网格生成可以使FDM应用于不规则区域，但是对区域的连续性等要求较严。使用FDM的好处在于易于编程，易于并行。 有限元方法：适合处理复杂区域，精度可选。缺憾在于内存和计算量巨大。并行不如FDM直观，FEM的并行计算是当前和将来应用的一个不错的方向。 3.3 常用的数值计算方法—3.3.3 有限元法 作业： 上机实验： * * 同时形核：同样的方法可用来处理三维晶球，这时把圆环确定的有效面积增量用球壳确定的有效体积增量 来代替，对于同时成核体系（N为单位体积的晶核数），则 （3.52） 3.2 材料科学与工程研究中数学模型—类比分析法 三维形核 不断成核体系：定义I为单位时间、单位体积中产生的晶核数，则 将上述情况进行归纳，可用一通式： 式中，k是同核密度及晶体一维生长速度有关的常数，称为结晶速度倍数；n是与成核方式及核结晶生长方式有关的常数。该式称为Avrami方程。 （3.53） （3.54） 模型检验： 图3.4为尼龙1010等温结晶体数据的Avrami处理结果，可见在结晶前期实验同理论相符，在结晶的最后部分同理论发生了偏离。 解释：因为生长着的球晶面相互接触后，接触区的增长即停止。在前期球晶尺寸较小，非晶部分多，球晶之间不致发生接触，随时间延长，球晶增长到满足相互接触的体积时，总体的结晶速度就要降低，Avrami方程将出现偏差。 3.2 材料科学与工程研究中数学模型—类比分析法 3.2 材料科学与工程研究中数学模型—数据分析法 4 数据分析法 当系统的结构性质不大清楚，无法从理论分析中得到系统的规律，也不便于类比分析，但有若干能表征系统规律、描述系统状态的