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

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

NONLINEAR VIBRATIONS OF AXIALLY ACCELERATING VISCOELASTIC TIMOSHENKO BEAMS

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

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

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

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

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

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

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

具有分数导数本构关系的非线性粘弹性Timoshenko梁动力学行为分析

本文利用分数导数型本构关系建立了在有限变形情况下Timoshenko梁的控制方程并利用Galerkin方法进行简化。然后利用一种存储部分历史数据的分数积分的计算方法对梁的控制方程进行求解。考察了载荷参数和分数导数参数对梁振动的影响,并采用非线性动力学中的各种数值方法,如时程曲线、功率谱、相图、Poincare截面等,揭示了非线性粘弹性Timoshenko梁丰富的动力学行为。     

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

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

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

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

分数导数型本构关系描述粘弹性梁的振动分析

本文研究粘弹性梁在周期激励作用下的受迫振动问题。梁的材料满足Kelvin-Volgt分数导数型本构关系。基于动力学方程、本构关系和应变－位移关系建立了小变形粘弹性梁的振动方程。采用分离变量法分析粘弹性梁的自由振动,导出模态坐标满足的常微分－积分方程和模态函数满足的常微分方程,对于两端简支的截面梁给出了固有频率和模态函数。对于简谐激励作用下粘弹性梁的受迫振动,利用模态叠加得到了稳态响应。最后给出数值算例说明本文方法的应用。     

轴向运动SMA层合梁的主共振分析

研究了均布横向载荷作用下轴向运动SMA(形状记忆合金)层合梁的横向非线性振动。考虑轴向运动效应、轴力等因素的综合影响,利用力平衡条件、变形协调方程及SMA多项式函数的本构关系,建立了SMA层合梁在均布横向载荷作用下的动力学方程。针对两端简支边界条件,通过伽辽金积分得到了轴向运动SMA层合梁横向振动微分方程。应用平均法得到了横向载荷作用下系统主共振幅频响应方程,对理论结果进行了数值验证;分析了轴向运动速度、温度、激励参数对系统稳态响应的影响。结果表明:轴向速度、轴向载荷的变化只对系统共振频率产生影响;在外激励幅值较大时,温度增加和SMA层增厚对系统产生了相同的减振效果。     

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.     

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.     

Effect of rotary inertia on stability of axially accelerating viscoelastic Rayleigh beams

The dynamic stability of axially moving viscoelastic Rayleigh beams is presented. The governing equation and simple support boundary condition are derived with the extended Hamilton’s principle. The viscoelastic material of the beams is described as the Kelvin constitutive relationship involving the total time derivative. The axial tension is considered to vary longitudinally. The natural frequencies and solvability condition are obtained in the multi-scale process. It is of interest to investigate the summation parametric resonance and principal parametric resonance by using the Routh-Hurwitz criterion to obtain the stability condition. Numerical examples show the effects of viscosity coefficients, mean speed, beam stiffness, and rotary inertia factor on the summation parametric resonance and principle parametric resonance. The differential quadrature method (DQM) is used to validate the value of the stability boundary in the principle parametric resonance for the first two modes.     

Stability of axially accelerating viscoelastic beams: asymptotic perturbation analysis and differential quadrature validation

An asymptotic perturbation method is proposed to investigate stability of an axially accelerating viscoelastic beam. The material time derivative is used in the viscoelastic constitutive relation. The axial speed is characterized as a simple harmonic variation about the constant mean speed. The stability condition can be determined via the asymptotic perturbation method. The differential quadrature scheme is developed to solve numerically the equation of axially accelerating viscoelastic beams with simple supports. The stability boundaries are numerically located in the summation parametric resonance and the principal parametric resonance. Numerical examples show the effects of the beam viscoelasticity and the mean axial speed. The numerical calculations validate the analytical results in the principal parametric resonance.     

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.     

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.     

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.     

Stability of axially accelerating viscoelastic beams: multi-scale analysis with numerical confirmations

Stability is investigated for an axially accelerating viscoelastic beam. The material time derivative is used in the viscoelastic constitutive relation, not simply the partial time derivative. The method of multiple scales is applied directly to the governing equation without discretization. When the axial speed is characterized as a simple harmonic variation about the constant mean speed, the instability conditions are presented for axially accelerating viscoelastic beams constrained by simple supports with rotational springs in parametric resonance. The finite difference schemes are developed to solve numerically the equation of axially accelerating viscoelastic beams with fixed supports for the instability regions in the principal parametric resonance. The numerical calculations confirm the analytical results. Numerical examples show the effects of the constraint stiffness, the mean axial speed, and the viscoelasticity.