Fatigue crack growth in high and low strength steel under torsional loading

During loading of a crack in mode III the crack surfaces in contact slide against each other giving rise to friction, abrasion and mutual support, thereby reducing the effective stress at the crack tip (“sliding mode crack closure”). This phenomenon was investigated in a high strength steel (AISI 4340) and in a low strength steel (AISI C1018) in circumferentially notched specimens under pure cyclic torsion and combined loading (cyclic torsion plus static axial load). The influence of sliding mode crack closure on fatigue crack propagation is shown and “true” crack growth values (without the sliding mode crack closure influence) are determined on the basis of an extrapolation procedure. Explanations are given for causes of the various fracture modes observed, such as “factory roof” fracture, macroscopically flat mode III fracture and “lamella” fracture. Finally the scientific and technical importance of sliding mode crack closure is demonstrated.     

Mechanics and fracture of crack tip deformable bi-material interface

Novel interface deformable bi-layer beam theory is developed to account for local effects at crack tip of bi-material interface by modeling a bi-layer composite beam as two separate shear deformable sub-layers with consideration of crack tip deformation. Unlike the sub-layer model in the literature in which the crack tip deformations under the interface peel and shear stresses are ignored and thus a “rigid” joint is used, the present study introduces two interface compliances to account for the effect of interface stresses on the crack tip deformation which is referred to as the elastic foundation effect; thus a flexible condition along the interface is considered. Closed-form solutions of resultant forces, deformations, and interface stresses are obtained for each sub-layer in the bi-layer beam, of which the local effects at the crack tip are demonstrated. In this study, an elastic deformable crack tip model is presented for the first time which can improve the split beam solution. The present model is in excellent agreements with analytical 2-D continuum solutions and finite element analyses. The resulting crack tip rotation is then used to calculate the energy release rate (ERR) and stress intensity factor (SIF) of interface fracture in bi-layer materials. Explicit closed-form solutions for ERR and SIF are obtained for which both the transverse shear and crack tip deformation effects are accounted. Compared to the full continuum elasticity analysis, such as finite element analysis, the present solutions are much explicit, more applicable, while comparable in accuracy. Further, the concept of deformable crack tip model can be applied to other bi-layer beam analyses (e.g., delamination buckling and vibration, etc.).     

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.     

Crack kinking from an interface crack with initial contact between the crack lips

The authors recently theoretically studied crack kinking and opening from an initially closed crack (without friction) in some homogeneous medium. The same problem, but for an interface crack, is considered here. Comninou has shown that the asymptotic stress field prior to kinking is governed by a single, mode II stress intensity factor (SIF). Using this result, plus a homogeneity property of the problems of elastic fracture mechanics with unilateral contact envisaged, a change of scale, and two reasonable hypotheses, we establish the expression of the SIF at the tip of the small, open crack extension. It is shown that whatever the geometry of the external boundary and the crack and whatever the loading, these SIF depend solely upon the initial (mode II) SIF (in a linear way), the kink angle and Dundurs' parameters α and β. Using this result and Goldstein and Salganik's “principle of local symmetry” to predict the kink angle, one finds that it is independent of the loading but does depend upon Dundurs' parameters α and β. This contrasts with the case of an ordinary (initially closed) crack in some homogeneous medium, for which the kink angle was not only independent of the loading but an absolute constant. However, it is numerically found that the influence of the mismatch of elastic properties upon the kink angle is rather weak.     

Gradient elasticity with surface energy: mode-III crack problem

In this paper an anisotropic strain-gradient dependent theory of elasticity is exploited, which contains both volumetric and surface energy gradient dependent terms. The theory is applied to the solution of the mode-III crack problem and is extending previous results by Aifantis and co-workers. The two boundary value problems corresponding to the “unclamped” and “clamped” crack tips, respectively, are solved analytically. It turns out that the first problem is physically questionable for some values of the surface energy parameter, whereas the second boundary value problem is leading to a cusping crack, which is consistent with Barenblatt's theory without the incorporation of artificial assumptions.     

Energy approach to crack growth in a fiber-reinforced brittle material modeled by an inclusion

Brittle materials randomly reinforced with a low volume fraction of strong, stiff and ductile fibers are considered, with specific reference to fiber-reinforced cements and concrete. Visible cracks in such materials are accompanied by a surrounding damage zone – together these constitute a very complex “crack system”. Enormous effort has been put into trying to understand the micromechanics of such systems. Almost all of these efforts do not deal with the “crack system” propagation behavior as a whole. The propagation process of such a “crack system” includes propagation of the visible crack and the growth of the damage zone. Propagation may take place by lengthening of the visible crack together with the concomitant lengthening of the surrounding damage zone, or simply by broadening of the damage zone while the visible crack length remains unchanged – or simultaneously by growth of both types. A phenomenological completely theoretical model (for an ideal material) is here proposed which can serve to examine the propagation process by means of energy principles, without recourse to the microscopic details of the process. An application of this theoretical approach is presented for the case of a damage zone evolving with a rectangular shape. This shape is chosen because it is expected that it will illustrate the nature of damage evolution and because the computational procedure necessary to follow the growth is the most straightforward.     

Crack propagation in wedged double cantilevered beam specimens

A closed form analytical solution of crack propagation in double cantilevered beam specimens opened at a constant rate has been found. Hamilton's principle for non-conservative systems was applied to describe the crack motion, under the assumption of a Bernoulli-Euler beam. The criterion of crack propagation is a critical bending moment at the crack tip. The calculations of beam motion take into account wave effects in the Bernoulli-Euler theory of elastic beams. The beam shape during the crack motion is found with a similarity transformation and expressed by Fresnel integrals. The boundary conditions satisfied are the fixed ones of zero bending moment and constant beam opening rate at the load end of the specimen and the moving ones of zero deflection and zero slope of the deflected beam at the tip of the moving crack. The fracture represents a moving critical bending moment. The analytical results show that the specific fracture surface energy is a unique function of the ratio of the crack length squared to the time subsequent to loading and this is computed from the recorded time-dependence of the crack length.     

Two general integrals of singular crack tip deformation fields

The Eshelby tensor E has vanishing divergence in a homogeneous elastic material, whereas the invariance of the crack tip J integral suggests, in accord with known solutions, that the product rE will have a finite limit at the tip. Here r is distance from the tip. These considerations are shown to lead to two general integrals of the equations governing singular crack tip deformation fields. Some of their consequences are discussed for analysis of crack tip fields in linear and nonlinear materials.     

Interaction integrals for fracture analysis of functionally graded piezoelectric materials

This paper presents domain form of the interaction integrals based on three independent formulations for computation of the stress intensity factors and electric displacement intensity factor for cracks in functionally graded piezoelectric materials. Conservation integrals of J-type are derived based on the governing equations for piezoelectric media and the crack tip asymptotic fields of homogeneous piezoelectric medium as auxiliary fields. Each of the formulation differs in the way auxiliary fields are imposed in the evaluation of interaction integral and each of them results in a consistent form of the interaction integral in the sense that extra terms naturally appears in their derivation to compensate for the difference in the chosen crack tip asymptotic fields of homogeneous and functionally graded piezoelectric medium. The additional terms play an important role of ensuring domain independence of the presented interaction integrals. Comparison of the numerically evaluated intensity factors through the three consistent formulations with those obtained using displacement extrapolation method is presented by means of two examples.     

J-integral and crack driving force in elastic–plastic materials

This paper discusses the crack driving force in elastic–plastic materials, with particular emphasis on incremental plasticity. Using the configurational forces approach we identify a “plasticity influence term” that describes crack tip shielding or anti-shielding due to plastic deformation in the body. Standard constitutive models for finite strain as well as small strain incremental plasticity are used to obtain explicit expressions for the plasticity influence term in a two-dimensional setting. The total dissipation in the body is related to the near-tip and far-field J-integrals and the plasticity influence term. In the special case of deformation plasticity the plasticity influence term vanishes identically whereas for rigid plasticity and elastic-ideal plasticity the crack driving force vanishes. For steady state crack growth in incremental elastic–plastic materials, the plasticity influence term is equal to the negative of the plastic work per unit crack extension and the total dissipation in the body due to crack propagation and plastic deformation is determined by the far-field J-integral. For non-steady state crack growth, the plasticity influence term can be evaluated by post-processing after a conventional finite element stress analysis. Theory and computations are applied to a stationary crack in a C(T)-specimen to examine the effects of contained, uncontained and general yielding. A novel method is proposed for evaluating J-integrals under incremental plasticity conditions through the configurational body force. The incremental plasticity near-tip and far-field J-integrals are compared to conventional deformational plasticity and experimental J-integrals.     

A generalized force measure of conditions at crack tips

A tentative measure of the forces tending to cause crack growth—the apparent crack extension force—is defined within the framework of continuum mechanics. By an associated fracture criterion initiation of growth may be predicted as well as the direction of preferred growth. The theory is specialized to elastic, viscoelastic and elastic-plastic materials. Under specified conditions the apparent crack extension force may be expressed by surface integrals over the boundary of an arbitrary part of the body for quasi-static deformation and for steady-state propagation of the crack. For plane deformation and for infinitesimal deformation under plane stress conditions these surface integrals reduce to path independent line integrals which include the J integral by Rice[1] and the G integral by Sih[2] as special cases.     

Calculation for path-domain independentJ integral with elasto-viscoplastic consistent tangent operator concept-based boundary element methods

This paper presents an elasto-viscoplastic consistent tangent operator (CTO) based boundary element formulation, and application for calculation of path-domain independentJ integrals (extension of the classicalJ integrals) in nonlinear crack analysis. When viscoplastic deformation happens, the effective stresses around the crack tip in the nonlinear region is allowed to exceed the loading surface, and the pure plastic theory is not suitable for this situation. The concept of consistency employed in the solution of increment viscoplastic problem, plays a crucial role in preserving the quadratic rate asymptotic convergence of iteractive schemes based on Newton's method. Therefore, this paper investigates the viscoplastic crack problem, and presents an implicit viscoplastic algorithm using the CTO concept in a boundary element framework for path-domain independentJ integrals. Applications are presented with two numerical examples for viscoplastic crack problems andJ integrals. The project supported by National Natural Science Foundation of China (9713008) and Zhejiang Natural Science Foundation Special Funds No. RC.9601     

Cracking characteristics of mixed mode dislocations near a lip-like mode crack

This work is concerned with the cracking characteristics of mixed mode dislocations near a lip-like mode crack, stress intensity and strain energy density factor are obtained by using conformal mapping, singularity analysis and Cauchy integrals. Shielding effect generated by screw dislocation near a lip-like mode crack decreases with the increment of the distance between screw dislocation and crack tip. Larger distance between two faces of the crack leads to the shielding effect waning. The strain energy density factor of mode III decreases with the increment of the distance between dislocation and crack tip. Larger distance between two faces of lip-like mode crack also leads to the strain energy density factor waning and encourages crack initiation; the shielding effects generated by edge dislocation near the crack decrease with the increment of the distance between edge dislocation and crack tip.     

Boundary element weight function analysis of a strip yield crack in a rotating disk

The boundary element method combined with subtration of Bueckner singular fields are used to obtain weight functions for an internal edge crack in a rotating annular disk. A previously developed, general representation of the weight function is used which leads to integrals that can be evaluated analytically to obtain the stress intensity factor and surface displacements of the crack. The determination of crack tip opening displacements for the strip yield crack is reduced to a non-singular integral which can be evaluated in closed form. The strip yield zone length and crack tip opening displacement are determined for an internal radial crack in a rotating annular disk for a range of crack lengths and rotational speeds.     

Logarithm strain singularity at tip of mode I growing crack in elasticand perfectly-plastic material

Reanalyzed in detail is the stress and strain distribution near the tip of a Mode I steadily growing crack in an elastic and perfectly-plastic material. The crack tip region is divided into five angular sectors, one of which is singular in character and represents a rapid transition zone that becomes a line of strain discontinuity in the limit as crack tip is approached. It is shown for an incompressible material (ν=0.5) under plane strain that the local strain in all the angular sectors possesses the same logarithm singularity, i.e., In r where r is the radial distance measured from the crack tip. This result also prevails for the compressible material ( v < 0.5) and resolves a long standing controversy concerning the strain singularity in the sector just ahead of the crack tip.     

SELF-SIMILAR CRACK EXPANSION METHOD FOR THE STRESS INTENSITY FACTOR ANALYSIS

The Self-Similar Crack Expansion (SSCE) method is proposed to evaluate stress intensi-ty factors at crack tips, whereby stress intensity factors of a crack can be determined by the crackopening displacement over the crack, not just by the local displacement around the crack tip. The crackexpansion rate is estimated by taking advantage of the crack self-similarity. Therefore, the accuracy ofthe calculation is improved. The singular integrals on crack tip elements are also analyzed and are pre-cisely evaluated in terms of a special integral analysis. Combination of these two techniques greatly in-creases the accuracy in estimating the stress distribution around the crack tip. A variety of two-dimen-sional cracks, such as subsurface cracks, edge cracks, and their interactions are calculated in terms ofthe self-similar expansion rate. Solutions are satisfied with errors less than 0.5% as compared with theanalytical solutions. Based on the calculations of the crack interactions, a theory for crack interactionsis proposed such that for a group of aligned cracks the summation of the square of SIFs at the right tipsof cracks is always equal to that at the left tips of cracks. This theory was proved by the mehtod ofSelf-Similar Crack Expansion in this paper.     

Plane-strain crack tip field and a dual-parameter fracture-criterion for non-linear fracture mechanics

Hancock and Cowling measured the critical crack tip opening displacements, δ_f, at fracture initiation in HY-80 steel specimens of six different configurations. δ_f varied from 90 μm in a deeply double-edge-cracked tensile panel to 900 μm in a single-edge-cracked tensile panel.McMeeking and Parks, and Shih and German have shown by their finite element calculations that the characteristics of the plane strain crack tip fields in both large scale yielding and general yielding are strongly dependent on specimen geometry and load level.In this study, the plane strain crack tip fields in the specimens tested by Hancock and Cowling were calculated using the finite element method. The crack tip triaxial tensile stress field is strongly affected by specimen geometric constraint, and the state of the triaxial tensile stress in a crack tip region is monitored by the ratio between the local tensile stress and the effective stress, i.e., ( ), at a distance x=2δ from the crack tip. The values of ( ) vary from 3.1 for the double-edge-cracked tensile panel to 1.7 for the single-edge-cracked tensile panel. The δ_f measured by Hancock and Cowling correlates very well with the ratio ( ). δ_f is a measure of the fracture ductility of the material ahead of the crack tip, and the ductility decreases with an increase in the triaxial tensile stress, i.e., the ratio ( ).     

A steady-state moving crack of longitudinal shear with an infinitely narrow plastic zone

A simple model of longitudinal shear with a plastic zone at the tip of the crack is considered in the paper. The plasticity zone is assumed to be narrow, and fulfillment of the Mises condition is stipulated on its boundary. The crack moves with a constant velocity without change in the length.Translated from Zhurnal Prikladnoi Mekhaniki i Teknieheskoi Fiziki, No. 2, pp. 121–126, March–April, 1973.     

Displacement discontinuity method for modeling axisymmetric cracks in an elastic half-space

This paper describes a displacement discontinuity method for modeling axisymmetric cracks in an elastic half-space or full space. The formulation is based on hypersingular integral equations that relate displacement jumps and tractions along the crack. The integral kernels, which represent stress influence functions for ring dislocation dipoles, are derived from available axisymmetric dislocation solutions. The crack is discretized into constant-strength displacement discontinuity elements, where each element represents a slice of a cone. The influence integrals are evaluated using a combination of numerical integration and a recursive procedure that allows for explicit integration of hyper- and Cauchy singularities. The accuracy of the solution at the crack tip is ensured by adding corrective stresses across the tip element. The method is validated by a comparison with analytical and numerical reference solutions.     

Retardation of a crack with connections between the faces using an induced thermoelastic stress field

This paper considers local temperature variations near the tip of a crack in the presence of regions in which the crack faces interact. It is assumed that these regions are adjacent to the crack tip and are comparable in size to the crack size. The problem of local temperature variations consists of delay or retardation of crack growth. For a crack with connections between the crack faces subjected to external tensile loads, an induced thermoelastic stress field, and the stresses at the connections preventing crack opening, the boundary-value problem of the equilibrium of the crack reduces to a system of nonlinear singular integrodifferential equations with a Cauchy kernel. The normal and tangential stresses at the connections are found by solving this system of equations. The stress intensity factors are calculated. The energy characteristics of cracks with tip regions are considered. The limiting equilibrium condition for cracks with tip regions is formulated using the criterion of limiting stretching of the connections.Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 1, pp. 133–143, January–February, 2005