数值流形方法中覆盖函数选用的建议

在数值流形方法中,常用的覆盖函数基并非是最佳选择。循着刚度矩阵的形成过程,分析了选用常用覆盖函数时,在非对角元上出现绝对值很大的元素之成因,且发现这会增加刚度矩阵的条件数,尤其是在用刚性弹簧约束位移的情况下,而这在数值流形方法中普遍而基本。建议采用局部化较好的覆盖函数,取代常用的关联于全局坐标的覆盖函数,可显著消除这一情况。建议方式简单明了,程序改动极小,对改善刚度矩阵性态却有很大作用。算例验证了这一建议的合理性,通过比较局部化的覆盖函数及全局性的覆盖函数所形成的刚度矩阵,表明前者形成了较小条件数的刚度矩阵。     

数值流形方法中覆盖位移函数的改进

通过对数值流形方法刚度矩阵形成的研究,提出一种改进的覆盖位移函数,以改善刚度矩阵计算。分析表明,改进的覆盖位移函数可使刚度矩阵的局部大数数量级明显的降低,提高计算效率,给矩阵求解提供适当的意见,并给出算例验证方法的正确性。     

数值流形方法在P型自适应分析中的初探

将数值流形方法应用于P型自适应分析，推导单元不同阶次的位移函数及其单元刚度矩阵表达式，引入了对应情况下的节理单元及其刚度矩阵，编制了计算程序，数值算例表明该算法可行。     

数值流形法的T样条局部加密算法

在岩土工程分析中求解精度控制常常是必需的,在数值流形法中可以通过控制数学覆盖网格的稀疏和覆盖位移的阶数来达到精度的要求。提出了基于等几何分析的数值流形方法,定义了相应的数学覆盖的构造形式,推导了基于二次B样条的9节点数值流形方法分析格式;针对基于Lagrange插值函数的4节点数值流形方法提出了基于T样条思想的数学覆盖网格的局部加密方法。算例计算结果表明,相对于4节点的数值流形方法,基于非均匀有理B样条的9节点数值流形方法具有更高的精度;基于T样条思想的加密网格在保持计算精度的前提下降低了自由度的数量,表明T样条加密是一种自然的局部加密算法。     

弹塑性数值流形方法在边坡稳定分析中的应用

针对现有数值流形方法只能进行弹性计算的不足,建立了一个能够反映完整岩块弹塑性变形特征的本构模型,并借助VC++开发了内置该本构模型的弹塑性数值流形程序。利用该程序模拟了含单节理岩样的室内压缩试验,分析得到了其强度和变形特性,计算结果符合实际的物理现象,表明程序是正确有效的。考虑到数值流形方法本身能够有效模拟材料的不连续变形,新增的弹塑性分析功能又可以反映岩石的强度特性,将弹塑性数值流形程序应用于某含有不连续面的岩石边坡的稳定性分析。并结合锚杆单元的使用,对比分析了不同锚固方案的加固效果。程序提供的变形、应力等计算结果表明:预应力锚杆不仅可以防止不连续面发生剪切破坏,增强坡体的稳定性,限制塑性变形的发展;而且可以使不连续软弱层面对岩体变形的消极影响得以减弱,起到提高岩体整体性的作用。     

三维数值流形法覆盖系统生成算法研究

覆盖系统生成效率低是当前制约三维数值流形法在大型实际工程中应用的瓶颈问题。通过改进三维覆盖系统的生成算法,提出了改进的布尔交运算生成三维覆盖系统的算法,并采用C++语言编写了相应的程序。根据数学网格与物理域的拓扑几何关系,确定采用传统布尔交运算还是由数学网格直接生成流形块体,并对生成的流形块体进行三维块体有效性检查,满足要求后即可生成新的三维流形单元,进而生成所有前处理三维覆盖系统。选取两个简单算例和一个边坡工程问题来说明所提方法的可行性和计算程序的正确性,并通过对比分析说明方法的计算效率。计算结果表明,改进的布尔交运算比传统布尔交运算更高效;随着网格密度的增加,生成一个单元的平均耗时逐渐减少;随着物理域边界复杂程度增加,所提方法计算效率逐渐降低。为今后采用数值流形法进行结构分析奠定基础,具有较强的实际应用价值。     

耦合渗流-变形的数值流形法裂隙岩质边坡稳定性分析

岩石裂隙渗流问题作为影响岩石边坡失稳的重要因素吸引了越来越多研究者采用不同的数值方法模拟其力学行为。提出了一种基于渗流-变形耦合模型的数值流形方法,并相应地讨论了岩体裂隙网络中地下水的流动、渗透压力和岩石变形的耦合效应。基于最小能量原理,分别推导渗流条件下单元边界流体压力和系统的整体平衡方程,以及块体接触算法求解任意岩石裂缝的局部安全系数。通过一系列验证实例对所提出的数模模型进行检验,验证了该数值模型的鲁棒性和有效性。利用该方法模拟了一个岩石边坡由于渗流作用而发生的坍塌事故,再现了岩石边坡的破坏过程。模拟结果表明,过大的水力压力导致垂直裂缝张开,增加岩体变形量,最终导致边坡破坏。研究表明,该方法在模拟大规模工程问题方面具有很大的潜力。     

基于加密物理片的数值流形法中局部网格加密

数值流形法的数学网格不需要适应求解域内各种不连续界面和边界,因此,总可以用规则的结构化网格建立数学覆盖。但对于大多数问题,在整个求解域上布置统一密度的网格显得浪费。因此,需要研究在结构化网格上实施局部加密,建立了加密物理片的方法用以解决这一问题。具体的实施过程中,确定了需要加密的区域后,先在这些区域布置规则的精细网格,然后找到这些区域中包含的原始网格中的物理片,用精细网格上建立的插值代替被加密物理片上的局部近似,从而提高了局部近似的阶次。数值算例结果表明,该方法收敛性良好。另外,如果所有物理片上的局部近似都采用0阶多项式(常数),那么将会得到正定的刚度矩阵。     

弹性静力学问题的边界积分形式的数值流形法

传统的数值流形法(NMM)一般均采用区域积分形式。结合边界单元法(BEM),提出了一种边界积分形式的数值流形法。该方法既能发挥NMM的可以灵活选取局部基的优势,又具有BEM降低问题求解维数的特点。针对二维的弹性静力学问题,对3个具有解析解的不同基准算例进行了数值应用,验证了所提方法的有效性和效率。计算结果表明,提高局部基的阶次可有效提高方法的计算精度。     

水文地质数值法存在的问题及其对策

本文提出了水文地质数值方法当今存在的两大带有全局性，根本性的问题，并针对这些问题指出了解决途径。     

三维数值流形法在裂纹扩展中的应用研究

运用三维数值流形法（3D NMM）进行三维裂纹扩展分析,并采用C++语言编写了相应的程序。充分利用三维数值流形法模拟裂纹扩展的优势,只需要更新裂纹尖端线附近的边界环路和流形单元,不需要使用阶跃函数。根据三维数值流形法计算得到的应力结果,应用非局部求迹方法分析每个裂纹尖端的破坏状态,如果发生破坏则沿垂直于其最大主应力方向扩展。针对裂纹扩展后的不同状况,采用四边形或三角形推进法。裂纹扩展后为了使变形后的面保持平面,必须对新生成的面进行三角化分割。对诸如单边裂纹、平行钱币型裂纹和倾斜钱币型裂纹扩展问题进行数值模拟。计算结果表明,采用三维数值流形法进行裂纹扩展模拟是可行的,文中方法对裂纹尖端线非闭合和闭合的情形均适用,且文中方法对于裂纹尖端线位于单元内部的非平面裂纹扩展也是有效的。     

基于改进数值流形法的接触裂纹问题研究

接触裂纹问题在工程结构中较为常见。结合数值流形法在裂纹处理上的优势,分析了压剪荷载作用下的接触裂纹问题,模拟了压剪裂纹渐进扩展过程。为了减少由于裂纹尖端位置不同而产生的误差,对裂纹尖端附近一定范围内的每一个物理覆盖附加奇异覆盖函数项,并根据裂纹尖端位置和单元含奇异物理覆盖的数目进行分区积分。选取一个压剪破坏算例,分析了法向接触力对应力强度因子计算结果的影响,并模拟了其渐进破坏过程。计算结果表明,所提方法在压剪裂纹问题方面的可行性,与未细化和覆盖细化方法得到的结果相比,更能准确地描述裂纹扩展路径。法向接触力对II型应力强度因子的贡献为0,对I型应力强度因子的影响较大,相对误差随网格密度变化明显,且法向接触力对I型应力强度因子的影响要比直接施加内压时的影响大。     

基于数值流形法的三维裂纹扩展研究

基于数值流形法（NMM）进行三维裂纹扩展分析研究,编写了相应的C++程序。充分发挥数值流形法在非连续变形分析领域的优势,不需要扩展有限元中的水平集和阶跃函数概念,应用数值流形法计算结果分析裂纹尖端线的破坏状况,对已有的非局部求迹法和三角形推进法进行简化和扩展,提出一种简化算法确定最终的新裂纹扩展面。应用简化算法对水平钱币型裂纹扩展问题进行数值模拟和对比分析研究。计算结果表明,基于数值流形法的三维裂纹扩展算法是可行的,采用简化处理之后,克服原有算法的不足,大大减少了新裂纹尖端线和新裂纹面的数目,降低了计算网格对新裂纹面的影响,提高了计算效率和扩大了应用范围。     

数值流形法中基于适合分析T样条的局部网格加密算法

数值流形方法中一般采用有限元网格或规则网格作为其数学覆盖系统,而规则的网格突出的优点是不需要适应求解域的边界和各种不连续面。采用规则的矩形网格作为数值流形方法中的数学网格,并借助适合分析的T样条实现了数值流形方法中的局部加密。适合分析的T样条定义在一个限制的T网格上,其基函数具有线性无关、单位分解、局部加密等许多重要性质,使得其非常适合用于工程设计及分析。当对一个适合分析的T网格加密后,所产生的新的网格往往不再是适合分析的T网格。基于此,提出了一种简单的数学网格加密算法,该算法能保证局部加密后的数学网格仍然是适合分析的。算例结果表明:在应力集中区和裂纹尖端等应力梯度较大区域,该算法均具有较强的适用性。     

改进的数值流形法在水力劈裂中的应用

数值流形法在非连续变形分析领域具有独特优势。结合裂纹尖端场函数的基本概念,分析了水力劈裂破坏问题,模拟了水力劈裂破坏过程,避免了扩展有限元中的阶跃函数和水平集概念。为了避免裂纹尖端在单元内部不同位置而产生误差,对裂纹尖端附近一定范围内的每一个物理覆盖附加奇异覆盖函数。选取一个算例比较分析了内水压力对应力强度因子的影响,当考虑裂纹面内压时,定量分析比较了各因素对应力强度因子的影响大小,并应用于分支裂纹水力劈裂破坏。计算结果表明,改进后的计算结果与解析解相吻合。未考虑裂纹面内压,误差往往较大。考虑裂纹面内压后,随着裂纹长度的增加,误差逐渐减小;随着网格密度的增加,误差也逐渐减小。分支裂纹的渐进破坏结果表明该改进方法的可行性,具有较大的实际应用价值。     

岩体黏弹性蠕变计算的高阶数值流形法研究

数值流形法（numerical manifold method）是一种新型的数值计算方法,已成功应用于岩土工程的诸多领域,但该方法尚未应用于岩土工程蠕变分析。近年来对高阶流形法的研究表明,对复杂的岩土工程问题,使用高阶覆盖函数可明显提高流形法的计算精度。为此,开展了用高阶流形法模拟蠕变的研究,在高阶流形法中引入“时步-初应变”法计算蠕变,以广义开尔文体为基础,推导了相关的计算公式,并编制了相应的计算程序,同时还通过算例,验证了方法的可行性和合理性。结果表明,高阶流形可以方便地与“时步-初应变”法结合用于蠕变计算,可较好地模拟蠕变变形。算例分析表明,在不改变网格密度情况下,仅通过采用高阶覆盖函数,高阶流形法可大幅提高传统流形法的计算精度。     

六面体覆盖的高阶数值流形方法的探讨

采用高阶近似位移覆盖函数,基于六面体数学覆盖网格建立了三维数值流形方法分析格式,给出了相应的子矩阵。利用MATLAB编制了与之对应的计算程序,对简单的地下洞室模型进行了计算,并将计算结果与其他数值分析方法结果进行了比较,证明了分析格式及相应程序的正确性和有效性。结果表明, 当数值流形方法的覆盖函数推广到高阶情况时,其求解精度会有相应的提高。最后,对该方法在隧道及地下工程的应用前景作了展望。     

基于数值流形法的重力坝水力劈裂研究

运用数值流形法在非连续介质领域的独特优势,并结合断裂力学的基本原理,分析了重力坝的水力劈裂破坏问题,实现了重力坝起裂、扩展和渐进破坏的全过程。充分发挥应力强度因子判别法和带抗拉强度的摩尔-库仑破坏准则各自的优势,根据不同情况选取其中某种方法确定断裂和扩展方向,且不需要预制初始裂纹。选取一个算例两种工况分析了重力坝在不考虑裂纹面内水压力和考虑裂纹面内水压力的水力劈裂破坏情况。计算结果表明,在重力坝的水力劈裂中,当不考虑裂纹面内水压力时,坝踵处裂纹向下游和深部扩展,上游折坡点处裂纹属于压剪型破坏;当考虑裂纹面内水压力时,坝踵处裂纹偏向深度方向扩展,折坡点处的裂纹逐渐由压剪破坏变成了拉剪破坏,并且考虑裂纹面内水压力后坝体破坏所需要的时步减少,故而考虑裂纹面内水后降低了重力坝的安全系数。该方法加深了对重力坝水力劈裂破坏的认识,具有较大的实际应用价值。     

Malvern Hills边坡溃曲破坏分析及数值流形法模拟

基于能量法推导了多层岩质边坡临界溃曲长度计算公式。考虑层间和交叉节理采用数值流形法对新西兰Malvern Hills边坡溃曲演变过程进行了模拟,并将理论计算、数值模拟结果与现场实测数据进行了对比。结果表明:通过预制层间和交叉节理,数值流形法可以准确模拟边坡溃曲破坏过程;边坡溃曲变形及失稳破坏过程依次为层间错动—轻微弯曲、坡脚牵引—剧烈隆起、加速滑移—形成滑坡;在自重长期作用下,边坡溃曲形成至破坏演化过程主要分为初始弯曲、急剧弯曲和滑坡形成3个阶段;β为交叉节理与斜坡法向夹角,在0°、15°、30°、45°这4种交叉节理中,β为45°时边坡最容易产生滑移弯曲变形,溃曲程度最大,并且经历滑移弯曲时步最少;β为30°～45°时,数值模拟结果与实际较为吻合。     

岩土力学数值分析方法研究

随着岩土工程领域的不断扩展与延伸,岩土力学数值分析方法得到了迅速发展,出现了各种各样的数值分析方法。将岩土力学常用的数值分析方法分为连续变形分析方法与非连续变形分析方法两大类,对各种方法的研究现状和最新进展进行了评述,着重介绍了各种方法的优缺点,并提出了解决问题的思路、方法和建议。