Integration of a continuum damage model for shale with the cutting plane algorithm

The stability of integration is essential to numerical simulations especially when solving nonlinear problems. In this work, a continuum damage mechanics model proposed by the first author is implemented with an integration method named cutting plane algorithm (CPA) to improve the robustness of the simulation. This integration method is one type of return mapping algorithm that bypasses the need for computing the gradients. We compare the current integration method with the previous direct method, and the result shows that the cutting plane algorithm exhibits excellent performance under large loading rate conditions. To enhance accuracy of the new method, a control procedure is utilized in the implementation of the algorithm based on error analysis. Thereafter, the theory of poromechanics is utilized with the damage model to account for the effects of fluid diffusion. Laboratory tests simulated with finite element method illustrate distinct behaviors of shale with different loading rates and indicate the development of microcrack propagation under triaxial compression. Copyright © 2016 John Wiley & Sons, Ltd.     

基于连续介质模型的隧道衬砌结构可靠性研究

基于连续介质力学模型,提出了一种简单易行的可靠性分析方法,该方法能利用现行通用的有限元程序,能得到显式的功能函数,但又不同于响应面法。以公路隧道圆形隧洞的开挖支护为例,采用此方法将围岩参数c (黏聚力)、E,? 和二次衬砌厚度D作为随机变量,采用有限元程序得到单一变量变化下的衬砌轴应力值,进而拟合出衬砌轴应力关于单一随机变量的关系式,并通过多元线性或非线性回归得到衬砌轴应力关于所有随机变量的显示表达式,再对圆形隧洞的衬砌结构进行了敏感性和可靠性分析,得到了一些有意义的结论。     

Computational microstructure modeling to estimate progressive moisture damage behavior of asphaltic paving mixtures

This paper presents a computational microstructure model to estimate the progressive moisture damage of multiphase asphaltic paving mixtures. Moisture damage because of water transport is incorporated with mechanical loading through a finite element method. To simulate nonlinear damage evolution in the mixtures, the model includes Fickian moisture diffusion, a cohesive zone model to simulate the gradual fracture process, and a degradation characteristic function to represent the reduction of material properties because of moisture infiltration. With the model developed, various parametric analyses are conducted to investigate how each model parameter affects the material‐specific moisture damage mechanism and damage resistance potential of the mixtures. Analysis results clearly demonstrate the significance of physical and mechanical properties of mixture components and geometric characteristics of microstructure for the better design of asphaltic paving mixtures and roadway structures. Copyright © 2012 John Wiley & Sons, Ltd.     

Finite element simulations of 3D planar hydraulic fracture propagation using a coupled hydro-mechanical interface element

Two-dimensional hydraulic fracturing simulations using the cohesive zone model (CZM) can be readily found in the literature; however, to our knowledge, verified 3D cohesive zone modeling is not available. We present the development of a 3D fully coupled hydro-mechanical finite element method (FEM) model (with parallel computation framework) and its application to hydraulic fracturing. A special zero-thickness interface element based on the CZM is developed for modeling fracture propagation and fluid flow. A local traction-separation law with strain softening is used to capture tensile cracking. The model is verified by considering penny-shaped hydraulic fracture and plain strain Kristianovich‑Geertsma‑de Klerk hydraulic fracture (in 3D) in the viscosity- and toughness-dominated regimes. Good agreement between numerical results and analytical solutions has been achieved. The model is used to investigate the influence of rock and fluid properties on hydraulic fracturing. Lower stiffness tip cohesive elements tend to yield a larger elastic deformation around the fracture tips before the tensile strength is reached, generating a larger fracture length and lower fracture pressure compared with higher stiffness elements. It is found that the energy release rate has almost no influence on hydraulic fracturing in the viscosity-dominated regime because the energy spent in creating new fractures is too small when compared with the total input energy. For the toughness-dominated regime, the released energy during fracturing should be accurately captured; relatively large tensile strength should be used in order to match numerical results to the asymptotic analytical solutions. It requires smaller elements when compared with those used in the viscosity-dominated regime.     

Simulating damage evolution and fracture propagation in sandstone containing a preexisting 3-D surface flaw under uniaxial compression

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.     

Modeling hydraulic fracture propagation using cohesive zone model equipped with frictional contact capability

We present a new pore pressure cohesive element for modeling the propagation of hydraulically induced fracture. The Park-Paulino-Roesler cohesive zone model has been employed to characterize the fracturing behavior. Coulomb’s frictional contact model has been incorporated into the element to model the possible shear reactivation of pre-existing natural fractures. The developed element has been validated through a series of single-element tests and an available analytical solution. Furthermore, intersection behaviors between the hydraulic fracture and the natural fracture under various conditions have been predicted using the present element, which shows good agreement with experimental results.     

A continuum model of layered rock masses with non-associative joint plasticity

Layered rock masses can be modelled either as standard, orthotropic continua if the layer bending can be neglected or as Cosserat continua if the influence of layer bending is essential. This paper presents a finite element smeared joint model based on the Cosserat theory. The layers are assumed to be elastic with equal thickness and equal mechanical properties. All the cosserat parameters are expressed through the elastic properties of layers, layer thickness and joint stiffness. Plastic-slip as well as tensile-opening of layer interface (joint) are accounted for in a manner similar to the conventional non-associative plasticity theory. As an application, the behaviour of an excavation in a layered rock mass is examined. The displacement and stress fields given by smeared joint models based on the Cosserat continuum and the conventional anisotropic continuum approaches are compared with those obtained from the discrete joint model. The conventional anisotropic continuum model is found to break-down completely when the effective shear modulus in the direction parallel to layering is low in comparison to the shear modulus of the intact layer, whereas the Cosserat model is found to be capable of accurately reproducing complex load–deflection patterns irrespective of the differences in shear moduli. © 1998 John Wiley & Sons, Ltd.     

Numerical simulation of fracture flow with a mixed-hybrid FEM stochastic discrete fracture network model

We introduce a discrete fracture network model of stationary Darcy flow in fractured rocks. We approximate the fractures by a network of planar circle disks, which is generated on the basis of statistical data obtained from field measurements. We then discretize this network into a mesh consisting of triangular elements placed in three-dimensional space. We use geometrical approximations in fracture planes, which allow for a significant simplification of the final triangular meshes. We consider two-dimensional Darcy flow in each fracture. In order to accurately simulate the channeling effect, we assign to each triangle an aperture defining its hydraulic permeability. For the discretization we use the lowest order Raviart-Thomas mixed finite element method. This method gives quite an accurate velocity field, which is computed directly and which satisfies the mass balance on each triangular element. We demonstrate the use of this method on a model problem with a known analytical solution and describe the generation and triangulation of the fracture network and the computation of fracture flow for a particular real situation.     

基于连续损伤的水平井射孔-近井筒三维破裂模拟

射孔作为井筒与储层之间的液流通道,是水力压裂过程中的重要可控性参数。为研究水平井射孔-近井筒破裂机制,采用岩层变形-流体渗流方程描述应力状态变化,应用连续损伤破裂单元表征三维破裂位置与形态演化,并开发有限元求解程序模拟分析了射孔对水平井初始破裂压力、破裂位置及近井筒裂缝复杂性的调控作用。通过与解析模型及射孔压裂物理模型试验结果对比,验证了模型及有限元程序的有效性;水平井破裂压力数值分析结果与现场测试数据吻合较好。研究表明：射孔可调控水平井破裂压力与初始破裂位置,同时对近井筒区域裂缝扩展形态影响显著。通过优化射孔参数可以引导初始破裂向最优破裂面扩展、有效降低破裂压力,减小由于螺旋射孔空间排布引起的水平井近井筒裂缝迂曲与复杂程度,提高致密油气藏压裂改造效果。     

Modeling hydraulic fracture propagation in permeable media with an embedded strong discontinuity approach

The paper presents an embedded strong discontinuity approach to simulate single hydraulic fracture propagation in the poroelastic medium under plane-strain conditions. The method enriches the strain field with the discontinuous deformation mode and allows the fracture to be modeled inside elements. The Mode-I fracture initiation and propagation are described by the trilinear cohesive law, which is implemented by the penalty method. The enhanced permeability inside the fractured elements is dependent on the fracture aperture. Hydraulic fracture propagation is driven by the high pressure gradient near the fracture. Fluid transfer between the fracture and bulk rock is automatically captured within the poroelastic framework. The numerical framework is verified by the comparisons with the asymptotic analytical solutions for single hydraulic fracture propagation.     

泥岩隧道施工过程中渗流场与应力场全耦合损伤模型研究

在连续损伤力学理论基础上,将塑性损伤演化及渗流相互耦合的概念引入Mohr-Coulomb 破坏准则,用于分析在孔隙压力和塑性损伤演化共同作用下岩石损伤演化机制,建立了相应的有限元损伤数值分析模型,并应用于比利时核废料库开挖过程中泥岩隧道附近围岩发生损伤演化、渗流场和应力场耦合过程分析中,得到了开挖引起的围岩损伤特性、孔隙压力以及渗透性的变化规律,为进一步研究隧道流变过程水力耦合特性合理的数值计算模型建立方法提供基础。     

Thermodynamic‐based model for coupling temperature‐dependent viscoelastic,viscoplastic, and viscodamage constitutive behavior of asphalt mixtures

Based on the continuum damage mechanics, a general and comprehensive thermodynamic‐based framework for coupling the temperature‐dependent viscoelastic, viscoplastic, and viscodamage behaviors of bituminous materials is presented. This general framework derives systematically Schapery‐type nonlinear viscoelasticity, Perzyna‐type viscoplasticity, and a viscodamage model analogous to the Perzyna‐type viscoplasticity. The resulting constitutive equations are implemented in the well‐known finite element code Abaqus via the user material subroutine UMAT. A systematic procedure for identifying the model parameters is discussed. Finally, the model is validated by comparing the model predictions with a comprehensive set of experimental data on hot mix asphalt that include creep‐recovery, creep, uniaxial constant strain rate, and repeated creep‐recovery tests in both tension and compression over a range of temperatures, stress levels, and strain rates. Comparisons between model predictions and experimental measurements show that the presented constitutive model is capable of predicting the nonlinear behavior of asphaltic mixes under different loading conditions. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical modeling of dynamic rock fracture with a combined 3D continuum viscodamage‐embedded discontinuity model

This paper deals with numerical modeling of dynamic failure phenomena in rate‐sensitive quasi‐brittle materials, such as rocks, with initial microcrack populations. To this end, a continuum viscodamage‐embedded discontinuity model is developed and tested in full 3D setting. The model describes the pre‐peak nonlinear and rate‐sensitive hardening response of the material behavior, representing the fracture‐process zone creation, by a rate‐dependent continuum damage model. The post‐peak response, involving the macrocrack creation accompanied by exponential softening, is formulated by using an embedded displacement discontinuity model. The finite element implementation of this model relies upon the linear tetrahedral element, which seems appropriate for explicit dynamic analyses involving stress wave propagation. The problems of crack locking and spreading typical of embedded discontinuity models are addressed in this paper. A combination of two remedies, the inclusion of viscosity in the spirit of Wang's viscoplastic consistency approach and introduction of isotropic damaging into the embedded discontinuity model, is shown to be effective in the present explicit dynamics setting. The model performance is illustrated by several numerical simulations. In particular, the dynamic Brazilian disc test and the Kalthoff–Winkler experiment show that the present model provides realistic predictions with the correct failure modes and rate‐dependent tensile strengths of rock at different loading rates. The ability of initial embedded discontinuity populations to model the initial microcrack populations in rocks is also successfully tested. Copyright © 2016 John Wiley & Sons, Ltd.     

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

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.     

多场耦合作用下岩体裂隙扩展演化关键问题研究

裂隙岩体热-水-应力（THM）耦合是目前研究的热点和难点。首先总结了裂隙岩体多场耦合的机制、模型、方法及研究内容,并通过分析裂隙对THM耦合的重要控制作用,提出了在THM耦合中考虑裂隙网络扩展演化及模拟的关键问题,同时指出了研究的3个关键点：（1）建立考虑裂隙网络演化的耦合模型;（2）裂隙扩展的数值模拟方法;（3）THM耦合及岩体变形、失稳全过程的数值模拟算法。随后通过对模拟多场耦合和裂隙扩展数值方法的归类比较,重点论述了目前适用于模拟多场耦合下裂隙扩展模拟的各种数值方法（包括有限单元法、无单元法、单位分解法、离散单元法、岩石破裂过程分析方法和数值流形方法）的优缺点,并通过对比研究,推荐采用数值流形方法（NMM）来实现对关键问题的模拟研究。最后,对研究思路和难点进行了初步探讨。     

基于扩展有限元法的岩体水力劈裂耦合模型

1. 浙江水利水电学院 水利与环境工程学院,浙江 杭州 310018;2. 上海交通大学 船舶海洋与建筑工程学院,上海 200240; 3. 浙江省钱塘江管理局勘测设计院,浙江 杭州 310016     

A rate‐dependent cohesive crack model based on anisotropic damage coupled to plasticity

In quasi‐brittle material the complex process of decohesion between particles in microcracks and localization of the displacement field into macrocracks is limited to a narrow fracture zone, and it is often modelled with cohesive crack models. Since the anisotropic nature of the decohesion process in separation and sliding is essential, it is particularly focused in this paper. Moreover, for cyclic and dynamic loading the unloading, load reversal (including crack closure) and rate dependency are essential features that are included in a new model. The modelling of degradation is based on a ‘localized’ version of anisotropic continuum damage coupled to inelasticity. The concept of strain energy equivalence between the states in the effective and nominal settings is adopted in order to define the free energy of the interface. The proposed fracture criterion is of the Mohr type, with a smooth transition of the failure and kinematics (slip and dilatation) characteristics between tension and shear. The chosen potential, of the Lemaitre‐type, for evolution of the dissipative processes is additively decomposed into plastic and damage parts, and non‐associative constitutive equations are obtained. The constitutive equations are integrated by applying the backward Euler rule and by using Newton iteration. The proposed model is assessed analytically and numerically and a typical calibration procedure for concrete is proposed. Copyright © 2006 John Wiley & Sons, Ltd.     

Thermo‐hydro‐mechanical modeling of damage in unsaturated porous media: Theoretical framework and numerical study of the EDZ

The damage model presented in this article (named ‘THHMD’ model) is dedicated to non‐isothermal unsaturated porous media. It is formulated by means of three independent strain state variables, which are the thermodynamic conjugates of net stress, suction and thermal stress. The damage variable is a second‐order tensor. Stress/strain relationships are derived from Helmholtz free energy, which is assumed to be the sum of damaged elastic potentials and ‘crack‐closure energies’. Damage is assumed to grow with tensile strains due to net stress, with pore shrinkage due to suction and with thermal dilatation. Specific conductivities are introduced to account for the effects of cracking on the intensification and on the orientation of liquid water and vapor flows. These conductivities depend on damage and internal length parameters. The mechanical aspects of the THHMD model are validated by comparing the results of a triaxial compression test with experimental measurements found in the literature. Parametric studies of damage are performed on three different heating problems related to nuclear waste disposals. Several types of loading and boundary conditions are investigated. The thermal damage potential is thoroughly studied. The THHMD model is expected to be a useful tool in the assessment of the Excavation Damaged Zone, especially in the vicinity of nuclear waste repositories. Copyright © 2011 John Wiley & Sons, Ltd.