Mode i and mixed mode fields with incomplete crack tip plasticity

Plane strain slip line fields, in which plasticity does not fully surround the crack tip have been developed for mode I and mixed mode I\II cracks under contained yielding. Analytical solutions have been assembled using slip line theory for the plastic sectors and semi-infinite wedge solutions for the elastic sectors. These solutions are compared with finite element solutions based on modified boundary layer formulations. The analytical solutions agree well with numerical solutions, and form a family of fields with incomplete plasticity around the crack tip.     

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.     

Numerical plane stress elastic–perfectly plastic crack analysis under Tresca yield condition with comparison to Dugdale plastic strip model

A number of plane stress numerical analyses of the mode I elastoplastic fracture mechanics problem have been performed in the past using the Huber–Mises yield criterion. This study employs instead the Tresca yield condition using an incremental theory of plasticity for a stationary crack. A commercial finite element program is used to solve the opening mode of fracture problem (mode I) for a square plate containing a central crack under generalized plane stress loading conditions. A biaxial uniform tensile traction is applied to the edges of a thin plate composed of a linear elastic non-work hardening material under small strain assumptions. The finite element results are compared with the analytical predictions of the Dugdale plastic strip model for a crack in an infinite plate subject to a biaxial uniform load at infinity.     

考虑材料循环塑性的疲劳裂纹扩展模拟

提出了一种考虑材料循环塑性性能的研究疲劳裂纹扩展与闭合行为的有限元模拟方法．对所选用的循环塑性本构关系进行了基本实验检验．探讨了在疲劳裂纹扩展有限元分析中网格尺寸的影响，给出了网格优化准则．研究了在循环硬化条件下考虑裂纹合效应时裂纹面张开廓形、裂纹尖端应力、应变场和正反向塑性区的演变规律．对于循环硬化和不同循环应力比Ｒ等因素对裂纹张开应力水平的影响也作了考察     

Finite deformation analysis of mechanism-based strain gradient plasticity: torsion and crack tip field

A finite deformation theory of mechanism-based strain gradient (MSG) plasticity is developed in this paper based on the Taylor dislocation model. The theory ensures the proper decomposition of deformation in order to exclude the volumetric deformation from the strain gradient tensor since the latter represents the density of geometrically necessary dislocations. The solution for a thin cylinder under large torsion is obtained. The numerical method is used to investigate the finite deformation crack tip field in MSG plasticity. It is established that the stress level around a crack tip in MSG plasticity is significantly higher than its counterpart (i.e. HRR field) in classical plasticity.     

Evaluation of crack initiation angle under mixed mode loading at diverse strain rates

Crack initiation angle, under mixed mode loading at several strain rates, is analysed using an experimental–numerical approach. The physical phenomenon for the problem at hand is influenced by the local and global conditions. One of such factors is the strain rate at the crack tip. For this purpose, PMMA plates with centred angled cracks under mixed mode loading were tested. The strain rate at the neighbourhood of the crack tip before crack propagation was evaluated. Considering that this material is strain rate sensitive, the numerical models were calibrated with the modulus of elasticity measured in tension tests at the observed strain rates. Numerical evaluations were performed with the finite element method in conjunction with the volume energy density criterion. An improvement in the evaluation of the crack propagation angle was observed. In order to complete the analysis, the crack initiation angle was also evaluated with the strain energy density factor S, considering the mechanical properties of PMMA, as evaluated at the observed strain rates, and the stress intensity factors k₁ and k₂. Results are in agreement with those observed experimentally.     

一种新的计算各向异性材料裂纹尖端应力强度因子杂交元法

利用有限元特征分析法研究了平面各向异性材料裂纹端部的奇性应力指数以及应力场和位移场的角分布函数,以此构造了一个新的裂纹尖端单元。文中利用该单元建立了研究裂纹尖端奇性场的杂交应力模型,并结合Hellinger-Reissner变分原理导出应力杂交元方程,建立了求解平面各向异性材料裂纹尖端问题的杂交元计算模型。与四节点单元相结合,由此提出了一种新的求解应力强度因子的杂交元法。最后给出了在平面应力和平面应变下求解裂纹尖端奇性场的算例。算例表明,本文所述方法不仅精度高,而且适应性强。     

Crack separation energy-rates for the DBCS model under biaxial modes of loading

A Modified version of the Dugdale-Bilby-Cottrell-Swinden (DBCS) model simulating the effect of plasticity at the tip of a crack in an infinite region was used by kfouri and rice (1978) to calculate the crack separation energy-rate G^Δ corresponding to a finite crack growth step Δa during plane strain mode I crack extension. The loading consisted of a remote uniaxial tension σp applied normally to the plane of the crack. Using Rice's path-independent integral J to characterize the applied load in the crack tip region, and assuming the length R of the crack tip plastic zone to be small compared with the length of the crack, an analytical expression was derived relating the ratios (J/G^Δ) and (2a/R) for small values of (2a/R). The analytical solution was incomplete in itself in that the value assumed in the plastic region of the DBCS model for the normal stress Y acting on the extending crack surfaces was the product of the yield stress in uniaxial tension σ_Y and an unknown parameter C, the value of which depends on the effect of the local hydrostatic stresses in the case of plane strain conditions. The analytical solution was compared with a numerical solution obtained from a plane strain elastic-plastic finite element analysis on a centre-cracked plate (CCP) of material obeying the von Mises yield criterion. The value used for the yield stress was 310 MN/m² and moderate isotropic linear hardening was applied with a tangent modulus of 4830 MN/m². A uniaxial tension σp was applied on the two appropriate sides of the plate. The comparisons showed that the analytical and finite element solutions were mutually consistent and they enabled the value of C to be established at 1.91. In the present paper similar comparisons are made between the analytical solution and the finite element solutions for the CCP of the same material under different biaxial modes of loading. By assuming further that the form of the analytical solution does not depend on the details of the geometry and of the loading at remote boundaries, a comparison has also been made with the results of a finite element analysis on a compact tension specimen (CTS) made of the same material as the CCP. The different values of C obtained in each case are consistent with investigations by other authors on the effect of load biaxiality on crack tip plasticity.     

The role of macroscopic hardening and individual length-scales on crack tip stress elevation from phenomenological strain gradient plasticity

This paper quantifies the effect of strain gradient plasticity (SGP) on crack tip stress elevation for a broad range of applied loading conditions and constitutive model parameters, including both macroscopic hardening parameters and individual material length-scales controlling gradient effects. Finite element simulations incorporating the Fleck-Hutchinson SGP theory are presented for an asymptotically sharp stationary crack. Results identify fundamental scaling relationships describing (i) the physical length-scales over which strain gradients are prominent, and (ii) the degree of stress elevation over conventional Hutchinson-Rice-Rosengren (HRR) fields. Results illustrate that the three length-scale theory predicts much larger SGP effects than the single length-scale theory. Critically, the first length-scale parameter dominates SGP stress elevation: this suggests that SGP effects in fracture can be predicted using the length-scales extracted from nanoindentation, which exhibits similar behavior. Transitional loading/material parameters are identified that establish regimes of SGP relevance: this provides the foundation for the rational application of SGP when developing new micromechanical models of crack tip damage mechanisms and associated subcritical crack propagation behavior in structural alloys.     

双材料界面裂纹应力强度因子的边界元分析

采用双材料基本解建立边界元法基本方程,计算双材料界面裂纹尖端附近的应用力和位移场。不离散界面,并设置面力奇异四分之一点裂尖单元以提高计算精度。数值结果表明,本文的方法具有较高的精度和效率。     

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.     

Interface crack between a compressible elastomer and a rigid substrate with finite slippage

We study the deformation of a crack between a soft elastomer and a rigid substrate with finite interfacial slippage. It is assumed that slippage occurs when the interfacial shear traction exceeds a threshold. This leads to a slip zone ahead of the crack tip where the shear traction is assumed to be equal to the constant threshold. We perform asymptotic analysis and determine closed-form solutions describing the near-tip crack opening displacement and the corresponding stress distributions. These solutions are consistent with numerical results based on finite element analysis. Our results reveal that slippage can significantly affect the deformation and stress fields near the tip of the interface crack. Specifically, depending on the direction of slippage, the crack opening profile may appear more blunted or sharpened than the parabola arising from for the case of zero interfacial shear traction or free slippage. The detailed crack opening profile is determined by the constant shear traction in the slip zone. More importantly, we find that the normal stress perpendicular to the interface can increase or decrease when slippage occurs, depending on the direction of slippage and the shear traction in the slip zone.     

裂尖大应变细观断裂研究

本文用反映空穴形核成长的Gurson本构方程来描述裂尖区域材料在大应变情形下的力学特性，并进一步考虑了空穴演变对材料杨氏模量的影响。文中用上述本构方程分别结合弹塑性大应变有限元方法对平面应变I型裂纹问题作了计算，分析了裂尖应力分布、裂尖形状变化和裂尖空穴演变过程，并与用Prandtl-Reuss本构方程教育处的结果作了比较。     

Crack tip field andJ-integral with strain gradient effect

The mode I plane strain crack tip field with strain gradient effects is presented in this paper based on a simplified strain gradient theory within the framework proposed by Acharya and Bassani. The theory retains the essential structure of the incremental version of the conventionalJ ₂ deformation theory. No higher-order stress is introduced and no extra boundary value conditions beyond the conventional ones are required. The strain gradient effects are considered in the constitutive relation only through the instantaneous tangent modulus. The strain gradient measures are included into the tangent modulus as internal parameters. Therefore the boundary value problem is the same as that in the conventional theory. Two typical crack problems are studied: (a) the crack tip field under the small scale yielding condition induced by a linear elastic mode-IK-field and (b) the complete field for a compact tension specimen. The calculated results clearly show that the stress level near the crack tip with strain gradient effects is considerable higher than that in the classical theory. The singularity of the strain field near the crack tip is nearly equal to the square-root singularity and the singularity of the stress field is slightly greater than it. Consequently, theJ-integral is no longer path independent and increases monotonically as the radius of the calculated circular contour decreases. The project supported by the National Natural Science Foundation of China (19704100 and 10202023) and the Natural Science Foundation of Chinese Academy of Sciences (KJ951-1-20)     

Analyses on micro-fracture of crack tip under large deformation conditions

In the present paper, Gurson's constitutive equation, which takes into account the development of voids, is used to study the behaviour of the material in the region near crack tip. Furthermore, the effect of void development on Young's modulus, which was not considered by Gurson, is taken into consideration. The analyses on void development, on stress distribution near crack tip, and on the variance of COD for the plane strain mode I problem are carried out with the large elastic-plastic deformation finite element method. The results are compared with those estimated from the Prandtl-Reuss constitutive equation.     

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.     

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.     

Mode I and mode II crack tip asymptotic fields with strain gradient effects

The strain gradient effect becomes significant when the size of fracture process zone around a crack tip is comparable to the intrinsic material lengthl, typically of the order of microns. Using the new strain gradient deformation theory given by Chen and Wang, the asymptotic fields near a crack tip in an elastic-plastic material with strain gradient effects are investigated. It is established that the dominant strain field is irrotational. For mode I plane stress crack tip asymptotic field, the stress asymptotic field and the couple stress asymptotic field can not exist simultaneously. In the stress dominated asymptotic field, the angular distributions of stresses are consistent with the classical plane stress HRR field; In the couple stress dominated asymptotic field, the angular distributions of couple stresses are consistent with that obtained by Huang et al. For mode II plane stress and plane strain crack tip asymptotic fields, only the stress-dominated asymptotic fields exist. The couple stress asymptotic field is less singular than the stress asymptotic fields. The stress asymptotic fields are the same as mode II plane stress and plane strain HRR fields, respectively. The increase in stresses is not observed in strain gradient plasticity for mode I and mode II, because the present theory is based only on the rotational gradient of deformation and the crack tip asymptotic fields are irrotational and dominated by the stretching gradient. The project supported by the National Natural Science Foundation of China (19704100), National Natural Science Foundation of Chinese Academy of Sciences (KJ951-1-20), CAS K.C. Wong Post-doctoral Research Award Fund and Post-doctoral Science Fund of China     

Asymptotic character of mixed mode in plane deformation of crack in rubber-like material

The asymptotic stress and strain distribution near a crack tip in rubber-like materials is determined by finite element for in-plane mixed mode loading. For large strain, the crack tip field is always in a state of uniaxial tension. The shear load affects only the orientation of the deformed near tip field in the space. A good agreement is obtained between the theoretical and numerical results.     

Stress and energy density fields near a blunted crack front in power hardening material

Using Jaumann and Dienes rates of Euler stress in elastic-plastic constitutive equations of finite deformation, plane strain finite element analysis for a compact tension specimen with a blunted crack front is made. The Euler stress, Kirchhoff stress and volume strain energy density near a blunted crack tip are computed. Constitutive relations with different deformation rates affect the the near crack tip solution in a region within an order of magnitude of the crack opening displacement. The results differed from the corresponding solution of deformation plasticity (or nonlinear elasticity) with increasing deformation. They are smaller in a local region of about 2 to 10 times of the crack opening distance.The volume energy density near the crack tip is computed, the stationary values of which determine the locations of extensive yielding and possible sites of crack initiation. It remained nearly constant with increasing deformation. Such a character tends to support the volume energy density criterion as a means for quantifying the ductile fracture behavior of metals.