Anisotropic damage mechanics as a novel approach to improve pre-and post-failure borehole stability analysis

Anisotropic damage distribution and evolution have a profound effect on borehole stress concentrations. Damage evolution is an irreversible process that is not adequately described within classical equilibrium thermodynamics. Therefore, we propose a constitutive model, based on non-equilibrium thermodynamics, that accounts for anisotropic damage distribution, anisotropic damage threshold and anisotropic damage evolution. We implemented this constitutive model numerically, using the finite element method, to calculate stress–strain curves and borehole stresses. The resulting stress–strain curves are distinctively different from linear elastic-brittle and linear elastic-ideal plastic constitutive models and realistically model experimental responses of brittle rocks. We show that the onset of damage evolution leads to an inhomogeneous redistribution of material properties and stresses along the borehole wall. The classical linear elastic-brittle approach to borehole stability analysis systematically overestimates the stress concentrations on the borehole wall, because dissipative strain-softening is underestimated. The proposed damage mechanics approach explicitly models dissipative behaviour and leads to non-conservative mud window estimations. Furthermore, anisotropic rocks with preferential planes of failure, like shales, can be addressed with our model.

Study Of The Damage-Induced Anisotropy Of Quasi-Brittle Materials Using The Component Assembling Model

Damage-induced anisotropy of quasi-brittle materials is investigated using component assembling model in this study. Damage-induced anisotropy is one significant character of quasi-brittle materials coupled with nonlinearity and strain softening. Formulation of such complicated phenomena is a difficult problem till now. The present model is based on the component assembling concept, where constitutive equations of materials are formed by means of assembling two kinds of components' response functions. These two kinds of components, orientational and volumetric ones, are abstracted based on pair-functional potentials and the Cauchy - Born rule. Moreover, macroscopic damage of quasi-brittle materials can be reflected by stiffness changing of orientational components, which represent grouped atomic bonds along discrete directions. Simultaneously, anisotropic characters are captured by the naturally directional property of the orientational component. Initial damage surface in the axial-shear stress space is calculated and analyzed. Furthermore, the anisotropic quasi-brittle damage behaviors of concrete under uniaxial, proportional, and nonproportional combined loading are analyzed to elucidate the utility and limitations of the present damage model. The numerical results show good agreement with the experimental data and predicted results of the classical anisotropic damage models.

Study of the damage-induced anisotropy of quasi-brittle materials using the component assembling model

Damage-induced anisotropy of quasi-brittle materials is investigated using component assembling model in this study. Damage-induced anisotropy is one significant character of quasi-brittle materials coupled with nonlinearity and strain softening. Formulation of such complicated phenomena is a difficult problem till now. The present model is based on the component assembling concept, where constitutive equations of materials are formed by means of assembling two kinds of components' response functions. These two kinds of components, orientational and volumetric ones, are abstracted based on pair-functional potentials and the Cauchy - Born rule. Moreover, macroscopic damage of quasi-brittle materials can be reflected by stiffness changing of orientational components, which represent grouped atomic bonds along discrete directions. Simultaneously, anisotropic characters are captured by the naturally directional property of the orientational component. Initial damage surface in the axial-shear stress space is calculated and analyzed. Furthermore, the anisotropic quasi-brittle damage behaviors of concrete under uniaxial, proportional, and nonproportional combined loading are analyzed to elucidate the utility and limitations of the present damage model. The numerical results show good agreement with the experimental data and predicted results of the classical anisotropic damage models.

Anisotropic characterization of crack growth in the tertiary flow of asphalt mixtures in compression

Asphalt mixtures exhibit primary, secondary, and tertiary stages in sequence during a rutting deterioration. Many field asphalt pavements are still in service even when the asphalt layer is in the tertiary stage, and rehabilitation is not performed until a significant amount of rutting accompanied by numerous macrocracks is observed. The objective of this study was to provide a mechanistic method to model the anisotropic cracking of the asphalt mixtures in compression during the tertiary stage of rutting. Laboratory tests including nondestructive and destructive tests were performed to obtain the viscoelastic and viscofracture properties of the asphalt mixtures. Each of the measured axial and radial total strains in the destructive tests were decomposed into elastic, plastic, viscoelastic, viscoplastic, and viscofracture strains using the pseudostrain method in an extended elastic-viscoelastic correspondence principle. The viscofracture strains are caused by the crack growth, which is primarily signaled by the increase of phase angle in the tertiary flow. The viscofracture properties are characterized using the anisotropic damage densities (i.e., the ratio of the lost area caused by cracks to the original total area in orthogonal directions). Using the decomposed axial and radial viscofracture strains, the axial and radial damage densities were determined by using a dissipated pseudostrain energy balance principle and a geometric analysis of the cracks, respectively. Anisotropic pseudo J-integral Paris' laws in terms of damage densities were used to characterize the evolution of the cracks in compression. The material constants in the Paris' law are determined and found to be highly correlated. These tests, analysis, and modeling were performed on different asphalt mixtures with two binders, two air void contents, and three aging periods. Consistent results were obtained; for instance, a stiffer asphalt mixture is demonstrated to have a higher modulus, a lower phase angle, a greater flow number, and a larger n1 value (exponent of Paris' law). The calculation of the orientation of cracks demonstrates that the asphalt mixture with 4% air voids has a brittle fracture and a splitting crack mode, whereas the asphalt mixture with 7% air voids tends to have a ductile fracture and a diagonal sliding crack mode. Cracks of the asphalt mixtures in compression are inclined to propagate along the direction of the external compressive load. © 2014 American Society of Civil Engineers.

Interpretation of the Superpave IDT strength test using a viscoelastic-damage constitutive model

This paper presents a new interpretation for the Superpave IDT strength test based on a viscoelastic-damage framework. The framework is based on continuum damage mechanics and the thermodynamics of irreversible processes with an anisotropic damage representation. The new approach introduces considerations for the viscoelastic effects and the damage accumulation that accompanies the fracture process in the interpretation of the Superpave IDT strength test for the identification of the Dissipated Creep Strain Energy (DCSE) limit from the test result. The viscoelastic model is implemented in a Finite Element Method (FEM) program for the simulation of the Superpave IDT strength test. The DCSE values obtained using the new approach is compared with the values obtained using the conventional approach to evaluate the validity of the assumptions made in the conventional interpretation of the test results. The result shows that the conventional approach over-estimates the DCSE value with increasing estimation error at higher deformation rates.

3D Anisotropic Elastoplastic-Damage Model and its Application in Simulating the Behavior of Rock Materials

A 3D anisotropic elastoplastic-damage model was presented based on continuum damage mechanics theory. In this model, the tensor decomposition technique is employed. Combined with the plastic yield rule and damage evolution, the stress tensor in incremental format is obtained. The derivate eigenmodes in the proposed model are assumed to be related with the uniaxial behavior of the rock material. Each eigenmode has a corresponding damage variable due to the fact that damage is a function of the magnitude of the eigenstrain. Within an eigenmodes, different damage evolution can be used for tensile and compressive loadings. This model was also developed into finite element code in explicit format, and the code was integrated into the well-known computational environment ABAQUS using the ABAQUS/Explicit Solver. Numerical simulation of an uniaxial compressive test for a rock sample is used to examine the performance of the proposed model, and the progressive failure process of the rock sample is unveiled.

An Anisotropic Model of Damage and Unilateral Effect for Brittle Materials

This work deals with the initial applications and formulation of an aniscitropic plastic-damage constitutive model proposed for non-linear analysis of reinforced concrete structures submitted to a loading with change of the sign. The original constitutive model is based on the fundamental hypothesis of energy equivalence between real and continuous medium following the concepts of the Continuum Damage Mechanics. The concrete is assumed as an initial elastic isotropic medium presenting anisotropy, permanent strains and bimodularity (distinct elastic responses whether traction or compression stress states prevail) induced by damage evolution. In order to take into account the bimodularity, two damage tensors governing the rigidity in tension or compression regimes are introduced. Then, some conditions are introduced in the original version of the model in order to simulate the damage unilateral effect. The three-dimensional version of the proposed model is analyzed in order to validate its formulation when compared to micromechanical theory. The one-dimensional version of the model is applied in the analyses of a reinforced concrete beam submitted to a loading with change of the sign. Despite the parametric identification problems, the initial applications show the good performance of the model.

An anisotropic elastic-viscoplastic damage model for bone tissue

A new anisotropic elastic-viscoplastic damage constitutive model for bone is proposed using an eccentric elliptical yield criterion and nonlinear isotropic hardening. A micromechanics-based multiscale homogenization scheme proposed by Reisinger et al. is used to obtain the effective elastic properties of lamellar bone. The dissipative process in bone is modeled as viscoplastic deformation coupled to damage. The model is based on an orthotropic ecuntric elliptical criterion in stress space. In order to simplify material identification, an eccentric elliptical isotropic yield surface was defined in strain space, which is transformed to a stress-based criterion by means of the damaged compliance tensor. Viscoplasticity is implemented by means of the continuous Perzyna formulation. Damage is modeled by a scalar function of the accumulated plastic strain D(κ) , reducing all element s of the stiffness matrix. A polynomial flow rule is proposed in order to capture the rate-dependent post-yield behavior of lamellar bone. A numerical algorithm to perform the back projection on the rate-dependent yield surface has been developed and implemented in the commercial finite element solver Abaqus/Standard as a user subroutine UMAT. A consistent tangent operator has been derived and implemented in order to ensure quadratic convergence. Correct implementation of the algorithm, convergence, and accuracy of the tangent operator was tested by means of strain- and stress-based single element tests. A finite element simulation of nano- indentation in lamellar bone was finally performed in order to show the abilities of the newly developed constitutive model.

Non-Destructive Evaluation of Defects and Damage in Composite Materials and Structures

The mechanical behaviour of composite materials differs from that of conventional structural materials owing to their heterogeneous and anisotropic nature. Different types of defects and anomalies get induced in these materials during the fabrication process. Further, during their service life, the components made of composite materials develop different types of damage. The performance and life of such components is governed by the combined effect of all these defects and damage. While porosity, voids, inclusions etc., are some defects those can get induced during the fabrication of composites, matrix cracks, interface debonds, delaminations and fiber breakage are major types of service induced damage which are of concern. During the service life of components made of composites, one type of damage can grow and initiate another type of damage. For example, matrix cracks can gradually grow to the interface and initiate debonds. Interface debonds in a particular plane can lead to delaminations. Consequently, the combined effect of different types of distributed damage causes the failure of the component. A set of non-destructive evaluation (NDE) methods is well established for testing conventional metallic materials. Some of them can also be utilized for composite materials as they are, and in some cases with a little different approach or modification. Ultrasonics, Radiography, Thermography, Fiber Optics, Acoustic Emision Techniques etc., to name a few. Detection, evaluation and characterization of different types of defects and damage encountered in composite materials and structures using different NDE tools is discussed briefly in this paper.

Closed-form nonlinear sectional analysis of pretwisted anisotropic beam with multiple delaminations

Fiber reinforced laminated composite open-section beams are widely used as bearingless rotor flex beams because of their high specific strength and stiffness as well as fatigue life. These laminated composite structures exhibit a number of different failure modes, including fiber-matrix debonding within individual layers, delamination or separation of the layers, transverse cracks through one or more layers and fiber fracture. Delamination is a predominant failure mode in continuous fiber reinforced laminated composites and often initiate near the free edges of the structure. The appearance of delaminations in the composite rotorcraft flexbeams can lead to deterioration of the mechanical properties and, in turn, the helicopter performance as well as safety. Understanding and predicting the influence of free-edge delamination on the overall behavior of the laminates will provide quantitative measures of the extent of the damage and help ensure their damage tolerance.

Relationship between DNA damage and sperm head birefringence

Birefringence or double refraction is the decomposition of a ray of light into two rays when it passes through an anisotropic material such as quartz. Sperm cells have been demonstrated to be optically anisotropic. The objective of this study was to evaluate the relationship between the pattern of human sperm head birefringence (SHBF) and DNA damage. A total of 26 patients with normal semen were included. DNA damage (fragmentation and denaturation) was evaluated in the sperm head in the context of birefringence, both total (SHBF-T) and partial (SHBF-P), by terminal deoxyribonucleotidyl transferase (TdT)-mediated dUDP nick-end labelling assay and acridine orange fluorescence, respectively. Positive DNA fragmentation in spermatozoa with SHBF-T (205/1053; 19.5%) was significantly higher (P < 0.0001) than in spermatozoa that presented SHBF-P (60/820; 7.3%). However, the percentage of denatured DNA in spermatozoa with SHBF-T (824/1256; 65.6%) was not significantly different from the ones with SHBF-P (666/1009; 66.0%). In conclusion, the data support a positive relationship between spermatozoa with total SHBF in their head and increased DNA fragmentation. (C) 2011, Reproductive Healthcare Ltd. Published by Elsevier Ltd. All rights reserved.

A generalized anisotropic quadric yield criterion and its application to bone tissue at multiple length scales

Nonlinear computational analysis of materials showing elasto-plasticity or damage relies on knowledge of their yield behavior and strengths under complex stress states. In this work, a generalized anisotropic quadric yield criterion is proposed that is homogeneous of degree one and takes a convex quadric shape with a smooth transition from ellipsoidal to cylindrical or conical surfaces. If in the case of material identification, the shape of the yield function is not known a priori, a minimization using the quadric criterion will result in the optimal shape among the convex quadrics. The convexity limits of the criterion and the transition points between the different shapes are identified. Several special cases of the criterion for distinct material symmetries such as isotropy, cubic symmetry, fabric-based orthotropy and general orthotropy are presented and discussed. The generality of the formulation is demonstrated by showing its degeneration to several classical yield surfaces like the von Mises, Drucker–Prager, Tsai–Wu, Liu, generalized Hill and classical Hill criteria under appropriate conditions. Applicability of the formulation for micromechanical analyses was shown by transformation of a criterion for porous cohesive-frictional materials by Maghous et al. In order to demonstrate the advantages of the generalized formulation, bone is chosen as an example material, since it features yield envelopes with different shapes depending on the considered length scale. A fabric- and density-based quadric criterion for the description of homogenized material behavior of trabecular bone is identified from uniaxial, multiaxial and torsional experimental data. Also, a fabric- and density-based Tsai–Wu yield criterion for homogenized trabecular bone from in silico data is converted to an equivalent quadric criterion by introduction of a transformation of the interaction parameters. Finally, a quadric yield criterion for lamellar bone at the microscale is identified from a nanoindentation study reported in the literature, thus demonstrating the applicability of the generalized formulation to the description of the yield envelope of bone at multiple length scales.

Damage tolerance of cold drawn ferritic-austenitic stainless steels wires for prestressed concrete

Damage tolerance of high strength cold-drawn ferritic–austenitic stainless steel wires is assessed by means of tensile fracture tests of cracked wires. The fatigue crack is transversally propagated from the wire surface. The damage tolerance curve of the wires results from the empirical failure load when given as a function of crack depth. As a consequence of cold drawing, the wire microstructure is orientated along its longitudinal axis and anisotropic fracture behaviour is found at macrostructural level at the tensile failure of the cracked specimens. An in situ optical technique known as video image correlation VIC-2D is used to get an insight into this failure mechanism by tensile testing transversally fatigue cracked plane specimens extracted from the cold-drawn wires. Finally, the experimentally obtained damage tolerance curve of the cold-drawn ferritic–austenitic stainless steel wires is compared with that of an elementary plastic collapse model and existing data of two types of high strength eutectoid steel currently used as prestressing steel for concrete.

A gradient-enhanced large-deformation continuum damage model for fibre-reinforced materials

A non-local gradient-based damage formulation within a geometrically non-linear setting is presented. The hyperelastic constitutive response at local material point level is governed by a strain energy which is additively composed of an isotropic matrix and of an anisotropic fibre-reinforced material, respectively. The inelastic constitutive response is governed by a scalar [1–d]-type damage formulation, where only the anisotropic elastic part is assumed to be affected by the damage. Following the concept in Dimitrijević and Hackl [28], the local free energy function is enhanced by a gradient-term. This term essentially contains the gradient of the non-local damage variable which, itself, is introduced as an additional independent variable. In order to guarantee the equivalence between the local and non-local damage variable, a penalisation term is incorporated within the free energy function. Based on the principle of minimum total potential energy, a coupled system of Euler–Lagrange equations, i.e., the balance of linear momentum and the balance of the non-local damage field, is obtained and solved in weak form. The resulting coupled, highly non-linear system of equations is symmetric and can conveniently be solved by a standard incremental-iterative Newton–Raphson-type solution scheme. Several three-dimensional displacement- and force-driven boundary value problems—partially motivated by biomechanical application—highlight the mesh-objective characteristics and constitutive properties of the model and illustratively underline the capabilities of the formulation proposed

A gradient-enhanced continuum damage model with application to fibre-reinforced tissues at finite strains

A non-local gradient-based damage formulation within a geometrically non-linear set- ting is presented. The hyperelastic constitutive response at local material point level is governed by a strain energy function which is additively composed by an isotropic neo-Hookean matrix and by an anisotropic fibre-reinforced material based on the model proposed by T. Gasser, R. Ogden, and G. Holzapfel.