Functionally graded penny-shaped cracks under dynamic loading

Dynamic load is applied to a functionally graded material with penny-shaped cracks. The materials are also transversely isotropic depending only on the axial coordinate z. The elastic region may be regarded to consist of many thin layers such that properties are constants within each layer, but they may vary from layer to layer. Laplace and Hankel transform are used in conjunction with the stiffness matrix approach. The Dual integral equations are then obtained by application of appropriate boundary and interface conditions. Stress intensity factors are then determined in the Laplace transform domain. Inversion yields the results in the time domain. Numerical examples show that multiple crack configurations in functionally graded materials can be treated where the continuously varied material properties can be divided into a finite number of layers with different properties.     

Elasticity solutions for a transversely isotropic functionally graded circular plate subject to an axisymmetric transverse load qrk

This paper considers the bending of transversely isotropic circular plates with elastic compliance coefficients being arbitrary functions of the thickness coordinate, subject to a transverse load in the form of qr^k (k is zero or a finite even number). The differential equations satisfied by stress functions for the particular problem are derived. An elaborate analysis procedure is then presented to derive these stress functions, from which the analytical expressions for the axial force, bending moment and displacements are obtained through integration. The method is then applied to the problem of transversely isotropic functionally graded circular plate subject to a uniform load, illustrating the procedure to determine the integral constants from the boundary conditions. Analytical elasticity solutions are presented for simply-supported and clamped plates, and, when degenerated, they coincide with the available solutions for an isotropic homogenous plate. Two numerical examples are finally presented to show the effect of material inhomogeneity on the elastic field in FGM plates.     

On asymmetric bending of functionally graded solid circular plates

Based on the three-dimensional elasticity equations, this paper studies the elastic bending response of a transversely isotropic functionally graded solid circular plate subject to transverse biharmonic forces applied on its top surface. The material properties can continuously and arbitrarily vary along the thickness direction. By virtue of the generalized England’s method, the problem can be solved by determining the expressions of four analytic functions. Expanding the transverse load in Fourier series along the circumferential direction eases the theoretical construction of the four analytic functions for a series of important biharmonic loads. Certain boundary conditions are then used to determine the unknown constants that are involved in the four constructed analytic functions. Numerical examples are presented to validate the proposed method. Then, we scrutinize the asymmetric bending behavior of a transversely isotropic functionally graded solid circular plate with different geometric and material parameters.     

Three-dimensional elasticity solution for bending of transversely isotropic functionally graded plates

This paper presents a three-dimensional elasticity solution for a simply supported, transversely isotropic functionally graded plate subjected to transverse loading, with Young’s moduli and the shear modulus varying exponentially through the thickness and Poisson’s ratios being constant. The approach makes use of the recently developed displacement functions for inhomogeneous transversely isotropic media. Dependence of stress and displacement fields in the plate on the inhomogeneity ratio, geometry and degree of anisotropy is examined and discussed. The developed three-dimensional solution for transversely isotropic functionally graded plate is validated through comparison with the available three-dimensional solutions for isotropic functionally graded plates, as well as the classical and higher-order plate theories.     

Revisiting displacement functions in three-dimensional elasticity of inhomogeneous media

The paper presents a three-dimensional solution to the equilibrium equations for linear elastic transversely isotropic inhomogeneous media. We assume that the material has constant Poisson’s ratios, and its Young’s and shear moduli have the same functional form of dependence on the co-ordinate normal to the plane of isotropy. We show, apparently for the first time, that stresses and displacements in such an inhomogeneous transversely isotropic elastic solid can be represented in terms of two displacement functions which satisfy the second- and fourth-order partial differential equations. We examine and discuss key aspects of the new representation; they include the relationship between the new displacement functions and Plevako’s solution for isotropic inhomogeneous material with constant Poisson’s ratio as well as the application of the new representation to some important classes of three-dimensional elasticity problems. As an example, the displacement function is derived that can be used to determine stresses and displacements in an inhomogeneous transversely isotropic half-space which is subjected to a concentrated force normal to a free surface and applied at the origin (Boussinesq’s problem).     

动态载荷下功能梯度复合材料的圆币形裂纹问题

研究了动态载荷下功能梯度材料中的圆币形裂纹问题．假设材料为横观各向同性,并且含有多个垂直于厚度方向的裂纹,材料参数沿轴向（与裂纹面垂直的方向）为变化的,沿该方向将材料划分为许多单层,各单层材料参数为常数,利用Hankel变换祛,在Laplace域内推导出了控制问题的对偶积分方程组．利用Laplace数值反演,得出了裂纹尖端的动态应力强度因子和能量释放率．研究了含两个裂纹的功能梯度接头结构,分析了材料非均匀性参数对应力强度因子和能量释放率的影响．     

Free axisymmetric vibration of FGM circular plates

Based on the three-dimensional theory, this work presents a direct displacement method to investigate the free axisymmetric vibration of transversely isotropic circular plates, whose material is functionally graded and properties obey the exponential law along the thickness direction of the plate. Under two boundary conditions, the solution satisfies all basic equations and the Corresponding boundary condition at every point. Thus, it is three-dimensional exact. Numerical examples are presented and compared with previous works. The present method can also be extended to the case of arbitrary distribution of the material properties along the thickness direction of the plate.     

Analytical solutions for plane problem of functionally graded magnetoelectric cantilever beam

In this paper, an exact analytical solution is presented for a transversely isotropic functionally graded magneto-electro-elastic (FGMEE) cantilever beam, which is subjected to a uniform load on its upper surface, as well as the concentrated force and moment at the free end. This solution can be applied for any form of gradient distribution. For the basic equations of plane problem, all the partial differential equations governing the stress field, electric, and magnetic potentials are derived. Then, the expressions of Airy stress, electric, and magnetic potential functions are assumed as quadratic polynomials of the longitudinal coordinate. Based on all the boundary conditions, the exact expressions of the three functions can be determined. As numerical examples, the material parameters are set as exponential and linear distributions in the thickness direction. The effects of the material parameters on the mechanical, electric, and magnetic fields of the cantilever beam are analyzed in detail.     

Elastic and viscoelastic solutions to rotating functionally graded hollow and solid cylinders

Analytical solutions to rotating functionally graded hollow and solid long cylinders are developed. Young＇s modulus and material density of the cylinder are assumed to vary exponentially in the radial direction, and Poisson＇s ratio is assumed to be constant. A unified governing equation is derived from the equilibrium equations, compatibility equation, deformation theory of elasticity and the stress-strain relationship. The governing second-order differential equation is solved in terms of a hypergeometric function for the elastic deformation of rotating functionally graded cylinders. Dependence of stresses in the cylinder on the inhomogeneous parameters, geometry and boundary conditions is examined and discussed. The proposed solution is validated by comparing the results for rotating functionally graded hollow and solid cylinders with the results for rotating homogeneous isotropic cylinders. In addition, a viscoelastic solution to the rotating viscoelastic cylinder is presented, and dependence of stresses in hollow and solid cylinders on the time parameter is examined.     

Several three-dimensional solutions for transversely isotropic functionally graded material plate welded with circular inclusion

The problem of a transversely isotropic functionally graded material (FGM) plate welded with a circular inclusion is considered. The analysis starts with the generalized England-Spencer plate theory for transversely isotropic FGM plates, which expresses a three-dimensional (3D) general solution in terms of four analytic functions. Several analytical solutions are then obtained for an infinite FGM plate welded with a circular inclusion and subjected to the loads at infinity. Three different cases are considered, i.e., a rigid circular inclusion fixed in the space, a rigid circular inclusion rotating about the x-, y-, and z-axes, and an elastic circular inclusion with different material constants from the plate itself. The static responses of the plate and/or the inclusion are investigated through numerical examples.     

Love solutions in the linear inhomogeneous transversely isotropic theory of elasticity

A general Love solution for the inhomogeneous transversely isotropic theory of elasticity with the elastic constants dependent on the coordinate z is proposed. This result may be considered as a generalization of the Love solutions we recently derived for the inhomogeneous isotropic theory of elasticity. The key steps of deriving the Love solution for the classical linear homogeneous transversely isotropic theory of elasticity are described for further use of the derivation procedure, which is then generalized to the inhomogeneous transversely isotropic case. Some particular cases of inhomogeneity traditionally used in the theory of elasticity are also examined. The significance of the derived solutions and their importance for the modeling of functionally graded materials are briefly discussed     

Green’s functions for transversely isotropic piezoelectric functionally graded multilayered half spaces

By virtue of two systems of vector functions and the propagator matrix method, Green’s functions for transversely isotropic, piezoelectric functionally graded (exponentially graded in the vertical direction), and multilayered half spaces are derived. It is observed that the homogeneous solution and propagator matrices for each functionally graded layer in the transformed domain are independent of the choice of the two systems of vector functions. For a point force and point charge density applied at any location of the functionally graded half space, Green’s functions are expressed in terms of one-dimensional infinite integrals. To carry out the numerical integral involving Bessel functions, an adaptive Gauss quadrature approach is introduced and modified. The piezoelectric functionally graded Green’s functions include those in the corresponding elastic functionally graded media as special results with the latter being also unavailable in the literature. Two piezoelectric functionally graded half-space models are analyzed numerically: one is a functionally graded PZT-4 half space, and the other a coated functionally graded PZT-4 layer over a homogeneous BaTiO₃ half space. The effects of different exponential factors on Green’s function components are clearly demonstrated, which could be useful in the design and manufacturing of piezoelectric functionally graded structures.     

Penny-shaped crack in transversely isotropic piezoelectric materials

Using a method of potential functions introduced successively to integrate the field equations of three-dimensional problems for transversely isotropic piezoelectric materials, we obtain the so-called general solution in which the displacement components and electric potential functions are represented by a singular function satisfying some special partial differential equations of 6th order. In order to analyse the mechanical-electric coupling behaviour of penny-shaped crack for above materials, another form of the general solution is obtained under cylindrical coordinate system by introducing three quasi-harmonic functions into the general equations obtained above. It is shown that both the two forms of the general solutions are complete. Furthermore, the mechanical-electric coupling behaviour of penny-shaped crack in transversely isotropic piezoelectric media is analysed under axisymmetric tensile loading case, and the crack-tip stress field and electric displacement field are obtained. The results show that the stress and the electric displacement components near the crack tip have (r ^−1/2) singularity. The project supported by the Natural Science Foundation of Shaanxi Province, China     

粘贴压电层功能梯度材料Timoshenko梁的热过屈曲分析

研究了上下表面粘贴压电层的功能梯度材料Timoshenko梁在升温及电场作用下的过屈曲行为。在精确考虑轴线伸长和一阶横向剪切变形的基础上,建立了压电功能梯度Timoshenko层合梁在热-电-机械载荷作用下的几何非线性控制方程。其中,假设功能梯度的材料性质沿厚度方向按照幂函数连续变化,压电层为各向同性均匀材料。采用打靶法数值求解所得强非线性边值问题,获得了在均匀电场和横向非均匀升温场内两端固定Timoshenko梁的静态非线性屈曲和过屈曲数值解。并给出了梁的变形随热、电载荷及材料梯度参数变化的特性曲线。结果表明,通过施加电压在压电层产生拉应力可以有效地提高梁的热屈曲临界载荷,延缓热过屈曲发生。由于材料在横向的非均匀性,即使在均匀升温和均匀电场作用下,也会产生拉-弯耦合效应。但是对于两端固定的压电-功能梯度材料梁,在横向非均匀升温下过屈曲变形仍然是分叉形的。     

Three-dimensional elastostatic solutions for transversely isotropic functionally graded material plates containing elastic inclusion

Based on the generalized England-Spencer plate theory, the equilibrium of a transversely isotropic functionally graded plate containing an elastic inclusion is studied. The general solutions of the governing equations are expressed by four analytic functions α(ζ), β(ζ), φ(ζ), and ψ(ζ) when no transverse forces are acting on the surfaces of the plate. Axisymmetric problems of a functionally graded circular plate and an infinite func-tionally graded plate containing a circular hole subject to loads applied on the cylindrical boundaries of the plate are firstly investigated. On this basis, the three-dimensional (3D) elasticity solutions are then obtained for a functionally graded infinite plate containing an elastic circular inclusion. When the material is degenerated into the homogeneous one, the present elasticity solutions are exactly the same as the ones obtained based on the plane stress elasticity, thus validating the present analysis in a certain sense.     

均布荷载作用下功能梯度简支梁弯曲的解析解

利用应力函数半逆解法,研究了均布载荷作用下、材料属性在厚度上任意变化的功能梯度简支梁弯曲的解析解,给出了各向应力应变与位移的解析显式表达式.首先根据平面应力状态的基本方程,得出了功能梯度梁的应力函数应满足的偏微分方程,并根据应力边界条件得出了各应力分布的表达式;进而根据功能梯度材料的本构方程和位移边界条件,得出了应变和位移的分布.最后,通过将本文的解退化到均质各向同性梁并与经典弹性解比较,证明了本文理论的正确性,并求解了材料组分呈幂律分布的功能梯度梁的应力和位移分布,分析了上下表层材料的弹性模量比λ与组分材料体积分数指数n对应力和位移分布的影响.     

Static analysis of simply supported functionally graded and layered magneto-electro-elastic plates

In this article, static analysis of functionally graded, anisotropic and linear magneto-electro-elastic plates have been carried out by semi-analytical finite element method. A series solution is assumed in the plane of the plate and finite element procedure is adopted across the thickness of the plate such a way that the three-dimensional character of the solution is preserved. The finite element model is derived based on constitutive equation of piezomagnetic material accounting for coupling between elasticity, electric and magnetic effect. The present finite element is modeled with displacement components, electric potential and magnetic potential as nodal degree of freedom. The other fields are calculated by post-computation through constitutive equation. The functionally graded material is assumed to be exponential in the thickness direction. The numerical results obtained by the present model are in good agreement with available functionally graded three-dimensional exact benchmark solutions given by Pan and Han [Pan, E., Han, F., in press. Green’s function for transversely isotropic piezoelectric functionally graded multilayered half spaces. Int. J. Solids Struct.]. Numerical study includes the influence of the different exponential factor, magneto-electro-elastic properties and effect of mechanical and electric type of loading on induced magneto-electro-elastic fields. In addition further study has been carried out on non-homogeneous transversely isotropic FGM magneto-electro-elastic plate available in the literature [Chen, W.Q., Lee, K.Y., Ding, H.J., 2005. On free vibration of non-homogeneous transversely isotropic magneto-electro-elastic plates].     

An exact solution for the steady-state thermoelastic response of functionally graded orthotropic cylindrical shells

We analyze the steady-state response of a functionally graded thick cylindrical shell subjected to thermal and mechanical loads. The functionally graded shell is simply supported at the edges and it is assumed to have an arbitrary variation of material properties in the radial direction. The three-dimensional steady-state heat conduction and thermoelasticity equations, simplified to the case of generalized plane strain deformations in the axial direction, are solved analytically. Suitable temperature and displacement functions that identically satisfy the boundary conditions at the simply supported edges are used to reduce the thermoelastic equilibrium equations to a set of coupled ordinary differential equations with variable coefficients, which are then solved by the power series method. In the present formulation, the cylindrical shell is assumed to be made of an orthotropic material, although the analytical solution is also valid for isotropic materials. Results are presented for two-constituent isotropic and fiber-reinforced functionally graded shells that have a smooth variation of material volume fractions, and/or in-plane fiber orientations, through the radial direction. The cylindrical shells are also analyzed using the Flügge and the Donnell shell theories. Displacements and stresses from the shell theories are compared with the three-dimensional exact solution to delineate the effects of transverse shear deformation, shell thickness and angular span.     

功能梯度压电、压磁柱对称问题的振动分析

本文在不考虑体力、体电流和体电荷的情况下, 假定压电、压磁柱壳的材料参数沿圆柱厚度方向呈幂函数分布时, 研究了径向载荷作用下功能梯度压电、压磁空心柱壳的空间柱对称径向振动问题. 首先在柱坐标系下, 由功能梯度材料的参数、本构、梯度和平衡方程推导得出外激励作用下以Bessel函数表示圆柱壳的应力、电势、磁势等物理量的稳态解, 进而对空间柱对称的压电、压磁功能梯度材料动力控制问题进行了理论分析. 可以看出, 当梯度参数时, 即完全退化为横观各项同性压电、压磁柱对称的振动问题, 与文献[16]的基本方程为柱坐标下得出的结果一致. 最后给出数值算例, 结果表明材料不均匀性对沿径向振动各物理量的显著影响, 且用一个特定不均匀性参数值可以优化电磁力耦合的性能, 这在现代工程设计中尤为重要.     

Exact Lévy-Type Solutions for Plate Bending Exist for Transversely Isotropic but Not for General Monoclinic Materials

An exact three-dimensional Lévy-type solution for the bending of an elastic slab is one in which there are edge loads only and the unknown displacement and stresses have very simple polynomial dependence on the thickness coordinate. Lévy obtained such a solution in 1877 for a linearly elastic, isotropic, plate-like body. The most general material that allows bending and stretching to be uncoupled is monoclinic (13 elastic constants). However, it is shown that only transversely isotropic materials (5 elastic constants) admit exact solutions having polynomial dependence in the thickness direction. Such solutions are listed explicitly.