共查询到20条相似文献，搜索用时 31 毫秒
1.
2.
A state-of-the-art review is conducted to highlight the fracture mechanism in rock blast and advantages and limitations of various methods in modelling it. A hybrid finite-discrete element method (FEM-DEM) is implemented to simulate rock fracture and resultant fragment muck-piling in various blasting scenarios. The modelled crushed, cracked and long radial crack zones are compared with those in literatures to calibrate the hybrid FEM-DEM. Moreover, the hybrid modelling reproduces the rock fragmentation process during blasting. It is concluded that the hybrid FEM-DEM is superior to continuous and discontinuous methods in terms of modelling dynamic fracture of rock under blast-induced impact load. 相似文献
3.
岩体中普遍存在着断层﹑节理和裂隙等结构面,这些结构面的存在和发展对岩体的整体强度﹑变形及稳定性有极大的影响。因此,研究岩体中原生结构面的萌生﹑发展以及贯通演化过程对评估岩体工程安全性和可靠性具有非常重要的理论与现实意义。扩展有限元法（XFEM）作为一种求解不连续问题的有效数值方法,模拟裂隙时独立于网格,因此,在模拟岩体裂隙扩展﹑水力劈裂等方面具有独特优势。针对扩展有限元法的基本理论及其在岩体裂隙扩展模拟中的应用展开了研究,建立了扩展有限元法求解岩体裂隙摩擦接触、岩体裂隙破坏等问题的数值模型,并将计算模型应用于岩质边坡稳定性分析和重力坝坝基断裂破坏等工程问题。 相似文献
4.
The paper deals with two-dimensional (2D) numerical modelling of hydro-fracking (hydraulic fracturing) in rocks at the meso-scale. A numerical model was developed to characterize the properties of fluid-driven fractures in rocks by combining the discrete element method (DEM) with computational fluid dynamics (CFD). The mechanical behaviour of the rock matrix was simulated with DEM and the behaviour of the fracturing fluid flow in newly developed and pre-existing fractures with CFD. The changes in the void geometry in the rock matrix were taken into account. The initial 2D hydro-fracking simulation tests were carried out for a rock segment under biaxial compression with one injection slot in order to validate the numerical model. The qualitative effect of several parameters on the propagation of a hydraulic fracture was studied: initial porosity of the rock matrix, dynamic viscosity of the fracking fluid, rock strength and pre-existing fracture. The characteristic features of a fractured rock mass due to a high-pressure injection of fluid were realistically modelled by the proposed coupled approach.
相似文献5.
运用ANSYS程序对水下的岩质边坡进行数值模拟分析,采用不同的材料强度来模拟岩石和软弱结构面。应用强度折减方法对边坡的强度参数进行折减,直至边坡出现破坏,分析边坡在达到破坏时内部的应力应变情况和位移变形情况。对于波浪力,本文采用简化的波压力方程模拟低频的水波对边坡的波压力影响,并用多个模型比较计算的方法得出最佳的模型尺寸和水下荷载的计算模式。 相似文献
6.
This paper proposes a three-dimensional coupled hydrothermal model for fractured rock based on the finite-discrete element method to simulate fluid flow and heat transport. The 3D coupled hydrothermal model is composed of three main parts: a heat conduction model for the rock matrix, a heat transfer model for the fluid in the fractures (including heat conduction and heat convection), and a heat exchange model between the rock matrix and the fluid in the fractures. Four examples with analytical solutions are provided to verify the model. A heat exchange experiment of circulating water in a cylindrical granite sample with one fracture is simulated. The simulation results agree well with the experimental results. The effects of the fracture aperture, fluid viscosity, and pressure difference on the heat exchange between the fluid and rock are studied. Finally, an application concerned with heat transport and fluid flow in fractured rock is presented. The simulation results indicate that the 3D fully coupled hydrothermal model can capture the fluid flow and temperature evolution of rocks and fluids. 相似文献
7.
水工隧洞内水外渗耦合分析 总被引：1，自引：0，他引：1
在水荷载体力理论的基础上,引入钢筋应变不均匀系数和混凝土应变不均匀系数。计算衬砌渗透系数,引入反映钢筋混凝土衬砌与围岩有条件联合承载特性的充水夹层单元模型,采用等效耦合分析方法,对水工隧洞内水外渗条件下的衬砌钢筋应力和围岩应力状态进行了计算分析。计算结果表明,与考虑衬砌与围岩完全联合承载的耦合分析方法相比,钢筋应力有了较大幅度的降低,在一定程度上反映了钢筋混凝土衬砌高压隧洞运行过程中钢筋应力水平普遍较低的这一工程特性。 相似文献
8.
9.
Timo Saksala 《国际地质力学数值与分析法杂志》2013,37(3):309-324
This paper considers numerical modelling of rock fracture induced by dynamic bit–rock interaction in percussive drilling. The work presented here extends the author's earlier research on the topic from the axisymmetric case to 3D case. The numerical method for modelling rock fracture includes a constitutive model for rock and a contact mechanics‐based technique to simulate the bit–rock interaction. The constitutive model is based on a combination of the recent viscoplastic consistency model, the isotropic damage concept and a parabolic compression cap. This model is improved here from its earlier state by calibrating the softening laws using fracture energies GIc and GIIc in tension and compression, respectively. Moreover, the viscosity modulus in tension is calibrated based on the dynamic Brazilian disc test. With these enhancements, the developed method is applied to 3D case of the bit–rock interaction problem assuming one symmetry plane. Single impact with single and multiple‐button bits is simulated. In the latter case, an initial borehole is modelled in order to simulate the usual in‐situ drilling conditions. The different failure types observed in the experiments as well as the interaction between the buttons resulting in chipping are realistically captured in the simulations. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
10.
块碎石作为道砟材料能够分散列车动荷载, 加之其内部含有大量空隙, 具有良好的对流换热特性, 作为冷却冻土路基的材料在青藏铁路建设中被广泛使用。但在列车荷载的作用下, 碎石集料会发生压密、 变形乃至破碎从而影响冷却效果。因此, 研究碎石层的变形过程及机理具有十分重要的现实意义。基于块体离散元法, 对块碎石集料的三轴剪切试验进行了数值分析, 将所得应力-应变曲线与室内试验结果进行了对比, 发现两者能较好地吻合, 说明块体离散元法能够较好地模拟块碎石集料的受力变形过程。结果表明： 增大围压或块体粒径, 块体单元受到的作用力加强, 集料的偏应力强度和抗剪强度值也随之增大。碎石块体在剪切作用下沿其接触面滑动分离, 形成X形剪切带是集料变形的主要形式, 此外在径向方向出现不同程度的扩张。基于试验和块体元研究路基碎石层的思路和方法可为今后评估青藏铁路碎石路基的热力稳定性提供理论依据和参数储备。 相似文献
11.
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. 相似文献
12.
A numerical model is described for coupled flow and mechanical deformation in fractured rock. The mechanical response of rock joints to changes in hydraulic pressure is strongly influenced by the geometric characteristics of the joint surfaces. The concept of this work is to combine straightforward finite element solutions with complex and realistic fracture surface geometry in order to reproduce the non-linear stress-deformation-permeability coupling that is commonly observed in fractures. Building on the numerous studies that have expanded the understanding of the key parameters needed to describe natural rough-walled fractures, new methods have been developed to generate a finite element mesh representing discrete fractures with realistic rough surface geometries embedded in a rock matrix. The finite element code GeoSys/Rockflow was then used to simulate the coupled effects of hydraulic stress, mechanical stress, and surface geometry on the evolving permeability of a single discrete fracture. The modeling concept was experimentally verified against examples from the literature. Modeling results were also compared to a simple interpenetration model. 相似文献
13.
Three dimensional numerical modelling was used to examine the fracture responses around cavities in rock masses experiencing the stress of excavation. In addition to the primary fracture zone in the near-field, numerical modelling generated a second fracture zone in the far-field and an elastic non-fracture zone between the two fields, i.e., fracture and non-fracture zones occurred alternately around a deep cavity. Further research illustrated that the dynamic load and static stress gradient are two necessary precursors for a far-field fracture in the excavation process. Neither quasi-static loading nor homogeneous stress conditions could induce a far-field fracture. A simple theory is introduced, suggesting that multiple fracture zones occur during excavation due to both the initial stress gradient and the dynamic load. This finding indicates that it may be possible to induce continuous rock fractures in deep underground rock masses by employing optimal excavation methods to generate multiple contiguous fracture zones. 相似文献
14.
不同条件下水压裂隙的发展特性对有效开采页岩气具有重要的指导作用。针对岩体在微观上为颗粒和孔隙的结构系统,提出离散元水力压裂数值模拟方法,离散元能量转化和能量守恒计算方法,建立了相应的三维离散元模型。采用自主研发的三维离散元模拟软件Mat DEM3D,通过控制模型的竖向应变与颗粒直径,来模拟地层中的应力与压裂速率的变化。模拟结果表明:(1)水力压裂产生裂隙的数量和方向受岩石的各向异性,压力状态和变化速率所影响。(2)裂隙在压缩波传播时发展,当水压力高速增加时,诱发的裂隙数量增多,并且有效能量(断裂热)百分比也随之增加,压裂作用也变得更明显。(3)当竖向应变为零时,50%的裂隙呈垂直状态,当竖向应变为-1×10-4时,裂隙趋于沿着最大压力方向发展,竖向裂隙的百分率增大。数值模拟和能量分析为定量地研究岩石水力压裂过程提供了一个新的方法。 相似文献
15.
Numerical modelling was applied to study the present-day state of stress and deformation under different tectonic loading conditions at the seismic BALTIC–SKJ profile in south-eastern Finland and in Estonia. The finite element method was used to solve the numerical problem. The two-dimensional model was constructed using the results from both seismic and thermal studies along the profile. The model is 700 km long and 200 km deep, and is roughly divided into an inhomogeneous, laterally layered crust and a homogeneous mantle lithosphere. Both the linear elastic and non-linear elasto-plastic rheologies were used. Elasto-plasticity was achieved by calculating a rheological strength as a function of depth along the profile. Different tectonic load cases were analysed with displacement, force and pressure type boundary conditions. Also, the effect of different strain rates was investigated. The results suggest that even with relatively low compressive stress levels the lower crust deforms in a plastic manner for a wet crustal rheology. When applying a dry crustal rheology, plastic yielding is attained only with much higher stress fields. 相似文献
16.
模拟岩石破裂过程的块体单元离散弹簧模型 总被引：2，自引：0，他引：2
在变形体离散元的基础上建立块体单元离散弹簧模型,并应用于岩石破裂过程的数值模拟研究。该模型以连续介质力学理论为基础,将块体单元离散为具有明确物理意义的弹簧系统,通过对弹簧系统的能量泛函求变分获得各弹簧的刚度系数,进而可以直接利用弹簧刚度求解单元的变形和应力,提高计算效率。以重力作用下的岩质边坡计算为例,通过与传统的有限元进行对比,验证该模型弹性计算结果的正确性。在该基础上,引入Mohr-Coulomb与最大拉应力的复合破坏准则,判断单元的破坏状态及破裂方向。当单元的内部破坏面确定后,则通过块体切割的方式实现单元破坏,并建立单元边界和单元内部的双重破裂机制,实现块体由连续到非连续的破裂过程,进而显示的模拟裂纹的形成和扩展。最后,以巴西圆盘劈裂、单轴压缩破裂以及三点弯曲梁等典型算例验证该方法,结果表明该方法可以较好地模拟拉伸、压剪等应力状态下裂纹的形成和扩展,从而可模拟岩石介质由连续到非连续的破裂过程。 相似文献
17.
加锚节理岩体有着复杂的细部结构和力学行为,常规的数值模型很难满足模拟这种力学行为的要求。进一步发展了笔者建立的加锚岩体的复合单元法,提出加锚节理岩体的复合单元法研究。在不考虑锚杆、节理等细部结构的情况下划分网格,分别定义了加锚岩体复合单元、不连续岩体复合单元和加锚节理岩体复合单元,并建立了其数值模型。同时,在复合单元中还定义了不同材料介质的子单元用来模拟加锚节理岩体内复杂的细部构造。对加锚节理复合单元法进行了理论推导和算例研究,通过与常规有限单元法计算结果的比较,表明了加锚节理岩体复合单元法的合理性,同时也突出体现了其前处理简便的优势。 相似文献
18.
The purpose of this paper is to investigate the application of the discrete element program UDEC for modelling underground explosions in a jointed rock mass. A field underground explosion test has been conducted and UDEC and AUTODYN were utilised in a coupled manner to predict the rock mass response due to the underground explosion. The dynamic equation for the state of the rock material obtained from impact tests was incorporated into the calibration modelling. The comparison of modelling results with test results and empirical formulas shows that UDEC is capable of modelling explosion wave propagation in jointed rock mass with high reliability. 相似文献
19.
R. Resende L. N. Lamas J. V. Lemos R. Calçada 《Rock Mechanics and Rock Engineering》2010,43(6):741-761
The goal of this paper is to simulate the interaction of stress waves and rock fractures in a particle micromechanical model. Stress waves travelling in fractured rock masses are slowed down and attenuated by natural heterogeneities, voids, microcracks and, above all, by faults and fractures. Considerable laboratory and theoretical investigation have uncovered the major aspects of this phenomenon, but models that cover the core mechanisms of the wave propagation in rock masses are necessary to investigate aspects of wave–fracture interaction, which are not completely clear, and in the future simulate full-scale real problems. The micromechanical model is based on the particle discrete element model that reproduces rock through a densely packed non-structured assembly of 2D disks with point contacts. The model of a hard rock core is developed and an irregular rock joint is generated at mid-height. A new contact constitutive model is applied to the particles in the joint walls. Numerical static joint compression tests are performed and a typical hyperbolic stress–displacement curve is obtained. Conditions for good quality wave transmission through non-jointed unorganized particulate media are determined, hybrid static–dynamic boundary conditions are established and plane waves are emitted into the compressed joint. The transmitted and reflected waves are extracted and analysed. Joint dynamic stiffness calculated according to the hypotheses of the Displacement Discontinuity Theory shows to increase with the static joint compression until the joint is completely closed. Still in its early stages of application, this rock micromechanical model enables the joint behaviour under static and dynamic loading to be analysed in detail. Its advantages are the reproduction of the real mechanics of contact creation, evolution and destruction and the possibility of visualizing in detail the joint geometry changes, which is hard to accomplish in the laboratory. 相似文献
20.
基于完整锚杆动测技术的围岩质量识别研究 总被引：1，自引：1，他引：0
完整锚杆系统是由完整的锚杆杆体、饱满的锚固介质以及围岩组成的有机体系。在描述锚杆低应变动力响应的数学力学模型中,将锚固介质与围岩对锚杆体的作用简化为沿杆长均匀分布的杆侧弹性系数与阻尼系数,以及均匀分布在杆底的刚度系数与阻尼系数,即锚杆系统的动力参数。它们不仅表征了锚杆锚固质量,而且体现出锚杆所处围岩的物理参数。通过建立围岩物理参数与完整锚杆系统动力参数的映射关系,可以对围岩质量进行识别。为此,建立了完整锚杆三维轴对称有限元计算模型,得到不同的围岩物理力学参数所对应的锚杆杆顶速度响应;根据完整锚杆低应变动力响应理论解和遗传算法建立了完整锚杆动力参数的识别方法;通过对不同围岩下完整锚杆的数值模拟结果的动力参数识别,建立了杆侧弹性系数与围岩变形模量之间的关系式,提出了一种基于完整锚杆动测技术的围岩质量识别方法。 相似文献