Overall moduli of solids with microcracks: Load-induced anisotropy

For a linearly elastic brittle solid containing microcracks that may be closed or may undergo frictional sliding, a general method is developed for estimating the overall instantaneous moduli which depend on the loading conditions. When the cracks are all open and when they are randomly distributed, then the overall response is isotropic. The moduli for this case have been obtained by Budiansky and O'C onnell (1976). On the other hand, when some cracks close, and when some closed cracks undergo frictional sliding, then the overall response becomes anisotropic and dependent on the loading conditions, as well as on the loading path. The self-consistent method is used to estimate the overall moduli. The effects of crack closure and loadinduced anisotropy are included. Several illustrative examples are worked out, showing the important influence of the load path on the overall response when crack closure and frictional sliding are involved.     

A micromechanics-based thermodynamic formulation of isotropic damage with unilateral and friction effects

This paper is devoted to micromechanical modeling of isotropic damage in brittle materials. The damaged materials will be considered as heterogeneous media composed of solid matrix weakened by isotropically distributed microcracks. The original contribution of the present work is to provide a proper micromechanical thermodynamic formulation for damage-friction modeling in brittle materials with the help of Eshelby’s solution to matrix-inclusion problems. The elastic and plastic strain energy involving unilateral effects will be fully determined. The condition of microcrack opening–closure transition will be determined in both strain-based and stress-based forms. The effect of spatial distribution of microcracks will also be taken into account. Further, the damage evolution law is formulated in a sound thermodynamic framework and inherently coupled with frictional sliding. As a first phase of validation, the proposed micromechanical model is finally applied to reproduce basic mechanical responses of ordinary concrete in compression tests.     

Dynamic damage constitutive relation of mesoscopic heterogenous brittle rock under rotation of principal stress axes

A micro-mechanics-based model is developed to investigate microcrack damage mechanism of four stages of brittle rock under rotation of the principal stress axes. They consist of linear elastic, non-linear hardening, rapid stress drop and strain softening. The frictional sliding crack model is applied to analyze microcracks nucleation, propagation and coalescence. The strain energy density factor approach is applied to determine the critical condition of microcrack nucleation, propagation and coalescence. The inelastic strain increments are formulated within the framework of thermodynamics with internal variables. Rotation of principal stress axes affect the dynamic damage constitutive relationship and the failure strength of brittle rock.     

An isotropic unilateral damage model coupled with frictional sliding for quasi-brittle materials

In this paper, we present an original extension of an isotropic damage model for quasi-brittle materials and assess its predictive capabilities. The proposed model accounts not only for unilateral behavior related to the opening and closure of microcracks but also for inelastic strains reflecting the frictional sliding along closed microcracks. More importantly, owing to its careful mathematical formulation with a particular attention paid to the continuous differentiability of the underlined thermodynamic potential, the model ensures the continuity of the inelastic stress–strain response. First applications show that it is able to predict the asymmetric behavior and hysteretic response of microcracked materials such as concrete and some rocks     

Micromechanical analysis of coupling between anisotropic damage and friction in quasi brittle materials: Role of the homogenization scheme

This paper deals with micromechanical analysis of anisotropic damage and its coupling with friction in quasi brittle materials. The anisotropic model is formulated in the framework of Eshelby-based homogenization methods. The emphasis is put on the study of effects of spatial distribution of microcracks and their interactions. Microcracks closure effects as well as coupling between damage evolution and frictional sliding on closed cracks lips are taken into account. The interaction of sliding and damage evolution is addressed by performing a global thermodynamic analysis on two macroscopic criteria established in the paper. The role of the homogenization scheme is discussed in detail through various applications.     

Analytical and numerical analysis of frictional damage in quasi brittle materials

Frictional sliding and crack growth are two main dissipation processes in quasi brittle materials. The frictional sliding along closed cracks is the origin of macroscopic plastic deformation while the crack growth induces a material damage. The main difficulty of modeling is to consider the inherent coupling between these two processes. Various models and associated numerical algorithms have been proposed. But there are so far no analytical solutions even for simple loading paths for the validation of such algorithms. In this paper, we first present a micro-mechanical model taking into account the damage-friction coupling for a large class of quasi brittle materials. The model is formulated by combining a linear homogenization procedure with the Mori–Tanaka scheme and the irreversible thermodynamics framework. As an original contribution, a series of analytical solutions of stress–strain relations are developed for various loading paths. Based on the micro-mechanical model, two numerical integration algorithms are exploited. The first one involves a coupled friction/damage correction scheme, which is consistent with the coupling nature of the constitutive model. The second one contains a friction/damage decoupling scheme with two consecutive steps: the friction correction followed by the damage correction. With the analytical solutions as reference results, the two algorithms are assessed through a series of numerical tests. It is found that the decoupling correction scheme is efficient to guarantee a systematic numerical convergence.     

Mechanical Behaviors of Saturated Sand under Complicated Loading

The different physical states of saturated sand, including shear elasticity, positive dilatancy, and negative dilatancy (preliminary negative dilatancy, secondary negative dilatancy and reversal negative dilatancy) are revealed based on the pore water pressure response of saturated sand in undrained dynamic torsional tests of thin cylinder samples and also checked by the drained cyclic triaxial tests under a given mean effective normal stress. According to the effective stress path of different physical states under the undrained cyclic torsional tests the physical state transformation surface, stress history boundary and yield surface are determined, and the state boundary surface is also determined by the range of effective frictional stress state movement. Based on the moving yield surface without rotation, and the expanding stress history boundary surface relevant to the stress path variations under different physical states in 3D stress space, a physical state model is proposed to provide a new app     

Partial sliding along a planar crack weakened by the presence of fluid

Partial sliding along a infinite planar crack that is locally weakened by the presence of fluid is analyzed. Under uniformly applied normal and shear loads the crack surface generates non-uniform frictional resistance that has a local minimum within a penny-shaped fluid-filled domain. After the applied shear load reaches the resistance minimum a circular sliding zone initiates and then spreads around as the applied shear load is gradually increased. The primary focus of this work is to analyze sliding evolution as a function of the applied loads and the induced fluid pressure. The growth of the sliding zone is studied based on the criterion that shear stress intensity factors are zero along the zone boundary. An analytical relation between the radius of the sliding zone and the applied shear load is determined.     

An interacting micro-crack damage model for failure of brittle materials under compression

A model is developed for brittle failure under compressive loading with an explicit accounting of micro-crack interactions. The model incorporates a pre-existing flaw distribution in the material. The macroscopic inelastic deformation is assumed to be due to the nucleation and growth of tensile “wing” micro-cracks associated with frictional sliding on these flaws. Interactions among the cracks are modeled by means of a crack-matrix-effective-medium approach in which each crack experiences a stress field different from that acting on isolated cracks. This yields an effective stress intensity factor at the crack tips which is utilized in the formulation of the crack growth dynamics. Load-induced damage in the material is defined in terms of a scalar crack density parameter, the evolution of which is a function of the existing flaw distribution and the crack growth dynamics. This methodology is applied for the case of uniaxial compression under constant strain rate loading. The model provides a natural prediction of a peak stress (defined as the compressive strength of the material) and also of a transition strain rate, beyond which the compressive strength increases dramatically with the imposed strain rate. The influences of the crack growth dynamics, the initial flaw distribution, and the imposed strain rate on the constitutive response and the damage evolution are studied. It is shown that different characteristics of the flaw distribution are dominant at different imposed strain rates: at low rates the spread of the distribution is critical, while at high strain rates the total flaw density is critical.     

Impact initiation of explosives and propellants via statistical crack mechanics

A statistical approach has been developed for modeling the dynamic response of brittle materials by superimposing the effects of a myriad of microcracks, including opening, shear, growth and coalescence, taking as a starting point the well-established theory of penny-shaped cracks. This paper discusses the general approach, but in particular an application to the sensitivity of explosives and propellants, which often contain brittle constituents. We examine the hypothesis that the intense heating by frictional sliding between the faces of a closed crack during unstable growth can form a hot spot, causing localized melting, ignition, and fast burn of the reactive material adjacent to the crack. Opening and growth of a closed crack due to the pressure of burned gases inside the crack and interactions of adjacent cracks can lead to violent reaction, with detonation as a possible consequence.This approach was used to model a multiple-shock experiment by Mulford et al. [1993. Initiation of preshocked high explosives PBX-9404, PBX-9502, PBX-9501, monitored with in-material magnetic gauging. In: Proceedings of the 10th International Detonation Symposium, pp. 459-467] involving initiation and subsequent quenching of chemical reactions in a slab of PBX 9501 impacted by a two-material flyer plate. We examine the effects of crack orientation and temperature dependence of viscosity of the melt on the response. Numerical results confirm our theoretical finding [Zuo, Q.H., Dienes, J.K., 2005. On the stability of penny-shaped cracks with friction: the five types of brittle behavior. Int. J. Solids Struct. 42, 1309-1326] that crack orientation has a significant effect on brittle behavior, especially under compressive loading where interfacial friction plays an important role. With a reasonable choice of crack orientation and a temperature-dependent viscosity obtained from molecular dynamics calculations, the calculated particle velocities compare well with those measured using embedded velocity gauges.     

A fracture damage model for Jointed rock masses and its application to the stability analysis of dam abutments

Based on continuum damage mechanics, for jointed rock masses, a fracture damage model is presented in this paper. First, the damage tensors are defined through the elastic-flexibility of intact rock and the equivalent elastic-damage flexibility for rock mass. Then, by the self-consistent principle of solid mechanics, the equivalent elastic-damage flexibility tensors involving the interaction between multicracks are deduced. The damage evolution law is proposed involving the mechanism of crack propagation process: frictional sliding, crack kinking, growing of branched tension cracks, interlinking of the microcracks near branched crack tips leading to the breakthrough of macro-cracks and finally the failure of rock mass. Thus the evolution of damage variables reasonably unified with the process of crack propagation is given. Finally, a plastic-brittle damage constitutive relation including brittle coupled strain rate, developed and applied to the stability analysis of complicated rock foundation of a dam in China, is described in this paper.     

Microscopic analysis of crack propagation for multiple cracks, inclusions and voids

The elastic crack interaction with internal defects, such as microcracks, voids and rigid inclusions, is investigated in this study for the purpose of analyzing crack propagation. The elastic stress field is obtained using linear theory of elasticity for isotropic materials. The cracks are modeled as pile-ups of edge dislocations resulting into a coupled set of integral equations, whose kernels are those of a dislocation in a medium with or without an inclusion or void. The numerical solution of these equations gives the stress intensity factors and the complete stress field in the given domain. The solution is valid for a general solid, however the propagation analysis is valid mostly for brittle materials. Among different propagation models the ones based on maximum circumferential stress and minimum strain energy density theories, are employed. A special emphasis is given to the estimation of the crack propagation direction that defines the direction of crack branching or kinking. Once a propagation direction is determined, an improved model dealing with kinked cracks must be employed to follow the propagation behavior.     

A microstructural elastoplastic model for unsaturated granular materials

The homogenization technique is used to obtain an elastoplastic stress–strain relationship for dry, saturated and unsaturated granular materials. Deformation of a representative volume of material is generated by mobilizing particle contacts in all orientations. In this way, the stress–strain relationship can be derived as an average of the mobilization behavior of these local contact planes. The local behavior is assumed to follow a Hertz–Mindlin’s elastic law and a Mohr–Coulomb’s plastic law. For the non-saturated state, capillary forces at the grain contacts are added to the contact forces created by an external load. They are calculated as a function of the degree of saturation, depending on the grain size distribution and on the void ratio of the granular assembly. Numerical simulations show that the model is capable of reproducing the major trends of a partially saturated granular assembly under various stress and water content conditions. The model predictions are compared to experimental results on saturated and unsaturated samples of silty sands under undrained triaxial loading condition. This comparison shows that the model is able to account for the influence of capillary forces on the stress–strain response of the granular materials and therefore, to reproduce the overall mechanical behavior of unsaturated granular materials.     

Micromechanical approach to the behavior of poroelastic materials

A micromechanical analysis of the behavior of saturated elastic porous media is presented. Pores are connected and saturated with a fluid at uniform pressure. The non-linear macroscopic response of the porous medium is assumed to be related to the existence of a network of microcracks in the solid phase. The progressive closure of cracks during loading is viewed to control the material non-linearity.Using averaging schemes based on both the classical and modified secant methods, the non-linearity of the overall response is shown to be sometimes controlled by an effective stress. The modified secant method accounts in a certain way for the non-uniform distribution of strains within the solid phase. Different definitions of the effective stress are found, each of them being a combination of the overall stress with the fluid pressure. The hypothesis of connexion between pores and microcracks is discussed in details. The consequences of this hypothesis on the non-linear overall response are illustrated and comparisons are made with experimental data.     

不同温度条件下含流体砂岩裂隙摩擦特性的试验研究

砂岩储层中由流体注入导致的地震活动与砂岩断层(裂隙)的摩擦行为有关.为了揭示不同温度条件下含流体砂岩裂隙的摩擦特性,在温度范围为25°C～140°C和有效法向应力范围为4～12 MPa的试验条件下,本文分别对干燥、水饱和以及注CO₂锯切砂岩裂隙进行了速度分级加载试验.试验结果表明:(1)对于干燥砂岩裂隙,增大有效法向应力和升高温度均能增大裂隙的初始摩擦系数,而改变有效法向应力对裂隙摩擦稳定性影响不明显,仅升高温度会略微降低其摩擦稳定性;(2)对于水饱和砂岩裂隙,裂隙的初始摩擦系数同样会随着有效法向应力的增大而增大,但会受到升温的弱化作用,而增大有效法向应力和升高温度均能降低裂隙的摩擦稳定性;(3)对于注CO₂砂岩裂隙,裂隙的初始摩擦系数受有效法向应力和温度变化的影响与水饱和砂岩裂隙相反,但裂隙的摩擦稳定性仅会随着温度的升高而降低,受有效法向应力的影响不明显.因此,砂岩裂隙的摩擦特性受有效法向应力、温度和注入流体类型的共同影响.该试验结果对理解流体注入诱发地震有一定的指示作用.     

Boundary-layer analysis of propagating mode II cracks in porous elastic media

The steady advance of semi-infinite cracks in porous elastic media is examined in the limiting cases of very slow and very rapid propagation. The low-speed limit is the elastic solution representing fully-drained behavior, with no boundary layers present. In the high-speed limit, two boundary layers are identified: one along the crack face and another centered at the crack tip. The zero-order outer solution is the elastic solution representing fully-undrained behavior, with both boundary layers acting as zones of transition to the effectively drained conditions at the crack itself. Closedform expressions are derived for the stresses and pore pressure within the crack-tip boundary layer. The true stress intensity factor is found to be reduced by an amount related to the disparity between the drained and undrained Poisson's ratios of the medium.     

Effective constitutive behavior of nonlinear solids containing penny-shaped cracks

This article studies the effective constitutive relations of brittle and ductile solids containing either aligned or randomly oriented penny-shaped cracks. The matrix material considered is linear before yielding and behaves with a pure power-law after yielding. An explicit yield function is obtained for the cracked body under general external loading conditions, and the effective stress potential is given in a simple optimization procedure. The results show that the presence of microcracks inside solids significantly lowers the strength of materials.     

Material instabilities in inelastic saturated porous media under dynamic loadings

Instabilities in inelastic saturated porous media are investigated here for general three-dimensional states under dynamic loadings using a perturbation approach.Under quasi-static conditions, unbounded growth of perturbations is related to the emergence of stationary discontinuities under drained or undrained conditions, while under dynamic conditions, unbounded growth is related either to the emergence of stationary discontinuities (and these are set by drained conditions) or to the appearance of the flutter phenomenon (acceleration waves).For associative behaviour unbounded growth always corresponds to localization under drained conditions and the onset of growth of perturbations occurs here only through divergence growth. It is only for non-associative flow that unbounded growth may correspond to undrained localization in quasi-static conditions and to flutter under dynamic conditions.     

芳纶纤维增强酚醛树脂摩擦材料的磨损机理研究

对芳纶增强酚醛树脂摩擦材料在模拟制动工况下的磨损机理进行了系统的实验研究,结果表明：在试验工况下,两种摩擦材料在与铸铁偶件对摩时均呈现粘着磨损和塑性变形特征;在高速和高压条件下塑性变形加剧,摩擦材料磨损表面可见熔融迹象;摩擦副接触表面发生材料的相互转移,两种摩擦材料均可在偶件表面形成转移膜,且在高速和高压条件下转移膜更易形成;Kevlar纤维作为增强相可以有效地提高摩擦材料的摩擦稳定性,并降低摩擦系数;摩擦材料的磨损机理主要为擦伤作用、粘着磨损和塑性变形。     

Thermal fracture of macrocrack with closure as influenced by microcracks

The fracture stability of macrocracks under uniform heat flux is analyzed to include the effect of a system of microcracks. The interaction of cracks leads to full or partial closure of crack surfaces. The boundary problem is stated and a solution is obtained. The domains where microcracks are closed and/or affect partial closing of the macrocrack are found. Evaluated is the macrocrack tip stress intensity factor accounting for closure.