两步型taylor-two - step taylor.docxVIP

摘要动态问题在科学技术和国民经 济的发展中有着广泛的应用背景 ，这类问题多 呈现材料非线性或几何非 线性且边界条件极为复杂，通常无法得到解析 解答，而 需要借助数值算法来进行 近似模拟和求解。目前求解动态问题的数值方 法主要为 有限单元法和无网格法， 有限单元法存在精度偏低以及处理网格畸变能 力不足的 问题，无网格法虽然处理 网格畸变能力较强，但其计算效率较低，很难 用于实际 工程计算。为了兼顾上述 算法的计算能力与计算效率，充分利用低阶单 元在解决 工程问 题时的巨 大优势， 本文从光 滑有限元 法和传统 Taylor-Galerkin 算法着 手， 提出一 种新型两 步型 Taylor-Galerkin 光滑有 限元法。 该方法基 于梯度光 滑技术， 有效提 高低阶单 元的精度 ，并利用 两步型 Taylor-Galerkin 算法的 动量及能 量守恒 性，在求解动态问题时能够展现良好的性能。针对不 同类型的 动态问题 ，本文将两 步型 Taylor-Galerkin 算法与 不同的梯 度 光滑算法相结合，并将它 们应用于弹性波传播问题和超弹性材料的大变 形问题的 研究中，通过数值算例的分析对它们的性质特点进行比较。本文的主要工作包含：(1)将两步型 Taylor-Galerkin 算法与边 光滑有限元法 结合，并将其用 于弹性波 在介质中的传播等问题的 分析中。该方法突破传统单元积分的限制，基 于三角形 单元的边形成基于边的光 滑积分域，有效地软化系统数值模型，提高计 算精度； 更值得注意的是：此算法 基于低阶线性三角形单元，计算流程简单，不 需要增加 任何内部变量或自由度，在保证高精度的同时具备较高的效率。(2)将两步型 Taylor-Galerkin 算法与基 于单元的光滑 有限单元法结 合，并运用 于超弹性材料大变形等问 题的求解中。在该算法中，数值积分在基于单 元形成的 光滑域中进行，克服了传 统四边形单元由于等参元的使用而对网格严重 畸变问题 产生局限性的缺点，对处 理网格严重畸变的问题具备独特的优势；此外 ，该方法 还具备了 两步型 Taylor- Galerkin 算法的能 量守恒性质 ，在复杂的 迭代过程中，能 量守恒可保证迭代的稳定 性，减少或避免高频振荡成分对计算的影响， 从而降低 计算结果的误差。(3)对 上 述 两 种 方 法 的 性 质 和 特 点 通 过 具 体 算 例 分 析 表 明 两 步 型 Taylor-Galerkin 边 光 滑 有 限 元 法 由 于 光 滑 域 的 形 成 方 式 以 及 稳 定 性 等 因 素 ， 较 适 用于弹 性波在介 质中的传 播问题； 两步型 Taylor-Galerkin 单元光 滑有限元 算法则 由于改善了传统等参元所存在的缺陷因此更适用于动态大变形问题的解答。关键词：两步型 Taylor-Galerkin 算法；梯度光滑有限元法；光滑域；能量守恒AbstractThe dynamic problems analyses are of great importance for science, technology and economy. As these problems are usually nonlinear or have extremely complex boundary conditions, the analytical solutions are not always available and numerical methods are employed. There are mainly two types of numerical methods, the finite element method and the meshless method. Due to low accuracy of lower order elements and poor capacity in dealing with elements distortion, the traditional finite element method has difficulty in solving complicated problems. The meshless methods can work perfectly to overcome disadvantages existed in finite element method, but they have problems of complex computational process and lower efficiency as well. In order to obtain hi