Energy dissipation and contour integral characterizing fracture behavior of incremental plasticity

Jep -integral is derived for characterizing the frac- ture behavior of elastic-plastic materials. The J ep -integral differs from Rice’s J-integral in that the free energy density rather than the stress working density is employed to define energy-momentum tensor. The J ep -integral is proved to be path-dependent regardless of incremental plasticity and deformation plasticity. The J epintegral possesses clearly clear physical meaning: (1) the value J ep tip evaluated on the infinitely small contour surrounding the crack tip represents the crack tip energy dissipation; (2) when the global steadystate crack growth condition is approached, the value of J ep farss calculated along the boundary contour equals to the sum of crack tip dissipation and bulk dissipation of plastic zone. The theoretical results are verified by simulating mode I crack problems.     

A novel method for measuring and characterizing dynamic fracture-initiation toughness of elastic-plastic materials

The characterization and testing methods of the dynamic fracture initiation toughness of elastic-plastic materials under tensile impact are studied. By using the self-designed bar-bar tensile impact apparatus, a novel test method for studying dynamic fracture-initiation has been proposed based on the one-dimensional test principle. The curve of average loadv. s. displacement is smooth until unstable crack propagation, and the kinetic energy which does not contribute to the crack growth can be removed from total work done by external-force to the specimen. The fracture-initiation point is determined by compliance-changing rate method. The results show that these methods are feasible and effective. Through the analysis of the conversion between work and energy of a fracture specimen, the dynamicJ-integral is adopted as a characteristic parameter for elastic-plastic materials under impact loading. TheJ-integral is calculated from and curves by using the formula proposed, by Rice. TheJ-integral at fracture initiation is employed to describe the dynamic fracture-initiation toughness of elastic-plastic materials and the experimental results indicate thatJ _ID can be regarded as a material constant.     

Computer simulated analyses on deformation and fracture of non-cracked and cracked specimens

A unified damage and fracture model, the combinatory work density model, which is suitable for either non-cracked body or cracked body has been suggested^[t−7]. In the present paper, the deformation and fracture of the two kinds of tensile spceimen and TPB specimen made of 40Cr steel have been simulated by using the new model together with the large elastic-plastic deformation finite element method. The results give a good picture of the whole deformation and fracture processes of the specimens in experiments; especially, the results on the TPB specimen can be used to obtain the relationship between load and displacement at the loading pointP-Δ, and between crack extension and displacement at the loading point Δa-Δ, the resistance curveJ _R -Δa and the fracture toughnessJ _1C . All the results are in remarkable agreement with those obtained by experiments. Therefore the model suggested here can be used to simulate crack initiation and propagation in non-cracked body and fracture initiation and crack stable propagation in cracked body. The project supported by National Natural Science Foundation of China     

In Plane Displacement Formulation for Finite Cracked Plates Under Mode I Using Grid Method and Finite Element Analysis

A grid method is used to experimentally determine the in-plane displacement fields around a crack tip in a Single-Edge-Notch (SEN) tensile polyurethane specimen. Horizontal displacement u _x-exp and vertical displacement u _y-exp are expressed as functions of circular coordinates centred on the crack tip. These are compared with the approximate solutions of linear elastic fracture mechanics with a view to studying the applicability to polymers. The results show that this solution is not in agreement with the experiments at the focused on the vicinity of a crack tip. Taking this into account, an FEA program is developed with CAST3M for the purpose of comparing the experimental displacements and the numerical data. New formulations of displacements u _x and u _y are then developed. These formulations are derived from the principle of superposition and based on Arakawa’s formulation. With the displacement gradients obtained from the FEA and the new formulations, the determination of J-integrals is found to be in very good agreement with those derived from numerical calculation. Consequently, the proposed formulations can give displacement fields compatible with the J-integral calculation for the region near the crack tip. An application based on an experimental test is proposed to evaluate the performances of the proposed formulations.     

Fracture toughness evaluation of natural rubber

A method for estimating the fracture toughness of rubber-like materials is presented. Experimental data of a notched natural rubber sheet is analysed by application of the path-independent J-integral. A finite element code for large elastic deformations is used to evaluate the deformed shape of the rubber at crack growth initiation.Discussed is a stability criterion based on the existence of a critical value of J in relation with the experimental results.     

Critical J-integral and tearing energies for fracture of reinforced natural rubber

The critical tearing energy and J-integral for initiation and rapid crack propagation at 22°C were determined for carbon black reinforced natural rubber. Tearing energies were 20 per cent and up to 100 per cent greater for initiation and rupture, respectively, and more scattered than the J values. Specimen geometry affects both fracture characterizing parameters because of the energy dissipated during deformation. However, the J-integral analysis can partition this energy between the crack tip region and the bulk of the specimen. Therefore, application of the J-integral concept for characterizing fracture of elastomers is promising.     

The surface-forming energy release rate based fracture criterion for elastic–plastic crack propagation

The J-integral based criterion is widely used in elastic–plastic fracture mechanics. However, it is not rigorously applicable when plastic unloading appears during crack propagation. One difficulty is that the energy density with plastic unloading in the J-integral cannot be defined unambiguously. In this paper, we alternatively start from the analysis on the power balance, and propose a surface-forming energy release rate (ERR), which represents the energy available for separating the crack surfaces during the crack propagation and excludes the loading-mode-dependent plastic dissipation. Therefore the surface-forming ERR based fracture criterion has wider applicability, including elastic–plastic crack propagation problems. Several formulae are derived for calculating the surface-forming ERR. From the most concise formula, it is interesting to note that the surface-forming ERR can be computed using only the stress and deformation of the current moment, and the definition of the energy density or work density is avoided. When an infinitesimal contour is chosen, the expression can be further simplified. For any fracture behaviors, the surface-forming ERR is proven to be path-independent, and the path-independence of its constituent term, so-called J_s-integral, is also investigated. The physical meanings and applicability of the proposed surface-forming ERR, traditional ERR, J_s-integral and J-integral are compared and discussed. Besides, we give an interpretation of Rice paradox by comparing the cohesive fracture model and the surface-forming ERR based fracture criterion.     

Energy methods for fracture-toughness determination in concrete

The concepts for evaluating the energetic fracture parameter,G _ic, experimentally for concrete have been disputed because of the inconsistency and diversity of the results obtained. In order to re-examine the validity of this fracture parameter, experimental investigations were conducted to compare various methods for determiningG _ic with considerations of effects of microcracking and slow crack growth prior to the onset of instability leading to failure. These studies were made using small beams in three-point bending and utilizing compliance calibration techniques to estimate crack lengths. The following approaches were used to determine the critical energy-release rate: (1) stress-intensity-factor method which uses linear-elastic fracture mechanics to relateG _ic toK _ic; (2) theJ-integral concept in which energy change per unit crack extension was measured experimentally; and (3) Petersson's Method (modified) in which the total energy absorbed to instability is divided by the uncracked cross-sectional area to obtain the energy per unit area during the fracture process. All of these approaches resulted in excellent agreement and also good consistency. Invariancy witha/W was obtained when extended crack lengths were considered. Paper was presented at V International Congress on Experimental Mechanics held in Montreal, Quebec, Canada on June 10–15, 1984.     

Modeling of crack growth in ductile solids: a three-dimensional analysis

A population of several spherical voids is included in a three-dimensional, small scale yielding model. Two distinct void growth mechanisms, put forth by [Int. J. Solids Struct. 39 (2002) 3581] for the case of a two-dimensional model containing cylindrical voids, are well contained in the model developed in this study for spherical voids. A material failure criterion, based on the occurrence of void coalescence in the unit cell model, is established. The critical ligament reduction ratio, which varies with stress triaxiality and initial porosity, is used to determine ligament failure between the crack tip and the nearest void. A comparison of crack initiation toughness of the model containing cylindrical voids with the model containing spherical voids reveals that the material having a sizeable fraction of spherical voids is tougher than the material having cylindrical voids. The proposed material failure determination method is then used to establish the fracture resistance curve (J–R curve) of the material. For a ductile material containing a small volume fraction of microscopic voids initially, the void by void growth mechanism prevails, which results in a J–R curve having steep slope. On the other hand, for a ductile material containing a large volume fraction of initial voids, the multiple voids interaction mechanism prevails, which results in a flat J–R curve. Next, the effect of T-stress on fracture resistance is examined. Finally, nucleation and growth of secondary microvoids and their effects on void coalescence are briefly discussed.     

Effect of specimen size and crack depth on 3D crack-front constraint for SENB specimens

Three-dimensional (3D) elastic–plastic finite element analyses (FEA) are performed to study constraint effect on the crack-front stress fields for single-edge notched bend (SENB) specimens. Both rectangular and square cross-section of the specimens with a deep crack of a/W=0.5 are considered to investigate the effect of specimen size. A square-cross-section specimen with a shallow crack of a/W=0.15 is also considered to examine the effect of crack depth. Stresses from FEA at the crack front on different planes of the specimen are compared with those determined by the J–A₂ three-term solution. Results show that in-plane stress fields can be characterized by the three-term solution throughout the thickness even in the region near the free surface. Cleavage fracture toughness data is compared to predict the effects of specimen size and crack depth on fracture behavior. It is found that the distributions of crack opening stress are nearly the same for the SENB specimens at the critical J which is consistent with the RKR model. Furthermore our results indicate that there is a distinct relationship between the crack-front constraint and the cleavage fracture toughness. By introducing the failure curves, the minimum fracture toughness and scatter band can be well captured using the J–A₂ approach.     

Yielding on inclined planes at the tip of a crack loaded in uniform tension

Plane-strain yielding from a crack in an infinite elastic body is represented here by a distribution of edge dislocations on two planes inclined at angles ±ga to the crack plane, and the equilibrium condition is solved numerically. Approximate analytical expressions are obtained for the plastic-zone length, the crack opening displacement, and the J-integral, as functions of the applied stress and α. A comparison with a co-planar model of the plastic zone gives very similar results for α ≈ 65°. It is shown that fracture criteria based either on a critical crack opening displacement (COD) or on a critical value of J are always different, and the use of the former may lead to critical defect-sizes which are twice as large as those given by the latter. Furthermore, COD appears not to be a well-defined material property. The critical J criterion gives a fracture stress which is α-dependent : this may be responsible for deviations towards results of linear elastic fracture mechanics when post-yield fracture mechanics is used to analyse extensive yielding. The changes in the stress field of the crack due to the existence of the plastic zone are also discussed.     

K Ic and J-resistance-curve measurements on Nevada tuff

Linear-elastic fracture mechanics and J-integral test methods were used to determine the static fracture behavior of a porous volcanic tuff. Notched and prefatigued specimens of two different sizes were tested in three-point bending. The fracture toughness (K _Ic ) and theJ-resistance (J vs. crack growth) curve for each specimen was determined. The results indicate that there is good agreement between the fracture parameters determined by the two methods; however, there is some dependence of the results on specimen size. Possible reasons for this specimen-size dependence are discussed.     

J-integral and crack driving force in elastic–plastic materials

This paper discusses the crack driving force in elastic–plastic materials, with particular emphasis on incremental plasticity. Using the configurational forces approach we identify a “plasticity influence term” that describes crack tip shielding or anti-shielding due to plastic deformation in the body. Standard constitutive models for finite strain as well as small strain incremental plasticity are used to obtain explicit expressions for the plasticity influence term in a two-dimensional setting. The total dissipation in the body is related to the near-tip and far-field J-integrals and the plasticity influence term. In the special case of deformation plasticity the plasticity influence term vanishes identically whereas for rigid plasticity and elastic-ideal plasticity the crack driving force vanishes. For steady state crack growth in incremental elastic–plastic materials, the plasticity influence term is equal to the negative of the plastic work per unit crack extension and the total dissipation in the body due to crack propagation and plastic deformation is determined by the far-field J-integral. For non-steady state crack growth, the plasticity influence term can be evaluated by post-processing after a conventional finite element stress analysis. Theory and computations are applied to a stationary crack in a C(T)-specimen to examine the effects of contained, uncontained and general yielding. A novel method is proposed for evaluating J-integrals under incremental plasticity conditions through the configurational body force. The incremental plasticity near-tip and far-field J-integrals are compared to conventional deformational plasticity and experimental J-integrals.     

A generalized force measure of conditions at crack tips

A tentative measure of the forces tending to cause crack growth—the apparent crack extension force—is defined within the framework of continuum mechanics. By an associated fracture criterion initiation of growth may be predicted as well as the direction of preferred growth. The theory is specialized to elastic, viscoelastic and elastic-plastic materials. Under specified conditions the apparent crack extension force may be expressed by surface integrals over the boundary of an arbitrary part of the body for quasi-static deformation and for steady-state propagation of the crack. For plane deformation and for infinitesimal deformation under plane stress conditions these surface integrals reduce to path independent line integrals which include the J integral by Rice[1] and the G integral by Sih[2] as special cases.     

Electro-mechanical crack systems: bound theorems,dual finite elements and error estimation

Because of the complexity of piezoelectric crack problems, it is hard to obtain closed-form solutions, and numerical methods are largely resorted to. Hence, the upper/lower bound estimation of piezoelectric fracture parameters is of theoretical and practical importance. in this paper, the path-independent integral I, which is the dual of the J-integral, for electro-mechanical coupling crack systems, is presented. The related bound theorems are established for J and I. Piezoelectric dual finite elements are presented for the numerical implementation of the bound analysis. Moreover, an error estimator is presented for the assessment of numerical accuracy of the piezoelectric fracture parameters.     

Global and local approaches to fracture normal to interfaces

The problem of a crack perpendicularly approaching a bimaterial interface is examined using both global and localapproaches to fracture. The global approach is based on the J-integral with a second parameter, Q, which scales the stress triaxiality ahead of the crack. The local approach is based on either brittle fracture(Beremin model ) or ductile fracture (Rice and Tracey model ). In the first case, the Weibull stress over the plasticzone is calculated. In the second case, the void growth rate is calculated at the tip of the crack over a representativevolume (generally associated with a characteristic length of the material ). After a brief summary of each approach,the results for a crack near an elastically homogeneous, plastically mismatched interface are presented. Thebehaviour of the bimaterial is expressed in relation to the behavior of the homogeneous material. It is shown thatthere is an effect on the crack behavior which depends on the direction of crack propagation, i.e. from the hardermaterial to the softer material or vice versa. This effect is examined as a function of change in yield strength ratioand hardening exponent, n. For the case of brittle fracture, the effect of changing the Weibull modulus, m, is also examined. The models based on the local approach show that both stress- and strain-controlledfracture mechanisms must be accounted for. This implies the necessity of using the two parameters J and Q in the global approach. This is due to the fact that the stress–strain fields ahead of the crack tip areaffected by the nature of the second material.     

Analyses of steady quasistatic crack growth in plane strain tension in elastic-plastic materials with non-isotropic hardening

Steady state crack propagation problems of elastic-plastic materials in Mode I, plane strain under small scale yielding conditions were investigated with the aid of the finite element method. The elastic-perfectly plastic solution shows that elastic unloading wedges subtended by the crack tip in the plastic wake region do exist and that the stress state around the crack tip is similar to the modified Prandtl fan solution. To demonstrate the effects of a vertex on the yield surface, the small strain version of a phenomenological J₂, corner theory of plasticity (Christoffersen, J. and Hutchinson, J. W. J. Mech. Phys. Solids,27, 465 C 1979) with a power law stress strain relation was used to govern the strain hardening of the material. The results are compared with the conventional J₂ incremental plasticity solution. To take account of Bauschinger like effects caused by the stress history near the crack tip, a simple kinematic hardening rule with a bilinear stress strain relation was also studied. The results are again compared with the smooth yield surface isotropic hardening solution for the same stress strain curve. There appears to be more potential for steady state crack growth in the conventional J₂ incremental plasticity material than in the other two plasticity laws considered here if a crack opening displacement fracture criterion is used. However, a fracture criterion dependent on both stress and strain could lead to a contrary prediction.     

In-Situ SEM J-Testing and Crack-Tip Micromechanics of Zircaloy-4 by Three-Point Bending

A three-point bend fixture has been designed, fabricated, and utilized to demonstrate the feasibility of performing in-situ J-testing at ambient and elevated temperatures inside a scanning electron microscope (SEM). Using the three-point bend test technique, in-situ SEM J-testing has been performed to measure the crack mouth opening displacement and crack extension as a function of the applied load in order to generate J-R curves for Zircaloy-4 at 25°C and 316°C. Once the J-R curve is determined, an equivalent K_J-resistance (K_J-R) curve is computed on the basis of a relationship between the J-integral (J) and the stress intensity factor (K). The J-R and K_J-R curves of Zircaloy-4 exhibit a rising R-curve behavior, while the elastic K-R curve underestimates the fracture resistance of Zircaloy-4 once substantial crack extension has occurred. For the specimen dimensions considered, the J-R curves generated by in-situ SEM J-tests are not sensitive to the specimen geometry and measure the actual fracture resistance of the material. Furthermore, the onset of crack extension is dictated by the emission of one or more slipbands from the crack tip, and a change in the crack-tip displacement field, followed by void formation along the slipband, and linkage of the voids with the main crack.     

三维计算断裂力学

断裂力学理论从1921年Griffith研究玻璃的脆性断裂问题开始,经历了从线弹性体系到弹塑性以及蠕变理论体系、从单参数到多参数体系和从理想的二维平面理论到实际三维含裂纹结构的三维断裂理论的发展历程。针对应力强度因子K和J积分以及C(t)积分的计算方法从理想化模型的理论计算发展到实际复杂工程结构裂纹体计算的各种商业软件平台以及专业的断裂理论分析平台。尤其是随着计算机技术的发展,对三维含裂结构的静态和扩展裂纹的计算模拟已经能够融入计算机辅助设计。结合本研究组近30年来在三维疲劳断裂理论和应用研究方面的体会,简述了三维计算断裂力学从裂纹体应力应变分析和断裂参数计算到三维蠕变断裂和疲劳裂纹扩展模拟的国内外进展,并对涉及的计算方法,包括原子尺度和跨尺度的计算模拟,以及目前面临的挑战性问题作了简要介绍和分析。     

Multiple specimen J-integral testing at intermediate rates

The effect of experimental technique on the measurement of elastic-plastic fracture toughness has been investigated. Multiple-specimenJ-integral measurements have been performed at an intermediate loading rate. One set of measurements was conducted completely under computer/servohydraulic control, while another set was performed with a fixture that mechanically limits the maximum displacement and load. For the single condition tested, a minimal effect was measured on the initiation fracture toughness. An increased effect was seen on the measured crack-growth-resistance (J-R) curve as the amount of crack growth increased. The consequences of such testing methods in terms of general (intermediate loading rate) testing will be discussed.