【工程力学】简谐荷载下斜拉索的超谐波共振分析

毕业设计（论文）题目简谐荷载下斜拉索的超谐波共振分析作者学院土木工程学院专业工程力学学号指导教师二一五年五月廿七日湖南科技大学毕业设计（论文）任务书土木工程学院院工程力学系（教研室）系（教研室）主任（签名）年月日学生姓名学号专业工程力学1设计（论文）题目及专题简谐荷载作用下斜拉索的3次超谐波共振分析2学生设计（论文）时间自2014年12月30日开始至2015年5月30日止3设计（论文）所用资源和参考资料相关专业的教材和专著；中文专业期刊上的相关论文；网上电子资源，例如外文资源库SPRINGERLINK全文数据库。4设计（论文）应完成的主要内容1查阅非线性振动理论中摄动法的相关文献，通过毕业论文的整个过程学会获取文献与追踪文献；2应用文献中的摄动法求解简谐荷载作用下斜拉索的3次超谐波共振的近似周期解与平均化方程；3算例分析，根据具体参数进行算例分析、讨论，得出论文的主要结论。5提交设计（论文）形式（设计说明与图纸或论文等）及要求1）毕业论文要求在5月30日以前间完成并提交；2）按照学校规定撰写毕业论文，包括论文总体格式、字体与字体大小、参考文献的书写格式、论文装订格式等。6发题时间2014年12月30日指导教师（签名）学生（签名）湖南科技大学毕业设计（论文）指导人评语指导人（签名）年月日湖南科技大学毕业设计（论文）评阅人评语评阅人（签名）年月日湖南科技大学毕业设计（论文）答辩记录日期学生学号班级题目提交毕业设计（论文）答辩委员会下列材料1设计（论文）说明书共页2设计（论文）图纸共页3指导人、评阅人评语共页毕业设计（论文）答辩委员会评语答辩委员会主任（签名）委员（签名）（签名）（签名）（签名）成绩摘要采用哈密顿变分原理，根据拟静态假设，建立斜拉索面内非线性运动方程，同时采用GALERKIN方法对斜拉索面内非线性运动微分方程进行离散，利用多尺度法分别对所得到的方程进行摄动分析，以引入新的时间尺度，将其代入非线性运动方程中消去久期项，设立新的方程，利用欧拉公式化简再经过分离实部虚部求解在稳态运动下方程的解，进而得到的斜拉索超谐波共振响应前三阶模态的的近似解。最后，利13用MATLAB和MATHMATIC对近似解进行数值计算，得到斜拉索面内非线性运动的相图曲线、时程曲线等并对其进行分析。关键词斜拉索；多尺度法；的超谐波共振；幅频响应曲线；时程响应曲线13ABSTRACTBYUSINGTHEHAMILTONSVARIATIONALPRINCIPLE,ACCORDINGTOTHEQUASISTATICASSUMPTION,THEESTABLISHMENTOFCABLESURFACEINNONLINEAREQUATIONSOFMOTIONATTHESAMETIME,THEGALERKINMETHODOFCABLEPLANENONLINEARMOTIONDIFFERENTIALEQUATIONSAREDISCRETIZED,RESPECTIVELYOFTHEOBTAINEDEQUATIONPERTURBATIONANALYSISUSINGTHEMETHODOFMULTIPLESCALES,BYINTRODUCINGANEWTIMESCALEWILLBESUBSTITUTEDINTOTHENONLINEAREQUATIONSOFMOTIONFORTHEELIMINATIONOFSECULARTERMS,THEESTABLISHMENTOFTHENEWEQUATIONANDEULERSFORMULAWASUSEDTOSIMPLIFYAGAINAFTERSEPARATIONREALPARTOFTHEIMAGINARYPARTOFTHESOLUTIONINSTEADYMOTIONEQUATIONSOLUTION,ACABLEANDTHESUPERHARMONICRESONANCERESPONSEANAPPROXIMATESOLUTIONFORTHETHIRDORDERMODESATLAST,THEPHASEDIAGRAMCURVESANDTIMEHISTORYCURVESOFTHENONLINEARMOTIONOFTHECABLEPLANEAREANALYZEDBYUSINGMATLABANDMATHMATICKEYWORDSSUSPENSIONCABLEMULTISCALEMETHODTHESUPERHARMONICRESONANCETHEAMPLITUDEFREQUENCYRESPONSECURVETIMEHISTORYRESPONSECURVE目录第一章前言1第二章斜拉索的运动微分方程321数学模型322离散化方程6第三章方程求解831多尺度法832代入方程求方程的通解933在稳态条件下求解方程1234得到方程的近似解12第四章数值计算分析1341基本数据1342幅频响应曲线1443时程曲线和相图曲线16第五章结论19参考文献20致谢21第一章前言本文对简谐荷载作用下斜拉索面内单自由度非线性超谐波共振进行了研究，利用哈密顿变分原理得到斜拉索面内单自由度位移振动微分方程，并采用多尺度法进行求解计算，通过数值分析，得到了斜拉索非线性超谐波共振下的多组曲线，包括相图，时程曲线，频率响应曲线等，在此基础上讨论简谐荷载下斜拉索超谐波共振振动的机制，利用这些图像数据相互补充说明，证明了本文所采用数据的合理性和精确性，并在此基础上得出相应结论。近年来，科学技术的进步以及实际施工（斜拉桥方面）的需要，促进了对斜拉索这种大幅振动的一系列实验和理论的研究，吴晓1在斜拉索非线性振动固有特性方面做了细致的研究，将斜拉索静平衡时索曲线纳入考虑范围,并在此基础上推导斜拉桥索大幅振动的非线性动力方程,利用傅立叶级数法研究与求解了斜拉索非线性运动的固有振动方程王连华2研究了斜拉索倾斜角对拉索周期运动的影响。对拉索位移离散空间动力学模型进行过深入研究，利用打靶法研究拉索的空间运动，探讨斜拉索在支撑运动作用下斜拉索大幅运动的原因，与此同时，在其论文中得出给出了相应结论，倍周期失稳和1L内共振强烈激发的模态非线性作用可能导致拉索的非平面运动。此外，倾角对斜拉索大幅振动有较大影响，小倾角斜拉索，竖向运动的惯性力是引起斜拉索大幅振动的主要因素；相反倾斜角较大时，拉索的大幅振动受轴向运动惯性力导致的参数激励控制。刘圣波3建立了拉索面内非线性振动模型，在充分考虑斜拉索重力现象分力垂度曲线的基础上，同时承受塔梁端部激励作用，并对其进行了分析研究。讨论了影响斜拉索非线性振动的诸多因素，其中主要包括初始拉力、激励幅值，并与按传统的结果进行了对比；熊蓉4在名为斜拉索非线性振动分析的论文中分析了斜拉索非线性振动的影响因素，讨论了非线性的振动机理，并且探讨了什么影响了斜拉索非线性振动，怎样影响的振动，做出了相当有意义的总结；袁从森5对斜拉索振动进行了研究，讨论了考虑索力变化对斜拉索振动的影响，解决了考虑弦向分布的重力分量时，计算以往计算精度不足问题；张剑6在九江长江公路大桥斜拉索振动特性研究中，基于实际的条件上，得出这样的结论阻尼比是否合适是是否能避免斜拉索在超谐波共振区产生大幅振动的关键因素，虽然此项研究得到的各斜拉索内共振和参数共振的共振区是基于九江长江大桥主梁形成的端部位移，但任具有一定普遍性，具有非凡意义；彭然7研究了桥面激励下斜拉索非线性耦合振动，他所建立的桥索塔振动模型模拟了桥面激励下斜拉索非线性振动，研究斜拉索在桥面激励下的非线性振动特性，采用伽辽金法离散控制方程，将理想条件下的系统频率比、桥面激励幅值、桥塔初始扰动、拉索初始索力、倾角等因素纳入考虑范围，探索它们对斜拉索面内振动的影响，在考虑阻尼的非理想条件下研究拉索面内振动特性；基于对斜拉桥全桥前30阶振动频率与全桥各拉索自振基频进行匹配的数据基础上，对可能发生共振的拉索及系统的振动特性进行分析。相应的，在国外，也有大量学者对斜拉索非线性做过大量研究，当然，由于斜拉索和斜拉索有很多的相似性，国内外对斜拉索的研究也从另一方面促进了斜拉索的发展，在此不一一列出。第二章斜拉索的运动微分方程21数学模型建立如图坐标系，分为直角坐标系，曲坐标系，分别用和来表IOX0VS、示，在这里我们可以把弹性斜拉索的变形分多个阶段进行描述，如下1自然无限伸长位置；0,QTQTQT引进一些新的自变量，因此关于T的导数变成为关于偏导数的展开式，即N201201,DTDDTTTT3600000121122232QDTQDTQDTTTTTTTTDDQ关于T的二阶偏导数22001102DT37210111012222300110QDQDQDQDQ将353637的结果代入33式得38222232200110101101223301COSQQAQAKT令两段，的系数分别相等得013902000CSDQT31013112DQA31122300002102010DQQA32代入方程求方程的通解式38为非其次线性方程，其通解为他所对应的一个特解和他所对应得齐次方程的通解的和是一个二阶齐次微分方程它关于的特征方程求200DQ0T解如下有共轭复根所以其通解为20R1,20RI3121020COSSINQATCIT其中是关于的任意函数。假定某复数使得1,AT102000COSSINEXPEXPCIIAIT则齐次方程的通解可以写成3130120120,QATITI求特解所以可以设特解将特解代入方程0COSFTK0COSQB可得200,DQT31422000CSCOSSBK解得则特解20KB020Q则通解为315001201202COS,EXP,EXPQATIATIT整理变化后得3160100100EXPEATIIII其中20K将代入310得0Q317210100222000EXPEXPDIDAITIITIAIC其中CC指前面各项的共轭复数315的解为消去久期项后各特解的叠加，先求的特解2010EXPDQIIT令代入上式可得10EXPBIT120IB31820120EXPQAIT令代入上式可得120EXPQBIT220AB33192012DQA令代入上式可得13Q32032020120EXPQIT令代入上式可得140EXP2QBI240B3212201DQ令代入上式可得150EXPQI2503222012EXPQAIT令代入上式可得160EXPQBIT260AB3232012EXPDQIT令代入上式可得170EXPQIT270将各特解叠加可得324222100020022002000EEXPEXP4AIIIAIITITC3将代入（311）消去中的久期项得Q12Q325202045628EXPIDAIAI其中32623402225300263014（316）式中我们设327121212,EXP,ATATI（这里的，都是实数,且根据欧拉公式）ASINCOS代入后3283230224563230224561EXPCOSSIDAIAAIII令化简并将实部和虚部分开我们得到2T32930063230204561SINCOSDAA33在稳态条件下求解方程当时得到将329两式平方后再相加可得20D262442426045050514AAAAA化简得22624222455540000001A22622540014A330226225400A34得到方程的近似解所以作为一次近似解33122022201COS3COSCOSCOS3461QTATTTATA式中33212245TAN/第四章数值计算分析41基本数据在这里，我们把研究对象取为湖南省洞庭湖大桥中的斜拉索之一作为研究对象进行数值分析。这个斜拉索的几何物理参数为索跨直径129LM19DM初始张力弹性模量350HKN520EMP倾斜角单位长度质量O8/KG重力加速度外部荷载振幅2981/GS1激励频率0根据资料，斜拉索的阻尼由粘性阻尼和气动阻尼，其中气动阻尼由于其受风速和方向影响，难于测量，所以不考虑它，且由于此斜拉索实测十阶模态阻尼比量级小，大约，所以为方便研究，本文阻尼比取值为00126；同时在此说明，本文1中的不稳定解不给出。与此同时王连华博士也曾测试与计算过这座桥的自振频率其结果是无量纲化后的结果，如下表所示表41斜拉索面内模态固有频率（摘自王连华博士的论文）面内模态摘录论文IRVINE理论实测值本文无量纲化前数据11015710157107PI220229202292142PI330346303463203PI440459404594234PI无量纲化后，本文各阶自振频率数值分别为1、2、3、4，可以看出，本文所用数据和理论计算出的数据是基本吻合的，而且与实测数据极为接近，同时，本文只研究一阶模态，故可取固有频率值为PI。根据第二章节的介绍及以上介绍的数据，利用MATHMATIC容易计算出运动微分方程中的关键参数二次项系数三次项系数2589716319842幅频响应曲线根据式3,52，并利用MATHMATIC计算数据后，通过软件ORINGIN画图,在这个频率响应曲线中，我们可以明显的看到，简谐荷载下斜拉索超谐波共振表现出硬弹簧性质，曲线向右弯曲，最终两条曲线会趋于一致。图41超谐波共振幅频响应曲线（K0001）0200204005010150202503AK08图42超谐波共振幅频响应曲线（K00008）0200204060050101502ABDK01002040608005010150202503AK02图43超谐波共振幅频响应曲线（K0002）下面两幅图K的值依次是00008和0002这两幅图是改变K的大小后得出的图，可以看出，外激励幅值并未改变其模态响应曲线的外形，所以外激励并非引起非线性的因素。与此同时，值得注意的是，外激励幅值却在一定程度上改变了调谐参数和振幅的数值关系，具体表现为随着K值得不断增大，左右两条曲线会在越来越远的地方交汇，偏转程度也有着一定程度的减缓，从而使得图像的偏转越来越集中调谐参数趋于零的区域，图像的跳跃显然在逐渐增加。所以，外激励幅值虽然不引起非线性，但外激励幅值确实影响非线性因素的重要因素。43时程曲线和相图曲线我们可以通过MATLAB轻松获得斜拉索面内单自由度的时程曲线和相图曲线，它们分别代表时间位移关系和速度位移曲线。0501001502002151050051152X103TIMEHISTORY21012X10354321012345X103PHASEDIAGRAM图44超谐波共振时程响应曲线（K0005）050100150200864202468X104TIMEHISTORY1050051X1032151050051152X103PHASEDIAGRAM图45超谐波共振时程响应曲线（K0002）05010015020054321012345X104TIMEHISTORY505X1041080604020020406081X103PHASEDIAGRAM图46超谐波共振时程响应曲线（K0001）05010015020015105005115X103TIMEHISTORY21012X1033210123X103PHASEDIAGRAM图47超谐波共振时程响应曲线（K0003）以上四幅图的K值分别是0005,0002,0001,0003，从这连续的四幅图中可以看出外激励幅值对斜拉索大幅振动的影响。从时程曲线看，尽管外激励幅值不断变化，振动趋向稳定的时间基本保持不变，或者说，变化幅度小，只是稳定阶段的振幅大小随着外激励的增加而增加；然而不稳定阶段的变化却很明显，随着外激励幅值的逐渐增大，不稳定阶段的振幅颠簸程度逐渐趋于缓和。从相图曲线看，变化也很明显，当外激励幅值上升时，相图曲线的收敛趋于逐渐外延，这是和时程曲线一致的，说明稳定阶段振幅逐渐增大了；同时曲线的外围逐渐紧凑，而内部逐渐变得系数，成向外扩散的环状，而且环形间距随着外激励幅值大而增大，这就表示，图像的收敛速度随着外激励幅值增大不断上升，这就给时程曲线做了个补充说明。第五章结论本文利用哈密顿变分原理，引入拟静态假设，建立了斜拉索面内振动微分方程接着利用GALERTKIN法对其进行离散；然后，利用多尺度法计算第一阶模态超谐波共振响应的近似解以及图像，通过数学分析和所对应的图像观察，通过计算分析可以到得出引起非线性振动的原因不是外激励荷载，调谐参数可以在一定的范围内对斜拉索的超谐波共振产生影响，其振动过程可以分为不稳定阶段和稳定阶段，而且外激励越大到达稳定阶段的时间就越长，在荷载作用下斜拉索的超谐波共振的响应振幅在一定的范围内不会因为初始条件的不同而发生变化，