轴向受力梁强非线性超谐波与次谐波共振的能量迭代法
周一峰1,唐进元2, 何旭辉1
(1.中南大学 土木建筑学院,湖南 长沙,410075;
2.中南大学 机电工程学院,湖南 长沙,410083)
摘要: 研究了在横向动载荷作用下存在轴向力的一类梁的强非线性振动;建立了梁振动的二阶强非线性非自治微分方程,并对求解强非线性自治系统的能量迭代法加以改进,用于求解梁的强非线性非自治系统,其方法是:由能量法得到主共振、超谐共振和次谐共振的一次近似解的表达式;引入牛顿迭代的思想和最小二乘法,得到高次近似解的表达式。研究结果表明:用改进后的能量迭代法求解强非线性非自治系统精度较高;分析这种非线性梁的振动时,除了要考虑其主共振外,还要考虑超谐共振和次谐共振。
关键词: 梁; 强非线性非自治系统; 超谐共振; 次谐共振; 能量迭代法
中图分类号:O322 文献标识码:A 文章编号: 1672-7207(2005)04-0698-06
Energy-iterative method for solution of super-harmonic and sub-harmonic resonance of strongly non-linear vibrations with an axially forcing beam
ZHOU Yi-feng1, TANG Jing-yuan2, HE Xu-hui1
(1.School of Civil and Architectural Engineering, Central South University, Changsha 410075, China;
2.School of Mechanical and Electrical Engineering, Central South University, Changsha 410083, China)
Abstract: Strongly non-linear vibrations of an axially forcing beam under transverse dynamic loads were investigated. The Duffing equation, the equation of strongly non-linear and non-autonomous oscillation system of the beam were set up. The energy-iteration method for strongly non-linear autonomous system was modified to deal with strongly non-linear non-autonomous system, and the modified new method was developed, the analytic expressions of solution of primary resonance and superharmonic and subharmonic resonance by the energy method were obtained, and the higher order analytic expression of solution was gotten by means of Newton iteration idea and least square technique principle.The results show that the modified energy iterative method for strongly non-linear and non-autonomous system is both effective and accurate, and that not only the primary resonance but also the superharmonic and subharmonic resonance must be considered in this sort of strongly non-linear beams.
Key words: beam; strongly non-linear and non-autonomous system; superharmonic resonance; subharmonic resonance; energy iterative method
在结构和机械振动中,许多模型可以归结为在横向动载荷作用下有轴向力的非线性振动,如带锯、输流管道、空中索道、纸张、皮带、纺织纤维、磁带和桥梁等可归结为弦或梁的模型[1]。A.H.Nayfeh等[2-4]归纳了弱非线性振动的各种摄动法,并研究了梁的主共振以及超谐、次谐共振和多个谐波的联合共振等各种弱非线性振动,用多尺度法研究了梁的弱非线性振动的内共振、组合谐波共振;L.Azrar等[5]用谐波平衡法研究梁的非线性自由振动及强迫振动;H.R.z[6]用多尺度法研究了轴向速度为变速度,且ε1时梁的弱非线性振动。
随着科学技术的迅速发展,梁的非线性振动尤其是强非线性振动问题日益突出。弱非线性振动分析的方法不能适用于强非线性振动问题的分析。为此,曾庆元等[7,8]建立了桥梁空间振动的有限元模型,用数值方法求解桥梁的非线性振动; 黄建亮等[9]用半解析的改进的L-P法,研究了轴向力为常量时梁的强非线性振动的超谐与次谐波解[9],并给出了非线性项系数ε〈1时的算例。能量—迭代法[10,11]在应用能量法 [12-14]求解的基础上进一步迭代,不仅能对系统进行定性分析,而且定量求解的结果精度较高。能量迭代法[10,11]仅能求解强非线性自治系统,在此,作者对其进行改进,推广应用于求解轴向力为变量时梁的强非线性非自治振动,求解当非线性项系数ε1时的主共振及超谐波与次谐波共振解。
1 梁的非线性振动微分方程
图1所示为受横向动载荷和轴向拉力的两端简支弹性弯曲梁结构图。设梁的长度为L,横截面积为A,横截面对中性轴的惯性矩为I,质量密度为ρ,抗弯刚度为EI,回转半径为r=[KF(]I/A[KF)],梁的横向载荷为q(x, t),其振动方程为:
图 1 受横向动载荷和轴向拉力的梁的结构图
Fig. 1 Beam under axially force and transverse dynamic loads
其中,S为轴力,其一般表达式为:
轴力包括2部分:EA·(u/x)和EA·(y/x)2/2。当u/x=o((y/x)2,即轴向力主要由梁横截面的转角位移引起时,可以略去梁轴向的应变对轴力的影响,即 ,故
将式(2)代入式(1),得:
采用基本模态 ,代入式(3),有
引入无量纲量:
方程(4)可以写为:
其中:非线性项3bx3/4即为轴力的影响项。取q(x, t)=hcosΩt·sin[SX(]πx[]L[SX)](其中,Ω为动载荷变化的频率),有
其中:P=h/(ρAr)。式(7)即梁的强非线性无阻尼受迫振动Duffing方程。此系统除有主共振响应外,还可能有超谐、次谐及超次谐共振响应[15-19]。
2 梁的主振动和超、次谐振动周期解的能量迭代解法
式(7)为二阶强非线性非自治系统,可以表达为一般形式:
对上式两端同乘dx后积分,得
式中:V(x)=∫g(x)dx,为系统的势能函数。式(9)表明系统的机械能变化等于非有势力所做的功。因为求周期解,系统在1个周期内平均能量守恒,应有:
取1个周期内的平均值x为其一次近似解,有
根据式(9)和式(11),有
设一次近似周期解为:
初始条件为:
故
将式(13)和式(14)代入式(11)得:
令
有
2.1 梁的主振动周期解
令
有
由此解得:
当Ω2>b时,a*为正实数,有周期解。
取一组数据:E=200 GPa,ρ=7855.7 kg/m, , L=3.46 m,P=1,则有b=17.232。对应于文献[9]中的ε,本例ε=3b/4=12.9241,此时取,存在主振动周期解。将Ω和b代入式(19)得:a*=2.56。为确定式(13)第2式中相应于a*的初相位φ*,令式(12)的被积函数中的非周期项为零,并将式(15)和式(17)代入式(12),令其非周期项为零,有
得φ*=0。故主振动一次近似周期解为:
为求二次近似解x2,将式(7)写为:
设
将式(21)中F(x2,[AKx·]2,Ωt)在(x1,[AKx·]1)邻域内展开为泰勒级数:
将式(22)和式(23)代入式(21),有
对照式(7), 。故式(24)即为
将x1=a*cosΩt代入式(25),有
这是关于w1的二阶周期系数的非齐次方程。按谐波平衡思想,令
由式(13)及φ*=0,有
将式(27)(取N=3)代入式(26),合并sinjΩt和cosjΩt各项,得:
采用最小二乘法,令
根据三角函数正交性,有
令
其中:i=1, 2, 3。
将式(32)写成矩阵方程形式:
其中:[M]为系数矩阵;{u}=(α0,α1,β1,…,α3,β3)T;{h}为常数列阵。解出{u},代入式(27)和式(22),即得二次近似的主振动周期解。算得
{u}=(0, -0.0385, 0, 0, 0, 0.0851, 0)T,
故二次近似主振动周期解表达式为:
用本方法求解所得结果与四阶Runge-kutta法数值解所得结果进行对比,结果见图2。
图 2 Duffing方程的主振动时间—位移图和相图
Fig. 2 Diagrams of time-displacement and phase of primary resonance of Duffing equation
2.2 梁的次谐振动周期解
令ω=Ω/3,求梁的1/3次谐振动周期解。据式(16),有
故当Ω2>9b时,a*有实数解。据前面的参数设置,Ω>12.45。取Ω=48,有
a*=4.963。
而式(7)中x3的系数不随ω的改变,ε=3b/4=12.9241。据式(12)被积函数中非周期项系数为零的条件,得sin3φ*=0,可取
φ*=π/3。
故一次近似1/3次谐解为
以下求二次近似的1/3次谐解。其方法与前面求二次近似解相同,只是此时
应用式(22)~(33),据式(36),将式(29)中的cosjθ和sinjθ展开,按cosjΩt和sinjΩt归并各项求得A0,Aj,Bj (j=1, 2, …, 5),最后求得{u}=(0,-0.1068,0,0,0,0.2015,0)T,代入式(22),得
x2=4.8292cos(Ωt/3+π/3)+0.2015(Ωt+π)。
该结果与四阶Runge-kutta法对照结果如图3所示。
图 3 Duffing方程1/3次谐振动的时间—位移图和相图
Fig. 3 Diagrams of time-displacement and phase of 1/3 subharmonic resonance of Duffing equation
2.3 梁的超谐振动周期解
令ω=3Ω, 。当Ω>1.384时,a*有实数解。取Ω=6,有a*=5.626,再由式(12)中被积函数的非周期项为零确定φ*,此时非周期项系数显然为零,故可取φ*=0。从而一次近似ω=3Ω超谐振动周期解为:
在求其二次近似解时,令θ=3Ωt。为使一次近似解谐波cos3Ωt得到修正,设w1为
按以上方法求A0,Aj和Bj (j=1, 2, …, 9)时,按jΩt的谐波归并各项,应用式(22)~(33),求得:
{u}=(0,0.0011,0,0,0,-0.0003,0)T。
即得
与四阶Rong-kutta法计算结果对照,结果见图4。
图 4 Duffing方程在ω=3Ω时的超谐振动时间—位移图和相图
Fig. 4 Diagrams of time-displacement and phase of sub harmonic resonance of Duffing equation at ω=3Ω
3 结 论
a. 应用能量原理与迭代方法求出强非线性轴向受力梁的主振动和超谐与次谐振动周期解。与四阶Runge-kutta法结果的对照图形可知,采用该方法可有效地求解该类梁的强非线性振动问题,在非线性项系数ε1时精度较高。
b. 由于考虑了轴力的影响,方程式(7)中出现强非线性项,使梁的振动不仅有主共振,还出现超谐共振、次谐共振,因此,在这类梁的设计中,不仅要考虑ω≈Ω的主共振,还应考虑各种超谐共振和次谐共振的影响。
参考文献:
[1]z H R, Pakdemirli M, Boyaci H. Non-linear vibrations and stability of an axially moving beam with time-dependent velocity[J]. Int J Non-linear Mechanics, 2001, 36: 107-115.
[2]Nayfeh A H, Mook D T. Nonlinear Oscillations[M]. New York: John Wily & Sons Interscience, 1979.
[3]Mook D T, Plant R H, Haquang N. The influence of an internal resonance on non-linear structural vibrations under subharmonic resonance conditions[J]. Journal of Sound and Vibration, 1985, 102: 437-492.
[4]Mook D T, Haquang N, Plant R H. The influence of an internal resonance on non-linear structural vibrations under combination resonance conditions[J]. Journal of Sound and Vibration, 1986, 104: 229-241.
[5]Azrar L, Benamar R, White R G. A semi-analytical approach to the non-linear dynamic response large vibration amplitudes (part I): general theory and application to the single mode approach to free and forced vibration analysis[J]. Journal of Sound and Vibration, 1999, 124(2): 183-207.
[6]z H R. Non-linear vibrations and stability analysis of tensioned pipes conveying fluid with variable velocity[J]. Int J Non-linear Mechanics, 2001, 36: 1031-1039.
[7]曾庆元,郭向荣.列车桥梁时变系统振动分析理论与应用[M].北京:中国铁道出版社,1999.
ZENG Qing-yuan, GU Xian-rong. Theory and application of vibration analysis with time variable for vehicle-beam system[M]. Beijing: China Railway Press,1999.
[8]郭文华,曾庆元.上承式钢板梁桥空间振动计算分析[J].振动与冲击, 2002, 21(1): 79-82.
GUO Wen-hua, ZENG Qing-yuan. Calculation and analysis spatial of vibration for deck steel plate girder bridges[J]. Vibration and Shock, 2002, 21(1): 79-82.
[9]黄建亮,黄惠仪,陈恒,等.梁的强非线性超、次谐波共振[J].振动与冲击,2004, 23(1): 1-3.
HUANG Jian-liang, HUANG Hui-yi, CHEN Heng, et al. Super-harmonic and sub-harmonic resonance of strongly non-linear vibrations of beam[J]. Vibration and Shock, 2004, 23(1): 1-3.
[10]周一峰,李骊. 一类强非线性系统周期解的能量迭代法[J].应用力学学报, 1997, 14(4): 26-30.
ZHOU Yi-feng, LI Li. The energy-iteration method for periodic solution of a class of strongly nonlinear system[J]. Chinese Journal of Applied Mechanics, 1997, 14(4): 26-30.
[11]周一峰. 强非线性振动系统周期解的能量迭代法[J]. 力学季刊, 2002, 23(3): 514-520.
ZHOU Yi-feng. The energy-iteration method for periodic solution of strongly nonlinear system[J].Chinese Quarterly of Mechanics, 2002, 23(4): 514-520.
[12]李骊. 强非线性振动系统的定性理论与定量方法[M]. 北京: 科学出版社,1997.
LI Li. The qualitative theory and quantitative analysis of nonlinear vibration systems[M]. Beijing: Science Press, 1997.
[13]LI Li. Energy method for computing periodic solutions of strongly nonlinear systems (Ⅰ)—autonomous systems[J]. Nonlinear Dynamics, 1996, 9(3): 223-247.
[14]LI Li. Energy method for approximate periodic solution of strongly nonlinear non-autonomous systems[J]. Nonlinear Dynamics, 1999, 12(3):237-261.
[15]黄安基.非线性振动[M].成都:西南交通大学出版社,1993.
HUANG An-ji. Nonlinear vibrations[M]. Chengdu: Southwest Jiaotong University Press, 1993.
[16]周纪卿,朱因远.非线性振动[M].西安:西安交通大学出版社,1998.
ZHOU Ji-qing, ZHU Ying-yuan. Nonlinear vibrations[M]. Xian: Xian Jiaotong University Press, 1998.
[17]刘延柱,陈立群.非线性振动[M]. 北京:高等教育出版社,2001.
LIU Yan-zhu, CHEN Li-qun. Nonlinear vibrations[M]. Beijing: Higher Education Press, 2001.
[18]张伟,杨绍普,徐鉴,等.非线性系统的周期振动和分岔[M].北京:科学出版社,2002.
ZHANG Wei, YANG Shao-pu, XU Jian, et al. The periodic solutions and their bifurcations of nonlinear vibration system[M]. Beijing: Science Press, 2002.
[19]刘延柱,陈立群.非线性动力学[M].上海:上海交通大学出版社,2000.
LIU Yan-zhu, CHEN Li-qun. Nonlinear dynamics[M]. Shanghai: Shanghai Jiaotong University Press, 2000.
收稿日期:2004-12-08
基金项目:国家自然科学基金资助项目(50475139);湖南省教委科研基金资助项目(湘教财字[1998]1号)
作者简介:周一峰(1946-),男,江苏东台人,副教授,从事非线性振动及一般力学研究
论文联系人: 周一峰,男,副教授; 电话:0731-2656176(H);E-mail:fyzhou@mail.csu.edu.cn