一种自适应影响域半径无网格Galerkin法的应用研究

文章在背景积分网格积分方式的基础上,采用基于最小移动二乘近似的一种自适应影响域半径无网格Galerkin法,运用线弹性断裂力学理论,对有限板单边裂纹的应力强度因子进行了分析.由于该方法仅需节点信息,而不需要节点的连接信息,从而避免了有限元方法中的网格重构,大大简化了裂纹扩展的分析过程.数值计算结果表明了该方法的有效性.     

轴对称问题中的无网格Galerkin法

采用无网格Galerkin法分析轴对称问题，得到弹性力学中的对称问题的无网格离散方程．将这一方法与有限元耦合，即在边界处布置有限单元，这样就可以用传统有限元方法方便地处理力学边界条件．算例考察表明：本文方法通过了分片检验，计算结果达到了较高的精度，最大误差不超过5％．     

基于S-R和分解定理的无网格Galerkin法求解几何非线性问题

针对几何非线性问题,基于S-R和分解定理与更新拖带坐标描述法,从运动变换的角度出发,由虚功率原理建立全局弱式积分方程.采用无网格Galerkin法(EFG)求得该问题的数值解.对于非线性问题的数值计算,通常采用增量和迭代法.S-R和分解定理给出了一种新的应变张量,该应变张量能合理地反映大位移大转动情况下物体的应变状态及转动情况,而且有利于求解几何非线性问题.数值算例验证了该方法的合理性和有效性.     

基于无网格局部Petrov—Galerkin方法的h型自适应分析

基于无网格局部Petrov-Galerkin方法进行了h型自适应分析.在进行自适应分析时以Von Mises等效应力作为应力高梯度判据,以最小节点允许距离作为应力高梯度区域加密方案.基于无网格局部Petrov-Galerkin方法实现了对二维线弹性平面问题的h型自适应分析.数值算例表明,基于无网格局部Petrov-Gaterkin方法的h型自适应分析具有较好的稳定性和收敛性.     

非线性Poisson-Boltzmann方程的改进无单元Galerkin法分析

【目的】利用改进无单元Galerkin法求解非线性Poisson-Boltzmann方程。【方法】将改进的移动最小二乘近似与非线性Poisson-Boltzmann方程的Galerkin弱形式耦合,建立了非线性Poisson-Boltzmann方程的改进无单元Galerkin法。基于改进移动最小二乘近似的误差结果下,推导了非线性Poisson-Boltzmann方程的改进无单元Galerkin法的误差估计。【结果】在Sobolev空间中获得了误差估计,并通过数值算例验证了理论结果。【结论】该方法具有较高的计算精度和较好的稳定性,误差随节点间距的减小而降低。     

粘性不可压流场数值模拟的无网格方法

采用无网格伽辽金法(elementfreeGalerkinmethod,EFG)对二维不可压粘性流场进行了数值模拟,主要计算了几个经典的流体力学问题:管道流、绕台阶流以及方腔驱动流。计算中采用标准的Galerkin离散,使用规则分割产生的背景网格来形成积分单元,并与有限元、有限差分法计算所得的速度剖面图进行了比较,结果相差很小。更重要的是,与有限元法相比,无网格虽然计算量比较大,但是它能够处理结点分布畸形甚至重合的情况,从而为解决多相流动直接数值模拟中网格重构的问题奠定基础。     

用径向基函数Galerkin法求两点边值问题

用径向基函数来构造Galerkin型无网格法是一种新型的数值算法,各种数值方法均得到广泛的应用.将其应用于求常微分方程中的两点边值问题,通过对算例中近似解与数值解产生的误差对比,说明用径向基函数Galerkin法稳定性好,且计算精度较高,是一种有效的求解方法.     

无网格法中幂权函数影响半径的优化

借鉴四次样条权函数最优影响半径的选取方法,对工程中常用的幂权函数建立了最优影响半径的计算模型。针对二维情形下的一次基、二次基求解了该模型。通过解二维弹性力学问题实例,分析了不同影响半径下数值解与解析解的误差,表明不同的幂权函数影响半径对数值解的影响是不可忽视的,同时也确认了本文给出的影响半径是可靠的、有效的、综合最优的。     

方形域上的全局同胚自适应网格生成法

给出了一种把任意维方形过渡域自适应映到最终计算域的映射;给出了此映射是全局同胚的充分性条件,以及对应此映射的数值网格生成迭代方法。用全局同胚自适应网格生成的例子进一步表明本文方法是实际可行的。     

紧积分算子特征值无网格Galerkin算法

将无网格Galerkin算法应用到紧积分算子特征值计算中．给出了MLS的基本过程以及无网格Galerkin算法的一般理论,并讨论了无网格Galerkin算法求解紧积分算子特征值问题的误差分析,最后给出特征值及特征向量在该算法下的收敛阶．     

无网格伽辽金法在裂纹扩展计算中的应用

采用了一种基于t-分布的新型权函数，提高了无网格伽辽金法的计算精度；采用完全变换法处理本质边界条件，实现了本质边界条件在节点处的精确施加；针对裂纹扩展中的实际情况，对动态影响半径法作了进一步的补充和改进．算例验证了方法的正确性和有效性．     

关于无单元法中的插值基函数选取的探讨

无单元法不需要单元信息，它采用了一种基于移动最小二乘（MLS）的插值函数。插值基函数对插值函数以及无单元法的计算精度影响很大。本文就不同的基函数对插值函数及无单元法的计算精度的影响作了分析比较，得出了一些有益的结论，并用算例说明了这些结论的正确性。     

无网格法与有限元法的耦合及其对功能梯度材料断裂计算的应用

提出了采用无网格Galerkin法与有限元耦合的方法来计算功能梯度材料中的J积分。这种耦合的方法不仅解决了无网格Galerkin法力学边界条件施加的难点,而且还克服了无网格Galerkin法耗时较多的缺点。此外,文中还给出了一个修正的适合功能梯度材料的J积分,它在功能梯度材料中是独立于积分路径的。在计算过程中,取积分网格中高斯点的材料常数来模拟材料特性的变化。算例的结果显示该方法具有较高的精度和效率。     

金属三维塑性成形过程无网格伽辽金法数值模拟技术

将无网格伽辽金法(EFGM)与三维刚(粘)塑性流动理论相结合,对EFGM在金属三维塑性成形过程数值模拟中的应用技术进行了研究.分别采用边界奇异权法和修正的罚函数法处理速度边界条件和体积不可压缩条件,采用反正切摩擦模型处理摩擦边界条件,推导了金属三维塑性成形过程EFGM法数值模拟的刚度方程,给出了关键算法.对长方体金属镦粗过程进行了数值模拟,并将数值结果与三维刚塑性有限元体积成形商品软件Deform3D计算结果作了比较.发现两者吻合良好,表明了本文方法的正确性和有效性.     

无网格伽辽金法(EFGM)求解接触问题

本文考虑到无网格伽辽金法只需节点信息,不需将节点连成单元,并且有精度高,后处理方便等优点,从而用它求解接触问题。这里采用的方法是将Ｋａｔｏｎａ界面单元引入ＥＦＧＭ,迭代求得两物体间的接触状态。算例表明,本文方法基本可行。     

伽辽金无网格法和有限元法的比较

无网格法在计算力学中成为一种区别于有限元法的新的数值计算方法。文章通过对无单元无网法、应变光滑稳定法、常规有限元法和杂交应力元法进行位移误差和应力误差比较分析;结果表明,无单元无网法和有限元完全积分法在许多问题上精度是可比较的,而有限元中杂交应力元法则明显优于其他方法。     

无网格伽辽金法在施工导流二维数值模拟中的应用

针对有限元法等传统数值计算方法存在受单元网格限制、前后处理工作复杂的问题,提出应用一种数值计算方法--无网格伽辽金法,对具有简单边界条件的水利水电工程施工导流的恒定二维浅水流动问题进行了分析、计算.同时利用有限元法进行了对比计算,从流速、水位的计算结果来看,两种计算方法结果相近、误差较小,表明采用无网格伽辽金法解决此类问题是可行的.     

无单元法应用于欧拉梁的动力特性分析

基于广义移动最小二乘法建立了同时考虑挠度和转角双变量的无单元法用于欧拉梁的动力计算.与传统有限元法相比,该方法只需输入节点信息无需定义单元,具有前处理简单的优势;与只考虑挠度的单变量无单元法相比,该方法具有更高的插值精度.运用双变量无单元法计算了4种不同边界条件欧拉梁的自振圆频率和振型,通过与理论解、有限元解、单变量无单元解的比较,表明无单元法在动力分析中的应用是可行的,欧拉梁的计算同时考虑挠度和转角双变量是必要的,该法在高阶振型计算中具有精度优势.