Finite-deformation analysis of the crack-tip fields under cyclic loading

The plane-strain crack subjected to mode I cyclic loading under small scale yielding was analysed. The influence of the load range, load ratio and overload on the near-tip deformation-, stress- and strain-fields was studied. Although the near-tip zones of appreciable cyclic plastic flow for all loading regimes matched closely one another, when scaled with (ΔK/σ_Y)², the activities of plastic flow within them manifested dependence on K_max and K_min, as well as on overload. Cyclic trajectories of the crack-tip opening displacement (CTOD) converged to stable self-similar loops of the sizes proportional to ΔK² and positions in CTOD-K plane dependent on the maximum K along the whole loading route, including an overload. Computed near-tip deformation evidenced plastic crack advance, this way visualising of the Laird–Smith concept of fatigue cracking. This crack growth by blunting-resharpening accelerated with rising ΔK and was halted by an overload. Crack closure upon unloading had no place. The affinities were revealed between computed near-tip stress–strain variables and the experimental trends of the fatigue crack growth rate, such as its dependence on K_max and K_min (or ΔK and K_max), and retardation by overload. Thus, the effects of loading parameters on fatigue cracking, hitherto associated with crack closure, are attributable to the stress–strain fields in front of it as the direct drives of the key fatigue constituents – damage accumulation and bond breaking.     

Finite deformation analysis of crack-tip opening in elastic-plastic materials and implications for fracture

Analyses of the stress and strain fields around smoothly-blunting crack tips in both non-hardening and hardening elastic-plastic materials, under contained plane-strain yielding and subject to mode I opening loads, have been carried out by use of a finite element method suitably formulated to admit large geometry changes. The results include the crack-tip shape and near-tip deformation field, and the crack-tip opening displacement has been related to a parameter of the applied load, the J-integral. The hydrostatic stresses near the crack tip are limited due to the lack of constraint on the blunted tip, limiting achievable stress levels except in a very small region around the crack tip in power-law hardening materials. The J-integral is found to be path-independent except very close to the crack tip in the region affected by the blunted tip. Models for fracture are discussed in the light of these results including one based on the growth of voids. The rate of void-growth near the tip in hardening materials seems to be little different from the rate in non-hardening ones when measured in terms of crack-tip opening displacement, which leads to a prediction of higher toughness in hardening materials. It is suggested that improvement of this model would follow from better understanding of void-void and void-crack coalescence and void nucleation, and some criteria and models for these effects are discussed. The implications of the finite element results for fracture criteria based on critical stress or strain, or both, is discussed with respect to transition of fracture mode and the angle of initial crack-growth. Localization of flow is discussed as a possible fracture model and as a model for void-crack coalescence.     

A reexamination of plasticity-induced crack closure in fatigue crack propagation

The crack closure concept is often used to consider the R-ratio and overload effects on fatigue crack growth. The presumption is that when the crack is closed, the external load produces negligible fatigue damage in the cracked component. The current investigation provides a reassessment of the frequently used concept with an emphasis on the plasticity-induced crack closure. A center cracked specimen made of 1070 steel was investigated. The specimen was subjected to plane-stress mode I loading. An elastic–plastic stress analysis was conducted for the cracked specimens using the finite element method. By applying the commonly used one-node-per-cycle debonding scheme for the crack closure simulations, it was shown that the predicted crack opening load did not stabilize when the extended crack was less than four times of the plastic zone size. The predicted opening load was strongly influenced by the plasticity model used. When the elastic–perfectly plastic (EPP) stress–strain relationship was used together with the kinematic hardening plasticity theory, the predicted crack opening load was found to be critically dependent on the element size of the finite element mesh model. For R = 0, the predicted crack opening load was greatly reduced when the finite element size became very fine. The kinematic hardening rule with the bilinear (BL) stress–strain relationship predicted crack closure with less dependence on the element size. When a recently developed cyclic plasticity model was used, the element size effect on the predicted crack opening level was insignificant. While crack closure may occur, it was demonstrated that cyclic plasticity persisted in the material near the crack tip. The cyclic plasticity was reduced but not negligible when the crack was closed. The traditional approaches may have overestimated the effect of crack closure in fatigue crack growth predictions.     

随机超载下疲劳裂纹扩展寿命的计算

针对在基本循环载荷上加入随机超载序列的疲劳裂纹扩展问题，应用裂纹闭合的概念考虑超载的迟滞效应，将延迟时间描述成纯离散型马尔可夫过程。对相应的柯尔莫哥洛夫－费勒微积分方程进行了分步求解，结合疲劳断裂分别出现在基本循环峰载作用时和超载作用时两种情况，计算出不同可靠度下裂纹的扩展寿命，并研究了超载大小和发生强度对扩展寿命的影响。     

Reverse shear on inclined planes at the tip of a sharp crack

Plane strain plastic yielding at a crack tip has been represented by edge dislocations with Burgers vectors parallel to symmetrical planes inclined at 70° and 45° to the plane of the crack. The plastic displacement and the stresses near the crack tip were calculated by a numerical method and the effect of a reduction in applied stress was determined. Removal of the whole or a part of the initial load produces reverse shear in regions of the slip band nearest the crack tip. The amount of reverse shear depends only on the reduction in the load and not on its initial value. The reverse shear is associated with the presence of negative dislocations and the stresses near the crack tip may become compressive even though the applied (remote) stress is still tensile. The degree and extent of compression depends on the reduction in applied stress and on its original value. It is argued that the residual compressive stresses produced under fluctuating loads may produce crack closure and crack arrest. The effect of residual plasticity in a slip band left behind a growing crack has been estimated. It is shown that after an overload the excess residual plasticity opposing crack opening rises to a maximum value when the crack tip has advanced some distance from the point where the overload was applied.     

Some possible alternative solutions of near-tip fields for steady crack growth in elastic perfectly-plastic media

Some possible alternative solutions of near-tip fields are studied for plane-strain Mode—I quasi-static steady crack growth in incompressible (ν=1/2) elastic perfectly-plastic media. A group of four-sector solutions and a three-sector solution in which the elastic-unloading region vanishes are given. Stress functions, plastic flow factors and plastic strains in each region are also given. Project supported by the State Education Commission under a funding program for Excellent University Young Faculties and National Natural Science Foundation of China.     

Calculation of near-tip non-singular stresses for pressurized cracks

When the crack surfaces are traction-free, there is only one constant term T in the near-tip stress field, which contributes uniformly to the stress component acting in the direction parallel to the crack flank. As to pressurized cracks, the non-singular part of the asymptotic stresses appears to be more complicated and is no longer characterized only by the constant T. In this work, an effective numerical approach is developed for calculation of the non-singular parts of the asymptotic near-tip stresses under the action of nonuniform crack surface pressures. With this approach, the near-tip non-singular stress field can be accurately evaluated by direct use of regular numerical methods such as finite elements.     

ASYMPTOTIC ANALYSIS OF PLANE-STRAIN MODE I STEADY-STATE CRACK GROWTH IN TRANSFORMATION TOUGHENING CERAMICS(Ⅱ)

Based on a constitutive law which includes the shear components of transformation plasticity, the asymptotic solutions to near-tip fields of plane-strain mode I steadity propagating cracks in transformed ceramics are obtained for the case of linear isotropic hardening. The stress singularity, the distributions of stresses and velocities at the crack tip are determined for various material parameters. The factors influencing the near-tip fields are discussed in detail.Project supported by the National Natural Science Foundation of China     

Prediction of fatigue crack growth life under variable-amplitude loading using finite element analysis

This article presents an elastic-plastic study aiming at predicting the fatigue crack growth (FCG) of 2024-T3 aluminum alloys under variable-amplitude loading. The proposed analysis needs the estimation of the residual stress distribution ahead of the crack tip during propagation. An elastic-plastic FE analysis has been implemented for modeling FCG using Chaboche's model. The FE study has been carried out through consideration of the loading history effect using the memory rules. Three different loading spectra have been applied in this work. The obtained results have been compared to the experimental ones and it has been proved that the suggested model has a better prediction of the FCG lives of cracked 2024-T3 aluminum alloy structures subjected to variable-amplitude loading.     

Numerical study of contact forces for crack closure analysis

Plasticity induced crack closure (PICC) has been widely studied using numerical models. Different numerical parameters can be considered to quantify the opening level, namely one based on the analysis of contact stresses at minimum load. A modified version of this parameter is proposed here, based on nodal contact forces instead of contact stresses. The predictions were found to be similar to those obtained from the contact status of 2nd node behind crack tip. The PICC_contact parameter was also found to be very consistent and adequate for parametric studies of the influence of different physical parameters. The contributions to the opening stress intensity factor of different points along crack flank were found to strongly decrease with distance to crack tip. The cumulative K_open between the crack tip and a distance of 0.1 mm was found to vary from 30% to 100%, increasing with stress ratio, R. Finally, a K solution was developed for punctual forces applied on crack flank and compared with a literature solution for infinite plates. A good agreement was found for plane strain state but significant differences of about 10% were found for plane stress state.     

A numerical study of crack tip constraint in ductile single crystals

In this work, the effect of crack tip constraint on near-tip stress and deformation fields in a ductile FCC single crystal is studied under mode I, plane strain conditions. To this end, modified boundary layer simulations within crystal plasticity framework are performed, neglecting elastic anisotropy. The first and second terms of the isotropic elastic crack tip field, which are governed by the stress intensity factor K and T-stress, are prescribed as remote boundary conditions and solutions pertaining to different levels of T-stress are generated. It is found that the near-tip deformation field, especially, the development of kink or slip shear bands, is sensitive to the constraint level. The stress distribution and the size and shape of the plastic zone near the crack tip are also strongly influenced by the level of T-stress, with progressive loss of crack tip constraint occurring as T-stress becomes more negative. A family of near-tip fields is obtained which are characterized by two terms (such as K and T or J and a constraint parameter Q) as in isotropic plastic solids.     

基于裂纹闭合效应的高载迟滞预测模型修正

裂纹闭合效应通常是导致I 型裂纹扩展在高载作用下发生迟滞效应的主要因素之一．本文采用汽车薄板QSTE340TM 材料，针对不同应力比，高载比条件下疲劳裂纹扩展行为进行了实验研究．论文通过断面分析，针对各参数对裂纹闭合效应的具体影响进行了分析讨论，认为裂纹作用区域随裂纹扩展而动态变化，从而提出了一种对有效应力强度因子幅的修正方法．通过在原有模型中引入幂函数形式的动态变量α，表征裂纹闭合效应的作用比例随裂纹长度的动态变化，取得了较好的效果．     

Creep deformation at crack tips in elastic-viscoplastic solids

The evaluation of crack growth tests under creep conditions must be based on the stress analysis of a cracked body taking into account elastic, plastic and creep deformation. In addition to the well-known analysis of a cracked body creeping in secondary (steady-state) creep, the stress field at the tip of a stationary crack is calculated for primary (strain-hardening) or tertiary (strain-softening) creep of the whole specimen. For the special hardening creep-law considered, a path-independent integral C¹_h, can be defined which correlates the near-tip field to the applied load.It is also shown how, after sudden load application, creep strains develop in the initially elastic or, for a higher load level, plastic body. Characteristic times are derived to distinguish between short times when the creep-zones, in which creep strains are concentrated, are still small, and long times when the whole specimen creeps extensively in primary and finally in secondary and tertiary creep. Comparing the creep-zone sizes with the specimen dimensions or comparing the characteristic times with the test duration, one can decide which deformation mechanism prevails in the bulk of the specimen and which load parameter enters into the near-tip stress field and determines crack growth behavior. The governing load parameter is the stress intensity factor K₁ if the bulk of the specimen is predominantly elastic and it is the J-integral in a fully-plastic situation when large creep strains are still confined to a small zone. The C¹_h-integral applies if the bulk of the specimen deforms in primary or tertiary creep, and C¹ is the relevant load parameter for predominantly secondary creep of the whole specimen.     

Dynamic steady crack growth in elastic-plastic solids—propagation of strong discontinuities

The near-tip field of a mode I crack growing steadily under plane strain conditions is studied. A key issue is whether strong discontinuities can propagate under dynamic conditions. Theories which impose rather restrictive assumptions on the structure of an admissible deformation path through a dynamically propagating discontinuity have been proposed recently. Asymptotic solutions for dynamic crack growth, based on such theories, do not contain any discontinuities. In the present work a broader family of deformation paths is considered and we show that a discontinuity can propagate dynamically without violating any of the mechanical constitutive relations of the material. The proposed theory for the propagation of strong discontinuities is corroborated by very detailed finite element calculations. The latter shows a plane of strong discontinuity emanating from the crack tip (with its normal pointing in the direction of crack advance) and moving with the tip. Elastic unloading ahead of and/or behind the plane of discontinuity and behind the crack tip have also been observed.The numerical investigation is performed within the framework of a boundary layer formulation whereby the remote loading is fully specified by the first two terms in the asymptotic solution of the elasto-dynamic crack tip field, characterized by K₁, and T. It is shown that the family of near-tip fields, associated with a given crack speed, can be arranged into a one-parameter field based on a characteristic length, L_g, which scales with the smallest dimension of the plastic zone. This extends a previous result for quasi-static crack growth.     

Finite element analysis for steady-state hydride-induced fracture in metals by composite model

Delayed hydride cracking, which is observed in hydride-forming metals, due to the precipitation of hydrides near the crack tip, is investigated under conditions of constant temperature and crack velocity, plane strain and small-scale hydride-precipitation. The coupling of the operating physical processes of hydrogen-diffusion, hydride precipitation and material deformation is taken into account. The material is assumed to be an elastic composite made of hydrides and solid solution, with properties depending locally on the volume fraction of the hydrides. In the present analysis, the composite elastic properties have been derived by a generalized self consistent model for particulate composites. With respect to hydride-precipitation, two cases have been considered: (i) precipitation in a homogeneous medium with elastic properties, equal to the effective properties of the composite and (ii) precipitation in an inhomogeneous medium, where the expanding hydride has different elastic properties than those of the surrounding solid solution. The differences between the near-tip field distributions, produced by the two precipitation models, are relatively small. The effect of the hydrogen concentration far from the crack tip, on the near-tip field is also studied. It is shown that for small crack growth velocities, near the threshold stress intensity factor, the remote hydrogen concentration weakly affects the normalized stress distribution in the hydride-precipitation zone, which is controlled by the thermodynamically required hydrostatic stress, under hydrogen chemical equilibrium. However, for values of the applied stress intensity factor and the crack tip velocity, away from the threshold stress intensity factor and crack arrest, the effect of remote hydrogen concentration on the normalized near-tip stress field is strong. Reduction of the remote hydrogen concentration generally leads to reduction of the hydride-precipitation zone and increase of the near-tip stresses. Also reduction of the remote hydrogen concentration leads to distributions closer to those under hydrogen chemical equilibrium.     

A model for predicting the retardation effect following a single overload

Many engineering components are subjected to variable amplitude loading history. It is well known that retardation in fatigue crack growth occurs due to application of single overloads in a constant amplitude loading block. Many models have been proposed to capture this counter intuitive phenomenon which has resulted in improved understanding of retardation effect following tensile overloads and consequently resulting in better life prediction models. The proposed study is focused on to evaluation of retardation in fatigue life due to application of a single overload. A model for prediction of crack growth and crack growth rate following single overloads is presented. Several modifications to Wheeler’s growth idea are proposed, which incorporate a consideration for effective stress intensity factor, based on Elber’s concept of crack closure, relationship between overload ratio and the Wheeler’s exponent, and fatigue growth rate calculations. The results presented here show that plastic zone interaction following overload and the consideration of crack closure explain retardation effect following a single overload. Correlation between analysis and experimental data obtained from several sources in literature show that the scheme, is robust and provides an insight into the nonlinear aspect of crack growth results. The model has been tested for 2024-T3 aluminum alloy and 6061-T6 aluminum alloy and thorough calibrations performed, established the fidelity of the program.     

Modeling quasi-static crack growth with the extended finite element method Part II: Numerical applications

In Part I [Int. J. Solids Struct., 2003], we described the implementation of the extended finite element method (X-FEM) within Dynaflow™, a standard finite element package. In our implementation, we focused on two-dimensional crack modeling in linear elasticity. For crack modeling in the X-FEM, a discontinuous function and the near-tip asymptotic functions are added to the finite element approximation using the framework of partition of unity. This permits the crack to be represented without explicitly meshing the crack surfaces and crack propagation simulations can be carried out without the need for any remeshing. In this paper, we present numerical solutions for the stress intensity factor for crack problems, and also conduct crack growth simulations with the X-FEM. Numerical examples are presented with a two-fold objective: first to show the efficacy of the X-FEM implementation in Dynaflow™; and second to demonstrate the accuracy and versatility of the method to solve challenging problems in computational failure mechanics.     

Plane-strain crack problems in microstructured solids governed by dipolar gradient elasticity

The present study aims at determining the elastic stress and displacement fields around the tips of a finite-length crack in a microstructured solid under remotely applied plane-strain loading (mode I and II cases). The material microstructure is modeled through the Toupin-Mindlin generalized continuum theory of dipolar gradient elasticity. According to this theory, the strain-energy density assumes the form of a positive-definite function of the strain tensor (as in classical elasticity) and the gradient of the strain tensor (additional term). A simple but yet rigorous version of the theory is employed here by considering an isotropic linear expression of the elastic strain-energy density that involves only three material constants (the two Lamé constants and the so-called gradient coefficient). First, a near-tip asymptotic solution is obtained by the Knein-Williams technique. Then, we attack the complete boundary value problem in an effort to obtain a full-field solution. Hypersingular integral equations with a cubic singularity are formulated with the aid of the Fourier transform. These equations are solved by analytical considerations on Hadamard finite-part integrals and a numerical treatment. The results show significant departure from the predictions of standard fracture mechanics. In view of these results, it seems that the classical theory of elasticity is inadequate to analyze crack problems in microstructured materials. Indeed, the present results indicate that the stress distribution ahead of the crack tip exhibits a local maximum that is bounded. Therefore, this maximum value may serve as a measure of the critical stress level at which further advancement of the crack may occur. Also, in the vicinity of the crack tip, the crack-face displacement closes more smoothly as compared to the standard result and the strain field is bounded. Finally, the J-integral (energy release rate) in gradient elasticity was evaluated. A decrease of its value is noticed in comparison with the classical theory. This shows that the gradient theory predicts a strengthening effect since a reduction of crack driving force takes place as the material microstructure becomes more pronounced.     

Crack-tip fields for fast fracture of an elastic-plastic material

An asymptotic analysis of the near-tip fields is given for transient crack propagation in an elastic-plastic material. The material is characterized by J₂ flow theory together with a bilinear effective stress-strain curve. Both plane stress and plane strain conditions have been considered. Explicit results are given for the order of the crack-tip singularity, the angular position at which unloading occurs, and the angular variation of the near-tip stresses, all as functions of the crack-tip speed and the ratio of the slopes of the two portions of the bilinear stress-strain relation. It was found that the results are much more sensitive to the elastic-plastic constitutive relation than to the crack speed. This result is important for numerical analyses of dynamic crack propagation problems.     

平面应变下紧凑拉伸试样的动态断裂韧性的实验研究

材料的动态断裂韧性是衡量材料在动载荷作用下抵杭裂纹扩展能力的重要指标,以往的材料动态断裂韧性测试多采用三点弯曲试样,而针对紧凑拉伸试样的动态断裂韧性研究很少.本文将紧凑拉伸试样(即CT试样)简化成等效弹簧质量模型,得到了CT试样动态应力强度因子的近似表达式.对Hopkinson压杆装置进行了改进,利用改进后的实验装置进...