溫度變化對懸索主共振響應影響分析

趙珧冰， 彭 劍

(1. 華僑大學 土木工程學院， 福建 廈門 361021； 2. 湖南科技大學 土木工程學院， 湖南 湘潭 411201；3. 湖南大學 土木工程學院， 長沙 410082)

溫度變化對懸索主共振響應影響分析

趙珧冰1， 彭 劍2, 3

(1. 華僑大學 土木工程學院， 福建 廈門 361021； 2. 湖南科技大學 土木工程學院， 湖南 湘潭 411201；3. 湖南大學 土木工程學院， 長沙 410082)

基于增量熱場理論，推導考慮均勻溫度變化影響下受諧波激勵的懸索非線性運動微分方程，利用Galerkin法將方程離散，并對離散后的方程進行線性分析。以懸索主共振響應為例，利用多尺度法求解其高階近似解及幅頻響應方程，通過算例研究溫度變化對不同Irvine參數的懸索前四階模態頻率以及幅頻響應曲線的影響。研究結果表明：懸索模態振型、頻率和主共振響應受溫度變化影響明顯，且與其Irvine參數密切相關；溫度變化有可能定性和定量地改變懸索非線性振動特性，其取決于線性項、平方和立方非線性項系數受溫度變化的影響程度；相同程度的升溫和降溫對懸索振動特性的影響不對稱。

懸索； 溫度變化； 多尺度法； 主共振

索結構是大跨度斜拉橋以及懸索橋重要的基本單元，在實際工程中，由于索承受的荷載極其復雜，周圍環境變化劇烈，因此不可避免會遭受損傷或性能退化，屬于易損結構[1]。為了保證橋梁結構運營安全，需要及時、正確地識別索結構的損傷特性，并采取有效的措施保證結構安全可靠[2]。而這些損傷特征信息往往淹沒在各種環境因素的影響中，其中以溫度變化影響較為突出[3]。為此研究索結構在溫度變化時的振動特性就變得非常重要。

無論是理論研究[4-10]還是工程實踐[11-12]，研究人員非常關注溫度效應和損傷效應對頻率、振型和張拉力等振動特征的影響。一方面，研究人員建立了考慮溫度效應影響下拉索靜、動力學分析模型；另一方面，基于拉索損傷程度、位置和范圍三個參數，建立了拉索損傷狀態下分析模型。然而懸索作為一類同時包含二次和三次非線性項的柔性結構，在荷載激勵下會產生豐富的非線性動力學現象。實驗結果表明[13]：溫度變化對懸索非線性振動中的分叉現象有顯著影響。因此基于考慮溫度效應影響下的非線性動力學模型，深入研究溫度變化對懸索非線性振動特性影響，既可以豐富拉索非線性動力學理論研究，還可以為結構損傷識別提供參考和依據。

本文以懸索為例，基于增量熱場理論，建立考慮溫度變化影響的動力學模型，推導非線性運動微分方程。采用多尺度法，求解其非線性主共振響應的高階近似解及幅頻響應方程。通過算例研究溫度變化對不同垂跨比懸索的非線性主共振響應影響。

1 數學模型

圖1 懸索構形及特性

1.1 運動微分方程

基于懸索靜平衡和熱應力平衡狀態時的平衡方程，受外激勵作用下懸索的運動方程可以表示為：

(1)

(2)

其中：cu和cv分別表示懸索在u和v方向阻尼系數。

根據鏈式法則，結合熱應力狀態時的平衡方程，可得[5]：

(3)

(4)

式中：Ξ和Ξd分別表示在熱應力平衡狀態下的靜拉力及在運動狀態下增大的動拉力沿懸索弦向分量，其表達式分別為：

和

由于溫度變化會影響懸索垂度和張拉力，因此引入了兩個無量綱參數κ2和χ2。基于擬靜態假設，懸索面內振動的運動方程為：

(5)

式中：Ω為外激勵的頻率。

為研究簡便，引入以下無量綱參數：

(6)

將式(6)代入式(5)中，略去表達式的上標星號，則相應的無量綱運動方程為：

(7)

式中：y(x)=4fx(1-x)和f=b/L。

1.2 Galerkin離散

假設懸索是由正對稱和反對稱模態組成的多自由度系統，利用Galerkin法對式(7)進行無窮維離散，面內運動的位移函數表示為：

φ(x)

(8)

式中：N為模態離散的數目；qn(t)是廣義坐標函數；φ(x)為n階模態函數。

將式(8)代入式(7)中，可得：

(9)

由此可見溫度變化并沒有改變懸索運動微分方程的形式，而是通過引入的兩個無量綱參數(χ2和κ2)來體現溫度變化的影響。

2 主共振響應分析

對于有阻尼系統，多尺度法處理起來更加方便，而且由于低階的多尺度近似解有可能產生定性的誤差，因此本文采用高階的多尺度法求解懸索非線性主共振響應的近似解[14]。

此時外激勵頻率表示為：

Ωn=ωn+ε2σn1+ε4σn2

(10)

式中：ωn為線性分析得到的模態頻率；ε為無量綱小參數；σn1和σn2為前兩階調諧參數。

利用多尺度法可得系統的幅頻響應方程：

(11)

(12)

而總的調諧參數參數表示為：

σn=σn1+σn2

(13)

因此系統的穩態響應解為：

(14)

3 算例分析

假設懸索彈性模量和阻尼系數均與溫度變化無關，其各項物理參數為：L=200 m、A=7.069×10-2m2、E=200 GPa、ρ=7 800 kg/m3、α=1.2×10-5℃以及g=9.81 m/s2。無量綱化后的激勵幅值以及阻尼系數分別為0.003和0.002 5。

圖2描繪了在整體均勻溫度變化影響下，懸索前四階模態頻率(前兩階正/反對稱模態)與Irvine參數λ2的關系曲線。由圖可知，前四階模態頻率受溫度變化的影響明顯，尤其是第一階正對稱模態頻率[15]。反對稱模態頻率雖與溫度變化呈現反比例關系，但是其模態振型不受溫度變化的影響。一階正/反對稱模態頻率之間存在交點，在交點附近容易發生模態間的內共振以及能量傳遞現象。然而在溫度變化的影響下，交點的位置發生了漂移，表明懸索內共振響應受溫度變化影響明顯。

圖2 溫度變化對Irvine參數與無量綱模態頻率關系曲線影響

Fig.2 Temperature effects on relationship curves between Irvine parameter and non-dimensional mode frequencies

圖3描繪了溫度變化對四種垂跨比懸索主共振幅頻響應曲線的影響。溫度變化和線性頻率之間的正、反比例關系與初始張拉力密切相關[6]。而懸索平方和立方非線性項系數的大小則和溫度變化呈現出明顯的正比例關系，即溫度上升，非線性項系數增大，反之則減小。

對于圖3(a)，此時懸索垂跨比約為0.006，對應Irvine參數約為1.45，從線性結果分析來看，溫度上升，線性頻率下降，溫度下降，線性頻率上升。此時從定性角度而言，懸索均表現出明顯的“硬”彈簧特性(曲線向右偏轉)；從定量角度來看，溫度變化對主共振響應影響比較有限。由此推斷，當懸索垂跨比非常小時，溫度變化對其非線性振動特性影響不明顯。

圖3(b)中，懸索的垂跨比為0.01，Irvine參數約為6.72。線性分析結果表明：無論溫度上升還是下降，線性頻率均會增大。而從非線性分析結果來看，在不考慮溫度變化(Δt=0 ℃)和溫度下降(Δt=-40 ℃)兩種情況下，懸索的幅頻響應曲線均向右偏轉。而當溫度上升時(Δt=+40 ℃)，懸索會表現出先偏“軟”再偏“硬”的振動特性，即曲線先向左然后向右偏轉。這一現象的原因在于當溫度上升時，懸索平方非線性項系數受溫度變化影響最為明顯，此時平方非線性在振動中占主導地位，從而使得曲線向左偏轉。

懸索的垂跨比進一步增大(f=0.015)，對應的Irvine參數也進一步增加(λ2=21.92)，此時溫度變化影響下的幅頻響應曲線如圖3(c)所示。當溫度上升時，懸索線性頻率也隨之增大。而對于主共振響應而言，無論溫度是上升還是下降，懸索均表現出先“軟”后“硬”的振動特性。但是從定量角度而言，當溫度上升時，升溫對立方非線性項系數的影響程度明顯大于降溫，因此隨著溫度上升，曲線向右偏轉程度明顯加劇。

最后懸索垂跨比進一步增大至0.02，對應的Irvine參數為60.04。由圖3(d)所示，溫度變化明顯改變了懸索主共振響應的振動特性。特別是當溫度上升時，此時懸索的“軟”彈簧特性幾乎消失。盡管平方非線性項的系數隨著溫度上升而增加，但是此時立方非線性項的系數隨著溫度的上升增大約為3倍。因此“硬”彈簧特性占據主導地位，此時曲線迅速向右偏轉。

(a)f=0.006

(b)f=0.01

(c)f=0.015

(d)f=0.02

圖3 溫度變化對懸索幅頻響應曲線影響：

Fig.3 Temperature effects on frequency response curves of the suspended cable

4 結 論

平方和立方非線性項的系數與溫度變化呈現出明顯的正比例關系，而線性項系數與溫度變化之間的比例關系則和懸索垂跨比密切相關；懸索非線性主共振響應受溫度變化影響較為明顯，與其垂跨比密切相關，但當垂跨比很小(初始張拉力很大)，溫度變化對懸索振動特性影響明顯減弱；懸索主共振響應特性在溫度變化影響下，有可能發生定性和定量地改變，取決于線性項、平方和立方非線性項受溫度變化的影響程度；相同程度的升溫和降溫對懸索振動特性影響不對稱。

[1] 李兆霞, 王瀅, 吳佰建, 等. 橋梁結構劣化與損傷過程的多尺度分析方法及其應用[J]. 固體力學學報, 2010, 21(6): 731-756.

LI Zhaoxia, WANG Ying, WU Baijian, et al. Multi-scale modeling and analyses on structural deterioration and damage in long-span bridges and its application[J]. Chinese Journal of Solid Mechanics, 2010, 21(6): 731-756.

[2] 李愛群, 丁幼亮, 王浩, 等. 橋梁健康監測海量數據分析與評估—“結構健康監測”研究進展[J]. 中國科學: 科學技術, 2012, 42(8): 972-984.

LI Aiqun, DING Youliang, WANG Hao, et al. Analysis and assessment of bridge health monitoring mass data-progress in research/development of Structural Health Monitoring[J]. Sci China Tech Sci, 2012, 42(8): 972-984.

[3] 陳波, 鄭瑾, 王建平. 橋梁結構溫度效應研究進展[J]. 武漢理工大學學報, 2010, 32(24): 79-83.

CHEN Bo, ZHENG Jin, WANG Jianping. State-of-the-art of the temperature effects of bridges[J]. Journal of Wuhan University of Technology, 2010, 32(24): 79-83.

[4] TREYSSEDE F. Free vibrations of cables under thermal stress[J]. Journal of Sound and Vibration, 2009, 327(1): 1-8.

[5] LEPIDI M, GATTULLI V. Static and dynamic response of elastic suspended cables with thermal effects[J]. International Journal of Solids and Structures, 2012, 49(9): 1103-1116.

[6] 趙珧冰, 孫測世, 彭劍, 等. 溫度變化對拉索頻率與索力的影響[J]. 應用力學學報, 2013, 30(6): 904-908.

ZHAO Yaobing, SUN Ceshi, PENG Jian, et al. The effects of temperature changes on the frequencies and tension forces of cable[J]. Chinese Journal of Applied Mechanics, 2013, 30(6): 904-908.

[7] 汪峰, 陳福青, 文曉旭, 等. 考慮溫度影響的斜拉索參數振動模型及響應分析[J] 重慶交通大學學報(自然科學版), 2016, 35(2): 1-5.

WANG Feng, CHEN Fuqing, WEN Xiaoxu, et al. Analysis of cable parametric vibration model and response with consideration of temperature effect[J]. Jounral of Chongqing Jiaotong University (Natural Science), 2016, 35(2): 1-5.

[8] LEPIDI M, GATTULLI V, VESTRONI F. Static and dynamic response of elastic suspended cables with damage[J]. International Journal of Solids and Structures, 2007, 44: 8194-8212.

[9] 王立彬, 王達, 吳勇. 損傷拉索的等效彈性模量及其參數分析[J]. 計算力學學報, 2015, 32(3): 339-345.

WANG Libin, WANG Da, WU Yong. The equivalent elastic modulus of damaged cables and parameter analysis[J]. Chinese Journal of Computational Mechanics, 2015, 32(3): 339-345.

[10] ZHU J, YE G R, XIANG Y Q, et al. Dynamic behavior of cable-stayed beam with localized damage[J]. Journal of Vibration and Control, 2011, 17(7): 1080-1089.

[11] 樊可清, 倪一清, 高贊明. 大跨度橋梁模態頻率識別中的溫度影響研究[J]. 中國公路學報, 2006, 19(2): 67-73.

FAN Keqing, NI Yiqing, GAO Zanming. Research on temperature Influences in long-span bridge[J]. China Journal of Highway and Transport, 2006, 19(2): 67-73.

[12] 李順龍, 李惠, 歐進萍, 等. 考慮溫度和風速影響下的橋梁結構模態參數分析[J]. 土木工程學報, 2009, 42(4): 100-106.

LI Shunlong, LI Hui, OU Jinping, et al. Identification of modal parameters of bridges considering temperature and wind effects[J]. China Civil Engineering Journal, 2009, 42(4): 100-106.

[13] REGA G, ALAGGIO R. Experimental unfolding of the nonlinear dynamics of a cable-mass suspended system around a divergence-hopf bifurcation[J]. Journal of Sound and Vibration, 2009, 322(3): 581-611.

[14] ZHAO Y B, SUN C S, WANG Z Q, et al. Analytical solutions for resonant response of suspended cables subjected to external excitation[J]. Nonlinear Dynamics, 2014, 78: 1017-1032.

[15] 趙珧冰, 孫測世, 彭劍, 等. 不同初拉力拉索對溫度變化的敏感性分析[J]. 中南大學學報(自然科學版), 2014, 45(5): 1680-1685.

ZHAO Yaobing, SUN Ceshi, PENG Jian, et al. Sensitivity analysis of different initial tension forces of suspended cable to temperature changes[J]. Journal of Central South University(Science and Technology), 2014, 45(5): 1680-1685.

Effects of temperature changes on primary resonances of suspended cables

ZHAO Yaobing1, PENG Jian2,3

(1. College of Civil Engineering, Huaqiao University, Xiamen 361021, China;2. School of Civil Engineering, Hunan University of Science and Technology, Xiangtan 411201, China;3. College of Civil Engineering, Hunan University, Changsha 410082, China)

Based on the incremental thermal field theory, the nonlinear vibration equations of a suspended cable under harmonic excitation considering uniform temperature variation were derived. Galerkin method was used to discretize the nonlinear vibration equations and the discretized equations were analyzed linearly. The higher order approximate solution and the amplitude-frequency response equation were obtained with the multi-scale method. The effects of temperature changes on the first 4 orders modal frequencies and amplitude-frequency response curves of the suspended cable with different Irvine parameters were studied with numerical examples. The results showed that modal shapes, modal frequencies and primary resonance responses of the suspended cable are obviously affected by temperature changes, and they are closely related to Irvine parameters of the suspended cable; the suspended cable’s nonlinear vibration characteristics may be changed quantitatively and qualitatively by temperature effects, they are dependent on the effects of temperature changes on coefficients of the linear term, the square one and the cubic one; the effects of temperature increase and its decrease at the same level on the vibration characteristics of the suspended cable are not symmetric.

suspended cable; temperature variations; multi-scale method; primary resonance

國家自然科學基金(11602089; 11402085); 福建省自然科學基金青年創新(2016J05011); 福建省中青年教師教育科研項目(JAT160025)；華僑大學高層次人才科研啟動費(15BS409)

2015-12-30 修改稿收到日期：2016-06-16

趙珧冰 男，博士，講師，1988年9月生

TH1

A

10.13465/j.cnki.jvs.2017.15.036