基于节点计算的数值分析方法在我国的起源.pdfVIP

基于节点计算的数值分析方法在我国的起源 孙海涛王元汉 (华中科技大学土木工程与力学学院，湖北武汉430074) 摘要基于节点计算的方法在解决传统有限元法处理较困难的领域 取得了明显的成效，受到普遍关注。中国学者在此类方法的发展中，从 B-Spline群条函数应用的角度，得到了，一种完全基于节点计算的数值分析方 法。本文证明了国内学者的方法与国外研究的等价性。国内学者所提方法 的计算结果在一些问题中显示出比国外学者所提方法有较大的优越性。 关键词数值分析，节点，控制方程，样条函数，逼近 1 引言 数值分析方法是许多工程和科学技术问题得以求解的重要方法，其 核心思想是用简单的数值插值或逼近代替复杂的真实场函数，通过一定 形式的方程控制形式，旌加给定的边界或初值条件后，将待求解问题转化 成线性方程组求解。按方程控制形式的不同，数值分析方法可分为强形 式和弱形式的数值分析方法。典型的强形式的数值分析方法是最早由 (FDM)，其他还有微分求积法E2I。强形式数值分析方法的基本特征是用 待求节点附近有限个节点上场函数值的线性组合，逼近场函数的各阶导 数(记函数自身为零阶导数)，然后，将此简化后的场函数逼近关系代入控 制微分方程，进而得到一组代数方程，对此代数方程组求解，即可得到给 定节点的函数值。FDM不需数值积分计算，是数值分析方法中最早得到 应用的方法，已经发展出各种函数及其导数插值格式，至今仍在计算流体 力学等领域中有着重要的应用，其求解过程为构造场函数及其各阶导数 逼近格式，施加边界条件，求解代数方程。但FDM在不规则域上边界条 ·226· 件难以精确满足，求解精度较低，满足函数及其各阶导数的插值逼近形式 是同一线性组合关系的插值格式很难找到，而且不同的问题，其程序编制 难于统一，因而其应用严重受限，在很多领域被弱形式的有限元法 在工程和科学技术中得到了广泛的应用，发展得最为活跃，影响最大的 FEM是典型的弱形式数值分析方法的代表。因为弱形式的求和特征， FEM需求解单元域上的积分计算问题，这点由于数值积分技术的发展， 总体上已得到很好的解决，其求解过程为按选定规则选择弱形式的控制 方程，构造场函数的逼近格式，施加边界条件，完成数值积分形成代数方 程，代数方程求解。但由于FEM分片插值特性，其数值解的连续性受到 分片插值函数的连续阶影响，最终数值解需作平滑处理(即所谓后处理)。 另外，FEM计算严重依赖于计算网格，计算网格的产生(所谓前处理)本 身就是一件十分困难的工作，特别当网格畸变严重时，FEM处理就会很 困难。再则，连续性问题导致的数值锁也是FEM的固有缺点。边界元法 (BEM)是另一种弱形式的数值分析方法，它具有将待求问题从全部求解 域向求解域边界转移的特点，而且所选取的基本解都能优先满足指定的 边界条件，因此在计算量和计算精度上较有优势，对一些无限域和半无限 域问题，比FEM有效。但边界元方法的核心工作——基本解的求取—— 十分困难，目前的结果也十分有限，对非线性问题还无能为力，程序的通 用性也差，还未产生大规模的应用结果。 弱形式的数值分析方法极大地扩张了场函数的选择空间，使得场函 数的选择更灵活多变。正是这一特点，弱形式的数值分析方法仍然继续 发展，其中基于节点计算的无网格法是过去十年中取得较大成功的数值 分析方法的代表。本文所指的基于节点的计算方法即是以弱形式表现的 无网格法数值方法。虽然许多国内学者在无网格法的其他方法中取得了 很好的成就，但本文仅从具有方法性原创意义角度探讨基于节点计算的 数值方法在我国的起源。 2基于节点计算的无网格法 数值分析方法最基本的工作是对场函数的数值逼近，基于节点计算的 ·227· 数值分析方法中，以无网格法最引人关注。无网格法的研究也是从场函数 逼近形式开始。无网格法有多种场函数逼近格式，基于散点数据逼近、函数 一性特点的场函数逼近格式等，在无网格法中都得到应用。每种插值格式 都形成一种类型的无网格法，若以这些插值形式的历史作为无网格法研究 起源，则有些牵强。较为一致的看法是：数值方法明确以不需要划分单元为 目的。据此标准，无网格法历史可以从20世纪70年代对非规则网格的 质点流体动力学(SPH)方法，成功地应用于天体物理、冲击波模拟、流体动