Large amplitude non-linear oscillations of a general conservative system

This paper presents new, approximate analytical solutions to large-amplitude oscillations of a general, inclusive of odd and non-odd non-linearity, conservative single-degree-of-freedom system. Based on the original general non-linear oscillating system, two new systems with odd non-linearity are to be addressed. Building on the approximate analytical solutions of odd non-linear systems developed by the authors earlier, we construct the new approximate analytical solutions to the original general non-linear system by combinatory piecing of the approximate solutions corresponding to, respectively, the two new systems introduced. These approximate solutions are valid for small as well as large amplitudes of oscillation for which the perturbation method either provides inaccurate solutions or is inapplicable. Two examples with excellent approximate analytical solutions are presented to illustrate the great accuracy and simplicity of the new formulation.     

An analytical approximate technique for a class of strongly non-linear oscillators

An analytical approximate technique for large amplitude oscillations of a class of conservative single degree-of-freedom systems with odd non-linearity is proposed. The method incorporates salient features of both Newton's method and the harmonic balance method. Unlike the classical harmonic balance method, accurate analytical approximate solutions are possible because linearization of the governing differential equation by Newton's method is conducted prior to harmonic balancing. The approach yields simple linear algebraic equations instead of non-linear algebraic equations without analytical solution. With carefully constructed iterations, only a few iterations can provide very accurate analytical approximate solutions for the whole range of oscillation amplitude beyond the domain of possible solution by the conventional perturbation methods or harmonic balance method. Three examples including cubic-quintic Duffing oscillators are presented to illustrate the usefulness and effectiveness of the proposed technique.     

多孔功能梯度材料Timoshenko梁的自由振动分析

基于Timoshenko梁理论研究多孔功能梯度材料梁(FGMs)的自由振动问题.首先,考虑多孔功能梯度材料梁的孔隙率模型,建立了两种类型的孔隙分布.其次,基于Timoshenko梁变形理论,给出位移场方程、几何方程和本构方程,利用Hamilton原理推导多孔功能梯度材料梁的自由振动控制微分方程,并进行无量纲化,然后应用微分变换法(DTM)对无量纲控制微分方程及其边界条件进行变换,得到含有固有频率的等价代数特征方程.最后,计算了固定-固定(C-C)、固定-简支(C-S)和简支-简支(S-S)三种不同边界下多孔功能梯度材料梁自由振动的无量纲固有频率.将其退化为均匀材料与已有文献数据结果对照,验证了正确性.讨论了孔隙率、细长比和梯度指数对多孔功能梯度材料梁无量纲固有频率的影响.     

Free vibration of functionally graded material beams with surface-bonded piezoelectric layers in thermal environment

Free vibration of statically thermal postbuckled functionally graded material （FGM） beams with surface-bonded piezoelectric layers subject to both temperature rise and voltage is studied. By accurately considering the axial extension and based on the Euler-Bernoulli beam theory, geometrically nonlinear dynamic governing equations for FGM beams with surface-bonded piezoelectric layers subject to thermo-electro- mechanical loadings are formulated. It is assumed that the material properties of the middle FGM layer vary continuously as a power law function of the thickness coordinate, and the piezoelectric layers are isotropic and homogenous. By assuming that the amplitude of the beam vibration is small and its response is harmonic, the above mentioned non-linear partial differential equations are reduced to two sets of coupled ordinary differential equations. One is for the postbuckling, and the other is for the linear vibration of the beam superimposed upon the postbuckled configuration. Using a shooting method to solve the two sets of ordinary differential equations with fixed-fixed boundary conditions numerically, the response of postbuckling and free vibration in the vicinity of the postbuckled configuration of the beam with fixed-fixed ends and subject to transversely nonuniform heating and uniform electric field is obtained. Thermo-electric postbuckling equilibrium paths and characteristic curves of the first three natural frequencies versus the temperature, the electricity, and the material gradient parameters are plotted. It is found that the three lowest frequencies of the prebuckled beam decrease with the increase of the temperature, but those of a buckled beam increase monotonically with the temperature rise. The results also show that the tensional force produced in the piezoelectric layers by the voltage can efficiently increase the critical buckling temperature and the natural frequency.     

Approximate solutions of axially moving viscoelastic beams subject to multi-frequency excitations

Approximate solutions of axially moving viscoelastic beams subject to multi-frequency excitations are presented. A non-linear partial-differential equation governing the transverse vibration of the beams is derived from Newton's second law, the Kelvin constitutive relationship, and the Lagrangian strain. Based on 1-term Galerkin's truncation, the governing equation is reduced to an ordinary differential equation. Three cases, including superharmonic resonance case, subharmonic resonance, and combination resonance are studied. The approximate solutions of the transverse vibration of the beams are obtained. Numerical results show that the approximate solutions are in good agreement with numerical results.     

Thermal-induced snap-through buckling of simply-supported functionally graded beams

The instability of functionally graded material (FGM) structures is one of the major threats to their service safety in engineering applications. This paper aims to clarify a long-standing controversy on the thermal instability type of simply-supported FGM beams. First, based on the Euler-Bernoulli beam theory and von Kármán geometric nonlinearity, a nonlinear governing equation of simply-supported FGM beams under uniform thermal loads by Zhang's two-variable method is formulated. Second, an approximate analytic solution to the nonlinear integro-differential boundary value problem under a thermal-induced inhomogeneous force boundary condition is obtained by using a semiinverse method when the coordinate axis is relocated to the bending axis (physical neutral plane), and then the analytical predictions are verified by the differential quadrature method (DQM). Finally, based on the free energy theorem, it is revealed that the symmetry breaking caused by the material inhomogeneity can make the simply-supported FGM beam under uniform thermal loads occur snap-through postbuckling only in odd modes; furthermore, the nonlinear critical load of thermal buckling varies non-monotonically with the functional gradient index due to the stretching-bending coupling effect. These results are expected to provide new ideas and references for the design and regulation of FGM structures.     

Numerical method for unsteady compressible boundary layers driven by compression or expansion wave

A numerical analysis is presented for the unsteady compressible laminar boundary layer driven by a compression or expansion wave. Approximate or series expansion methods have been used for the problems because of the characteristics of the governing equations, such as non-linearity, coupling with the thermal boundary layer equation and initial conditions. Here a transformation of the governing equations and the numerical linearization technique are introduced to deal with the difficulties. First, the governing equations are transformed for the initial conditions by Howarth and semisimilarity variables. These transformations reduce the number of independent variables from three to two and the governing equations from partial to ordinary differential equations at the initial point. Next, the numerical linearization technique is introduced for the non-linearity and the coupling with the thermal boundary layer equation. Because the non-linear terms are linearized without sacrifice of numerical accuracy, the solutions can be obtained without numerical iterations. Therefore the exact numerical solution, not approximate or series expansion, can be obtained. Compared with the approximate or series expansion method, this method is much improved. Results are compared with the series expansion solutions.     

Two-dimensional elasticity solution for bending of functionally graded beams with variable thickness

This paper studies the stress and displacement distributions of functionally graded beam with continuously varying thickness, which is simply supported at two ends. The Young’s modulus is graded through the thickness following the exponential-law and the Poisson’s ratio keeps constant. On the basis of two-dimensional elasticity theory, the general expressions for the displacements and stresses of the beam under static loads, which exactly satisfy the governing differential equations and the simply supported boundary conditions at two ends, are analytically derived out. The unknown coefficients in the solutions are approximately determined by using the Fourier sinusoidal series expansions to the boundary conditions on the upper and lower surfaces of the beams. The effect of Young’s modulus varying rules on the displacements and stresses of functionally graded beams is investigated in detail. The two-dimensional elasticity solution obtained can be used to assess the validity of various approximate solutions and numerical methods for the aforementioned functionally graded beams.     

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.     

Nonlinear vibration analysis of functionally graded beams considering the influences of the rotary inertia of the cross section and neutral surface position

In the paper work, the nonlinear vibration response of functionally graded (FG) Euler–Bernoulli beam resting on elastic foundation is studied. Based on von Kármán’s geometric nonlinearity, the partial differential governing equations describing the nonlinear vibration of FG Euler–Bernoulli beam are derived from Hamilton’s principle and are reduced to an ordinary nonlinear differential equation with quadratic and cubic nonlinear terms via Galerkin’s procedure. Due to unsymmetrical material variation along the thickness of FG beam, the neutral surface concept is proposed to remove the stretching and bending coupling effect and the rotary inertia of the cross section is incorporated to obtain an analytical solution. Numerical results are presented to show the effects of the nonlocal parameters and vibration amplitude on the frequency responses. This results may be useful in design and engineering applications.     

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

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

Thermoelastically coupled axisymmetric nonlinear vibration of shallow spherical and conical shells

The problem of axisymmetric nonlinear vibration for shallow thin spherical and conical shells when temperature and strain fields are coupled is studied. Based on the large deflection theories of yon Ktirrntin and the theory of thermoelusticity, the whole governing equations and their simplified type are derived. The time-spatial variables are separated by Galerkin ‘ s technique, thus reducing the governing equations to a system of time-dependent nonlinear ordinary differential equation. By means of regular perturbation method and multiple-scales method, the first-order approximate analytical solution for characteristic relation of frequency vs amplitude parameters along with the decay rate of amplitude are obtained, and the effects of different geometric parameters and coupling factors us well us boundary conditions on thermoelustically coupled nonlinear vibration behaviors are discussed.     

Free vibration of non-uniform axially functionally graded beams using the asymptotic development method

The asymptotic development method is applied to analyze the free vibration of non-uniform axially functionally graded(AFG) beams, of which the governing equations are differential equations with variable coefficients. By decomposing the variable flexural stiffness and mass per unit length into reference invariant and variant parts, the perturbation theory is introduced to obtain an approximate analytical formula of the natural frequencies of the non-uniform AFG beams with different boundary conditions.Furthermore, assuming polynomial distributions of Young's modulus and the mass density, the numerical results of the AFG beams with various taper ratios are obtained and compared with the published literature results. The discussion results illustrate that the proposed method yields an effective estimate of the first three order natural frequencies for the AFG tapered beams. However, the errors increase with the increase in the mode orders especially for the cases with variable heights. In brief, the asymptotic development method is verified to be simple and efficient to analytically study the free vibration of non-uniform AFG beams, and it could be used to analyze any tapered beams with an arbitrary varying cross width.     

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.     

An analytical study on the nonlinear vibration of functionally graded beams

Nonlinear vibration of beams made of functionally graded materials (FGMs) is studied in this paper based on Euler-Bernoulli beam theory and von Kármán geometric nonlinearity. It is assumed that material properties follow either exponential or power law distributions through thickness direction. Galerkin procedure is used to obtain a second order nonlinear ordinary equation with quadratic and cubic nonlinear terms. The direct numerical integration method and Runge-Kutta method are employed to find the nonlinear vibration response of FGM beams with different end supports. The effects of material property distribution and end supports on the nonlinear dynamic behavior of FGM beams are discussed. It is found that unlike homogeneous beams, FGM beams show different vibration behavior at positive and negative amplitudes due to the presence of quadratic nonlinear term arising from bending-stretching coupling effect.     

Construction of approximate analytical solutions to strongly nonlinear damped oscillators

An analytical approximate method for strongly nonlinear damped oscillators is proposed. By introducing phase and amplitude of oscillation as well as a bookkeeping parameter, we rewrite the governing equation into a partial differential equation with solution being a periodic function of the phase. Based on combination of the Newton’s method with the harmonic balance method, the partial differential equation is transformed into a set of linear ordinary differential equations in terms of harmonic coefficients, which can further be converted into systems of linear algebraic equations by using the bookkeeping parameter expansion. Only a few iterations can provide very accurate approximate analytical solutions even if the nonlinearity and damping are significant. The method can be applied to general oscillators with odd nonlinearities as well as even ones even without linear restoring force. Three examples are presented to illustrate the usefulness and effectiveness of the proposed method.     

功能梯度板的非线性动力分析

非线性材料功能梯度板件的动力分析是属于在数学方程上同时具有变系数、非线性、非定常特征的固体力学问题.文中首先将问题的变系数非线性偏微分方程组转化为各向异性常系数非线性常微分方程,然后用小参数法求得解析解,适用于各种形状、边界及功能梯度分布的板件非线性弹性振动分析.     

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.     

热荷载作用下Timoshenko功能梯度夹层梁的静态响应

在精确考虑轴线伸长和一阶横向剪切变形的基础上建立了Timoshenko功能梯度夹层梁在热载荷作用下的几何非线性控制方程.采用打靶法数值求解所得强非线性边值问题,获得了两端固支功能梯度夹层梁在横向非均匀升温作用下的静态热过屈曲和热弯曲变形数值解.分析了功能梯度材料参数变化、不同表层厚度和升温参数对夹层梁弯曲变形、拉-弯耦...     

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.