Stress analysis of functionally graded shells using a 7-parameter shell element

In this paper, finite element stress analysis of functionally graded structures using a high-order spectral/hp shell finite element is presented. The shell element is based on a seven-parameter first-order shear deformation theory in which the seventh parameter, in addition to the usual six degrees of freedom, is the thickness stretch. The continuum shell element is utilized for the numerical simulations of the fully geometrically nonlinear response of functionally graded elastic shell structures. Several nontrivial shell problems are considered to report deflections and stresses, the latter being the main focus of the current paper. It is found that the stresses computed in the current study agree only in some cases with those of ANSYS and/or ABAQUS and thus requires additional study to determine the cause of the disagreement.     

Coupled thermoelasticity of functionally graded cylindrical shells

The coupled thermoelastic response of a functionally graded circular cylindrical shell is studied. The coupled thermoelastic and the energy equations are simultaneously solved for a functionally graded axisymmetric cylindrical shell subjected to thermal shock load. A second-order shear deformation shell theory that accounts for the transverse shear strains and rotations is considered. Including the thermo-mechanical coupling and rotary inertia, a Galerkin finite element formulation in space domain and the Laplace transform in time domain are used to formulate the problem. The inverse Laplace transform is obtained using a numerical algorithm. The shell is graded through the thickness assuming a volume fraction of metal and ceramic, using a power law distribution. The results are validated with the known data in the literature.     

Finite elements based on a first-order shear deformation moderate rotation shell theory with applications to the analysis of composite structures

The first-order shear deformation moderate rotation shell theory of Schmidt and Reddy [R. Schmidt and J. N. Reddy, J. Appl. Mech. 55, 611–617 (1988)] is used as a basis for the development of finite element models for the analysis of the static, geometrically non-linear response of anisotropic and laminated structures. The incremental, total Lagrangian formulation of the theory is developed, and numerical solutions are obtained by using the isoparametric Lagrangian 9-node and Serendipity 8-node shell finite elements. Various integration schemes (full, selective reduced, and uniformly reduced integration) are applied in order to detect and to overcome the effects of shear and membrane locking on the predicted structural response. A number of sample problems of isotropic, orthotropic, and multi-layered structures are presented to show the accuracy of the present theory. The von Kármán-type first-order shear deformation shell theory and continuum 2D theory are used for comparative analyses.     

Postbuckling of FGM cylindrical shells under combined axial and radial mechanical loads in thermal environments

A postbuckling analysis is presented for a shear deformable functionally graded cylindrical shell of finite length subjected to combined axial and radial loads in thermal environments. Heat conduction and temperature-dependent material properties are both taken into account. The temperature field considered is assumed to be a uniform distribution over the shell surface and varied in the thickness direction only. Material properties are assumed to be temperature-dependent, and graded in the thickness direction according to a simple power law distribution in terms of the volume fractions of the constituents. The formulations are based on a higher order shear deformation shell theory with von Kármán–Donnell-type of kinematic nonlinearity. A boundary layer theory of shell buckling, which includes the effects of nonlinear prebuckling deformations, large deflections in the postbuckling range, and initial geometric imperfections of the shell, is extended to the case of functionally graded cylindrical shells. A singular perturbation technique is employed to determine the interactive buckling loads and postbuckling equilibrium paths. The numerical illustrations concern the postbuckling response of perfect and imperfect cylindrical shells with two constituent materials subjected to combined axial and radial mechanical loads and under different sets of thermal environments. The results reveal that the temperature field and volume fraction distribution have a significant effect on the postbuckling behavior, but they have a small effect on the imperfection sensitivity of the functionally graded shell.     

Postbuckling analysis of axially loaded functionally graded cylindrical panels in thermal environments

A postbuckling analysis is presented for a functionally graded cylindrical panel of finite length subjected to axial compression in thermal environments. Material properties are assumed to be temperature dependent, and graded in the thickness direction according to a simple power law distribution in terms of the volume fractions of the constituents. The governing equations of a functionally graded cylindrical panel are based on Reddy’s higher order shear deformation shell theory with a von Kármán–Donnell-type of kinematic nonlinearity and including thermal effects. Two cases of the in-plane boundary conditions are considered. The nonlinear prebuckling deformations and initial geometric imperfections of the panel are both taken into account. A boundary layer theory of shell buckling, which includes the effects of nonlinear prebuckling deformations, large deflections in the postbuckling range, and initial geometric imperfections of the shell, is extended to the case of functionally graded cylindrical panels under axial compression. A singular perturbation technique is employed to determine the buckling loads and postbuckling equilibrium paths. The numerical illustrations concern the postbuckling behavior of axially loaded, perfect and imperfect, functional graded cylindrical panels with two constituent materials and under different sets of thermal environments. The influences played by temperature rise, volume fraction distributions, the character of in-plane boundary conditions, transverse shear deformation, panel geometric parameters, as well as initial geometric imperfections are studied.     

Nonlinear vibration of functionally graded graphene platelet-reinforced composite truncated conical shell using first-order shear deformation theory

In this study, the first-order shear deformation theory(FSDT) is used to establish a nonlinear dynamic model for a conical shell truncated by a functionally graded graphene platelet-reinforced composite(FG-GPLRC). The vibration analyses of the FG-GPLRC truncated conical shell are presented. Considering the graphene platelets(GPLs) of the FG-GPLRC truncated conical shell with three different distribution patterns, the modified Halpin-Tsai model is used to calculate the effective Young's modulus. Hamilton's principle, the FSDT, and the von-Karman type nonlinear geometric relationships are used to derive a system of partial differential governing equations of the FG-GPLRC truncated conical shell. The Galerkin method is used to obtain the ordinary differential equations of the truncated conical shell. Then, the analytical nonlinear frequencies of the FG-GPLRC truncated conical shell are solved by the harmonic balance method. The effects of the weight fraction and distribution pattern of the GPLs, the ratio of the length to the radius as well as the ratio of the radius to the thickness of the FG-GPLRC truncated conical shell on the nonlinear natural frequency characteristics are discussed. This study culminates in the discovery of the periodic motion and chaotic motion of the FG-GPLRC truncated conical shell.     

功能梯度板动力分析的光滑有限元法研究

为提高功能梯度板动响应问题求解精度,基于一阶剪切变形板理论,提出了求解功能梯度板自由振动问题的Cell-Based 光滑有限元格式．功能梯度板材料属性沿厚度方向呈梯度连续变化,计算系统刚度矩阵时在光滑域内进行光滑梯度操作,可提高求解精度．采用Cell-Based 光滑有限元法,讨论了长厚比、形状因子和边界条件对两种典型功能梯度方板无量纲频率参数的影响,并与FEM法的解和文献中的解做了对比．数值算例的结果表明,光滑梯度操作可改善有限元系统的刚度,Cell-Based 光滑有限元法的计算结果更加逼近真实解,从而为功能梯度材料的进一步应用打下基础．     

Free vibration analysis of functionally graded laminated sandwich cylindrical shells integrated with piezoelectric layer

An analytical method for the three-dimensional vibration analysis of a functionally graded cylindrical shell integrated by two thin functionally graded piezoelectric (FGP) layers is presented. The first-order shear deformation theory is used to model the electromechanical system. Nonlinear equations of motion are derived by considering the von Karman nonlinear strain-displacement relations using Hamilton’s principle. The piezoelectric layers on the inner and outer surfaces of the core can be considered as a sensor and an actuator for controlling characteristic vibration of the system. The equations of motion are derived as partial differential equations and then discretized by the Navier method. Numerical simulation is performed to investigate the effect of different parameters of material and geometry on characteristic vibration of the cylinder. The results of this study show that the natural frequency of the system decreases by increasing the non-homogeneous index of FGP layers and decreases by increasing the non-homogeneous index of the functionally graded core. Furthermore, it is concluded that by increasing the ratio of core thickness to cylinder length, the natural frequencies of the cylinder increase considerably.     

An efficient and simple refined theory for nonlinear bending analysis of functionally graded sandwich plates

In this paper, an efficient and simple refined theory is presented for nonlinear bending analysis of functionally graded sandwich plates. The theory presented is variationally consistent, does not require the shear correction factor, and gives rise to transverse shear stress variations such that the transverse shear stresses vary parabolically across the plate thickness, satisfying shear-stress-free surface conditions. The energy concept along with the present theory and the first- and third-order shear deformation theories is used to predict the large deflection and the stress distribution across the thickness of functionally graded sandwich plates.     

R-functions theory applied to investigation of nonlinear free vibrations of functionally graded shallow shells

Nonlinear free vibration of functionally graded shallow shells with complex planform is investigated using the R-functions method and variational Ritz method. The proposed method is developed in the framework of the first-order shear deformation shallow shell theory. Effect of transverse shear strains and rotary inertia is taken into account. The properties of functionally graded materials are assumed to be varying continuously through the thickness according to a power law distribution. The Rayleigh–Ritz procedure is applied to obtain the frequency equation. Admissible functions are constructed by the R-functions theory. To implement the proposed approach, the corresponding software has been developed. Comprehensive numerical results for three types of shallow shells with positive, zero and negative curvature with complex planform are presented in tabular and graphical forms. The convergence of the natural frequencies with increasing number of admissible functions has been checked out. Effect of volume fraction exponent, geometry of a shape and boundary conditions on the natural and nonlinear frequencies is brought out. For simply supported rectangular FG shallow shells, the results obtained are compared with those available in the literature. Comparison demonstrates a good accuracy of the approach proposed.     

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

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

基于物理中面FGM中厚板自由振动的有限元分析

基于物理中面和一阶剪切变形板理论,研究了不同边界条件下功能梯度材料(FGM)中厚板的自由振动问题．假设功能梯度板的材料性质沿厚度方向按幂函数规律连续变化．根据哈密顿原理建立了FGM板有限元形式的自由振动方程,利用MATLAB软件编写程序进行了计算．通过数值算例,讨论了不同边界条件下FGM中厚板的无量纲频率随材料梯度指数和厚宽比的变化情况,并与经典板理论下的频率进行了比较．     

Thermoelastic buckling response of thick functionally graded plates

The thermoelastic buckling behavior of a thick plate made of a functionally graded material is investigated in this paper by using an exponential shear deformation plate theory. A simple power law based on the rule of mixtures is used to estimate the effective material properties as functions of the plate thickness. The neutral surface position for such functionally graded plates is determined on the basis of the nonlinear strain-displacement relations. Uniform, linear, and nonlinear temperature distributions across the plate are considered. An analytical approach is presented to find the critical buckling temperature, which can be used in engineering calculations. A numerical solution of the problem with the use of an exponential dependence for shear strains is presented. The results obtained are compared with available data.     

Thermomechanical vibration analysis of a functionally graded shell with flowing fluid

This paper reports the results of an investigation into the vibration of functionally graded cylindrical shells with flowing fluid, embedded in an elastic medium, under mechanical and thermal loads. By considering rotary inertia, the first-order shear deformation theory (FSDT) and the fluid velocity potential, the dynamic equation of functionally graded cylindrical shells with flowing fluid is derived. Here, heat conduction equation along the thickness of the shell is applied to determine the temperature distribution and material properties are assumed to be graded distribution along the thickness direction according to a power-law in terms of the volume fractions of the constituents. The equations of eigenvalue problem are obtained by using a modal expansion method. In numerical examples, effects of material composition, thermal loading, static axial loading, flow velocity, medium stiffness and shell geometry parameters on the free vibration characteristics are described. The new features in this paper are helpful for the application and the design of functionally graded cylindrical shells containing fluid flow.     

Thermal instability of functionally graded deep spherical shell

Thermal instability of deep spherical shells made of functionally graded material (FGM) is studied in this paper. The governing equations are based on the first-order theory of shells and the Sanders nonlinear kinematics equations. It is assumed that the mechanical properties are linear functions of thickness coordinate. The constituent material of the functionally graded shell is assumed to be a mixture of ceramic and metal. The analytical solutions are obtained for three types of thermal loadings including the uniform temperature rise (UTR), the linear radial temperature (LRT), and the nonlinear radial temperature (NRT). Results are validated with the known data in literature.     

基于高阶剪切变形理论的四边形求积元板单元及其应用

近年来由各类新型复合材料或功能梯度材料构成的板结构在工程领域得到了广泛应用,其显著特点是材料性能沿板厚变化.为合理考虑横向剪切应变,许多学者基于Reddy高阶剪切变形理论,构建了不同的有限元单元对该类板结构进行分析,但其中满足$C^{1}$连续条件的单元相对较少.本文基于Reddy高阶剪切变形理论,采用求积元方法,建立了$C^{1}$连续的四边形板单元.利用该单元对均质材料、复合材料、功能梯度材料构成的等厚度矩形板、变厚度矩形板及等厚度斜板的线弹性弯曲和自由振动问题进行了计算分析,并与现有文献中的相应计算结果进行了对比.研究表明:基于高阶剪切变形理论的四边形求积元板单元具有较高的计算效率和良好的适应性,文中各类材料构成的等变厚度矩形板及等厚度斜板均只需1个单元即可得到理想的计算结果.对于等/变厚度矩形板,可仅使用9$\times$9个积分点,而对于等厚度斜板,随着斜角的增大,所需积分点的数目逐渐增多至15$\times $15.该四边形求积元板单元可进一步用于新型复合材料板的非线性分析.      

Nonlinear Transient Response of Functionally Graded Material Sandwich Doubly Curved Shallow Shell Using New Displacement Field

In this paper, the nonlinear transient dynamic response of functionally graded material(FGM) sandwich doubly curved shell with homogenous isotropic material core and functionally graded face sheet is analyzed using a new displacement field on the basis of Reddy's third-order shear theory for the first time. The equivalent material properties for the FGM face sheet are assumed to obey the rule of simple power law function in the thickness direction.Based on Reddy'stheory of higher shear deformation, a new displacement field is developed by introducing the secant function into transverse displacement. Four coupled nonlinear differential equations are obtained by applying Hamilton's principle and Galerkin method. It is assumed that the FGM sandwich doubly curved shell is subjected to step loading, air-blast loading,triangular loading, and sinusoidal loading, respectively. On the basis of double-precision variablecoefficient ordinary differential equation solver, a new program code in FORTRAN software is developed to solve the nonlinear transient dynamics of the system. The influences of core thickness, volume fraction, core-to-face sheet thickness ratio, width-to-thickness ratio and blast type on the transient response of the shell are discussed in detail through numerical simulation.     

纵横载荷作用下功能梯度梁的弯曲和过屈曲

基于一阶非线性梁理论和物理中面概念,导出了纵横向载荷作用下功能梯度材料(FGM)梁非线性弯曲和过屈曲问题的控制方程,并获得了该问题的精确解;据此解研究了梯度材料性质、外载荷、横向剪切变形以及边界条件等因素对功能梯度材料梁非线性力学行为的影响,分析中假设功能梯度材料性质只沿梁厚度方向,并按成分含量的幂指数函数形式变化。结果表明:纵横载荷共同作用下,功能梯度梁的弯曲构形将有无限多个;随着梯度指数的增大,梁的变形减小,临界载荷升高;随着长高比的增大,横向剪切变形的影响减小。     

Rotating sandwich cylindrical shells with an FGM core and two FGPM layers: free vibration analysis

The free vibration analysis of a rotating cylindrical shell with an analytical method is investigated. The shell is considered as a sandwich structure, where the middle layer is a functionally graded material(FGM) shell, and it is surrounded by two piezoelectric layers. Considering piezoelectric materials to be functionally graded(FG),the material properties vary along the thickness direction as one innovation of this study.Applying the first-order shear deformation theory(FSDT), the equations of motion of this electromechanical system are derived as the partial differential equations(PDEs) using Hamilton's principle. Then, the Galerkin procedure is used to discretize the governing equations, and the present results are compared with the previously published results for both isotropic and FGM shells to verify the analytical method. Finally, the effects of FGM and functionally graded piezoelectric material(FGPM) properties as well as the thickness ratio and the axial and circumferential wave numbers on the natural frequencies are studied. Moreover, the Campbell diagram is plotted and discussed through the governing equations. The present results show that increasing the non-homogeneous index of the FGM decreases the natural frequencies on the contrary of the effect of non-homogeneous index of the FGPM.