FRACTURE ANALYSIS OF A FUNCTIONALLY GRADED STRIP UNDER PLANE DEFORMATION

In this paper the plane elasticity problem for a functionally graded strip containing a crack is considered. It is assumed that the reciprocal of the shear modulus is a linear function of the thickness-coordinate, while the Possion's ratio keeps constant. By utilizing the Fourier transformation technique and the transfer matrix method, the mixed boundary problem is reduced to a system of singular integral equations that are solved numerically. The influences of the geometric parameters and the graded parameter on the stress intensity factors and the strain energy release rate are investigated. The numerical results show that the graded parameters, the thickness of the strip and the crack size have significant effects on the stress intensity factors and the strain energy release rate.     

断裂力学判据的评述

从Inglis 和Griffith 的著名论文到Irwin 和Rice 等的奠基性贡献,对断裂力学中的线弹性断裂力学的K判据,界面断裂力学的G判据,和弹塑性断裂力学的J 判据作了扼要的综述. 介绍了在界面断裂力学G判据的基础上提出的界面断裂力学的K判据,以说明断裂力学的判据存在改进的可能性. 在综述中归纳出断裂力学判据中目前还没有较好解决的几个问题. 在总结以往断裂力学研究经验的基础上,指出裂纹端应力奇异性的源是对断裂力学判据存在的问题作进一步研究的切入点. 探讨了裂纹端应变间断的奇点是裂纹端应力奇异性的源的问题,从而对裂纹端应力强度因子的物理意义进行了讨论. 最后,阐述了进行可靠的裂纹端应力场的弹塑性分析是改进弹塑性断裂力学判据的关键,而进行可靠的裂纹端应力场的弹塑性分析的前提是要通过裂纹端应力奇异性的源的研究来获得作用在裂纹端的造成裂纹端应变间断的有限值应力.      

A MOVING CRACK IN A NONHOMOGENEOUS MATERIAL STRIP

This paper considers an anti-plane moving crack in a nonhomogeneous material strip of finite thickness. The shear modulus and the mass density of the strip are considered for a class of functional forms for which the equilibrium equation has analytical solutions. The problem is solved by means of the singular integral equation technique. The stress field near the crack tip is obtained. The results are plotted to show the effect of the material non-homogeneity and crack moving velocity on the crack tip field. Crack bifurcation behaviour is also discussed. The paper points out that use of an appropriate fracture criterion is essential for studying the stability of a moving crack in nonhomogeneous materials. The prediction whether the unstable crack growth will be enhanced or retarded is strongly dependent on the type of the fracture criterion used. Based on the analysis, it seems that the maximum 'anti-plane shear' stress around the crack tip is a suitable failure criterion for moving cracks in nonhomogeneous materials.     

Pre-kinking of a moving crack in a magnetoelectroelastic material under in-plane loading

A constant moving crack in a magnetoelectroelastic material under in-plane mechanical, electric and magnetic loading is studied for impermeable crack surface boundary conditions. Fourier transform is employed to reduce the mixed boundary value problem of the crack to dual integral equations, which are solved exactly. Steady-state asymptotic fields near the crack tip are obtained in closed form and the corresponding field intensity factors are expressed explicitly. The crack speed influences the singular field distribution around the crack tip and the effects of electric and magnetic loading on the crack tip fields are discussed. The crack kinking phenomena is investigated using the maximum hoop stress intensity factor criterion. The magnitude of the maximum hoop stress intensity factor tends to increase as the crack speed increases.     

Mixed mode fatigue crack growth in test coupons made from 2024-T3 aluminum

In this paper, two different fracture criteria are applied to determine the crack trajectory or angle of crack propagation in test specimens containing inclined cracks emanating from open holes. Also, different crack growth rate models are assumed for each criterion. The maximum principal stress criterion is used with a crack growth-rate equation based on an effective stress intensity factor. The strain energy density criterion is used with a crack growth-rate equation corresponding to an effective strain energy density factor. The crack growth-rate models for each criterion were constructed using unpublished fatigue crack growth data for 2024-T3 aluminum.     

A Mode- II Griffith Crack in Decagonal Quasicrystals

By using the method of stress functions, the problem of mode-II Griffith crack in decagonal quasicrystals was solved. First, the crack problem of two-dimensional quasicrystals was decomposed into a plane strain state problem superposed on anti-plane state problem and secondly, by introducing stress functions, the18 basic elasticity equations on coupling phonon-phason field of decagonal quasicrystals were reduced to a single higher-order partial differential equations. The solution of this equation under mixed boundary conditions of mode-II Griffith crack was obtained in terms of Fourier transform and dual integral equations methods. All components of stresses and displacements can be expressed by elemental functions and the stress intensity factor and the strain energy release rate were determined. Biography: GUO Yu-cui (1962-), Associate professor, Doctor     

Mode-I crack projection procedure using the strain energy density criterion

A practice used in linear elastic fracture mechanics is the projection of a crack onto a plane normal to the principal tensile stress axes for computing the stress intensity factor K_I. The minimum strain-energy criterion is applied for different crack configurations with the introduction of a safety factor S_i which is the ratio of the strain energy density factor of the projected crack and that of the original crack. Numerous crack configurations are investigated to illustrate the degree of conservativeness of the crack projection procedure.     

A long circular cylinder with a circumferential edge crack subjected to a uniform shearing stress

This paper contains an analysis of the stress distribution in a long circular cylinder of isotropic elastic material with a circumferential edge crack when it is deformed by the application of a uniform shearing stress. The crack with its center on the axis of the cylinder lies on the plane perpendicular to that axis, and the cylindrical surface is stress-free. By making a suitable representation of the stress function for the problem, the problem is reduced to the solution of a pair of singular integral equations. This pair of singular integral equations is solved numerically, and the stress intensity factor due to the effect of the crack size is tabulated.     

Mixed mode axisymmetric annular cracks in a finite layer: Off angle crack initiation

The solutions of axisymmetric Volterra type climb and glide edge dislocations are obtained in a layer by means of the Hankel transforms. Utilizing the same procedure, Green’s function solution is obtained for a layer under self-equilibration normal ring traction. The distributed dislocation technique is used to construct integral equations for a system of co-axial annular cracks where the layer is under axisymmetric normal loads. These equations are solved numerically to obtain dislocation density on the cracks surfaces. The results are employed to determine stress intensity factors for annular and penny-shaped cracks and the interaction between two co-axial penny-shaped cracks is studied. Moreover, the stress intensity factors of the interacting cracks are determined such that they can be further used in conjunction with strain energy density (SED) failure criterion to obtain the possible direction of crack initiation that may not be apparent under mixed mode conditions.     

A Zener–Stroh crack interacting with a coated inclusion with generalized Irwin plastic zone correction

The elastic–plastic fracture behavior of a Zener–Stroh crack interacting with a coated inclusion in composite materials has been investigated with crack tip plastic zone corrections. With the distributed dislocation method, the crack problem is formulated into a set of singular integral equations which are solved numerically. The plastic zone sizes at the both crack tips are determined by a generalized Irwin model where Von Mises stress yielding criterion is used. The stress intensity factor (SIF), the plastic zone size (PZS), the crack tip opening displacement (CTOD) and the effective stress intensity factor have been evaluated. In the numerical examples, the influence of the inclusion shear modulus, the coating-layer thickness and shear modulus, as well as the distance between the crack and inclusion, on the SIF, the PZS and the CTOD are discussed in detail. Numerical examples show that increasing the shear modulus or the thickness of the coating phase, the influence of the inclusion on the normalized SIF and the normalized PZS will be shielded.     

Thermally induced fracture of a smart functionally graded composite structure

Discussed is the fracture behavior of a cracked smart actuator on a substrate under thermal load. The actuator is made of piezoelectric material with functionally graded material (FGM) properties. Integral transform method is used to reduce the problem to the solution of a set of singular integral equations and is solved numerically. This paper is completed by including graphical plots of the thermal flow, stress and electric displacement intensity factors around the crack for different crack positions and material gradients. Directions of crack initiation are also predicted by using the energy density criterion.     

Stable crack growth along a brittle\ductile interface—II. Small scale yielding solutions and interfacial toughness predictions

The problem of a crack growing steadily and quasi-statically along a brittle\ductile interface under plane strain, mixed mode, and small scale yielding conditions is considered. The ductile material is assumed to be characterized by the J₂-flow theory of plasticity with linear strain hardening, while the brittle material is assumed to be linear elastic. A displacement-based finite element method, exploiting the convective nature of the problem, is utilized to solve the relevant boundary value problem. In Part I of this work, the corresponding asymptotic problem was solved. This paper addresses the full-field problem in order to validate the asymptotic solutions, and to explore the physical implications of the results. The numerical full-field results are found to be in good agreement with the analytical asymptotic solutions. In particular, the full-field results strongly suggest that the stress fields in the vicinity of the crack tip are variable-separable of the power singular type; and also that the mode mix of the near-tip stress fields is, to a large extent, independent of the applied elastic mode mix. The amplitude (the plastic stress intensity factor) and the regions of validity of the asymptotic fields are estimated from the full-field results, and are observed to be strongly dependent on the applied mode mix. The remote elastic loading fields appear to influence the near-tip fields, primarily, through the plastic stress intensity factor. The present work also explores the suggestion made by Bose and Ponte Castaneda, 1992 that the solutions to the small scale yielding problem may be used in the context of a standard crack growth criterion, requiring that continued growth take place with a fixed near-tip crack opening profile, to obtain theoretical predictions for the dependence of interfacial toughness on the applied mode mix. Based on the numerical results, predictions for mixed mode toughness of the brittle\ductile interface are reported. The results, which are in qualitative agreement with available experimental data and also with some recent theoretical results, predict a strong dependence of interfacial toughness on mode mix. This suggests that ductility provides the main operating mechanism for explaining the dependence of interfacial toughness on the mode mix of the applied loading fields, during steady crack growth.     

弹性T项对裂尖参数和裂纹扩展路径稳定性的作用

研究了弹性T项在主裂纹与近裂尖微空洞干涉问题中对裂尖参数和裂纹扩展路径稳定性的影响作用.利用“伪力”法,并通过解决主裂纹与近裂尖微空洞干涉问题,对远场纯1型载荷和远场弹性T项载荷下,该干涉问题中弹性T项的影响作用进行了分析从数值结果可以看出：由于空洞的存在;释放了弹性T项所引起的应力,弹性T项对裂尖参数;应力强度因子和J积分都有直接显著的影响,因而,它对该载荷下的裂纹扩展路径的稳定性有控制作用。     

Development of a multiaxial fatigue theory by considering constraint effects on small mixed-mode cracks

The effects of the transverse strain (the normal strain in the crack-line direction) on the near-tip fields of small shallow surface cracks (Case A cracks) in power-law hardening materials are investigated by finite element analyses. The small Case A cracks are under plane stress, general yielding, and mixed mode I and II conditions. Constant effective stress contours representing the intense straining zones near the tip, deformed crack-tip profiles and near-tip mode mixity factors are presented for different transverse strains in the crack-line direction. Based on the concept of characterization of fatigue crack growth by the cyclic J-integral, the effects of the transverse strain on J are investigated. The results suggest that the fatigue life prediction based on multiaxial fatigue theories and the critical plane approach should include the constraint effects due to the transverse strain. Consequently, the concept of constant fatigue life contour on the Γ-plane in multiaxial fatigue theories is generalized to the constant fatigue life surface in the Γ-space where the shear strain and the two normal strains are the three axes. Finally, a damage parameter as a function of the shear strain and the two normal strains is proposed for evaluation of fatigue damage under multiaxial loading conditions.     

On anti-plane shear behavior of a Griffith permeable crack in piezoelectric materials by use of the non-local theory

In this paper, the non-local theory of elasticity is applied to obtain the behavior of a Griffith crack in the piezoelectric materials under anti-plane shear loading for permeable crack surface conditions. By means of the Fourier transform the problem can be solved with the help of a pair of dual integral equations with the unknown variable being the jump of the displacement across the crack surfaces. These equations are solved by the Schmidt method. Numerical examples are provided. Unlike the classical elasticity solutions, it is found that no stress and electric displacement singularity is present at the crack tip. The non-local elastic solutions yield a finite hoop stress at the crack tip, thus allowing for a fracture criterion based on the maximum stress hypothesis. The finite hoop stress at the crack tip depends on the crack length and the lattice parameter of the materials, respectively. The project supported by the National Natural Science Foundation of China (50232030 and 10172030)     

Thermal response of a partially insulated interface crack in a graded coating-substrate structure under thermo-mechanical disturbance: Energy release and density

Consider the thermal fracture problem of a functionally graded coating-substrate structure of finite thickness with a partially insulated interface crack subjected to thermal-mechanical supply. A new model is proposed that the heat conduction through the crack region occurs and the temperature drop across the crack surfaces is the result of the thermal resistance. For the first time, real fundamental solutions are derived for the fracture analysis of functionally graded materials. The complicated mixed boundary problems of equations of heat conduction and elasticity are converted analytically into singular integral equations, which are solved numerically. The asymptotic expressions with higher order terms for the singular integral kernels are considered to improve the accuracy and efficiency of the numerical integration. Explicit expressions of various failure modes including stress intensity factors, energy release rate and strain energy density, are provided. Numerical results are presented to illustrate the effects of non-homogeneity parameters and the dimensionless thermal resistance on the temperature distribution along the crack surfaces and extended crack line, the thermal stress intensity factors and minimum strain energy density.     

Bonded half planes containing an arbitrarily oriented crack

The plane elastostatic problem for two bonded half planes containing an arbitrarily oriented crack in the neighborhood of the interface is considered. Using Mellin Transforms the problem is formulated as a system of singular integral equations. The equations are solved for various crack orientations, material combinations, and external loads. The numerical results given in the paper include the stress intensity factors, the strain energy release rates, and the probable cleavage angles giving the direction of crack propagation.     

Crack emanating from an open notch

This paper is concerned with the plane elastostatic problem of a crack emanating from an open notch where it is assumed that the crack is directed parallel to one side of the notch. The problem is reduced to two simultaneous Fredholm integral equations which are then solved numerically for three loading conditions of the crack faces. In particular, it is shown that the crack tip intensity factors are almost independent of crack length if the singular stress field associated with the uncracked notch is used as the loading of the crack faces.     

导电薄板内电流密度分布与反平面剪切的比拟

定量分析电流密度在含裂纹载流薄板内的分布是当前利用电流热效应止裂技术中一个首先要解决的问题．由于裂纹的存在,电流密度在裂尖形成带奇异性分布的高度密集．现有的分析方法往往比较复杂或局限于特殊布置形式的裂纹．通过电流密度分布与弹性力学里反平面剪切问题的比拟,把分析含裂纹载流薄板内电流密度的分布等效于考虑相应的III型裂纹问题,并比照III型裂纹的应力强度因子来定义电流密度因子．而对于裂纹问题的处理可采用分布位错法这一断裂力学里便利有效的分析手段．由给出的算例可见,所提出的比拟解法可以方便精确地求解电流密度在裂尖附近的奇异分布,并有助于对这一奇异性在概念上的直观理解．     

SCATTERING OF HARMONIC ANTI-PLANE SHEAR STRESS WAVES BY A CRACK IN FUNCTIONALLY GRADED PIEZOELECTRIC/PIEZOMAGNETIC MATERIALS

In this paper, the dynamic behavior of a permeable crack in functionally graded piezoelectric/piezomagnetic materials is investigated. To make the analysis tractable, it is assumed that the material properties vary exponentially with the coordinate parallel to the crack. By using the Fourier transform, the problem can be solved with the help of a pair of dual integral equations in which the unknown is the jump of displacements across the crack surfaces. These equations are solved to obtain the relations between the electric filed, the magnetic flux field and the dynamic stress field near the crack tips using the Schmidt method. Numerical examples are provided to show the effect of the functionally graded parameter and the circular frequency of the incident waves upon the stress, the electric displacement and the magnetic flux intensity factors of the crack.