Determination of the dynamic stress intensity factorsK I d andK II d,for a mixed-mode propagating crack

In this paper, the dynamic propagation problem of a mixed-mode crack was studied by means of the experimental method of caustics. The initial curve and caustic equations were derived under the mixed-mode dynamic condition. A multi-point measurement method for determining the dynamic stress intensity factors,K _I ^d , andK _II ^d , and the position of the crack tip was developed. Several other methods were adopted to check this method, and showed that it has a good precision. Finally, the dynamic propagating process of a mixed-mode crack in the three-point bending beam specimen was investigated with our method.     

Thermoelastic investigations for fatigue life assessment

An investigation is presented on the suitability and accuracy of a thermoelastic technique for the analysis of fatigue cracks. The stress intensity factor ranges ΔK _I and ΔK _II are determined from thermoelastic data recorded from around the tip of a sharp slot in a steel specimen under biaxial load, in order to assess the accuracy of the technique. ΔK _I and ΔK _II are determined to within 4% and 9% of a theoretical prediction, respectively. The results from a similar test on a fatigue crack under biaxial load are also presented. These show that thermoelastic stress analysis is a rapid and accurate way of analyzing mixed-mode fatigue cracks. A discussion is given on the potential of thermoelastic stress analysis of propagating cracks.     

FEM analysis on mixed-mode fracture of CSM-GRP

A FEM analysis for studying mixed-mode fracture problem of chopped strand mat glass fibre reinforced polyester laminate is presented. The analysis is formulated on the basis of 8-node quadrilateral isoparametric element. The collapsed triangular quarter-point singular elements were used for calculating stress intensity factors K_Ι and K_Ⅱ.The crack propagation process was computed by implementing constraint release technique. Three different approaches to the solution of stress intensity factors K_Ι and K_Ⅱ were compared. The effect of constraint condition imposed upon the displacement of the three collapsed nodes of the crack tip elements on the K_Ι and K_Ⅱ results was evaluated. The mixed-mode critical stress intensity factors K_ΙC and K_ⅡC were estimated for CSM-GRP through the consideration of K_Ι and K_Ⅱ calculated and the measured failure load and critical crack length in the experiment.     

Plasticity aspects of interfacial fracture for polymer/glass specimens

Experimental results suggest that the interfacial fracture resistance is minimal for approximate near tip Mode I accompanied by positive and negative near tip Mode II. Finite-strain FE analysis is made for an elastic–plastic medium bonded to an ideally elastic medium with an interface crack. Small-scale plasticity conditions are invoked and examined in relation to the elastic–plastic stress distribution along the bond line. Plasticity engenders a tendency to turn near tip biaxiality towards pure Mode I regardless of the mixed-mode loading. High levels of hydrostatic stress are attained. For different mode mixities of the applied load, the dependence of the elastic–plastic normal bond stress on load level is examined. It is found that under positive Mode II loading, the normal bond stress σ_yy tends to saturate as the load level rises. This does not occur for Mode I and negative Mode II loading. In addition, deformation patterns inside the plastic zone are examined for mixed-mode situations. A displacement criterion based on the normal bond crack opening suggests a dependence of the critical load level on the extent of mixed mode. Under positive mode II fracture, traces of the ductile material are found at the top of the elastic substrate. Some of these conclusions appear to be consistent with the fracture patterns observed for LD-polyethylene/glass interfacial mixed-mode fracture.     

Resolved shear stress intensity coefficient and fatigue crack growth in large crystals

Fatigue crack growth is caused primarily by shear decohesion due to dislocation motion in the crack tip region. The resolved shear stress, which drives dislocation in a crystal, is strongly orientation dependent, and therefore, the cyclic plastic deformation of the shear decohesion process is highly anisotropic.The crack planes are often inclined to the loading axis both in the inplane orientation and in the thickness direction. This inclination induces all three modes of the crack tip stress field, K_I, K_II, and K_III.Fatigue crack growth in large-grain Al 7029 aluminum alloy was studied. The crack tip stress fields of the test specimens are calculated with the finite element method. The values of K_I, K_II, and K_III are evaluated. The orientation of the crystal at a crack tip was determined with the Laue X-ray method. The crystal orientation and the calculated crack tip stress fields are used to compute the resolved shear stress intensity of each of the twelve slip systems of the crystal at the crack tip. The resolved shear stress field of a slip system is linearly proportional to the resolved shear stress intensity coefficient, RSSIC.The values of RSSIC thus evaluated are used to analyze the orientations of the crack plane and to correlate with the shear fatigue crack growth rate.     

On the higher order asymptotic analysis of a non-uniformly propagating dynamic crack along an arbitrary path

Transient mixed-mode elastodynamic crack growth along arbitrary smoothly varying paths is considered. Asymptotically, the crack tip stress field is square root singular with the angular variation of the singular term depending weakly on the instantaneous values of the crack tip speed and on the mode-I and mode-II stress intensity factors. However, for a material particle at a small distance away from the moving crack tip, the local stress field will depend not only on the instantaneous values of the crack tip speed and stress intensity factors, but also on the past history of these time dependent quantities. In addition, for cracks propagating along curved paths the stress field is also expected to depend on the nature of the curved crack path. Here, a representation of the crack tip fields in the form of an expansion about the crack tip is obtained in powers of radial distance from the tip. The higher order coefficients of this expansion are found to depend on the time derivative of crack tip speed, the time derivatives of the two stress intensity factors as well as on the instantaneous value of the local curvature of the crack path. It is also demonstrated that even if cracks follow a curved path dictated by the criterion K ₁₁ ^d =0, the stress field may still retain higher order asymmetric components related to non-zero local curvature of the crack path.     

Plane-strain state of an elastoplastic body with a crack under mixed-mode loading

A mixed-mode (I + II) crack model with a plastic strip on its continuation under plane strain is proposed. The stress components within the strip are determined from the yield conditions, stress limitation, and relationship between the normal stress components defined via the principal stress state. The crack parameters are analyzed for the Mises yield condition. In the quasibrittle case, the governing system of equations includes stress intensity factors K _I, K _II, and T-stresses     

Effect of residual surface stress on the fracture of nanoscale materials

The problem of a mode I crack in nanomaterials under a remote mechanical load is investigated. The effect of the residual surface stress on the crack surface is considered and the solutions to the crack opening displacement (COD) and the stress intensity factor (K_I) are obtained. The results show that the surface effect on the crack deformation and crack tip field are prominent at nanoscale. Moreover, COD and K_I are influenced by the residual surface stress not only on the surface near the crack tip region but also on the entire crack surface.     

On Cauchy singular integrals and stress intensity factors for 2D mode I cracks in an infinite solid

This paper comments on some of the different numerical techniques commonly employed in evaluating Cauchy singular integrals of the first kind; e.g. as pertaining to 2D through cracks in a brittle material undergoing Mode I loading. In addition, a different more direct method is proposed here. Also, two different ways to calculate the stress intensity factor (K_I) are contrasted. The accuracy attained by the different methods in calculating K_I, and the factors affecting the calculation, are compared. Finally, comments on calculating the stress field of a 2D crack and important considerations are presented.     

A study ofJ-integral of the orthotropic composite material

The relation between J-integral near model I crack tip in the orthotropic plate and displacement derivative is derived in this paper. Meanwhile, the relation between stress intensity factor K _I and displacement is also given in this paper. With sticking film moire interferometry method, the three-point bending beam is tested, thus the values of J-integral and K _I can be obtained from the displacement field, and then the truth of relation formula between J-integral and K _I in the orthotropic composite materials is experimentally verified.     

Stress intensity factors of an oblique edge crack subjected to normal and shear tractions

Stress intensity factors (SIFs) were obtained for an oblique crack under normal and shear traction and remote extension loads. The oblique crack was modeled as the pseudodislocation. The stress field due to tractions was solved by the Flamant solution. The SIR of Mode I and Mode II (K_Iand K_II) were then obtained. Finite element analysis was performed with ABAQUS and compared with the analytical solutions. The analytical solutions were in good agreement with the results of FEM. From investigating SIFs and their ranges, the following results were obtained. The growth rate of an oblique edge crack decreased due to the reduction in the SIF ranges. The crack driving force depended on the obliquity, the normal traction and the ratio of crack to traction length. The peak value of shear traction was found as a key parameter to accelerate the crack growth.     

Mixed-mode crack growth in plates under three-point bending

A combined theoretical and experimental study of the problem of crack growth in a plate subjected to unsymmetrical three-point bending was undertaken. The opening-modeK _I and sliding-modeK _II stress-intensity factors describing the local stress field around the crack tip were determined by a finite-element computer program. The crack growth was analyzed by the maximum circumferential stress and the minimum strain-energy density criteria. The critical loads for crack growth and the crack trajectories were determined both by theory and experiment. The experimental results corroborated the theoretical predictions.     

An experimental analysis of displacement field of mixed mode crack

In this paper, the characteristic properties ofv (y-direction displacement) field surrounding the tip of a mixed mode crack are studied. These properties can be used to evaluate the rigid body rotation of the crack tip, theK _I SIF and the ratio ofK _II SIF toK _I.The authors employ a film to record the displacement information based on the technique of moire interferometry with sticking films. By using the data taken from the moire pattern and treating them with the damping least square method, all of the parameters of the crack can be obtained accurately.     

Anti-plane shear of kinked interface crack in bonded dissimilar anisotropic solid

Complex potentials are derived to describe the anti-plane singular shear stress fields around a kinked crack, the main portion of which is embedded along the interface of two dissimilar anisotropic elastic media. This is accomplished by formulating the problem as singular integral equations with generalized Cauchy kernels. The shear stress singularity at the kink differs from the familiar inverse square root of the local distance; it is found to influence the magnitude of the Mode III crack tip stress intensity factor, K₃. Numerical results of K₃ are obtained and displayed in graphical forms for different degree of material anisotropy and crack dimensions.     

Numerical simulation of stress intensity factors via the thermoelastic technique

The purpose of this study is to investigate the accuracy of the least squares method for finding the in-plane stress intensity factorsK _I andK _II using thermoelastic data from isotropic materials. To fully understand the idealized condition ofK _I andK _II calculated from thermoelastic experiments, the total stress field calculated from finite element analysis is used to take the place of data obtained from real thermoelastic experiments. In the finite element analysis, theJ-integral is also calculated to compare with (K _I ² +K _II ² )/E evaluated by the least squares method. The stress fields near the crack tip are dominated by the two stress intensity factors; however, the edge effect will cause inaccuracy of the thermoelastic data near the crack tip. Furthermore, the scan area of thermoelastic experiments cannot be too small. Therefore, we suggest that three or four terms of stress function be included in the least squares method for evaluating stress intensity factors via the thermoelastic technique. In the idealized condition, the error can be smaller than 3 percent from our numerical simulations. If only ther ^–1/2 term (K _I andK _II ) is included in the least squares method, even in the idealized case the error can be up to 20 percent.     

Axial crack propagation and arrest in a pressurized cylinder: An experimental-numerical analysis

The feasibility of using a previously developed crack-kinking criterion to predict crack arrest at a tear strap in a pressurized fuselage was studied with instrumented axial rupture tests of 21 models of an idealized fuselage. A rapidly propagating axial crack, which was initiated from a precrack, kinked immediately upon extension and propagated diagonally until it turned circumferentially and propagated along the tear straps. An elastodynamic finite element analysis of the rupturing model fuselage yielded the mixed-mode stress intensity factors,K _I andK _II , and the remote stress component, σ _OX . This numerical procedure was also used to predict the crack trajectories in full-scale fuselage rupture tests. All numerical results agreed well with their measured counterparts regardless of size.     

Crack front curvature effect on stress intensity factor of through and part-through thickness cracks

The character of the local stresses and displacements are determined for a through crack with finite radius of curvature in a finite thickness plate. Numerical results obtained from the boundary element method show that the solutions are sufficiently accurate for /a ≤ 0.03 and 0.03 ≤ /a ≤ 0.1, where and a represent, respectively, the crack front radius of curvature and crack dimension such that a is the width of a through thickness crack and the depth of a part-through crack. For /a ≤ 0.1, the asymptotic singular stress field dominates such that the Mode I stress intensity factor K₁ can be evaluated. As the crack border radius of curvature is increased for /a ≥ 0.1, the non-singular terms become significant such that K_I would no longer dominate. Other failure criteria would have to be invoked to address fracture initiation.     

A method for measuring mode I crack tip constraint under static and dynamic loading conditions

A novel experimental technique for measuring crack tipT-stress, and hence in-plane crack tip constraint, in elastic materials has been developed. The method exploits optimal positioning of stacked strain gage rosette near a mode I crack tip such that the influence of dominant singular strains is negated in order to determineT-stress accurately. The method is demonstrated for quasi-static and low-velocity impact loading conditions and two values of crack length to plate width ratios (a/W). By coupling this new method with the Dally-Sanford single strain gage method for measuring the mode I stress intensity factorK _I , the crack tip biaxiality parameter is also measured experimentally. Complementary small strain, static and dynamic finite element simulations are carried out under plane stress conditions. Time histories ofK _I andT-stress are computed by regression analysis of the displacement and stress fields, respectively. The experimental results are in good agreement with those obtained from numerical simulations. Preliminary data for critical values ofK _I and β for dynamic experiments involving epoxy specimens are reported. Dynamic crack initiation toughness shows an increasing trend as β becomes more negative at higher impact velocities.     

钢轨表面剥离掉块路径预测研究

轮轨滚动接触下,钢轨表面会产生典型的鱼钩形剥离掉块,其形成机理目前暂未明确.为了探究轮轨滚动接触下钢轨表面裂纹扩展机理,基于最大周向拉应力准则,建立轮轨滚动接触疲劳计算模型,提出裂尖扩展路径预测方法,并对不同初始角度裂纹的扩展路径进行预测.结果表明,钢轨表面微裂纹为Ⅰ-Ⅱ复合型裂纹,随着裂纹长度增加,KⅠ先增加后减小,...     

Strain-gage methods for measuring the opening-mode stress-intensity factor,K I

Measurements of strain near a crak tip with electrical-resistance strain gages do not usually provide a reliable measure ofK _I because of local yielding, three-dimensional effects and limited regions for strain-gage placement. This paper develops expressions for the strains in a valid region removed from the crack tip, and indicates procedures for locating and orienting the gages to accurately determineK _I from one or more strain-gage readings.