首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 46 毫秒
1.
数值流形方法对岩土工程开挖卸荷问题的模拟   总被引:9,自引:0,他引:9  
基于数值流形方法数学覆盖与物理网格的相对独立性,考虑任一时刻体系中真实存在的各平衡力项对体系平衡的作用,提出用数值流形方法模拟岩土工程开挖卸荷的方法。相对于传统数值方法,数值流形方法在模拟开挖时不用计算开挖面上的释放荷载,不用为卸除的材料部分准备专门的单元,简单的数学网格可以适应任意的开挖过程,使其对开挖的模拟更加简便和有效,同时高阶的流形方法对开挖问题具有较高的求解精度。  相似文献   

2.
朱爱军  邓安福  曾祥勇 《岩土力学》2006,27(Z2):146-150
基于流形方法的有限覆盖技术,确定流形元覆盖系统下锚杆的位移函数,由此建立全长粘结锚杆的数值流形方法模型,定义了锚杆流形单元,并推导了锚杆流形单元的数值计算格式。该模型符合全粘结的物理意义,是一种局部解析的锚杆数值方法,具有较高的精度。同时锚杆与数学覆盖相对独立,只要求锚杆完全处于数学网格的覆盖之中,对其相对位置没有要求,所以一般网格可以适应锚杆的复杂布置。该模型可以应用于岩土锚杆模拟及钢筋与混凝土共同作用分析等问题。  相似文献   

3.
水是影响岩土工程安全性的活跃因素,开挖是岩土工程中常见的工程作用,作为一种较新的数值方法,数值流形方法在这两方面的研究仍然较少,影响了数值流形方法在实际工程中的应用。针对结构面可能承受三种不同类型的水压力作用,分别给出了相应的计算公式以计算水压力对岩体结构的影响。数值方法中常删除开挖区域块体和单元的做法模拟开挖,这种处理方法易于描述且物理意义明确,相应的算例表明将它引入数值流形方法之后,数值流形方法能够合理地模拟开挖后的围岩变形,扩展了数值流形方法的应用范围并为更深入研究流形方法奠定了基础。  相似文献   

4.
柯锦福  王水林 《岩土力学》2020,41(10):3473-3480
运用布尔碎片运算实现了简单的三维块体切割功能,形成覆盖整个求解域的六面体网格,再将六面体单元拆分生成48个四面体单元,从而生成四面体数学覆盖。运用布尔交运算将四面体单元与求解域求交集得到流形块体,再根据三维拓扑有向性原理和三维单纯形积分理论,形成了有向边、有向环、有向面和有向壳4种有向几何数据结构,用来构造生成封闭的有向三维流形单元。定义了有向流形单元的连通内面对和连通的有向流形单元的概念,利用有向流形单元的连通内面对搜索生成物理覆盖体系。概括总结了基于修正对称和反对称分解的三维数值流形元法的求解计算要点,在不考虑三维接触、三维裂纹尖端奇异场和三维裂纹扩展的假设下,模拟了三维节理面的有限塑性变形张拉过程,得到了比较合理的数值模拟结果,验证了前处理和计算求解算法的正确性。  相似文献   

5.
运用布尔碎片运算实现了简单的三维块体切割功能,形成覆盖整个求解域的六面体网格,再将六面体单元拆分生成48个四面体单元,从而生成四面体数学覆盖。运用布尔交运算将四面体单元与求解域求交集得到流形块体,再根据三维拓扑有向性原理和三维单纯形积分理论,形成了有向边、有向环、有向面和有向壳四种有向几何数据结构,用来构造生成封闭的有向三维流形单元。定义了有向流形单元的连通内面对和连通的有向流形单元的概念,利用有向流形单元的连通内面对搜索生成物理覆盖体系。概括总结了基于修正对称和反对称分解的三维数值流形元法的求解计算要点,在不考虑三维接触、三维裂纹尖端奇异场和三维裂纹扩展的假设下,模拟了三维节理面的有限塑性变形张拉过程,得到了比较合理的数值模拟结果,验证了前处理和计算求解算法的正确性。  相似文献   

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

7.
杨石扣  任旭华  张继勋 《岩土力学》2016,37(9):2706-2711
三维流形单元的生成是进行三维数值流形分析的首要问题之一。详细研究了三维流形单元的生成过程,并采用C++语言编写了相应的程序。借鉴二维流形单元的形成技术,基于拓扑学的“有向性”原理,将点、有向边、有向环、有向面和有向壳等作为三维块体的基本数据结构。将材料体和数学网格进行布尔交运算,并对形成的流形块体进行有效性检测,满足要求后即形成新的三维流形单元。每个数学网格的顶点作为新流形单元的数学覆盖,再对数学覆盖进行细分,形成流形单元的物理覆盖。分别选取凹形体、空心体和包含有限结构面的材料体与数学网格进行布尔交运算,并选取一个典型工程来检查该方法和程序的可行性。计算结果表明,该方法可以对复杂块体(凹形体、空心体和包含有限结构面的体)进行处理,为今后进行复杂结构计算和分析奠定基础,具有较强的适应性和可靠性。  相似文献   

8.
提出了一种改进的四边形网格生成算法。该法在行波法生成平面三角形网格的基础上,将一个三角形单元分解为三个四边形单元。经过优化处理,获得优良的计算域四边形网格。针对岩土工程结构问题,还提出了网格剖分中尖灭不连续面和锚杆及锚索的处理方案。地下洞室和边坡算例表明了这种方法的可行性和可靠性。  相似文献   

9.
三维地层的三棱柱剖分与土方计算   总被引:5,自引:0,他引:5  
周翠英  董立国 《岩土力学》2006,27(2):204-208
三维地层剖分算法是实现地层三维可视化与土方计算的重要基础。采用三棱柱体建模方法构造三维地层,分析任意平面切割三棱柱的不同情形,得出了针对三棱柱地层的剖分算法;提出了以切割面为分界面,将单个三棱柱剖分成多个小三棱柱的三角投影分割法--将三棱柱顶面、底面及切割截面按一定规则进行合理三角化划分,使划分后三角形在水平面上的投影相重合;将该3层三角形组合成以切割面为分界面的2层小三棱柱,从而完成剖分。土方计算中首先给出开挖体的边界控制方程,根据该控制方程所限定的内边界得到开挖区域,结合上述剖分算法,对该区域内所有小三棱柱体积求和,得到开挖土方量。将以上算法应用到广东省某重大水利工程中,实现了基坑模拟开挖及其土方计算。  相似文献   

10.
Wilson非协调数值流形方法   总被引:1,自引:0,他引:1  
魏高峰  冯伟 《岩土力学》2006,27(2):189-192
三维数值流形方法中,当数学覆盖取六面体体单元时流形单元总体位移函数中所包含的多项式并不是完全的,非完全的高次项非但对改善精度不起作用,而且还可能起相反的作用。为此,基于Wilson非协调元理论,推导了附加非协调位移基本项的流形元通用公式,通过内参静力凝聚处理,导出了消除单元内参后的单元应变矩阵、单元刚度矩阵,建立了非协调数值流形方法。数值试验表明,在规则数学网格覆盖下它们能够保证收敛,有较高的精度,从而证明所建方法的可行性。  相似文献   

11.
A major challenge in seepage analysis is to locate the phreatic surface in an unconfined aquifer. The phreatic surface is unknown and assumed as a discontinuity separating the seepage domain into dry and wet parts, thus should be determined iteratively with special schemes. In this study, we systematically developed a new numerical manifold method (NMM) model for unconfined seepage analysis. The NMM is a general numerical method for modeling continuous and discontinuous deformation in a unified mathematical form. The novelty of our NMM model is rooted in the NMM two‐cover‐mesh system: the mathematical covers are fixed and the physical covers are adjusted with iterations to account for the discontinuity feature of the phreatic surface. We developed an energy‐work seepage model, which accommodates flexible approaches for boundary conditions and provides a form consistent with that in mechanical analysis with clarified physical meaning of the potential energy. In the framework of this energy‐work seepage model, we proposed a physical concept model (a pipe model) for constructing the penalty function used in the penalty method to uniformly deal with Dirichlet, Neumann, and material boundaries. The new NMM model was applied to study four example problems of unconfined seepage with varying geometric shape, boundary conditions, and material domains. The comparison of our simulation results to those of existing numerical models for these examples indicates that our NMM model can achieve a high accuracy and faster convergence speed with relatively coarse meshes. This NMM seepage model will be a key component of our future coupled hydro‐mechanical NMM model. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

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

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

14.
Summary  Numerical Manifold Method (NMM) is an extension of Discontinuous Deformation Analysis which contains and combines FEM and joint/block oriented DDA in a unified form. At present, application of NMM is mainly limited to two-dimensional analysis which is not realistic for many underground excavation works. In the present paper, three-dimensional NMM based on tetrahedron and hexahedron elements is proposed. The corresponding three-dimensional cover contact detection algorithms which have some fundamental properties different from the corresponding two-dimensional analysis are also developed. Based on the present formulations, engineering problems with arbitrarily oriented discontinuities can be modeled and realistic modeling of construction works becomes possible. Author’s address: Y. M. Cheng, Department of Civil and Structural Engineering, Hong Kong Polytechnic University, Room TV707, Hung Hom Kowloon, Hong Kong, P.R. China  相似文献   

15.
徐栋栋  郑宏 《岩土力学》2014,35(8):2385-2393
数值流形方法(NMM)的最大优势在于可以统一地处理岩土力学中的连续和非连续变形问题。它在求解断裂力学问题时无需强制裂纹与数学网格保持一致,非常适合应用于岩土工程中由连续到非连续的破坏过程模拟。在裂纹扩展过程中,裂纹与数学网格的相对位置将会是任意的,如裂纹尖端可能落在网格内部、网格节点上或网格边上等。因此,对同一条裂纹,通过旋转和移动数学网格构造了它们之间的这种相对位置关系以及一些可能对计算结果产生影响的极端情况,并以应力强度因子作为衡量标准,研究了NMM在处理线弹性断裂力学问题时的网格依赖性。研究表明,NMM即使在处理强奇异性问题时依然有着很好的网格无关性,进一步证实了它在模拟裂纹扩展问题时的鲁棒性。  相似文献   

16.
In this paper, the numerical manifold method (NMM) is extended to study wave propagation across rock masses. First, improvements to the system equations, contact treatment, and boundary conditions of the NMM are performed, where new system equations are derived based on the Newmark assumption of the space–time relationship, the edge‐to‐edge contact treatment is further developed for the NMM to handle stress wave propagation across discontinuities, and the viscous non‐reflection boundary condition is derived based on the energy minimisation principle. After the modification, numerical comparisons between the original and improved NMM are presented. The results show that the original system equations result in artificial numerical damping, which can be overcome by the Newmark system equations. Meanwhile, the original contact scheme suffers some calculation problems when modelling stress wave propagation across a discontinuity, which can be solved by the proposed edge‐to‐edge contact scheme. Subsequently, the influence of the mesh size and time step on the improved NMM for stress wave propagation is studied. Finally, 2D wave propagation is modelled, and the model's results are in good agreement with the analytical solution. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

17.
徐栋栋  杨永涛  郑宏  邬爱清 《岩土力学》2016,37(10):2984-2992
数值流形方法(NMM)基于两套覆盖(数学和物理覆盖)和接触环路而建立,能够统一地处理岩土工程中的连续和非连续变形分析问题。与其他基于单位分解理论的数值方法一样,NMM可以自由地提高物理片上局部位移函数(多项式)的阶次,从而在不加密网格的情况下显著地提高计算精度,但有可能会使总体刚度矩阵奇异,产生线性相关问题。针对这种情况,引入了一种新的高次多项式形式的局部位移函数,在此基础上,建立了新的NMM求解体系,并应用于求解一般的弹性力学问题。结果表明:它有效地消除了线性相关问题;较之传统局部位移函数取一次多项式的NMM,达到了更高的精度;节点应力是连续的;定义在物理片上的所有自由度都具有明确的物理含义,其中第3~5个刚好是物理片所对应插值点处的应变分量,因此,直接获得此处的应力状态。该方法可以很容易地推广到其他基于单位分解的数值方法中。  相似文献   

18.
徐栋栋  孙冠华  郑宏 《岩土力学》2013,34(2):526-534
数值流形法(NMM)在接触处理过程包括:接触方式判断,包括点-点和点-线接触;开-闭迭代,确定块体系统的约束状态;以及接触传递,将有效的接触对传递到下一个时步。过程稍显繁琐和耗时,且与真实的物理接触状态有异。由Munjiza所提出的NBS(no binary search)接触检测算法,将单元映射到规则格子中,以链表结构将其有效地连接在一起,只在单元所在格子以及周围格子内部进行接触判断,接触检测效率大为提高,计算量仅随单元数线性增长,内存需求也很低。在计算接触力时以所定义的势为媒介,用重叠面积来衡量接触力的大小,属于分布式接触力,更接近于实际,避免了原NMM接触处理过程的繁琐。因此,将其作为一种平行的接触处理方法引入到NMM中,并以算例验证,证实了算法的可行性。  相似文献   

19.
The present study focuses on the hydraulic behaviour of joints, and, specially, on its numerical implementation in terms of the FEM analysis using a discrete fracture flow approach. Fluid flow through discontinuities has traditionally been modelled using special elements of zero-thickness, which we can classify into single, double and triple-nodded. Single node elements are the simplest and consist of ‘line’ or ‘pipe’ elements which are superimposed onto the standard continuum mesh and that can only model the longitudinal conductivity through the discontinuity. On the other hand, some authors have included the influence of a transversal conductivity, and the subsequent localized potential drop, by using triple node interface elements. In those, the two nodes of the adjacent continuum elements represent the potentials in the pore system on each side of the interface, and a third node in the middle represents the average potential of the fluid in the channel represented by the discontinuity. Finally, double node interface elements have also been proposed, which have the advantage of making it possible to use the same FE mesh for both mechanical and flow analysis. In some cases the influence of a transversal conductivity is not considered and, therefore, although geometrically double-nodded, these elements belong to the single node type and when time comes to solve the system the two nodes must have the same potential, which can only be obtained by the ‘trick’ of prescribing the equivalence of these two d.o.f. before solving the global system of equations. This limitation may, however, be avoided by assuming that the potential in the channel is the average of the two sides of the interface. Based in this simple assumption, an alternative flow interface model has been recently developed and implemented, which preserves both longitudinal and transversal conductivities. An application example is developed and solved with the three types of interfaces described. The results offer useful information regarding the range of applicability and limitations of the new double-nodded interface element proposed. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

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

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