en
×

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。
通讯作者:

严巧赟,E-mail:qyyan@shu.edu.cn

中图分类号:O32

文献标识码:A

文章编号:1672-6553-2023-21(9)-041-009

DOI:10.6052/1672-6553-2023-061

参考文献 1
GHAYESH M H,KAZEMIRAD S,AMABILI M.Coupled longitudinal-transverse dynamics of an axially moving beam with an internal resonance [J].Mechanism and Machine Theory,2012,52:18-34.
参考文献 2
HU W P,YE J,DENG Z C.Internal resonance of a flexible beam in a spatial tethered system [J].Journal of Sound and Vibration,2020,475:115286.
参考文献 3
DING H,HUANG L L,MAO X Y,et al.Primary resonance of traveling viscoelastic beam under internal resonance [J].Applied Mathematics and Mechanics-English Edition,2017,38(1):1-14.
参考文献 4
WANG Y R,WEI Y H.Internal resonance analysis of a fluid-conveying tube resting on a nonlinear elastic foundation [J].European Physical Journal Plus,2020,135(4):364.
参考文献 5
WANG G X,DING H,CHEN L Q.Dynamic effect of internal resonance caused by gravity on the nonlinear vibration of vertical cantilever beams [J].Journal of Sound and Vibration,2020,474:115265.
参考文献 6
GUILLOT V,GIVOIS A,COLIN M,et al.Theoretical and experimental investigation of a 1:3 internal resonance in a beam with piezoelectric patches [J].Journal of Vibration and Control,2020,26(13-14):1119-1132.
参考文献 7
DING H,HUANG L L,DOWELL E,et al.Stress distribution and fatigue life of nonlinear vibration of an axially moving beam [J].Science China Technological Sciences,2019,62(7):1123-1133.
参考文献 8
YOUNESIAN D,SADRI M,ESMAILZADEH E.Primary and secondary resonance analyses of clamped-clamped micro-beams [J].Nonlinear Dynamics,2014,76(4):1867-1884.
参考文献 9
顾伟,张博,丁虎,等.2:1内共振条件下变转速预变形叶片的非线性动力学响应 [J].力学学报,2020,52(7):1131-1142.GU W,ZHANG B,DING HU,et al.Nonlinear dynamic response of pre-deformed blade with variable rotational speed under 2:1 internal resonance [J].Chinese Journal of Theoretical and Applied Mechanics,2020,52(7):1131-1142.(in Chinese)
参考文献 10
毕勤胜,陈予恕.双摆内共振分岔分析 [J].应用数学和力学,2000,21(3):226-234.BI Q S,CHEN Y S.Bifurcation analysis of a double pendulum with internal resonance [J].Applied Mathematics and Mechanics,2000,21(3):226-234.(in Chinese)
参考文献 11
魏明海,肖仪清.内外共振联合作用下索-梁组合结构非线性振动分析 [J].振动与冲击,2012,31(7):79-84.WEI H M,XIAO Y Q.Nonlinear vibration analysis for a cable-beam coupled system under simultaneous internal and external resonances [J].Vibration and Shock,2012,31(7):79-84.(in Chinese)
参考文献 12
LI S H,YANG S P,CHEN L Q,et al.Effects of parameters on dynamic responses for a heavy vehicle-pavement-foundation coupled system [J].International Journal of Heavy Vehicle Systems,2012,19(2):207-224.
参考文献 13
DING H,CHEN L Q.Approximate and numerical analysis of nonlinear forced vibration of axially moving viscoelastic beams [J].Acta Mechanica Sinica,2011,27(3):426-437.
参考文献 14
DONG Z J,MA X Y.Analytical solutions of asphalt pavement responses under moving loads with arbitrary non-uniform tire contact pressure and irregular tire imprint [J].Road Materials and Pavement Design,2018,19(8):1887-1903.
参考文献 15
杨予,滕念管,黄醒春,等.承受移动均布质量的简支梁振动反应分析 [J].振动与冲击,2005,24(3):19-22+6.YANG Y,TENG N G,HUANG X C,et al.Vibration analysis of a simply supported beam traversed by uniform distributed moving mass [J].Vibration and Shock,2005,24(3):19-22+6.(in Chinese)
参考文献 16
刘涛,周洋忻,胡伟鹏.轴向运动功能梯度梁横向振动问题的保结构分析 [J].动力学与控制学报,2022,20(6):101-105.LIU T,ZHOU Y X,HU W P.Structure-preserving analysis on transverse vibration of functionally graded beam with an axial velocity [J].Journal of Dynamics and Control,2022,20(6):101-105.(in Chinese)
参考文献 17
MARTINEZ-CASTRO A E,MUSEROS P,CASTILLO-LINARES A.Semi-analytic solution in the time domain for non-uniform multi-span Bernoulli-Euler beams traversed by moving loads [J].Journal of Sound and Vibration,2006,294(1-2):278-297.
参考文献 18
张弛,毛晓晔,丁虎,等.受轴向激励弹性支承梁的稳定性分析 [J].动力学与控制学报,2022,20(3):66-76.ZHANG C,MAO X Y,DING H,et al.Stability analysis of axially excited beam with elastic boundary [J].Journal Dynamics and Control,2022,20(3):66-76.(in Chinese)
参考文献 19
YOUNESIAN D,MARJANI S R,ESMAILZADEH E.Nonlinear vibration analysis of harmonically excited cracked beams on viscoelastic foundations [J].Nonlinear Dynamics,2013,71(1-2):109-120.
参考文献 20
DING H,YANG Y,CHEN L Q,et al.Vibration of vehicle-pavement coupled system based on a Timoshenko beam on a nonlinear foundation [J].Journal of Sound and Vibration,2014,333(24):6623-6636.
参考文献 21
ANSARI M,ESMAILZADEH E,YOUNESIAN D.Frequency analysis of finite beams on nonlinear Kelvin-Voight foundation under moving loads [J].Journal of Sound and Vibration,2011,330(7):1455-1471.
参考文献 22
SAPOUNTZAKIS E J,KAMPITSIS A E.Nonlinear response of shear deformable beams on tensionless nonlinear viscoelastic foundation under moving loads [J].Journal of Sound and Vibration,2011,330(22):5410-5426.
目录contents

    摘要

    研究了有限长弹性基础上梁在移动载荷作用下的内共振响应.建立了移动集中力激励的非线性粘弹性基础支承的有限长Euler-Bernoulli梁模型,并对非线性偏微分方程进行离散,在第三阶固有频率与第一阶固有频率成三倍关系时,采用多尺度方法导出了3:1内共振的可解性条件,研究了有无移动载荷时基础阻尼和非线性刚度对梁内共振条件下自由振动响应和受迫振动响应的影响规律.在此基础上,应用Lyapunov第一方法确定了系统的稳定性条件.

    Abstract

    The internal external resonance response of a beam on a finite length elastic foundation under moving loads is studied. A finite length Euler-Bernoulli beam model supported on a nonlinear viscoelastic foundation excited by a moving concentrated force is established, and the nonlinear partial differential equation is discretized. When the third-order natural frequency and the first-order natural frequency are three times related, the solvability condition for 3:1 internal resonance is derived by using a multi-scale method. The effects of foundation damping and nonlinear stiffness on the free vibration response and forced vibration response of a beam under internal resonance conditions with or without moving loads are studied. On this basis, the stability conditions of the system are determined by using the first Lyapunov method.

  • 引言

  • 移动荷载下基础梁能够作为很多工程结构的力学模型,例如,受到来自高速、重载车辆的作用的非线性基础路面结构.研究移动载荷作用下基础梁的振动分析是结构动力学中的一个重要课题,也是机械、桥梁和铁路工程等不同学科的研究热点.

  • Ghayesh等[1]以3:1内共振为例,采用数值的方法研究了梁轴向运动引起的纵向和横向耦合非线性振动的动力学特性,根据不同的参数分析其表现出的动态特性.Hu等[2]利用多尺度法,得到了空间轨道系留系统中两弹簧悬挂的空间柔性梁的内共振条件,研究了内共振对姿态稳定性的影响以及构件间能量传递趋势的影响.Ding等[3]在3:1内共振条件下,研究了粘弹性移动梁受迫振动的稳态周期响应,采用数值算例分析了粘弹性行为对稳态周期响应的影响.Wang等[4]采用输送流体线性梁和非线性弹簧来模拟在弹性基础上放置的流体输送管的振动,用多尺度法求出了各模态在稳态下的频率响应,并用各模态的振幅来检验内共振,通过改变流体的流速来分析系统的稳定性.Wang和Ding[5]研究了悬臂梁的非线性自由振动、1/3超谐共振以及重力引起的3:1内共振,采用时间多尺度方法,获得其对应振动的模态响应,运用谐波平衡法求解超谐波共振中的激励分量和模型的挠度,研究了重力引起的内部共振对垂直悬臂梁非线性振动的动态影响.Guillot等[6]对压电片梁中1:3内共振现象,以实验和理论的方法进行了研究.Ding等[7]采用直接多尺度方法研究了非线性振动对轴向运动梁的应力分布和疲劳寿命的影响,揭示了内共振对V型带稳态响应和疲劳寿命的影响.Younesian等[8]采用Galerkin方法建立了微梁的无量纲运动控制方程,结合多尺度方法求解非线性方程组,并确定了主共振和次共振条件,研究了不同谐振情况下的谐波响应.顾伟等[9]研究了预变形叶片在变速条件下的非线性动力学行为,详细研究了温度梯度、阻尼、转速扰动幅度等系统参数对叶片动态响应的影响,在2:1内共振条件下,研究了立方项对方程的影响.毕勤胜和陈予恕[10]分析了复摆自治系统在1:1内共振时的混沌和分岔特性.魏明海等[11]研究了内共振和外共振联合作用下的索-梁组合结构非线性振动问题.在众多关于飞线振动的研究成果中,与道路相关的非线性系统内-外联合共振的研究鲜有报道.

  • 目前对于连续体的共振问题多是将连续体以梁或板的形式进行研究[12-14].杨予等[15]将桥简化为简支梁,研究其在均布移动载荷作用下的振动响应,指出梁的动态响应与其自振频率、载荷的行驶频率和荷载与梁质量的振动相关.刘涛等[16]分析了轴向运动速度和材料的非均匀性对轴向运动功能梯度梁发生共振时的影响.Martínez-Castro[17]以时变模态方程得到了移动荷载作用下非均匀多跨Euler-Bernoulli梁的时域半解析解,并证实了该方法具有很高的精度和鲁棒性.张弛等[18]针对轴向移动梁结构,讨论了边界支撑参数对振动的影响.Younesian等[19]采用Galerkin方法结合多尺度方法求解非线性运动控制方程,研究了由非线性粘弹性基础支承的裂纹梁的频率响应,分析了不同参数对频率响应解的影响.Ding等[20]着眼于车辆与路面的耦合非线性振动问题,将路面建模为非线性基础上的Timoshenko梁,研究了汽车-路面耦合系统的动力学响应.

  • 本文建立了在移动载荷作用下的非线性粘弹性基础上的有限长Euler-Bernoulli梁,且系统的第三阶固有频率与第一阶固有频率成三倍关系,用多尺度法得到内共振条件下非线性系统的自由振动和强迫振动响应.通过参数分析,研究了不同参数对系统动态响应的影响规律.

  • 1 非线性粘弹性基础梁模型

  • 考虑在移动载荷作用下基于非线性粘弹性Kelvin基础的有限长Euler-Bernoulli梁模型,如图1所示.基础梁密度为ρ,梁的宽度和厚度分别为bh,惯性矩为I,路面长度为L,弹性模量是E.非线性Kelvin基础模型可以如下表示[2122]

  • 图1 非线性粘弹性基础上的有限Euler-Bernoulli梁示意图

  • Fig.1 Finite Euler-Bernoulli beam on nonlinear viscoelastic foundation

  • P=k1Y+k3Y3+μYT
    (1)
  • 其中,P是在梁上任意点发生位移的力,k1k3分别是基础的线性刚度和非线性刚度系数,μ是基础的阻尼系数.根据Euler-Bernoulli梁理论,Kelvin基础系统的控制微分方程可以表示为

  • EI4wX4+ρA2wT2+k1w+k3w3+μwT=

  • Fδ(X-VT)
    (2)
  • 假定梁两端均为简支,因此,边界条件如下:

  • w (X, T) X=0=w (X, T) X=L=0,

  • EIψX(X,T)X=0=EIψX(X,T)X=L=0
    (3)
  • 根据边界条件,梁的位移函数的表达式如下:

  • w(X,T)=i=1 qi(T)φi(X)
    (4)
  • 其中,φiX)表示第i阶模态函数,qiT)表示关于时间T的第i阶广义位移.梁的模态函数表示为

  • φi(X)=siniπXL
    (5)
  • 将式(4)代入式(2)得,

  • i=1n q¨i (T) +2εμ^q˙i (T) +ωi2qi (T) siniπXL+

  • k^3i=1n qi(T)siniπXL3=εLfδ(X-VT)
    (6)
  • 其中,μρA=2εμ^k1ρA+i4π4EIρAL4=ωi2k3ρA=k^3FρA=L2εfε为小参数(0<ε1).

  • 根据三角函数的正交性,上式等号两边同时乘以φmX=sinmπXL,然后在区间长度[0,L]上进行积分,取Galerkin截断前三阶来求解,可得:

  • q¨1 (T) +ω12q1 (T) =-2εμ^q˙1 (T) +ε-α1q13-

  • α2q12q3-α3q1q22-α4q1q32-α5q22q3+

  • εfsinΩ1T
    (7)
  • q¨2 (T) +ω22q2 (T) =-2εμ^q˙2 (T) +

  • ε-α6q12q2-α7q1q2q3-α8q23-

  • α9q2q32+εfsinΩ2T
    (8)
  • q¨3 (T) +ω32q3 (T) =-2εμ^q˙3 (T) +

  • ε-α10q13-α11q12q3-α12q1q22-α13q22q3-

  • α14k^3q33+εfsinΩ3T
    (9)
  • 其中,Ωi=iπV/Lα2=-3/4k^3α1=α5=α8=α12=α14=3/4k^3α3=α4=α6=α7=α9=α11=α13=3/2k^3α10=-1/4k^3.

  • 运用多尺度法,将方程的解设为

  • qi=qi0T0,T1+εqi1T0,T1+
    (10)
  • 其中,Tn=εntn=1,23,将式(10)代入式(7)~(9)可得:

  • ε0阶:

  • D02q10+ω12q10=0
    (11)
  • D02q20+ω22q20=0
    (12)
  • D02q30+ω32q30=0
    (13)
  • ε1阶:

  • D02q11+ω12q11=-2D0D1q10-2μ^D0q10-

  • α1q103-α2q102q30-α3q10q202-α4q10q302-

  • α5q202q30+fsinΩ1T
    (14)
  • D02q21+ω22q21=-2D0D1q20-2μ^D0q20-

  • α6q102q20-α7q10q20q30-α8q203-α9q20q302+

  • 2fsinΩ2T
    (15)
  • D02q31+ω32q31=-2D0D1q30-2μ^D0q30-

  • α10k^3q103-α11q102q30-α12q10q202-

  • α13q202q30-α14q303+2fsinΩ3T
    (16)
  • 由式(11)~(13)可以得到下列解的形式.

  • q10=A1T1ejω1T0+cc
    (17)
  • q20=A2T1ejω2T0+cc
    (18)
  • q30=A3T1ejω2T0+cc
    (19)
  • 其中,为上式中所对应的共轭项.

  • 将式(17)~(19)代入到式(14)~(16)可得:

  • D02q11+ω12q11=-2jω1A1'ejω1T0+2jω1A-1'e-jω1T0-

  • 2jμ^ω1A1ejω1T0+2jμ^ω1A-1e-jω1T0+

  • -α1A13e3jω1T0-3α1A12A-1ejω1T0-

  • 3α1A1A-12e-jω1T0-α1A-13e-3jω1T0-

  • α2A12A3ej2ω1+ω3T0-2α2A1A-1A3ejω3T0-

  • α2A-12A3e-j2ω1-ω3T0-α2A12A-3ej2ω1-ω3T0-

  • 2α2A1A-1A-3e-jω3T0-α2A-12A-3e-j2ω1+ω3T0-

  • α3A1A22ejω1+2ω2T0-2α3A1A2A-2ejω1T0-

  • α3A1A-12ejω1-2ω2T0-α3A-1A22e-jω1-2ω2T0-

  • 2α3A-1A2A-2e-jω1T0-α3A-1A-22e-jω1+2ω2T0-

  • α4A1A32ejω1+2ω3T0-2α4A1A3A-3ejω1T0-

  • α4A1A-32ejω1-2ω3T0-α4A-1A32e-jω1-2ω3T0-

  • 2α4A-1A3A-3e-jω1T0-α4A-1A-32e-jω1+2ω3T0-

  • α5A22A3ej2ω2+ω3T0-2α5A2A-2A3ejω3T0-

  • α5A-22A3e-j2ω2-ω3T0-α5A22A-3ej2ω2-ω3T0-

  • 2α5A2A-2A-3e-jω3T0-α5A-22A-3e-j2ω2+ω3T0+

  • 1/2fejΩ1T-π2+e-jΩ1T-π2
    (20)
  • D02q21+ω22q21=-2jω2A2'ejω2T0+2jω2A-2'e-jω2T0-

  • 2jμ^ω2A2ejω2T0+2jμ^ω2A-2e-jω2T0+

  • -α6k^3A12A2ej2ω1+ω2T0-2α6A-1A1A2ejω2T0-

  • α6A-12A2e-j2ω1-ω2T0-α6A12A-2ej2ω1-ω2T0-

  • 2α6A-1A1A-2e-jω2T0-α6A-12A-2e-j2ω1+ω3T0

  • α7A1A2A3ejω1+ω2+ω3T0-α7A-1A2A3e-jω1-ω2-ω3T0-

  • α7A1A-12ejω1-2ω2T0-α3A-1A22e-jω1-2ω2T0-

  • α7A1A2A-3ejω1+ω2-ω3T0-α7A-1A2A-3e-jω1-ω2+ω3T0-

  • α7A1A-2A-3ejω1-ω2-ω3T0-α7A-1A-2A-3e-jω1+ω2+ω3T0-

  • α8A23e3jω2T0-3α8A22A-2ejω2T0-

  • 3α8A2A-22e-jω2T0-α8A-23e-3jω2T0-

  • α9A2A32ejω2+2ω3T0-2α9A2A-3A3A3ejω22T0-

  • α9A2A-32ejω2-2ω3T0-α9A-2A32e-jω2-2ω3T0-

  • 2α9A-2A-3A3e-jω2T0-α9A-2A-32e-jω2+2ω3T0+

  • fejΩ2T-π2+e-jΩ2T-π2
    (21)
  • D02q31+ω32q31=-2jω3A3'ejω3T0+2jω3A-3'e-jω3T0

  • 2jμ^ω3A3ejω3T0+2jμ^ω3A-3e-jω3T0+-α10A13e3jω1T0-3α10A12A-1ejω1T0-3α10A1A-12e-jω1T0-α10A-13e-3jω1T0-α11A12A3ej2ω1+ω3T0-2α11A1A-1A3ejω3T0-α11A-12A3e-j2ω1-ω3T0-α11A12A-3ej2ω1-ω3T0-2α11A1A-1A-3e-jω3T0-α11A-12A-3e-j2ω1+ω3T0-α12A1A22ejω1+2ω2T0-2α12A1A2A-2ejω1T0-α12A1A-22ejω1-2ω2T0-α12A-1A22e-jω1-2ω2T0-2α12A-1A2A-2e-jω11T0-α12A-1A-22e-jω1+2ω2T0-α13A22A3ej2ω2+ω3T0-2α13A2A-2A3ejω3T0-α13A-22A3e-j2ω2-ω3T0-α13A22A-3ej2ω2-ω3T0-2α13A2A-2A-3e-jω3T0-α13A-22A-3e-j2ω2+ω3T0-

  • α14A33e3jω3T0-3α14A32A-3ejω3T0-3α14A3A-32e-jω3T0-α14A-33e-3jω3T0-

  • fejΩ3T-π2+e-jΩ3T-π2
    (22)
  • 2 3:1内共振条件下的稳定性分析

  • 为研究3:1内共振下基础梁的自由振动和强迫振动响应,将系统第3阶模态的共振频率与第1阶模态共振频率的3倍设为相近,即,进而可以确定基础的线性刚度.其中,以及下文提到的均为引入的调谐因子.

  • 2.1 3:1内共振下的自由振动

  • 分析自由振动,即考虑在无移动载荷作用下的振动响应.消除久期项的条件为:

  • -2jω1A1'-2jμ^ω1A1-3α1A12A-1-α2A-12A3ejσ1T1-2α3A1A2A-2-

  • 2α4A1A3A-3=0
    (23)
  • -2jω2A2'-2jμ^ω2A2-2α6A-1A1A2-

  • 3α8A22A-2-2α9A2A-3A3=0
    (24)
  • -2jω3A3'-2jμ^ω3A3-α10A13e-jσ1T1-

  • 2α11A1A-1A3-2α13A2A-2A3-3α14A32A-3=0
    (25)
  • A1=1/2a1ejβ1A2=1/2a2ejβ2A3=1/2a3ejβ3,并代入到上式,再分离实部和虚部,可得:

  • ω1β1'a1-3/8α1a13-1/8α2a12a3cosγ-

  • 1/4α4a1a22-1/4α4a1a32=0
    (26)
  • -ω1a1'-a1μ^ω1-1/8α2a12a3sinγ=0
    (27)
  • ω2β2'a2-1/4α6a12a2-1/4α9a2a32-

  • 3/8α8a23=0
    (28)
  • -ω2a2'-μ^ω2a2=0
    (29)
  • ω3β3'a3-1/8α10a13cosγ-1/4α11a12a3-

  • 1/4α13a22a3-3/8α14a33=0
    (30)
  • -ω3a3'-μ^ω3a3+1/8α10a13sinγ=0
    (31)
  • 其中,γ=σ1T1+β3-3β1.

  • 式(29)在稳态条件下的响应为

  • a2=0
    (32)
  • 由式(27)、式(31)和式(32)可得:

  • a1a1'+ω3α2ω1α10a3a3'=-μ^a12-ω3α2ω1α10μ^a3
    (33)
  • f13T1=a12+ω3α2ω1α10a32,将其代入式(33),可得:

  • df13T1dT1=-2μ^f13T1
    (34)
  • 由式(34)可得:

  • f13T1=Ce-2μ^T1
    (35)
  • 其中C为大于0的常数.

  • 从式(35)可以看出,当T1→+∞时,f13(+∞)=0,从而理论分析上可以说明系统是随着时间逐渐衰减为0.同样,在时域下变量a1a3的数值解可由式(26)~式(31)得到,如图2~图7所示.图2通过数值仿真给出了系统在自由振动下的衰减过程.

  • 阻尼和非线性刚度对系统自由振动的影响如图3所示.从图3(a)和3(c)可以看出随着阻尼增加,系统的衰减速度就会越快;图3(b)和3(d)可以看出,非线性刚度会增加系统的振动频率,但对其衰减速度基本无影响.

  • 图2 自由振动下幅值的衰减过程

  • Fig.2 Attenuation process of amplitude under free vibration

  • 图3 自由振动下参数对系统响应的影响

  • Fig.3 Influence of parameters on system response under free vibration

  • 2.2 3:1内共振下的受迫振动分析

  • 在3:1内共振(即ω3≈3ω1时)下分别分析Ω1ω1Ω2ω2Ω3ω3时系统发生强迫共振条件.

  • (1)情况I:Ω1ω1

  • 假定

  • Ω1=ω1+εσ2
    (36)
  • 消除久期项,并分离实部和虚部,可以得到:

  • ω1β1'a1-3/8α1a13-1/4/α4a1a22-1/4α4a1a32-

  • 1/8α2a12a3cosγ1+f/2cosγ2=0
    (37)
  • -ω1a1'-μ^ω1a1-1/8α2a12a3sinγ1+

  • f/2sinγ2=0
    (38)
  • ω2β2'a2-1/4α6a12a2-1/4α9a2a32-

  • 3/8α8a23=0
    (39)
  • -ω2a2'-μ^ω2a2=0
    (40)
  • ω3β3'a3-1/8α10a13cosγ1-1/4α11a12a3-

  • 1/4α13a22a3-3/8α14a33=0
    (41)
  • -ω3a3'-μ^ω3a3+1/8α10a13sinγ1=0
    (42)
  • 其中,γ1=σ1T1+β3-3β1γ2=σ2T1-β1-π/2.

  • (2)情况II:Ω2ω2

  • 假定

  • Ω2=ω2+εσ2
    (43)
  • 消除久期项,并分离实部和虚部,可以得到:

  • ω1β1'a1-3/8α1a13-1/4α4a1a22-1/4α4a1a32-

  • 1/8α2a12a3cosγ1=0
    (44)
  • -μ^ω1a1-1/8α2a12a3sinγ1=0
    (45)
  • ω2σ2a2-1/4/α6a12a2-1/4α9a2a32-

  • 3/8α8a23+f/2cosγ2=0
    (46)
  • -μ^ω2a2+f/2sinγ2=0
    (47)
  • ω33β1'-σ1a3-1/8α10a13cosγ1-1/4α11a12a3-

  • 1/4α13a22a3-3/8α14a33=0
    (48)
  • -μ^ω3a3+1/8α10a13sinγ1=0
    (49)
  • 其中,γ1=σ1T1+β3-3β1γ2=σ2T1-β2-π/2.

  • (3)情况III:Ω3ω3

  • 假定

  • Ω3=ω3+εσ3
    (50)
  • 消除久期项,并分离方程的实部和虚部可以得到式(51)~式(56).

  • ω1β1'a1-3/8α1a13-1/8α2a12a3cosγ-

  • 1/4α4a1a22-1/4α4a1a32=0
    (51)
  • -ω1a1'-a1μ^ω1-1/8α2a12a3sinγ=0
    (52)
  • ω2β2'a2-1/4α6a12a2-1/4α9a2a32-

  • 3/8α8a23=0
    (53)
  • -ω2a2'-μ^ω2a2=0
    (54)
  • ω33β1'-σ1a3-1/8α10a13cosγ1-

  • 1/4α11a12a3-1/4α13a22a3-

  • 3/8α14a33+1/2fcosγ2=0
    (55)
  • -μ^ω3a3+1/8α10a13sinγ1+1/2fsinγ2=0
    (56)
  • 其中,γ1=σ1T1+β3-3β1γ2=σ2T1-β3-π/2.

  • (4)稳定性分析

  • 有限长梁和非线性基础的物理几何参数如表1所示.以情况I为例分析系统的稳定性,观察方程组式(37)~式(42)可以看出稳态下a2=a2=0,则方程组可以写为

  • ω1β1'a1-3/8α1a13-1/4α4a1a32-

  • 1/8α2a12a3cosγ1+f/2cosγ2=0
    (57)
  • -ω1a1'-μ^ω1a1+f/2sinγ2=0
    (58)
  • ω3β3'a3-1/8α10a13cosγ1-1/4α11a12a3-

  • 3/8α14a33=0
    (59)
  • -ω3a3'-μ^ω3a3+1/8α10a13sinγ1=0
    (60)
  • 又因为

  • β3'=γ1'-3γ2'+3σ2-σ1
    (61)
  • β1'=σ2-γ2'
    (62)
  • 表1 基础梁和移动载荷的物理几何参数[1]

  • Table1 Physical geometric parameters of foundation beams and moving load

  • 将式(61)和式(62)代入式(57)~式(60)可得,

  • a1'=-μ^a1-α28ω1a12a3sinγ1+f2ω1sinγ2
    (63)
  • a3'=-μ^a3+α108ω3k^3a13sinγ1
    (64)
  • γ1'=σ1+α114ω3-9α18ω1k^3a12+3α148ω3-

  • 3α44ω1a32+α108ω3a3a13-3α28ω1a1a3cosγ1+

  • 3f2a1ω1cosγ2
    (65)
  • γ2'=σ2-3α18ω1a12-α44ω1a32-α28ω1a1a3cosγ1+

  • f2a1ω1cosγ2
    (66)
  • a1a3γ1γ2T为状态向量,通过线性化式(63)~式(66)可以得到Jacobian矩阵

  • (67)
  • 其中

  • J11=-μ^-α24ω1a1a3sinγ1,

  • J12=-α28ω1a12sinγ1,

  • J13=-α28ω1a12a3cosγ1

  • J14=f2ω1cosγ2, J21=3α108ω3k^3a1sinγ1,

  • J22=-μ^, J23=α108ω3k^3a13cosγ1, J24=0,

  • J31=α112ω3-9α14ω1k^3a1+

  • 2α108ω3a3a12-3α28ω1a3cosγ1-3f2a12ω1cosγ2,

  • J32=3α144ω3-3α42ω1a3+

  • -α108ω3a32a13-3α28ω1a1cosγ1,

  • J33=-α108ω3a3a13-3α28ω1a1a3sinγ1,

  • J34=-3f2a1ω1sinγ2

  • J41=-3α14ω1a1-α28ω1a3cosγ1-f2a12ω1cosγ2

  • J42=-α42ω1a3-α28ω1a1cosγ1,

  • J43=α28ω1a1a3sinγ1, J44=-f2a1ω1sinγ2

  • 根据Lyapunov第一方法,可以判定:若矩阵J的特征值的实部全为负则系统稳定,反之系统不稳定.本文根据以上所述的方法,分别给出了在情况I、情况II和情况III下当k3=8×109.9N/m4时的稳定性分析图,如图4~图6所示.

  • 图4 k3=8×109.9N/m4时,幅频响应及其稳定性(情况Ⅰ)

  • Fig.4 Amplitude-frequency response and its stability at k3=8×109.9N/m4 (Case Ⅰ)

  • 图4(c)画出了情况Ⅰ下特征值实部的变化情况,可以看出在区间(σ2Bσ2A),即38.45<σ2<132.15时,系统存在正实特征值,则该部分存在不稳定解,其它部分为1个稳定解.与图4(c)相对应的,图4(a)和图4(b)分别给出了第一阶模态和第三阶模态的幅频相应图.同时,结合表2可知由A和C点计算得到的Jacobian矩阵特征值都具有负实部,可以判定A和C点所在的曲线为稳定解用实线表示;由B点计算得到的Jacobian矩阵特征值实部中有两个具有正实部,可以判定B点所在的曲线为不稳定解用虚线表示.

  • 表2 图4标注各点所对应的Jacobian矩阵特征值

  • Table2 Figure4 shows the eigenvalues of Jacobian matrix corresponding to each point

  • 运用同样的方法分别判断了系统在情况Ⅱ和情况Ⅲ下的稳定性.从图5(a)和图5(b)可以看出在区间(σ2Bσ2A),即38.16<σ2<131.67时,系统存在2个稳定解和1个不稳定解,其它部分为稳定解.从图6(a)和图6(b)可以看出在区间(σ2Aσ2C)和(σ2Fσ2B),系统存在2个稳定解和1个不稳定解;在区间(σ2Cσ2E),系统存在3个稳定解和2个不稳定解;在区间(σ2Eσ2D),系统存在3个稳定解和4个不稳定解;在区间(σ2Dσ2F),系统存在2个稳定解和3个不稳定解,其它区间只存在1个稳定解.

  • 图5 k3=8×109.9N/m4时,幅频响应及其稳定性(情况Ⅱ)

  • Fig.5 Amplitude-frequency response and its stability at k3=8×109.9N/m4 (Case Ⅱ)

  • 图6 k3=8×109.9N/m4时,幅频响应及其稳定性(情况Ⅲ)

  • Fig.6 Amplitude-frequency response and its stability at k3=8×109.9N/m4 (Case Ⅲ)

  • 3 结论

  • 研究了在移动荷载作用下非线性粘弹性基础梁的振动问题.采用Galerkin方法对非线性偏微分运动方程进行离散,并利用多尺度法求解了非线性运动控制方程.在此基础上,研究了非线性基础梁在3:1 内共振下的自由振动,以及在内共振作用下梁受外激励作用的第一、第二和第三阶模态主共振的非线性动力学行为及其动力学响应的稳定性问题.研究揭示了基础参数对内共振下弹性梁振动响应的影响规律,并发现在发生第三阶主共振时,幅频响应中会出现7个解支,包括3个稳定解支和4个不稳定解支.

  • 参考文献

    • [1] GHAYESH M H,KAZEMIRAD S,AMABILI M.Coupled longitudinal-transverse dynamics of an axially moving beam with an internal resonance [J].Mechanism and Machine Theory,2012,52:18-34.

    • [2] HU W P,YE J,DENG Z C.Internal resonance of a flexible beam in a spatial tethered system [J].Journal of Sound and Vibration,2020,475:115286.

    • [3] DING H,HUANG L L,MAO X Y,et al.Primary resonance of traveling viscoelastic beam under internal resonance [J].Applied Mathematics and Mechanics-English Edition,2017,38(1):1-14.

    • [4] WANG Y R,WEI Y H.Internal resonance analysis of a fluid-conveying tube resting on a nonlinear elastic foundation [J].European Physical Journal Plus,2020,135(4):364.

    • [5] WANG G X,DING H,CHEN L Q.Dynamic effect of internal resonance caused by gravity on the nonlinear vibration of vertical cantilever beams [J].Journal of Sound and Vibration,2020,474:115265.

    • [6] GUILLOT V,GIVOIS A,COLIN M,et al.Theoretical and experimental investigation of a 1:3 internal resonance in a beam with piezoelectric patches [J].Journal of Vibration and Control,2020,26(13-14):1119-1132.

    • [7] DING H,HUANG L L,DOWELL E,et al.Stress distribution and fatigue life of nonlinear vibration of an axially moving beam [J].Science China Technological Sciences,2019,62(7):1123-1133.

    • [8] YOUNESIAN D,SADRI M,ESMAILZADEH E.Primary and secondary resonance analyses of clamped-clamped micro-beams [J].Nonlinear Dynamics,2014,76(4):1867-1884.

    • [9] 顾伟,张博,丁虎,等.2:1内共振条件下变转速预变形叶片的非线性动力学响应 [J].力学学报,2020,52(7):1131-1142.GU W,ZHANG B,DING HU,et al.Nonlinear dynamic response of pre-deformed blade with variable rotational speed under 2:1 internal resonance [J].Chinese Journal of Theoretical and Applied Mechanics,2020,52(7):1131-1142.(in Chinese)

    • [10] 毕勤胜,陈予恕.双摆内共振分岔分析 [J].应用数学和力学,2000,21(3):226-234.BI Q S,CHEN Y S.Bifurcation analysis of a double pendulum with internal resonance [J].Applied Mathematics and Mechanics,2000,21(3):226-234.(in Chinese)

    • [11] 魏明海,肖仪清.内外共振联合作用下索-梁组合结构非线性振动分析 [J].振动与冲击,2012,31(7):79-84.WEI H M,XIAO Y Q.Nonlinear vibration analysis for a cable-beam coupled system under simultaneous internal and external resonances [J].Vibration and Shock,2012,31(7):79-84.(in Chinese)

    • [12] LI S H,YANG S P,CHEN L Q,et al.Effects of parameters on dynamic responses for a heavy vehicle-pavement-foundation coupled system [J].International Journal of Heavy Vehicle Systems,2012,19(2):207-224.

    • [13] DING H,CHEN L Q.Approximate and numerical analysis of nonlinear forced vibration of axially moving viscoelastic beams [J].Acta Mechanica Sinica,2011,27(3):426-437.

    • [14] DONG Z J,MA X Y.Analytical solutions of asphalt pavement responses under moving loads with arbitrary non-uniform tire contact pressure and irregular tire imprint [J].Road Materials and Pavement Design,2018,19(8):1887-1903.

    • [15] 杨予,滕念管,黄醒春,等.承受移动均布质量的简支梁振动反应分析 [J].振动与冲击,2005,24(3):19-22+6.YANG Y,TENG N G,HUANG X C,et al.Vibration analysis of a simply supported beam traversed by uniform distributed moving mass [J].Vibration and Shock,2005,24(3):19-22+6.(in Chinese)

    • [16] 刘涛,周洋忻,胡伟鹏.轴向运动功能梯度梁横向振动问题的保结构分析 [J].动力学与控制学报,2022,20(6):101-105.LIU T,ZHOU Y X,HU W P.Structure-preserving analysis on transverse vibration of functionally graded beam with an axial velocity [J].Journal of Dynamics and Control,2022,20(6):101-105.(in Chinese)

    • [17] MARTINEZ-CASTRO A E,MUSEROS P,CASTILLO-LINARES A.Semi-analytic solution in the time domain for non-uniform multi-span Bernoulli-Euler beams traversed by moving loads [J].Journal of Sound and Vibration,2006,294(1-2):278-297.

    • [18] 张弛,毛晓晔,丁虎,等.受轴向激励弹性支承梁的稳定性分析 [J].动力学与控制学报,2022,20(3):66-76.ZHANG C,MAO X Y,DING H,et al.Stability analysis of axially excited beam with elastic boundary [J].Journal Dynamics and Control,2022,20(3):66-76.(in Chinese)

    • [19] YOUNESIAN D,MARJANI S R,ESMAILZADEH E.Nonlinear vibration analysis of harmonically excited cracked beams on viscoelastic foundations [J].Nonlinear Dynamics,2013,71(1-2):109-120.

    • [20] DING H,YANG Y,CHEN L Q,et al.Vibration of vehicle-pavement coupled system based on a Timoshenko beam on a nonlinear foundation [J].Journal of Sound and Vibration,2014,333(24):6623-6636.

    • [21] ANSARI M,ESMAILZADEH E,YOUNESIAN D.Frequency analysis of finite beams on nonlinear Kelvin-Voight foundation under moving loads [J].Journal of Sound and Vibration,2011,330(7):1455-1471.

    • [22] SAPOUNTZAKIS E J,KAMPITSIS A E.Nonlinear response of shear deformable beams on tensionless nonlinear viscoelastic foundation under moving loads [J].Journal of Sound and Vibration,2011,330(22):5410-5426.

  • 参考文献

    • [1] GHAYESH M H,KAZEMIRAD S,AMABILI M.Coupled longitudinal-transverse dynamics of an axially moving beam with an internal resonance [J].Mechanism and Machine Theory,2012,52:18-34.

    • [2] HU W P,YE J,DENG Z C.Internal resonance of a flexible beam in a spatial tethered system [J].Journal of Sound and Vibration,2020,475:115286.

    • [3] DING H,HUANG L L,MAO X Y,et al.Primary resonance of traveling viscoelastic beam under internal resonance [J].Applied Mathematics and Mechanics-English Edition,2017,38(1):1-14.

    • [4] WANG Y R,WEI Y H.Internal resonance analysis of a fluid-conveying tube resting on a nonlinear elastic foundation [J].European Physical Journal Plus,2020,135(4):364.

    • [5] WANG G X,DING H,CHEN L Q.Dynamic effect of internal resonance caused by gravity on the nonlinear vibration of vertical cantilever beams [J].Journal of Sound and Vibration,2020,474:115265.

    • [6] GUILLOT V,GIVOIS A,COLIN M,et al.Theoretical and experimental investigation of a 1:3 internal resonance in a beam with piezoelectric patches [J].Journal of Vibration and Control,2020,26(13-14):1119-1132.

    • [7] DING H,HUANG L L,DOWELL E,et al.Stress distribution and fatigue life of nonlinear vibration of an axially moving beam [J].Science China Technological Sciences,2019,62(7):1123-1133.

    • [8] YOUNESIAN D,SADRI M,ESMAILZADEH E.Primary and secondary resonance analyses of clamped-clamped micro-beams [J].Nonlinear Dynamics,2014,76(4):1867-1884.

    • [9] 顾伟,张博,丁虎,等.2:1内共振条件下变转速预变形叶片的非线性动力学响应 [J].力学学报,2020,52(7):1131-1142.GU W,ZHANG B,DING HU,et al.Nonlinear dynamic response of pre-deformed blade with variable rotational speed under 2:1 internal resonance [J].Chinese Journal of Theoretical and Applied Mechanics,2020,52(7):1131-1142.(in Chinese)

    • [10] 毕勤胜,陈予恕.双摆内共振分岔分析 [J].应用数学和力学,2000,21(3):226-234.BI Q S,CHEN Y S.Bifurcation analysis of a double pendulum with internal resonance [J].Applied Mathematics and Mechanics,2000,21(3):226-234.(in Chinese)

    • [11] 魏明海,肖仪清.内外共振联合作用下索-梁组合结构非线性振动分析 [J].振动与冲击,2012,31(7):79-84.WEI H M,XIAO Y Q.Nonlinear vibration analysis for a cable-beam coupled system under simultaneous internal and external resonances [J].Vibration and Shock,2012,31(7):79-84.(in Chinese)

    • [12] LI S H,YANG S P,CHEN L Q,et al.Effects of parameters on dynamic responses for a heavy vehicle-pavement-foundation coupled system [J].International Journal of Heavy Vehicle Systems,2012,19(2):207-224.

    • [13] DING H,CHEN L Q.Approximate and numerical analysis of nonlinear forced vibration of axially moving viscoelastic beams [J].Acta Mechanica Sinica,2011,27(3):426-437.

    • [14] DONG Z J,MA X Y.Analytical solutions of asphalt pavement responses under moving loads with arbitrary non-uniform tire contact pressure and irregular tire imprint [J].Road Materials and Pavement Design,2018,19(8):1887-1903.

    • [15] 杨予,滕念管,黄醒春,等.承受移动均布质量的简支梁振动反应分析 [J].振动与冲击,2005,24(3):19-22+6.YANG Y,TENG N G,HUANG X C,et al.Vibration analysis of a simply supported beam traversed by uniform distributed moving mass [J].Vibration and Shock,2005,24(3):19-22+6.(in Chinese)

    • [16] 刘涛,周洋忻,胡伟鹏.轴向运动功能梯度梁横向振动问题的保结构分析 [J].动力学与控制学报,2022,20(6):101-105.LIU T,ZHOU Y X,HU W P.Structure-preserving analysis on transverse vibration of functionally graded beam with an axial velocity [J].Journal of Dynamics and Control,2022,20(6):101-105.(in Chinese)

    • [17] MARTINEZ-CASTRO A E,MUSEROS P,CASTILLO-LINARES A.Semi-analytic solution in the time domain for non-uniform multi-span Bernoulli-Euler beams traversed by moving loads [J].Journal of Sound and Vibration,2006,294(1-2):278-297.

    • [18] 张弛,毛晓晔,丁虎,等.受轴向激励弹性支承梁的稳定性分析 [J].动力学与控制学报,2022,20(3):66-76.ZHANG C,MAO X Y,DING H,et al.Stability analysis of axially excited beam with elastic boundary [J].Journal Dynamics and Control,2022,20(3):66-76.(in Chinese)

    • [19] YOUNESIAN D,MARJANI S R,ESMAILZADEH E.Nonlinear vibration analysis of harmonically excited cracked beams on viscoelastic foundations [J].Nonlinear Dynamics,2013,71(1-2):109-120.

    • [20] DING H,YANG Y,CHEN L Q,et al.Vibration of vehicle-pavement coupled system based on a Timoshenko beam on a nonlinear foundation [J].Journal of Sound and Vibration,2014,333(24):6623-6636.

    • [21] ANSARI M,ESMAILZADEH E,YOUNESIAN D.Frequency analysis of finite beams on nonlinear Kelvin-Voight foundation under moving loads [J].Journal of Sound and Vibration,2011,330(7):1455-1471.

    • [22] SAPOUNTZAKIS E J,KAMPITSIS A E.Nonlinear response of shear deformable beams on tensionless nonlinear viscoelastic foundation under moving loads [J].Journal of Sound and Vibration,2011,330(22):5410-5426.

  • 微信公众号二维码

    手机版网站二维码