DYNAMIC STRESS FIELD AROUND THE MODE Ⅲ CRACK TIP IN AN ORTHOTROPIC FUNCTIONALLY GRADED MATERIAL

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｒｅｃｅｎｔｙｅａｒｓ,ｇｒｅａｔａｔｔｅｎｔｉｏｎｓｈａｖｅｂｅｅｎｐａｉｄｔｏｔｈｅｒｅｓｅａｒｃｈｏｆＦｕｎｃｔｉｏｎａｌｌｙＧｒａｄｅｄＭａｔｅｒｉａｌｓ(ＦＧＭ).Ｆｒｏｍｔｈｅｖｉｅｗｐｏｉｎｔｓｏｆａｐｐｌｉｅｄｍｅｃｈａｎｉｃｓ,ＦＧＭａｒｅｎｏｎ＿ｈｏｍｏｇｅｎｅｏｕｓｓｏｌｉｄｓ.Ｔｈｅｎｏｎ＿ｈｏｍｏｇｅｎｅｉｔｙｏｆＦＧＭｈａｓａｇｒｅａｔｉｎｆｌｕｅｎｃｅｏｎｔｈｅｉｒｍｅｃｈａｎｉｃａｌｂｅｈａｖｉｏｒ,ｅｓｐｅｃｉａｌｌｙｗｈｅｎｔｈｅｃｏｍｐ…     

Impact failure prediction of Mode III crack in orthotropic functionally graded strip

Mode III impact of a crack in an orthotropic functionally graded strip is investigated. The shear moduli in two directions of the material are assumed to vary proportionately with gradient. Laplace transform and Fourier cosine transform are used to reduce the problem to solving a Fredholm integral equation. The crack tip stress field is obtained by considering the asymptotic behavior of Bessel function. Energy density factor criterion is applied to obtain the maximum of minimum energy density and direction of crack initiation. Numerical results are given graphically. The effects of orthotropy, nonhomogeneity and height of the strip on the energy density factor are discussed.     

DYNAMIC STRESS FIELD AROUND THE MODE Ⅲ CRACK TIP IN AN ORTHOTROPIC FUNCTIONALLY GRADED MATERIAL

The problem of a Griffith crack in an unbounded orthotropic functionally graded material subjected to antipole shear impact was studied. The shear moduli in two directions of the functionally graded material were assumed to vary proportionately as definite gradient. By using integral transforms and dual integral equations, the local dynamic stress field was obtained. The results of dynamic stress intensity factor show that increasing shear moduli’s gradient of FGM or increasing the shear modulus in direction perpendicular to crack surface can restrain the magnitude of dynamic stress intensity factor.     

Damage accumulation and crack growth in bilinear materials with softening: Application of strain energy density theory

A pseudo-elastic damage-accumulation model is developed by application of the strain energy density theory. The three-point bending specimen is analyzed to illustrate the crack growth characteristics according to a linear elastic softening constitutive law that is typical of concrete materials. Damage accumulation is accounted for by the decrease of elastic modulus and fracture toughness. Both of these effects are assessed by means of the strain energy density functions in the elements around a slowly moving crack. The rate of change of the strain energy density factor S with crack growth as expressed by the relation dS/da = constant is shown to describe the failure behavior of concrete. Results are obtained for different loading steps that yield different slopes of lines in an S versus a (crack length) plot. The lines rotate about the common intersect in an anti-clockwise direction as the load steps are increased. The intersect shifts upward according to increase in the specimen size. In this way, the combined interaction of material properties, load steps and specimen geometry and size are easily analyzed in terms of the failure mode or behavior that can change from the very brittle to the ductile involving stable crack growth. An upper limit on specimen or structural size is established beyond which stable crack growth ceases to occur and failure corresponds to unstable crack propagation or catastrophic fracture. The parameters that control the failure mode are the threshold values of the strain energy density function (dW/dV)_c and the strain energy density factor S_c.     

Failure instability of a debonded interface in the presence of a main crack

The problem of fracture initiating from an edge crack in a nonhomogeneous beam made of two dissimilar linear elastic materials that are partially bonded along a common interface is studied by the strain energy density theory. The beam is subjected to three-point bending and the unbonded part of the interface is symmetrically located with regard to the applied loading. The applied load acts on the stiffer material, while the edge crack lies in the softer material. Fracture initiation from the tip of the edge crack and global instability of the composite beam are studied by considering both the local and global stationary values of the strain energy density function, dW/dV. A length parameter l defined by the relative distance between the maximum of the local and global minima of dW/dV is determined for evaluating the stability of failure initiation by fracture. Predictions on critical loads for fracture initiation from the tip of the edge crack, crack trajectories and fracture instability are made. In the analysis the load, the length of the edge crack and the length and position of the interfacial crack remained unchanged. The influence of the ratio of the moduli of elasticity of the two materials, the position of the edge crack and the width of the stiffer material on the local and global instability of the beam was examined. A general trend is that the critical load for crack initiation and fracture instability is enhanced as the width and the modulus of elasticity of the stiffer material increase. Thus, the stiffer material acts as a barrier in load transfer.     

Growth of crack caused by temperature gradients with change in surface insulation

This work is concerned with thermoelastic stress and failure analysis of a centrally cracked panel subjected to temperature gradients while the insulation on the crack surface is varied. The corresponding temperature and thermoelastic stress fields are obtained by application of the finite element method. According to the strain energy density criterion, the crack grows incrementally when the maximum of the minimum strain energy density function reaches a critical value for a given material. Crack growth resistance curves involving plots of the strain energy density factor S versus the half crack length a are developed for crack surfaces with varying degree of heat resistance. The resulting curves are straight lines satisfying the condition dS/da = const. and useful for determining combined influence of thermal loading and structural geometry that lead to global instability.     

Crack initiation direction from interface of bonded dissimilar media

Debonded region of an interface between two dissimilar materials are modeled as a line crack that tends to enhance the initiation of failure by fracture. Depending on the load that interacts with dissimilar materials, no a priori knowledge of how failure would initiate from an existing interface crack is assumed. By application of the strain energy density criterion, potential crack initiation sites are obtained for different biaxial loading states and materials with dissimilar properties.Numerical results are obtained for an epoxy/aluminum medium. In each case, a finite line segment of debonding is assumed. Uniform stresses are applied normal and parallel to the interface so that a biaxial load factor k determines the relative magnitude of biaxiality. Positive and negative k correspond, respectively, to applied tension and compression parallel to the interface. For a fixed ratio of the elastic moduli, crack initiation angles measured from the interface would increase with positive k and decrease with an increase of negative k. These findings are presented for different values of k. The direction of maximum yield initiation could also be determined from the stationary values of the strain density function. These locations are identified with elements that undergo excessive distortion while the possible fracture sites are assumed to coincide with regions where dilatational effects would dominate.     

A pseudo-elastic material damage model, part II: Application to the scaling of specimen size, material behavior and loading rate

A failure criterion is presented which relates the strain energy density of the material to both yielding and fracture. Cumulative material damage throughout a structural component may be monitored and the relative influence of yielding and stable crack growth assessed. The criterion is demonstrated, using finite element analysis, for center cracked panel specimens differing by material toughness values. From crack growth increment predictions using the uniaxial stress-strain behavior of the material, the criterion predicts the critical value of the strain energy density factor S_c governing crack instability.     

ANALYSIS OF A CRACK IN A FUNCTIONALLY GRADED STRIP WITH A POWER FORM SHEAR MODULUS

The plane strain problem of a crack in a functionally graded strip with a power form shear modulus is studied. The governing equation in terms of Airy＇s stress function is solved exactly by means of Fourier transform. The mixed boundary problem is then reduced to a system of singular integral equations and is solved numerically to obtain the stress intensity factor at crack-tip. The maximum circumferential stress criterion and the strain energy density criterion are both employed to predict the direction of crack initiation. Numerical examples are given to show the influence of the material gradation models and the crack sizes on the mode-I and mode-II stress intensity factors. The dependence of the critical kink-angle on the crack size is examined and it is found that the crack kink-angle decreases with the increase of the normalized crack length, indicating that a longer crack tends to follow the original crack-line while it is much easier for a shorter crack to deviate from the original crack-line.     

Influence of indenter geometry on half-plane with edge crack subjected to fretting condition

An edge crack is analyzed to study fretting failure. A flat punch with rounded corners and a half-plane are regarded as an indenter and a substrate, respectively. Plane strain condition is considered. Contact shear traction in the case of partial slip is evaluated numerically. It is assumed that an initial crack is extended to the point of minimum strain energy density in the half-plane from the trailing edge of contact. Dislocation density function method is used to evaluate K_I and K_II. The variations of K_I and K_II during crack growth are examined in the case of indentation by a punch with different ratio of the flat region (l) to the punch width (L). Sih's minimum strain energy density theory [1] is also applied to predict the propagation direction of the initial crack. The direction evaluated is similar to that found in the experiment. Stress intensity factor ranges (ΔK_I and ΔK_II) are examined during cyclic shear on the contact. For the design of contacting bodies, a suggestible geometry of punch for alleviating cracking failure is studied.     

Stable growth of a crack interacting with a circular inclusion

The problem of failure of a plate containing a circular inclusion and a crack is studied. The crack is oriented along a diameter of the inclusion and the plate is subjected to a remote uniaxial stress perpendicular to the crack axis. The process of slow stable crack growth from initiation to termination is studied by the strain energy density theory. The crack growth is simulated by predicting finite increments of crack extension when material elements near the crack tip absorb a critical amount of strain energy density level, . Unstable crack growth occurs when the strain energy density factor S reaches a critical value where r_c is the critical size of the final crack increment prior to instability. The stress at crack initiation and the critical stress and crack length at failure are determined. The influence of the mechanical properties of the plate and the inclusion, the relative position of the inclusion and the crack and the crack length on the characteristic quantities of stable crack growth is analyzed. The dependence of the stable crack growth process on the loading rate is also investigated. Results are displayed in graphical form.     

Strain energy density prediction of crack initiation and direction in cracked T-beams and pipes

In this paper, the S-theory is applied to determine crack initiation and direction for cracked T-beams and circumferentially cracked pipes. It makes use of a parameter called strain energy density factor, S, which is a function of the stress intensity factors. The strain energy density theory provides a more general treatment of fracture mechanics problems by virtue of its ability in describing the multiscale feature of material damage and in dealing with mixed mode crack propagation problem. A simple method for obtaining approximate stress intensity factors is also applied. It takes into account the elastic crack tip stress singularity while using the elementary beam theory. Some basic loading conditions in beams and pipes are studied.     

Initiation and growth characterization of corner cracks near circular hole

The finite element analysis of crack problems often incorporates the asymptotic character of the local solution into the formulation. Embedment of stress or strain singularities can impose serious restrictions on the outcome and inconsistencies in predicting crack and/or growth. These restrictions are discussed in connection with the problem of two diametrically opposite corner cracks near a circular hole subjected to remote uniform tension. Enforced in the numerical treatment is the 1/r character of the strain energy density function local to the corner crack border where r is the radial distance measured from the crack front. The tendency for the corner crack to become a through crack is predicted by assuming that each point of the crack border extends by an amount proportional to the strain energy density factor. The path would correspond to the loci of minimum strain energy density function. Numerical results are displayed graphically and discussed in connection with crack initiation and non-self-similar crack growth.     

Crack growth resistance in non-perfect plasticity: Isotropic versus kinematic hardening

The Strain Energy Density Theory is applied for analyzing energy dissipation and crack growth in the three-point bending specimen when the material behavior follows a multilinear strain-hardening stress-strain relationship. The problem is solved through the application of incremental theory of plasticity and finite element method.The rate of change of the strain energy density factor S with crack length a is verified to be governed by the relation . Results are obtained for isotropic and kinematic hardening. Moreover, the effects of loading step and specimen size are pointed out.     

Crack-tip stress fields in functionally graded materials with linearly varying properties

Crack-tip stress fields for a stationary crack along or inclined to the direction of property gradation in functionally graded materials (FGMs) are obtained through an asymptotic analysis coupled with Westergaard’s stress function approach. The elastic modulus of the FGM is assumed to vary linearly along the gradation direction. The first six terms for a crack along the direction of property variation and first four terms for a crack inclined to the direction of property variation in the expansion of the stress field are derived to explicitly bring out the influence of nonhomogeneity on the structure of the stress field. Using these stress fields, contours of constant maximum shear stress and constant out of plane displacement are generated and the effect of inclination of property gradation direction on these contours is discussed. The strain energy density criterion is applied to obtain critical conditions for crack initiation and the effect of property gradation is discussed. It is shown that the materials with varying properties can offer more resistance to crack propagation and will suppress crack growth in some situations.     

Evaluation of crack initiation angle under mixed mode loading at diverse strain rates

Crack initiation angle, under mixed mode loading at several strain rates, is analysed using an experimental–numerical approach. The physical phenomenon for the problem at hand is influenced by the local and global conditions. One of such factors is the strain rate at the crack tip. For this purpose, PMMA plates with centred angled cracks under mixed mode loading were tested. The strain rate at the neighbourhood of the crack tip before crack propagation was evaluated. Considering that this material is strain rate sensitive, the numerical models were calibrated with the modulus of elasticity measured in tension tests at the observed strain rates. Numerical evaluations were performed with the finite element method in conjunction with the volume energy density criterion. An improvement in the evaluation of the crack propagation angle was observed. In order to complete the analysis, the crack initiation angle was also evaluated with the strain energy density factor S, considering the mechanical properties of PMMA, as evaluated at the observed strain rates, and the stress intensity factors k₁ and k₂. Results are in agreement with those observed experimentally.     

An extended strain energy density failure criterion by differentiating volumetric and distortional deformation

We extend Sih’s strain energy density criterion (Sih, 1974) for crack kinks and material failure by weighting differently the volumetric and distortional parts in the extended strain energy density factor. The work is inspired by the factor that failure by microscopic shearing governed by distortion and microscopic separation controlled by hydrostatic tension represent distinct deformation processes, and should be treated differently as we count their influences to material failure. With the weight parameter introduced to the extended strain energy density factor criterion, we explain satisfactorily several critical experiments which reported crack kink in samples subjected to mixed-mode loading. The extended strain energy density idea is also used to derive a generalized pressure-dependent yielding criterion, which supplies a theoretical basis for those novel strength criteria for materials like bulk metallic glasses. Corresponding methods to determine the two material parameters, the critical strain energy density factor and the weight parameter quantifying the relative contribution by distortion over volumetric deformation, are discussed.     

Non-self-similar crack growth in elastic-plastic finite thickness plate

This work is concerned with non-self-similar crack growth in medium strength metal plates while the loading step, plate thickness and material properties are altered. The three-dimensional elastic-plastic finite element stress analysis is combined with the strain energy density criterion for modeling the material damage process from crack initiation to final global instability including the intervening stage of slow crack growth. Both inelastic deformation and crack growth are accounted for each increment of loading such that the redistribution of stresses and strains are made for each new crack profile. Numerical results are obtained for the center cracked plate configuration under uniform extension with twenty-seven (27) different combinations of specimen thickness, loading step and material type. The fracture toughness S_c being related to K_1c for three different materials are predicted analytically from the corresponding uniaxial tensile test data. Effective strain energy density factor and half crack length are defined so that the results can be compared with their two-dimensional counterparts. Crack growth resistance curves (R-curves) are constructed by plotting as a function of . The condition is found to prevail during slow crack growth. Translation and/or rotation of the lines can yield results other than those calculated and serve a useful purpose for scaling component size and test time. The minimum thickness requirement for the ASTM valid K_1c test is also discussed in connection with predictions based on the strain energy density criterion. The corresponding K_1c for smaller specimens that exhibit moderate ductility and nonlinearity can also be obtained analytically. In such cases, the influence of loading step can be significant and should not be neglected. Notwithstanding the shortcomings of the theory of plasticity, the qualitative features of non-self-similar crack growth are predicted by the strain energy density criterion. Any refinements on the analytical modeling of the material damage process would only affect the results qualitatively, a subject that is left for future investigation.     

Crack-extension resistance of polycarbonate material

Crack-extension resistance for the polycarbonate material is examined by application of the strain energy density criterion and the incremental theory of plasticity. The energy state ahead of a slow moving crack in a three-point bend specimen is obtained for each load increment and used to determine the crack growth characteristics. The analytical results are displayed by plotting the strain energy density factor S as a function of crack length and compared with available experimental data on the polycarbonate material. Standard deviations and mean errors are computed for the experimentally measured and analytically determined values of S and are shown to be much lower than those based on the J-integral parameter. Modeling of the polycarbonate material by the theory of plasticity still remains much to be desired. Crack growth calculations are performed for a strain hardening parameter α = 0.85 that controls the proportion of isotropic and kinematic hardening. Nevertheless, the criterion dS/da = const. is shown to collate well with the experimental crack growth data.     

Torsional impact of transversely isotropic solid with functionally graded shear moduli and a penny-shaped crack

The torsional impact response of a penny-shaped crack in an unbounded transversely isotropic solid is considered. The shear moduli are assumed to be functionally graded such that the mathematics is tractable. Laplace transform and Hankel transform are used to reduce the problem to solving a Fredholm integral equation. The crack tip stress fields are obtained. Investigated are the influence of material nonhomogeneity and orthotropy on the dynamic stress intensity factor. The peak value of the dynamic stress intensity factor can be suppressed by increasing the shear moduli's gradient and/or increasing the shear modulus in a direction perpendicular to the crack surface.