Bifurcation and chaos of an axially accelerating viscoelastic beam

This paper investigates bifurcation and chaos of an axially accelerating viscoelastic beam. The Kelvin–Voigt model is adopted to constitute the material of the beam. Lagrangian strain is used to account for the beam's geometric nonlinearity. The nonlinear partial–differential equation governing transverse motion of the beam is derived from the Newton second law. The Galerkin method is applied to truncate the governing equation into a set of ordinary differential equations. By use of the Poincaré map, the dynamical behavior is identified based on the numerical solutions of the ordinary differential equations. The bifurcation diagrams are presented in the case that the mean axial speed, the amplitude of speed fluctuation and the dynamic viscoelasticity is respectively varied while other parameters are fixed. The Lyapunov exponent is calculated to identify chaos. From numerical simulations, it is indicated that the periodic, quasi-periodic and chaotic motions occur in the transverse vibrations of the axially accelerating viscoelastic beam.     

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).     

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

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

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.     

Nonlinear analysis of a parametrically excited beam with intermediate support by using Multi-dimensional incremental harmonic balance method

In this paper, a nonlinear Euler-Bernoulli beam under a concentrated harmonic excitation with intermediate nonlinear support is investigated. Continuous expression for the kinetic energy, potential energy and dissipation function are constructed. An energy method based on the Lagrange equation combined with the Galerkin truncation is used for discretizing the governing equation. The Multi-dimensional incremental harmonic balance method (MIHBM) is derived, and the comparisons between the numerical results and the approximate analytical solutions based on the MIHBM verify the excellent accuracy of the MIHBM. The steady state dynamic of the beam is investigated by MIHBM. In order to investigate the energy transmission and understand the vibration response of the Euler-Bernoulli beam, the effects of the key parameters on the dynamic behaviors are studied and discussed, individually. The results show that the amplitude-frequency curves exhibits softening nonlinear behavior in the super-harmonic resonance region, and near resonant region the hardening nonlinear behavior is observed depending on the different parameters. Nonlinear dynamic analysis, such as bifurcation, 3-D frequency spectrum, waveform, frequency spectrum, phase diagram and Poincaré map, are also presented in order to study the influences of the key parameters on the vibration behaviors for the beam in a more accurate manner. In addition, the path to chaotic motion is observed to be through a sequence of the periodic motion and quasi-periodic motion.     

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

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

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

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

非线性粘弹性柱的稳定性和混沌运动

研究了受轴向周期力作用的各向同性简支柱的动力学稳定性。假定粘弹性材料满足Lea-derman非线性本构关系。导出运动方程为非线性偏微分-积分方程,并利用Galerkin方法简化为非线性微分-积分方程。应用平均法进行了稳定性分析,并用数值结果进行验证。数值结果还表明系统可能存在混沌运动。     

Non-linear responses of a one-sided constrained beam with base excitation

The paper discussed the non-linear responses of a buckled beamunder base excitation and constrained by a one-sided motionrestraint. The geometric non-linearity due to axial extensionis taken into account. We apply the Galerkin method to the governingpartial differential equation of the transverse motion to obtaina general model of n degrees of freedom (nDOF). The resultsof dynamic response for the 1DOF mode of a pinned-pinned beamare presented for the system both with and without one-sidedmotion restraint. The regular and irregular motions of the 1DOFmodel for the beam are represented in the forms of time trace,phase plot, bifurcation diagram and power spectra. An in-depthstudy based on an energy approach is done to illustrate thenon-linear responses resulting from the multiplicity of resonantsolutions. It is shown that the effect of the motion restrainton the dynamics of the system is significant.     

轴向变速运动弦线的非线性振动的稳态响应及其稳定性

研究具有几何非线性的轴向运动弦线的稳态横向振动及其稳定性．轴向运动速度为常平均速度与小简谐涨落的叠加．应用Hamilton原理导出了描述弦线横向振动的非线性偏微分方程．直接应用于多尺度方法求解该方程．建立了避免出现长期项的可解性条件．得到了近倍频共振时非平凡稳态响应及其存在条件．给出数值例子说明了平均轴向速度、轴向速度涨落的幅值和频率的影响．应用Liapunov线性化稳定性理论,导出倍频参数共振时平凡解和非平凡解的不稳定条件．给出数值算例说明相关参数对不稳定条件的影响．     

车路耦合非线性振动高阶Galerkin截断研究

将移动车辆模型化为运动的两自由度质量-弹簧-阻尼系统,道路模型化为立方非线性黏弹性地基上的弹性梁,并将路面不平度设定为简谐函数．通过受力分析,建立车路非线性耦合振动高阶偏微分方程．采用高阶Galerkin截断结合数值方法求解耦合系统的动态响应．首次研究不同截断阶数对车路耦合非线性振动动态响应的影响,确定Galerkin截断研究车路耦合振动的收敛性．研究结果表明,对于软土地基的沥青路面,耦合振动的动态响应,需要150阶以上的截断才能达到收敛效果．并通过高阶收敛的Galerkin截断研究了系统参数对车路耦合非线性振动动态响应的影响.     

不可压饱和多孔弹性梁的大挠度非线性数学模型

在孔隙流体仅存在沿梁轴线方向扩散的假定下,建立了微观不可压饱和多孔弹性梁大挠度问题的非线性数学模型．利用Galerkin截断法,研究了固定端不可渗透、自由端可渗透的饱和多孔弹性悬臂梁在自由端突加集中载荷作用下的非线性弯曲,得到了梁骨架的挠度、弯矩以及孔隙流体压力等效力偶等的时间响应和沿轴线的分布．比较了大挠度非线性和小挠度线性理论的结果,揭示了两者间的差异．研究发现大挠度理论的结果小于相应的小挠度理论结果,并且,大挠度理论的结果趋于其稳态值的时间小于相应的小挠度理论结果趋于其稳态值的时间．     

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

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

非线性粘弹性梁的动力学行为

建立了描述受周期荷载作用的均匀粘弹性梁动力学行为的非线性偏微分-积分方程,梁的材料满足Leaderman非线性本构关系,对于两端简支的情形用Galerkin方法进行了2阶截断后,简化为常微分-积分方程,进一步简化为便于进行数值实验的常微分方程,最后用数值方法比较了1阶和2阶截断系统的动力学行为。     

Analytical study on dynamic responses of a curved beam subjected to three-directional moving loads

Considering the warping resistance, inertia force and moving three-directional loads, a more comprehensive set of governing equations for vertical, torsional, radial and axial motions of the curved beam are derived. The analytical solutions for vertical, torsional, radial and axial responses of the curved beam subjected to three-directional moving loads are obtained, using the Galerkin method to discretize the partial differential equations and the modal superposition method to decouple the ordinary differential equations. The analytical results are compared with the numerical integration and a published work to verify the validity of the proposed solutions. Effects of Galerkin truncation terms and damping ratio on solution convergence are also discussed. Considering first-mode and higher-mode truncation respectively, the conditions of resonance and cancellation are analyzed for vertical, torsional, radial and axial motions of the curved beam. Taking a curved bridge under passage of a vehicle as an example, the influences of system parameters, such as vehicle speed, braking acceleration, bridge curve radius, bridge span and bridge deck elastic modulus, on bridge midpoint vibration are explored. The proposed approach and results may be beneficial to enhance understanding the three-directional vehicle-induced dynamic responses of curved bridges. It is shown that when the axial motion, or the multiple moving loads are involved, the first-order truncation are not accurate enough and one should use higher-mode truncation to study the responses of curved beams. In addition, it is necessary to consider damping in the vibration study of curved beams.     

Stability analysis of cantilever carbon nanotubes subjected to partially distributed tangential force and viscoelastic foundation

Dynamic instability of cantilever carbon nanotubes conveying fluid embedded in viscoelastic foundation under a partially distributed tangential force is investigated based on nonlocal elasticity theory and Euler–Bernouli beam theory. The present study has incorporated the effects of nonlocal parameter, Knudsen number, surface effects and magnetic field. And two main parameters have also considered, namely partially distributed tangential force and foundation. It is assumed that viscoelastic foundation has modeled as Kelvin–Voigt, Maxwell and Standard linear solid types. The size-dependent governing equation of transverse vibration is derived using Hamilton’s variational principle and discretized by the Galerkin truncation method. A detailed parameter study is carried out, indicating the stability behavior of the nanotubes. In the light of numerical results, it is shown that variables considered in nondimensional equations have significant effects on natural frequencies and flutter velocities, especially for the foundation distribution length and model as well as the partially distributed tangential force.     

Post-buckling and nonlinear free vibration analysis of geometrically imperfect functionally graded beams resting on nonlinear elastic foundation

In this paper, post-buckling and nonlinear vibration analysis of geometrically imperfect beams made of functionally graded materials (FGMs) resting on nonlinear elastic foundation subjected to axial force are studied. The material properties of FGMs are assumed to be graded in the thickness direction according to a simple power law distribution in terms of the volume fractions of the constituents. The assumptions of a small strain and moderate deformation are used. Based on Euler–Bernoulli beam theory and von-Karman geometric nonlinearity, the integral partial differential equation of motion is derived. Then this partial differential equation (PDE) problem, which has quadratic and cubic nonlinearities, is simplified into an ordinary differential equation (ODE) problem by using the Galerkin method. Finally, the governing equation is solved analytically using the variational iteration method (VIM). Some new results for the nonlinear natural frequencies and buckling load of the imperfect functionally graded (FG) beams such as the effects of vibration amplitude, elastic coefficients of foundation, axial force, end supports and material inhomogeneity are presented for future references. Results show that the imperfection has a significant effect on the post-buckling and vibration response of FG beams.     

粘弹性运动带动力响应分析

基于Kelvin粘弹性材料本构模型及带运动方程,建立了运动带非线性动力学分析模型.基于该模型和Lie群分析方法推导了匀速运动及简谐运动带线性问题的解析解;基于该非线性模型的数值仿真讨论了运动带材料参数、带稳态运动速度、扰动速度对系统动态响应的影响.结果表明:1)当带匀速运动时,无论系统是线性还是非线性,运动带横向振动"频率"都随着带运动稳态速度增加而减小.2)随着材料粘性增加,系统耗散能力逐渐增强,动态响应逐渐减小.3)当带运动速度简谐波动时,系统动态响应随扰动速度增大而增大.扰动频率对带横向振动影响较大.     

微曲输流管道振动固有频率分析与仿真北大核心CSCD

首次建立了基于Timoshenko梁理论的微曲输流管道横向振动的动力学模型,并分析了流体流动影响下微曲管道横向自由振动的固有特征.采用广义Hamilton原理,导出了考虑流体影响的微曲管道横向振动的控制方程,通过Galerkin截断对控制方程离散化,再由广义本征值问题得到管道横向振动的固有频率,并研究了液体流速和弯曲幅度对管道横向固有振动特征的影响.发展了基于等效刚度和等效阻尼方法的考虑流体影响的微曲管道振动分析的有限元仿真计算方法,并通过有限元软件实现数值仿真,验证了Galerkin截断的分析结果以及所建立的Timoshenko微曲管道动力学模型的有效性.研究表明,流体的流速以及管道的弯曲幅度对管道横向振动固有频率均有显著影响.     

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.