Dynamical analysis of axially moving plate by finite difference method

The complex natural frequencies for linear free vibrations and bifurcation and chaos for forced nonlinear vibration of axially moving viscoelastic plate are investigated in this paper. The governing partial differential equation of out-of-plane motion of the plate is derived by Newton’s second law. The finite difference method in spatial field is applied to the differential equation to study the instability due to flutter and divergence. The finite difference method in both spatial and temporal field is used in the analysis of a nonlinear partial differential equation to detect bifurcations and chaos of a nonlinear forced vibration of the system. Numerical results show that, with the increasing axially moving speed, the increasing excitation amplitude, and the decreasing viscosity coefficient, the equilibrium loses its stability and bifurcates into periodic motion, and then the periodic motion becomes chaotic motion by period-doubling bifurcation.     

Dynamical Stability of Viscoelastic Column with Fractional Derivative Constitutive Relation

ＩｎｔｒｏｄｕｃｔｉｏｎＭｏｒｅａｎｄｍｏｒｅｐｒｏｂｌｅｍｓｉｎｓｃｉｅｎｃｅａｎｄｅｎｇｉｎｅｅｒｉｎｇｈａｖｅｃｏｎｃｅｒｎｅｄｔｈｅｓｔａｂｉｌｉｔｙａｎａｌｙｓｉｓｏｆｔｈｅｖｉｓｃｏｅｌａｓｔｉｃｓｔｒｕｃｔｕｒｅｓｄｕｅｔｏｔｈｅｉｒｂｒｏａｄａｐｐｌｉｃａｔｉｏｎｓ .Ｉｎ [1]ＣｅｄｅｒｂａｕｍａｎｄＭｏｎｄａｐｐｌｉｅｄｔｈｅｍｕｌｔｉｐｌｅ＿ｓｃａｌｅｓｍｅｔｈｏｄｔｏｔｒｅａｔｓｔａｂｉｌｉｔｙｏｆａｖｉｓｃｏｅｌａｓｔｉｃｃｏｌｕｍｎｕｎｄｅｒａｐｅｒｉｏｄｉｃａｘｉａｌｌ…     

Complex-mode Galerkin approach in transverse vibration of an axially accelerating viscoelastic string

Under the consideration of harmonic fluctuations of initial tension and axially velocity, a nonlinear governing equation for transverse vibration of an axially accelerating string is set up by using the equation of motion for a 3-dimensional deformable body with initial stresses. The Kelvin model is used to describe viscoelastic behaviors of the material. The basis function of the complex-mode Galerkin method for axially accelerating nonlinear strings is constructed by using the modal function of linear moving strings with constant axially transport velocity. By the constructed basis functions, the application of the complex-mode Galerkin method in nonlinear vibration analysis of an axially accelerating viscoelastic string is investigated. Numerical results show that the convergence velocity of the complex-mode Galerkin method is higher than that of the real-mode Galerkin method for a variable coefficient gyroscopic system.     

非线性粘弹性板的失稳条件

研究了给定面内周期激励作用下简支各向同性均匀粘弹性板平衡失稀问题,板的材料特性由Leaderman非线性本构关系描述,将板的动力学方程进行（Galerkin截断得到简化数学模型为弱非线性系统,采用平均法得到系统的平均化方程,对平均化方程进行稳定性分析得到了板平衡失稳的解析条件,对原系统用数值仿真进行研究,数值结果表明,随着激励幅值的增加或粘弹性材料系数的减少,系统平衡点推失稳,激励幅值和粘弹性材料系数的临界值均与解析结果接近。     

多自由度参激系统稳定性的数值分析

提出多自由度周期参激系统稳定性的数值直接法。通过将扰动方程表示成状态方程形式,再根据Flo-quet理论将扰动解表示成指数特征分量与周期分量之积,并将其周期分量与系统周期系数展成Fourier级数,导出一系列代数方程,建立矩阵特征值问题,从而由数值求解特征值可直接确定参激系统的稳定性。该方法可用于一般周期参激阻尼系统,特征值矩阵不含逆子阵。应用于斜拉索在支座周期运动激励下的参激振动不稳定性分析,数值结果表明该方法的有效性。     

GENERALIZED FINITE SPECTRAL METHOD FOR 1D BURGERS AND KDV EQUATIONS

A generalized finite spectral method is proposed. The method is of high-order accuracy. To attain high accuracy in time discretization, the fourth-order Adams-Bashforth-Moulton predictor and corrector scheme was used. To avoid numerical oscillations caused by the dispersion term in the KdV equation, two numerical techniques were introduced to improve the numerical stability. The Legendre, Chebyshev and Her-mite polynomials were used as the basis functions. The proposed numerical scheme is validated by applications to the Burgers equation (nonlinear convection- diffusion problem) and KdV equation (single solitary and 2-solitary wave problems), where analytical solutions are available for comparison. Numerical results agree very well with the corresponding analytical solutions in all cases.     

Mathematical modeling for dynamic stability of sandwich beam with variable mechanical properties of core

The paper is devoted to mathematical modelling of static and dynamic stability of a simply supported three-layered beam with a metal foam core. Mechanical properties of the core vary along the vertical direction. The field of displacements is formulated using the classical broken line hypothesis and the proposed nonlinear hypothesis that generalizes the classical one. Using both hypotheses, the strains are determined as well as the stresses of each layer. The kinetic energy, the elastic strain energy, and the work of load are also determined. The system of equations of motion is derived using Hamilton’s principle. Finally, the system of three equations is reduced to one equation of motion, in particular, the Mathieu equation. The Bubnov-Galerkin method is used to solve the system of equations of motion, and the Runge-Kutta method is used to solve the second-order differential equation. Numerical calculations are done for the chosen family of beams. The critical loads, unstable regions, angular frequencies of the beam, and the static and dynamic equilibrium paths are calculated analytically and verified numerically. The results of this study are presented in the forms of figures and tables.     

矩形脉冲作用下复合材料层合扁球壳的冲击屈曲分析

研究了计及横向剪切的复合材料层合扁球壳在矩形脉冲载荷作用下的非线性动力屈曲问题;采用Galerkin方法得到以顶点挠度表达的动力响应方程,并用Runge-Kutta方法进行数值求解,应用Budiansky-Roth准则(简称B-R准则)确定冲击屈曲的临界荷载;讨论了壳体几何尺寸和物理参数对复合材料层合扁球壳冲击屈曲的影响;数值算例表明,该方法是可行的.     

扁球面网壳的混沌运动研究

在圆形三向网架非线性动力学基本方程的基础上,用拟壳法给出了圆底扁球面三向网壳的非线性动力学基本方程．在固定边界条件下,引入了异于等厚度壳的无量纲量,对基本方程和边界条件进行无量纲化,通过Galerkin作用得到了一个含二次、三次的非线性动力学方程．为求Melnikov函数,对一类非线性动力系统的自由振动方程进行了求解,得到了此类问题的准确解．在无激励情况下,讨论了稳定性问题．在外激励情况下,通过求Melnikov函数,给出了可能发生混沌运动的条件．通过数字仿真绘出了平面相图,证实了混沌运动的存在．     

Combination,principal parametric and internal resonances of an accelerating beam under two frequency parametric excitation

This study analyses the nonlinear transverse vibration of an axially moving beam subject to two frequency excitation. Focus has been made on simultaneous resonant cases i.e. principal parametric resonance of first mode and combination parametric resonance of additive type involving first two modes in presence of internal resonance. By adopting the direct method of multiple scales, the governing nonlinear integro-partial differential equation for transverse motion is reduced to a set of nonlinear first order ordinary partial differential equations which are solved either by means of continuation algorithm or via direct time integration. Specifically, the frequency response plots and amplitude curves, their stability and bifurcation are obtained using continuation algorithm. Numerical results reveal the rich and interesting nonlinear phenomena that have not been presented in the existent literature on the nonlinear dynamics of axially moving systems.     

Peristaltic flow of a Johnson-Segalman fluid through a deformable tube

To understand theoretically the flow properties of physiological fluids we have considered as a model the peristaltic motion of a Johnson–Segalman fluid in a tube with a sinusoidal wave traveling down its wall. The perturbation solution for the stream function is obtained for large wavelength and small Weissenberg number. The expressions for the axial velocity, pressure gradient, and pressure rise per wavelength are also constructed. The general solution of the governing nonlinear partial differential equation is given using a transformation method. The numerical solution is also obtained and is compared with the perturbation solution. Numerical results are demonstrated for various values of the physical parameters of interest.      

粘弹性桩的混沌运动分析

研究了在轴向载荷和周期性横向载荷共同作用下非线性粘弹性嵌岩桩的混沌运动情况。假定桩和土体分别满足Leaderman非线性粘弹性和线性粘弹性本构关系，得到的运动方程为非线性偏微分．积分方程；利用Galerkin方法将方程简化为非线性常微分一积分方程，同时利用非线性动力系统中的数值方法，进行了数值计算，得到了不同载荷参数、几何参数、材料参数时粘弹性桩发生周期运动、多周期运动及混沌运动的时程曲线、相图、功率谱、Poincare截面图，同时得到了挠度-载荷、挠度-几何参数、挠度-材料参数等分叉图，考察了各种参数的影响。数值结果表明非线性粘弹性桩在一定的条件下可以通过倍周期分叉的方式进入混沌运动状态，且桩的载荷参数、几何参数、材料参数对其运动状态有较大的影响。     

爆炸冲击下复合材料层合扁球壳的动力屈曲

研究计及横向剪切的复合材料层合扁球壳在爆炸冲击载荷作用下的非线性轴对称动力屈曲问题。通过在复合材料层合扁球壳非线性稳定性的基本方程中增加横向转动惯量项并引入R.H.Cole理论的爆炸冲击力,得到爆炸冲击下复合材料层合扁球壳的动力控制方程,应用Galerkin方法得到用顶点挠度表达的爆炸冲击动力响应方程,并采用Runge-Kutta方法进行数值求解,采用Budiansky-Roth准则确定冲击屈曲的临界载荷,讨论了壳体几何尺寸对复合材料层合扁球壳冲击屈曲的影响;数值算例表明,此方法是可行的。     

A Newton type iterative method for heat-conduction inverse problems

An inverse problem for identification of the coefficient in heat-conduction equation is considered. After reducing the problem to a nonlinear ill-posed operator equation, Newton type iterative methods are considered. The implicit iterative method is applied to the linearized Newton equation, and the key step in the process is that a new reasonable a posteriori stopping rule for the inner iteration is presented. Numerical experiments for the new method as well as for Tikhonov method and Bakushikskii method are given, and these results show the obvious advantages of the new method over the other ones.     

A Newton type iterative inverse method for heat-conduction problems

An inverse problem for identification of the coefficient in heat-conduction equation is considered. After reducing the problem to a nonlinear ill-posed operator equation, Newton type iterative methods are considered. The implicit iterative method is applied to the linearized Newton equation, and the key step in the process is that a new reasonable a posteriori stopping rule for the inner iteration is presented. Numerical experiments for the new method as well as for Tikhonov method and Bakushikskii method are given, and these results show the obvious advantages of the new method over the other ones.     

Effects of thermo-magneto-electro nonlinearity characteristics on the stability of functionally graded piezoelectric beam

Due to the increasing interests in using functionally graded piezoelectric materials(FGPMs) in the design of advanced micro-electro-mechanical systems, it is important to understand the stability behaviors of the FGPM beams. In this study, considering the effects of geometrical nonlinearity, temperature, and electricity in the constitutive relations and the effect of the magnetic field on the FGPM beam, the Euler-Bernoulli beam model is adopted, and the nonlinear governing equation of motion is derived via Hamilton's principle. A perturbation method, which can decompose the deflection into static and dynamic components, is utilized to linearize the nonlinear governing equation. Then,a dynamic stability analysis is carried out, and the approximate analytical solutions for the nonlinear frequency and boundary frequencies of the unstable region are obtained.Numerical examples are performed to verify the present analysis. The effects of the static deflection, the static load factor, the temperature change, and the magnetic field flux on the stability behaviors of the FGPM beam are discussed. From the proposed analytical solutions and numerical results, one can easily and clearly find the effects of various controlled parameters, such as geometric and physical properties of the system, on the mechanical behaviors of structures, and the conclusions are very important and useful for the design of micro-devices.     

Subharmonic resonance of a clamped-clamped buckled beam with 1:1 internal resonance under base harmonic excitations

The subharmonic resonance and bifurcations of a clamped-clamped buckled beam under base harmonic excitations are investigated. The nonlinear partial integrodifferential equation of the motion of the buckled beam with both quadratic and cubic nonlinearities is given by using Hamilton's principle. A set of second-order nonlinear ordinary differential equations are obtained by spatial discretization with the Galerkin method. A high-dimensional model of the buckled beam is derived, concerning nonlinear coupling. The incremental harmonic balance (IHB) method is used to achieve the periodic solutions of the high-dimensional model of the buckled beam to observe the nonlinear frequency response curve and the nonlinear amplitude response curve, and the Floquet theory is used to analyze the stability of the periodic solutions. Attention is focused on the subharmonic resonance caused by the internal resonance as the excitation frequency near twice of the first natural frequency of the buckled beam with/without the antisymmetric modes being excited. Bifurcations including the saddle-node, Hopf, perioddoubling, and symmetry-breaking bifurcations are observed. Furthermore, quasi-periodic motion is observed by using the fourth-order Runge-Kutta method, which results from the Hopf bifurcation of the response of the buckled beam with the anti-symmetric modes being excited.     

Symplectic analysis for wave propagation in one-dimensional nonlinear periodic structures

The wave propagation problem in the nonlinear periodic mass-spring structure chain is analyzed using the symplectic mathematical method. The energy method is used to construct the dynamic equation, and the nonlinear dynamic equation is linearized using the small parameter perturbation method. Eigen-solutions of the symplectic matrix are used to analyze the wave propagation problem in nonlinear periodic lattices. Nonlinearity in the mass-spring chain, arising from the nonlinear spring stiffness effect, has profound effects on the overall transmission of the chain. The wave propagation characteristics are altered due to nonlinearity, and related to the incident wave intensity, which is a genuine nonlinear effect not present in the corresponding linear model. Numerical results show how the increase of nonlinearity or incident wave amplitude leads to closing of transmitting gaps. Comparison with the normal recursive approach shows effectiveness and superiority of the symplectic method for the wave propagation problem in nonlinear periodic structures.     

The vibration and stability analysis of moderate thick plates by the method of lines

The method of lines based on Hu Hai-chang's theory for the vibration and stability of moderate thick plates is developed. The standard nonlinear ordinary differential equation (ODE) system for natural frequencies and critical load is given by use of ODE techniques, and then any indicated eigenvalue could be obtained directly from ODE solver by employing the so-called initial eigenfunction technique instead of the mode orthogonality condition.Numerical examples show that the present method is very effective and reliable.     

SYMPLECTIC STRUCTURE OF POISSON SYSTEM

When the Poisson matrix of Poisson system is non-constant, classical symplectic methods, such as symplectic Runge-Kutta method, generating function method, cannot preserve the Poisson structure. The non-constant Poisson structure was transformed into the symplectic structure by the nonlinear transform.Arbitrary order symplectic method was applied to the transformed Poisson system. The Euler equation of the free rigid body problem was transformed into the symplectic structure and computed by the midpoint scheme. Numerical results show the effectiveness of the nonlinear transform.