非齐次边界条件下轴向运动梁的非线性振动

轴向运动系统的横向非线性振动一直是国内外研究的热点课题之一.目前相关研究大都是针对齐次边界条件的.但是在工程实际中,非齐次边界条件更为常见,而针对非齐次边界条件的研究相对较少.为深入研究非齐次边界条件对轴向运动系统横向非线性振动的影响,本文以轴向变速运动黏弹性Euler梁为例,引入由黏弹性引起的非齐次边界条件,同时还引入由轴向加速度引起的径向变化张力,建立梁横向振动的积分-偏微分型运动方程,并导出了相应的非齐次边界条件.采用直接多尺度法分析了梁的次谐波参数共振.由可解性条件得到了梁的稳态响应,并根据Routh-Hurvitz判据确定了系统稳态响应的稳定性.通过数值例子讨论了黏弹性系数,轴向运动速度,轴向速度脉动幅值和非线性系数对幅频响应的影响,并详细对比分析了非齐次边界条件和齐次边界条件对幅频响应的影响.结果表明:随着黏弹性系数的增大,非齐次边界条件下的零解失稳区域和稳态响应幅值比齐次边界条件下的失稳区域和幅值大,非齐次边界条件对高阶次谐波参数共振的影响更加显著.最后,引入微分求积法来验证直接多尺度法的近似解结果.     

黏弹性轴向运动梁非线性受迫振动稳态幅频响应

本文研究了黏弹性轴向运动梁横向受迫振动稳态幅频响应问题.在控制方程的推导中,对黏弹性本构关系采用物质导数.把多尺度法直接应用于梁横向振动的非线性控制方程,利用可解性条件消除长期项,得到系统稳态的幅频响应曲线.运用Lyapunov一次近似理论分析幅频响应曲线的稳定性.通过算例研究了黏性系数,外部激励幅值以及非线性项系数对稳态幅频响应曲线及其稳定性的影响.运用数值方法对两端固定边界下黏弹性轴向运动梁的控制方程直接数值解,分析梁横向非线性振动的稳态幅频响应,通过数值算例验证直接多尺度法的结论.     

轴向变速粘弹性Rayleigh梁参数共振稳定性

运用近似解析方法和数值方法研究轴向变速运动黏弹性Rayleigh梁的次谐波共振和组合共振的稳定性区域。基于变分原理,考虑梁断面旋转惯性的影响,推导轴向速度有周期波动的微变形梁横向振动的数学模型;采用多尺度方法建立前两阶次谐波共振和组合共振范围内的参数振动的可解性条件;进而确定梁两端简支边界条件下,因共振而产生的失稳区域;通过微分求积方法求解表征细长Rayleigh梁横向振动的运动微分方程。数值算例分析了黏弹性系数和扭转系数对梁振动失稳区域的影响,将数值仿真结果与近似解析方法的结论进行比较。算例表明:近似解析解的精度较高,第一、第二阶主共振的最大误差分别为3.206%、4.213%。     

NONLINEAR VIBRATIONS OF AXIALLY ACCELERATING VISCOELASTIC TIMOSHENKO BEAMS

研究了轴向加速黏弹性Timoshenko梁的非线性参数振动。参数激励是由径向变化张力和轴向速度波动引起的。引入了取决于轴向加速度的径向变化张力,同时还考虑了有限支撑刚度对张力的影响。应用广义哈密尔顿原理建立了Timoshenko梁耦合平面运动的控制方程和相关的边界条件。黏弹性本构关系采用Kelvin模型并引入物质时间导数。耦合方程简化为具有随时间和空间变化系数的积分-偏微分型非线性方程。采用直接多尺度法分析了Timoshenko梁的组合参数共振。根据可解性条件得到了Timoshenko梁的稳态响应,并应用Routh-Hurvitz判据确定了稳态响应的稳定性。最后通过一系列数值例子描述了黏弹性系数、平均轴向速度、剪切变形系数、转动惯量系数、速度脉动幅值、有限支撑刚度参数以及非线性系数对稳态响应的影响。     

变速轴向运动三参数模型黏弹性梁的稳定性

本文研究了速度变化的轴向运动三参数模型黏弹性梁在主参数共振以及组合参数共振范围内的稳定性.轴向运动梁的黏弹性本构关系采用三参数模型并引入了物质时间导数.运用渐进摄动法,直接求解梁的控制微分方程并导出了当运动参数激励频率接近某一阶固有频率2倍或接近某两阶固有频率之和时主参数共振和组合参数共振的稳定性条件.在解谐参数和激励振幅平面上,可以找出由于共振而产生的失稳区域.数值结果给出了梁的刚度系数、黏弹性系数及轴向平均速度对失稳区域的影响.在发生组合共振和主共振时,随着刚度系数E1的变大,失稳区域变小;刚度系数E2的变大,失稳区域变大.随着黏弹性系数的变大,失稳区域变小.发生组合共振时,随着平均速度的变大,失稳区域变小;发生主共振时,随着平均速度的变大,失稳区域变大.     

时变张力作用下轴向运动梁的分岔与混沌

轴向运动结构的横向参激振动一直是非线性动力学领域的研究热点之一. 目前研究较多的是轴向速度摄动的动力学模型,参数激励由速度的简谐波动产生. 但在工程应用中,存在轴向张力波动的运动结构较为广泛,而针对轴向张力摄动的模型研究较少. 本文研究了时变张力作用下轴向变速运动黏弹性梁的分岔与混沌. 考虑随着时间周期性变化的轴向张力,计入线性黏性阻尼,采用Kelvin模型的黏弹性本构关系,给出了梁横向非线性 振动的积分--偏微分控制方程. 首先应用四阶Galerkin截断方法将控制方程离散化,然后采用四阶Runge-Kutta方法计算系统的数值解,进而确定其动力学行为. 基于梁中点的横向位移和速度的数值结果,仿真了梁沿平均轴速、张力摄动幅值、张力摄动频率以及黏弹性系数变化的倍周期分岔与混 沌运动,并且通过计算系统的最大李雅普诺夫指数来识别其混沌行为. 结果表明:较小的平均轴速有助于梁的周期运动,梁在临界速度附近容易发生倍周期分岔与混沌行为. 随着张力摄动幅值的增大,梁的振动幅值的混沌区间不断增大. 较小的黏弹性系数和张力摄动频率更容易使梁发生混沌运动. 最后,给出时程图、频谱图、相图以及Poincaré 映射图来确定梁的混沌运动.      

时变张力作用下轴向运动梁的分岔与混沌

轴向运动结构的横向参激振动一直是非线性动力学领域的研究热点之一.目前研究较多的是轴向速度摄动的动力学模型,参数激励由速度的简谐波动产生.但在工程应用中,存在轴向张力波动的运动结构较为广泛,而针对轴向张力摄动的模型研究较少.本文研究了时变张力作用下轴向变速运动黏弹性梁的分岔与混沌.考虑随着时间周期性变化的轴向张力,计入线性黏性阻尼,采用Kelvin模型的黏弹性本构关系,给出了梁横向非线性振动的积分—偏微分控制方程.首先应用四阶Galerkin截断方法将控制方程离散化,然后采用四阶Runge-Kutta方法计算系统的数值解,进而确定其动力学行为.基于梁中点的横向位移和速度的数值结果,仿真了梁沿平均轴速、张力摄动幅值、张力摄动频率以及黏弹性系数变化的倍周期分岔与混沌运动,并且通过计算系统的最大李雅普诺夫指数来识别其混沌行为.结果表明:较小的平均轴速有助于梁的周期运动,梁在临界速度附近容易发生倍周期分岔与混沌行为.随着张力摄动幅值的增大,梁的振动幅值的混沌区间不断增大.较小的黏弹性系数和张力摄动频率更容易使梁发生混沌运动.最后,给出时程图、频谱图、相图以及Poincaré映射图来确定梁的混沌运动.     

轴向变速运动黏弹性梁稳定性渐近分析和数值验证

研究了轴向变速运动黏弹性梁参数振动的稳定性.对黏弹性本构关系采用物质时间导数,轴向速度用关于恒定平均速度的简单谐波变化来描述.发展浙近摄动法确定稳定性条件.应用微分求积法数值求解简支边界条件下的轴向变速运动黏弹性梁方程,并进而确定次谐波参数共振的稳定性边界.数值结果显示了梁的黏性阻尼和轴向平均速度的影响并验证了次谐波共振的解析结果.     

数值仿真轴向运动黏弹性梁非线性参激振动

分别通过两种直接数值方法研究速度变化的经典边界条件下轴向运动黏弹性梁参数振动的稳定性。在控制方程的推导中,采用物质导数黏弹性本构关系和只对时间取偏导数的黏弹性本构关系;分别运用有限差分法和微分求积法对两种经典边界下轴向变速运动黏弹性梁的非线性控制方程求数值解,计算得到梁中点非线性参数振动的稳定稳态响应。数值结果表明,两种黏弹性本构关系对应的稳态响应存在明显差别,同时发现两种直接数值方法的仿真结果基本吻合,证明数值仿真具有较高精度。     

简谐荷载作用下弹性地基上不可伸长梁的2次超谐共振响应研究

基于已建立的弹性地基上不可伸长梁的非线性动力学模型,利用梁的量纲归一化运动方程和多尺度方法求得梁2次超谐共振的幅频响应方程和位移的二次近似解。进而,运用梁的幅频响应曲线对其超谐共振响应特性进行研究,同时分析了弹性地基模型、Winkler参数、外激励幅值、边界条件等对该共振响应的影响效应。结果表明:弹性地基模型中剪切参数的引入增大了梁2次超谐共振响应的幅值和多值区域;弹性地基Winkler参数的增加会抑制系统的共振响应,但同时会增加系统动力响应的软弹簧特性;在外激励幅值较小的情况下,系统共振响应未展现出明显的非线性特征;边界约束对弹性地基剪切参数作用于梁2次超谐共振响应的效应有显著影响,可在一定程度上改变系统响应幅值及多值区域。     

Steady-state response of axially moving viscoelastic beams with pulsating speed: comparison of two nonlinear models

Principal parametric resonance in transverse vibration is investigated for viscoelastic beams moving with axial pulsating speed. A nonlinear partial-differential equation governing the transverse vibration is derived from the dynamical, constitutive, and geometrical relations. Under certain assumption, the partial-differential reduces to an integro-partial-differential equation for transverse vibration of axially accelerating viscoelastic nonlinear beams. The method of multiple scales is applied to two equations to calculate the steady-state response. Closed form solutions for the amplitude of the vibration are derived from the solvability condition of eliminating secular terms. The stability of straight equilibrium and nontrivial steady-state response are analyzed by use of the Lyapunov linearized stability theory. Numerical examples are presented to highlight the effects of speed pulsation, viscoelascity, and nonlinearity and to compare results obtained from two equations.     

Periodic responses and chaotic behaviors of an axially accelerating viscoelastic Timoshenko beam

This paper investigates the steady-state periodic response and the chaos and bifurcation of an axially accelerating viscoelastic Timoshenko beam. For the first time, the nonlinear dynamic behaviors in the transverse parametric vibration of an axially moving Timoshenko beam are studied. The axial speed of the system is assumed as a harmonic variation over a constant mean speed. The transverse motion of the beam is governed by nonlinear integro-partial-differential equations, including the finite axial support rigidity and the longitudinally varying tension due to the axial acceleration. The Galerkin truncation is applied to discretize the governing equations into a set of nonlinear ordinary differential equations. Based on the solutions obtained by the fourth-order Runge–Kutta algorithm, the stable steady-state periodic response is examined. Besides, the bifurcation diagrams of different bifurcation parameters are presented in the subcritical and supercritical regime. Furthermore, the nonlinear dynamical behaviors are identified in the forms of time histories, phase portraits, Poincaré maps, amplitude spectra, and sensitivity to initial conditions. Moreover, numerical examples reveal the effects of various terms Galerkin truncation on the amplitude–frequency responses, as well as bifurcation diagrams.     

Asymptotic analysis on weakly forced vibration of axially moving viscoelastic beam constituted by standard linear solid model

The weakly forced vibration of an axially moving viscoelastic beam is investigated.The viscoelastic material of the beam is constituted by the standard linear solid model with the material time derivative involved.The nonlinear equations governing the transverse vibration are derived from the dynamical,constitutive,and geometrical relations.The method of multiple scales is used to determine the steady-state response.The modulation equation is derived from the solvability condition of eliminating secular terms.Closed-form expressions of the amplitude and existence condition of nontrivial steady-state response are derived from the modulation equation.The stability of nontrivial steady-state response is examined via the Routh-Hurwitz criterion.     

Internal resonance of an axially transporting beam with a two-frequency parametric excitation

This paper investigates the transverse 3:1 internal resonance of an axially transporting nonlinear viscoelastic Euler-Bernoulli beam with a two-frequency parametric excitation caused by a speed perturbation. The Kelvin-Voigt model is introduced to describe the viscoelastic characteristics of the axially transporting beam. The governing equation and the associated boundary conditions are obtained by Newton’s second law. The method of multiple scales is utilized to obtain the steady-state responses. The Routh-Hurwitz criterion is used to determine the stabilities and bifurcations of the steady-state responses. The effects of the material viscoelastic coefficient on the dynamics of the transporting beam are studied in detail by a series of numerical demonstrations. Interesting phenomena of the steady-state responses are revealed in the 3:1 internal resonance and two-frequency parametric excitation. The approximate analytical method is validated via a differential quadrature method.     

Using Fourier differential quadrature method to analyze transverse nonlinear vibrations of an axially accelerating viscoelastic beam

In this paper, a Fourier expansion-based differential quadrature (FDQ) method is developed to analyze numerically the transverse nonlinear vibrations of an axially accelerating viscoelastic beam. The partial differential nonlinear governing equation is discretized in space region and in time domain using FDQ and Runge–Kutta–Fehlberg methods, respectively. The accuracy of the proposed method is represented by two numerical examples. The nonlinear dynamical behaviors, such as the bifurcations and chaotic motions of the axially accelerating viscoelastic beam, are investigated using the bifurcation diagrams, Lyapunov exponents, Poincare maps, and three-dimensional phase portraits. The bifurcation diagrams for the in-plane responses to the mean axial velocity, the amplitude of velocity fluctuation, and the frequency of velocity fluctuation are, respectively, presented when other parameters are fixed. The Lyapunov exponents are calculated to further identify the existence of the periodic and chaotic motions in the transverse nonlinear vibrations of the axially accelerating viscoelastic beam. The conclusion is drawn from numerical simulation results that the FDQ method is a simple and efficient method for the analysis of the nonlinear dynamics of the axially accelerating viscoelastic beam.     

NONLINEAR DYNAMICS OF AXIALLY ACCELERATING VISCOELASTIC BEAMS BASED ON DIFFERENTIAL QUADRATURE

This paper investigates nonlinear dynamical behaviors in transverse motion of an axially accelerating viscoelastic beam via the differential quadrature method. The governing equation, a nonlinear partial-differential equation, is derived from the viscoelastic constitution relation using the material derivative. The differential quadrature scheme is developed to solve numerically the governing equation. Based on the numerical solutions, the nonlinear dynamical behaviors are identified by use of the Poincare map and the phase portrait. The bifurcation diagrams are presented in the case that the mean axial speed and the amplitude of the speed fluctuation are respectively varied while other parameters are fixed. The Lyapunov exponent and the initial value sensitivity of the different points of the beam, calculated from the time series based on the numerical solutions, are used to indicate periodic motions or chaotic motions occurring in the transverse motion of the axially accelerating viscoelastic beam.     

Primary resonance in forced vibrations of in-plane translating viscoelastic plates with 3:1 internal resonance

Nonlinear forced vibrations of in-plane translating viscoelastic plates subjected to plane stresses are analytically and numerically investigated on the steady-state responses in external and internal resonances. A nonlinear partial-differential equation with the associated boundary conditions governing the transverse motion is derived from the generalized Hamilton principle and the Kelvin relation. The method of multiple scales is directly applied to establish the solvability conditions in the primary resonance and the 3:1 internal resonance. The steady-state responses are predicted in two patterns: single-mode and two-mode solutions. The Routh?CHurvitz criterion is used to determine the stabilities of the steady-state responses. The effects of the in-plane translating speed, the viscosity coefficient, and the excitation amplitude on the steady-state responses are examined. The differential quadrature scheme is developed to solve the nonlinear governing equations numerically. The numerical calculations confirm the approximate analytical results regarding the single-mode solutions of the steady-state responses.