共查询到20条相似文献，搜索用时 27 毫秒
1.
In the absence of initial cracks, the material behavior is limited by its strength, usually defined in homogeneous conditions (of stress and strain). Beyond this limit, in quasi‐brittle case, cracks may propagate and the material behavior tends to be well described by fracture mechanics. Discrete element approaches show consistent results dealing with this transition during rupture. However, the calibration of the parameters of the numerical models (i.e., stiffness, strength, and toughness) may be quite complex and sometimes only approximative. Based on a brittle rupture criterion, we analyze the biaxial response of uncracked samples. Thus, tensile and compressive strengths are analytically identified and become direct parameters of our discrete model. Furthermore, a physically reliable crack initiation (and subsequent propagation) is shown to be induced during rupture and verified by the simulation of three‐point bending and diametral compression tests. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
2.
We investigate the initial rupture process of microearthquakes to reveal relationships between nucleation region sizes and eventual earthquake sizes. In order to obtain high quality waveform data, we installed a trigger recording system with a sampling frequency of 10 kHz at the base of a deep borehole at the Nojima Fault, Japan. We analyze waveform data of 31 events around the borehole, with seismic moment ranging from 4.2 × 109 Nm to 7.1 × 1011 Nm. We use both a circular crack model with an accelerating rupture velocity (SK model) [Sato, T., Kanamori, H., 1999. Beginning of earthquakes modeled with the Griffith's fracture criterion, Bull. Seism. Soc. Am., 89, 80-93.], which generates a slow initial phase of velocity pulse, and a circular crack model with a constant rupture velocity (SH model) [Sato, T, Hirasawa, T., 1973. Body wave spectra from propagating shear cracks, J. Phys. Earth, 21, 415-431.], which generates a ramp-like velocity pulse. Source parameters of these two models are estimated by waveform inversion of the first half cycle of the observed velocity pulse applying both a grid search and a non-linear least squares method. 14 of 31 events are never reproduced by the SH model with a constant Q operator. But SK model with a constant Q operator provides a size of the pre-existing crack, corresponding to the size of the nucleation regions, and a size of the eventual crack. We recognize that (i) the eventual seismic moment is approximately scaled as the cube of the size of pre-existing cracks, (ii) the eventual seismic moment is scaled as the cube of the size of eventual cracks, and (iii) the size of eventual cracks is roughly proportional to the size of pre-existing cracks. We, thus, conclude that the size of eventual earthquakes is controlled by the size of the nucleation regions. 相似文献
3.
An energy approach is proposed as a complement to the stress approach commonly considered for investigating soil desiccation cracking. The elastic strain energies before and after crack initiation are estimated by both numerical and analytical solutions. The energy released by cracking is then compared with the fracture energy to discuss crack initiation conditions. This leads to combined energy and stress conditions for crack initiation following Leguillon's theory. An approximate analytical solution is derived from a variational formulation of the porous elastic body equations. A cohesive zone model and finite element code are used to simulate crack propagation in an unsaturated porous body. This analysis shows that the energy criterion is reached before the stress criterion, and this can explain unstable crack propagation at the beginning. The approximate analytical solution allows predicting correctly the crack depth and opening in its initiation stage. 相似文献
4.
A relationship between tensile strength and loading stress governing the onset of mode I crack propagation obtained via numerical investigations using a bonded particle model 
下载免费PDF全文
下载免费PDF全文 Numerical models based on the discrete element method are used to study the fracturing process in brittle rock‐like materials under direct and indirect tension. The results demonstrate the capacity of the model to capture the essential characteristics of fracture including the onset of crack propagation, stable and unstable crack growth, arrest and reinitiation of fracturing, and crack branching. Simulations of Brazilian indirect tension tests serve to calibrate the numerical model, relating macroscopic tensile strength of specimens to their micromechanical breakage parameters. A second suite of simulations reveals a linear relationship between the tensile strength of specimens and the loading stress for which mode I tensile crack propagation ensues. Based on these results, a crack initiation criterion for brittle materials is proposed, prescribing the stressing conditions required to induce tensile failure. Such a criterion, if broadly applicable, provides a practical means to rapidly assess the failure potential of brittle materials under tensile loads. 相似文献
5.
This article is devoted to numerical modeling of anisotropic damage and plasticity in saturated quasi‐brittle materials such as rocks and concrete. The damaged materials are represented by an isotropic poroelastic matrix containing a number of families of microcracks. Based on previous works, a discrete thermodynamic approach is proposed. Each family of microcracks exhibits frictional sliding along crack surfaces as well as crack propagation. The frictional sliding is described by a Coulomb–Mohr‐type plastic criterion by taking into account the effect of fluid pressure through a generalized effective stress concept. The damage evolution is entirely controlled by and coupled with the frictional sliding. The effective elastic properties as well as Biot's coefficients of cracked porous materials are determined as functions of induced damage. The inelastic deformation due to frictional sliding is also taken into account. The procedure for the identification of the model's parameters is presented. The proposed model is finally applied to study both mechanical and poromechanical responses of a typical porous brittle rock in drained and undrained compression tests as well as in interstitial pressure controlled tests. The main features of material behaviors are well reproduced by the model. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
6.
基于Drucker-Prager（下简称D-P）准则,建立压缩载荷作用下的非贯通节理岩石的弹塑性断裂模型。针对节理岩石小范围屈服翼裂纹尖端塑性区,推导了D-P屈服准则的纯I、纯II及I、II复合型3种翼裂纹无量纲塑性区径长函数,并与Mises准则的塑性区进行对比;结果表明,D-P准则的I型和复合型塑性区较Mises屈服准则的塑性区大,且其II型及I、II复合型塑性区在翼裂纹上下表面不连续。进一步,引入断裂软化因子以表征节理岩石裂隙断裂扩展后的断裂软化规律,考虑非贯通节理岩石复合型断裂软化,是由于节理尖端翼裂纹应变能密度超过最小应变能密度导致其成核扩展引起的,提出用应变能密度的指数函数形式表征断裂软化变量的演化;塑性屈服函数采用Borja等的应力张量3个不变量的硬化/软化函数,反映塑性内变量及应力状态对硬化函数的影响;建立节理岩石的弹塑性断裂本构关系及其数值算法,并用回映隐式积分算法编制了弹塑性断裂模型的程序。以单轴压缩下非贯通节理岩石为例,分析岩石断裂韧度、节理摩擦系数和节理倾角等参数的影响,结果表明,所提出的弹塑性断裂模型与数值和试验结果比较吻合。 相似文献
7.
基于完全流-固耦合的弹塑性理论给出了水力压裂数值计算的弥散裂缝模型,其中材料的弹性部分采用线弹性本构关系,塑性部分采用摩尔-库仑破坏准则及强化准则。依据当前的有效应力状态修正渗透系数来模拟压裂液在裂缝中的流动。渗透系数的修改使用双曲正切函数,并采用平均有效应力作为水力裂缝的起裂判据。在ABAQUS软件中通过用户自定义程序添加了该模型。根据岩石的切面照片建立了含有硬包裹体分布的非均质岩石的有限元计算模型,模拟了中心点注水条件下的水力压裂传播过程,讨论了在常应力状态下非均质岩石中开裂区域、典型位置的应力路径变化和裂缝传播范围随时间变化的特点。进行了多种条件下含有硬包裹体分布的岩石材料的数值试验,得出了基岩材料的弹性模量、凝聚力和渗透系数以及注水速率对峰值注水压力、平均注水压力和裂缝开度的影响规律。 相似文献
8.
Inglis [1] has solved the problem of distribution of stress in an elastic plate around an elliptical hole. His works clarify the role of cracks in the failure of an elastic material. However, his solution cannot be applied to saturated clay because he considers only total stresses, while, in saturated clay, the criterion of rupture should be expressed in terms of effective and not total stresses. The solution of Atkinson and Craster [2] using Biot's poroelasticity theory, shows that there is no high pore pressure in the vicinity of the crack tips for saturated clay. The major difference between this approach and the Biot's theory of is that, in saturated clay, strain is a function of the variation of the effective stress [3], while, in poroelastic media, strain is only a function of the variation of the total stress [4, Equation 2.2]. Also in their solution there is continuity between the pore fluid and the inner fluid in the crack. Their solution is valid for poroelastic media involving a movement of the pore fluid. In our solution there is no movement of the pore fluid (Undrained condition). In this paper we have solved the same problem as Inglis [1], but for the particular case of saturated clay obeying elastic law. By solving this problem we obtained the expressions for pore pressure, effective stress, total stress and displacements. The results show that not only the total stress but also the pore pressure and the effective stress are also high in the vicinity of the crack tips. A new failure criterion, based on Griffith's strain energy principle [5] and maximum tensile stress [6], valid for saturated clay is developed in this paper. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
9.
硬岩裂纹起裂强度和损伤强度取值方法探讨 总被引：1，自引：0，他引：1
岩石裂纹起裂强度σci和损伤强度σcd的研究对探讨岩石渐进破坏过程和脆性破坏机制,建立岩石的起裂强度和长期强度准则具有重要意义。鉴于此,系统总结了岩石起裂强度和损伤强度的4种计算方法,包括裂纹应变模型计算法,声发射参数取值法,侧向应变、体积应变曲线观察法和移动点回归法,并以大理岩和花岗岩在单轴压缩条件下为例,详细描述了用不同方法计算σci和σcd时存在的问题。研究结果表明,裂纹应变模型法物理意义明确,可以方便和准确地确定岩石的σci和σcd;声发射撞击率可以定性或半定量了解岩石的渐进破坏和损伤过程,但无法准确定量确定σci和σcd;移动点回归法计算体积应变刚度时受取值点间隔影响较大,而且加载过程中不同阶段的特征并不明显,难以确定起裂强度值。研究结论可为加深对岩石裂纹起裂强度和损伤强度的重要意义的认识,了解两个强度值的计算方法并合理选取其计算方法提供有益的参考。 相似文献
10.
将闭合裂纹表面的有效剪应力引入裂尖应力强度因子的计算,获得了无限大板和有限宽板含中心闭合裂纹在不同裂纹长度、倾角以及摩擦系数下裂纹尖端的应力强度因子值。引入等径向剪应力线 这一概念,建立了闭合裂纹断裂的有效剪应力准则：（1）岩板内闭合裂纹将沿着等 线上双剪应力的和最小的方向扩展;（2）裂纹尖端的应力强度因子KⅡ达到材料的临界值KⅡC,裂纹将开始扩展。该准则成功预测了闭合裂纹的临界起裂角 ,与各种经典复合型断裂准则计算Ⅱ型裂纹起裂角比较,结果较为接近。将其应用于闭合裂纹的断裂判定是安全的。 相似文献
11.
在研究双向压缩条件下压剪复合型裂纹应力分布特征及断裂破坏机制基础上,考虑渗透压对初始裂隙面上有效正应力的影响,提出高低渗透压环境的判定准则,并基于滑动裂纹模型理论及最大周向拉应力破坏准则,得到不同渗透压环境下初始裂隙尖端微裂纹起裂特征与规律。研究结果表明：压剪复合应力条件下,初始裂隙尖端发育微裂纹的最优倾角与裂隙面摩擦系数直接相关,随裂隙面摩擦系数的增大,最优初始裂隙倾角由45°起逐渐增大;低渗透压条件下,渗流场的存在使裂纹面摩擦系数发生弱化,进而使得最优初始裂隙倾角向45°靠近,而渗透压直接降低裂隙面上有效正应力且与裂隙倾角无关,其仅仅影响裂隙体材料的初裂强度;高渗透压条件下,初始裂隙面由压剪复合应力状态转化为拉剪复合应力状态,并在拉剪复合应力场作用下,尖端微裂纹起裂角随KI/KII的不断增大,由70.5°逐渐趋近于0°。 相似文献
12.
考虑岩石闭合裂纹壁面间存在的摩擦力对裂纹尖端应力场的影响,应用最大周向应力理论得到压剪复合裂纹的断裂角。在此基础上,依据岩石裂纹尖端双向受力时的破坏特征,结合最大周向应力准则与修正的格里菲斯（Griffith）强度理论,建立了考虑摩擦效应的闭合裂纹失稳扩展的岩石压剪断裂判据。研究结果表明：断裂角受裂纹和荷载方向的夹角、裂纹壁面之间的摩擦系数、侧压力系数的影响;当压剪裂纹的断裂角是某个定值时,纯II型裂纹的断裂韧度与纯I型裂纹的断裂韧度的比值只与岩石裂纹表面的摩擦系数取值有关,而与其他岩石力学参数无关。此研究成果可为压剪应力作用下裂隙岩体的失稳破坏提供参考。 相似文献
13.
Computational model coupling mode II discrete fracture propagation with continuum damage zone evolution 下载免费PDF全文
下载免费PDF全文 We propose a numerical method that couples a cohesive zone model (CZM) and a finite element‐based continuum damage mechanics (CDM) model. The CZM represents a mode II macro‐fracture, and CDM finite elements (FE) represent the damage zone of the CZM. The coupled CZM/CDM model can capture the flow of energy that takes place between the bulk material that forms the matrix and the macroscopic fracture surfaces. The CDM model, which does not account for micro‐crack interaction, is calibrated against triaxial compression tests performed on Bakken shale, so as to reproduce the stress/strain curve before the failure peak. Based on a comparison with Kachanov's micro‐mechanical model, we confirm that the critical micro‐crack density value equal to 0.3 reflects the point at which crack interaction cannot be neglected. The CZM is assigned a pure mode II cohesive law that accounts for the dependence of the shear strength and energy release rate on confining pressure. The cohesive shear strength of the CZM is calibrated by calculating the shear stress necessary to reach a CDM damage of 0.3 during a direct shear test. We find that the shear cohesive strength of the CZM depends linearly on the confining pressure. Triaxial compression tests are simulated, in which the shale sample is modeled as an FE CDM continuum that contains a predefined thin cohesive zone representing the idealized shear fracture plane. The shear energy release rate of the CZM is fitted in order to match to the post‐peak stress/strain curves obtained during experimental tests performed on Bakken shale. We find that the energy release rate depends linearly on the shear cohesive strength. We then use the calibrated shale rheology to simulate the propagation of a meter‐scale mode II fracture. Under low confining pressure, the macroscopic crack (CZM) and its damaged zone (CDM) propagate simultaneously (i.e., during the same loading increments). Under high confining pressure, the fracture propagates in slip‐friction, that is, the debonding of the cohesive zone alternates with the propagation of continuum damage. The computational method is applicable to a range of geological injection problems including hydraulic fracturing and fluid storage and should be further enhanced by the addition of mode I and mixed mode (I+II+III) propagation. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
14.
Analytical simulation of the dynamic compressive strength of a granite using the sliding crack model
A sliding crack model is employed to simulate rock strength under dynamic compression. It is assumed that the growth and nucleation of a sliding crack array presented results in the shear fault failure and dominate the mechanical properties of rock material. The pseudo‐tractions method is used to calculate the stress intensity factor of the sliding crack array under compression. With the utilization of a dynamic crack growth criterion, the growth of the sliding crack array is studied and the simulated strengths of a granite under dynamic compression are correspondingly obtained. It is concluded that the simulated rock strengths increase with increasing strain rates at different confining pressures, and the rising rates have a trend to decrease with increasing confining pressures. It is also indicated that the simulated rock strengths increase with increment of confining pressure at different strain rates, and the rising rates are almost identical at different strain rates. The simulation results are validated by the experimental data for the granite. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
15.
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. 相似文献
16.
We presented a finite‐element‐based algorithm to simulate plane‐strain, straight hydraulic fractures in an impermeable elastic medium. The algorithm accounts for the nonlinear coupling between the fluid pressure and the crack opening and separately tracks the evolution of the crack tip and the fluid front. It therefore allows the existence of a fluid lag. The fluid front is advanced explicitly in time, but an implicit strategy is needed for the crack tip to guarantee the satisfaction of Griffith's criterion at each time step. We enforced the coupling between the fluid and the rock by simultaneously solving for the pressure field in the fluid and the crack opening at each time step. We provided verification of our algorithm by performing sample simulations and comparing them with two known similarity solutions. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
17.
在冻融条件下岩石微裂隙中的水发生相变,体积膨胀,对微裂隙产生很大的冻胀力,当冻胀力超过岩石的抗拉强度时,微裂隙扩展。温度升高时,水又进入新的微裂隙,如此反复循环造成了岩石的损伤。据此,将岩石中的微裂隙等效为扁平状椭圆裂隙,基于断裂力学建立了单条微裂隙下裂隙扩展长度与冻胀力的关系,考虑岩石中微裂隙的分布,将岩石冻融条件下的应变分解为初始损伤应变、附加损伤应变和塑性应变,建立了弹塑性冻融损伤本构模型。最后,通过岩石冻融试验对该模型的合理性进行了验证,结果表明,该模型能够较好地模拟岩石在不同冻融次数下的应力-应变关系曲线。 相似文献
18.
In a cracked material, the stress intensity factors (SIFs) at the crack tips, which govern the crack propagation and are associated with the strength of the material, are strongly affected by the crack inclination angle and the orientation with respect to the principal direction of anisotropy. In this paper, a formulation of the boundary element method (BEM), based on the relative displacements of the crack tip, is used to determine the mixed‐mode SIFs of isotropic and anisotropic rocks. Numerical examples of the application of the formulation for different crack inclination angles, crack lengths, and degree of material anisotropy are presented. Furthermore, the BEM formulation combined with the maximum circumferential stress criterion is adopted to predict the crack initiation angles and simulate the crack propagation paths. The propagation path in cracked straight through Brazilian disc specimen is numerically predicted and the results of numerical and experimental data compared with the actual laboratory observations. Good agreement is found between the two approaches. The proposed BEM formulation is therefore suitable to simulate the process of crack propagation. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
19.
In order to apply the mechanical properties (measured on material specimens or laboratory‐sized models) to large structures (such as concrete dams), a non‐linear theory able to predict the size‐scale effect has to be used. One of these theories was first proposed by Hillerborg and co‐workers (fictitious crack model) and is based on the earlier works by Barenblatt and Dugdale for metals (cohesive crack model). It is based on the existence of a fracture process zone (FPZ), where the material undergoes strain softening. The behaviour of the material outside the FPZ is linear elastic. A large number of short‐time laboratory tests were executed, by varying the load, under crack mouth opening displacement control. Since concrete exhibits a time‐dependent behaviour, an interaction between creep and micro‐crack growth occurs in the FPZ. Therefore, different testing conditions can be applied: rupture can be achieved by keeping the load constant before peak value (pre‐peak tests), or after peak value and after an unloading and reloading procedure (post‐peak tests). The crack propagation rate is shown to be small enough to neglect inertial forces and large enough to keep the time‐dependent behaviour of the process zone as dominant compared to the behaviour of the undamaged and viscoelastic zone. Due to the variability in material microstructure from one specimen to another, experimental data show large ranges of scatter. Well established methods in probability theory require sufficient experimental data in order to assume a probability density distribution. The objective of this study is to investigate the ranges of variation of the time response under constant load in simple structural elements associated with pre‐selected variation (fuzziness) in the main material parameters. For situations where the values of the material parameters are of a non‐stochastic nature, the fuzzy set approach to modelling variability has been proposed as a better and more natural approach. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
20.
Influence of liquid bridges on the mechanical behaviour of polydisperse granular materials 总被引：1，自引：0，他引：1
We investigate a polydisperse granular material in which the particle interactions are governed by a capillary force law. The cohesion force for a grain‐pair with unequal diameters is expressed as an explicit function of the inter‐particle distance and the volume of the liquid bridge. This analytical relation is validated by experiments on a reference material. Then, it is completed by a rupture criterion and cast in the form of a force law that accounts for solid contact, capillary force and rupture characteristics of a grain‐pair. Finally, in order to evaluate the influence of capillary cohesion on the macroscopic behaviour, radial and axial compression tests on cylindrical assemblies of wet particles are simulated using a 3D distinct element method. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献