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.     

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.     

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.     

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

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

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.     

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.     

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.     

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.     

Photoelastic determination ofK I andK II : A numerical study on experimental data

The most efficient photoelastic methods to obtain stress intensity factors are those based on stress functions series expansions. The coefficients of these expansions are fitted to the experimental isochromatic pattern using an overdeterministic Newton-Raphson least squares method. In this paper, a study has been carried out to analyze the influence on the results of several numerical and experimental factors. It is shown that accurate values of the stress intensity factorsK _I andK _II can be obtained by following some recommendations given in the text and summarized in the conclusions at the end of the paper.     

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.     

Classification of stress-intensity factors from isochromatic-fringe patterns

The stresses in the local neighborhood of a crack tip have been used to develop a relation between the isochromatic-fringe orderN, its position parametersr and θ and the stress field expressed in terms of stress intensities,K _I ,K _II , and a far-field stress σ _ox . This relation was programmed and a plotting routine was developed to map isochromatic (σ₁ ? σ₂) fields in the neighborhood of the crack tip. The stress intensitiesK _I andK _II and the far-field stress σ _ox were varied and isochromatic fields were constructed for each combination. As bothK _II and σ _ox influence the size, shape and orientation of the isochromatics loops in a systematic manner, the pictorial representation of the isochromatic fields can be used to classify the state of stress (K _I ,K _II and σ _ox ) at the crack tip. Isochromatics which classify six different states of stress have been illustrated and methods used to determineK _I ,K _II , andσ _ox in five of the six states are given.     

An Experimental Study of Mixed Mode Crack Initiation and Growth in Functionally Graded Materials

Quasi-static mixed mode crack initiation and growth in functionally graded materials (FGMs) was studied through fracture experiments on polymer-based FGMs manufactured by selective ultraviolet irradiation poly(ethylene carbon monoxide)—a photo-sensitive copolymer that becomes more brittle and stiffer under ultraviolet irradiation. The objective of the study was to determine whether crack kinking criteria for homogeneous materials, e.g., maximum hoop stress criterion, also hold for FGMs. Single edge notched tension specimens with different spatial variations of Young's modulus, failure stress and failure strain, were tested. Near tip mode mixity was introduced either by inclining the crack to the remote loading direction, as in the case of homogeneous materials, or to the direction of material gradient, or both. A full-field digital image correlation technique was used to measure in real-time the displacement field around the crack tip while it propagated through the graded material, and to extract the fracture parameters of stress intensity factor K _I and K _II , and the T-stress. It was found that the nonsingular T-stress term in the asymptotic expansion for stresses plays a very important role in accurately measuring fracture parameters. It was also found that the maximum tangential stress criterion can be applied to the case of FGMs to predict crack kinking provided that the effect of the T-stress is accounted for and the process zone size is small compared to the intrinsic material gradient length scale. However, for accurate crack path prediction at a length scale comparable to the material gradient, detailed material property information is required. In general, the crack will propagate towards a region that exhibits less fracture toughness, but, unlike the case of homogeneous materials, along a path where K _II is not necessarily equal to zero.     

A crack perpendicular to and terminating at a bimaterial interface

Using dislocation simulation approach, the basic equation for a finite crack perpendicular to and terminating at a bimaterial interface is formulated. A novel expansion method is proposed for solving the problem. The complete solution to the problem, including the explicit formulae for theT stresses ahead of the crack tip and the stress intensity factors are presented. The stress field characteristics are analysed in detail. It is found that normal stresses {ie27-1} and {ie27-2} ahead of the crack tip, are characterised byQ fields if the crack is within a stiff material and the parameters |p _T | and |q _T | are very small, whereQ is a generalised stress intensity factor for a crack normal to and terminating at the interface. If the crack is within a weak material, the normal stresses {ie27-3} and {ie27-4} are dominated by theQ field plusT stress. This work was supported by the Swedish Research Council for Engineering Sciences.     

Analysis of dynamic mixed-mode isochromatics

The mixed-mode, elastodynamic state of stress in the neighborhood of a constant-velocity crack tip is used to generate numerically unsymmetric isochromatics. Unsymmetry associated with the third-order terms of a mixed-mode stress field, with and without the Mode II singular stress term, is also investigated. In extractingK _I from an unsymmetric isochromatic pattern, errors in the Mode I fracture parameters due to the assumed presence ofK _II in aK _I stress field were found to be significant when data are taken more than 4 mm from the crack tip. Paper was presented at V International Congress on Experimental Mechanics held in Montreal, Quebec, Canada on June 10–15, 1984.     

Stress state in the vicinity of an inclined elliptical defect and stress intensity factors for biaxial loading of a plate

The problem of determining the stress state of a plate with an inclined elliptical notch under biaxial loading is considered. The Kolosov-Muskhelishvili method is used to obtain an expression for the stress near the vertex of an inclined ellipse, whose particular case are expressions for the stress in the case of an inclined crack. The stress intensity factors K _I and K _II were determined experimentally by holographic interferometry in the case of extension of a plate with an inclined crack-like defect. The calculation results are compared with experimental data. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 1, pp. 118–127, January–February, 2009.     

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.     

The effect of contact stress intensity factors on fatigue crack propagation

Using the technique of Dimensional Analysis the phenomenon of crack closure is modelled using the concept of a contact stress intensity factor Kc. For constant amplitude loading, a simple expression, Kc_max = g(R) ΔK, is obtained without making idealized assumptions concerning crack tip behaviour. Further, by assuming that crack closure arises from the interaction of residual plasticity in the wake of the crack and crack tip compressive stresses, the function g(R) is shown to be constant for non-workhardening materials. This implies that any dependency of Kc_max on R must be attributed to the workhardening characteristic of the material. With Kc known, an “effective” stress intensity factor Ke may be calculated and incorporated into a crack growth law of the form da/dn = f(ΔKe). From analysis, it can be deduced that for a workhardening material, Kc_max will decrease as R increases and the effective stress intensity factor will increase. This means that the fatigue crack propagation rate will increase with R, in accordance with experimental observations.     

Simultaneous measurements of stress intensity and toughness for fast-running cracks in steel

Simultaneous measurements of the dynamicstress-intensity factorK _I ^dyn and the dynamic-fracture toughnessK _ID were made in a high-strength steel to investigate the relation between energy delivered to and energy absorbed by rapidly propagating cracks. Values ofK _I ^dyn were obtained intermittently during the propagation history by the shadow optical method of caustics from high-speed photographs of the moving crack tips. Values ofK _ID were calculated from temperature maxima recorded by thermocouples near the crack path. The results indicate that for fast-running cracks, the change in energy available at the crack tip can be significantly less than the energy absorbed in crack extension, suggesting that currently used dynamic-energy-balance methods for determining dynamic-fracture toughnesses may provide erroneous values.     

Semi-weight function method on computation of stress intensity factors in dissimilar materials

Semi-weight function method is developed to solve the plane problem of two bonded dissimilar materials containing a crack along the bond. From equilibrium equation, stress and strain relationship, conditions of continuity across interface and free crack surface, the stress and displacement fields were obtained. The eigenvalue of these fields is lambda. Semi-weight functions were obtained as virtual displacement and stress fields with eigenvalue-lambda. Integral expression of fracture parameters, KⅠ and KⅡ, were obtained from reciprocal work theorem with semi-weight functions and approximate displacement and stress values on any integral path around crack tip. The calculation results of applications show that the semi-weight function method is a simple, convenient and high precision calculation method.     

Study of pure and mixed-mode fracture of a brittle material

In order to study the behavior of a crack in a linear-elastic material in plane mixed mode (modes I+II), a specimen's shape and loading have been specially adapted. The specimens are first precracked in mode I and then subjected to monotonic loading until instability is reached by an original device which makes it possible to control the nature of the mixed mode applied and which is adjustable from pure mode I to pure mode II. After the specimens are fractured, the lengths of the initial cracks and the kinking angles were measured in the plane-strain area. Then the stress-intensity factorsK _I andK _II at the moment when crack instability appeared were calculated. This made it possible to apply two criteria: maximum principal stress and maximum stressintensity factorK ₁ ^* at the onset of kinking. From comparing the calculated values with the experimental values we may note that there is good agreement with respect to the crackkinking directions. However, for the limit load values considerable divergences have been recorded which are analyzed.