Free vibration analysis of functionally graded material beams based on Levinson beam theory

Free vibration response of functionally graded material (FGM) beams is studied based on the Levinson beam theory (LBT). Equations of motion of an FGM beam are derived by directly integrating the stress-form equations of elasticity along the beam depth with the inertial resultant forces related to the included coupling and higherorder shear strain. Assuming harmonic response, governing equations of the free vibration of the FGM beam are reduced to a standard system of second-order ordinary differential equations associated with boundary conditions in terms of shape functions related to axial and transverse displacements and the rotational angle. By a shooting method to solve the two-point boundary value problem of the three coupled ordinary differential equations, free vibration response of thick FGM beams is obtained numerically. Particularly, for a beam with simply supported edges, the natural frequency of an FGM Levinson beam is analytically derived in terms of the natural frequency of a corresponding homogenous Euler-Bernoulli beam. As the material properties are assumed to vary through the depth according to the power-law functions, the numerical results of frequencies are presented to examine the effects of the material gradient parameter, the length-to-depth ratio, and the boundary conditions on the vibration response.     

A higher-order theory for static and dynamic analyses of functionally graded beams

The higher-order theory is extended to functionally graded beams (FGBs) with continuously varying material properties. For FGBs with shear deformation taken into account, a single governing equation for an auxiliary function F is derived from the basic equations of elasticity. It can be used to deal with forced and free vibrations as well as static behaviors of FGBs. A general solution is constructed, and all physical quantities including transverse deflection, longitudinal warping, bending moment, shear force, and internal stresses can be represented in terms of the derivatives of F. The static solution can be determined for different end conditions. Explicit expressions for cantilever, simply supported, and clamped-clamped FGBs for typical loading cases are given. A comparison of the present static solution with existing elasticity solutions indicates that the method is simple and efficient. Moreover, the gradient variation of Young’s modulus and Poisson’s ratio may be arbitrary functions of the thickness direction. Functionally graded Rayleigh and Euler–Bernoulli beams are two special cases when the shear modulus is sufficiently high. Moreover, the classical Levinson beam theory is recovered from the present theory when the material constants are unchanged. Numerical computations are performed for a functionally graded cantilever beam with a gradient index obeying power law and the results are displayed graphically to show the effects of the gradient index on the deflection and stress distribution, indicating that both stresses and deflection are sensitive to the gradient variation of material properties.     

Free vibration of functionally graded beams based on both classical and first-order shear deformation beam theories

The free vibration of functionally graded material （FGM） beams is studied based on both the classical and the first-order shear deformation beam theories. The equations of motion for the FGM beams are derived by considering the shear deforma- tion and the axial, transversal, rotational, and axial-rotational coupling inertia forces on the assumption that the material properties vary arbitrarily in the thickness direction. By using the numerical shooting method to solve the eigenvalue problem of the coupled ordinary differential equations with different boundary conditions, the natural frequen- cies of the FGM Timoshenko beams are obtained numerically. In a special case of the classical beam theory, a proportional transformation between the natural frequencies of the FGM and the reference homogenous beams is obtained by using the mathematical similarity between the mathematical formulations. This formula provides a simple and useful approach to evaluate the natural frequencies of the FGM beams without dealing with the tension-bending coupling problem. Approximately, this analogous transition can also be extended to predict the frequencies of the FGM Timoshenko beams. The numerical results obtained by the shooting method and those obtained by the analogous transformation are presented to show the effects of the material gradient, the slenderness ratio, and the boundary conditions on the natural frequencies in detail.     

CLASSICAL AND HOMOGENIZED EXPRESSIONS FOR BUCKLING SOLUTIONS OF FUNCTIONALLY GRADED MATERIAL LEVINSON BEAMS

The relationship between the critical buckling loads of functionally graded material(FGM) Levinson beams(LBs) and those of the corresponding homogeneous Euler-Bernoulli beams(HEBBs) is investigated. Properties of the beam are assumed to vary continuously in the depth direction. The governing equations of the FGM beam are derived based on the Levinson beam theory, in which a quadratic variation of the transverse shear strain through the depth is included.By eliminating the axial displacement as well as the rotational angle in the governing equations,an ordinary differential equation in terms of the deflection of the FGM LBs is derived, the form of which is the same as that of HEBBs except for the definition of the load parameter. By solving the eigenvalue problem of ordinary differential equations under different boundary conditions clamped(C), simply-supported(S), roller(R) and free(F) edges combined, a uniform analytical formulation of buckling loads of FGM LBs with S-S, C-C, C-F, C-R and S-R edges is presented for those of HEBBs with the same boundary conditions. For the C-S beam the above-mentioned equation does not hold. Instead, a transcendental equation is derived to find the critical buckling load for the FGM LB which is similar to that for HEBB with the same ends. The significance of this work lies in that the solution of the critical buckling load of a FGM LB can be reduced to that of the HEBB and calculation of three constants whose values only depend upon the throughthe-depth gradient of the material properties and the geometry of the beam. So, a homogeneous and classical expression for the buckling solution of FGM LBs is accomplished.     

基于一种广义梁理论的弹性地基FGM简支梁自由振动解析解

基于一种扩展的n阶广义剪切变形梁理论(n-GBT),应用Hamilton原理,建立了以轴向位移、横向位移及转角为未知函数的Winkler-Pasternak弹性地基功能梯度材料(FGM)梁的自由振动方程,采用Navier法获得了弹性地基FGM简支梁自由振动的精确解。与多种梁理论预测结果进行比较,讨论并给出了GBT阶次n的理想取值;分析了梯度指标、跨厚比及地基刚度对FGM梁频率的影响。结果表明:本文方法有效且适用范围广,若采用高阶剪切梁理论模型,宜取n≥3的奇数;FGM梁的自振频率随材料梯度指标的增大而减小;随跨厚比的增加而增大,但当跨厚比大于20,跨厚比增加对频率的影响很小;随地基刚度的增加而增大,地基刚度足够大时,频率趋于收敛。     

功能梯度梁与均匀梁静动态解间的相似转换

基于Euler-Bernoulli 梁理论, 研究了功能梯度材料梁的弯曲、屈曲和自由振动. 通过分析和比较功能梯度材料梁 和均匀梁的控制方程, 得到了功能梯度材料梁与均匀梁的解之间的相似转换关系. 在给定功 能梯度材料梁的材料性质在横向按幂函数分布的情况下, 导出了解之间的相似转换系数的解 析表达式. 该系数集中反映了功能梯度梁的材料非均匀性. 因此, 可将功能梯度材料梁的静 动态问题的求解转换为同样载荷和边界条件下均匀梁的静动态问题求解以及相似转换系数的 计算.     

功能梯度与均匀圆板弯曲解的线性转换关系

基于一阶剪切理论, 研究了功能梯度材料圆板与均匀圆板轴对称弯曲解之间的线性转换关系. 通过理论分析和比较 功能梯度材料圆板和均匀圆板在一阶剪切理论下的位移形式的轴对称弯曲控制方程, 发现了功能梯度材料圆板的转角与均匀圆板的转角之间的相似转换关系. 在假设材料性质沿板厚连续变化的情况下, 给出了相似转换系数的解析表达式. 在此基础上, 进一步导出了一阶剪切理论下功能梯度圆板的挠度与经典理论下, 均匀圆板的挠度之间的线性关系. 从而, 可将功能梯度材料圆板在一阶剪切理论下的弯曲问题求解, 转化为相应均匀薄圆板在经典理论下的弯曲问题求解, 以及转换系数的计算问题. 这一方法为功能梯度非均匀中厚度圆板的求解提供了简捷途径, 而且更便于工程应用. 作为例子, 采用上述方法分别求得了周边简支和夹紧条件下, 梯度圆板在均布横向载荷作用下的弯曲解析解, 该解答与Reddy得到的结果完全吻合.      

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

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

Bending and free vibrational analysis of bi-directional functionally graded beams with circular cross-section

The bending and free vibrational behaviors of functionally graded (FG) cylindrical beams with radially and axially varying material inhomogeneities are investigated. Based on a high-order cylindrical beam model, where the shear deformation and rotary inertia are both considered, the two coupled governing differential motion equations for the deflection and rotation are established. The analytical bending solutions for various boundary conditions are derived. In the vibrational analysis of FG cylindrical beams, the two governing equations are firstly changed to a single equation by means of an auxiliary function, and then the vibration mode is expanded into shifted Chebyshev polynomials. Numerical examples are given to investigate the effects of the material gradient indices on the deflections, the stress distributions, and the eigenfrequencies of the cylindrical beams, respectively. By comparing the obtained numerical results with those obtained by the three-dimensional (3D) elasticity theory and the Timoshenko beam theory, the effectiveness of the present approach is verified.     

A microstructure-dependent Timoshenko beam model based on a modified couple stress theory

A microstructure-dependent Timoshenko beam model is developed using a variational formulation. It is based on a modified couple stress theory and Hamilton's principle. The new model contains a material length scale parameter and can capture the size effect, unlike the classical Timoshenko beam theory. Moreover, both bending and axial deformations are considered, and the Poisson effect is incorporated in the current model, which differ from existing Timoshenko beam models. The newly developed non-classical beam model recovers the classical Timoshenko beam model when the material length scale parameter and Poisson's ratio are both set to be zero. In addition, the current Timoshenko beam model reduces to a microstructure-dependent Bernoulli-Euler beam model when the normality assumption is reinstated, which also incorporates the Poisson effect and can be further reduced to the classical Bernoulli-Euler beam model. To illustrate the new Timoshenko beam model, the static bending and free vibration problems of a simply supported beam are solved by directly applying the formulas derived. The numerical results for the static bending problem reveal that both the deflection and rotation of the simply supported beam predicted by the new model are smaller than those predicted by the classical Timoshenko beam model. Also, the differences in both the deflection and rotation predicted by the two models are very large when the beam thickness is small, but they are diminishing with the increase of the beam thickness. Similar trends are observed for the free vibration problem, where it is shown that the natural frequency predicted by the new model is higher than that by the classical model, with the difference between them being significantly large only for very thin beams. These predicted trends of the size effect in beam bending at the micron scale agree with those observed experimentally. Finally, the Poisson effect on the beam deflection, rotation and natural frequency is found to be significant, which is especially true when the classical Timoshenko beam model is used. This indicates that the assumption of Poisson's effect being negligible, which is commonly used in existing beam theories, is inadequate and should be individually verified or simply abandoned in order to obtain more accurate and reliable results.     

Free vibration characteristics of sectioned unidirectional/bidirectional functionally graded material cantilever beams based on finite element analysis

Advancements in manufacturing technology, including the rapid development of additive manufacturing (AM), allow the fabrication of complex functionally graded material (FGM) sectioned beams. Portions of these beams may be made from different materials with possibly different gradients of material properties. The present work proposes models to investigate the free vibration of FGM sectioned beams based on onedimensional (1D) finite element analysis. For this purpose, a sample beam is divided into discrete elements, and the total energy stored in each element during vibration is computed by considering either the Timoshenko or Euler-Bernoulli beam theory. Then, Hamilton's principle is used to derive the equations of motion for the beam. The effects of material properties and dimensions of FGM sections on the beam's natural frequencies and their corresponding mode shapes are then investigated based on a dynamic Timoshenko model (TM). The presented model is validated by comparison with three-dimensional (3D) finite element simulations of the first three mode shapes of the beam.     

功能梯度与均匀圆板静动态解之间的相似转换关系

基于经典板理论,研究了功能梯度材料圆板的轴对称弯曲、屈曲和自由振动解与相应的均匀材料圆板解之间的转换关系.通过消去拉-弯耦合项得到了以挠度函数表示的功能梯度圆板的弯曲、屈曲和自由振动控制方程.分析功能梯度圆板与均匀圆板的控制方程之间的相似性,得到了功能梯度材料圆板与均匀圆板的解之间解的相似转换关系,在假定FGM圆板的材料性质沿厚分别以幂函数和指数函数的度变规律后,给出了相应的转换系数的解析表达式.该系数集中反映了功能梯度圆板的材料非均匀性.在已知均匀材料圆板轴对称解的条件下,可将功能梯度材料圆板轴对称问题的求解转化为相似转换系数的计算问题.这一方法可为非均匀板的求解提供了十分便捷有效的途径,而且便于工程应用.     

Bending of functionally graded cantilever beam with power-law non-linearity subjected to an end force

A realistic beam structure often exhibits material and geometrical non-linearity, in particular for those made of metals. The mechanical behaviors of a non-linear functionally graded-material (FGM) cantilever beam subjected to an end force are investigated by using large and small deformation theories. Young's modulus is assumed to be depth-dependent. For an FGM beam of power-law hardening, the location of the neutral axis is determined. The effects of depth-dependent Young's modulus and non-linearity parameter on the deflections and rotations of the FGM beams are analyzed. Our results show that different gradient indexes may change the bending stiffness of the beam so that an FGM beam may bear larger applied load than a homogeneous beam when choosing appropriate gradients. Moreover, the bending stress distribution in an FGM beam is completely different from that in a homogeneous beam. The bending stress arrives at the maximum tensile stress at an internal position rather than at the surface. Obtained results are useful in safety design of linear and non-linear beams.     

饱和多孔弹性Timoshenko梁的大挠度分析

基于微观不可压饱和多孔介质理论和弹性梁的大挠度变形假设,考虑梁剪切变形效应,在梁轴线不可伸长和孔隙流体仅沿轴向扩散的限定下,建立了饱和多孔弹性Timoshenko梁大挠度弯曲变形的非线性数学模型.在此基础上,利用Galerkin截断法,研究了两端可渗透简支饱和多孔Timoshenko梁在突加均布横向载荷作用下的拟静态弯曲,给出了饱和多孔 Timoshenko梁弯曲变形时固相挠度、弯矩和孔隙流体压力等效力偶等随时间的响应.比较了饱和多孔Timoshenko梁非线性大挠度和线性小挠度理论以及饱和多孔 Euler-Bernoulli梁非线性大挠度理论的结果,揭示了他们间的差异,指出当无量纲载荷参数q＞l0时,应采用饱和多孔Timoshenko梁或Euler-Bernoulli梁的大挠度数学模型进行分析,特别的,当梁长细比λ＜30时,应采用饱和多孔Timoshenko梁大挠度数学模型进行分析.     

功能梯度中厚圆/环板轴对称弯曲问题的解析解

基于一阶剪切变形板理论，导出了热/机载荷作用下，位移形式的功能梯度 中厚圆/环板轴对称弯曲问题的控制方程，获得了问题的位移和内力的一般解析解. 作为特 例，分别研究了边缘径向固定和可动的夹紧和简支的4种实心功能梯度圆板，给出了它们的 解，并分析了热/机载荷作用下解的形态，讨论了横向剪切变形、材料梯度常数和边界条件， 对板的轴对称弯曲行为的影响.     

考虑剪切效应的旋转FGM楔形梁刚柔耦合动力学建模与仿真

本文对一类中心刚体-柔性梁系统在大范围转动下的刚柔耦合动力学问题进行了研究. 柔性梁为功能梯度材料(functionally graded materials, FGM)楔形变截面梁,材料体积分数在梁轴向呈幂律分布变化. 以弧长坐标来描述柔性FGM梁的几何位移关系,分别使用倾角和拉伸应变变量描述柔性梁的横向弯曲和纵向拉伸变形,并计及剪切效应. 采用假设模态法离散变形场,运用第二类拉格朗日方程进行方程推导,得到系统考虑剪切效应的刚柔耦合动力学模型. 基于全新的刚柔耦合动力学建模理论,研究不同轴向材料梯度分布的FGM楔形梁,通过数值仿真计算,分析讨论不同的转速、梯度分布规律以及变截面参数对系统动力学特性的影响. 结果表明,剪切效应对大高跨比的FGM楔形梁的变形影响较为明显,不容忽略;材料梯度分布规律和截面参数的选取均会对旋转FGM楔形梁的动力学响应和频率产生较大影响. 本文提出的考虑剪切效应的倾角刚柔耦合动力学模型是对以往非剪切模型的进一步完善,可应用于工程中的 Timoshenko梁结构的动力学问题求解.      

非保守功能梯度弹性组合曲梁的精确模型及其数值解

基于直法线假设,采用可伸长梁的几何非线性理论,建立了功能梯度材料弹性组合曲梁受切线均布随从力作用下的静态大变形数学模型。该模型不仅计及了轴线伸长,同时也精确地考虑了梁的初始曲率对变形的影响以及轴向变形与弯曲变形之间的耦合效应。用打靶法数值求解了由金属和陶瓷两相材料所构成的一种FGM组合曲梁在沿轴线均布切向随动载荷作用下的非线性平面弯曲问题,给出了不同梯度指标下FGM弹性曲梁随载荷参数大范围变化的平衡路径,并与金属和陶瓷两种单相材料曲梁的相应特性进行了比较。     

A size-dependent Reddy–Levinson beam model based on a strain gradient elasticity theory

A size-dependent Reddy–Levinson beam model is developed based on a strain gradient elasticity theory. Governing equations and boundary conditions are derived by using Hamilton’s principle. The model contains three material length scale parameters, which may effectively capture the size effect in micron or sub-micron. This model can degenerate into the modified couple stress model or even the classical model if two or all material length scale parameters are taken to be zero respectively. In addition, the present model recovers the micro scale Timoshenko and Bernoulli–Euler beam models based on the same strain gradient elasticity theory. To illustrate the new model, the static bending and free vibration problems of a simply supported micro scale Reddy–Levinson beam are solved respectively; the results are compared with the reduced models. Numerical results reveal that the differences in the deflection, rotation and natural frequency predicted by the present model and the other two reduced Reddy–Levinson models are getting larger as the beam thickness is comparable to the material length scale parameters. These differences, however, are decreasing or even diminishing with the increase of the beam thickness. This study may be helpful to characterize the mechanical properties of small scale beam-like structures for a wide range of potential applications.     

考虑前屈曲耦合变形时功能梯度简支梁的稳定性分析

在某些边界条件下,功能梯度材料（FGM）梁会由于拉弯耦合产生前屈曲耦合变形,而该变形对FGM梁的稳定性有影响。本文假设FGM梁的材料性质只沿厚度方向进行变化,基于经典非线性梁理论和物理中面概念,推导出FGM梁的平衡方程以及包含前屈曲耦合变形影响的屈曲控制方程,并用打靶法进行数值求解。讨论了前屈曲耦合变形、梯度指数以及材料性质的温度依赖等因素对FGM梁非线性变形和稳定性的影响。     

经典理论与一阶理论之间简支梁特征值的解析关系

利用Euler-Bernoulli梁理论(EBT)和Timoshenko梁理论(一阶理论,TBT)之间,梁的特征值问题在数学上的相似性,研究了不同梁理论之间特征值的关系。将特征值问题的求解转化为一个代数方程的求解,并导出了不同梁理论之间梁的特征值之间的精确解析关系。因此,只要已知梁的经典结果(临界载荷和固有频率),便很容易从这些关系中获得一阶梁理论下的相应结果。这些解析结果清楚地显示了横向剪切变形对经典结果影响的本质特点。另外,从这些关系中获得的含有剪切变形影响的结果,可以用于检验一阶理论下梁特征值数值结果的有效性、收敛性以及精确性等问题。