功能梯度材料亚表面缺陷的热波多重散射问题

基于热传导波动模型,采用波函数展开法,研究了半无限功能梯度材料亚表面球形缺陷的热波多重散射.给出了热波散射的一般解.温度波由调制光束在材料表面激发,球形缺陷表面的边界条件为绝热,非均匀参数为指数函数变化.分析了结构几何参数和物理参数对温度分布的影响,并给出了温度变化的数值结果.本研究可为功能梯度材料的分析研究、物理反问题和红外热波成像等提供理论基础和参考数据.     

Scattering of the harmonic anti-plane shear waves by a crack in functionally graded piezoelectric materials

In the present paper, the dynamic behavior of a Griffth crack in the functionally graded piezoelectric material (FGPM) is investigated. It is assumed that the elastic stiffness, piezoelectric constant, dielectric permittivity and mass density of the FGPM vary continuously as an exponential function, and that FGPM is under the anti-plane mechanical loading and in-plane electrical loading. By using the Fourier transform and defining the jumps of displacement and electric potential components across the crack surface as the unknown functions, two pairs of dual integral equations are derived. To solve the dual integral equations, the jumps of the displacement and electric potential components across the crack surface are expanded in a series of Jacobi polynomial. Numerical examples are provided to show the effects of material properties on the stress and the electric displacement intensity factors.     

Size effect on dynamic stability of functionally graded microbeams based on a modified couple stress theory

Dynamic stability of microbeams made of functionally graded materials (FGMs) is investigated in this paper based on the modified couple stress theory and Timoshenko beam theory. This non-classical Timoshenko beam model contains a material length scale parameter and can interpret the size effect. The material properties of FGM microbeams are assumed to vary in the thickness direction and are estimated though Mori–Tanaka homogenization technique. The higher-order governing equations and boundary conditions are derived by using the Hamilton’s principle. The differential quadrature (DQ) method is employed to convert the governing differential equations into a linear system of Mathieu–Hill equations from which the boundary points on the unstable regions are determined by Bolotin’s method. Free vibration and static buckling are also discussed as subset problems. A parametric study is conducted to investigate the influences of the length scale parameter, gradient index and length-to-thickness ratio on the dynamic stability characteristics of FGM microbeams with hinged–hinged and clamped–clamped end supports. Results show that the size effect on the dynamic stability characteristics is significant only when the thickness of beam has a similar value to the material length scale parameter.     

Dynamic stress intensity factors around two parallel cracks in a functionally graded layer bonded to dissimilar half-planes subjected to anti-plane incident harmonic stress waves

The time-harmonic problem of determining the stress field around two parallel cracks in functionally graded materials (FGMs) is studied. The Fourier transform technique is used to reduce the boundary conditions to four simultaneous integral equations which are then solved by expanding the differences of crack surface displacements in a series. The unknown coefficients in the series are obtained by the Schmidt method. Numerical calculations are carried out for dynamic stress intensity factors (DSIF) in FGMs.     

Free vibration and stability of functionally graded circular cylindrical shells according to a 2D higher-order deformation theory

A two-dimensional (2D) higher-order deformation theory is presented for vibration and buckling problems of circular cylindrical shells made of functionally graded materials (FGMs). The modulus of elasticity of functionally graded (FG) shells is assumed to vary according to a power law distribution in terms of the volume fractions of the constituents. By using the method of power series expansion of continuous displacement components, a set of fundamental governing equations which can take into account the effects of both transverse shear and normal deformations, and rotatory inertia is derived through Hamilton’s principle. Several sets of truncated Mth order approximate theories are applied to solve the eigenvalue problems of simply supported FG circular cylindrical shells. In order to assure the accuracy of the present theory, convergence properties of the fundamental natural frequency for the fundamental mode r=s=1 are examined in detail. A comparison of the present natural frequencies of isotropic and FG shells is also made with previously published results. Critical buckling stresses of simply supported FG circular cylindrical shells subjected to axial stress are also obtained and a relation between the buckling stress and natural frequency is presented. The internal and external works are calculated and compared to prove the numerical accuracy of solutions. Modal transverse shear and normal stresses are calculated by integrating the three-dimensional (3D) equations of motion in the thickness direction satisfying the stress boundary conditions at the outer and inner surfaces. The 2D higher-order deformation theory has an advantage in the analysis of vibration and buckling problems of FG circular cylindrical shells.     

半无限功能梯度材料结构中圆孔对弹性波的多重散射

基于弹性波多体散射理论,采用波函数展开法,研究了半无限指数梯度材料中圆孔对弹性波的多重散射和动应力集中,得到了问题的解析解,给出了圆孔动应力集中系数的数值解,分析了圆孔与边界的距离、入射波波数以及材料的非均匀参数对圆孔周围动应力集中系数的影响。分析表明:梯度材料的非均匀参数小于零时对最大动应力影响较小,但是对动应力在圆周的分布有较大影响;大于零时对最大动应力和动应力在圆周的分布影响都很大,特别是在圆孔与边界的距离较小时影响更大。     

Investigation of anti-plane shear behavior of a Griffith permeable crack in functionally graded 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 functionally graded piezoelectric materials under the anti-plane shear loading for the permeable electric boundary conditions. To make the analysis tractable, it is assumed that the material properties vary exponentially with coordinate vertical to the crack. By means of the Fourier transform, the problem can be solved with the help of a pair of dual-integral equations that the unknown variable is the jump of the displacement across the crack surfaces. These equations are solved by use of the Schmidt method. Numerical examples are provided. Unlike the classical elasticity solutions, it is found that no stress and electric displacement singularities are present near the crack tips. The non-local elastic solutions yield a finite hoop stress at the crack tips, thus allows us to using the maximum stress as a fracture criterion. The finite hoop stresses at the crack tips depend on the crack length, the functionally graded parameter and the lattice parameter of the materials, respectively.     

Influence of microstructure on fracture toughness distribution in ceramic-metal functionally graded materials

This paper deals with the influence of microstructure on fracture toughness distribution in functionally graded materials (FGMs) consisting of partially stabilized zirconia (PSZ) and austenitic stainless steel SUS 304. FGMs and non-graded composites (non-FGMs) with fine and coarse microstructures are fabricated by powder metallurgy using PSZ and two kinds of SUS 304 powders. The fracture toughness is determined by conventional tests for several non-FGMs with each material composition and by a method utilizing stable crack growth in FGMs. The obtained results on the fracture toughness are as follows: (1) The fracture toughness increases with an increase in a content of SUS 304 on both FGMs and non-FGMs. (2) On the fracture toughness of the non-FGMs, the influence of microstructure is negligible. (3) On the FGMs, the fracture toughness is higher in the FGM with fine microstructure than in the FGM with coarse microstructure. (4) The fracture toughness of the FGMs is higher than that of the non-FGMs especially in the case of fine microstructure. Finally, the residual stress in the FGMs created in a fabrication process is estimated from the difference in fracture toughness between the FGMs and non-FGMs.     

梯度材料中任意形孔对反平面剪切波散射与动应力

该文基于弹性动力学理论,采用复变函数与保角映射方法,研究了指数梯度材料中任意形孔洞对弹性波的散射与动应力集中,给出了问题的解析解.并以求解椭圆孔动应力集中系数为例,分析了入射波数和材料非均匀参数等对椭圆孔动应力分布的影响.     

Finite element evaluation of mixed mode stress intensity factors in functionally graded materials

This paper is directed towards finite element computation of fracture parameters in functionally graded material (FGM) assemblages of arbitrary geometry with stationary cracks. Graded finite elements are developed where the elastic moduli are smooth functions of spatial co‐ordinates which are integrated into the element stiffness matrix. In particular, stress intensity factors for mode I and mixed‐mode two‐dimensional problems are evaluated and compared through three different approaches tailored for FGMs: path‐independent J^*_k‐integral, modified crack‐closure integral method, and displacement correlation technique. The accuracy of these methods is discussed based on comparison with available theoretical, experimental or numerical solutions. Copyright © 2001 John Wiley & Sons, Ltd.     

Effect of a graded interlayer on the stress intensity factor of cracks in a joint under thermal loading

In a bimaterial joint with and without a graded interlayer, the stress intensity factor of cracks perpendicular to the interface was calculated for a thermal loading by a homogeneous change in the temperature. In joints without an interlayer, the stress intensity factor increases to infinity as the crack approaches the interface for the case of the Young’s moduli E₁/E₂>1 (crack in material 1). Introducing a graded interlayer with a continuous transition in the material properties between the two joined materials leads to a continuous change in the stress intensity factor if the crack propagates from material 1 into material 2. Results are presented for different transition functions of the material properties and for different thickness ratios of the layers. The possible beneficial effect of a graded interlayer is discussed.     

Consideration of spatial variation of the length scale parameter in static and dynamic analyses of functionally graded annular and circular micro-plates

This article introduces new methods for static and free vibration analyses of functionally graded annular and circular micro-plates, which can take into account spatial variation of the length scale parameter. The underlying higher order continuum theory behind the proposed approaches is the modified couple stress theory. A unified way of expressing the displacement field is adopted so as to produce numerical results for three different plate theories, which are Kirchhoff plate theory (KPT), Mindlin plate theory (MPT), and third-order shear deformation theory (TSDT). Governing partial differential equations and corresponding boundary conditions are obtained following the variational approach and the Hamilton's principle. Derived systems of differential equations are solved numerically by utilizing the differential quadrature method (DQM). Comparisons to the results available in the literature demonstrate the high level of accuracy of the numerical results generated through the developed methods. Extensive analyses are presented in order to illustrate the influences of various geometric and material parameters upon static deformation profiles, stresses, and natural vibration frequencies. In particular, the length scale parameter ratio -which defines the length scale parameter variation profile-is shown to possess a profound impact on both static and dynamic behaviors of functionally graded annular and circular micro-plates.     

Analysis of a transiently propagating crack in functionally graded materials under mode I and II

Crack tip stress and displacement fields for a transiently propagating crack along gradient in functionally graded materials (FGMs) with a linear variation of shear modulus are developed. The higher order terms of the transient stress and displacement fields at crack tip were obtained by transforming the general partial differential equations of the dynamic equilibrium into Laplace’s equations whose solutions have harmonic functions. Thus, the fields can be expressed very simply. Using these stress components, isochromatics and the first invariant at crack tip are generated.The results show that the isochromatics (constant maximum shear stress) for mode I crack tilt backward around the crack tip with an increase of crack tip acceleration , and tilt forward around the crack tip with an increase of rate of change of dynamic mode I stress intensity factor . The isochromatics for mixed mode crack move to upper direction with an increases of and , and lower direction with an increase of . Contours of the first stress invariant for mode I crack enlarge around the crack tip with an increase of , and decrease around the crack tip with an increase of . As decreases at crack initiation, the predicted kinking angles increase. As increases, the predicted kinking angles also increase.     

Dynamic stress field for torsional impact of a penny-shaped crack in a transversely isotropic functional graded strip

The torsional impact response of a penny-shaped crack in a transversely isotropic strip is considered. The shear moduli are assumed to be functionally graded such that the mathematics is tractable. Laplace and Hankel transforms 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. Investigated are the effects of material nonhomogeneity and orthotropy and strip’s highness on the dynamic stress intensity factor. The peak 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. The dynamic behavior varies little with the increasing of the strip’s highness.     

带形域内圆柱形夹杂对SH型导波的散射

本文采用波函数展开法建立了带形域内圆形夹杂对稳态SH型导波散射的级数解答。该问题采用镜像法求解,利用圆柱形夹杂边界位移和应力的连续性条件,并借助内域型Graf加法公式将解答归结为一组无穷代数方程组的求解,截断该无穷代数方程组可求得该问题的数值结果。为了便于问题的分析,算例给出了带形域中最低阶对称波型的SH型导波入射时圆柱形夹杂周边的动应力分布情况。结果表明,带形域自由边界的存在以及夹杂的几何位置对动应力集中分布有较大的影响。     

Dynamic stress intensity factors due to scattering of harmonic waves by a crack in an infinite medium

An analytical method for calculating dynamic stress intensity factors in the mixed mode (combination of opening and sliding modes) using complex functions theory is presented. The crack is in infinite medium and subjected to the plane harmonic waves. The basis of the method is grounded on solving the two‐dimensional wave equations in the frequency domain and complex plane using mapping technique. In this domain, solution of the resulting partial differential equations is found in the series of the Hankel functions with unknown coefficients. Applying the boundary conditions of the crack, these coefficients are calculated. After solving the wave equations, the stress and displacement fields, also the J‐integrals are obtained. Finally using the J‐integrals, dynamic stress intensity factors are calculated. Numerical results including the values of dynamic stress intensity factors for a crack in an infinite medium subjected to the dilatation and shear harmonic waves are presented.     

Determination of thermal stress intensity factors for an interface crack in a graded orthotropic coating-substrate structure

The thermal fracture problem of an interface crack between a graded orthotropic coating and the homogeneous substrate is investigated by two different approaches. For the case that most of the material properties in the graded orthotropic coating are assumed to vary as an exponential function, the integral transform and singular integral equation technique is used to obtain some analytical results. In order to analyze the case with more complex material distribution, an interaction integral is presented to evaluate the thermal stress intensity factors of cracked functionally graded materials (FGMs), and then the element-free Galerkin method (EFGM) is developed to obtain the final numerical results. The good agreement is obtained between the numerical results and the analytical ones. In addition, the influence of material gradient parameters and material distribution on the thermal fracture behavior is also presented.