A comprehensive analysis of functionally graded sandwich plates: Part 1—Deflection and stresses

A two-dimensional solution is presented for bending analysis of simply supported functionally graded ceramic–metal sandwich plates. The sandwich plate faces are assumed to have isotropic, two-constituent material distribution through the thickness, and the modulus of elasticity and Poisson’s ratio of the faces are assumed to vary according to a power-law distribution in terms of the volume fractions of the constituents. The core layer is still homogeneous and made of an isotropic ceramic material. Several kinds of sandwich plates are used taking into account the symmetry of the plate and the thickness of each layer. We derive field equations for functionally graded sandwich plates whose deformations are governed by either the shear deformation theories or the classical theory. Displacement functions that identically satisfy boundary conditions are used to reduce the governing equations to a set of coupled ordinary differential equations with variable coefficients. Numerical results of the sinusoidal, third-order, first-order and classical theories are presented to show the effect of material distribution on the deflections and stresses.     

Dynamic fracture behaviors of cracks in a functionally graded magneto-electro-elastic plate

In this paper the dynamic anti-plane problem for a functionally graded magneto-electro-elastic plate containing an internal or an edge crack parallel to the graded direction is investigated. The crack is assumed to be magneto-electrically impermeable. Integral transforms and dislocation density functions are employed to reduce the problem to Cauchy singular integral equations. Field intensity factors and energy release rate are derived, analyzed and partially calculated numerically. The effects of material graded index, loading combination parameter (including size and direction) and geometry criterion of the plate on the dynamic energy release rate are shown graphically. Numerical results indicate that increasing the graded index can all retard the crack extension, and that both the applied magnetic field loadings and electric field loadings play a dominant role in the dynamic fracture behaviors of crack tips.     

Postbuckling of FGM plates with piezoelectric actuators under thermo-electro-mechanical loadings

A postbuckling analysis is presented for a simply supported, shear deformable functionally graded plate with piezoelectric actuators subjected to the combined action of mechanical, electrical and thermal loads. The temperature field considered is assumed to be of uniform distribution over the plate surface and through the plate thickness and the electric field considered only has non-zero-valued component E_Z. The material properties of functionally graded materials (FGMs) are assumed to be graded in the thickness direction according to a simple power law distribution in terms of the volume fractions of the constituents, and the material properties of both FGM and piezoelectric layers are assumed to be temperature-dependent. The governing equations are based on a higher order shear deformation plate theory that includes thermo-piezoelectric effects. The initial geometric imperfection of the plate is taken into account. Two cases of the in-plane boundary conditions are considered. A two step perturbation technique is employed to determine buckling loads and postbuckling equilibrium paths. The numerical illustrations concern the postbuckling behavior of perfect and imperfect, geometrically mid-plane symmetric FGM plates with fully covered or embedded piezoelectric actuators under different sets of thermal and electric loading conditions. The effects played by temperature rise, volume fraction distribution, applied voltage, the character of in-plane boundary conditions, as well as initial geometric imperfections are studied.     

Two new refined shear displacement models for functionally graded sandwich plates

Two refined displacement models, RSDT1 and RSDT2, are developed for a bending analysis of functionally graded sandwich plates. Unlike any other theory, the number of unknown functions involved is only four, as against five in case of other shear deformation theories. The developed models are variationally consistent, have strong similarity with classical plate theory in many aspects, do not require shear correction factor, and give rise to transverse shear stress variation such that the transverse shear stresses vary parabolically across the thickness satisfying shear stress-free surface conditions. The accuracy of the analysis presented is demonstrated by comparing the results with solutions derived from other higher-order models. The functionally graded layers are assumed to have isotropic, two-constituent material distribution through the thickness, and the modulus of elasticity, Poisson’s ratio of the faces, and thermal expansion coefficients are assumed to vary according to a power-law distribution in terms of the volume fractions of the constituents. The core layer is still homogeneous and made of an isotropic ceramic material. Numerical results for deflections and stresses of functionally graded metal–ceramic plates are investigated. It can be concluded that the proposed models are accurate and simple in solving the bending behavior of functionally graded plates.     

Analytical solutions for single- and multi-span functionally graded plates in cylindrical bending

A recently developed plate theory using the concept of shape function of the transverse coordinate parameter is extended to determine the stress distribution in an orthotropic functionally graded plate subjected to cylindrical bending. The transfer matrix method is presented to derive the shape function. The equations governing the plate deformation are then solved analytically using the transfer matrix method for arbitrary boundary conditions. For a simply supported functionally graded plate, a comparison of the present solution with the exact elasticity solution, the first- and third-order shear deformation plate theories is presented and discussed. It is demonstrated that the present method yields more accurate stresses than the first- and third-order shear deformation theories. The effect of boundary conditions and inhomogeneity of material on the displacements and stresses in functionally graded plates are investigated. A multi-span functionally graded plate with arbitrary boundary conditions is further considered to demonstrate the efficiency of the present method.     

Nonlinear dynamic response of a functionally graded plate with a through-width surface crack

A nonlinear vibration analysis of a simply supported functionally graded rectangular plate with a through-width surface crack is presented in this paper. The plate is subjected to a transverse excitation force. Material properties are graded in the thickness direction according to exponential distributions. The cracked plate is treated as an assembly of two sub-plates connected by a rotational spring at the cracked section whose stiffness is calculated through stress intensity factor. Based on Reddy’s third-order shear deformation plate theory, the nonlinear governing equations of motion for the FGM plate are derived by using the Hamilton’s principle. The deflection of each sub-plate is assumed to be a combination of the first two mode shape functions with unknown constants to be determined from boundary and compatibility conditions. The Galerkin’s method is then utilized to convert the governing equations to a two-degree-of-freedom nonlinear system including quadratic and cubic nonlinear terms under the external excitation, which is numerically solved to obtain the nonlinear responses of cracked FGM rectangular plates. The influences of material property gradient, crack depth, crack location and plate thickness ratio on the vibration frequencies and transient response of the surface-racked FGM plate are discussed in detail through a parametric study.     

Higher-order crack tip fields for functionally graded material plate with transverse shear deformation

The crack tip fields are investigated for a cracked functionally graded material (FGM) plate by Reissner’s linear plate theory with the consideration of the transverse shear deformation generated by bending. The elastic modulus and Poisson’s ratio of the functionally graded plates are assumed to vary continuously through the coordinate y, according to a linear law and a constant, respectively. The governing equations, i.e., the 6th-order partial differential equations with variable coefficients, are derived in the polar coordinate system based on Reissner’s plate theory. Furthermore, the generalized displacements are treated in a separation-of-variable form, and the higher-order crack tip fields of the cracked FGM plate are obtained by the eigen-expansion method. It is found that the analytic solutions degenerate to the corresponding fields of the isotropic homogeneous plate with Reissner’s effect when the in-homogeneity parameter approaches zero.     

Thermoelastic buckling behavior of thick functionally graded rectangular plates

Thermoelastic buckling behavior of thick rectangular plate made of functionally graded materials is investigated in this article. The material properties of the plate are assumed to vary continuously through the thickness of the plate according to a power-law distribution. Three types of thermal loading as uniform temperature raise, nonlinear and linear temperature distribution through the thickness of plate are considered. The coupled governing stability equations are derived based on the Reddy’s higher-order shear deformation plate theory using the energy method. The resulted stability equations are decoupled and solved analytically for the functionally graded rectangular plates with two opposite edges simply supported subjected to different types of thermal loading. A comparison of the present results with those available in the literature is carried out to establish the accuracy of the presented analytical method. The influences of power of functionally graded material, plate thickness, aspect ratio, thermal loading conditions and boundary conditions on the critical buckling temperature of aluminum/alumina functionally graded rectangular plates are investigated and discussed in detail. The critical buckling temperatures of thick functionally graded rectangular plates with various boundary conditions are reported for the first time and can be served as benchmark results for researchers to validate their numerical and analytical methods in the future.     

Effects of multiple cracking on the residual strength behavior of thermally shocked functionally graded ceramics

When subjected to severe thermal shocks a functionally graded ceramic (FGC) suffers strength degradation due to the thermally-induced damages in the material. Multiple surface cracking has been observed as one of the dominant defects/damages affecting the thermal shock behavior of ceramics. This paper presents a thermo-fracture mechanics model to investigate the thermal shock residual strength behavior of elastically homogeneous but thermally graded FGCs undergoing multiple surface cracking. We consider an FGC plate with an array of parallel edge cracks at the thermally shocked surface. A Fourier transform/superposition method is used to derive the singular integral equation of the thermal shock crack problem. The critical thermal shock that causes crack propagation and thermal shock damage are determined using linear elastic fracture mechanics. The thermal shock residual strength of the FGC as a function of thermal shock severity and crack density (crack spacing) is subsequently evaluated. Numerical calculations are carried out for two FGC materials, i.e., Al₂O₃/Si₃N₄ and TiC/SiC FGCs, to illustrate the effects of crack density (crack spacing) and material gradation on the thermal shock strength behavior of FGCs. It is found that a higher crack density (lower crack spacing) together with appropriately graded material properties significantly enhances the residual strength of the thermally shocked FGCs.     

黏弹性功能梯度材料裂纹问题的有限元方法

针对组分材料体积含量任意分布的黏弹性功能梯度材料裂纹问题建立有限元分析途径. 通过Laplace变换,将黏弹性问题转化到象空间中求解,基于反映材料非均匀的梯度单元和裂纹尖端奇异特性的奇异单元计算象空间中的位移、应力和应变场,应用虚拟裂纹闭合方法得到应变能释放率,分别由应力和应变能释放率确定应力强度因子. 给出这些断裂参量在物理空间和象空间之间的对应关系,由数值逆变换求出其在物理空间的相应值. 文中分析两端均匀受拉的黏弹性边裂纹板条,首先针对松弛模量表示为空间函数和时间函数乘积的特殊梯度材料进行计算,结合对应原理验证方法的有效性. 然后分析组分材料体积含量具有任意梯度分布的情形,由Mori-Tanaka方法预测象空间中的等效松弛模量. 计算结果表明,蠕变加载条件下,应变能释放率随时间增加,其增大程度与黏弹性组分材料体积含量相关. 由于梯度材料的非均匀黏弹性性质,产生应力重新分布,导致应力强度因子随时间变化,其变化范围与组分材料的体积含量分布方式有关.     

热环境中旋转运动功能梯度圆板的强非线性固有振动

研究热环境中旋转运动功能梯度圆板的非线性固有振动问题．针对金属－陶瓷功能梯度圆板，考虑几何非线性、材料物理属性参数随温度变化以及材料组分沿厚度方向按幂律分布的情况，应用哈密顿原理推得热环境中旋转运动功能梯度圆板的非线性振动微分方程．考虑周边夹支边界条件，利用伽辽金法得到了横向非线性固有振动方程，并确定了静载荷引起的静挠度．用改进的多尺度法求解强非线性方程，得出非线性固有频率表达式．通过算例，分析了旋转运动功能梯度圆板固有频率随转速、温度等参量的变化情况．结果表明，非线性固有频率随金属含量的增加而降低；随转速和圆板厚度的增大而升高；随功能梯度圆板表面温度的升高而降低．     

热环境中功能梯度圆形薄板的混沌运动

针对陶瓷-金属功能梯度圆板,同时考虑几何非线性、材料物性参数随温度变化且材料组分沿厚度方向按幂律分布的情况,应用虚功原理给出了热载荷与横向简谐载荷共同作用下的非线性振动偏微分方程。在固支无滑动的边界条件下,通过引入位移函数,利用伽辽金方法得到了达芬型非线性动力学方程。利用Melnikov方法,给出了热环境中功能梯度圆板可能发生混沌运动的临界条件。通过数值算例,给出了不同体积分数指数和温度的同宿分岔曲线,平面相图和庞加莱映射图,讨论其对临界条件的影响,证实了系统混沌运动的存在。通过分岔图和与其相对应的最大李雅普诺夫指数图,分析了激励频率和激励幅值对倍周期分岔的影响及变化规律,发现系统可出现周期、倍周期和混沌等复杂动力学响应。     

功能梯度材料板中Lamb波传播特性研究

对材料性能参数沿厚度连续变化的横观各向同性热应力缓和型功能梯度材料板中Lamb波的传播问题,采用幂级数法,求得其相速度方程.借助数值算例,分析了参数梯度变化对Lamb波频散曲线的影响,并与相应陶瓷板和金属板中的频散曲线进行了对比.进一步研究了参数梯度变化对波结构的影响,揭示了Lamb波在这种非均质板中的传播行为,所得结果可以为功能梯度材料及结构的超声表征与检测提供理论依据.     

Thermal buckling analysis of functionally graded beam with longitudinal crack

The thermal buckling problem of functionally graded beam with longitudinal crack is presented in the paper. The whole beam is divided into four sub-beams and each one is modeled as a Timoshenko beam. The buckling governing equation of each sub-beam in thermal environment is established by using Hamilton Principle. Combining with the boundary conditions, the continuous conditions of the displacements and the forces, the buckling governing equations are solved by both the analytical and numerical methods. The buckling modes and critical buckling temperatures are obtained, and the effects of the functionally graded index, crack length, crack depth, and crack longitudinal location on the buckling characteristics of beams are discussed in numerical examples.     

弹性地基加肋功能梯度板自由振动分析的无网格法

基于一阶剪切变形理论和移动最小二乘近似研究Winkler弹性地基上加肋功能梯度板的固有频率。假设功能梯度板的材料性质沿厚度方向按幂函数连续变化,基于物理中面和移动最小二乘近似分别推导功能梯度板和肋条的动能和势能,再通过引入位移协调条件,建立板和肋条节点参数转换关系,叠加两者的总能量,然后利用Hamilton原理推导加肋功能梯度板自由振动控制方程。采用完全转换法施加边界条件。通过将本文的计算结果与有限元以及文献的结果对比,验证方法的收敛性以及准确性。     

Stability analysis of functionally graded rectangular plates under nonlinearly varying in-plane loading resting on elastic foundation

In this research work, an exact analytical solution for buckling of functionally graded rectangular plates subjected to non-uniformly distributed in-plane loading acting on two opposite simply supported edges is developed. It is assumed that the plate rests on two-parameter elastic foundation and its material properties vary through the thickness of the plate as a power function. The neutral surface position for such plate is determined, and the classical plate theory based on exact neutral surface position is employed to derive the governing stability equations. Considering Levy-type solution, the buckling equation reduces to an ordinary differential equation with variable coefficients. An exact analytical solution is obtained for this equation in the form of power series using the method of Frobenius. By considering sufficient terms in power series, the critical buckling load of functionally graded plate with different boundary conditions is determined. The accuracy of presented results is verified by appropriate convergence study, and the results are checked with those available in related literature. Furthermore, the effects of power of functionally graded material, aspect ratio, foundation stiffness coefficients and in-plane loading configuration together with different combinations of boundary conditions on the critical buckling load of functionally graded rectangular thin plate are studied.     

Nonlinear bending and post-buckling of a functionally graded circular plate under mechanical and thermal loadings

Based on the classical nonlinear von Karman plate theory, axisymmetric large deflection bending of a functionally graded circular plate is investigated under mechanical, thermal and combined thermal–mechanical loadings, respectively, and axisymmetric thermal post-buckling behavior of a functionally graded circular plate is also investigated. The mechanical and thermal properties of functionally graded material (FGM) are assumed to vary continuously through the thickness of the plate, and obey a simple power law of the volume fraction of the constituents. Governing equations for the problem are derived, and then a shooting method is employed to numerically solve the equations. Effects of material constant n and boundary conditions on the temperature distribution, nonlinear bending, critical buckling temperature and thermal post-buckling behavior of the FGM plate are discussed in details.     

Non-linear analysis of functionally graded fiber reinforced composite laminated plates,Part II: Numerical results

In this Part, the extensive parametric studies performed are reported and numerical results are presented for the non-linear vibration, non-linear bending and compressive postbuckling of uniformly distributed and functionally graded fiber reinforced unsymmetric cross-ply and/or antisymmetric angle-ply laminated plates resting on Pasternak elastic foundations under different hygrothermal environmental conditions. The numerical results show that the functionally graded fiber reinforcement has a significant effect on the postbuckling response and load-bending moment curves of plate bending, whereas this effect is less pronounced on the load-deflection curves of plate bending and the linear and non-linear frequencies of the same plate.     

平行于功能梯度材料夹层的币型裂纹起裂条件

分析了功能梯度材料中币型裂纹的扩展问题．裂纹平行于无限域中功能梯度材料夹层,受有与裂纹面成任意角度的拉应力．假定功能梯度材料夹层与两个半无限域均匀介质完全粘合,其弹性模量沿厚度方向变化．采用基于层状材料广义Kelvin基本解的边界元方法分析裂纹问题,给出了均布正应力和剪应力作用下裂纹的应力强度因子、将应力强度因子耦合于应变能密度断裂判据,讨论了裂纹体在拉伸应力作用下的起裂条件．     

面内预应力功能梯度材料矩形板的横向弯曲

功能梯度材料（FGM)是一类具有广阔应用前景的新型复合材料。本文考虑材料物性参数随坐标和温度变化的特性，研究横向荷载和面内预加荷载作用下FGM矩形板在各种边界条件下的弯曲问题。给出了基于一维微分求积格式的Galerkin技术的半解析方法，并以ZrO2/Ti-6Al-4V板为例考察了材料组分，温度相关性，面内预加荷载，边界约束条件等对FGM板弯曲行为的影响。结果表明，FGM板的弯曲变形介于各向同性陶瓷板与各向同性金属板之间，且随板抗弯刚度的增大而逐步减小，在高温下条件下必须考虑材料物性和温度的相关性。