Dynamic stress-intensity factors for unsymmetric dynamic isochromatics

The mixed mode, near-field state of stresses sourrounding a crack propagating at constant velocity is used to derive a relation between the dynamic stress-intensity factorsK _I,K _II, the remote stress component σ _ox and the dynamic isochromatics. This relation, together with an over-deterministic least-square method, form the basis of a datareduction procedure for extracting dynamic,K _I,K _II and σ _ox from the recorded dynamic photoelastic pattern surrounding a running crack. The overdeterministic least-square method is also used to fit static isochromatics to the numerically generated dynamic isochromatics. The resultant staticK _I,K _II and σ _ox are compared with the corresponding dynamic values and estimats of errors involved in using static analysis to process dynamic isochromatic data are obtained. The data-reduction procedure is then used to evaluate the branching stress-intensity factor associated with crack branching and the mixed-mode stress-intensity factors associated with crack curving.     

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.     

A Novel Experimental Approach for Calculating Stress Intensity Factors from Isochromatic Data

In the current paper a new method for calculating the stress intensity factor from isochromatic data is proposed. The method is based on the calculation of a map of the retardation using only the light- and dark-field images of a crack from a circular polariscope. The experimental retardation is compared with the theoretical one predicted by Westergaard’s model using between 700 and 800 data points in an array around the crack tip region. In this way, the necessity of unwrapping and calibrating the isochromatic fringe order distribution is avoided and a large potential source of error is eliminated. Subsequently, K_I, K_II and σ_ox are inferred by minimizing an error function. To show the potential and efficacy of the method, K values have been obtained from photoelastic images captured during fatigue tests conducted in a polycarbonate Center-Cracked-Tension (CCT) specimen. Results show an excellent level of agreement with those predicted from theory, highlighting the potential of the proposed methodology.     

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.     

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.     

Photoelastic determination of mixed-mode stress-intensity factorsK I,KII andK III

On the basis of existing photoelastic methods for the determination ofK _I andK _II, this paper presents an experimental method for determiningK _III with photoelastic data, and a photoelastic method for comprehensively determiningK _I,K _II andK _III under the complex stress condition. A frozen three-dimensional photoelastic model is first used to determineK _I andK _II from the slice perpendicular to the flaw edge. Then, from that slice, a sub-slice is taken to determine the factorK _III. This method is examined by comparison with two test models.     

Theoretical and experimental studies in fracture mechanics of crack-like defects under biaxial loading

It is a common point of view in fracture mechanics that, for any geometry of the body with a crack and any boundary conditions for the loading acting in the body plane, the stress and displacement components near the crack tip can be approximated in the framework of the theory of elasticity by a one-parameter or one-term representation, i.e., strictly in terms of the stress intensity coefficients K _I and K _II for an arbitrary failure crack [1, 2]. The authors of [2] specified the Westergaard function of the singular solution for a central crack under the biaxial loading of a plate. This approximate two-component solution has satisfactory accuracy. It is clear from [2] that this method cannot be admitted as a general statement [1], although it has long been assumed to be correct. The cause is that one cannot reasonably justify neglecting the second term in the Williams representation of the stress components in the plane case in the form of eigenfunction series; the contribution of this term in the rectangular coordinate system x, y is independent of the distance from the crack tip. This method may result in a serious mistake, from both the qualitative and quantitative viewpoints, in the prediction of local stresses, displacements, and related variables that are of interest. Apparently, this can best be demonstrated by an example of biaxial loading of a plate with a crack [1]. The unfounded neglect of the second term (whose contribution is independent of the distance from the crack tip) in the series representing the stress components is the source of the above-mentioned difficulties. In this problem, the influence of the load applied in the direction parallel to the crack plane manifests itself only in the second term of the series [3]. Therefore, this term should be clearly determined and studied in detail in the case of technological welding defects (faulty fusions, incomplete fusions, undercuts, and slag inclusions) and crack-like defects (scratches and cuts) in the base metal. The influence of the stress σ _OX along the crack axis on the stress tensor σ _x , σ _y , τ _xy and on the displacements u _x and u _y is confirmed by experimental studies of cracks by the photoelasticity method [4].     

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.     

A new criterion for mixed mode fracture initiation based on the crack tip plastic core region

In this work, we propose a new criterion for mixed mode I-II crack initiation angles based on the characteristics of the plastic core region surrounding the crack tip. The shape and size of the plastic core region are thoroughly analyzed under different loading conditions and a new formulation for the non-dimensional variable radius of the core region is presented for mixed mode (K_I−K_II) fracture. The proposed criterion states that the crack extends in the direction of the local or global minimum of the plastic core region boundary depending on the resultant stress state at the crack tip. The results show a well-defined correlation between the plastic core region characteristics and crack extension angles predicted by other criteria. The proposed criterion is formulated for various loading conditions and is compared with other available criteria against the limited available experimental data. It is shown that the proposed criterion provides a better agreement with the experimental data.     

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.     

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.     

On the relationship between critical tensile stress and fracture toughness in mild steel

An analysis is presented which relates the critical value of tensile stress (σ_f) for unstable cleavage fracture to the fracture toughness (K_Ic) for a high-nitrogen mild steel under plane strain conditions. The correlation is based on (i) the model for cleavage cracking developed by E. Smith and (ii) accurate plastic-elastic solutions for the stress distributions ahead of a sharp crack derived by J. R. Rice and co-workers. Unstable fracture is found to be consistent with the attainment of a stress intensification close to the tip such that the maximum principal stress σ_yy exceeds σ_f over a characteristic distance, determined as twice the grain size. The model is seen to predict the experimentally determined variation of K_Ic with temperature over the range -150 to -75°C from a knowledge of the yield stress and hardening properties. It is further shown that the onset of fibrous fracture ahead of the tip can be deduced from the position of the maximum achievable stress intensification. The relationship between the model for fracture ahead of a sharp crack, and that ahead of a rounded notch, is discussed in detail.     

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.     

A photoelastic study of crack repair

Birefringent coatings have been employed to study the effectiveness of an adhesively bonded repair of a center-cracked tension panel. The repair was one sided, with photoelastic coatings applied to the opposite side. Photoelastic coatings were also applied over the patch. Analysis methods are presented to permit the stress intensity factor to be determined from the isochromatic fringe patterns recorded from both continuous andX- andY-edged coatings. The results showed that the one-sided adhesively bonded patch reduced the stress intensity factor; however, the repair did not markedly change the character of stress distributions. Fringe loops formed near the crack tips for both the cracked and repaired tension panels. The primary difference was in the size of the loops. The reduction inK _I due to repair was smaller than anticipated, but even small improvements in ΔK _I markedly enhance the life of a repaired panel. The Paris power law is used to show the relation between the reduction in ΔK _I and the improvement in the crack growth rateda/dN. Fringe patterns from the birefringent coatings applied to the patch provided a means not only to investigate the stresses in the patch but also to detect the initiation of the local debonding of the adhesive in the neighborhood of the crack. The birefringent coating on the patch is an approach for producing an optically “smart” repair.     

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.     

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.     

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

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

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.     

A stress analytical solution for Mode III crack within modified gradient elasticity

The strain gradient exists near a crack tip may significantly influence the near-tip stress field. In this paper, the strain gradient and the internal length scales are introduced into the basic equations of mode III crack by the modified gradient elasticity (MGE). By using a complex function approach, the analytical solution of stress fields for mode III crack problem is derived within MGE. When the internal length scales vanish, the stress fields can be simplified to the stress fields of classical linear elastic fracture mechanics. The results show that the singularity of the shear stress is made up of two parts, r^−1/2 part and r^−3/2 part, and the sign of the stress σ_yz changes. With the increase of l_x, the peak value of σ_yz decrease and its location moves farther from the fracture vertex. The influence of strain gradient for mode III crack problem cannot be ignored.     

An analysis on three-dimensional effects near the tip of a crack in an elastic plate

An infinite plate containing a finite through crack under tensile loading is analysed by Fourier transform based on the Kane-Mindlin kinematic assumptions for the quasi-three-dimensional deformation of plates in extension. The asymptotic expressions of stress and displacement fields near the crack tip, the variation of the stress intensity factor with the plate-thickness and the three-dimensional deformation zone near the crack tip are investigated. The results of the analysis show that, (a) the crack-tip stress and displacement fields accounting for the plate-thickness effects are different from the plane stress solutions and this is true even for extremely small parameter (=1–vh/6 a). In a very small region near the crack tip, plane strain solutions prevail; (b) the ratio of the stress intensity factor K_I to the corresponding plane stress one K_I, K_I/K _I ^o , approaches 1/(1–v²) as tends to zero; (c) plane stress solutions can give satisfactory results for points a distance from the crack tip greater than about three-fourths of the plate-thickness; (d) the linear elastic result for the zone of three-dimensional effects is approximately valid for an elasto-plastic material with linear strain-hardening when the plastic tangential mudulus E_t is not very small.The Project Supported by National Natural Science Foundation of China.