共查询到20条相似文献，搜索用时 31 毫秒
1.
Particle manifold method (PMM) is a new extension of the numerical manifold method (NMM). PMM uses a mathematical cover system to describe the motion and deformation of a particle‐based physical domain. By introducing the concept of particle into NMM, PMM takes the advantages of easy topological and contact operations with particles. In this article, the methodology, formulations and implementation of the method are presented, together with modelling examples for validation. It is found that good solutions for both continuous and discontinuous problems are obtained by the new developed PMM. Due to the underlying coupled continuum‐discontinuum property of PMM, it has great potential for modelling of geomechanical problems. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
2.
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. 相似文献
3.
裂隙岩体热-水-应力（THM）耦合是目前研究的热点和难点。首先总结了裂隙岩体多场耦合的机制、模型、方法及研究内容,并通过分析裂隙对THM耦合的重要控制作用,提出了在THM耦合中考虑裂隙网络扩展演化及模拟的关键问题,同时指出了研究的3个关键点：（1）建立考虑裂隙网络演化的耦合模型;（2）裂隙扩展的数值模拟方法;（3）THM耦合及岩体变形、失稳全过程的数值模拟算法。随后通过对模拟多场耦合和裂隙扩展数值方法的归类比较,重点论述了目前适用于模拟多场耦合下裂隙扩展模拟的各种数值方法（包括有限单元法、无单元法、单位分解法、离散单元法、岩石破裂过程分析方法和数值流形方法）的优缺点,并通过对比研究,推荐采用数值流形方法（NMM）来实现对关键问题的模拟研究。最后,对研究思路和难点进行了初步探讨。 相似文献
4.
Coupled hydro-mechanical (HM) processes are significant in geological engineering such as oil and gas extraction, geothermal energy, nuclear waste disposal and for the safety assessment of dam foundations and rock slopes, where the geological media usually consist of fractured rock masses. In this study, we developed a model for the analysis of coupled hydro-mechanical processes in porous rock containing dominant fractures, by using the numerical manifold method (NMM). In the current model, the fractures are regarded as different material domains from surrounding rock, i.e., finite-thickness fracture zones as porous media. Compared with the rock matrix, these fractured porous media are characterized with nonlinear behavior of hydraulic and mechanical properties, involving not only direct (poroelastic) coupling but also indirect (property change) coupling. By combining the potential energy associated with mechanical responses, fluid flow and solid–fluid interactions, a new formulation for direct HM coupling in porous media is established. For indirect coupling associated with fracture opening/closure, we developed a new approach implicitly considering the nonlinear properties by directly assembling the corresponding strain energy. Compared with traditional methods with approximation of the nonlinear constitutive equations, this new formulation achieves a more accurate representation of the nonlinear behavior. We implemented the new model for coupled HM analysis in NMM, which has fixed mathematical grid and accurate integration, and developed a new computer code. We tested the code for direct coupling on two classical poroelastic problems with coarse mesh and compared the results with the analytical solutions, achieving excellent agreement, respectively. Finally, we tested for indirect coupling on models with a single dominant fracture and obtained reasonable results. The current poroelastic NNM model with a continuous finite-thickness fracture zone will be further developed considering thin fractures in a discontinuous approach for a comprehensive model for HM analysis in fractured porous rock masses. 相似文献
5.
6.
在岩土工程分析中求解精度控制常常是必需的,在数值流形法中可以通过控制数学覆盖网格的稀疏和覆盖位移的阶数来达到精度的要求。提出了基于等几何分析的数值流形方法,定义了相应的数学覆盖的构造形式,推导了基于二次B样条的9节点数值流形方法分析格式;针对基于Lagrange插值函数的4节点数值流形方法提出了基于T样条思想的数学覆盖网格的局部加密方法。算例计算结果表明,相对于4节点的数值流形方法,基于非均匀有理B样条的9节点数值流形方法具有更高的精度;基于T样条思想的加密网格在保持计算精度的前提下降低了自由度的数量,表明T样条加密是一种自然的局部加密算法。 相似文献
7.
It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
8.
9.
基于数值流形法（NMM）进行三维裂纹扩展分析研究,编写了相应的C++程序。充分发挥数值流形法在非连续变形分析领域的优势,不需要扩展有限元中的水平集和阶跃函数概念,应用数值流形法计算结果分析裂纹尖端线的破坏状况,对已有的非局部求迹法和三角形推进法进行简化和扩展,提出一种简化算法确定最终的新裂纹扩展面。应用简化算法对水平钱币型裂纹扩展问题进行数值模拟和对比分析研究。计算结果表明,基于数值流形法的三维裂纹扩展算法是可行的,采用简化处理之后,克服原有算法的不足,大大减少了新裂纹尖端线和新裂纹面的数目,降低了计算网格对新裂纹面的影响,提高了计算效率和扩大了应用范围。 相似文献
10.
运用三维数值流形法（3D NMM）进行三维裂纹扩展分析,并采用C++语言编写了相应的程序。充分利用三维数值流形法模拟裂纹扩展的优势,只需要更新裂纹尖端线附近的边界环路和流形单元,不需要使用阶跃函数。根据三维数值流形法计算得到的应力结果,应用非局部求迹方法分析每个裂纹尖端的破坏状态,如果发生破坏则沿垂直于其最大主应力方向扩展。针对裂纹扩展后的不同状况,采用四边形或三角形推进法。裂纹扩展后为了使变形后的面保持平面,必须对新生成的面进行三角化分割。对诸如单边裂纹、平行钱币型裂纹和倾斜钱币型裂纹扩展问题进行数值模拟。计算结果表明,采用三维数值流形法进行裂纹扩展模拟是可行的,文中方法对裂纹尖端线非闭合和闭合的情形均适用,且文中方法对于裂纹尖端线位于单元内部的非平面裂纹扩展也是有效的。 相似文献
11.
Development of a discontinuous approach for modeling fluid flow in heterogeneous media using the numerical manifold method 
下载免费PDF全文
下载免费PDF全文 In the numerical modeling of fluid flow in heterogeneous geological media, large material contrasts associated with complexly intersected material interfaces are challenging, not only related to mesh discretization but also for the accurate realization of the corresponding boundary constraints. To address these challenges, we developed a discontinuous approach for modeling fluid flow in heterogeneous media using the numerical manifold method (NMM) and the Lagrange multiplier method (LMM) for modeling boundary constraints. The advantages of NMM include meshing efficiency with fixed mathematical grids (covers), the convenience of increasing the approximation precision, and the high integration precision provided by simplex integration. In this discontinuous approach, the elements intersected by material interfaces are divided into different elements and linked together using the LMM. We derive and compare different forms of LMMs and arrive at a new LMM that is efficient in terms of not requiring additional Lagrange multiplier topology, yet stringently derived by physical principles, and accurate in numerical performance. To demonstrate the accuracy and efficiency of the NMM with the developed LMM for boundary constraints, we simulate a number of verification and demonstration examples, involving a Dirichlet boundary condition and dense and intersected material interfaces. Last, we applied the developed model for modeling fluid flow in heterogeneous media with several material zones containing a fault and an opening. We show that the developed discontinuous approach is very suitable for modeling fluid flow in strongly heterogeneous media with good accuracy for large material contrasts, complex Dirichlet boundary conditions, or complexly intersected material interfaces. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
12.
The concurrent multiscale method, which couples the discrete element method (DEM) for predicting the local micro‐scale evolution of the soil particle skeleton with the finite element method (FEM) for estimating the remaining macro‐scale continuum deformation, is a versatile tool for modeling the failure process of soil masses. This paper presents the separate edge coupling method, which is degenerated from the generalized bridging domain method and is good at eliminating spurious reflections that are induced by coupling models of different scales, to capture the granular behavior in the domain of interest and to coarsen the mesh to save computational cost in the remaining domain. Cundall non‐viscous damping was used as numerical damping to dissipate the kinetic energy for simulating static failure problems. The proposed coupled DEM–FEM scheme was adopted to model the wave propagation in a 1D steel bar, a soil slope because of the effect of a shallow foundation and a plane‐strain cone penetration test (CPT). The numerical results show that the separate edge coupling method is effective when it is adopted for a problem with Cundall non‐viscous damping; it qualitatively reproduces the failure process of the soil masses and is consistent with the full micro‐scale discrete element model. Stress discontinuity is found in the coupling domain. Copyright © 2017 John Wiley & Sons, Ltd. 相似文献
13.
为了保障岩土工程结构能长期正常使用,需要对其蠕变变形进行分析。“时步-初应变”法是一种常用的计算岩石蠕变的方法。数值流形方法是一种新兴的数值计算方法,常用于计算节理岩体的变形,但尚未被试用于计算蠕变变形。在原数值流形方法的程序中增加了基于“时步-初应变”法的计算模块,通过对广义开尔文模型进行的模拟,显示新程序可以正确反映岩石的黏弹性蠕变趋势,并能够计算包含节理的岩体的蠕变变形,改进后的数值流形方法不但能够模拟岩石的线弹性变形,而且可以模拟岩石的黏弹性蠕变,比原流形方法更能全面地模拟岩石的变形,扩展了数值流形方法在岩土工程中的使用范围。 相似文献
14.
This paper presents a second-order work analysis in application to geotechnical problems by using a novel effective multiscale approach. To abandon complicated equations involved in conventional phenomenological models, this multiscale approach employs a micromechanically-based formulation, in which only four parameters are involved. The multiscale approach makes it possible a coupling of the finite element method (FEM) and the micromechanically-based model. The FEM is used to solve the boundary value problem (BVP) while the micromechanically-based model is utilized at the Gauss point of the FEM. Then, the multiscale approach is used to simulate a three-dimensional triaxial test and a plain-strain footing. On the basis of the simulations, material instabilities are analyzed at both mesoscale and global scale. The second-order work criterion is then used to analyze the numerical results. It opens a road to interpret and understand the micromechanisms hiding behind the occurrence of failure in geotechnical issues. 相似文献
15.
16.
A simple three-dimensional heat transfer model is developed to consider the hindering effect of cracks on heat transfer. The 3D heat transfer model can also be applied to numerical methods such as the combined finite-discrete element method (FDEM), discrete element method (DEM), discontinuous deformation analysis (DDA), the numerical manifold method (NMM), and the finite element method (FEM) to construct thermo-mechanical coupling models that allow these methods to solve thermal cracking problems and dynamically consider the hindering effect of cracks on heat transfer. In the 3D heat transfer model, the continuous-discontinuous medium is discretized into independent tetrahedral elements, and joint elements are inserted between adjacent tetrahedral elements. Heat transfer calculations for continuous-discontinuous media are converted to heat conduction in tetrahedral elements and the heat exchange between the adjacent tetrahedral elements through the joint element. If the joint element between adjacent tetrahedral elements breaks (ie, a crack generates), the heat exchange coefficient of the joint element is reduced to account for the hindering effect of cracks on heat conduction. Then the model and the FDEM are combined to build a thermo-mechanical coupling model to simulate thermal cracking. The thermally induced deformation, stress, and cracking are investigated by the thermo-mechanical coupling model, and the numerical results are compared with analytical solutions or experimental results. The 3D heat transfer model and thermo-mechanical model can provide a powerful tool for simulating heat transfer and thermal cracking in a continuous-discontinuous medium. 相似文献
17.
《Geomechanics and Geoengineering》2013,8(2):147-153
With high-order numerical manifold method (NMM) or high-order discontinuous deformation analysis method (DDA), computational accuracy of structure deformation can be improved greatly. However, poor accuracy is obtained and even computation is not convergent while treating large deformation problems, due to inaccurate or incorrect high-order initial stress formulae. Based on 2D triangular mathematical meshes and polynomial cover functions in high-order NMM, exact formulae for high-order initial stresses are deduced to depict configuration change of structures under large deformations. The formulae are expressed in polynomial forms so as to be used in simplex integrations. The approach is also extended to high-order DDA. Comparing with analytical solutions, accurate results for large deformation of a cantilever beam prove the validity of the proposed formulae. 相似文献
18.
We pay a revisit to some classical geomechanics problems using a novel computational multiscale modelling approach. The multiscale approach employs a hierarchical coupling of the finite element method (FEM) and the discrete element method. It solves a boundary value problem at the continuum scale by FEM and derives the material point response from the discrete element method simulation attached to each Gauss point of the FEM mesh. The multiscale modelling framework not only helps successfully bypass phenomenological constitutive assumptions as required in conventional modelling approaches but also facilitates effective cross‐scale interpretation and understanding of soil behaviour. We examine the classical retaining wall and footing problems by this method and demonstrate that the simulating results can be well validated and verified by their analytical solutions. Furthermore, the study sheds novel multiscale insights into these classical problems and offers a new tool for geotechnical engineers to design and analyse geotechnical applications based directly upon particle‐level information of soils. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
19.
针对数值流形方法特有的覆盖剖分方式,提出了一种模拟岩土工程中开挖过程的算法。该算法采取某种措施,在覆盖剖分过程中将开挖面视为特殊的不连续面,这种不连续面将其所在的数学网格剖分成不同的流形单元,但却不对所在的数学覆盖作剖分。这样,开挖面两侧虽分属不同的流形单元,但开挖面两侧同一数学网格内的流形单元却具有相同的物理覆盖。采用该算法,无需对开挖面处的单元进行特殊处理,可在整个分析域采用统一的网格形式;同时,打破了原有数值流形方法的限制,将开挖面的位置完全当作连续介质来处理,避免了因将其视为不连续面而产生的误差。验证了算法的可靠性后,将其应用于某假想隧道的开挖模拟,计算结果表明该算法具有一定的应用前景。 相似文献
20.
弹塑性数值流形方法在边坡稳定分析中的应用 总被引：2，自引：0，他引：2
针对现有数值流形方法只能进行弹性计算的不足,建立了一个能够反映完整岩块弹塑性变形特征的本构模型,并借助VC++开发了内置该本构模型的弹塑性数值流形程序。利用该程序模拟了含单节理岩样的室内压缩试验,分析得到了其强度和变形特性,计算结果符合实际的物理现象,表明程序是正确有效的。考虑到数值流形方法本身能够有效模拟材料的不连续变形,新增的弹塑性分析功能又可以反映岩石的强度特性,将弹塑性数值流形程序应用于某含有不连续面的岩石边坡的稳定性分析。并结合锚杆单元的使用,对比分析了不同锚固方案的加固效果。程序提供的变形、应力等计算结果表明:预应力锚杆不仅可以防止不连续面发生剪切破坏,增强坡体的稳定性,限制塑性变形的发展;而且可以使不连续软弱层面对岩体变形的消极影响得以减弱,起到提高岩体整体性的作用。 相似文献