CFRP索斜拉梁面内自由振动建模及参数分析

康厚军 解维东++郭铁丁��

摘要：利用张紧弦和欧拉梁振动理论分别描述斜拉梁结构中索与梁的振动，通过索梁连接处的动态平衡条件，建立斜拉梁平面内自由振动理论.利用传递矩阵法和边界条件对斜拉梁结构平面内自由振动的特征值问题进行求解.同时，建立斜拉梁的有限元模型，有限元法所得结果与本文理论研究非常吻合，证明了本文理论和方法的正确性.最后对CFRP索斜拉梁平面内自由振动进行参数分析.研究表明，CFRP索斜拉梁的基本动力学性能优于传统钢索斜拉梁.

关键词：CFRP索；斜拉梁；传递矩阵法；振动分析；频率；振型

中图分类号： O343.9 文献标识码：A

碳纤维增强复合材料（Carbon Fiber Reinforced Polymer，简称CFRP）是由多股连续有机纤维丝在惰性气体中经高温炭化，并经拉挤成型技术和必要的表面处理而形成的一种新型复合材料.采用CFRP制成的拉索具有耐腐蚀性强、自重轻（仅为钢材的1/5左右）、强度高（钢材的8～10倍，弹性模量最高可达1 000 GPa，抗拉强度可达2 700 MPa[1]）、抗疲劳性能好等优点，相比传统钢拉索优势明显， 因此，CFRP斜拉索将有很好的应用前景.目前，国内外学者已从理论上证明了CFRP索相对于钢索的静动力特性有不同程度的改善[2-4]，CFRP索也已投入实际应用[5-6].截至目前国内外已建成CFRP索斜拉桥6座，其进一步的应用研究和基础研究已成为国内外研究的一个热点.我国已成功采用CFRP拉索替换钢拉索建造试验性质的人行斜拉桥[5]，未来斜拉桥也有采用CFRP拉索的趋势，尤其是对于特大跨径桥梁，CFRP索将具有足够的优势.然而，我国对于CFRP的研究还主要集中在应用加固方面，作为大跨度柔性结构，其动力学问题比较突出，相关研究却很少见到.

斜拉梁结构由于其良好的受力性能和优美的外观被广泛应用于土木工程和海洋工程，如斜拉桥、房屋建筑中的雨棚、塔吊以及桅杆结构等.由于斜拉梁中索和梁2种结构单元有着很大的力学差异，特别是索跟梁的耦合，历来是国内外学者研究的重点和难点.Fung[7]通过Hamilton原理和有限元法推导出的非线性时变微分方程研究了斜拉梁中索的长度和张力随时间变化的振动问题.Gattulli等人[8-9]通过经典变分公式得到了斜拉梁横向振动的运动控制方程，将其与有限元方法和试验进行对比，并考虑了面内和面外的振动；赵跃宇等人[10]利用索梁组合结构的连接条件和边界条件，建立了索梁组合结构的约化运动学控制方程，利用Galerkin模态截断得到了该系统的多模态离散动力学方程；Wang等人[11]通过Halmilton原理得到索梁组合结构的动力学运动方程，通过边界和连续性条件以及分离变量法，得到结构的频率方程和相应的振型表达式，并对固有频率进行了讨论.这些研究工作都只考虑了梁的横向振动，没有考虑纵向振动问题，并且在索梁连接条件的处理上各不相同，存在较大的局限性.

传递矩阵法（Transfer Matrix Method，简称TMM）是20世纪20年代建立起来的一种用矩阵来描述多输入多输出的线性系统的输出与输入之间关系的方法.相比于有限元方法，该方法计算精度不随划分段数而改变，许多学者和工程技术人员将传递矩阵法应用于解决工程实际问题，例如Kang和Wang等人[12-14]用传递矩阵法来研究索拱结构和悬索桥的动力学问题.

针对以上问题和方法，本文将同时考虑索和梁的纵横向振动，利用张紧弦和欧拉梁振动微分方程，在索梁结合处考虑它们的动态平衡并将索端和梁端内力和纵横向位移进行耦合，利用传递矩阵法求解系统振动的特征值问题.为了验证本文中索梁理论和传递矩阵法运用的正确性，我们将建立斜拉梁的有限元模型，对本文理论研究和有限元法结果进行对比，对本文的理论和求解方法进行验证.最后将对CFRP索斜拉梁的特征值问题进行参数分析，同时和传统钢索斜拉梁进行对比研究.

3特征值分析

为研究CFRP索斜拉梁的特征值问题，即固有频率和模态，选取如下物理参数：索为CFRP索，单位长度质量为10.4 kg/m，横截面积为6.273×10-3 m2，弹性模量为210 GPa，初始索力为1 MN，倾斜角度为30°；梁为钢筋混凝土箱梁，长100 m，单位长度质量为4.4×104 kg/m，横截面面积为16.3 m2，截面惯性矩为9.8 m4，弹性模量为34.5 GPa.

为了验证本文理论方法在斜拉梁结构中运用的正确性，我们用有限元软件ANSYS12.0建立了同样参数的斜拉梁有限元模型，其中索用Link1单元，梁用Beam3单元，划分单元数为200，然后比较本文理论和有限元法得到的频率和振型.表1分别列出了通过有限元法和本文理论研究两种情况下（左端梁固支和简支）的斜拉梁的前5阶频率.图3给出了第一种情况（左端梁固支）的前5阶振型.可以发现，两种方法所得的结果几乎完全吻合.因此，表1和图3不仅可以说明本文理论的正确性，还为下面的CFRP索斜拉梁面内自由振动的研究作了铺垫.考虑到工程实际中第一种情况（梁左端固支）的斜拉梁更常见，下面的研究只考虑梁左端固支情况的斜拉梁.

图4给出了不同索力和拉索倾斜角度对CFRP索斜拉梁面内自由振动的各阶频率的影响.一阶频率几乎不随索力大小而改变，倾角的变化有一定的影响，各高阶频率随索力的增大而增大，随拉索倾斜角度的增大而减小，变化较明显.斜拉梁一阶频率对索力和拉索倾斜角度的变化不敏感，原因主要为斜拉梁结构的第一阶振动以梁的振动为主，而索的振动主要是由梁的振动拖动产生.这时，索对于悬臂梁相当于起一个弹性支承的作用，弹性支承主要由索的轴向刚度和倾斜角度决定，索力的改变对弹性支承的影响相对较小.对于2，3，4和5阶的振动，可从振型看出，除二阶振型为索与梁的联合振动外，主要为索的振动，索力和拉索倾斜角度变化时，索的参数发生变化，直接影响到索的振动，因此这几阶频率变化较明显.当索力增大时，斜拉梁整个系统刚度增大，而拉索倾斜角度增加时，拉索变长，其质量也跟着增大，刚度却减小，根据等效频率公式ωeq=keqmeq，频率也就相应地增大和减小了.另外，仔细观察会发现所有相邻两阶频率随索力和拉索倾角的变化发生靠近而又分离的现象，并非两个频率变化曲线交叉，而是两条频率变化曲线转向了（Veering现象），这时两阶振型会发生快速且连续的交换[17]，并且系统两个模态之间发生能量传递，很容易发生内共振现象，这对指导斜拉梁设计，特别是其振动控制具有重要参考价值.

图5给出了斜拉索在不同索力、材料和弹性模量下对斜拉梁结构一阶频率的影响.其中，Ecc中下标第二个c表示CFRP索， Ecg中下标g表示钢索.从中可发现，当采用CFRP索时，索力对一阶频率的影响微乎其微；当采用钢索且索力小于0.5 MN时，一阶频率随索力的增大而增大，当索力大于0.5 MN时，CFRP索和钢索斜拉梁的一阶频率随索力变化的曲线几乎是重合的.这是因为CFRP索斜拉梁不论是大索力下还是小索力下其一阶振型均如图3（a）所示，这样一种模态是梁拖动索振动的模态，所以随着索力的增加其频率基本不变.当采用钢索时，由于其质量要比CFRP索质量大，受其影响振型随索力的变化如图6所示.可看到一阶振型的变化过程是由索振动为主到索梁整体振动再到梁振动为主.因此其一阶频率变化曲线是先增大后持平的变化过程.另外，CFRP索斜拉梁一阶频率随拉索弹性模量的增大而增大，说明可以通过提高拉索弹性模量来提高斜拉梁整体结构的刚度，这是因为4种弹性模量下斜拉梁的振型均如图3（a）所示，此时斜拉梁可以看成是一端固支一端弹簧支撑的梁模型，其振动频率与弹簧刚度有关，弹簧刚度越大，振动频率越大，反之越小.

图7反映了斜拉索在不同材料、索力和弹性模量下对斜拉梁结构二阶和三阶频率的影响.可以看出CFRP索斜拉梁的4条曲线均有一个上升段，之后持平，持平段曲线特征与图6类似.因此我们猜测，上升段的振型是渐变的过程，当到达持平段后，振型基本不再变化.为了验证我们的猜测，我们提取出弹性模量为210 GPa的CFRP索斜拉梁索力在0.3 MN，0.6 MN和1 MN的二阶振型和索力在1 MN，5 MN和10 MN的三阶模态如图8所示.从图8可看出随着索力的增加，第二、三阶振型均是从拉索振动为主到斜拉梁整体振动再到梁振动为主的变化过程，证明我们的猜测是正确的.另外，可以发现使用钢索的斜拉梁要相比于使用CFRP索的斜拉梁随着索力的增加较慢进入持平状态，说明振动阶数越高，拉索质量对其影响越明显.

索力/MN

综合分析图6和图8，可发现索力对斜拉梁结构的动力学特性的影响，主要体现在索与梁刚度相对变化.当索力较小时，拉索振动明显，随着索力的增大，索振动慢慢地弱化，最后变为随梁振动的“摆动”.这是因为索力增大使拉索的横向刚度显著增大（应力刚化），最后拉索所表现出的性质就类似于刚度很大的弹簧.

4结论

本文建立了不考虑垂度影响的CFRP索斜拉梁面内自由振动的力学模型，利用简单的张紧弦和欧拉梁振动理论，采用分离变量法得到它们的振型函数，通过考虑索梁连接处的动态平衡条件，将索和梁的振动耦合到一起，利用传递矩阵法得到斜拉梁面内自由振动的各阶频率方程，从而求得各阶频率值.最后讨论了斜拉梁面内自由振动在不同索力、拉索倾角和拉索材料的变化情况.这种研究方法不仅将复杂的问题简单化，而且能反映实际工程中斜拉梁应有的振动特性，并由此得到以下结论：

1） CFRP斜拉梁结构的面内第一阶自振频率几乎不受索力变化的影响，但随着拉索倾角的改变有不同程度的变化，而钢索斜拉梁第一阶频率则随索力和倾角变化较大.这说明CFRP索斜拉梁的刚度相对稳定.

2） 斜拉梁结构的面内二阶以上振动模态表现出受索力和倾角变化的敏感性，都可能出现频率变化曲线转向（veering）现象，因此为了避免内共振对结构产生不利影响，设计或建造斜拉梁时应该避免使用这些可能产生内共振的参数.

3） CFRP索斜拉梁基本动力学性能优于钢索斜拉梁，特别是在较低索力下和高阶频率上尤为突出，并且弹性模量的增大，对结构的一阶频率的影响较大，振动阶数越高，影响越小.由于工程实际中，高阶振动出现的概率要远小于低阶振动，所以高弹性模量的CFRP索在斜拉梁结构中有着更广阔的应用前景.

4） 随着索力的增加，各阶振动的振型均经历从索振动为主到索梁全局振动再到梁振动为主的变化过程，拉索表现出的性质越来越像一根弹簧，这对拉索振动控制具有重要参考意义.

参考文献

[1]ACI 440.4R04Prestressing concrete structures with FRP tendons[S]. Farmington Hills， USA： American Concrete Institute， 2004： 440.4R-10.

[2]梅葵花， 吕志涛. CFRP斜拉索的静力特性分析[J]. 中国公路学报， 2004，17（2）： 43-45.

MEI Kuihua， LV Zhitao. Static characteristic analysis of CFRP cables[J]. China Journal of Highway and Transport， 2004，17（2）： 43-45. （In Chinese）

[3]梅葵花， 吕志涛， 孙胜江. CFRP拉索的非线性参数振动特性[J]. 中国公路学报， 2007， 20（1）： 52-57.

MEI Kuihua， LV Zhitao， SUN Shengjiang. Property of nonlinear parametric vibration of CFRP cables[J]. China Journal of Highway and Transport， 2007， 20（1）： 52-57. （In Chinese）

[4]康厚军，赵跃宇，朱志辉，等. 强迫激励下CFRP斜拉索面内分叉特性[J]. 湖南大学学报：自然科学版， 2014， 41（9）： 8-13.

KANG Houjun， ZHAO Yueyu， ZHU Zhihui， et al. Inplane bifurcation behavior of inclined CFRP cables subject to external excitation[J]. Journal of Hunan University： Natural Sciences， 2014， 41（9）： 8-13. （In Chinese）

[5]吕志涛， 梅葵花.国内首座CFRP索斜拉桥的研究[J]. 土木工程学报， 2007，40（1）： 54-59.

LV Zhitao， MEI Kuihua. First application of CFRP cables for a cablestayed bridge in China[J]. China Civil Engineering Journal， 2007， 40（1）： 54-59. （In Chinese）

[6]GRACE N F， NAVARRE F C， NACEY R B， et al. Designconstruction of bridge street bridgefirst CFRP bridge in the United States[J]. PCI Journal， 2002，47（5）： 20-35.

[7]FUNG R， LU L， HUANG S. Dynamic modeling and vibration analysis of a flexible cablestayed beam structure[J]. Journal of Sound and Vibration， 2002， 254（4）： 717-726.

[8]GATTULLI V， MORANDINI M， PAOLONE A. A parametric analytical model for nonlinear dynamics in cablestayed beam[J].Earthquake Engineering and Structural Dynamics， 2002， 31：1281-1300.

[9]GATTULLI V， LEPIDI M. Nonlinear interactions in the planar dynamics of cablestayed beam[J]. International Journal of Solids and Structures， 2003， 40： 4729-4748.

[10]赵跃宇， 蒋丽忠， 王连华， 等. 索梁组合结构的动力学建模理论及其内共振分析[J]. 土木工程学报， 2004， 37（3）： 69-72.

ZHAO Yueyu， JIANG Lizhong， WANG Lianhua， et al. The dynamical modelling theory and internal resonance of cablebeam composite structure[J].China Civil Engineering Journal， 2004， 37（3）： 69-72. （In Chinese）

[11]WANG Z Q， YI Z P， LUO Y S. Modeling and nonlinear modal characteristics of the cablestayed beam[J].European Journal of Mechanics A/Solids， 2014， 47： 58-69.

[12]ZHAO Y Y， KANG H J. Inplane free vibration analysis of cablearch structure[J]. Journal of Sound and Vibration， 2008， 172： 363-379.

[13]KANG H J， ZHAO Y Y， ZHU H P. Outofplane free vibration analysis of a cablearch structure[J]. Journal of Sound and Vibration， 2013， 332： 907-921.

[14]WANG Z Q， KANG H J， SUN C S， et al. Modeling and parameter analysis of inplane dynamics of a suspension bridge with transfer matrix method[J]. Acta Mechanic， 2014， 225： 3423-3435.

[15]CLOUGH R W， PENZIEN J. Dynamics of structures[M]. 3rd ed. Berkeley，CA，USA： Computers & Structures， Inc， 2003： 368-375.

[16]HUANH Z H， JONES N P. Damping of tautcable systems： effects of linear elastic spring support[J]. Journal of Engineering Mechanics， 2011， 137： 512-518.

[17]GATTULLI V， LEPIDI M. Localization and veering in the dynamics of cablestayed bridges[J]. Computers & Structures， 2007， 85： 1661-1678.