Non-uniform rational B-spline based free vibration analysis of axially functionally graded tapered Timoshenko curved beams

The free vibration of axially functionally graded (FG) tapered Timoshenko curved beams is studied with the numerical approach. By using the non-uniform rational B-spline (NURBS) basis functions, the exact geometry and the generalized displacement field are formulated. Variable geometric parameters and material properties, including the curvature, cross-sectional area, area moment of inertia, mass density, and Young’s modulus, are expanded as functions of the coordinate in a parametric domain. Based on Hamilton’s principle, the weak formulation is derived by applying a refined constitutive relation which considers the thickness effect. Natural frequencies and mode shapes are obtained from the eigenvalue equation. Circular, elliptic, and parabolic curved beams are considered in numerical examples. The obtained results are in good agreement with those in the existing studies and those calculated by the finite element software ANSYS. Moreover, the effects of the material gradient, taper ratio, slenderness ratio, and heightspan ratio on vibration behaviors are discussed.     

Initial value method for free vibration of axially loaded functionally graded Timoshenko beams with nonuniform cross section

Free vibration of nonuniform axially functionally graded Timoshenko beams subjected to combined axially tensile or compressive loading is studied. An emphasis is placed on the effect of tip and distributed axial loads on the natural frequencies and mode shapes for an inhomogeneous cantilever beam including material inhomogeneity and geometric non-uniform cross section. The initial value method is developed to determine the natural frequencies. The method’s effectiveness is verified by comparing our results with previous ones for special cases. Natural frequencies of standing/hanging Timoshenko beams are calculated for four different cross sections. The influences of shear rigidity, taper ratio, gradient index, tip force, and axially distributed loading on the natural frequencies of clamped-free beams are discussed. Material inhomogeneity and geometric non-uniform cross-section strongly affect higher-order vibration frequencies and mode shapes.     

New shape functions for non-uniform curved Timoshenko beams with arbitrarily varying curvature using basic displacement functions

Basic Displacement Functions (BDFs) are introduced and derived using the basic principles of structural mechanics. BDFs are then utilized to obtain new shape functions for arbitrarily curved non-uniform beams, including the effects of shear deformation and extensibility of neutral axis. The main virtue of the proposed shape functions is their susceptibility to the variations in cross-sectional area, moment of inertia, and curvature along the axis of the beam element. Competence of the present method in static and free vibration analyses, as well as its applicability to the special case of straight Timoshenko beam has been verified through several numerical examples.     

Out-of-plane responses of a circular curved Timoshenko beam due to a moving load

For the cases of using the finite curved beam elements and taking the effects of both the shear deformation and rotary inertias into consideration, the literature regarding either free or forced vibration analysis of the curved beams is rare. Thus, this paper tries to determine the dynamic responses of a circular curved Timoshenko beam due to a moving load using the curved beam elements. By taking account of the effect of shear deformation and that of rotary inertias due to bending and torsional vibrations, the stiffness matrix and the mass matrix of the curved beam element were obtained from the force–displacement relations and the kinetic energy equations, respectively. Since all the element property matrices for the curved beam element are derived based on the local polar coordinate system (rather than the local Cartesian one), their coefficients are invariant for any curved beam element with constant radius of curvature and subtended angle and one does not need to transform the property matrices of each curved beam element from the local coordinate system to the global one to achieve the overall property matrices for the entire curved beam structure before they are assembled. The availability of the presented approach has been verified by both the existing analytical solutions for the entire continuum curved beam and the numerical solutions for the entire discretized curved beam composed of the conventional straight beam elements based on either the consistent-mass model or the lumped-mass model. In addition to the typical circular curved beams, a hybrid curved beam composed of one curved-beam segment and two identical straight-beam segments subjected to a moving load was also studied. Influence on the dynamic responses of the curved beams of the slenderness ratio, moving-load speed, shear deformation and rotary inertias was investigated.     

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

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

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

基于一阶非线性梁理论和物理中面概念,导出了纵横向载荷作用下功能梯度材料(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.     

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.     

Stochastic response of an axially loaded composite Timoshenko beam exhibiting bending–torsion coupling

A spectral finite element method is proposed to investigate the stochastic response of an axially loaded composite Timoshenko beam with solid or thin-walled closed section exhibiting bending–torsion materially coupling under the stochastic excitations with stationary and ergodic properties. The effects of axial force, shear deformation (SD) and rotary inertia (RI) as well as bending–torsion coupling are considered in the present study. First, the damped general governing differential equations of motion of an axially loaded composite Timoshenko beam are derived. Then, the spectral finite element formulation is developed in the frequency domain using the dynamic shape functions based on the exact solutions of the governing equations in undamped free vibration, which is used to compute the mean square displacement response of axially loaded composite Timoshenko beams. Finally, the proposed method is illustrated by its application to a specific example to investigate the effects of bending–torsion coupling, axial force, SD and RI on the stochastic response of the composite beam.     

Natural frequencies of composite beams with a variable fiber volume fraction including rotary inertia and shear deformation

The present study is concerned with the vibration analysis of symmetric composite beams with a variable fiber volume fraction through thickness. First-order shear deformation and rotary inertia have been included in the analysis. The solution procedure is applicable to arbitrary boundary conditions. Continuous gradation of the fiber volume fraction is modeled in the form of an m-th power polynomial of the coordinate axis in the thickness direction of the beam. By varying the fiber volume fraction within the symmetric composite beam to create a functionally graded material （FGM）, certain vibration characteristics are affected. Results are presented to demonstrate the effects of shear deformation, fiber volume fraction, and boundary conditions on the natural frequencies and mode shapes of composite beams.     

Size-dependent vibration of functionally graded curved microbeams based on the modified strain gradient elasticity theory

On the basis of the modified strain gradient elasticity theory, the free vibration characteristics of curved microbeams made of functionally graded materials (FGMs) whose material properties vary in the thickness direction are investigated. A size-dependent first-order shear deformation beam model is developed containing three internal material length scale parameters to incorporate small-scale effect. Through Hamilton’s principle, the higher-order governing equations of motion and boundary conditions are derived. Natural frequencies of FGM curved microbeams corresponding to different mode numbers are evaluated for over a wide range of material property gradient index, dimensionless length scale parameter and aspect ratio. Moreover, the results obtained via the present non-classical first-order shear deformation beam model are compared with those of degenerated beam models based on the modified couple stress and the classical theories. It is found that the difference between the natural frequencies predicted by the various beam models is more significant for lower values of dimensionless length scale parameter and higher values of mode number.     

Buckling and vibration of axially functionally graded nonuniform beams using differential transformation based dynamic stiffness approach

The free vibration of axially functionally graded (FG) non-uniform beams with different boundary conditions is studied using Differential Transformation (DT) based Dynamic Stiffness approach. This method is capable of modeling any beam (Timoshenko or Euler, centrifugally stiffened or not) whose cross sectional area, moment of Inertia and material properties vary along the beam. The effectiveness of the method is confirmed by comparing the present results with existing closed form solutions and numerical results. In FG beams, flexural rigidity and mass density may take majority of functions including polynomials, trigonometric and exponential functions (converted to polynomial expressions). DT based Dynamic stiffness approach is proved to be a versatile and simple approach compared to many other methods already proposed.     

高效的三维曲梁单元

三维井眼中延伸数千米的三维细长圆截面钢钻柱应力分析问题是一个复杂的力学问题,通常使用有限元数值分析方法对其进行受力分析。而在进行有限元分析时,现有的圆弧曲粱单元和空间直粱单元在几何上都不能很好地模拟三维曲线形状的钻柱。为了确保计算精度．其单元划分势必不能过大,结果是计算时间长,收敛性差。为了解决这一问题,显然必须构建一种新的较有效的曲梁单元。基于自然坐标系,依据圆截面空间曲粱单元节点有6个自由度——3个线位移和3个角位移,利用包含全部刚体位移模式和常应变的形函数,忽略剪切变形,假设变形后的梁轴线的弯曲曲率改变为线性变化,建立起了保证收敛性的具有12个自由度的有初始曲率和挠率的圆截面空间曲梁的有限元模型。为了证明给出的有限元模型的高效性,分析了几个静态问题,并与现有文献中的解析解或数值结果进行了比较。基于所给出的结果,可望该有限元模型可以作为分析三维空间曲粱结构的有效工具。     

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.     

Large deflection of curved elastic beams made of Ludwick type material

The large deflection of an axially extensible curved beam with a rectangular cross-section is investigated. The elastic beam is assumed to satisfy the Euler-Bernoulli postulation and be made of the Ludwick type material. Through reasonably simplified integration, the strain and curvature of the axis of the beam are presented in implicit formulations. The governing equations involving both geometric and material nonlinearities of the curved beam are derived and solved by the shooting method. When the initial curvature of the beam is zero, the curved beam is degenerated into a straight beam,and the predicted results obtained by the present model are consistent with those in the open literature. Numerical examples are further given for curved cantilever and simply supported beams, and the couplings between elongation and bending are found for the curved beams.     

Modified shooting approach to the non-linear periodic forced response of isotropic/composite curved beams

Large amplitude periodic forced vibration of curved beams under periodic excitation is investigated using a three-noded beam element. The element is based on the higher-order shear deformation theory satisfying interlayer continuity of displacements and transverse shear stress, and top-bottom conditions on the latter. The periodic responses are obtained using shooting technique coupled with Newmark time marching and arc length continuation algorithm developed. The second order governing differential equations of motion are solved without transforming to the first order differential equations thereby resulting in a computationally more efficient algorithm. The effects of excitation amplitude, support conditions and beam curvature on the frequency versus response amplitude relation are highlighted. The typical frequency response curves for isotropic and cross-ply laminated curved beams are presented. Phenomenon of strong modal interactions is observed.     

Elasticity solutions for plane anisotropic functionally graded beams

This paper considers the plane stress problem of generally anisotropic beams with elastic compliance parameters being arbitrary functions of the thickness coordinate. Firstly, the partial differential equation, which is satisfied by the Airy stress function for the plane problem of anisotropic functionally graded materials and involves the effect of body force, is derived. Secondly, a unified method is developed to obtain the stress function. The analytical expressions of axial force, bending moment, shear force and displacements are then deduced through integration. Thirdly, the stress function is employed to solve problems of anisotropic functionally graded plane beams, with the integral constants completely determined from boundary conditions. A series of elasticity solutions are thus obtained, including the solution for beams under tension and pure bending, the solution for cantilever beams subjected to shear force applied at the free end, the solution for cantilever beams or simply supported beams subjected to uniform load, the solution for fixed–fixed beams subjected to uniform load, and the one for beams subjected to body force, etc. These solutions can be easily degenerated into the elasticity solutions for homogeneous beams. Some of them are absolutely new to literature, and some coincide with the available solutions. It is also found that there are certain errors in several available solutions. A numerical example is finally presented to show the effect of material inhomogeneity on the elastic field in a functionally graded anisotropic cantilever beam.     

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

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

正交各向异性功能梯度材料开口圆柱壳的自由振动分析

区别于一般圆柱壳,开口圆柱壳沿周向是不封闭的,因此具有四个边界,本文根据轴向梁式振动和轴向曲拱振动特性对各种端部与侧边边界条件下的壳体提出统一的位移振型函数,并根据哈密顿原理建立了材料参数与空间坐标相关的正交各向异性开口圆柱壳的动力变分方程,求出了不同材料属性下开敞圆柱壳固有频率与振型解的一般解析表达式,适用于任意边界条件下不同材料的开敞圆柱壳自由振动分析．