Functionally graded materials (FGMs) are composites with continuous material property variations. The distinct interfaces between the reinforcement and the matrix in classical composites are potential damage initiation sites. The concept of FGM aims at avoiding the material mismatch at the interfaces. Functionally graded materials originated from the need for a material that has high-toughness at very high operating temperatures that occur in rocket nozzles and aeroplane engines. One of the early applications of graded materials can be thus found in thermal barrier coatings of gas turbine blades. Recent applications of FGMs include optoelectronics, ballistic impact resistance structures, wear resistant coatings and others. Although the manufacturing and applications of FGMs are well developed the basic mechanics of failure is not well understood, which is important in developing engineering design methodologies. Modern day design practice uses the concepts of fracture mechanics and the fracture properties of graded materials is not well understood. Most studies in the literature have assumed that the material response of the bulk functionally graded material to be elastic even though the constituents are nominally ductile. Some asymptotic analysis available in the literature have described the effect of ductility on the fracture parameters. However, these analysis are not complete in the sense that they have some undetermined constants. The present thesis aims at performing whole-field finite element (FE) simulations of a crack embedded in a ductile functionally graded material subjected to an anti-plane shear (mode-3) loading. A J2-deformation theory based power-law hardening nonlinear material response is assumed. The material property variation is assumed to be in the radial-direction (r-FGM), tangential to the crack (x-FGM), normal to the crack plane (y-FGM) and also at an arbitrary angle to the crack-plane (xy-FGM). Yet another power law described the material property variation. The competition between the indices of the hardening and material property variation is understood by performing a parametric analysis by varying both systematically. Our results indicate that the first most singular term of the asymptotic series remains unaffected. For some values of the material property variation index, the second asymptotic term is affected. The semi-closed form solutions available in the literature were unable to decipher the relative range of dominance of the first and second terms. From the present whole-field FEM analysis were able to extract this relative range of dominance. Our results indicate the range of dominance of the first term is least for FGMs when the material property variation is in the direction to the crack (x-FGM), and it is more for y-FGM.

No description available.

Functionally graded materials (FGMs) are composites that have continuously varying material properties, which eliminate undesirable stress concentrations that might otherwise occur in layered composites. The concept of inhomogeneously varying properties is observed in nature; examples include bones, teeth, shells and timber. Modern engineering applications of FGMs include thermal barrier coatings, wear-resistant coatings, biomedical implants and MEMS devices. Syntactic foams, particle filled nano-composites are examples of inhomogeneous materials of current interest. Analyses and experiments available in the literature have focused on characterizing the inhomogeneous material modulus and density variations. Common techniques employed are nano-indentation and wave propagation studies. There are a few fracture mechanics analyses and experiments available in the literature; most of which are devoted to measuring the fracture toughness of graded materials. A few fracture analyses of graded materials are devoted to deriving asymptotic stress, strain and displacement fields around stationary and steadily growing cracks in inhomogeneous materials. Only a few studies exist that deal with understanding the effect of material property inhomogeneity on the spontaneous crack propagation. In the present thesis the effect of material property inhomogeneity on the dynamic fracture mechanics of cracks in FGMs is described. Numerical analysis of the elastodynamic initial boundary value problem is performed using a spectral scheme. Spectral scheme is a special numerical technique developed to simulate spontaneous, planar crack propagation in a variety of materials. The method is numerically efficient as it can be implemented on parallel machines with ease. The numerical scheme is versatile and can handle any state-and rate-dependent traction-separation laws (cohesive zone models) or frictional laws. Spectral scheme has successfully been used in simulating intersonic crack propagation, earthquake slip dynamics and also direct silicon wafer bonding process used in realizing 3D MEMS structures. In the present work, the spectral formulation accounts for the inhomogeneous variation in the material wave speeds in the medium. The effect of inhomogeneity on spontaneous crack propagation due to in-plane mixed-mode loading is also addressed here. A parametric study has been performed by varying the inhomogeneity length scales independently in the top and bottom half-spaces. The effect of inhomogeneity in shear wave speed on the dynamic stress intensity factors (SIFs) of a crack propagating in a quasi-steady-state along the interface between the two functionally graded half-spaces is studied. A symmetric hardening FGM offers the maximum fracture resistance, while the fracture resistance is minimum for a symmetric softening FGM. Our simulation shows that increasing the inhomogeneity in the wave speed leads to eliminate the overshoot in the dynamic stress intensity factor. The magnitude of the steady-state (long-time) SIF increases indicating an increase in the fracture resistance. The effect of the inhomogeneous wave speed on the mode-3 crack propagation characteristics is demonstrated by taking snapshots of the crack opening at a time interval. The magnitude of the crack sliding displacement is found to increase with increase in the inhomogeneity. The effect of the material property inhomogeneity on the mode-1 crack propagation is simulated to track the crack opening displacements. The inhomogeneity is assumed to be symmetric about the weak-plane. Our spectral scheme developed here for functionally graded material with exponential variation in the material properties is capable of simulating independent bimaterial combinations. When the graded material becomes progressively stiffer and denser (hardening), the crack opening displacement in reduced, indicating an increase in the fracture resistance. On the other hand, for the softening FGMs the crack opening displacement increases indicating a reduction in fracture toughness. It is noted that the cohesive fracture resistance on the weak-plane remains same in all the FGMs.

Design based on damage tolerance concepts has become mandatory in high technology structures. These concepts are also essential for evaluating life extension of aged structures which are in service beyond originally stipulated life. Fracture analysis of such structures in the presence of single or multiple three-dimensional flaws is essential for this approach. Surface cracks are the most commonly occurring flaws and development of accurate methods of analysis for such cracks is essential for structural integrity evaluation of newly designed or aged structures. The crack fronts of these surface flaws are usually approximated mathematically to be of either part-elliptical or part-circular in geometry. In this thesis, some of the issues related to fatigue crack growth of single and multiple surface cracks are studied in detail. Here emphasis is given to the development of simple and accurate post-processing techniques to estimate stress intensity factors for surface cracks, development and/or implementation of simple numerical methods to simulate three-dimensional single and multiple cracks in fatigue and their experimental verification. Modified virtual crack closure integral (MVCCI) technique for estimation of strain energy release rates has been improved (chapter II) to deal with curved crack front and unequal elements across the crack front. The accuracy of this method is evaluated and presented in this chapter for certain benchmark surface flaw problems. The improved MVCCI is used in the investigation of interaction between multiple surface cracks in three-dimensional solids. The interaction effects are studied for both interacting and coalescing phases as observed to occur in the growth of multiple surface cracks. Extensive numerical work is performed to study the effects of various parameters such as aspect ratio, thickness ratio, interspacing on the interaction factors. These solutions are used in formulating empirical equations to estimate interaction factors. This facilitated the development of a simple semi-analytical method to study fatigue crack growth of multiple cracks. The growth of surface cracks under fatigue loading in the finite width specimens of an aero-engine superalloy has been studied experimentally (presented in chapter III). Four configurations for single semi-elliptical cracks are considered. Fatigue crack growth is simulated by two models viz. two degrees of freedom and "multi degrees of freedom with ellipse fit'. These models are sometimes referred to as semi-analytical models as the crack growth is predicted by numerical integration combining Paris equation with an empirical form of stress intensity factor solution. In order to use two degrees of freedom model for fatigue crack growth prediction of semi-elliptical cracks, empirical solution for the Ml range of geometric parameters for stress intensity factor is required for the considered configurations. The available Newman-Raju solution is useful for this purpose within a limited range of surface crack length to width (c/W) of the specimen. Based on the present finite element results, the empirical equations are developed for extended values of c/W. It is well understood that the fatigue prediction for two-dimensional crack can be improved by inclusion of crack closure effects. Usually, in semi-analytical models for growth of surface cracks under fatigue loading, the crack closure is included as a ratio of crack closure factor at surface and depth locations of semi-elliptical crack. In the present work, this ratio for the considered material of specimens is obtained by an experimental study. The difference in characteristics of preferred propagation path between semi-elliptical crack in a finite width plate and a wide plate is clearly brought out. Current crack growth predictions for most of the structures are based on the presence of only a single crack. However, in structures several cracks may initiate simultaneously within a stress critical zone and may interact depending upon their geometry, spatial location, structure geometry and mode of loading. In this work various configurations of twin semi-elliptical cracks have been studied by experiments. The beachmarks created on the specimens during experiments are used in the investigation of crack shape progression during fatigue. A three degrees of freedom crack growth model for interacting and coalescing cracks has been proposed. The experimentally determined crack shape and lives have been compared with the corresponding values from numerical simulation. The correlation of experimental results with numerical predictions was carried out through improved MVCCI for eight-noded brick elements. This has worked well in the configurations analysed. However, it is known in literature that there are benefits of using 20-noded singular elements. There could be special situations where the regular elements could fail, and singular elements could be essential. For this purpose, further development of MVCCI were carried out using 20-noded quarter-point elements (presented in chapter IV). Also a novel technique of decomposed crack closure integral (DCCI) was developed (presented in chapter V) for both regular and singular elements to represent the variation of MVCCI more accurately along the crack front. It is well known that quarter-point elements at crack front produce the required singularity at the crack tip and give accurate stress distribution with fewer degrees of freedom than conventional elements. Thus to develop more efficient post-processing tools, the MVCCI expressions are formulated for 20-noded singular quarter-point element for various assumptions regarding stress and displacement distributions in the elements across the crack front. A comprehensive study is presented (chapter IV) on MVCCI for 20-noded singular brick element including various simplified expressions for three-dimensional part-through cracks in pure and mixed-mode state of deformation of fracture. The developed MVCCI expressions are also valid for 15-noded quarter-point Penta elements. The reduction in model size can further be obtained if 12-noded three-dimensional singular element is employed at the crack front and eight-noded elements are used away from the crack front. The MVCCI expressions are also developed for 12-noded singular element and their accuracy is evaluated by numerical solutions. Presently, MVCCI, estimates the average stress intensity factor at the center of each element along the crack front. In this thesis, a Decomposed Crack Closure Integral (DCCI) is formulated to represent an assumed variation of stress intensity factor along the crack front in each element. The DCCI is formulated for 8-noded brick, 20-noded conventional brick and 20-noded singular brick elements. The numerical examples presented here deal with three-dimensional problems of patch repair technology and part-through cracks. The technique showed a major advantage for the patch repair problems where SIF variations along the crack front are of significance and large mesh sizes are computationally expensive. This along with MVCCI for 12-noded and 20-noded singular elements formed a part of the work on development of accurate and effective post-processing tools. It is expected that the present work will be helpful in damage tolerance design and assessment of aerospace structures and the experimental work performed as a part of this thesis will enhance confidence in the damage tolerance analysis. The thesis is concluded in chapter VI presenting the contributions of this thesis and projecting future lines of work possible in this area.