(新)计算材料学——精品.docVIP

得 分 评卷人 沈阳师范大学 硕士研究生期末考试答题册 考试科目： 计 算 材 料 学 姓 名： 刘芳慧 学 号： 培养单位： 沈阳师范大学 专 业： 材料物理与化学 考试时间： 2013. 6. 2 沈阳师范大学研究生处印制 有限元法及其应用 1. 引言 有限元法是计算力学中的一种高效能、常用的计算方法，兴起于20世纪50年代末60年代初，主要应用在数学、现代力学以及与计算机科学相互渗透、综合利用的边缘科学。所谓有限元，通俗的讲就是对一个真实的系统用有限个单元来描述。而有限元法，则是把求解区域看作由许多小的在节点处相互连接的单元所构成，其模型给出基本方程的子域近似解，由于单元可以被分割成各种形状和大小不同的尺寸，所以它能很好地适应复杂的几何形状、复杂的材料特性和复杂的边界条件。最初应用在工程科学技术中，用于模拟并且解决工程力学、热学、电磁学等物理问题。由于有成熟而强大的软件系统的支持，有限元法已成为一种非常受欢迎的、应用极广的数值计算方法[1]。 有限元法的基本思想是将连续的求解区域离散为一组有限个、且按一定方式相互联结在一起的单元组合体。由于单元能按不同的联结方式进行组合，且单元本身可以有不同形状，因此可以模型化几何形状复杂的求解区域。有限元法作为数值分析方法的一个重要特点是利用在每一个单元内假设的近似函数，分片地表示全求解域上待求的未知场函数。单元内的近似函数通常由未知场函数或其导数在单元的各个节点的数值和其插值函数表达。这样，一个问题的有限元分析中， 未知场函数或其导数在各个节点上的数值就成为新的未知量即自由度，从而使一个连续的无限自由度问题变成离散的有限自由度问题。一经求解出这些未知量，就可通过插值函数计算出各个单元内场函数的近似值，从而得到整个求解域上的近似解。显然，随着单元数目的增加，即单元尺寸的缩小，或者随着单元自由度的增加及插值函数精度的提高，解的近似程度将不断改进。如果单元是满足收敛要求的，近似解最后将收敛于精确解。 2. 有限元法的起源与发展 有限元法基本思想的提出，可以追溯到Courant在1943年的工作。他首先尝试应用在一系列三角形区域上定义的分片连续函数和最小能原理相结合，用来求解St.Venant扭转问题。此后，不少应用数学家、物理学家和工程师分别从不同角度对有限元法的离散理论、离散方法及应用进行了研究。有限元法的实际应用是随着电子计算机的出现而开始的。首先是Turner，Clough等人于1956年将刚架分析中的位移法推广到弹性力学平面问题，并用于飞机结构的分析。他们首先给出了采用三角形单元法求解平面应力问题的正确解答。三角形单元的刚性矩阵和结构的求解方程是由弹性理论的方程通过直接刚度法确定的。他们的研究工作开始了利用电子计算机求解复杂弹性力学问题的新阶段。1960年Clough进一步求解了平面弹性问题，并第一次提出了“有限元法”名称，使人们更清楚地认识到有限元法的特性和功能。 近四十多年来，伴随着电子计算机科学和技术的快速发展，有限元法作为工程分析的有效方法，在理论、方法的研究、计算机程序的开发以及应用领域的开拓诸方面均取得了根本性的发展。它不仅成功地运用于固体静力学，近年来也应用于解决动力学问题。在瞬态动力学数值计算机里，通常采用集中质量法代替一致质量矩阵，将求解过程化为显式计算，这样可避开复杂的矩阵运算，从而大大减少计算工作量[2]。 目前在进行大型复杂工程结构中的物理场分析时，为了估计并控制误差，常用基于后验误差估计的自适应有限元法。基于后处理法计算误差，与传统算法不同，将网格自适应过程分成均匀化和变密度化2个迭代过程。在均匀化迭代过程中，采用均匀网格尺寸对整体区域进行网格划分，以便得到一个合适的起始均匀网格；在变密度化迭代过程中只进行网格的细化操作，并充分利用上一次迭代的结果，在单元所在的曲边三角形区域内部进行局部网格细化，保证了全局网格尺寸分布的合理性，使得不同尺寸的网格能光滑衔接，从而提高网格质量。整个方案简单易行，稳定可靠，数次迭代即可快速收敛，生成的网格布局合理，质量高。 3. 有限元法的分类 有限元法可以分为两类，即线弹性有限元法和非线性有限元法。其中线弹性有限元法是非线性有限元法的基础，二者不但在分析方法和研究步骤上有类似之处，而且后者常常要引用前者的某些结果。 线弹性有限元是以理想弹性体为研究对象的，所考虑的变形建立在小变形假设的基础上。在这类问题中，材料的应力与应变呈线性关