Effect of constraint on void growth near a blunt crack tip

The growth of a single cylindrical hole ahead of a blunt crack tip was studied using large deformation finite element analysis in three-point bend specimens with different precrack depth. The effect of small second phase particles was taken into account by incorporating Gurson’s constitutive equation. The effects of strain hardening and the initial distance from the hole to the crack tip were also investigated. The results show that the variation of crack tip opening displacement with load is not sensitive to constraint level. The effects of constraint on the growth of hole and ductile initiation toughness are diminished with decreasing initial distance from the hole to the blunt crack tip.

Numerical investigation on stable crack growth in plane stress

Large deformation finite element analysis has been carried out to investigate the stress-strain fields ahead of a growing crack for compact tension .a=W D 0:5/ and three-point bend .a=W D 0:1 and 0:5/ specimens under plane stress condition. The crack growth is controlled by the experimental J -integral resistance curves measured by Sun et al. The results indicate that the distributions of opening stress, equivalent stress and equivalent strain ahead of a growing crack are not sensitive to specimen geometry. For both stationary and growing cracks, similar distributions of opening stress and triaxiality can be found along the ligament. During stable crack growth, the crack-tip opening displacement (CTOD) resistance curve and the cohesive fracture energy in the fracture process zone are independent of specimen geometry and may be suitable criteria for characterizing stable crack growth in plane stress.

Effect of constraint on ductile crack growth and ductile-brittle fracture transition of a carbon steel

Ductile-brittle fracture transition was investigated using compact tension (CT) specimens from -70oC to 40oC for a carbon steel. Large deformation finite element analysis has been carried out to simulate the stable crack growth in the compact tension (CT, a/W=0.6), three point-point bend (SE(B), a/W=0.1) and centre-cracked tension (M(T), a/W=0.5) specimens. Experimental crack tip opening displacement (CTOD) resistance curve was employed as the crack growth criterion. Ductile tearing is sensitive to constraint and tearing modulus increases with reduced constraint level. The finite element analysis shows that path-dependence of J-integral occurs from the very beginning of crack growth and ductile crack growth elevates the opening stress on the remaining ligament. Cleavage may occur after some ductile crack growth due to the increase of opening stress. For both stationary and growing cracks, the magnitude of opening stress increases with increasing in-plane constraint. The ductile-brittle transition takes place when the opening stress ahead of the crack tip reaches the local cleavage stress as the in-plane constraint of the specimen increases.

Effect of crack depth and specimen width on fracture toughness of a carbon steel in the ductile-brittle transition region

The effects of crack depth (a/W) and specimen width W on the fracture toughness and ductile±brittle transition have been investigated using three-point bend specimens. Finite element analysis is employed to obtain the stress-strain fields ahead of the crack tip. The results show that both normalized crack depth (a/W) and specimen width (W) affect the fracture toughness and ductile±brittle fracture transition. The measured crack tip opening displacement decreases and ductile±brittle transition occurs with increasing crack depth (a/W) from 0.1 to 0.2 and 0.3. At a fixed a/W (0.2 or 0.3), all specimens fail by cleavage prior to ductile tearing when specimen width W increases from 25 to 40 and 50 mm. The lower bound fracture toughness is not sensitive to crack depth and specimen width. Finite element analysis shows that the opening stress in the remaining ligament is elevated with increasing crack depth or specimen width due to the increase of in-plane constraint. The average local cleavage stress is dependent on both crack depth and specimen width but its lower bound value is not sensitive to constraint level. No fixed distance can be found from the cleavage initiation site to the crack tip and this distance increases gradually with decreasing inplane constraint.

An enriched radial point interpolation method (e-RPIM) for analysis of crack tip fields

In this paper, an enriched radial point interpolation method (e-RPIM) is developed the for the determination of crack tip fields. In e-RPIM, the conventional RBF interpolation is novelly augmented by the suitable trigonometric basis functions to reflect the properties of stresses for the crack tip fields. The performance of the enriched RBF meshfree shape functions is firstly investigated to fit different surfaces. The surface fitting results have proven that, comparing with the conventional RBF shape function, the enriched RBF shape function has: (1) a similar accuracy to fit a polynomial surface; (2) a much better accuracy to fit a trigonometric surface; and (3) a similar interpolation stability without increase of the condition number of the RBF interpolation matrix. Therefore, it has proven that the enriched RBF shape function will not only possess all advantages of the conventional RBF shape function, but also can accurately reflect the properties of stresses for the crack tip fields. The system of equations for the crack analysis is then derived based on the enriched RBF meshfree shape function and the meshfree weak-form. Several problems of linear fracture mechanics are simulated using this newlydeveloped e-RPIM method. It has demonstrated that the present e-RPIM is very accurate and stable, and it has a good potential to develop a practical simulation tool for fracture mechanics problems.

An enriched radial point interpolation method (e-RPIM) for analysis of crack tip fields

In this paper, an enriched radial point interpolation method (e-RPIM) is developed the for the determination of crack tip fields. In e-RPIM, the conventional RBF interpolation is novelly augmented by the suitable trigonometric basis functions to reflect the properties of stresses for the crack tip fields. The performance of the enriched RBF meshfree shape functions is firstly investigated to fit different surfaces. The surface fitting results have proven that, comparing with the conventional RBF shape function, the enriched RBF shape function has: (1) a similar accuracy to fit a polynomial surface; (2) a much better accuracy to fit a trigonometric surface; and (3) a similar interpolation stability without increase of the condition number of the RBF interpolation matrix. Therefore, it has proven that the enriched RBF shape function will not only possess all advantages of the conventional RBF shape function, but also can accurately reflect the properties of stresses for the crack tip fields. The system of equations for the crack analysis is then derived based on the enriched RBF meshfree shape function and the meshfree weak-form. Several problems of linear fracture mechanics are simulated using this newlydeveloped e-RPIM method. It has demonstrated that the present e-RPIM is very accurate and stable, and it has a good potential to develop a practical simulation tool for fracture mechanics problems.

Numerical solution of stress intensity factors of multiple cracks in great number with eigen COD boundary integral equations

Based on the eigen crack opening displacement (COD) boundary integral equations, a newly developed computational approach is proposed for the analysis of multiple crack problems. The eigen COD particularly refers to a crack in an infinite domain under fictitious traction acting on the crack surface. With the concept of eigen COD, the multiple cracks in great number can be solved by using the conventional displacement discontinuity boundary integral equations in an iterative fashion with a small size of system matrix. The interactions among cracks are dealt with by two parts according to the distances of cracks to the current crack. The strong effects of cracks in adjacent group are treated with the aid of the local Eshelby matrix derived from the traction BIEs in discrete form. While the relatively week effects of cracks in far-field group are treated in the iteration procedures. Numerical examples are provided for the stress intensity factors of multiple cracks, up to several thousands in number, with the proposed approach. By comparing with the analytical solutions in the literature as well as solutions of the dual boundary integral equations, the effectiveness and the efficiencies of the proposed approach are verified.

Efficient solution of multiple cracks in great number using eigen COD boundary integral equations with iteration procedure

A newly developed computational approach is proposed in the paper for the analysis of multiple crack problems based on the eigen crack opening displacement (COD) boundary integral equations. The eigen COD particularly refers to a crack in an infinite domain under fictitious traction acting on the crack surface. With the concept of eigen COD, the multiple cracks in great number can be solved by using the conventional displacement discontinuity boundary integral equations in an iterative fashion with a small size of system matrix to determine all the unknown CODs step by step. To deal with the interactions among cracks for multiple crack problems, all cracks in the problem are divided into two groups, namely the adjacent group and the far-field group, according to the distance to the current crack in consideration. The adjacent group contains cracks with relatively small distances but strong effects to the current crack, while the others, the cracks of far-field group are composed of those with relatively large distances. Correspondingly, the eigen COD of the current crack is computed in two parts. The first part is computed by using the fictitious tractions of adjacent cracks via the local Eshelby matrix derived from the traction boundary integral equations in discretized form, while the second part is computed by using those of far-field cracks so that the high computational efficiency can be achieved in the proposed approach. The numerical results of the proposed approach are compared not only with those using the dual boundary integral equations (D-BIE) and the BIE with numerical Green's functions (NGF) but also with those of the analytical solutions in literature. The effectiveness and the efficiency of the proposed approach is verified. Numerical examples are provided for the stress intensity factors of cracks, up to several thousands in number, in both the finite and infinite plates.

Fatigue crack propagation analysis of plaque rupture

Rupture of atheromatous plaque is the major cause of stroke or heart attack. Considering that the cardiovascular system is a classic fatigue environment, plaque rupture was treated as a chronic fatigue crack growth process in this study. Fracture mechanics theory was introduced to describe the stress status at the crack tip and Paris' law was used to calculate the crack growth rate. The effect of anatomical variation of an idealized plaque cross-section model was investigated. The crack initiation was considered to be either at the maximum circumferential stress location or at any other possible locations around the lumen. Although the crack automatically initialized at the maximum circumferential stress location usually propagated faster than others, it was not necessarily the most critical location where the fatigue life reached its minimum. We found that the fatigue life was minimum for cracks initialized in the following three regions: the midcap zone, the shoulder zone, and the backside zone. The anatomical variation has a significant influence on the fatigue life. Either a decrease in cap thickness or an increase in lipid pool size resulted in a significant decrease in fatigue life. Comparing to the previously used stress analysis, this fatigue model provides some possible explanations of plaque rupture at a low stress level in a pulsatile cardiovascular environment, and the method proposed here may be useful for further investigation of the mechanism of plaque rupture based on in vivo patient data.

Numerical analysis for the interlayer toughening of piezoelectric smart composite plate

The piezoelectric composite material could engender stress concentration resulting from small cracks during layers easily, as the cracks growth will lead to the failure of the whole structure. In this paper, a finite element model for piezoelectric composite materials by ABAQUS including interlayer crack was established, and the J integral and crack tip stress of different types PZT patches were calculated by using the equivalent integral method. Then, the J integral for adhesive layers with different thickness, elastic modulus considering and not considering piezoelectricity was investigated. The results show that the J integral of mode I, II reduces with thicker adhesive layer and lower elastic modules, and the J integral of mode II decreases more sharply than that of mode I.

A material model for flexural crack simulation in reinforced concrete elements using ABAQUS

This paper presents a material model to simulate load induced cracking in Reinforced Concrete (RC) elements in ABAQUS finite element package. Two numerical material models are used and combined to simulate complete stress-strain behaviour of concrete under compression and tension including damage properties. Both numerical techniques used in the present material model are capable of developing the stress-strain curves including strain softening regimes only using ultimate compressive strength of concrete, which is easily and practically obtainable for many of the existing RC structures or those to be built. Therefore, the method proposed in this paper is valuable in assessing existing RC structures in the absence of more detailed test results. The numerical models are slightly modified from the original versions to be comparable with the damaged plasticity model used in ABAQUS. The model is validated using different experiment results for RC beam elements presented in the literature. The results indicate a good agreement with load vs. displacement curve and observed crack patterns.

Fatigue crack growth and fatigue life evaluation for welded steel bridge members with initial crack

Initial crack widely exists in the welded members of steel bridge induced by the welding procedure or by the fatigue damage crack initiation. The behavior of crack growth with a view to fatigue damage accumulation on the tip of cracks is discussed. Fatigue life of welded components with initial crack in bridges under traffic loading is investigated. Based on existing fatigue experiment results of welded members with initial crack and the fatigue experiment results of welded bridge members under constant stress cycles, the crack would keep semi-elliptical shape with variable ratio of a/c during the crack propagation. Based on the concept of continuum damage accumulated on the tip of fatigue cracks,the fatigue damage law suitable for steel bridge members under traffic loading is modified to consider the crack growth.The virtual crack growth method and the semi-elliptical crack shape assumption are proposed in this paper to deduce a new model of fatigue crack growth rate for welded bridge members under traffic loading. And the calculated method of the stress intensity factor necessary for evaluation of the fatigue life of welded bridge members with cracks is discussed.The proposed fatigue crack growth model is then applied to calculate the crack growth and the fatigue life of existing welded members with fatigue experimental results. The fatigue crack propagation computation results show that the ratio of crack depth to the half crack surface length a/c is variable during crack propagation process and the stress cycle increases with the increase of a0/c0 with certain a0/t0 .The calculated and measured fatigue lives are generally in good agreement,at some initial conditions of cracking, for welded members widely used in steel bridges.

Experimental Study of crack propagation in carbon steel using acoustic emission

Acoustic emission technique has become a significant and powerful structural health monitoring tool for structures. Researches to date have been done on crack location, fatigue crack propagation in materials and severity assessment of failure using acoustic emission technique. Determining severity of failure in steel structures using acoustic emission technique is still a challenge to accurately determine the relationship between the severity of crack propagation and acoustic emission activities. In this study three point bending test on low carbon steel samples along with acoustic emission technique have been used to determine crack propagation and severity. A notch is introduced at the tension face of the loading point to the samples to initiate the crack. The results show that the percentage of load drop of the steel specimen has a reciprocal relationship with the crack opening i.e. crack opening zones are influenced by the loading rate. In post yielding region, common acoustic emission signal parameters such as, signal strength, energy and amplitudes are found to be higher than those at pre-yielding and at yielding.

Evaluation of material crack using acoustic emission technique

Cracks in civil structures can result in premature failure due to material degradation and can result in both financial loss and environmental consequences. This thesis reports an effective technique using Acoustic Emission (AE) technique to assess the severity of the crack propagation in steel structures. The outcome of this work confirms that combination of AE parametric analysis and signal processing techniques can be used to evaluate crack propagation under different loading configurations. The technique has potential application to assess and monitor the condition of civil structures.