共查询到20条相似文献，搜索用时 31 毫秒
1.
2.
Damage‐viscoplastic consistency model for rock fracture in heterogeneous rocks under dynamic loading
Timo Saksala 《国际地质力学数值与分析法杂志》2010,34(10):1041-1062
This paper presents a damage‐viscoplastic consistency model for numerical simulation of brittle fracture in heterogeneous rocks. The model is based on a combination of the recent viscoplastic consistency model by Wang and the isotropic damage concept with separate damage variables in tension and compression. This approach does not suffer from ill‐posedness, caused by strain softening, of the underlying boundary/initial value problem since viscoplasticity provides the regularization by introducing a length scale effect under dynamic loading conditions. The model uses the Mohr–Coulomb yield criterion with the Rankine criterion as a tensile cut‐off. The damage law in compression is calibrated via the degradation index concept of Fang and Harrison. Thereby, the model is able to capture the brittle‐to‐ductile transition occurring in confined compression at a certain level of confinement. The heterogeneity of rock is accounted for by the statistical approach based on the Weibull distribution. Numerical simulations of confined compression test in plane strain conditions demonstrate a good agreement with the experiments at both the material point and structural levels as the fracture modes are realistically predicted. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
3.
The solutions of stress and displacement of a circular opening excavated in brittle and strain-softening rock mass incorporating rockbolts effectiveness and seepage force are presented in this study. The evolution equation is reconstructed for the strength parameters that incorporate these factors. Based on the evolution equation, an improved numerical method and stepwise procedure are presented which are compatible with the Mohr–Coulomb (M–C) and the generalized Hoek–Brown (H–B) failure criteria, respectively. Then given three interaction mechanisms between rockbolts and surrounding rock, solutions for stress and displacement are proposed in line with the improved numerical method and numerical stepwise procedure. The proposed approach can be reduced to Fahimifar and Soroush’s (Tunn Undergr Space Technol 20:333–343, 2005) solutions for special cases. The proposed method was validated by field monitoring data and FLAC results of Yanzidong tunnel. Examples under the M–C and generalized H–B failure criteria for rock mass are generated through MATLAB programming. Moreover, parametric studies are conducted to highlight the influence of rockbolts effectiveness in combination with seepage force on the stress and displacement of very good, average, and very poor surrounding rock. Results show that in this case, stress confinement is higher and tunnel convergences are lower than the corresponding stresses and displacements obtained in non-reinforced tunnels. Displacement and plastic radius are also higher than those without considering seepage force. 相似文献
4.
The action of tunnel excavation reduces the in-situ stresses along the excavated circumference and can therefore be simulated by unloading of cavities from the in-situ stress state. Increasing evidence suggests that soil behavior in the plane perpendicular to the tunnel axis can be modelled reasonably by a contracting cylindrical cavity, while movements ahead of an advancing tunnel heading can be better predicted by spherical cavity contraction theory. In the past, solutions for unloading of cavities from in-situ stresses in cohesive-frictional soils have mainly concentrated on the small strain, cylindrical cavity model. Large strain spherical cavity contraction solutions with a non-associated Mohr–Coulomb model do not seem to be widely available for tunnel applications. Also, cavity unloading solutions in undrained clays have been developed only in terms of total stresses with a linear elastic-perfectly plastic soil model. The total stress analyses do not account for the effects of strain hardening/softening, variable soil stiffness, and soil stress history (OCR). The effect of these simplifying assumptions on the predicted soil behavior around tunnels is not known. In this paper, analytical and semi-analytical solutions are presented for unloading of both cylindrical and spherical cavities from in-situ state of stresses under both drained and undrained conditions. The non-associated Mohr-Coulomb model and various critical state theories are used respectively to describe the drained and undrained stress-strain behaviors of the soils. The analytical solutions presented in this paper are developed in terms of large strain formulations. These solutions can be used to serve two main purposes: (1) to provide models for predicting soil behavior around tunnels; (2) to provide valuable benchmark solutions for verifying various numerical methods involving both Mohr–Coulomb and critical state plasticity models. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
5.
Enriched finite element procedures for analyzing decoupled bolts installed in rock mass 总被引：1，自引：0，他引：1
Numerical procedures are developed to analyze interaction between fully grouted bolts and rock mass using ‘enriched finite element method (EFEM)’. A solid element intersected by a rock bolt along any arbitrary direction is termed as ‘enriched’ element. The nodes of an enriched element have additional degrees of freedom for determining displacements, stresses developed in the bolt rod. The stiffness of the enriched element is formulated based on properties of rock mass, bolt rod and grout, orientation of the bolt and borehole diameter. Decoupling at grout–bolt interface and elasto‐plastic behavior of rock mass have also been incorporated into the EFEM procedures. The results of this method are compared with analytical pull‐out test results presented by Li and Stillborg (Int. J. Rock Mech. Min. Sci. 1999; 36 :1013–1029). In addition, a numerical example of a bolted tunnel is provided to demonstrate the efficacy of the proposed method for practical applications. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
6.
A large strain analysis of undrained expansion of a spherical/cylindrical cavity in a soil modelled as non‐linear elastic modified Cam clay material is presented. The stress–strain response of the soil is assumed to obey non‐linear elasticity until yielding. A power‐law characteristic or a hyperbolic stress–strain curve is used to describe the gradual reduction of soil stiffness with shear strain. It is assumed that, after yielding, the elasto‐plastic behaviour of the soil can be described by the modified Cam clay model. Based on a closed‐form stress–strain response in undrained condition, a numerical solution is obtained with the aid of simple numerical integration technique. The results show that the stresses and the pore pressure in the soil around an expanded cavity are significantly affected by the non‐linear elasticity, especially if the soil is overconsolidated. The difference between large strain and small strain solutions in the elastic zone is not significant. The stresses and the pore pressure at the cavity wall can be expressed as an approximate closed‐form solution. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
7.
In this article we present closed‐form solutions for the undrained variations in stress, pore pressure, deformation and displacement inside hollow cylinders and hollow spheres subjected to uniform mechanical pressure instantaneously applied to their external and internal boundary surfaces. The material is assumed to be a saturated porous medium obeying a Mohr–Coulomb model failure criterion, exhibiting dilatant plastic deformation according to a non‐associated flow rule which accounts for isotropically strain hardening or softening. The instantaneous response of a porous medium submitted to an instantaneous loading is undrained, i.e. without any fluid mass exchange. The short‐term equilibrium problem to be solved is now formally identical to a problem of elastoplasticity where the constitutive equations involve the undrained elastic moduli and particular equivalent plastic parameters. The response of the model is presented (i) for extension and compression undrained triaxial tests, and (ii) for unloading problems of hollow cylinders and spheres through the use of appropriately developed closed‐form solutions. Numerical results are presented for a plastic clay stone with strain hardening and an argilite with strain softening. The effects of plastic dilation, of the strain softening law and also of geometry of the cavity on the behaviour of the porous medium have been underlined. Analytical solutions provide valuable benchmarks enabling various numerical methods in undrained conditions with a finite boundary to be verified. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
8.
基于Hoek-Brown准则的隧洞围岩变形研究 总被引：1，自引：0，他引：1
由于在许多实际条件下,比如节理岩体中,线性的M-C准则不太适用,非线性的Hoek-Brown比较适用,因此,可以尝试使用这一非线性屈服准则对洞室变形进行研究。研究隧洞变形时,将围岩分为弹性区、应变软化区、塑性流动区。采用Hoek-Brown准则和非关联流动法则对洞室变形进行了理论推导;软化区域围岩参数随着塑性变形增加而变化,解析法难以求得应力,采用龙格－库塔方法进行数值计算,求解得到塑性软化区和流动区半径,并最终求得洞室变形。通过算例计算表明,在不考虑软化区和流动区时,方法和Carranza-Torres计算结果相差甚小;随着原岩应力的增加,膨胀角对洞室变形的影响增大。 相似文献
9.
考虑塑性应变率梯度的单轴压缩岩样轴向响应 总被引：9，自引：7，他引：2
基于剪切应变率梯度格式,采用解析方式研究了岩石材料在单轴压缩条件下的应变软化的结构响应。根据非局部连续介质模型,提出了一维二阶剪切应变率梯度格式。非局部剪切应变率与局部剪切应变率及其二阶梯度有关。将经典塑性理论中的局部剪切应变率替换为非局部剪切应变率,可以直接得到局部剪切应变率的封闭解析解,而不必通过将局部剪切应变对时间求导获得。通过对局部剪切应变率积分,得到了沿剪切带方向的相对剪切速度。试件峰值强度后的端部速度由弹性及塑性两部分构成。前一部分由虎克定律描述;后一部分与相对剪切速度有关。对弹性及塑性两部分速度求和,得到了单轴压缩岩样剪切破坏问题轴向响应的解析式。研究表明：试样高度越大、内部长度越小、剪切软化模量越大及泊松比越小,则岩样的轴向响应倾向于脆性。根据岩样与矿柱的相似性,岩样响应倾向于脆性,意味着矿柱将失去稳定性,发生矿柱岩爆。目前的基于剪切应变率梯度格式的主要优点是简洁。 相似文献
10.
Reaction softening by dissolution–precipitation creep in a retrograde greenschist facies ductile shear zone,New Hampshire,USA 
下载免费PDF全文
下载免费PDF全文 R. J. McAleer D. L. Bish M. J. Kunk K. R. Sicard P. M. Valley G. J. Walsh B. A. Wathen R. P. Wintsch 《Journal of Metamorphic Geology》2017,35(1):95-119
We describe strain localization by a mixed process of reaction and microstructural softening in a lower greenschist facies ductile fault zone that transposes and replaces middle to upper amphibolite facies fabrics and mineral assemblages in the host schist of the Littleton Formation near Claremont, New Hampshire. Here, Na‐poor muscovite and chlorite progressively replace first staurolite, then garnet, and finally biotite porphyroblasts as the core of the fault zone is approached. Across the transect, higher grade fabric‐forming Na‐rich muscovite is also progressively replaced by fabric‐forming Na‐poor muscovite. The mineralogy of the new phyllonitic fault‐rock produced is dominated by Na‐poor muscovite and chlorite together with late albite porphyroblasts. The replacement of the amphibolite facies porphyroblasts by muscovite and chlorite is pseudomorphic in some samples and shows that the chemical metastability of the porphyroblasts is sufficient to drive replacement. In contrast, element mapping shows that fabric‐forming Na‐rich muscovite is selectively replaced at high‐strain microstructural sites, indicating that strain energy played an important role in activating the dissolution of the compositionally metastable muscovite. The replacement of strong, high‐grade porphyroblasts by weaker Na‐poor muscovite and chlorite constitutes reaction softening. The crystallization of parallel and contiguous mica in the retrograde foliation at the expense of the earlier and locally crenulated Na‐rich muscovite‐defined foliation destroys not only the metastable high‐grade mineralogy, but also its stronger geometry. This process constitutes both reaction and microstructural softening. The deformation mechanism here was thus one of dissolution–precipitation creep, activated at considerably lower stresses than might be predicted in quartzofeldspathic rocks at the same lower greenschist facies conditions. 相似文献
11.
Timo Saksala 《国际地质力学数值与分析法杂志》2010,34(13):1362-1386
This paper presents a damage–viscoplastic cap model for rocks with brittle and ductile behavior under low‐velocity impact loading, which occurs, e.g. in percussive drilling. The model is based on a combination of the recent viscoplastic consistency model by Wang and the isotropic damage concept. This approach does not suffer from ill posedness—caused by strain softening—of the underlying boundary/initial‐value problem since viscoplasticity provides a regularization under dynamic loading by introducing an internal length scale. The model uses the Drucker–Prager (DP) yield function with the modified Rankine criterion as a tension cut‐off and a parabolic cap surface as a compression cut‐off. The parabolic cap is smoothly fitted to the DP cone. The strain softening law in compression is calibrated with the degradation index concept of Fang and Harrison. Thereby, the model is able to capture the brittle‐to‐ductile transition and hardening behavior of geomaterials under highly confined compression, which is the prevailing stress state under a bit‐button in percussive drilling. Rock strength heterogeneity is characterized statistically at the structural level using the Weibull distribution. An explicit time integrator is chosen for solving the FE‐discretized equations of motion. The contact constraints due to the impact of an indenter are imposed with the forward increment Lagrange multiplier method that is compatible with explicit time integrators. The model is tested at the material point level with various uniaxial and triaxial tests. At the structural level confined compression, uniaxial tension tests and a rock sample under low‐velocity impact are simulated. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
12.
13.
The essential difference in the formation of conjugate shear zones in brittle and ductile deformation is that the intersection angle between brittle conjugate faults in the contractional quadrants is acute(usually ~60°) whereas the angle between conjugate ductile shear zones is obtuse(usually 110°). The Mohr-Coulomb failure criterion, an experimentally validated empirical relationship, is commonly applied for interpreting the stress directions based on the orientation of the brittle shear fractures. However, the Mohr-Coulomb failure criterion fails to explain the formation of the low-angle normal fault, high-angle reverse fault, and the conjugate strike-slip fault with an obtuse angle in the σ1 direction. Although it is ten years since the Maximum-Effective-Moment(MEM) criterion was first proposed, and increasingly solid evidence in support of it has been obtained from both observed examples in nature and laboratory experiments, it is not yet a commonly accepted model to use to interpret these antiMohr-Coulomb features that are widely observed in the natural world. The deformational behavior of rock depends on its intrinsic mechanical properties and external factors such as applied stresses, strain rates, and temperature conditions related to crustal depths. The occurrence of conjugate shear features with obtuse angles of ~110° in the contractional direction on different scales and at different crustal levels are consistent with the prediction of the MEM criterion, therefore ~110° is a reliable indicator for deformation localization that occurred at medium-low strain rates at any crustal levels. Since the strain–rate is variable through time in nature, brittle, ductile, and plastic features may appear within the same rock. 相似文献
14.
E. Alonso L. R. Alejano F. Varas G. Fdez‐Manin C. Carranza‐Torres 《国际地质力学数值与分析法杂志》2003,27(13):1153-1185
A literature review has shown that there exist adequate techniques to obtain ground reaction curves for tunnels excavated in elastic‐brittle and perfectly plastic materials. However, for strain‐softening materials it seems that the problem has not been sufficiently analysed. In this paper, a one‐dimensional numerical solution to obtain the ground reaction curve (GRC) for circular tunnels excavated in strain‐softening materials is presented. The problem is formulated in a very general form and leads to a system of ordinary differential equations. By adequately defining a fictitious ‘time’ variable and re‐scaling some variables the problem is converted into an initial value one, which can be solved numerically by a Runge–Kutta–Fehlberg method, which is implemented in MATLAB environment. The method has been developed for various common particular behaviour models including Tresca, Mohr–Coulomb and Hoek–Brown failure criteria, in all cases with non‐associative flow rules and two‐segment piecewise linear functions related to a principal strain‐dependent plastic parameter to model the transition between peak and residual failure criteria. Some particular examples for the different failure criteria have been run, which agree well with closed‐form solutions—if existing—or with FDM‐based code results. Parametric studies and specific charts are created to highlight the influence of different parameters. The proposed methodology intends to be a wider and general numerical basis where standard and newly featured behaviour modes focusing on obtaining GRC for tunnels excavated in strain‐softening materials can be implemented. This way of solving such problems has proved to be more efficient and less time consuming than using FEM‐ or FDM‐based numerical 2D codes. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
15.
基于Drucker-Prager（下简称D-P）准则,建立压缩载荷作用下的非贯通节理岩石的弹塑性断裂模型。针对节理岩石小范围屈服翼裂纹尖端塑性区,推导了D-P屈服准则的纯I、纯II及I、II复合型3种翼裂纹无量纲塑性区径长函数,并与Mises准则的塑性区进行对比;结果表明,D-P准则的I型和复合型塑性区较Mises屈服准则的塑性区大,且其II型及I、II复合型塑性区在翼裂纹上下表面不连续。进一步,引入断裂软化因子以表征节理岩石裂隙断裂扩展后的断裂软化规律,考虑非贯通节理岩石复合型断裂软化,是由于节理尖端翼裂纹应变能密度超过最小应变能密度导致其成核扩展引起的,提出用应变能密度的指数函数形式表征断裂软化变量的演化;塑性屈服函数采用Borja等的应力张量3个不变量的硬化/软化函数,反映塑性内变量及应力状态对硬化函数的影响;建立节理岩石的弹塑性断裂本构关系及其数值算法,并用回映隐式积分算法编制了弹塑性断裂模型的程序。以单轴压缩下非贯通节理岩石为例,分析岩石断裂韧度、节理摩擦系数和节理倾角等参数的影响,结果表明,所提出的弹塑性断裂模型与数值和试验结果比较吻合。 相似文献
16.
随着地下工程开挖深度的增加,深部岩体将处于高应力和复杂的地质环境中,产生与浅埋洞室破坏模式迥异的分区破裂现象。深部洞室在动力卸荷作用下,基于应变梯度理论和损伤软化模型,建立了弹塑性损伤软化动力模型,推导了含有应变梯度项的运动方程、平衡方程和边界条件,提出相应的破坏判据,采用Runge-Kutta方法和Matlab数值软件求得不同卸载时刻围岩附加位移场、应力场和开挖后围岩总位移场、应力场,得到深部洞室围岩分区破裂的动态形成过程和发展规律。由理论计算值与地质力学模型试验实测值对比分析得知,围岩的径向位移、径向应力和切向应力出现波峰和波谷交替振荡的变化规律,理论计算得到的破裂区和非破裂区的宽度和数量与试验实测值有很好的一致性,证实了该模型分析分区破裂现象的适用性,对以后深部地下工程围岩变形破坏和支护设计提供理论支持。 相似文献
17.
We consider discontinuous bifurcations as the indicator of a localized failure for a class of composites that are characterized by elastic fibres reinforcing an elastic–plastic matrix. A macroscopic tangent stiffness tensor for the fibre‐reinforced composite is developed by consistently homogenizing the contribution of fibres in a spherical representative volume element. Analytical solutions are derived for the critical hardening modulus and corresponding bifurcation directions for the case of plane strain loading. Properties of the solutions are further illustrated on the example of the non‐associated Drucker–Prager model at onset of yielding. Results show that presence of fibres decreases the critical hardening modulus, thus inhibiting the onset of strain localization. The rate of decrease in the critical hardening modulus is the highest for pure shear, followed by uniaxial tension, uniaxial compression, biaxial tension and biaxial compression. The main fibre parameters that control the onset of strain localization are their volumetric content and their stiffness modulus whereby very stiff fibres can produce the most significant decrease in the critical hardening modulus, especially for the state of biaxial tension. The critical hardening modulus for the non‐associated Drucker–Prager model exhibits a full range of localization modes including compaction bands, dilation bands, and transition in the form of shear bands regardless of the presence of fibres. Presence of fibres affects bifurcation directions, except in the case when Poisson's ratio of the matrix is equal to 0.25. The results demonstrate stabilizing effects of fibres by which they provide the control against the onset of strain localization. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
18.
Large strain similarity solution for a spherical or circular opening excavated in elastic‐perfectly plastic media 下载免费PDF全文
下载免费PDF全文 K.H. Park 《国际地质力学数值与分析法杂志》2015,39(7):724-737
This paper deals with the unloading problem of a spherical or circular opening excavated in elastic‐perfectly plastic media with a nonassociated Mohr–Coulomb yield criterion. A large strain similarity solution, using incremental velocity approach, is presented by replacing partial differential equations from stress equilibrium, constitutive law, consistency condition, and displacement equation with first‐order ordinary differential equations. The classical Runge–Kutta method is used to solve the first‐order ordinary differential equations. Comparisons among small and large strain solutions are made using some data sets of soil and rock. The results show that the displacements by large strain similarity solution are smaller than those by exact small strain solution and somewhat larger than those by large strain solution using total strain approach. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
19.
A plastic deviatoric model with hardening is developed on the basis of geomechanical tests performed in the saturated case on low permeable porous material such as argillite. This model is a generalized Mohr–Coulomb plastic criterion combined with a Drucker–Prager plastic potential and the hardening parameter is the plastic distortion. Three different hardening functions have been introduced on the basis of triaxial tests: an increase of friction angle, a decrease of cohesion after a threshold and a contractancy to dilatancy transition for volumetric plastic strain. This plastic model has been adapted to the partially saturated case. The effective stress is expressed thanks to the equivalent interstitial pressure π. Numerical results are presented for the excavation and monotonous ventilation of a deep cylindrical cavity. A first plastification due to excavation is followed by a second one due to desaturation. The extent of the non-saturated zone provokes an extent of a plastic zone in the rock mass. Analysis shows that the origin of the plastification can be found in the deviatoric stresses because mean effective stresses are compressive during drying. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
20.
Torsional tests are frequently used in material testing to determine shear stress–strain relationships. In the tests, torque is applied to a solid or hollow cylindrical specimen, whose twist angle is recorded as the response. While the raw test results are readily available in terms of torque–twist sets, converting them into intrinsic stress–strain relationships is not always straightforward, owing to the nonlinear stress–strain behavior and the non‐uniform radial distribution of the stress and strain fields. It follows that the approaches currently used to interpret torsional data are essentially empirical based. The present paper addresses this data interpretation issue from a more rigorous perspective. The paper includes two parts: the first part concerns hollow cylindrical specimens under torsion, where an analytical‐based recursive procedure is developed; and the second part deals with solid torsional specimens, for which an analytical solution given by Nadai is directly employed. Examples pertaining to these two cases are examined and the resulting shear stress–strain curves are compared with those based on a number of popular empirical methods. The proposed analytical approaches can serve as a benchmark for checking the accuracy and validity of those empirical propositions. It is worth noting that, for brittle materials exhibiting softening behavior, only the analytical method gives the correct interpretation. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献