Biot平面固结分析的流形单元法

基于数值流形方法和有限覆盖技术，提出了适用于Biot固结分析的三节点平面协调流形元。由于土骨架位移和孔隙水压力的节点覆盖函数(Lagrange插值函数)阶次可分别任意选择，该单元是一组满足位移和孔压插值阶次不同且所有节点具有相同自由度数的新型u-p混合模式单元，并且更加方便编程。数值分析表明，位移和孔压的节点覆盖函数阶次分别取一次和零次的流形单元(T1-0)是该组单元中最为有效的。与等价四边形等参元相比，T1-0流形元能给出精度更高的初期孔压和位移。     

弹性力学中的一种非协调数值流形方法

通过引入数学和物理双重网格，将插值域与 积分域分别定义在不同的覆盖上，即在数学网格上进行插值函数的构造，物理网格上完成 系统能量泛函积分运算，最后通过覆盖权函数将二者联结在一起. 它的优点是单元网格划 分随意，不受复杂边界形状和二相材料界面的限制，单元可以是任意形状，是较之于有限 元方法更一般的数值模拟方法. 在4节点四边形数值流形方法中，由于单元总体位移函数 包含的完全多项式不完全，使得计算精度不够精确，为此，在单元总体位移函数上附 加非协调位移基本项，使之趋于完全，提出了弹性力学问题的一种改进的数值流形 方法------非协调数值流形方法. 通过内部自由度静力凝聚处理，导出了消除内参后的单元应变矩阵 和单元刚度矩阵，使得在不增加广义节点自由度的前提下，大大提高了数值流形方法的计 算精度和计算效率. 同时对非协调项进行了显式处理，可以对工程实践起到更切实的帮助. 数值试验表明，它们能够保证收敛，有较高的精度，对畸变不敏感，从而证明了该方法的 可行性.     

基于单位分解法的无网格数值流形方法

在数值流形方法和单位分解法的基础上，提出了无网格数值流形方法. 无网格数值流形 方法在分析时采用了双重覆盖系统，即数学覆盖和物理覆盖. 数学覆盖提供的节点形成求解 域的有限覆盖和单位分解函数；而物理覆盖描述问题的几何区域及其域内不连续性. 与原有 的数值流形方法相比，无网格数值流形方法的数学覆盖形状更加灵活，可以用一系列节点的 影响域来建立数学覆盖和单位分解函数，具有无网格方法的特性，从而摆脱了传统的数值流 形方法中网格所带来的困难. 与无网格方法相比，由于采用了有限覆盖技术，试函数的构造 不受域内不连续的影响，克服了原有的无网格方法在处理不连续问题时所遇到的困难. 详细推导了无网格数值流形方法的试函数和求解方程，最后给出了算例，验证了该方法的正 确性.     

热传导问题的非协调数值流形方法

数值流形方法通过引入数学与物理双重网格，将插值域与积分域分别定义在两个不同的覆盖上，其优点是网格划分随意，不受复杂边界形状和材料界面的限制，是较之于有限元方法更一般化的数值模拟方法。在计算精度方面，数值流形方法远远高于有限元法。但它的精度还是不够理想。为此本文在单元总体位移场上附加非协调位移基本项，使单元位移函数趋于完全，构造了非协调流形单元来改善流形单元的计算精度和计算效率，并将其应用于热传导问题，推导了势问题的非协调数值流形方法。     

高精度四节点四边形流形单元

基于数值流形方法的四点四边形流形单元在物理覆盖上使用高阶覆盖位移函数能够得到高精度的数值结果,且可以在求解区域的不同地方混合使用各阶覆盖函数来提高求解效率,具有编程 和前后处理简单等优点,弥补了有限元的不足,计算结果表明,数值解与理论解吻合。     

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

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

有限覆盖径向点插值方法理论及其应用

数值流形方法能够统一地处理连续与非连续变形问题,有限覆盖技术是这种方法的核心。无网格方法前处理过程比较简单,径向点插值法是其中的一种计算格式。本文将有限覆盖技术与径向点插值方法相结合发展了有限覆盖径向点插值无网格方法,综合了数值流形方法与点插值方法的各自优点,能够有效地处理连续与非连续性问题,由此所构造的形函数具有Kronecker δ-函数属性,能够有效地处理位移边界条件。本文在阐述了这种方法基本原理的基础上,通过算例分析与数值计算论证了本文所建议方法的可靠性及其有效性。     

配点型广义节点无网格法的基本原理及其应用

借鉴流形方法思想,引入广义节点的概念,对传统的无网格法进行了改进,建立了可具有任意高阶多项式插值函数的广义节点无网格方法。同时采用径向插值函数构造具有插值特性的逼近函数;采用配点法建立系统的离散方程。在阐述了这种方法基本原理的同时,针对线弹性力学问题给出了这种方法的数值计算列式。与传统无网格方法相比,这种方法更具有一般性;同时由于采用了配点法而不需要背景积分网格,所以可以认为这种方法是某种真正意义上的无网格法。当选取0阶广义节点位移插值函数时便可得到传统的无网格法;在不增加支持域内节点数目的条件下,通过选取高阶广义节点位移插值函数可以提高计算精度。最后通过算例分析,对0阶、1阶及2阶广义节点无网格法与现有的有关解答进行了对比,论证了其合理性。     

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

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

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

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

A novel four-node quadrilateral element with continuous nodal stress

Formulation and numerical evaluation of a novel four-node quadrilateral element with continuous nodal stress（Q4-CNS）are presented.Q4-CNS can be regarded as an improved hybrid FE-meshless four-node quadrilateral element（FE-LSPIM QUAD4）, which is a hybrid FE-meshless method.Derivatives of Q4-CNS are continuous at nodes, so the continuous nodal stress can be obtained without any smoothing operation.It is found that,compared with the standard four-node quadrilateral element（QUAD4）,Q4- CNS can achieve significantly better accuracy and higher convergence rate.It is also found that Q4-CNS exhibits high tolerance to mesh distortion.Moreover,since derivatives of Q4-CNS shape functions are continuous at nodes,Q4-CNS is potentially useful for the problem of bending plate and shell models.     

基于重叠划分的自由网格四边形单元计算方法

提出了一种基于重叠划分的自由网格四边形单元计算方法。这一方法将四边形单元引入到自由网格计算方法中，不仅提高了计算的精度，同时还保留了自由网格计算方法的特点。方法首先对分析域内自动生成的每一个节点建立一套临时三角形单元，利用这些临时三角形单元组合生成四边形单元，以节点为单位进行计算。由于各矩阵的计算与组集均以节点为中心进行处理，因而特别适合于并行计算环境。在详细介绍自由网格四边形单元计算方法的基础上，利用数值算例证实了这一方法改善计算精度方面的有效性。     

平面广义四节点等参元GQ4及其性能探讨

广义节点有限元是将传统有限元方法中的节点广义化,在不增加节点个数的前提下,仅通过提高广义节点的插值函数的阶次,从而达到提高有限元解精度的目的.与现有的p型和hp型有限元不同,在这种新的有限元中,节点自由度全部定义在节点处,在理论与程序实现上与传统有限元方法具有很好的相容性,传统有限元方法是这种新方法的广义节点退化为0阶时的特殊情形.文中主要讨论了这一新方法的四节点等参元（记为GQ4）的形式.对GQ4进行的各种数值试验表明,所发展的广义四节点等参单元具有精度高且无剪切自锁与体积自锁等的特点.     

黏弹性人工边界等效荷载计算的改进方法

黏弹性人工边界在场地地震反应和结构-地基动力相互作用等问题的计算中已得到了广泛的应用.地震波在黏弹性人工边界中的输入是通过将地震波转化为作用于人工边界处的等效载荷来实现的.计算等效节点载荷的常规方法默认边界节点对应区域内的应力为均布力,但实际上该节点对应区域内的应力分布通常是不均匀的.本文在有限元方法结合黏弹性局部人工边界的显式时域波动方法的基础上,建立了无限域散射问题地震波等效载荷计算的一种改进方法.该方法采用细化网格与应力积分相结合的方法计算人工边界等效节点力,有效地降低了人工边界上等效节点力的计算误差.以不同角度入射地震波的二维算例为例,算例给出的波场位移云图和节点位移时程曲线验证了本文方法的有效性,其计算精度与网格尺寸和地震波入射角度密切相关,且网格越小、入射角度越小,计算精度越高.对于相同的网格尺寸,本文采用方法的计算精度明显高于常规方法,尤其是对于斜入射问题优势更为明显.     

平面弹性力学问题的离散元法

根据离散元的基本原理,基于变形体的理论提出了适用于平面弹性力学问题的界面位移、应变和应力模式,建立了求解平面弹性力学问题的离散元方程和相应的迭代求解方法.通过界面位移可以简洁地将位移和力的边界条件引入离散系统的控制方程,也可以方便地求解节点位移.数值算例表明,与具有相同网格的有限元结果相比,离散元能同时给出精度相对较高的应力解和精度相当的位移解.     

裂纹问题的一致性高阶无网格法

一致性高阶无网格法能高效精确地求解连续体问题，尤其是能得到高精度的应力场。本文将该方法拓展到应力解析精度至关重要的裂纹问题（即非连续体问题）的数值分析。采用背景积分网格描述裂纹几何，基于无需增加节点额外自由度的虚拟节点法描述裂纹处位移场的间断，提出了虚拟节点的引入算法和断裂单元的数值积分方法。为进一步模拟裂纹扩展，采用相互作用积分方法计算应力强度因子，裂纹的扩展方向由最大周向应力准则确定。数值结果表明，本文发展方法能够精确地通过间断分片试验；相较于标准的高阶无网格法和低阶一致性无网格法，本文的一致性高阶无网格法显著改善了应力强度因子的计算精度，能够准确预测裂纹扩展路径。     

Numerical solutions of two-dimensional anisotropic crack problems

A complete set of series form solutions of stress and displacement functions, including all higher order terms, around the crack tip for anisotropic crack problems have been newly derived by eigenfunction expansion approach. The analytical solutions of displacement functions were classified into four cases with respect to different types of complex parameters and different corresponding physical meanings. By employing these displacement functions as global interpolation functions, fractal two-level finite element method (F2LFEM) was applied to evaluate the stress intensity factors (SIFs) for various kinds of anisotropic crack problems. In the method of F2LFEM, the infinite number of nodal displacements was transformed to a small set of generalized coordinates by fractal transformation technique. New element matrices need not be generated and the singular numerical integration was avoided completely. Numerical examples of the four cases were studied and high accurate results of SIFs were obtained.