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.     

Mixed Mode Dynamic Fracture in Particulate Reinforced Functionally Graded Materials

A detailed analytical and experimental investigation is presented to understand the dynamic fracture behavior of functionally graded materials (FGMs) under mode I and mixed mode loading conditions. Crack-tip stress, strain and displacement fields for a mixed mode crack propagating at an angle from the direction of property gradation were obtained through an asymptotic analysis coupled with a displacement potential approach. This was followed by a comprehensive series of experiments to gain further insight into the behavior of propagating cracks in FGMs. Dynamic photoelasticity coupled with high-speed photography was used to obtain crack tip velocities and dynamic stress fields around the propagating cracks. Birefringent coatings were used to conduct the photoelastic study due to the opaqueness of the FGMs. Dynamic fracture experiments were performed using different specimen geometries to develop a dynamic constitutive fracture relationship between the mode I dynamic stress intensity factor (K _ID ) and crack-tip velocity ( ) for FGMs with the crack moving in the direction of increasing fracture toughness. A similar -K _ID relation was also obtained for matrix material (polyester) for comparison purposes. The results obtained show that crack propagation velocities in FGMs were about 80% higher than the polyester matrix. Crack arrest toughness was found to be about 10% lower than the value of local fracture toughness in FGMs.     

A novel method in determination of dynamic fracture toughness under mixed mode I/II impact loading

The objective of this paper is to propose a novel methodology for determining dynamic fracture toughness (DFT) of materials under mixed mode I/II impact loading. Previous experimental investigations on mixed mode fracture have been largely limited to qusi-static conditions, due to difficulties in the generation of mixed mode dynamic loading and the precise control of mode mixity at crack tip, in absence of sophisticated experimental techniques. In this study, a hybrid experimental–numerical approach is employed to measure mixed mode DFT of 40Cr high strength steel, with the aid of the split Hopkinson tension bar (SHTB) apparatus and finite element analysis (FEA). A fixture device and a series of tensile specimens with an inclined center crack are designed for the tests to generate the components of mode I and mode II dynamic stress intensity factors (DSIF). Through the change of the crack inclination angle β (=90°, 60°, 45°, and 30°), the K_II/K_I ratio is successfully controlled in the range from 0 to 1.14. A mixed mode I/II dynamic fracture plane, which can also exhibit the information of crack inclination angle and loading rate at the same time, is obtained based on the experimental results. A safety zone is determined in this plane according to the characteristic line. Through observation of the fracture surfaces, different fracture mechanisms are found for pure mode I and mixed mode fractures.     

Stress intensity factor range and propagation mode of surface cracks under rolling–sliding contact

Finite element analyses were conducted in order to evaluate the mode I and mode II stress intensity factors for inclined edge cracks under cyclic contact load under rolling and rolling–sliding condition. The SIF range depends on crack orientation, crack length to Hertzian contact zone half-width ratio, friction between the crack faces and friction on the contact surface. The results were combined in two compact functions that determine the ΔK_I and ΔK_II values. The crack propagation mode and direction were investigated using both the maximum stress criterion and the minimum strain energy density criterion. The results are displayed in graph form, which allows a fast evaluation of the crack growth condition.     

平面应变下紧凑拉伸试样的动态断裂韧性的实验研究

材料的动态断裂韧性是衡量材料在动载荷作用下抵杭裂纹扩展能力的重要指标,以往的材料动态断裂韧性测试多采用三点弯曲试样,而针对紧凑拉伸试样的动态断裂韧性研究很少.本文将紧凑拉伸试样(即CT试样)简化成等效弹簧质量模型,得到了CT试样动态应力强度因子的近似表达式.对Hopkinson压杆装置进行了改进,利用改进后的实验装置进...     

Stress intensification in cracked shank of tightened bolt

Defects or cracks in the shank of bolts can degrade their load carrying capacity. The ways with which loading and residual stress intensify the crack border stress field can be reflected through the stress intensity factor quantity as defined in the linear elastic fracture mechanics theory. Use is made of the stiffness derivative method where quarter-point singular finite elements are used in the numerical calculation. Improved accuracy is achieved by considering the displacements not only of the main nodes but also of those quarter-point nodes in plane normal and adjacent to the crack.Numerical results are obtained for a semi-elliptical shaped crack in the bolt shank owing to tension, bending, residual stress and stress caused by tightening of the bolt. Maximum value of the Mode I stress intensity factor K_a due to tension or bending could prevail either at the deepest point on the crack border or at the root of the shank where the crack border terminates depending on the aspect ratio of the ellipse. In general, K₁ at the deepest point of crack penetration is larger than that at the free surface for tension and bending for a fixed crack depth with reference to the bolt diameter. Tightening of the bolt tends to increase K₁ at the free surface if the crack depth is small. The opposite is obtained for deeper cracks. Assumed residual stress effect obtained from experimental data is found to have negligible influence on the stress intensity factor when compared with that arising from tensile load.     

Dynamic stress intensity factor K III and dynamic crack propagation characteristics of anisotropic materials

Based on the mechanics of anisotropic materials, the dynamic propagation problem of a mode Ⅲ crack in an infinite anisotropic body is investigated. Stress, strain and displacement around the crack tip are expressed as an analytical complex function, which can be represented in power series. Constant coefficients of series are determined by boundary conditions. Expressions of dynamic stress intensity factors for a mode Ⅲ crack are obtained. Components of dynamic stress, dynamic strain and dynamic displacement around the crack tip are derived. Crack propagation characteristics are represented by the mechanical properties of the anisotropic materials, i.e., crack propagation velocity M and the parameter ~. The faster the crack velocity is, the greater the maximums of stress components and dynamic displacement components around the crack tip are. In particular, the parameter α affects stress and dynamic displacement around the crack tip.     

Effect of temperature on the fracture modes of BHW-35 steel under mode-II loading

Finite element method (FEM) has been used to analyze the stress and strain fields and the stress tri-axial levels around the tip of the crack under mode- II loading. The results show that: under mode- II loading, the direction of the maximum tensile stress and that of the maximum tri-axial levels (R _o ) exist at an angle of –75. 3° from the original crack plane; the maximum shear stress andR _o = 0 exist along the original crack plane.Mode- II loading experiment using BHW-35 steel at different temperatures show that there are two kinds of fracture mode, opening mode (or tensile mode) and sliding mode (or shear mode). A decrease in temperature causes the fracture mode to change from shear mode to tensile mode. For BHW-35 steel, this critical temperature is about –90 C. Actually, under any kind of loading mode (mode I . mode II , mode III or mixed mode), there always exist several kinds of potenital fracture modes (for example, opening mode, sliding mode, tearing mode or mixed mode). The effect of temperature under mode- II loading is actually related to the change of the elastic-plastic properties of the material.     

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.     

Brittle fracture: Variation of fracture toughness with constraint and crack curving under mode I conditions

The effect of constraint on brittle fracture of solids under predominantly elastic deformation and mode I loading conditions is studied. Using different cracked specimen geometry, the variation of constraint is achieved in this work. Fracture tests of polymethyl methacrylate were performed using single edge notch, compact tension and double cantilever beam specimens to cover a bread range of constraint. The test data demonstrate that the apparent fracture toughness of the material varies with the specimen geometry or the constraint level. Theory is developed using the critical stress (strain) as the fracture criterion to show that this variation can be interpreted using the critical stress intensity factorK _Cand a second parameterT orA ₃,whereT andA ₃are the amplitudes of the second and the third term in the Williams series solution, respectively. The implication of this constraint effect to the ASTM fracture toughness value, crack tip opening displacement fracture criterion and energy release rateG _Cis discussed. Using the same critical stress (strain) as the fracture criterion, the theory further predicts crack curving or instability under mode I loading conditions. Experimental data are presented and compared with the theory.     

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.     

Unsteady Dynamic Crack Propagation in a Brittle Polymer

Dynamic crack propagation in a brittle polymer, poly(methyl-methacrylate) (PMMA), was studied using the method of caustics in combination with a Cranz–Schardin high-speed camera. Four different types of specimen geometry and loading method were employed to achieve the crack acceleration, deceleration, and/or reacceleration processes in one fracture event. The dynamic stress intensity factor K _ID and crack velocity were obtained in the course of the crack propagation and the corresponding relationship was determined. The effect of the crack acceleration and deceleration on the K _ID-velocity relationships was as follows: (1) the variations of K _ID and the velocity were strongly influenced by the specimen geometry and loading method; (2) the velocity change was qualitatively in accord with K _ID; (3) K _ID for a constant crack velocity was larger when the crack decelerated than it was when the crack accelerated or reaccelerated; (4) K _ID for an acceleration-free crack was uniquely related to the velocity; and (5) K _ID could be expressed as two parametric functions of the velocity and acceleration.     

Dynamic crack tip fields and dynamic crack propagation characteristics of anisotropic material

Based on mechanics of anisotropic material, the dynamic crack propagation problem of I/II mixed mode crack in an infinite anisotropic body is investigated. Expressions of dynamic stress intensity factors for modes I and II crack are obtained. Components of dynamic stress and dynamic displacements around the crack tip are derived. The strain energy density theory is used to predict the dynamic crack extension angle. The critical strain energy density is determined by the strength parameters of anisotropic materials. The obtained dynamic crack tip fields are unified and applicable to the analysis of the crack tip fields of anisotropic material, orthotropic material and isotropic material under dynamic or static load. The obtained results show Crack propagation characteristics are represented by the mechanical properties of anisotropic material, i.e., crack propagation velocity M and fiber direction α. In particular, the fiber direction α and the crack propagation velocity M give greater influence on the variations of the stress fields and displacement fields. Fracture angle is found to depend not only on the crack propagation but also on the anisotropic character of the material.     

Dynamic mode I perturbation solution for a moving crack unsteadily

Rice et al. (Journal of Mechanics and Physics of Solids42, 813–843) analyze the propagation of a planar crack with a nominally straight front in a model elastic solid with a single displacement component. Using the form of Willis et al. (Journal of the Mechanics and Physics of Solids43, 319–341), of dynamic mode I weight functions for a moving crack, we address that problem solved by Rice et al. in the 3D context of elastodynamic theory. Oscillatory crack tip motion results from constructive-destructive interference of stress intensity waves. Those waves, including system of the dilatational, shear and Rayleigh waves, interact on each other and with moving edge of crack, can lead to continuing fluctuations of the crack front and propagation velocity.     

General solutions of plane problem for power function curved cracks

A new exact and universal conformal mapping is proposed. Using Muskhelishvili's complex potential method, the plane elasticity problem of power function curved cracks is investigated with an arbitrary power of a natural number, and the general solutions of the stress intensity factors (SIFs) for mode I and mode II at the crack tip are obtained under the remotely uniform tensile loads. The present results can be reduced to the well-known solutions when the power of the function takes different natural numbers. Numerical examples are conducted to reveal the effects of the coefficient, the power, and the projected length along the x-axis of the power function curved crack on the SIFs for mode I and mode II.     

Analytical and experimental determination of dynamic impact stress intensity factor for 40 Cr steel

A linear elastic three-dimensional finite element analysis is made to analyze the near field stress behavior of an edge cracked rectangular bar simply supported and subjected to central impact at the back side of the crack. The material is made of 40 Cr steel. Determined numerically are the local time histories of the stress wave, displacement near load point, crack tip strain, and dynamic stress intensity factor K_(d)¹. The above quantities were also measured experimentally by performing impact tests; they agreed well with the analytical results and determine the load at fracture initiation and hence the critical dynamic stress intensity factor K_(d)^1c. The interaction effect between the loading bar and specimen appears to be negligible.     

An improved experimental method for the determination of the critical stress intensity factor KIc of brittle materials

The critical stress intensity factor K_Ic is determined by a simple and accurate method, using small test specimens and a simple procedure in this paper.Single edge V-notched tension specimens made of PMMA are subjected to a load which is slowly increased until the crack begins to move from the notch tip. During the crack propagation event shadow patterns at the tip of the crack are recorded in a video recorder. Under these loading conditions, the creating real crack propagate slowly until the crack propagation velocity take an abrupt increase and the entire fracture of the specimen takes place. The stress intensity factor which correspond to the transition from the slow to fast crack speed, is the critical stress intensity factor K_Ic and it can be the fracture toughness of the material.The results are accurate and in good agreement with those values of K_Ic which are calculated by approximate theoretical expressions.The purpose of this paper is to introduce an improved, simple and accurate experimental method for the determination of fracture toughness of brittle materials.     

Finite element analysis of mode III steady crack growth in anisotropic hardening materials

In this paper, the steady crack growth of mode III under small scale yielding conditions is investigated for anisotropic hardening materials by the finite element method. The elastic-plastic stiffness matrix for anisotropic materials is given. The results show the significant influences of anisotropic hardening behaviour on the shape and size of plastic zone and deformation field near the crack tip. With a COD fracture criterion, the ratio of stress intensity factorsk _ss/k_c varies appreciably with the anisotropic hardening parameterM and the hardening exponentN.     

Fracture dynamics problems of mode I semi-infinite crack for anisotropic orthotropic body

By means of the theory of complex functions, fracture dynamics problems of mode I semi- infinite crack for anisotropic orthotropic body were researched. Analytical solutions of stress, displacement, and dynamic stress intensity factor under the action of moving increasing loads Px ³/t ³, Pt ⁴/x ³, respectively, are very easily obtained utilizing the approaches of self-similar functions. In the light of relevant material’s coefficients, the alterable rule of dynamic stress intensity factor was depicted very well. The correlative closed solutions are attained based on the Riemann–Hilbert problems. After those analytical solutions were applied by the superposition principle, the solutions of discretional complex problems could be attained.