弹性力学的复变量数值流形方法

数值流形方法通过引入数学和物理双重网格,将插值域和积分域分别定义在两个不同的覆盖上来完成系统能量泛函积分运算. 当采用高阶函数构造位移函数时,广义节点自由度将大大增加. 在求解系统的平衡方程中,运算量是与自由度的三次方成正比的,因此数值流形方法的计算量是较大的. 为此,在复变量理论的基础上,采用一维基函数建立二维问题的逼近试函数,然后将其应用于弹性力学的数值流形方法,提出了复变量数值流形方法,推导了弹性力学的复变量数值流形方法的公式. 与传统的数值流形方法相比,复变量数值流形方法具有计算量小、精度高的优点.      

冲击载荷作用下岩体破坏规律的数值流形方法模拟研究

为了更好地模拟岩体这种典型的连续与非连续介质共存的复杂结构体在冲击载荷作用下的破坏规律,在近年来新出现的数值流形方法基本理论的基础上,利用其算法程序对岩体中存在的不连续面如节理、裂隙等对岩体破裂效果的影响进行了数值模拟分析,计算结果充分显示了数值流形方法在模拟这类问题中的准确性和有效性。     

数值流形方法研究及应用进展

基于有限覆盖技术的数值流形方法是一种新的广义的数值方法.该方法的场函数近似原理和有限元、无网格、单位分解等方法相似,但在网格划分、覆盖形式、近似函数等方面有其自身的特点和优势.对该方法近年来在理论研究和应用方面取得的重要进展进行了综述.在理论研究方面, 目前已对不同形式物理覆盖流形单元的性能进行了研究,结果表明流形单元的精度较有限单元高,且提高覆盖函数的阶次能提高单元的精度;同时理论研究已由二维低阶流形方法推广到三维高阶流形方法,由线性流形方法推广到非线性流形方法,由基于能量原理的流形方法推广到基于加权余量的流形方法,非协调流形方法、无网格流形方法等也已开展了研究; 此外,覆盖系统的自动生成、覆盖函数的形式以及边界条件的处理方法等流形方法相关理论的研究也取得了进展.在应用方面,开展了有关岩石破坏和裂纹扩展等非连续变形分析更深入的研究,并已逐步推广到金属塑性变形分析、多孔介质变形分析以及温度场的数值分析等多个领域.针对目前流形方法的研究和应用现状,该文展望了流形方法理论及实现方法的研究方向、及其在计算流体力学、金属成形等大变形问题、多物理场分析等领域的应用前景.      

数值流形方法的对象设计

研究了数值流形方法的对象设计方法和组织方式,为流形方法的理论研究和向三维问题的扩展打下良好的基础,研究发现数值流形方法具有编程和前后处理简单的特点,且仅用三角形流形单元和一阶近似的覆盖位移函数就可达到有限元多结点等参元的求解精度,具有深远的工程意义,计算结果表明,数值解与理论解吻合。     

爆炸荷载下岩石破坏的数值流形方法模拟

为了更好地利用数值流形方法对动力学问题进行分析,在对原数值流形方法中的动力学问题求解思想进行分析的基础上,采用动力有限元方法中的Newmark法对该算法进行了改进。改进后的数值流形方法与原来相比具有三个明显的优势：（1）当选择合适的参数后,该方法能够保证解的无条件收敛;（2）可以采用比原算法大得多的时间步长;（3）充分考虑了动力学问题中的阻尼效应。最后通过一个算例说明了改进后的数值流形方法能够很好地模拟岩石在冲击载荷作用下破坏的全过程,克服了有限元法不能模拟岩石破坏后块体运动情况的不足。     

数值流形方法的粘性边界问题初探

在实际工程数值流形方法分析中,采用固定约束边界的方法处理无限域或者半无限域的情况,边界处应力波的反射造成模拟结果与实际情况不符.本文基于Lysmer等人提出的粘性边界理论,在边界上设置阻尼器,推导相应粘性边界条件下流形单元刚度矩阵的数值计算格式,经岩石长条中弹性波传播算例,并与有限元结果对比,验证了该粘性边界的有效性,有利于数值流形方法的工程中推广应用.     

基于有限变形理论的数值流形方法研究

原有数值流形方法通过积累每一时步的小变形而得到结构最终的大变形,然而,当结构发生大变形、大转动时往往产生较大计算误差. 针对该问题,从动量守恒方程以及应力边界条件的积分弱形式出发,引入流形方法的插值函数,建立了基于有限变形理论的数值流形方法. 通过对比改进前后流形方法的计算迭代格式,指出了原有流形方法计算大变形问题时的误差来源. 最后,通过大变形悬臂梁和旋转块体算例对有限变形流形方法进行了验证. 数值结果表明,改进后的流形方法能够很好地处理大变形大转动问题,消除了转动所带来的计算误差,其计算结果与解析解及ABAQUS 软件求得的数值解相吻合.      

数值流形方法及其在岩石力学中的应用

数值流形方法是目前岩石力学分析的主要方法之一.该方法起源于不连续变形分析,主要用于统一求解连续和非连续问题,其核心技术是在分析时采用了双重网格:数学网格提供的节点形成求解域的有限覆盖和权函数;而物理网格为求解的积分域.数学网格被用来建立数学覆盖,数学覆盖与物理网格的交集定义为物理覆盖,由物理覆盖的交集形成流形单元.流形方法的优点在于它使用了独立的数学和物理网格,具有和有限元明显不同的定义形式,且数学网格对于同一问题不同的求解精度的需求可以很方便地细化.由于该方法考虑了块体运动学,可以模拟节理岩体裂隙的开裂和闭合过程,因而在岩石力学中得到了广泛应用,近年来许多学者对该方法进行了研究.本文简要叙述了节理岩体的数值方法从连续到非连续的发展过程,详细地介绍了数值流形方法的组成和数值流形方法在岩石力学及其相关领域的研究和发展概况,最后就作者所关心的一些问题,如三维问题的数值流形方法、数值流形方法在物理非线性问题和裂纹扩展问题中的应用、相关的耦合方法等进行了探讨.      

弹性力学的复变量无网格流形方法

为克服无网格流形方法配点过多、计算速度慢、容易形成病态方程组等缺点,将复变量移动最小二乘法与无网格流形方法相结合,提出了弹性力学的复变量无网格流形方法。分别采用线性基本与二次基进行计算,并与无网格流形方法相比。研究表明该方法计算量小、精度高。     

数值流形单元法数学网格自适应

基于数值流形方法和有限覆盖技术,将有限元法的后验误差估计理论及h型网格自适应技术推广应用到数值流形单元法中,提出了数值流形单元法的后验误差估计方法和数学网格自适应技术,并编制了相应的程序。数值算例表明,经过网格自适应,可以在粗糙的初始网格基础上得到质量比较理想的网格,计算结果可达到用户要求的精度。     

高次数值流形方法及其在弹性地基板中的应用

根据数值流形方法的思想,用Hermite插值函数表示三角形板单元位移模式中的权函数,并在此基础上,在不增加额外覆盖数的条件下,构建了一种新的更高次痊移模式,各种条件的计算结果表明这种位移模式具有精度高的,优点,很有实用意义。     

数值流形方法在连续体数值分析中的应用

在平面连续作的数值计算中使用数值流形方法的各阶覆盖位移函数和统一的三节点三角形流形单元,能够得到高精度的数值结果,且具有编程和前后处理简单等优点,克服了有限元法复杂理论的不足．文中还引入了三角形单元Ｈａｍｍｅｒ数值积分,使流形方法的程序编制通用化．     

基于有限变形理论的数值流形方法研究简

原有数值流形方法通过积累每一时步的小变形而得到结构最终的大变形,然而,当结构发生大变形、大转动时往往产生较大计算误差. 针对该问题,从动量守恒方程以及应力边界条件的积分弱形式出发,引入流形方法的插值函数,建立了基于有限变形理论的数值流形方法. 通过对比改进前后流形方法的计算迭代格式,指出了原有流形方法计算大变形问题时的误差来源. 最后,通过大变形悬臂梁和旋转块体算例对有限变形流形方法进行了验证. 数值结果表明,改进后的流形方法能够很好地处理大变形大转动问题,消除了转动所带来的计算误差,其计算结果与解析解及ABAQUS 软件求得的数值解相吻合.     

基于适合分析T样条的高阶数值流形方法

数值流形方法是一种非常灵活的数值计算方法,连续体的有限单元方法和块体系统的非连续变形分析方法只是这一数值方法的特例.数值流形方法中高阶位移函数的构造可通过提高权函数的阶次来实现,这种方法往往需要沿单元边界配置适当的边内节点,这些结点的出现增加了前处理的复杂性,特别是对于大型复杂的空间问题.另一方面,在数值流形方法中可通过缩小单元尺寸(h加密)来提高求解精度.当模拟裂纹扩展时,这种细化策略可用来克服裂纹尖端的奇异性.一个传统的解决方案是细化整个网格,但这会导致计算效率的显著降低.将适合分析的T样条(analysis-suitable T-spline,AST)引入数值流形方法中来建立高阶数值流形方法的分析格式,有效的避免了该问题的出现.AST样条基函数具有线性无关,单位分解,局部加密等许多重要性质,使得其非常适合用于工程设计及分析.在引入AST样条后,可通过改变数学覆盖的构造形式建立不同阶次的数值流形方法分析格式;AST样条自身的局部加密性质也使得数值流形方法中的数学网格局部加密更容易实现.算例结果表明:随着AST样条基函数阶次的提高,数值流形方法的计算结果有了明显的改善;基于AST样条基函数的数值流形方法在保持计算精度的前提下降低了自由度的数量.     

固定网格的数值流形方法研究

针对单纯几何非线性的材料大变形问题, 提出一种新的研究思路------固定数学网格的数值流形方法, 简称固定网格流形法, 可以看作是采用了固定网格的拉格朗日方法. 它充分利用数值流形方法的数学网格与材料物理边界分离的特性, 具备拉格朗日法和欧拉法各自的优势, 避免了原始拉格朗日法的网格扭曲问题以及欧拉法对移动边界难以精确描述和迁移项较难处理的问题. 采用数值流形方法的大变形分步计算格式, 使得固定网格流形法实现起来并不复杂, 仅需要每步切割网格形成新的流形单元, 以及对初应力载荷进行适当的处理, 而后者是固定网格流形法的关键. 针对固定的矩形数学网格开展研究, 采用一阶多项式覆盖函数的高阶流形法, 给出了两种初应力计算方法, 并用悬臂梁大变形算例验证了固定网格流形法的可行性, 将来需要进一步解决初应力载荷所带来的计算稳定问题.      

多裂纹扩展的数值流形法

数值流形法的求解体系建立在两套覆盖(包括数学覆盖和物理覆盖) 和接触环路的基础之上,实现了对连续和非连续问题的统一求解. 在处理裂纹问题时,数学覆盖无需与裂纹重合,方便岩体破坏过程的模拟. 通过在裂纹尖端影响区域内的物理片上增加用于模拟应力奇异性的增强位移函数,发展了扩展的数值流形法. 在此基础上,提出一种多裂纹扩展的控制算法,并给出了裂纹扩展过程中材料体的整体响应. 针对典型的线弹性断裂力学问题, 给出的数值算例表明所建议的方法是正确有效的.      

EVALUATION ON STRESS INTENSITY FACTOR OF CRACK UNDER DYNAMIC LOAD USING NUMERICAL MANIFOLD METHOD

相较于传统有限元,数值流形方法（numerical manifold method, NMM） 的一个显著优点是在处理裂纹问题时网格无需与裂纹重合,这就方便了岩体破坏过程的模拟. 基于包含裂尖增强函数的NMM,采用Newmark 隐式动力学算法进行时间积分,重点研究了动力载荷条件下裂纹动态应力强度因子（dynamic stress intensity factor,DSIF） 的求解方法. 针对典型的线弹性动力裂纹问题,给出了NMM 的数值算例. 结果表明NMM 能够准确计算动载荷条件下裂纹的DSIF,并且具有较好的收敛性.