Determination of critical stress triaxiality along yield locus of isotropic ductile materials under plane strain condition

It is widely accepted that failure due to plastic deformation in metals greatly depends on the stress triaxiality factor (TF). This article investigates the variation of stress triaxiality along the yield locus of ductile materials. Von Mises yield criteria and triaxiality factor have been used to determine the critical limits of stress triaxiality for the materials under plane strain condition. A generalized mathematical model for triaxiality factor has been formulated and a constrained optimization has been carried out using genetic algorithm. Finite element analysis of a two dimensional square plate has been carried out to verify the results obtained by the mathematical model. It is found that the set of values of the first and the second principal stresses on the yield locus, which results in maximum stress triaxiality, can be used to determine the location at which crack initiation may occur. Thus, the results indicate that while designing a certain component, such combination of stresses which leads the stress triaxiality to its critical value, should be avoided.     

焊接接头中的裂端应力三轴性

采用弹塑性有限元方法研究了平面应变条件下裂位于焊缝中心不同强度匹配和裂纹深度焊接接头的裂端应力三轴性。结果发现,与均匀材料相比,同等载荷水平下,高匹配接头裂端应力三轴性降低,低匹配接头裂端应力三轴性升高;随裂纹变浅,应力三轴性降低,强度匹配的影响也更为显著。从焊缝和母材塑性变形与应变硬化交互作用方面进行了解释。     

Comparison of cohesive zone parameters and crack tip stress states between two different specimen types

The cohesive zone parameters (separation energy and cohesive strength) and the crack tip triaxialities are compared between a compact tension (CT) and a double edge notched tension (DENT) specimen with smooth side-surfaces. The material is a pressure vessel steel 20MnMoNi55. The cohesive zone parameters are determined by fitting the simulated crack extensions near the midsection to the experimental data. The purpose of the study is to understand the relationship between the cohesive zone parameters and the crack tip stress triaxiality. The results show that for the same cohesive zone parameters the crack tip triaxiality near the midsection is lower in DENT specimens than in CT specimens. When the separation energy is set constant for CT and DENT specimens, the cohesive strength for the DENT specimens should be significantly lower than that for the CT specimens in order to make the simulated crack extensions near the midsection fit to the experimental data. Near the midsection, the cohesive strength and crack tip triaxiality influence each other: the specimen with a higher stress triaxiality has a higher cohesive strength; an increase of cohesive strength results in an increase of the crack tip triaxiality.     

Crack growth prediction and non-linear analysis for an elasto-plastic solid

In this paper, a variable radius for the plastic zone is introduced and a maximum principal stress criterion is proposed for the prediction of crack initiation and growth. It is assumed that the direction of crack initiation coincides with the direction of the maximum principal stress. The von Mises yield criterion is applied to define the plastic zone, instead of assuming a plastic zone with a constant distance r from the crack tip. An improvement is made to this fracture criterion, and the criterion is extended to study the crack growth characteristics of mixed mode cracks. Based on the concept of frictional stress intensity factor, k_f, the rate of fatigue crack propagation, db/dN, is postulated to be a function of the effective stress intensity factor range, Δk_eff. Subsequently, this concept is applied to predict crack growth due to fatigue loads. The proposed crack growth model is discussed by comparing the experimental results with those obtained using the maximum principal stress criterion.     

Investigation of fracture characterizations of functionally graded materials

In this paper, the mode I crack problem of functionally gradient materials (FGMs) with the gradient direction parallel to the crack is discussed, and the differences of stress distribution between the gradient materials and the homogeneous materials are analyzed. It is shown that a mode I crack problem of FGMs with the gradient direction parallel to the crack direction can become a mixed‐mode crack problem. In FGMs, the crack initiation angles are determined by the fracture toughness gradient, elastic modulus and crack mode. If the gradient coefficients are small, the crack initiation angles in FGMs are the same as those in homogeneous materials. If the elastic modulus gradient is large, the principal stress terms without the gradient coefficients can be ignored in obtaining the crack initiation angle. In this study, all the above results are generalized to the mixed‐mode crack problems with arbitrary angle between the gradient direction and the crack direction.     

混凝土结构I型裂纹裂尖塑性区研究

应用双剪统一强度理论,研究了I型裂纹的塑性变形问题。给出了包含反映材料拉压性能差异的参数拉压比及反映中间主应力效应的参数b的I型裂纹裂尖塑性区形状和大小的统一解。已有的Tresca准则、Mises准则和Mohr-Coulomb准则解均是本文的特例或线性逼近。针对混凝土结构,画出了不同参数b情况下的裂尖塑性区半径变化图。得出了材料拉压比对I型裂纹裂尖塑性区影响很大。b对I型裂纹裂尖塑性区影响随拉压比的不同而不同,拉压比较大时,b对塑性区影响大,拉压比较小时,b对塑性区影响小的结论。该统一解可以适应于各种不同材料,能充分发挥材料潜力,具有普遍性和广泛的适应性,有一定的工程应用价值。结论对于研究各种材料的断裂问题有参考作用。     

Transferability of specimen J–R curve to straight pipe with circumferential surface flaw

Specimen J–R curve is extensively used for structural integrity of large components. It is well known that J–R curve heavily depends on constraint level ahead of crack tip in remaining ligament. In earlier work, it was demonstrated that J–R curve from Three Point Bending (TPB) specimen is transferable to straight pipe with circumferential through wall crack. In this paper, the transferability of J–R curve is investigated from TPB specimen to pipe with circumferential surface crack. A 16 in. diameter pipe with circumferential surface crack and TPB specimen machined from same piping material (SA333Gr6 Steel) are tested. Consequently, 3D finite element analysis (FEA) has been performed on surface cracked pipe and TPB specimen. Crack‐initiation load is also predicted for surface cracked pipe by FEA and compared with experimental result. J–R curve is calculated for the pipe using experimental data, that is, load, load line displacement and crack growth. J–R curve of pipe is compared with TPB specimen and it is found that the pipe is predicting much higher J–R curve than TPB. This difference of J–R curve is investigated by evaluating stress triaxiality in remaining ligament for both cases. Stress triaxiality is quantified using triaxiality factor (h) ahead of crack tip for pipe and TPB specimen. It is found that the TPB specimen has considerably higher constraint level than pipe with surface crack, which is well supported by trend of J–R curves for specimen and pipe. A study has also been carried out to investigate the effect of internal pressure on the stress triaxiality. It is found that there is negligible difference in stress triaxiality because of internal pressure. The stress triaxiality is re‐established as a qualitative parameter to assess the transferability of J–R curve from specimen to component.     

Effect of strength mis-match and dynamic loading on ductile fracture initiation

It has been well known that ductile fracture of steels is accelerated by triaxial stresses. The characteristics of ductile crack initiation in steels are evaluated quantitatively using a two-parameter criterion based on equivalent plastic strain and stress triaxiality.The present study focuses on the effects of geometrical discontinuity, strength mis-match, which can elevate plastic constraint due to heterogeneous plastic straining, and loading rate on the critical condition for ductile fracture initiation using a two-parameter criterion. Fracture initiation testing has been conducted under static and dynamic loading using circumferentially notched round-bar specimens. In order to evaluate the stress/strain state in the specimens, especially under dynamic loading, a thermal elastic-plastic dynamic finite element (FE) analysis considering the temperature rise due to plastic deformation has been carried out.The tensile tests on specimens with an undermatching interlayer showed that the relationship between the critical equivalent plastic strain to initiate ductile fracture and stress triaxiality was equivalent to that obtained on homogeneous specimens under static loading. Moreover, the two-parameter criterion for ductile fracture initiation is shown to be independent of the loading rate. It was demonstrated that the critical global strain to initiate ductile fracture in specimens with strength mis-match under various loading rate can be estimated based on the local criterion, that is two-parameter criterion obtained on homogeneous specimens under static tension, by mean of FE-analysis taken into account accurately both strength mis-match and dynamic loading effects on stress/strain behaviors.     

Stress state and stress‐induced microstructural evolution around the crack tip of G115 steel after dwell‐fatigue crack propagation

Many components operate under creep‐fatigue loading causing an increasing need to learn the knowledge about creep‐fatigue crack growth (CFCG). The stress and microstructural evolution around the crack tip of G115 steel after CFCG were investigated. According to the finite element simulations, the variations of equivalent stress and stress triaxiality in the crack tip zone are presented. Furthermore, the evolutions of martensitic laths, dislocations, and precipitates were systematically studied through electron backscatter diffraction and transmission electron microscopy observations. The laths at the crack tip or under a larger hold time are wider than those remote from the crack tip or under a shorter dwell time. Meanwhile, the dislocation densities reduced significantly at the crack tip or under a larger hold time. The different variations of laths' width and dislocation densities resulted from the different stress triaxiality and creep strain. W‐rich Laves and Cu‐rich phases appeared during CFCG. The Laves phase coarsened rapidly because of the stress‐accelerated diffusion.     

Effect of stress triaxiality levels in crack tip regions on the characteristics of void growth and fracture criteria

The behaviour of void growth in the crack tip regions of four specimen geometries with different stress triaxiality levels have been investigated by the FEM method and experimental observations in plane strain and plane stress cases respectively. It was found that the shape change of growing voids, the configurations of a blunting crack tip and the sizes of decreasing ligament between the void and the crack tip are strongly dependent upon the stress triaxiality levels. Under the condition of plane stress, the stress triaxiality on the ligaments of cracks are nearly the same for different specimen geometries, also the void growth, crack tip blunting asnd decreasing of ligament size are identical for various specimens with increasing load levels, which lead to the conclusion that the J_i-value is independent of specimen geometries. However, in the plain strain case, different void growth, crack tip blunting and decreasing in ligament size for various stress triaxiality levels directly caused the J_i-value to be dependent on the specimen geometries. It was found that when the void is linked to the blunting crack tip by the extrapolation to the zero ligament from FEM calculations, the J_i values, measured experimentally, are underestimated slightly.     

Finite element simulation of fracture behaviour for aged duplex stainless steels

In this paper, the local approach model developed by Gurson–Tvergaard has been applied to simulate both the crack initiation and the crack growth of aged duplex stainless steel. The parameters of the Gurson–Tvergaard model have been obtained, from axisymmetric notched specimen testing, as a function of the ageing time at 400°C, the ferrite content of the steel and the stress triaxiality. After that, to simulate the fracture of CT specimens, finite element (FE) calculations have been effected in order to obtain the stress triaxiality value at each point on the process zone ahead of the crack tip of these specimens. The adequate damage parameters concerning triaxiality are determined from the ones obtained at the notched specimens, in order to be used in FE simulations of fracture behaviour. With them, the corresponding J−Δa curves have been simulated as representative of both the crack initiation and crack propagation stages, and compared with experimental results in order to validate the methodology proposed.     

Off-axis fatigue cracking behaviour in notched fibre metal laminates

The off‐axis fatigue cracking behaviour of notched fibre metal laminates under constant amplitude loading conditions was investigated experimentally and numerically. It was found that the off‐axis fatigue crack initiation life decreased as the off‐axis angles increased. This indicated that the off‐axis laminates raised the applied stress level in the aluminium (Al) layer and subsequently resulted in earlier cracking in the Al layer. The off‐axis fatigue crack initiation lives of notched fibre metal laminates were predicted using lamination theory and an energy‐based critical plane fatigue damage analysis from the literature. After a crack initiated in the Al layer, it was observed that the crack path angles of the off‐axis specimens were neither perpendicular to the fibre nor to the loading direction. A finite‐element model was established for predicting the crack path angles.     

Mechanical heterogeneity and validity of J-dominance in welded joints

Fracture criterion of the J-integral finds wide application in the integrity evaluation of welded components, but there exist some confused problems such as the dependence of the fracture toughness on the strength mis-matching and specimen geometry which need to be clarified. It is rough and unsuitable to attribute the variation of J-integral fracture parameter simply to the effect of mechanical heterogeneity. In the present paper, a two-dimensional finite element method is employed to analyze the distribution and variation of crack tip field of welded joints with different strength mis-matching in four kinds of specimen geometry, and then the validity of J-dominance in welded joints is investigated. It is found that the crack tip field of mis-matched joint is different from that of either the weld metal or base metal of which the joint is composed, but it is situated between those of weld metal and base metal. Under the plane strain, there is obvious difference in stress triaxiality for different strength mis-matched joints. The validity of J-dominance in welded joint can not be obtained by comparing whether the stress triaxiality meets that required by the HRR solution because of the existence of mechanical inhomogeneity. By ascertaining if the stress triaxiality of welded joint near the crack tip is dependent of specimen geometry, the conclusion can be arrived at: for plane stress the validity of J-dominance is valid, whilst for plane strain the validity of J-dominance is lost. Based on the above, attempt has been made to point out that the influence of mechanical heterogeneity on the fracture toughness of weldment arises from the variation of constraint intensity-crack tip stress triaxiality. Compared with the effect of mechanical heterogeneity on the stress triaxiality, the losing of validity of J-dominance in mis-matched joint under plane strain may play a more critical role in the variation of J-integral fracture parameter of weldment.     

Mesoscopic mechanisms in fatigue crack initiation in an aluminium alloy

The mechanistic aspects of process of initiation of a mode‐I fatigue crack in an aluminium alloy (AA 2219‐T87) are studied in detail, both computationally as well as experimentally. Simulations are carried out under plane strain conditions with fatigue process zone modelled as stress‐state–dependent cohesive elements along the expected mode‐I failure path. An irreversible damage parameter that accounts for the progressive microstructural damage due to fatigue is employed to degrade cohesive properties. The simulations predict the location of initiation of the fatigue crack to be subsurface where the triaxiality and the opening tensile stresses are higher in comparison with that at the notch surface. Examination of the fracture surface profile of fracture test specimens near notch tip reveals a few types of regions and existence of a mesoscopic length scale that is the distance of the location of highest roughness from the notch root. A discussion is developed on the physical significance of the experimentally observed length scale.     

Three-dimensional modeling of ductile crack growth: Cohesive zone parameters and crack tip triaxiality

For 10 mm thick smooth-sided compact tension specimens made of a pressure vessel steel 20MnMoNi55, the interrelations between the cohesive zone parameters (the cohesive strength, T_max, and the separation energy, Γ) and the crack tip triaxiality are investigated. The slant shear-lip fracture near the side-surfaces is modeled as a normal fracture along the symmetry plane of the specimen. The cohesive zone parameters are determined by fitting the simulated crack extensions to the experimental data of a multi-specimen test. It is found that for constant cohesive zone parameters, the simulated crack extension curves show a strong tunneling effect. For a good fit between simulated and experimental crack growth, both the cohesive strength and the separation energy near the side-surface should be considerably lower than near the midsection. When the same cohesive zone parameters are applied to the 3D model and a plane strain model, the stress triaxiality in the midsection of the 3D model is much lower, the von-Mises equivalent stress is distinctly higher, and the crack growth rate is significantly lower than in the plane strain model. Therefore, the specimen must be considered as a thin specimen. The stress triaxiality varies dramatically during the initial stages of crack growth, but varies only smoothly during the subsequent stable crack growth. In the midsection region, the decrease of the cohesive strength results in a decrease of the stress triaxiality, while the decrease of the separation energy results in an increase of the triaxiality.     

ANALYSIS AND MODELISATION OF SHORT CRACK GROWTH BY DUCTILE FRACTURE MICROMECHANISMS

Abstract— A ductile medium strength steel has been modelled by means of the Gurson model, and been used to investigate the effect of crack tip constraint in several fracture mechanics specimens. Both numerical and experimental results have been obtained, in the course of the crack extension process, for single edge notch bending specimens with different crack length-to-width ratios. The geometries with the shorter cracks always exhibited higher J values at initiation and steeper J crack growth resistance curves, and these results have been explained in terms of the stress and strain fields and damage development in the region ahead of the crack tip.     

裂纹尖端形变功与韧性启裂

研究了高强度结构钢典型强化组织板条马氏体、粒状贝氏体和板条马氏体加粒状贝氏体混合组织的韧性启裂.结果表明;高强度结构钢的韧性启裂依赖于强度、塑性和形变硬化能力,表现为裂纹尖端应力应变场中消耗的弹性应变能和塑性形变功不同.其相对量可用无量纲参量H表示,并与材料的流变应力有较好的对应关系.     

Effects of short cracks produced at blunted pre-crack tip on stress and strain distributions

The present work investigates problems: (1) How are the plastic strain and the stress (triaxiality) re-distributed after a short crack initiated, extended and blunted at the pre-crack tip? (2) How do the above changes put a crucial effect on the triggering of the cleavage fracture? Based on the previous observations of configuration changes and fracture surfaces of pre-crack tips, Finite element method (FEM) simulations of a short crack initiated, extended and blunted at a pre-crack tip and calculations of distributions of stress, strain and triaxiality are carried out for 3PB pre-cracked HSLA steel specimens tested at -130^°C. The results reveal that: as long as the fatigue pre-crack is only blunted, in its vicinity a region where the accumulated strain is sufficient to nucleate a crack, and a region where the stress (triaxiality) is sufficient to propagate a crack nucleus are separated by a distance. The nucleated crack cannot be propagated and the cleavage fracture cannot be triggered. While a short crack produced at the fully blunted fatigue pre-crack, the strain retains, the stress (triaxiality) is rebuilt. An initiated and significantly extended and then blunted short crack makes a tip configuration, which on one hand is much sharper than that of the fully blunted original pre-crack tip, on other hand is wide enough to spread its effects into the high stress covered region. This sharpened crack tip configuration re-builds a ‘sharper’ distribution of stress (triaxiality) and makes two regions metioned above closer. Finally the two regions overlap each other and a cleavage crack can be initiated and propagated at a distance ahead of the blunted fatigue pre-crack.     

Fatigue crack growth law of API X80 pipeline steel under various stress ratios based on J‐integral

Load‐controlled three‐point bending fatigue tests were conducted on API X80 pipeline steel to investigate the effects of stress ratio and specimen orientation on the fatigue crack growth behaviour. Because of the high strength and toughness of X80 steel, crack growth rate was measured and plotted versus ΔJ with stress ratio. The fatigue crack length is longer in the transverse direction, whereas the fatigue crack growth rates are nearly the same in different orientations. Finally, a new fatigue crack growth model was proposed. The effective J‐integral range was modified by ΔJ_p in order to correlate crack closure effect due to large‐scale yield of crack tip. The model was proved to fit well for fatigue crack growth rate of API X80 at various stress ratios of R > 0.     

Dependence of ductile crack formation in tensile tests on stress triaxiality, stress and strain ratios

An experimental and numerical study on ductile crack formation in tensile tests was conducted. Five different specimens including flat specimens, smooth round bars, notched bars (two types) and flat-grooved plates were investigated. Von Mises equivalent strain to crack formation, stress triaxiality, and stress and strain ratios at critical locations, were obtained. Accuracy of the Bridgman formulas for stresses in necked round bars, and McClintock's model for flat-grooved plates, were studied. A relationship between the stress triaxiality and equivalent strain to crack formation was determined in a high stress triaxiality range for Al 2024-T351. More importantly, it was found that equivalent strain and stress triaxiality are the two most important factors governing crack formation, while stress and strain ratios cause secondary effects. It appears possible to make a good prediction of crack formation with equivalent strain and stress triaxiality.