RMVT- and PVD-based finite cylindrical layer methods for the three-dimensional buckling analysis of multilayered FGM cylinders under axial compression

The unified formulations of finite cylindrical layer methods (FCLMs) based on the Reissner mixed variational theorem (RMVT) and the principle of virtual displacements (PVD) are developed for the three-dimensional (3D) linear buckling analysis of simply-supported, multilayered functionally graded material (FGM) circular hollow cylinders and laminated composite ones under axial compression. The material properties of the FGM layer are assumed to obey the power-law distributions of the volume fraction of the constituents through the thickness coordinate. In these formulations, the cylinder is divided into a number of finite cylindrical layers, in which the trigonometric functions and Lagrange polynomials are used to interpolate the in- and out-of-surface variations of the primary variables of each individual layer, respectively, as well as the related order of each primary variable can be freely chosen, such as the layerwise linear, quadratic or cubic function distribution through the thickness coordinate. The accuracy and convergence of the RMVT- and PVD-based FCLMs developed in this article are assessed by comparing their solutions with the exact 3D solutions available in the literature.     

Approximation of functionally graded plates with non-conforming finite elements

In this paper rectangular plates made of functionally graded materials (FGMs) are studied. A two-constituent material distribution through the thickness is considered, varying with a simple power rule of mixture. The equations governing the FGM plates are determined using a variational formulation arising from the Reissner–Mindlin theory. To approximate the problem a simple locking-free Discontinuous Galerkin finite element of non-conforming type is used, choosing a piecewise linear non-conforming approximation for both rotations and transversal displacement. Several numerical simulations are carried out in order to show the capability of the proposed element to capture the properties of plates of various gradings, subjected to thermo-mechanical loads.     

Functionally graded cylinders : Asymptotically exact analytical formulations

This paper uses Variational Asymptotic Method (VAM) to obtain asymptotically exact analytical solutions for cylinders made up of Functionally Graded Materials (FGM). VAM splits a three-dimensional elasticity problem into a two-dimensional linear cross-sectional problem and a one-dimensional beam problem. It ensures the asymptotic correctness since the method makes no ad hoc assumptions. This is accomplished by taking advantage of certain small parameters inherent to beam-like structures. This technique has been successfully used for a variety of problems but it has never been implemented on cylinders made up of FGM. Starting with the variable geometry of the cylindrical cross-section made up of radially non-homogeneous material properties controlled by a volume fraction variable, the 3-D problem has been formulated and solved analytically despite the presence of geometrical non-linearities. Closed form analytical solutions are obtained to predict the exact response of functionally graded cylinders. The influence of inner and outer material properties, radius values and the variation trend of material composition on the mechanical behavior is highlighted. Results obtained from the present theory are successfully validated using a commercially available 3-D FEM solver ABAQUS. Analytical solutions obtained can analyze the behavior of any cylinder with FGM quickly and accurately.     

Nonlinear dynamic response for functionally graded shallow spherical shell under low velocity impact in thermal environment

Based on Giannakopoulos’s 2-D functionally graded material (FGM) contact model, a modified contact model is put forward to deal with impact problem of the functionally graded shallow spherical shell in thermal environment. The FGM shallow spherical shell, having temperature dependent material property, is subjected to a temperature field uniform over the shell surface but varying along the thickness direction due to steady-state heat conduction. The displacement field and geometrical relations of the FGM shallow spherical shell are established on the basis of Timoshenko–Midlin theory. And the nonlinear motion equations of the FGM shallow spherical shell under low velocity impact in thermal environment are founded in terms of displacement variable functions. Using the orthogonal collocation point method and the Newmark method to discretize the unknown variable functions in space and in time domain, the whole problem is solved by the iterative method. In numerical examples, the contact force and nonlinear dynamic response of the FGM shallow spherical shell under low velocity impact are investigated and effects of temperature field, material and geometrical parameters on contact force and dynamic response of the FGM shallow spherical shell are discussed.     

Variational formulation of the material failure process in solids by embedded discontinuities model

This work presents a variational formulation of the material failure process, idealized as strain or displacement discontinuities, by weak, strong, or discrete embedded discontinuities into a continuum. It is shown that the solution of the proposed variational formulation may be approximated by different types of finite elements with embedded discontinuities. The developed displacement approximation of a finite element split by the discontinuity leads to a symmetric stiffness matrix, which considers not only the continuity of tractions but also the rigid body relative motions of the portions in which the element is split. The variational formulation of a continuum with more than one discontinuity in its interior is developed. It is shown that this formulation may lead to finite elements with embedded discontinuities that can be classified as displacement, force, mixed, and hybrid models. To show the effectiveness of the proposed formulation, the classical example of a bar under tension is solved using one and 2D finite element approximations. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Post buckling response of functionally graded materials plate subjected to mechanical and thermal loadings with random material properties

In this paper, the second order statistics of post buckling response of functionally graded materials plate (FGM) subjected to mechanical and thermal loading with nonuniform temperature changes subjected to temperature independent (TID) and dependent (TD) material properties is examined. Material properties such as material properties of each constituent’s materials, volume fraction index are taken as independent random input variables. The basic formulation is based on higher order shear deformation theory (HSDT) with von-Karman nonlinear kinematic using modified C⁰ continuity. A direct iterative based C⁰ nonlinear finite element method (FEM) combined with mean centered first order perturbation technique (FOPT) proposed by last two authors for the composite plate is extended for Functionally Graded Materials (FGMs) plate with reasonable accuracy to compute the second order statistics (mean and coefficient of variation) of the post buckling load response of the FGM plates. The effect of random material properties with amplitude ratios, volume fraction index, plate thickness ratios, aspect ratios, boundary conditions and types of loadings subjected to TID and TD material properties are presented through numerical examples. The performance of outlined present approach is validated with the results available in literatures and independent Monte Carlo simulation (MCS).     

Static response and free vibration analysis of FGM plates using higher order shear deformation theory

Free vibration and static analysis of functionally graded material (FGM) plates are studied using higher order shear deformation theory with a special modification in the transverse displacement in conjunction with finite element models. The mechanical properties of the plate are assumed to vary continuously in the thickness direction by a simple power-law distribution in terms of the volume fractions of the constituents. The fundamental equations for FGM plates are derived using variational approach by considering traction free boundary conditions on the top and bottom faces of the plate. Results have been obtained by employing a continuous isoparametric Lagrangian finite element with 13 degrees of freedom per node. Convergence tests and comparison studies have been carried out to demonstrate the efficiency of the present model. Numerical results for different thickness ratios, aspect ratios and volume fraction index with different boundary conditions have been presented. It is observed that the natural frequency parameter increases for plate aspect ratio, lower volume fraction index n and smaller thickness ratios. It is also observed that the effect of thickness ratio on the frequency of a plate is independent of the volume fraction index. For a given thickness ratio non-dimensional deflection increases as the volume fraction index increases. It is concluded that the gradient in the material properties plays a vital role in determining the response of the FGM plates.     

Exact solutions of functionally graded piezoelectric material sandwich cylinders by a modified Pagano method

The three-dimensional (3D) coupled analysis of simply-supported, functionally graded piezoelectric material (FGPM) circular hollow sandwich cylinders under electro-mechanical loads is presented. The material properties of each FGPM layer are regarded as heterogeneous through the thickness coordinate, and obey an exponent-law dependent on this. The Pagano method is modified to be feasible for the study of FGPM sandwich cylinders. The modifications are as follows: a displacement-based formulation is replaced by a mixed formulation; a set of the complex-valued solutions of the system equations is transferred to the corresponding set of real-valued solutions; a successive approximation method is adopted to approximately transform each FGPM layer into a multilayered piezoelectric one with an equal and small thickness for each layer in comparison with the mid-surface radius, and with the homogeneous material properties determined in an average thickness sense; and a transfer matrix method is developed, so that the general solutions of the system equations can be obtained layer-by-layer, which is significantly less time-consuming than the usual approach. A parametric study is undertaken of the influence of the aspect ratio, open- and closed-circuit surface conditions, and material-property gradient index on the assorted field variables induced in the FGPM sandwich cylinders.     

弹性边界约束旋转功能梯度圆柱壳结构自由振动行波特性分析

对旋转功能梯度圆柱壳自由振动行波特性及边界约束影响进行了分析研究.将功能梯度材料的物理特性表示成沿壳体厚度方向指数变化的函数,基于Love壳体理论,将圆柱壳3个方向的振动位移场采用改进Fourier(傅立叶)级数方法展开, 进而改善位移函数在边界位置求导连续性,结合旋转圆柱壳结构能量原理描述与Rayleigh Ritz法,推导旋转功能梯度圆柱壳自由振动特征方程.通过将计算结果与现有文献结果对比验证了该文模型的正确性与收敛性.随后,通过算例讨论分析了功能梯度材料特性参数、几何参数、边界条件及约束弹簧刚度对旋转功能梯度圆柱壳自由振动行波振动特性的影响.结果表明:边界条件在环向波数n较小或长径比L/R较小的情况下对行波特性影响较为明显;随着厚径比H/R的增大,边界条件的影响逐渐减小;边界约束弹簧对行波特性影响程度取决于模态阶数情况;功能梯度材料特性参数对前后行波频率的影响随着模态序数的增大而逐渐增大.     

Free edge stress prediction in thick laminated cylindrical shell panel subjected to bending moment

A thick composite cylindrical shell panel with general layer stacking is studied to investigate the free edge and 3D stresses in the panel which is subjected to pure bending moment. To this aim, a Galerkin based layerwise formulation is presented to discretize the governing equation of the panel to ordinary differential equations. Employing a reduced displacement field for the cylindrical panel, the governing equations for thick panel are developed in terms of displacements and a set of coupled ordinary differential equations is obtained. The governing equations are solved analytically for free edge boundary conditions and applied pure bending moment. The accuracy of numerical results is examined and the distribution of interlaminar and in-plane stresses is studied. The free edge stresses are studied and the effect of radius to thickness ratio, width to thickness ratio and layer stacking on the distribution of stresses is investigated. The focus of numerical results is on the prediction of boundary layer and free edge stress distribution.     

The dynamic analysis of the functionally graded piezoelectric (FGP) shell panel based on three-dimensional elasticity theory

In this paper, three-dimensional elasticity solution is extended to investigate a FGPM finite length, simply supported shell panel under dynamic pressure excitation. The host panel is assumed to be of some functionally graded piezoelectric material (FGPM). The ordinary differential equations (o.d.e.) are derived from the highly coupled partial differential equations (p.d.e.) using series expansions of mechanical and electrical displacements. The resulting system of ordinary differential equations is solved by means of Galerkin finite element method. At last, numerical examples are presented for a FGPM shell panel. To verify the validity of code and formulation, the results of a FGM panel and a FGM plate are compared with the published results.     

Semi-analytic modelling of transversely isotropic magneto-electro-elastic materials under frictional sliding contact

Functionally graded magneto-electro-elastic (FGMEE) materials has been increasingly used in engineering applications, particularly in smart material or intelligent structure systems. This paper proposes a semi-analytical approach for sliding frictional contact problem between a rigid insulating sphere and a transversely isotropic FGMEE film and half-space based on frequency response functions (FRFs). Multilayered approximation is used to model the functionally graded material (FGM), and the FRFs for each MEE layer are derived explicitly. The unknown coefficients in FRFs are formulated by two matrix equations, and their efficient solution process is proposed. Based on the obtained FRFs, a highly efficient semi-analytical model (SAM) is developed which is able to solve the three-dimensional frictional contact of FGMEE materials with arbitrary layer designs. The model is validated with finite element method and the literature. Furthermore, the pressure/stress distribution and electric/magnetic potential are studied in different FGM designs to investigate the influence of material layout.     

Natural frequency analysis of functionally graded material beams with axially varying stochastic properties

The functionally graded material (FGM) has a potential to replace ordinary ones in engineering reality due to its superior thermal and dynamical characteristics. In this regard, the paper presents an effective approach for uncertain natural frequency analysis of composite beams with axially varying material properties. Rather than simply assuming the material model as a deterministic function, we further extend the FGM property as a random field, which is able to account for spatial variability in laboratory observations and in-field data. Due to the axially varying input uncertainty, natural frequencies of the stochastically FGM (S-FGM) beam become random variables. To this end, the Karhunen–Loève expansion is first introduced to represent the composite material random field as the summation of a finite number of random variables. Then, a generalized eigenvalue function is derived for stochastic natural frequency analysis of the composite beam. Once the mechanistic model is available, the brutal Monte-Carlo simulation (MCS) similar to the design of experiment can be used to estimate statistical characteristics of the uncertain natural frequency response. To alleviate the computational cost of the MCS method, a generalized polynomial chaos expansion model developed based on a rather small number of training samples is used to mimic the true natural frequency function. Case studies have demonstrated the effectiveness of the proposed approach for uncertain natural frequency analysis of functionally graded material beams with axially varying stochastic properties.     

Full 3D finite element Cosserat formulation with application in layered structures

In this paper, a full three-dimensional (3D) finite element Cosserat formulation is developed within the principles of continuum mechanics in the small deformation framework. The developed finite element formulation is general; however, the proposed constitutive laws incorporate the effect of the internal length parameter of 3D layered continua. The extension of the existing two-dimensional (2D) Cosserat formulation to the 3D framework is novel and is consistent with plate theory which can be considered as the 3D version of beam theory. The results demonstrate a high level of consistency with the analytical solutions predicted by plate theory as well as predictions by alternative numerical techniques such as the discrete element method.     

Transient piezothermoelastic analysis of a cross-ply laminated cylindrical panel bonded to a piezoelectric actuator

In this study, the theoretical treatment of transient piezothermoelastic problem is developed for a cross-ply laminated cylindrical panel bonded to a piezoelectric actuator due to nonuniform heat supply. By using the exact solutions for cross-ply laminate and piezoelectric layer of crystal class mm2, the theoretical analysis of a transient piezothermoelasticity is developed for a simple supported cylindrical composite panel under the state of plane strain. Analysis of a piezothermoelastic problem leads to an appropriate electric potential applied to the piezoelectric layer which suppresses the induced thermoelastic displacement in the radial direction at the midpoint on the free surface of the cross-ply laminate. Some numerical results for the temperature change, the displacement, the stress in a transient state when the transient thermoelastic displacement is controlled are shown in figures.     

Nonlinear free vibration analysis of SSMFG cylindrical shells resting on nonlinear viscoelastic foundation in thermal environment

In this paper, a semi-analytical method for the free vibration behavior of spiral stiffened multilayer functionally graded (SSMFG) cylindrical shells under the thermal environment is investigated. The distribution of linear and uniform temperature along the direction of thickness is assumed. The structure is embedded within a generalized nonlinear viscoelastic foundation, which is composed of a two-parameter Winkler-Pasternak foundation augmented by a Kelvin-Voigt viscoelastic model with a nonlinear cubic stiffness. The cylindrical shell has three layers consist of ceramic, FGM, and metal in two cases. In the first model i.e. Ceramic-FGM-Metal (CFM), the exterior layer of the cylindrical shell is rich ceramic while the interior layer is rich metal and the functionally graded material is located between these layers and the material distribution is in reverse order in the second model i.e. Metal-FGM-Ceramic (MFC). The material constitutive of the stiffeners is continuously changed through the thickness. Using the Galerkin method based on the von Kármán equations and the smeared stiffeners technique, the problem of nonlinear vibration has been solved. In order to find the nonlinear vibration responses, the fourth order Runge–Kutta method is utilized. The results show that the different angles of stiffeners and nonlinear elastic foundation parameters have a strong effect on the vibration behaviors of the SSMFG cylindrical shells. Also, the results illustrate that the vibration amplitude and the natural frequency for CFM and MFC shells with the first longitudinal and third transversal modes (m = 1, n = 3) with the stiffeners angle θ = 30°, β = 60° and θ = β = 30° is less than and more than others, respectively.     

An analytical method for magnetothermoelastic analysis of functionally graded hollow cylinders

This paper considers magnetothermoelastic behavior of a functionally graded material (FGM) hollow cylinder, placed in a uniform magnetic field, subjected to thermal and mechanical loads. Exact solutions for stresses and perturbations of the magnetic field vector in FGM hollow cylinders is determined by using the infinitesimal theory of magnetothermoelasticity. Numerical results indicate that the inhomogeneous constants presented in the present study are useful parameters from a design point of view in that it can be tailored for specific applications to control the stress and perturbation of magnetic field vector distributions. This research is helpful for the optimum design annular cylindrical FGM sensors/actuators.     

粘弹性方程全离散化有限体积元格式及数值模拟

本文利用有限体积元方法研究二维粘弹性方程, 给出一种时间二阶精度的全离散化有限体积元格式, 并给出这种全离散化有限体积元解的误差估计, 最后用数值例子验证数值结果与理论结果是相吻合的. 通过与有限元方法和有限差分方法相比较, 进一步说明了全离散化有限体积元格式是求解二维粘弹性方程数值解的最有效方法之一.     

Asymmetric Green's functions for a functionally graded transversely isotropic tri-material

This paper considers the elastic analysis of a functionally graded transversely isotropic tri-material solid under the arbitrary distribution of applied static loads. Using two displacement potential functions, for three-dimensional point-load and patch-load configurations, Green's functions for displacement and stress components are generated in the form of infinite line-integrals. These solutions are shown to be analytically reducible to the special cases of exponentially graded bi-material, exponentially graded half-space and homogeneous tri-material Green's functions. It also encompasses a functionally graded finite layer on a rigid base with surface loading with two cases of interfacial conditions, rigid-bonded and rigid-frictionless. Finally, for the special case of a functionally graded layer sandwiched between two homogeneous layers, using several numerical displays, the effect of material inhomogeneity on the responses is studied and the accuracy of numerical scheme is verified.     

Time-harmonic response of transversely isotropic and layered poroelastic half-spaces under general buried loads

A semi-analytical method is developed for solving the dynamic response of transversely isotropic, multilayered, and poroelastic half-spaces with different surface hydraulic conditions and subjected to time-harmonic vertical and horizontal loads buried in the layered half-space. The coupled governing equations of motion are presented in details in terms of the Biot's poroelastodynamic theory via the (u,p) formulation. The cylindrical system of vector functions is introduced to express the unknown primary quantities so that the coupled governing partial differential equations can be reduced and separated into two sets of first-order ordinary differential equations (i.e., the LM- and N-types). A recursive relation for the expansion coefficients among different layers is established by virtue of the stable and efficient dual variable and position method. Making use of the boundary and interface conditions, the fundamental solutions are obtained in terms of the vector-function system. The corresponding physical-domain solutions are then derived via an accurate semi-infinite integral algorithm. The developed fundamental solutions are carefully checked with existing solutions, and numerical examples are further presented to demonstrate the effect of material anisotropy, loading depth, material layering, and surface hydraulic condition on the dynamic response, which should be useful to design engineers. These solutions could be further served as benchmarks for future numerical methods.