共查询到20条相似文献，搜索用时 46 毫秒
1.
将由室内试验总结得到的岩石微观胶结模型嵌入离散元软件,对Lac du Bonnet花岗岩石进行预制单裂隙单轴压缩试验DEM数值模拟,分析了压缩过程中裂隙试样中应力的分布,并与理论计算结果进行对比分析,同时对各种断裂判据中裂纹起裂角的预测值进行了适用性的对比分析。结果表明,离散元模拟试样破坏形态与试验结果相近;离散元分析得到的应力分布与理论解在定性上相似;当预制角度较小时,侧向应力都处于拉压状态;由于裂隙左右两端压应变的集中造成了裂隙上下面拉应变的产生,造成了裂隙周围特殊的应力分布;当裂隙角度较大时,应力集中现象已不明显,因而,理论值与试验值有偏差;在断裂判据中最大周应力准则和最大能量释放率准则得到的裂纹扩展角与室内试验与DEM结果中的数值较为吻合。 相似文献
2.
A contact model for rock is established and imbedded into a DEM software by summarizing the bond granule tests. DEM simulation of uniaxial compression test on the pre-cracked Lac du Bonnet granite is performed, and then stress distributions are further analyzed and compared with the theoretical results. Different fracture criteria are employed to predict the crack initiation angles that are compared with theoretical ones. The results show that the failure modes obtained from DEM simulation are similar to experimental results, and stress distributions in DEM simulation are qualitatively similar to theoretical values. When the angle of pre-crack is small, the lateral stresses are compressive and tensile. The compressive strains concentrate at two edges, resulting in the tensile strains in the up-and downward cracks. When the angle of the pre-crack is large enough, the stress concentration is unobvious, leading to a discrepancy between the DEM and theoretical results. The crack extension angle resulting from uniaxial compression measured from DEM tests are in good agreement with those acquired from experimental tests. These angles are consistent with theoretical predictions by the maximum circumferential stress criterion and the maximum energy release rate criterion. 相似文献
3.
Lekan Olatayo Afolagboye Jianming He Sijing Wang 《Geotechnical and Geological Engineering》2018,36(1):105-133
Crack initiation and coalescence behavior of rock or rock-like specimens containing artificial flaws under uniaxial compression have been subjects of intensive investigation in the past. Most of these investigations however focused on crack initiation and coalescence between two or more parallel flaws. Although there have been few experimental studies on non-parallel flaws, these studies did not address the influence of geometrical factors such as ligament length and ligament angle on the crack initiation and coalescence behavior of non-parallel flaws. To investigate whether the individual geometrical factors have similar effects on the crack initiation and coalescence behavior of both parallel and non-parallel flaws, we conducted uniaxial compression tests to investigate crack cracking and coalescence processes in rock like material containing two non-parallel flaws. The paper presents the influence of individual geometrical factors on the crack initiation process and coalescence pattern of non-parallel flaws. Initiation of primary first cracks from all the tips of the two flaws did not occur simultaneously in all the flaw configurations. The flaw configuration of the non-parallel flaws influences the crack initiation, crack trajectories and coalescence behavior. The crack coalescence pattern changes with an increasing ligament angle from indirect to shear crack or mixed tensile-shear crack to tensile crack coalescence. The chance of direct coalescence is reduced with an increase in ligament length. In conclusion, the crack initiation and coalescence behavior of prismatic rock-like specimens with non-parallel flaws, as influenced by the geometrical factors, are analogous to the cracking and coalescence pattern observed in specimens with parallel flaws. 相似文献
4.
Strength,fragmentation and fractal properties of mixed flaws 总被引：1，自引：1，他引：0
Experiments on Portland cement samples containing mixed flaws are conducted to investigate the strength, fragmentation and fractal properties. Flaw geometry is a new combination of two edge-notched flaws and an imbedded flaw, which is different from those in the previous studies, where parallel or coplanar flaws are used. The physical implications of the shear-box test applied to result to rock slopes are studied. The physical and analytical fragmentation characteristics of preflawed samples are analyzed through the sieve test and fractal theory, respectively. Three different patterns of tensile cracks and shear cracks are observed. A sliding crack model is presented to elucidate the brittle failure flaws. In all of the cases of the shear-box tests, the coalescence is produced by the linkage of shear cracks, and two types of coalescence (Type C1 and Type C2) have been classified, which tend to confirm the observations from the numerical model and field of jointed rock slopes. The shear strength is a function of the flaw geometry and the shear–normal stress ratio. The result of sieve tests indicates that the fragment size distribution of fragments has the fractal property, providing a physical understanding of the fragmentation mechanism. The fragments under the shear-box test have fractal dimensions between 2.2 and 2.6, which are larger than those under the compression test but similar to those in the fault cores. The fragmentation in the case of Type C2 has a smaller fractal dimension, corresponding to a larger shear strength. 相似文献
5.
An Experimental Study of Crack Coalescence Behaviour in Rock-Like Materials Containing Multiple Flaws Under Uniaxial Compression 总被引：4，自引：1，他引：3
Experiments on man-made flawed rock-like materials are applied extensively to study the mechanical behaviour of rock masses as well as crack initiation modes and crack coalescence types. A large number of experiments on specimens containing two or three pre-existing flaws were previously conducted. In the present work, experiments on rock-like materials (formed from a mixture of sand, plaster, limestone and water at mass ratio of 126:9:9:16) containing multiple flaws subjected to uniaxial compression were conducted to further research the effects of the layout of pre-existing flaws on mechanical properties, crack initiation modes and crack coalescence types. Compared with previous experiments in which only three types of cracks were found, the present experiments on specimens containing multiple flaws under uniaxial compression revealed five types of cracks, including wing cracks, quasi-coplanar secondary cracks, oblique secondary cracks, out-of-plane tensile cracks and out-of-plane shear cracks. Ten types of crack coalescence occurred through linkage among wing cracks, quasi-coplanar secondary cracks, oblique secondary cracks, out-of-plane shear cracks and out-of-plane tensile cracks. Moreover, the effects of the non-overlapping length and flaw angle on the complete stress–strain curves, the stress of crack initiation, the peak strength, the peak strain and the elastic modulus were also investigated in detail. 相似文献
6.
Cracking and coalescence behavior in a rectangular rock-like specimen containing two parallel (stepped and coplanar) pre-existing open flaws under uniaxial compression load has been numerically studied by a parallel bonded-particle model, which is a type of bonded-particle model. Crack initiation and propagation from two flaws replicate most of the phenomena observed in prior physical experiments, such as the type (tensile/shear) and the initiation stress of the first crack, as well as the coalescence pattern. Eight crack coalescence categories representing different crack types and trajectories are identified. New coalescence categories namely “New 1” and “New 2”, which are first observed in the present simulation, are incorporated into categories 3 and 4, and category 5 previously proposed by the MIT Rock Mechanics Research Group, respectively. The flaw inclination angle (β), the ligament length (L) (spacing between two flaws) and the bridging angle (α) (inclination of a line linking up the inner flaw tips, between two flaws) have different effects on the coalescence patterns, coalescence stresses (before, at or post the peak stress) as well as peak strength of specimens. Some insights on the coalescence processes, such as the initiation of cracks in the intact part of specimens at a distance away from the flaw tips, and coalescence due to the development and linkage of a number of steeply inclined to vertical macro-tensile cracks are revealed by the present numerical study. 相似文献
7.
In this paper, a coupled constitutive model is proposed for anisotropic damage and permeability variation in brittle rocks under deviatoric compressive stresses. The formulation of the model is based on experimental evidences and main physical mechanisms involved in the scale of microcracks are taken into account. The proposed model is expressed in the macroscopic framework and can be easily implemented for engineering application. The macroscopic free enthalpy of cracked solid is first determined by approximating crack distribution by a second‐order damage tensor. The effective elastic properties of damaged material are then derived from the free enthalpy function. The damage evolution is related to the crack growth in multiple orientations. A pragmatic approach inspired from fracture mechanics is used for the formulation of the crack propagation criterion. Compressive stress induced crack opening is taken into account and leads to macroscopic volumetric dilatancy and permeability variation. The overall permeability tensor of cracked material is determined using a micro–macro averaging procedure. Darcy's law is used for fluid flow at the macroscopic scale whereas laminar flow is assumed at the microcrack scale. Hydraulic connectivity of cracks increases with crack growth. The proposed model is applied to the Lac du Bonnet granite. Generally, good agreement is observed between numerical simulations and experimental data. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
8.
The failure process of brittle rock submitted to a compression state of stress with different confining pressures is investigated in this paper based on discrete element method (DEM) simulations. In the DEM model, the rock sample is represented by bonding rigid particles at their contact points. The numerical model is first calibrated by comparing the macroscopic response with the macroscopic response of Beishan granite obtained from laboratory tests. After the validation of numerical model in terms of macroscopic responses, the failure process of the DEM model under unconfined and confined compression is studied in micro‐scale in detail. The contact force network and its relation to the development of micro‐cracks and evolution of major fractures are studied. Confining pressure will prohibit the development of tensile cracks and hence alter the failure patterns. An in‐depth analysis of micro‐scale response is carried out, including the orientation distribution and probability density of stress acting on parallel bonds, the effect of particle size heterogeneity on bond breakage and the evolution of fabric tensor and coordination number of parallel bond. The proposed micromechanical analysis will allow us to extract innovative features emerged from the stresses and crack evolution in brittle rock failure process. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
9.
The phase field model represents sharp cracks by diffusive mushy-zone and can simulate crack propagation automatically. Propagation and coalescence of quasi-static cracks in Brazilian disks are investigated by a phase field model. The phase field modeling is implemented in Comsol Multiphysics and initially verified by a benchmark of three-point bending test. The Brazilian disk specimens containing no initial crack, a single and two pre-existing cracks subjected to compression are then tested by the phase field model. Crack propagation patterns along with the load–displacement curves are fully discussed. Meanwhile, the effects of length scale parameter and critical energy release rate on crack propagation are evaluated. In addition, the effect of crack inclination angle on the pre-cracked Brazilian disk specimens is also investigated. The numerical results obtained by the phase field model are in good agreement with previous experimental and numerical results.
相似文献10.
Crack evolution is initiated by the occurrence of tensile wing cracks and is then further promoted due to the crack coalescence caused by the extension of a central tensile crack segment between two relatively adjacent flaws. To understand such progressive failures in rock, a parallelized peridynamics coupled with a finite element method is utilized. Through this method, the initiation position of tensile wing cracks is observed with respect to varying inclination angles of a flaw, and then its corresponding shifting mechanism is investigated. In addition, the phenomenon of the position shifting being sensitive to various flaw shapes is discussed. Moreover, it is observed that the inclination angle of a central flaw affects the initiation position of other flaws; therefore, the initiation positions of tensile wing crack emanating from other neighboring flaws are analyzed with their angles. Following tensile wing cracks, a central tensile crack segment occurs in the bridging region between a central flaw and other neighboring flaws; the developmental patterns caused by the crack segment are discussed as well. Finally, the role a central tensile crack segment plays in the formation of crack coalescence and specimen failure is investigated in detail. The numerical results in this paper demonstrate good fidelity with established physical test results and complement them, thereby expanding the understanding of fracturing morphology in rock specimens with various flaws. 相似文献
11.
高水压下岩体裂纹扩展的渗流-断裂耦合机制与数值实现 总被引：1，自引：0，他引：1
采用渗流力学、断裂力学理论结合Monte Carlo方法描述岩体裂纹的随机分布,研究高水压作用下岩体原生裂纹的变形和翼形裂纹的萌生、扩展、贯通的渗流-断裂耦合作用机制,建立高水压作用下岩体裂纹的渗流-断裂耦合数学模型,给出该数学模型的求解策略与方法,在Fortran95平台下开发高水压下岩体裂纹扩展的渗流-断裂耦合分析程序HWFSC.for。高水压下岩体裂纹扩展的渗流-断裂耦合体现在岩体裂纹网络和渗流初始条件都随渗流时步变化。对高压注水岩体裂纹扩展过程进行渗流-断裂耦合分析。结果表明,高压注水条件下,岩体裂纹扩展存在起动水压力,当水压力大于起动水压力时,裂纹尖端开始萌生翼形裂纹,随着裂纹水压力的增加,翼形裂纹扩展,进而与其他裂纹搭接贯通,停止扩展。渗流-断裂耦合分析考虑了裂纹动、静水压力对裂纹产生的法向扩张效应及翼形裂纹的扩展而形成新的渗流通道两方面的影响,连通裂纹数随渗流的发展而增加。岩体裂纹的渗流-断裂耦合分析,能较真实地再现岩体裂纹的水力劈裂现象,描述岩体裂纹的扩展、贯通过程及与之相耦合的渗流响应。 相似文献
12.
含双裂隙岩石裂纹演化机理的离散元数值分析 总被引：2，自引：0，他引：2
采用离散单元法探讨了预制双裂隙岩石的裂纹演化机理。用近期从试验资料提取的无胶结厚度含抗转动能力的岩石微观力学模型和相应的离散单元法商业软件,模拟了含不同预制倾角的双裂隙岩石试样在单轴压缩作用下裂纹的扩展与贯通规律,揭示了裂纹演化的宏微观机理。同时,将离散元法DEM岩石试样的裂纹的扩展和贯通规律以及强度特性与室内试验结果进行了比较分析。结果表明,预制裂隙之间以及端点处的拉应力集中是导致裂隙岩石破坏的主要原因,且DEM数值试验得到裂纹的演化规律与室内试验结果较为一致。含30°的预制裂隙的岩石试样最容易起裂,含75°的预制裂隙的岩石试样最困难起裂,造成此种现象的原因可能是裂纹在垂直于主应力方向上的长度不同导致试样受拉区域大小不同。 相似文献
13.
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.
相似文献14.
微裂隙对于岩石的物理性质有重要影响.通过光学显微镜对显微裂隙进行的观察,难以将薄片的制作过程中产生的裂隙同自然裂隙区分;其次,观测的裂隙大小限制在0.1mm或以上数量级,一些重要细节观察不到.扫描电镜的高放大倍数和三维分辨率使其成为适合于裂隙研究的理想工具.利用扫描电镜观测了浙江花岗岩在室温下由51.6MPa压力产生的裂隙的发育过程.观察花岗岩的表面以研究其微裂隙和矿物的解理、晶形及破裂作用.微裂隙分为3种类型:晶体内裂隙(完全发育在颗粒内部),晶体间裂隙(穿过颗粒边界进入其他颗粒中),颗粒边界裂隙(沿颗粒边界发育或与边界重合).本研究中的花岗岩为中—粗粒(1～6mm),经过蚀变,约含40%斜长石,25%钾长石,25%石英,10%镁铁矿物(主要为黑云母和白云母).(1)石英裂隙表面的石英颗粒具有边界破裂和贯穿颗粒破裂.贯穿石英颗粒的破裂常呈贝壳状断口.变形集中带与滑动面相伴,其上发育有擦痕和沟槽,这是石英典型的脆性破裂类型.石英矿物的贝壳状破裂产状指示其变形作用早期的碎裂作用,在变形作用中发生交代作用.一些平滑的破裂面上没有石英颗粒的微裂隙构造,但在石英表面却有明显的擦痕面,可指示应力的方向.由于石英的各向异性,这些面也是由贯穿颗粒破裂所形成.一些开放的裂隙集中于石英晶体之中或沿 相似文献
15.
微裂隙对于岩石的物理性质有重要影响.通过光学显微镜对显微裂隙进行的观察,难以将薄片的制作过程中产生的裂隙同自然裂隙区分;其次,观测的裂隙大小限制在0.1mm或以上数量级,一些重要细节观察不到.扫描电镜的高放大倍数和三维分辨率使其成为适合于裂隙研究的理想工具.利用扫描电镜观测了浙江花岗岩在室温下由51.6MPa压力产生的裂隙的发育过程.观察花岗岩的表面以研究其微裂隙和矿物的解理、晶形及破裂作用.微裂隙分为3种类型:晶体内裂隙(完全发育在颗粒内部),晶体间裂隙(穿过颗粒边界进入其他颗粒中),颗粒边界裂隙(沿颗粒边界发育或与边界重合).本研究中的花岗岩为中-粗粒(1~6mm),经过蚀变,约含40%斜长石,25%钾长石,25%石英,10%镁铁矿物(主要为黑云母和白云母). 相似文献
16.
Summary ?Cracks that initiate from pre-existing discontinuities can link with other cracks or with other discontinuities and produce
failure in a rock mass. The Displacement Discontinuity Method (DDM), FROCK, is used in this investigation to model experimental
observations on pre-cracked specimens of gypsum. In these experiments two fractures, which were either both open or closed,
were placed through the thickness of the specimens, and detailed observations of the cracking process were performed as the
specimens were loaded in uniaxial compression. The following aspects are studied for both open and closed fractures: 1) crack
initiation stress; 2) direction and propagation of the new cracks; 3) type of coalescence and stress at which it occurs. Modeling
is done considering the actual size of the specimens. Relations between the direction of initiation for each type of crack,
the orientation of the initial fractures, and the type and coalescence are established. In addition, comparisons between results
from experiments and predictions from the model are presented. The numerical results are in agreement with the experiments. 相似文献
17.
Peng-Zhi Pan Jonny Rutqvist Xia-Ting Feng Fei Yan Quan Jiang 《Rock Mechanics and Rock Engineering》2014,47(6):2183-2198
We present a formulation of a discontinuous cellular automaton method for modeling of rock fluid pressure induced fracture propagation and coalescence without the need for remeshing. Using this method, modelers discretize a discontinuous rock-mass domain into a system composed of cell elements in which the numerical grid and crack geometry are independent of each other. The level set method, which defines the relationship between cracks and the numerical grid, is used for tracking the crack location and its propagation path. As a result, no explicit meshing for crack surfaces and no remeshing for crack growth are needed. Discontinuous displacement functions, i.e., the Heaviside functions for crack surfaces and asymptotic crack-tip displacement fields, are introduced to represent complex discontinuities. When two cracks intersect, the tip enrichment of the approaching crack is annihilated and is replaced by a Heaviside enrichment. We use the “partition of unity” concept to improve the integral precision for elements, including crack surfaces and crack tips. From this, we develop a cellular automaton updating rule to calculate the stress field induced by fluid pressure. Then, the stress is substituted into a mixed-mode fracture criterion. The cracking direction is determined from the stress analysis around the crack tips, where fracture fluid is assumed to penetrate into the newly developed crack, leading to a continuous crack propagation. Finally, we performed verification against independent numerical models and analytic solutions and conducted a number of simulations with different crack geometries and crack arrangements to show the robustness and applicability of this method. 相似文献
18.
E.G. Bombolakis 《Tectonophysics》1973,18(3-4)
Static friction along inclined cracks in photoelastic models increases with both the loading and frictional displacement before the larger crack-wall asperities are broken or completely over-ridden. Elastic shocks resulting from successive stick-slip during this stage can be repeated more or less indefinitely with repeated loading and unloading. Locked-in residual stresses, especially around crack tips, result when the model is unloaded, because of frictional coupling between crack walls. Favorable crack arrays for initial growth in the model include particular sets of en-echelon cracks inclined 45° to the stress axis, but the critical orientation may be smaller in brittle rock if crack-wall friction in rock is greater than in the model. Axial growth of en-echelon cleavage cracks, inclined at angles smaller than 45°, was observed in feldspar during deformation of a pegmatite having a mineralogy and texture similar to granite. Their growth follows predictions derived from photoelastic model studies.Crack growth in the pegmatite begins between half and two-thirds of the ultimate strength. The first flaws to grow, however, also include pre-existing axially oriented cleavage cracks in the feldspar. Crack growth occurs randomly throughout the specimen as stress is increased, without much evidence that grain boundaries are activated for crack growth. But when the applied stress approaches the ultimate strength, two new features are observed. There is an abrupt development of finite frictional slip along favorably inclined flaws and grain boundaries, beginning with displacements of the order of the dimensions of grain-boundary asperities. Crack growth still occurs at various locations throughout the specimen at this stage, but there is also a detectable concentration of growth along potential shear zones. Flaw-wall friction appears to be one of the critical factors that determine the pegmatite's ultimate strength and the instability of through-going fracture. 相似文献
19.
Thierry Reuschié 《地学学报》1992,4(5):591-597
Fracture in a heterogeneous medium is simulated by a network model. Heterogeneity is introduced by randomly distributing an initial population of cracks on a triangular network of bonds. Each ruptured bond represents a crack. External stresses are applied onto the network and the evolution of the crack population is analysed when the stresses are increased. A crack will propagate to an intact bond if the critical crack extension force is reached. This propagation leads to the coalescence of bonds: crack clusters are formed. When their shapes are complex we approximate them by simpler ones to compute a crack extension force. Crack interactions are introduced by adding for a given bond the contributions of all adjacent dusters. By using this kind of approach combining fracture mechanics and network modelling, we are able to simulate the rupture of a rock specimen under compression/ traction or shear stresses and to give some characteristics of the macroscopic fracture. 相似文献
20.
The present numerical study, which is an extension of our previous numerical analysis on cracking processes of a single pre-existing flaw, focuses on the coalescence of two pre-existing parallel open flaws in rock subjected to a uniaxial compressive loading. To facilitate a systematic investigation, the arrangements of the flaw pair are classified into 11 categories. Simulations engaging AUTODYN are conducted on each category. The numerical results are compared with some published physical experimental test results. Eleven typical coalescence patterns are obtained, which are in good agreement with the experimental results, which include two coalescence patterns obtained in flaw pair arrangements (II) and (VIII″) not being reported in previous studies. The information gathered in the simulations helps identify the type (tensile/shear) of each crack segment involved in the coalescence. Most of the coalescence cracks initiate at or around the flaw tips, except those in flaw pair arrangements (II) and (IX′) with a very short ligament length, in which the coalescence cracks initiate on the flaw surfaces away from the flaw tip regions. Based on the numerical simulation results, the properties of the 11 coalescence patterns are obtained. Except those in flaw pair arrangements (II) and (IX′), the other coalescence patterns can be interpreted with respect to the basic crack types—tensile wing crack, horsetail crack and anti-wing crack. In addition, based on the type of crack segments involved in coalescence, namely tensile and shear, the coalescence can be classified into T mode (tensile mode), S mode (shear mode) and TS mode (mixed tensile–shear mode). 相似文献