共查询到19条相似文献，搜索用时 111 毫秒
1.
2.
模拟三维裂纹问题的扩展有限元法 总被引：4，自引：1，他引：3
扩展有限元法是一种在常规有限元框架内求解强和弱不连续问题的新型数值方法,其计算网格与不连续面相互独立,因此模拟移动不连续面时无需对网格进行重新剖分。给出了模拟三维裂纹问题的扩展有限元法。在常规有限元位移模式中,基于单位分解的思想加进一个阶跃函数和二维渐近裂尖位移场,反映裂纹处位移的不连续性。用两个水平集函数表示裂纹。采用线性互补法求解裂纹面非线性接触条件,不需要迭代,提高了计算效率。采用两点位移外推法计算裂纹前缘应力强度因子。给出了3个三维弹性静力问题算例,其结果显示了所提方法能获得高精度的应力强度因子,并能有效地处理裂纹面间的接触问题,同时表明扩展有限元结合线性互补法求解不连续问题具有较好的前景。 相似文献
3.
岩体中普遍存在着断层﹑节理和裂隙等结构面,这些结构面的存在和发展对岩体的整体强度﹑变形及稳定性有极大的影响。因此,研究岩体中原生结构面的萌生﹑发展以及贯通演化过程对评估岩体工程安全性和可靠性具有非常重要的理论与现实意义。扩展有限元法（XFEM）作为一种求解不连续问题的有效数值方法,模拟裂隙时独立于网格,因此,在模拟岩体裂隙扩展﹑水力劈裂等方面具有独特优势。针对扩展有限元法的基本理论及其在岩体裂隙扩展模拟中的应用展开了研究,建立了扩展有限元法求解岩体裂隙摩擦接触、岩体裂隙破坏等问题的数值模型,并将计算模型应用于岩质边坡稳定性分析和重力坝坝基断裂破坏等工程问题。 相似文献
4.
扩展有限元法是一种在常规有限元框架内求解强和弱不连续问题的新型数值方法,其原理是在裂尖附近用一些奇异函数和沿裂纹面用阶跃函数加强传统有限元的基,以考虑跨过裂纹的位移场的不连续,该加强策略允许计算网格独立于不连续体几何。讨论了扩展有限元法的一些数值方面,主要包括：水平集法确定界面和加强节点与加强方式、裂尖加强范围的选择、J积分区域的确定和积分方案等。 相似文献
5.
《物探化探计算技术》2015,(3)
无单元Galerkin法作为较成熟的一种无网格方法,已成功应用于有限元法触及的领域,还解决了如大变形、裂纹扩展及高速冲击等网格方法较难处理的问题,但其最大的缺陷在于系统方程的离散需借助背景网格,因此该方法并非真正意义上的无网格方法。无网格局部径向基点插值法采用子域法构造系统方程,加权残量只要求在局部积分域消除,大大降低了对背景网格的依赖,向真正的无网格迈进了一大步.这里将此方法用于大地电磁二维正演,介绍了该方法的基本原理;从大地电磁二维边值问题出发,利用子域法推导了与之对应的无网格局部弱式系统方程,并用高斯积分将其离散化;论述了局部径向基点插值法较无单元Galerkin法及有限元法的优缺点;最后通过二维模型的计算验证了算法的有效性。 相似文献
6.
针对数值流形方法特有的覆盖剖分方式,提出了一种模拟岩土工程中开挖过程的算法。该算法采取某种措施,在覆盖剖分过程中将开挖面视为特殊的不连续面,这种不连续面将其所在的数学网格剖分成不同的流形单元,但却不对所在的数学覆盖作剖分。这样,开挖面两侧虽分属不同的流形单元,但开挖面两侧同一数学网格内的流形单元却具有相同的物理覆盖。采用该算法,无需对开挖面处的单元进行特殊处理,可在整个分析域采用统一的网格形式;同时,打破了原有数值流形方法的限制,将开挖面的位置完全当作连续介质来处理,避免了因将其视为不连续面而产生的误差。验证了算法的可靠性后,将其应用于某假想隧道的开挖模拟,计算结果表明该算法具有一定的应用前景。 相似文献
7.
裂隙岩体因含有发育程度不同的裂隙、节理和断层等不连续面,致其渗透性具有各向异性、不连续性等特点,因此传统的有限元法对分布密集的裂隙岩体渗流场求解有一定的难度。本文提出了采用无单元Glaerkin法求解有自由面裂隙渗流问题,并推导了无单元法求解渗流场的基本方程和积分格式,给出了应用罚函数法处理渗流边界条件和自由面处理方法。采用IDL语言编制了二维无单元法计算软件LIDAREFM。文中以北京怀柔桥梓镇某裂隙岩体边坡渗流场计算为例,研究了复杂裂隙共同作用下渗流场特性和自由面分布,讨论了不同开度、不同连通程度的裂隙对渗流场的影响。研究结果表明:无单元法可以较好地解决有密集裂隙的岩体渗流场的求解问题,实现了裂隙处结点任意加密以及积分网格的独立布置,避免了对有自由面和裂隙穿越的子域的重新处理,简化了渗流问题的求解过程。 相似文献
8.
基于复合单元法,结合三维热传导-对流方程和“充填模型”,提出了裂隙岩体不稳定温度场的复合单元模型。该模型前处理简便快捷,计算网格生成时无需考虑裂隙的存在,网格剖分不受限制,随后利用复合单元前处理程序,依据裂隙的位置和方位将其自动离散在单元内。对常规热传导-对流方程进行自伴随性调整,应用变分原理,推导出裂隙岩体不稳定温度场的复合单元算法,该算法可分别计算出岩块子单元和裂隙的温度值,且可真实反映裂隙中水流与相邻岩块间的热能量交换规律。将复合单元数值模型计算的不稳定温度场结果与相应的实测数据进行对比分析可知,数值计算结果与实测数据基本一致,验证了裂隙岩体不稳定温度场复合单元算法的可靠性与有效性。算例分析表明,裂隙中水流与相邻岩块间有明显的热传导和热对流作用。 相似文献
9.
将无网格局部Petrov-Galerkin算法用于大地电磁二维正演。介绍了该方法的基本原理;从大地电磁二维边值问题出发,利用子域法详细推导了与之对应的局部Petrov-Galerkin弱式方程,并用高斯积分法将其离散化。论述了无网格局部Petrov-Galerkin法较无单元Galerkin法及有限元法的优缺点,最后通过二维模型的计算验证了算法的有效性。 相似文献
10.
运用布尔碎片运算实现了简单的三维块体切割功能，形成覆盖整个求解域的六面体网格，再将六面体单元拆分生成48个四面体单元，从而生成四面体数学覆盖。运用布尔交运算将四面体单元与求解域求交集得到流形块体，再根据三维拓扑有向性原理和三维单纯形积分理论，形成了有向边、有向环、有向面和有向壳4种有向几何数据结构，用来构造生成封闭的有向三维流形单元。定义了有向流形单元的连通内面对和连通的有向流形单元的概念，利用有向流形单元的连通内面对搜索生成物理覆盖体系。概括总结了基于修正对称和反对称分解的三维数值流形元法的求解计算要点，在不考虑三维接触、三维裂纹尖端奇异场和三维裂纹扩展的假设下，模拟了三维节理面的有限塑性变形张拉过程，得到了比较合理的数值模拟结果，验证了前处理和计算求解算法的正确性。 相似文献
11.
岩石裂纹的扩展是一个经典的不连续问题,常规有限元方法难以实现裂纹扩展过程的仿真模拟。扩展有限元法(XFEM)实现了计算网格与不连续面相互独立,因此模拟移动的不连续面时无需对网格进行重新剖分。本文介绍了XFEM基本原理和岩石断裂力学常用判据,尝试对岩石类材料单缝Ⅰ型三点弯曲、单缝剪切和双缝平板实验进行模拟。分析结果表明:扩展有限元模拟岩石类材料断裂问题不受网格划分限制,裂纹以实际应力场分布随机扩展;直观地给出岩样的微裂纹产生、演化,直至完全破坏的全过程,并与实验结果吻合。该方法能够应用到岩石断裂力学方面的研究,模拟岩石类材料的宏细观破坏过程,为解决复杂问题提供了方便的途径。 相似文献
12.
This paper presents a fracture mapping (FM) approach combined with the extended finite element method (XFEM) to simulate coupled deformation and fluid flow in fractured porous media. Specifically, the method accurately represents the impact of discrete fractures on flow and deformation, although the individual fractures are not part of the finite element mesh. A key feature of FM‐XFEM is its ability to model discontinuities in the domain independently of the computational mesh. The proposed FM approach is a continuum‐based approach that is used to model the flow interaction between the porous matrix and existing fractures via a transfer function. Fracture geometry is defined using the level set method. Therefore, in contrast to the discrete fracture flow model, the fracture representation is not meshed along with the computational domain. Consequently, the method is able to determine the influence of fractures on fluid flow within a fractured domain without the complexity of meshing the fractures within the domain. The XFEM component of the scheme addresses the discontinuous displacement field within elements that are intersected by existing fractures. In XFEM, enrichment functions are added to the standard finite element approximation to adequately resolve discontinuous fields within the simulation domain. Numerical tests illustrate the ability of the method to adequately describe the displacement and fluid pressure fields within a fractured domain at significantly less computational expense than explicitly resolving the fracture within the finite element mesh. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
13.
The strength and deformability of rock mass primarily depend on the condition of joints and their spacing and partially on
the engineering properties of rock matrix. Till today, numerical analysis of discontinuities e.g. joint, fault, shear plane
and others is conducted placing an interface element in between two adjacent rock matrix elements. However, the applicability
of interface elements is limited in rock mechanics problems having multiple discontinuities due to its inherent numerical
difficulties often leading to non-convergent solution. Recent developments in extended finite element method (XFEM) having
strong discontinuity imbedded within a regular element provide an opportunity to analyze discrete discontinuities in rock
masses without any numerical difficulties. This concept is based on partition of unity principle and can be used for cohesive
rock joints. This paper summarizes the mathematical frameworks for the implementation of strong discontinuities in 3 and 6
nodded triangular elements and also provides numerical examples of the application of XFEM in one and two dimensional problems
with single and multiple discontinuities. 相似文献
14.
The high computational costs associated with the implicit formulation of discontinuous deformation analysis (DDA) have been one of the major obstacles for its implementation to engineering problems involving jointed rock masses with large numbers of blocks. In this paper, the Newmark-based predictor-corrector solution (NPC) approach was modified to improve the performance of the original DDA solution module in modeling discontinuous problems. The equation of motion for a discrete block system is first established with emphasis on the consideration of contact constraints. A family of modified Newmark-based predictor-corrector integration (MNPC) scheme is then proposed and implemented into a unified analysis framework. Comparisons are made between the proposed approach and the widely used constant acceleration (CA) integration approach and central difference (CD) approach, regarding the stability and numerical damping features for a single-degree-of-freedom model, where the implications of the proposed approach on open-close iteration are also discussed. The validity of the proposed approach is verified by several benchmarking examples, and it is then applied to two typical problems with different numbers of blocks. The results show that the original CA approach in DDA is efficient for the simulation of quasi-static deformation of jointed rock masses, while the proposed MNPC approach leads to improved computational efficiency for dynamic analysis of large-scale jointed rock masses. The MNPC approach therefore provides an additional option for efficient DDA of jointed rock masses. 相似文献
15.
节理岩体中隧道开挖与支护的数值模拟 总被引：11，自引：0，他引：11
由于节理、断层等不连续面的存在造成岩体变形的不连续性,并且这些不连续面对岩体变形、应力等力学行为造成重要影响。对已有的非连续变形分析程序进行了两点改进,应用改进的程序模拟了节理岩体中隧洞开挖与支护的过程,研究了不同节理倾角岩体的应力分布特性、开挖后隧洞围岩的变形和应力分布规律以及支护后衬砌的变形与应力特点。研究结果表明,节理倾角对节理岩体系统的应力和变形有着直接的影响,并且影响着隧洞开挖后的内空收敛效应;隧洞开挖顺序对围岩的变形和稳定也有重要影响。改进后的非连续变形分析方法和程序可以用来分析复杂节理岩体系统中隧洞的开挖与支护。 相似文献
16.
To improve the computational efficiency of the numerical manifold method for discontinuous deformation simulations, a spatial-domain coupled explicit-implicit time integration algorithm is proposed. A subdomain partition algorithm based on a super manifold element is developed for the numerical manifold method to simulate dynamic motions of blocky rock mass. In different subdomains, explicit or implicit time integration method is employed respectively based on its contact and motion status. These subdomains interact through assembling the corresponding explicit or implicit time integration-based matrices of different rock blocks. The computational efficiency of the discontinuity system under dynamic loading is improved by partially diagonalizing the global matrices. Two verification examples of a sliding block along an inclined plane under a horizontal acceleration excitation and a multiblock system acted on by dynamic forces are studied to examine the accuracy of the proposed numerical method, respectively. A highly fractured rock mass situated on an inclined slope subjected to seismic excitations is then studied to show the computational efficiency of the developed algorithm. The simulated results are in good agreement with those from the versions using purely implicit or explicit time integration algorithm for the numerical manifold method. The computational efficiency is shown to be higher using the proposed algorithm, which demonstrates its potential for application in dynamic analysis of highly fractured rock masses. 相似文献
17.
Preexisting flaws and rock heterogeneity have important ramifications on the process of rock fracturing and on rock stability in many applications. Therefore, there is great interest in numerical modelling of rock fracture and the underlying mechanisms. We simulated damage evolution and fracture propagation in sandstone specimens containing a preexisting 3-D surface flaw under uniaxial compression. We applied the linear elastic damage model based on the unified strength theory following the rock failure process analysis code. However, in contrast to the rock failure process analysis code, we used the finite element method with tetrahedron elements on unstructured meshes. It provided higher geometrical flexibility and allowed for a more accurate representation of the disk-shaped flaw with various flaw depths, angles, and lengths through locally adapted meshes. The rock heterogeneity was modelled by sampling the initial local Young's modulus from a Weibull distribution over a cubic grid. The values were then interpolated to the computational finite element method mesh. This method introduced an additional length scale for the rock heterogeneity represented by the cell size in the sampling grid. The generation of three typical surface cracking patterns, called wing cracks, anti-wing cracks, and far-field cracks, were identified in the simulation results. These depend on the geometry of the preexisting surface flaw. The simulated fracture propagation, coalescence types, and failure modes for the specimens with preexisting surface flaw show good agreement with recent experimental studies. 相似文献
18.
In natural rock masses, the shapes of three‐dimensional (3‐D) blocks cut by arbitrary fracture networks may be very complex. Owing to the geometric complexity and difficulty of mesh discretization of 3‐D blocks and fracture facets, explicit consideration of fracture networks in flow analysis of fractured porous medium (FPM) is very challenging. Using the numerical manifold method based on independent covers (NMMIC), an FPM flow model was proposed that can deal with very complex 3‐D fracture networks. In this paper, the convergence of NMMIC was first demonstrated. The theoretical basis of the arbitrary refinement of computational meshes was proven. Moreover, three peculiarities of NMMIC meshes, that is, arbitrary shape, arbitrary connection, and arbitrary refinement of independent covers, were concluded. Finally, some two‐dimensional (2‐D) tunnel flow examples were analyzed and the numerical results were compared with the analytical results. 3‐D examples with complex fracture distributions were also analyzed. In addition, the computational scale of the developed program was tested by increasing the number of computational elements. The results show that our model can accurately analyze the groundwater flow of rocks surrounding tunnels with complex fracture distributions. 相似文献
19.
Savithru Jayasinghe David L. Darmofal Nicholas K. Burgess Marshall C. Galbraith Steven R. Allmaras 《Computational Geosciences》2018,22(1):107-123
This paper presents a space-time adaptive framework for solving porous media flow problems, with specific application to reservoir simulation. A fully unstructured mesh discretization of space and time is used instead of a conventional time-marching approach. A space-time discontinuous Galerkin finite element method is employed to achieve a high-order discretization on the anisotropic, unstructured meshes. Anisotropic mesh adaptation is performed to reduce the error of a specified output of interest, by using a posteriori error estimates from the dual-weighted residual method to drive a metric-based mesh optimization algorithm. The space-time adaptive method is tested on a one-dimensional two-phase flow problem, and is found to be more efficient in terms of computational cost (degrees-of-freedom and total runtime) required to achieve a specified output error level, when compared to a conventional first-order time-marching finite volume method and the space-time discontinuous Galerkin method on structured meshes. 相似文献