Transverse nonlinear dynamics of axially accelerating viscoelastic beams based on 4-term Galerkin truncation

This paper investigates bifurcation and chaos in transverse motion of axially accelerating viscoelastic beams. The Kelvin model is used to describe the viscoelastic property of the beam material, and the Lagrangian strain is used to account for geometric nonlinearity due to small but finite stretching of the beam. The transverse motion is governed by a nonlinear partial-differential equation. The Galerkin method is applied to truncate the partial-differential equation into a set of ordinary differential equations. When the Galerkin truncation is based on the eigenfunctions of a linear non-translating beam subjected to the same boundary constraints, a computation technique is proposed by regrouping nonlinear terms. The scheme can be easily implemented in practical computations. When the transport speed is assumed to be a constant mean speed with small harmonic variations, the Poincaré map is numerically calculated based on 4-term Galerkin truncation to identify dynamical behaviors. The bifurcation diagrams are present for varying one of the following parameter: the axial speed fluctuation amplitude, the mean axial speed and the beam viscosity coefficient, while other parameters are unchanged.     

Two-dimensional nonlinear dynamics of an axially moving viscoelastic beam with time-dependent axial speed

In the present study, the coupled nonlinear dynamics of an axially moving viscoelastic beam with time-dependent axial speed is investigated employing a numerical technique. The equations of motion for both the transverse and longitudinal motions are obtained using Newton’s second law of motion and the constitutive relations. A two-parameter rheological model of the Kelvin–Voigt energy dissipation mechanism is employed in the modelling of the viscoelastic beam material, in which the material time derivative is used in the viscoelastic constitutive relation. The Galerkin method is then applied to the coupled nonlinear equations, which are in the form of partial differential equations, resulting in a set of nonlinear ordinary differential equations (ODEs) with time-dependent coefficients due to the axial acceleration. A change of variables is then introduced to this set of ODEs to transform them into a set of first-order ordinary differential equations. A variable step-size modified Rosenbrock method is used to conduct direct time integration upon this new set of first-order nonlinear ODEs. The mean axial speed and the amplitude of the speed variations, which are taken as bifurcation parameters, are varied, resulting in the bifurcation diagrams of Poincaré maps of the system. The dynamical characteristics of the system are examined more precisely via plotting time histories, phase-plane portraits, Poincaré sections, and fast Fourier transforms (FFTs).     

Study on vibration and stability of an axially translating viscoelastic Timoshenko beam: Non-transforming spectral element analysis

Engineering systems, such as rolled steel beams, chain and belt drives and high-speed paper, can be modeled as axially translating beams. This article scrutinizes vibration and stability of an axially translating viscoelastic Timoshenko beam constrained by simple supports and subjected to axial pretension. The viscoelastic form of general rheological model is adopted to constitute the material of the beam. The partial differential equations governing transverse motion of the beam are derived from the extended form of Hamilton's principle. The non-transforming spectral element method (NTSEM) is applied to transform the governing equations into a set of ordinary differential equations. The formulation is similar to conventional FFT-based spectral element model except that Daubechies wavelet basis functions are used for temporal discretization. Influences of translating velocities, axial tensile force, viscoelastic parameter, shear deformation, beam model and boundary condition types are investigated on the underlying dynamic response and stability via the NTSEM and demonstrated via numerical simulations.     

Non-linearly parametric resonances of an axially moving viscoelastic sandwich beam with time-dependent velocity

Non-linearly parametric resonances of an axially moving viscoelastic sandwich beam are investigated in this paper. The beam is moving with a time-dependent velocity, namely a harmonically varied velocity about the mean velocity. The partial differential equation is discretized into nonlinear ordinary differential equations via the method of Galerkin truncation and then the steady-state response is obtained using the method of multiple scales, an approximate analytical method. The tuning equations are obtained by eliminating secular terms and the amplitude of the vibration is derived from the tuning equations expressed in polar form, and two bifurcation points are obtained as well. Additionally, the stability conditions of trivial and nontrivial solutions are analyzed using the Routh–Hurwitz criterion. Eventually, the effects of various parameters such as the thickness of core layer, mean velocity, initial tension, and the amplitude of axially moving velocity on amplitude–frequency response curves and unstable regions are investigated.     

Transverse vibration of viscoelastic sandwich beam with time-dependent axial tension and axially varying moving velocity

Nonlinearly parametric resonances of axially accelerating moving viscoelastic sandwich beams with time-dependent tension are investigated in this paper. Based on the Kelvin differential constitutive equation, the controlling equation of the transverse vibration of a beam with large deflection is established. The system has been subjected to a time varying velocity and a harmonic axial tension. Here the governing equation of motion contains linear parametric terms and two frequencies, one is the frequency of axially moving velocity and the other one is the frequency of varying tension. The method of multiple scales is applied directly to the governing equation to obtain the complex eigenfunctions and natural frequencies of the system. The elimination of secular terms leads to the steady-state response and amplitude of vibrations. The influence of various parameters such as initial tension on natural frequencies and the amplitude of axial fluctuation, the phase angle between the two frequencies on response curves has been investigated for two different resonance conditions. With the help of numerical results, it has been shown that by using suitable initial tension, the amplitude of axial fluctuation, the phase angle, the vibration of the sandwich beam can be significantly controlled.     

小垂度粘弹性索非线性响应及振动主动控制

研究具有初始应力的小垂度粘弹性索的非线性动态响应及振动主动控制。在假定索材料的本构关系为一般微分本构类型的基础上,建立小垂度粘弹性索的运动微分方程;应用Galerkin方法将其转化为可用Runge-Kutta数值积分方法求解的一系列三阶非线性常微分方程。在仅考虑面内的横向振动及忽略非线性的情况下得到了连续状态空间中的状态方程,将状态方程离散为差分方程形式,并用矩阵指数来逐步近似状态转移矩阵;基于二次性能指标的最小化得到了最优的控制力与状态向量。最后通过数值仿真研究说明了粘性参数对索动态响应的影响。     

轴向变速运动粘弹性弦线横向振动的复模态Galerkin方法

在考虑初始张力和轴向速度简谐涨落的情况下,利用含预应力三维变形体的运动方程,建立了轴向变速运动弦线横向振动的非线性控制方程,材料的粘弹性行为由Kelvin模型描述．利用匀速运动线性弦线的模态函数构造了变速运动非线性弦线复模态Galerkin方法的基底函数,并借助构造出来的基底函数研究了复模态Galerkin方法在轴向变速运动粘弹性弦线非线性振动分析中的应用．数值结果表明,复模态Galerkin方法相比实模态Galerkin方法对变系数陀螺系统有较高的收敛速度．     

轴向运动三参数黏弹性梁弱受迫振动的渐近分析

研究了轴向运动三参数黏弹性梁的弱受迫振动．建立了轴向运动三参数黏弹性梁受迫振动的控制方程．使用多尺度法渐近分析了运动梁的稳态响应,导出了解稳定性边界方程、稳态振幅的表达式以及稳态响应非零解的存在条件．依据Routh-Hurwitz定律决定了非线性稳态响应非零解的稳定性．     

轴向运动粘弹性板的横向振动特性

研究了轴向运动粘弹性矩形薄板的动力特性和稳定性问题．从二维粘弹性微分型本构关系出发,建立了轴向运动粘弹性板的运动微分方程．采用微分求积法,对四边简支、一对边简支一对边固支两种边界条件下粘弹性板的无量纲复频率进行了数值计算．分析了薄板的长宽比、无量纲运动速度及材料的无量纲延滞时间对其横向振动及稳定性的影响．     

Control and exponential stabilization for the equation of an axially moving viscoelastic strip

The aim through this work is to suppress the transverse vibrations of an axially moving viscoelastic strip. A controller mechanism (dynamic actuator) is attached at the right boundary to control the undesirable vibrations. The moving strip is modeled as a moving beam pulled at a constant speed through 2 eyelets. The left eyelet is fixed in the sense that there is no transverse displacement (see Figure 1 ). The mathematical model of this system consists of an integro‐partial differential equation describing the dynamic of the strip and an integro‐differential equation describing the dynamic of the actuator. The multiplier method is used to design a boundary control law ensuring an exponential stabilization result.     

非线性弹性梁中的混沌带现象

研究了非线性弹性梁的混沌运动,梁受到轴向载荷的作用。非线性弹性梁的本构方程可用三次多项式表示。计及材料非线性和几何非线性,建立了系统的非线性控制方程。利用非线性Galerkin法,得到微分动力系统。采用Melnikov方法对系统进行分析后发现,当载荷P₀和f满足一定条件时,系统将发生混沌运动,且混沌运动的区域呈现带状。还详尽分析了从次谐分岔到混沌的路径,确定了混沌发生的临界条件。     

Linear and nonlinear vibrations of fractional viscoelastic Timoshenko nanobeams considering surface energy effects

In this paper, the linear and nonlinear vibrations of fractional viscoelastic Timoshenko nanobeams are studied based on the Gurtin–Murdoch surface stress theory. Firstly, the constitutive equations of fractional viscoelasticity theory are considered, and based on the Gurtin–Murdoch model, stress components on the surface of the nanobeam are incorporated into the axial stress tensor. Afterward, using Hamilton's principle, equations governing the two-dimensional vibrations of fractional viscoelastic nanobeams are derived. Finally, two solution procedures are utilized to describe the time responses of nanobeams. In the first method, which is fully numerical, the generalized differential quadrature and finite difference methods are used to discretize the linear part of the governing equations in spatial and time domains. In the second method, which is semi-analytical, the Galerkin approach is first used to discretize nonlinear partial differential governing equations in the spatial domain, and the obtained set of fractional-order ordinary differential equations are then solved by the predictor–corrector method. The accuracy of the results for the linear and nonlinear vibrations of fractional viscoelastic nanobeams with different boundary conditions is shown. Also, by comparing obtained results for different values of some parameters such as viscoelasticity coefficient, order of fractional derivative and parameters of surface stress model, their effects on the frequency and damping of vibrations of the fractional viscoelastic nanobeams are investigated.     

Mixed quadratic-cubic autocatalytic reaction–diffusion equations: Semi-analytical solutions

Semi-analytical solutions for autocatalytic reactions with mixed quadratic and cubic terms are considered. The kinetic model is combined with diffusion and considered in a one-dimensional reactor. The spatial structure of the reactant and autocatalyst concentrations are approximated by trial functions and averaging is used to obtain a lower-order ordinary differential equation model, as an approximation to the governing partial differential equations. This allows semi-analytical results to be obtained for the reaction–diffusion cell, using theoretical methods developed for ordinary differential equations. Singularity theory is used to investigate the static multiplicity of the system and obtain a parameter map, in which the different types of steady-state bifurcation diagrams occur. Hopf bifurcations are also found by a local stability analysis of the semi-analytical model. The transitions in the number and types of bifurcation diagrams and the changes to the parameter regions, in which Hopf bifurcations occur, as the relative importance of the cubic and quadratic terms vary, is explored in great detail. A key outcome of the study is that the static and dynamic stability of the mixed system exhibits more complexity than either the cubic or quadratic autocatalytic systems alone. In addition it is found that varying the diffusivity ratio, of the reactant and autocatalyst, causes dramatic changes to the dynamic stability. The semi-analytical results are show to be highly accurate, in comparison to numerical solutions of the governing partial differential equations.     

Transverse free vibration and stability of axially moving nanoplates based on nonlocal elasticity theory

Transverse dynamical behaviors of axially moving nanoplates which could be used to model the graphene nanosheets or other plate-like nanostructures with axial motion are examined based on the nonlocal elasticity theory. The Hamilton's principle is employed to derive the multivariable coupling partial differential equations governing the transverse motion of the axially moving nanoplates. Subsequently, the equations are transformed into a set of ordinary differential equations by the method of separation of variables. The effects of dimensionless small-scale parameter, axial speed and boundary conditions on the natural frequencies in sub-critical region are discussed by the method of complex mode. Then the Galerkin method is employed to analyze the effects of small-scale parameter on divergent instability and coupled-mode flutter in super-critical region. It is shown that the existence of small-scale parameter contributes to strengthen the stability in the super-critical region, but the stability of the sub-critical region is weakened. The regions of divergent instability and coupled-mode flutter decrease even disappear with an increase in the small-scale parameter. The natural frequencies in sub-critical region show different tendencies with different boundary effects, while the natural frequencies in super-critical region keep constants with the increase of axial speed.     

Dynamical properties of discrete Lotka–Volterra equations

A discrete version of the Lotka–Volterra differential equations for competing population species is analyzed in detail in much the same way as the discrete form of the logistic equation has been investigated as a source of bifurcation phenomena and chaotic dynamics. It is found that in addition to the logistic dynamics – ranging from very simple to manifestly chaotic regimes in terms of governing parameters – the discrete Lotka–Volterra equations exhibit their own brands of bifurcation and chaos that are essentially two-dimensional in nature. In particular, it is shown that the system exhibits “twisted horseshoe” dynamics associated with a strange invariant set for certain parameter ranges.     

Vibration of axially moving hyperelastic beam with finite deformation

In this paper, we study the vibration of an axially moving hyperelastic beam under simply supported condition. The kinematic of the axially moving beam have been described by Eulerian-Lagrangian formulation. In continuum mechanics frame, the finite deformation formula and a higher order shear deformation beam theory are applied to describe the deformation of the axially moving hyperelastic beam. In these formulas the material parameter, shear deformation and the geometric non-linearity have been taken into account. Through the Hamilton principle, the governing equations of nonlinear vibration are obtained, where the transverse vibration is coupled with the longitudinal vibration. When the velocity is a constant, the critical speed and natural frequencies are determined by solving the corresponding linear equations. Meantime, effects of the geometrical and material parameters on the critical speed and natural frequencies have been investigated. Comparisons among the critical velocities of the hyperelastic and Euler linear beam are also made. The results show that the critical velocity of hyperelastic beam is larger than that of linear Euler–Bernoulli beam. For the natural frequencies, we have the same conclusions. Lastly, by the multiple scales method, the leading order analytical solutions of the equilibrium state of axially moving hyperelastic beam in the supercritical regime are obtained. Furthermore the amplitudes of analytical solutions of the hyperelastic beam have been compared with that of linear Euler–Bernoulli beam. The effects of the material and geometrical parameters on the asymptotic solutions and the amplitude has been analyzed.     

微分本构粘弹性轴向运动弦线横向振动分析的差分法

给出了微分本构粘弹性轴向运动弦线横向振动数值仿真的一种差分法．文中建立了具有微分本构的粘弹性运动弦线的横向振动模型;通过对系统的控制方程和本构方程在不同的分数节点离散,得到一种新的差分方法．利用这一方法,弦线振动方程的数值计算过程可以交替地显式进行,且有较小的截断误差和好的数值稳定性．与通用的方法比较,新的方法计算简单、方便．文中利用方程的不变量检验了数值结果的可靠性,并利用这一方法给出了一类弦线模型的参数振动分析．     

Free vibration and stability of tapered Euler–Bernoulli beams made of axially functionally graded materials

The free vibration and stability of axially functionally graded tapered Euler–Bernoulli beams are studied through solving the governing differential equations of motion. Observing the fact that the conventional differential transform method (DTM) does not necessarily converge to satisfactory results, a new approach based on DTM called differential transform element method (DTEM) is introduced which considerably improves the convergence rate of the method. In addition to DTEM, differential quadrature element method of lowest-order (DQEL) is used to solve the governing differential equation, as well. Carrying out several numerical examples, the competency of DQEL and DTEM in determination of free longitudinal and free transverse frequencies and critical buckling load of tapered Euler–Bernoulli beams made of axially functionally graded materials is verified.     

Chaotic vibration analysis of rotating,flexible, continuous shaft-disk system with a rub-impact between the disk and the stator

The chaotic vibration analysis of a rotating flexible continuous shaft-disk system with rub-impact is studied. The system is modeled as a continuous shaft with a rigid disk in its mid-section with Coriolis and centrifugal effects included. The governing partial differential equations of motion are extracted based on the Euler–Bernoulli beam theory. The assumed modes method is used to discretize partial differential equations and the resulting equations are solved via numerical methods. Time series, phase plane portrait, power spectra, Poincaré map, bifurcation diagrams, and Lyapunov exponents are used to analyze the vibration behavior of the system. Initially, the case is investigated in which no Coriolis or centrifugal effects are considered. Then, another case is studied in which these effects are considered. The results confirm the claim that the rub-impact occurs at lower speed ratios due to the Coriolis and centrifugal forcing effects, and that the dynamic behaviors of the system for the two cases are much different as a result of the rub-impact in the second case. Periodic, quasi-periodic, sub-harmonic, and chaotic states can be observed while the appearance or disappearance of the chaos is different. The centrifugal forcing effect plays a greater role than that of the Coriolis force on the incidence of the rub-impact. These results can be useful in identifying the undesirable behaviors in these types of rotating systems.     

Radiation effects on the MHD flow near the stagnation point of a stretching sheet

The radiation effects using the Rosseland approximation on the flow of an incompressible viscous electrically conducting fluid over a stretching sheet near the stagnation point in the presence of uniform transverse magnetic field is studied. The governing equations transform to ordinary differential equation by using suitable similarity transformation and then by a perturbation technique the numerical results for temperature distribution were obtained and discussed graphically.