首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到16条相似文献,搜索用时 125 毫秒
1.
数值流形方法及其在岩石力学中的应用   总被引:9,自引:0,他引:9  
李树忱  程玉民 《力学进展》2004,34(4):446-454
数值流形方法是目前岩石力学分析的主要方法之一.该方法起源于不连续变形分析,主要用于统一求解连续和非连续问题,其核心技术是在分析时采用了双重网格:数学网格提供的节点形成求解域的有限覆盖和权函数;而物理网格为求解的积分域.数学网格被用来建立数学覆盖,数学覆盖与物理网格的交集定义为物理覆盖,由物理覆盖的交集形成流形单元.流形方法的优点在于它使用了独立的数学和物理网格,具有和有限元明显不同的定义形式,且数学网格对于同一问题不同的求解精度的需求可以很方便地细化.由于该方法考虑了块体运动学,可以模拟节理岩体裂隙的开裂和闭合过程,因而在岩石力学中得到了广泛应用,近年来许多学者对该方法进行了研究.本文简要叙述了节理岩体的数值方法从连续到非连续的发展过程,详细地介绍了数值流形方法的组成和数值流形方法在岩石力学及其相关领域的研究和发展概况,最后就作者所关心的一些问题,如三维问题的数值流形方法、数值流形方法在物理非线性问题和裂纹扩展问题中的应用、相关的耦合方法等进行了探讨.  相似文献   

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

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

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

5.
基于流形覆盖思想的无网格方法的研究   总被引:20,自引:3,他引:17  
本语言基于流形思想,利用有限覆盖,单位分解等概念,引入建立在覆盖上的覆盖函数和具有紧支撑特性的单位分解函数,建立场逼近的近似表达,由弱形式的Galerkin变分得到数值分析模型,结合边界条件用于边值问题的求解,由此建立了一类新的无网格数值方法,论文采用这种方法分析了平面弹性问题,分析了体积闭锁现象,h、p型收敛性等,提出了一种选择覆盖大小的方案,且对狭长城采用了椭圆覆盖形式,取得了比较好的效果。  相似文献   

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

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

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

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

10.
复杂裂纹问题的多边形数值流形方法求解   总被引:1,自引:0,他引:1  
数值流形方法是一种能统一处理连续和不连续问题的有效数值方法。该方法采用的数学覆盖系统可完全独立于物理域,能很好地求解各类裂纹问题,而n边形单元(n>4)则具有网格划分灵活,求解精度高等优点。本文基于数值流形方法,采用正六边形数学单元求解线弹性复杂裂纹问题。在导出相关方程的基础上对典型裂纹问题进行了分析,通过互能积分法得到了裂尖的应力强度因子,计算结果与参考解吻合得较好。除此之外,文中还对不同单元上的求解精度进行了比较,结果表明采用正六边形单元的求解精度较正四边形单元和正三角形单元上的精度均更高。  相似文献   

11.
The meshless manifold method is based on the partition of unity method and the finite cover approximation theory which provides a unified framework for solving problems dealing with both continuum with and without discontinuities. The meshless manifold method employs two cover systems. The mathematical cover system provides the nodes for forming finite covers of the solution domain and the partition of unity functions. And the physical cover system describes geometry of the domain and the discontinuous surfaces in the domain. The shape functions are derived by the partition of unity and the finite covers approximation theory. In meshless manifold method, the mathematical finite cover approximation theory is used to model cracks that lead to interior discontinuities in the displacement. Therefore, the discontinuity is treated mathematically instead of empirically by the existing methods. However, one cover of a node is divided into two irregular sub-covers when the meshless manifold method is used to model the discontinuity. As a result, the method sometimes causes numerical errors at the tip of a crack. To improve the precision of the meshless manifold method, the enriched methods are introduced in this work for crack problems.  相似文献   

12.
单位分解扩展无网格法(PUEM)是一种求解不连续问题的新型无网格方法.其基于单位分解思想,通过在传统无网格法的近似函数中加入扩展项来反映由裂纹所产生的不连续位移场.详细描述了水平集方法,PUEM不连续近似函数的构造及控制方程的离散.针对裂纹扩展问题,提出了一种十分简单的水平集更新算法;讨论了不同的节点数、高斯积分阶次以及围线积分区域对应力强度因子计算结果的影响,并给出了合理的参数;模拟了边裂纹和中心裂纹的扩展问题,并与XFEM的数值结果进行了比较.数值算例表明,本文方法具有较高的计算精度,是模拟裂纹扩展非常有效的方法,具有广阔的应用前景.  相似文献   

13.
平面裂纹问题的h, p, hp型自适应无网格方法的研究   总被引:17,自引:0,他引:17  
刘欣  朱德懋  陆明万  张雄 《力学学报》2000,32(3):308-318
无网格方法以其独特的优点:不需“网格”(即节点间的连接信息)划分,特别适合自适应的分析,在分析中只需要高梯度域简单地插入离散点(h型)或保持模型节点数、分布、覆盖大小均不变,中增加高误差覆盖上的函数的多项式阶次(p型),便可以得到更高精度的数值模型。针对平面弹性问题发展和推导一种显式后验误差指示公式,对平面裂纹实例进行了h型,p型,hp型三种不同类型的无网格自适应分析,数值分析结果表明了这种自适应  相似文献   

14.
The incompatible numerical manifold method (INMM) is based on the finite cover approximation theory, which provides a unified framework for problems dealing with continuum and discontinuities. The incompatible numerical manifold method employs two cover systems as follows. The mathematical cover system provides the nodes for forming finite covers of the solution domain and the weighted functions, and the physical cover system describes geometry of the domain and the discontinuous surfaces therein. In INMM, the mathematical finite cover approximation theory is used to model cracks that lead to interior discontinuities in the process of displacement. Therefore, the discontinuity is treated mathematically instead of empirically by the existing methods. However, one cover of a node is divided into two irregular sub-covers when the INMM is used to model the discontinuity. As a result, the method sometimes causes numerical errors at the tip of a crack. To improve the precision of the INMM, the analytical solution is used at the tip of a crack, and thus the cover displacement functions are extended with higher precision and computational efficiency. Some numerical examples are given.  相似文献   

15.
扩展有限元法(XFEM)及其应用   总被引:46,自引:3,他引:43  
扩展有限元法(extended finite element method, XFEM)是1999年提出的一种求解不连续力学问题的数值方法, 它继承了常规有限元法(CFEM) 的所有优点, 在模拟界面、裂纹生长、复杂流体等不连续问题时特别有效, 短短几年间得到 了快速发展与应用. XFEM与CFEM的最根本区别在于, 它所使用的网格与结构内部的几何或 物理界面无关, 从而克服了在诸如裂纹尖端等高应力和变形集中区进行高密度网格剖分所带 来的困难, 模拟裂纹生长时也无需对网格进行重新剖分. 重点介绍XFEM的基本原理、 实施步骤及应用实例等, 并进行必要的评述. 单位分解概念保证了XFEM的收敛, 基于此, XFEM 通过改进单元的形状函数使之包含问题不连续性的基本成分, 从而放松对网格密度的过分要 求. 水平集法是XFEM中常用的确定内部界面位置和跟踪其生长的数值技术, 任何内部界面 可用它的零水平集函数表示. 第2和第3节分别简要介绍单位分解法和水平集法; 第4节和第5节介绍XFEM的基本思想、详细实施步骤和若干应用实例, 同时修正了以往文 献中的一些不妥之处; 最后, 初步展望了该领域尚需进一步研究的课题.  相似文献   

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

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号