四节点四边形数值流形方法及其改进

本文对四节点四边形流形元提出了改进措施,将覆盖位移函数用自然坐标表示,使得在一般非规则有限数学覆盖网格下,数值积分变得比较容易,克服了现有四节点四边形流形单元数值积分困难的缺点。数值算例将其应用于复合材料数值模拟,计算结果表明,当覆盖位移函数采用完全一阶等参多项式时,计算精度较传统有限元法有很大改进。在力应集中或应力突变的区域,无需网格加密,只需提高覆盖位移函数的阶次即可。     

基于余弦样条函数的有限点阵法

应用节点影响域的概念,提出了基于余弦样条函数的有限点阵方法.利用余弦样条函数的性质,通过张量积的形式构造余弦样条函数正则解空间,用于逼近场函数,余弦样条基的线性组合使得边界条件处理如同有限元法一样方便.余弦样条与边界型方法结合,可用于求解不规则域问题.数值实例的计算结果表明,文中方法避免了高阶多项式构造形函数所带来的数值振荡,解的连续性不受限制,为改进计算精度而加密点阵所导致的计算量增加量较少,计算收敛快.     

一种建筑材料细观力学数值模拟的新方法

数值流形方法通过数学和物理双重网格,分析连续和非连续问题,已应用于模拟节理岩体裂隙的开裂与闭合问题.但对于裂纹尖端的局部化问题,数值流形方法需要像有限元那样在裂纹尖端设置细密单元.本文利用裂纹尖端解析解将数值流形方法的基函数进行扩展,推导了相应的试函数.从最小势能原理出发提出了断裂力学的数值流形方法,推导了相应的求解方程,将其应用于建筑材料细观力学数值模拟.最后给出两个数值算例,将计算结果与解析解对比,说明该方法的正确性和可行性.     

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

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

Kirchhoff-Love壳的等几何分析

论文利用等几何分析研究了基于Kirchhoff-Love理论的薄壳的静态问题.等几何分析采用等参思想,将精确描述几何形状的NURBS基函数同时作为场变量的插值函数,保证了在分析和网格优化过程中模型的几何精确性,并可以轻易地构造任意高阶连续的单元.该方法具有很高的数值精度.计算结果表明,在等几何分析中,NURBS单元的阶次越高,网格数越多,计算结果越精确.     

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

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

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

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

基于三角网格的四阶Hermite型对流守恒格式

提出一种基于三角网格的求解双曲对流方程的高阶守恒型格式.该格式首先在每个三角单元上重构二元三次Hermite插值多项式,以当前时刻单元节点处解的函数值、一阶空间导数值和该单元的积分平均值为插值条件.然后,利用Semi-Lagrange方法得到单元节点处的下一时刻解的函数值及导数值,而下一时刻的解的单元积分平均值由有限体积方法得到.本文所提出的格式将原始CIP方法从结构网格推广到非结构网格上,使得CIP方法能灵活地用于处理复杂边界问题.该格式为显式紧致格式,计算简单且易于实现.数值实验表明,该格式对于光滑解问题能达到四阶空间精度,而对于非光滑解问题能准确地捕捉激波的位置,改进了原始CIP格式的不守恒性.     

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

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

高阶数值流形方法的初应力公式

高阶数值流形方法和高阶DDA方法可以显著提高结构变形的计算精度,但目前涉及几何非线性问题的研究成果大都计算精度差甚至不收敛,这是由高阶初应力公式的不准确或不正确引起的。本文介绍数值流形方法的大变形计算格式,基于平面三角形数学网格和多项式覆盖函数,提出高阶流形法的两种初应力处理方法,首次导出了高阶初应力的准确公式。该公式在分步计算的初应力累加中考虑了大变形结构的构形变化,并将初应力表示成多项式函数形式以满足单纯形积分的要求。文中给出的悬臂梁大变形数值算例与理论解的对比结果证明了方法的正确性。本文的方法和公式也适用于三维四面体数学网格,稍加修改后将可应用于高阶DDA方法和常规的有限元方法。     

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

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

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

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

基于NMM的EFG方法及其裂纹扩展模拟

数值流形方法（numerucal manifold method,NMM）通过引入数学覆盖和物理覆盖两套系统来统一处理连续和非连续问题. 通过用移动最小二乘插值（moving least squares interpolation,MLS）中的节点影响域构造数学覆盖,得到了基于数值流形方法的无网格伽辽金法（element free Galerkin,EFG）. 该方法在保证前处理简单的同时,又能方便处理如裂纹等不连续问题. 建立了适用于小变形和大变形的裂纹扩展计算格式,并通过对曲折裂纹（kinked crack）的处理,在不加密的情况下实现了任意小步长的裂纹扩展,大大提高了在固定网格中模拟裂纹扩展的实用性. 大小变形的结果对比表明,按照不考虑构型变化的小变形计算,结果可能偏于危险.      

THE NMM-BASED EFG METHOD AND SIMULATION OF CRACK PROPAGATION

数值流形方法（numerucal manifold method,NMM）通过引入数学覆盖和物理覆盖两套系统来统一处理连续和非连续问题. 通过用移动最小二乘插值（moving least squares interpolation,MLS）中的节点影响域构造数学覆盖,得到了基于数值流形方法的无网格伽辽金法（element free Galerkin,EFG）. 该方法在保证前处理简单的同时,又能方便处理如裂纹等不连续问题. 建立了适用于小变形和大变形的裂纹扩展计算格式,并通过对曲折裂纹（kinked crack）的处理,在不加密的情况下实现了任意小步长的裂纹扩展,大大提高了在固定网格中模拟裂纹扩展的实用性. 大小变形的结果对比表明,按照不考虑构型变化的小变形计算,结果可能偏于危险.     

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

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

基于T样条的变网格等几何薄板动力学分析

具有大位移、大变形的薄板在接触碰撞等工况下, 其局部应变会产生剧烈变化. 为了保证对其进行动力学分析的精度和计算效率, 本文整合计算机辅助设计(CAD)与计算机辅助工程(CAE)系统, 提出了一种基于T样条曲面的变网格柔性系统等几何分析方法. 首先, 建立基于T样条曲面单元的基尔霍夫薄板运动学模型, 并根据非线性格林?拉格朗日应变建立由T样条曲面单元离散的薄板弹性模型. 其次, 通过在T网格中的局部区域插入节点的方式, 达到T样条曲面网格局部更新的目的. 利用T样条混合函数细化算法得到计算新广义变量的转换矩阵, 并结合广义α法创建了变自由度系统动力学方程的求解算法, 形成了系统的T样条单元局部细化算法. 最后, 静力学算例与柔性单摆模型分别验证了T样条薄板弹性模型的正确性, 以及T样条薄板单元在动力学分析上的精度和收敛性. 通过对受冲击柔性薄板的动力学分析表明, 本文所提出T样条单元及局部细化算法可以只在接触碰撞等应变剧烈变化的区域实现局部网格细化, 从而控制系统自由度数, 提高计算效率.      

2D numerical manifold method based on quartic uniform B-spline interpolation and its application in thin plate bending

A new numerical manifold （NMM） method is derived on the basis of quartic uniform B-spline interpolation. The analysis shows that the new interpolation function possesses higher-order continuity and polynomial consistency compared with the conven- tional NMM. The stiffness matrix of the new element is well-conditioned. The proposed method is applied for the numerical example of thin plate bending. Based on the prin- ciple of minimum potential energy, the manifold matrices and equilibrium equation are deduced. Numerical results reveal that the NMM has high interpolation accuracy and rapid convergence for the global cover function and its higher-order partial derivatives.     

基于最佳质量网格的薄板问题的非协调流形方法

对于大部分非协调板单元,使用规则网格能得到很好的效果。但是,当网格不规则时,非协调元的数值特性将变得很差,甚至收敛性得不到保证。为解决网格依赖性问题,许多专家学者提出了改造单元,如拟协调元法和广义协调元法,这些方法能解决收敛性问题,但是数值实践证明没有一种单元能在所有情况下都具有良好的数值特性。考虑到流形方法采用两套完全独立的覆盖系统,可以用规则的数学网格来作为数学覆盖进行插值,取得最佳的插值效果,单元收敛性便能得到保证。再结合适用于流形方法的变分提法,建立起流形方法处理非规则物理边界非协调板单元的一般格式。以ACM薄板单元为例,与ANSYS、拟协调元法和广义协调元法进行了对比,证明本文方法在处理具有曲线边界的薄板弯曲问题时具有收敛快和精度高等优势。