基于頻率的索力測量方法研究

趙震寅

(浙江理工大學(xué)建筑工程學(xué)院，浙江 杭州 310018)

拉索索力測量方法主要包括：千斤頂油壓表讀數(shù)法、壓力傳感器測試法、磁通量測試法及振動頻率測試法。

在各類測量方法中，千斤頂油壓表讀數(shù)法和壓力傳感器測試法僅可用于施工階段拉索張拉時(shí)的監(jiān)測，無法用于成橋階段的索力檢測。而相較于前述測量方法，頻率法能夠靈活運(yùn)用于各階段的橋梁索力檢測，在操作便捷的同時(shí)也保證了精度，擁有較為完備的測量體系，目前絕大多數(shù)工程都傾向于采用頻率法對橋梁進(jìn)行索力檢測。

但在實(shí)際測量索力的過程中，頻率法使用的因素有多方面影響，包括拉索的抗彎剛度、邊界條件、計(jì)算長度、垂度以及測量溫度，頻率法發(fā)展至今誕生了諸多基于影響因素的索力計(jì)算理論和修正方法。鑒于上述情況，本文梳理了現(xiàn)有頻率法測量索力的研究進(jìn)展，著重論述了各研究成果的優(yōu)化原理，指明了各自的優(yōu)點(diǎn)和局限性，歸納了不同方法的適用條件，并展望了頻率法在今后的研究趨勢。

1 頻率法索力測量原理

處于工作狀態(tài)的拉索會在諸如環(huán)境因素或人工激勵(lì)作用下發(fā)生微小振動，通過使用高分辨率的頻率傳感器以及對應(yīng)的數(shù)據(jù)采集設(shè)備，可獲得拉索振動時(shí)的若干階自振頻率。再結(jié)合索力與拉索線密度、索長、自振頻率、邊界條件、抗彎剛度等的關(guān)系式得出拉索索力[1]。拉索索力計(jì)算公式基于以下假設(shè)：

(1)拉索垂跨比很小，不考慮垂度的影響；

(2)拉索只在豎向平面內(nèi)振動，忽略水平方向的移動，且豎直方向的運(yùn)動也是微小的；

(3)拉索的幾何形狀近似用拋物線表示[2]。

如圖1所示，取拉索任意一段微小長度作為微元體，其平衡方程如下：

圖1 拉索微元體受力圖

(1)

(2)

式中，Q—剪力，N；M—彎矩，N·mm—拉力，N；u—y方向上的位移，mm。

根據(jù)梁彎矩和中和軸彎曲的曲率關(guān)系整理公式(1)和(2)可得到振動微分方程：

(3)

方程的解為：

u(x,t)=u(x)·y(t)

(4)

通過采取離散變量法最終可得到如下關(guān)系式：

y(t)=A1sinωt+A2cosωt

(5)

u(x)=Asinαx+Bcosαx+Csinhβx+Dcoshβx

(6)

式中，y(t)—時(shí)間函數(shù)；u(x)—振型函數(shù)；A、B、C、D—根據(jù)邊界條件確定的待定系數(shù)；ω—拉索的自振圓頻率，rad/s。其中：

(7)

式中，m—拉索單位長度質(zhì)量，kg/m；T—拉索索力，N；EI—拉索的抗彎剛度，N·mm2。

上述推導(dǎo)過程采用了兩端鉸接的拉索模型，為滿足一般邊界條件的設(shè)定，在現(xiàn)有鉸接的基礎(chǔ)上，再在兩端支座處各添加一個(gè)彈簧約束，以此作為最終的力學(xué)模型。假定兩端彈簧鉸的約束剛度相同，k為彈性嵌固系數(shù)。此時(shí)拉索的邊界條件為：

(8)

將(8)代入方程(6)中，通過行列式運(yùn)算和簡化最終得到了一般邊界條件下的拉索索力關(guān)于自振頻率的表達(dá)式：

k2[2αβ(1-cosαLcoshβL)+(α2-β2)
sinαLsinhβL]+2kEI(α3cosαLsinhβL-
α2βsinαLcoshβL+αβ2cosαLsinhβL-β3sinαLcoshβL)
-(EL)2(α2+β2)sinαLsinβL=0

(9)

由于上式的求解較為困難，在實(shí)測中缺乏便捷性，因此理論上只考慮簡支邊界的情況，并逐步發(fā)展成了最為通用的拉索索力計(jì)算理論，即張拉弦理論。

張拉弦理論在支座處僅考慮鉸接情況，拉索不受轉(zhuǎn)動約束，此時(shí)k=0，代入方程(9)中，并結(jié)合式(3)可得到索力的表達(dá)式：

(10)

式中，L—拉索的計(jì)算長度，m；fn—拉索第n階自振頻率，Hz；n—振動階數(shù)。

當(dāng)忽略抗彎剛度時(shí)，公式則轉(zhuǎn)變?yōu)椋?/p>

(11)

式(10)最先由Irvin[3]提出，是基于簡支理論推導(dǎo)的拉索索力計(jì)算公式，考慮抗彎剛度EI的影響，多用于短索索力計(jì)算，也是最基礎(chǔ)的索力計(jì)算模型；式(11)是實(shí)際運(yùn)用最廣泛的經(jīng)典索力計(jì)算公式，因其不考慮抗彎剛度，該狀態(tài)下的拉索為柔性，所以往往用于細(xì)長索的索力計(jì)算，結(jié)果也較為準(zhǔn)確，計(jì)算過程多采用拉索的低階振頻。

從推導(dǎo)過程能夠看出，經(jīng)典索力計(jì)算公式是建立在簡支邊界條件下的，但是實(shí)際拉索的邊界條件多處于簡支和固支之間，亦或是完全固支，直接將該公式運(yùn)用于所有類型拉索的索力計(jì)算中難免會造成誤差，另外計(jì)算長度L的取值準(zhǔn)確程度也會影響到計(jì)算結(jié)果，拉索自身的抗彎剛度、垂度以及所處環(huán)境的溫度都會直接影響到索的振動頻率fn，進(jìn)而影響索力，因此有多項(xiàng)研究探討了上述因素在頻率法中的影響。

2 應(yīng)用中的修正

2.1 邊界條件影響

在索體系的橋梁當(dāng)中，拉索的連接方式主要有銷栓式和承壓式，兩類連接方式的共同點(diǎn)為各自端部至少有一側(cè)可等效為固支，因此必須考慮邊界條件對索力的影響，固支梁理論便是對這一構(gòu)造特性的論證。

Zui[2]提出的固支梁理論就是在方程(9)中考慮k=∞的情況，此時(shí)拉索的轉(zhuǎn)動約束被視作無窮大，可將邊界條件視作固支。

為方便計(jì)算，在此過程中引入了參數(shù)ξ和C。

(12)

(13)

ζ為無量綱參數(shù)，也稱作拉彎比，用于表示拉索的動力特性與張緊弦的接近程度，取值越大接近程度越高，抗彎剛度影響也越小。

根據(jù)ξ取值范圍的不同和所測得的拉索不同階振動頻率，Zui等[2]通過曲線擬合的方法得出了拉索索力-振頻的關(guān)系式：

(1)采用一階振型頻率計(jì)算時(shí)：

(14)

(2)采用二階振型頻率計(jì)算時(shí)：

(15)

(3)采用高階(三階及以上)振型頻率計(jì)算時(shí)：

(16)

上述固支梁理論為后續(xù)索力公式的改進(jìn)提供了依據(jù)，國內(nèi)學(xué)者也提出了眾多基于固支梁理論下考慮拉索邊界條件的索力計(jì)算公式：

任偉新等[4]采用能量法和曲線擬合，建立了分別考慮索垂度和抗彎剛度影響由基頻計(jì)算索力的公式：

(17)

陳淮等[5]采用數(shù)值計(jì)算和曲線擬合，推導(dǎo)出了考慮吊索抗彎剛度和吊索兩端邊界條件影響的索力計(jì)算公式：

(1)采用一階振型頻率計(jì)算時(shí)：

(18)

(2)采用三階振型頻率計(jì)算時(shí)：

(19)

安振源等[6]使用迭代法和曲線擬合方法，建立了索力與一階振動頻率的關(guān)系式：

(20)

以上代表性理論皆是基于兩端固支的邊界條件考慮的，一定程度上彌補(bǔ)了早期國內(nèi)索力計(jì)算領(lǐng)域的空白，雖然顧及了邊界條件但缺陷也是明顯的，例如邊界條件的定義存仍在理想化情況，采用高階頻率計(jì)算時(shí)精度難以把控，公式運(yùn)用也并非十分方便。

為了提高頻率法公式的便捷性，甘泉[7]采用了將固支歐拉梁的振型函數(shù)作為兩端固接的振型函數(shù)：

φn(x)=[cosh(ax)-cos(ax)]-
En[sinh(ax)-sin(ax)]

(21)

(22)

根據(jù)上述條件研究建立了一套較為便捷的可用于固支和簡支邊界的拉索索力計(jì)算公式，相較于經(jīng)典固支梁理論的振型函數(shù)，采用歐拉梁振型函數(shù)可降低公式的復(fù)雜程度，但所基于的邊界條件仍為固支。

對于實(shí)際服役的橋梁而言，拉索邊界條件多介于鉸支和固支之間，拉索端部或多或少會受到抗彎剛度影響，而技術(shù)人員在計(jì)算拉索索力時(shí)習(xí)慣于考慮單一的邊界狀態(tài)，由此必將造成檢測結(jié)果的偏差。

劉文峰等[8]從拉索實(shí)際的邊界條件出發(fā)，提出了固支和簡支耦合形式的振型函數(shù)：

(23)

(24)

式中，r—振型的階數(shù)；l—拉索計(jì)算長度，m；k—固支振動和簡支振動的一個(gè)比值，其大小代表了固支在邊界中所占的比重，可在后續(xù)求解索力的過程中求得。研究最終基于耦合邊界形式的振型函數(shù)給出了拉索索力-剛度關(guān)系式，方法較好地考慮了拉索實(shí)際的端部情況，精度優(yōu)良，但計(jì)算過程仍略顯復(fù)雜。

孫永明等[9]在指出工程實(shí)際拉索與理想勻質(zhì)拉索模型區(qū)別的基礎(chǔ)上，研究了端部性質(zhì)對頻率法的影響，并最終建立了可計(jì)入端部性質(zhì)影響的修正索力計(jì)算公式：

(25)

2.2 抗彎剛度影響

抗彎剛度對索力的影響主要體現(xiàn)在索長較短或索徑較大的拉索上。從構(gòu)造角度出發(fā)，拉索是由數(shù)根平行的鋼筋、鋼絲或鋼絞線經(jīng)小角度扭絞形成。若拉索的鋼絲之間完全滑移，則拉索慣性矩就是各根鋼絲的慣性矩之和；若拉索的鋼絲之間完全粘結(jié)，即各根鋼絲形成了一個(gè)整體，則索的慣性矩為等截面鋼棒的慣性矩，兩種假設(shè)的結(jié)果相差較大。而實(shí)際工程中拉索內(nèi)部構(gòu)造往往介于彼此分離和完全粘結(jié)之間，嚴(yán)格意義上來說拉索的抗彎剛度是一個(gè)未知量。

蘇成等[10]采用有限元法及樣條擬合技術(shù)獲取了拉索索力與頻率的對應(yīng)關(guān)系，提出一種基于多階頻率測試結(jié)果的斜拉索抗彎剛度識別方法，通過實(shí)測數(shù)據(jù)分析得出了拉索的實(shí)際抗彎剛度為0.37(EI)max,其中(EI)max為索按全截面計(jì)算時(shí)的抗彎剛度。

但是實(shí)際拉索的抗彎剛度并不是一個(gè)恒定的值，會因其他參數(shù)的改變而發(fā)生細(xì)微的變化，例如索力和溫度，因此在探尋抗彎剛度對索力影響的過程中，不少學(xué)者另辟蹊徑，提出了一種避開抗彎剛度求算索力的思路，如孟少平等[11]運(yùn)用能量方法建立了索力、抗彎剛度和固有圓頻率之間的關(guān)系：

(26)

式中，φ(x)—振型函數(shù)，通過建立一端鉸支一端固支的振型函數(shù)，配合聯(lián)立1階和2階固有圓頻率的表達(dá)式ω1和ω2，可消去抗彎剛度EI，并最終得到索力關(guān)于前兩階固有振頻的計(jì)算式：

(27)

上式將抗彎剛度作為隱式參數(shù)計(jì)算，回避了工程中拉索抗彎剛度難以識別的問題，在無法精確獲取抗彎剛度的情況下優(yōu)勢顯著；同樣，袁博等[12]也利用能量法，提出了在拉索抗彎剛度未知的情況下，用于計(jì)算拉索索力的計(jì)算方法；Seyed等[13]在推導(dǎo)考慮剛度影響的索力計(jì)算公式時(shí)也選擇了規(guī)避抗彎剛度的方法，但采用的原理并非能量法，首先計(jì)算拉索實(shí)測頻率和計(jì)算頻率間的差值，將該差值定義為關(guān)于剛度EI和索力H的誤差函數(shù)E，再通過E求得一階固有振動頻率后，即可計(jì)算出索力。

當(dāng)不考慮規(guī)避抗彎剛度時(shí)，Chen等[14]通過對一系列纜索進(jìn)行抗彎剛度實(shí)驗(yàn)后得出結(jié)論，纜索的抗彎剛度會隨著索兩端張拉力的增大而增大；而張永平等[15]進(jìn)一步提出纜索抗彎剛度與拉力的非線性關(guān)系模型：

ω2=n2π2T/mL2+n4π4EI/mL2

(28)

式中，ω—纜索的振動頻率，Hz；n—振頻階數(shù)。上式可用于頻率法中拉索抗彎剛度的修正，但仍需針對不同的邊界條件進(jìn)行更深入的研究，完善抗彎剛度的修正模型。

燕啟清等[16]建立了一種考慮短吊桿抗彎剛度的索力計(jì)算公式：

(29)

(30)

(31)

式中，f1—吊桿一階振動頻率，Hz；TQ—標(biāo)定吊桿索力比例系數(shù)時(shí)的千斤頂張拉力，N；Tb—第一次張拉階段吊桿張拉力設(shè)計(jì)值，N；Ts—成橋后吊桿索力設(shè)計(jì)值，N；λ—吊桿邊界條件修正系數(shù)；C—吊桿橫向抗彎剛度修正值；K—吊桿索力比例系數(shù)；η—吊桿索力偏差系數(shù)。

η=1-Tb(Ts+C)/Ts(Tb+C)

(32)

當(dāng)|η|≤Δ時(shí)，可不考慮吊桿橫向抗彎剛度對索力的影響；當(dāng)|η|≥Δ時(shí)，必須用式(32)進(jìn)行索力修正，此時(shí)將被修正的吊桿定義為短吊桿，橫向抗彎剛度對索力的影響不可忽略，Δ的值則根據(jù)相應(yīng)的技術(shù)規(guī)范取定。

2.3 計(jì)算長度

有關(guān)拉索計(jì)算長度的研究本質(zhì)上是由邊界條件引申而出的一類子問題，由于邊界剛度的存在，在計(jì)算拉索索力時(shí)無法取整段索長作為實(shí)際的計(jì)算長度，取值的關(guān)鍵在于精確判別端部抗彎剛度限制拉索有效振動的范圍。

拉索通常是由：套筒、錨杯、中部自由振動的索體3部分組成，如圖2所示。錨固區(qū)域具備較大的剛度，能夠有效限制振動，但是其限制自由部分振動的能力不顯著，在索力計(jì)算時(shí)索長可取拉索能振動部分的長度，設(shè)計(jì)文件中一般為拉索兩端錨墊板之間的距離[17]。陳陽[18]通過采用牛頓切線法迭代計(jì)算有效長度的方法，得出吊索計(jì)算長度可取為上下錨固點(diǎn)之間的理論長度，驗(yàn)證了設(shè)計(jì)規(guī)范中的內(nèi)容，但是其實(shí)驗(yàn)對象為9～38m的吊桿，且實(shí)測數(shù)據(jù)表明索力誤差可達(dá)8%，精度尚待提高。而且在長度更大或是具有一定垂度的拉索有效長度取值中該方法或許并不適用。

圖2 拉索構(gòu)造示意圖

經(jīng)有關(guān)學(xué)者實(shí)驗(yàn)分析得出，拉索有效長度取值的精確程度對長索的影響較小，但對短索的影響較大，有數(shù)據(jù)表明誤差最高可達(dá)66%[19]，這就意味著在面對短索索力的計(jì)算時(shí)要更為嚴(yán)格地取定拉索有效振動長度。

宋一凡和賀拴海[20]引入了動力計(jì)算長度的概念，解決了用振動法測估索力時(shí)面臨的有效索長取值難題，提出了基于二階固有頻率的動力計(jì)算長度lde的實(shí)用公式：

(33)

為減緩拉索振動，保證橋梁整體的穩(wěn)定性，在實(shí)際布置拉索時(shí)一般會在套筒內(nèi)安裝減振片等阻尼裝置，或是在外部安裝減振器。減振器的設(shè)置會縮短拉索的自由振動長度因此需重新定義有效長度，胡利平和林育洪[21]采用的方法是分別測取拉索安裝減振器前后的頻率，再利用2次頻率值的不同來換算索長：

(34)

式中，l′—換算后的索長，m；l—實(shí)際索長，m；fn、f′n—安裝減振器前后的某一階拉索振動頻率，Hz，一般可取基頻。假如事先沒有進(jìn)行該項(xiàng)操作，則可先取拉索上下減振器之間的長度作為計(jì)算長度再適當(dāng)估計(jì)結(jié)果誤差范圍進(jìn)行調(diào)整的辦法。相較于其他理論，采用此方法的分析結(jié)果較快捷，具備一定的工程推廣性，但并沒有針對不同邊界條件的拉索進(jìn)行公式的細(xì)化，在使用過程中也難免會產(chǎn)生誤差；林天然[22]在采用張拉弦理論計(jì)算索力的基礎(chǔ)上提出用吊桿計(jì)算長度的修正值對索力進(jìn)行修正，通過對某一大跨度拱橋吊桿進(jìn)行索力測試，最終確定修正值為減振器到錨固端距離的0.9～1.1倍之間，該法同樣是一種快捷的計(jì)算長度修正手段，但在特殊情況下也同樣面臨著精度不足的問題。

Li[23]根據(jù)能量法發(fā)展了一種基于附加質(zhì)量塊的有效長度確定方法，得到了固有頻率和有效長度間存在的解析關(guān)系：

(35)

式中，ω1、ωm1—施加質(zhì)量塊前后拉索的一階固有振頻，Hz；m—拉索的單位長度質(zhì)量，kg；M—質(zhì)量塊的質(zhì)量，kg。Li將該公式運(yùn)用于某吊桿拱橋的索力檢測中，通過與實(shí)測值進(jìn)行比對顯示此方法可將索力誤差控制在5%，在短索中的預(yù)測精度表現(xiàn)較為出色。

當(dāng)面對具有一定垂度的長拉索時(shí)，基于拋物線理論的長度計(jì)算法是設(shè)計(jì)中普遍采用并且認(rèn)可度較高的方法，但隨著拉索長度的增加，垂度的非線性影響會不斷增大，傳統(tǒng)拋物線理論的計(jì)算精度會因此受到影響，考慮到這一問題，孟慶成等[24]提出采用懸鏈線公式來計(jì)算斜拉索的無應(yīng)力長度，長度的計(jì)算公式為：

l0=l-Δl

(36)

(37)

(38)

式中，l0—拉索無應(yīng)力長度，m；H1—拉索下錨點(diǎn)拉力的水平分量，N；gcb—拉索的單位長度重量，kg/m；a—上下錨點(diǎn)之間的水平距離，m；λ1—下錨點(diǎn)處拉索切線與水平的夾角。實(shí)測結(jié)果顯示，該公式下的各階頻率估算的索力波動較小，改進(jìn)了用拋物線法近似求解索長的傳統(tǒng)方法，可用于大跨度橋梁的拉索計(jì)算中。

除上述理論外，學(xué)者們在研究過的程中也建立了眾多從不同原理出發(fā)的，關(guān)于計(jì)算拉索有效長度的實(shí)用公式。周正茂[25]提出了振型節(jié)點(diǎn)法，利用3階以上模態(tài)振型的中間節(jié)點(diǎn)間距作為計(jì)算長度，該方法的優(yōu)點(diǎn)在于回避了復(fù)雜的邊界條件，原理簡潔思路新穎，可直接運(yùn)用于公式(10)中，以三階模態(tài)為例，選取三階模態(tài)頻率f3作為計(jì)算頻率，取中間相鄰2個(gè)模態(tài)節(jié)點(diǎn)間的距離作為計(jì)算長度L，因?yàn)榇藭r(shí)只取了一段距離，所以階數(shù)n取為1，在更高階模態(tài)中n可取2或更大，L也要進(jìn)行相應(yīng)調(diào)整。方法的不足之處在于其精度易受索長、振型階次和邊界條件的影響，仍有很大的改進(jìn)空間。

史杰[26]通過建立靜載試驗(yàn)前后吊桿索力差和頻率差的解析式，反推得出了索力計(jì)算長度的表達(dá)式：

(39)

其中：

(40)

(41)

式中，n—振動階數(shù)；ρ—拉索線密度，kg/m；δT—加載前后索力的增量，N；fn,1、fn,2—加載前后對應(yīng)的n階振動頻率，Hz。由公式(41)可知索力增量和頻率平方差成正比，因此Kn可通過拉索加載前后的索力增量和頻率數(shù)據(jù)擬合得到。此方法原理簡單，操作步驟簡潔，可運(yùn)用于新舊橋梁的索力檢測。

2.4 垂度影響

實(shí)際許多的橋梁拉索會呈現(xiàn)出下垂的狀態(tài)，呈現(xiàn)出很強(qiáng)的非線性性質(zhì)[27]。拉索垂度的成因包括：索力、索長、自重、傾角等，垂度的存在對拉索的最直接影響就是其振動頻率。在拉索長度確定的前提下，拉索垂度對基頻的影響較大，但是對高階頻率(四階及以上)的影響較小。

早期的研究例如島田忠幸[28]通過引入?yún)?shù)Γ來量化拉索垂度對固有頻率的影響：

(42)

Irvine[3]通過推導(dǎo)微元體平衡方程，提出了用無量綱參數(shù)λ2考慮拉索垂度對索力的影響，該參數(shù)也成為了拉索弦振動理論中的基本參數(shù)之一：

(43)

(44)

參數(shù)λ2反映的是索的垂度和彈性模量對基頻的影響，具體表現(xiàn)為當(dāng)λ2較大時(shí)，索力較小，拉索的垂度較大，對基頻的影響也較大，反之則越小。

由于直接采用λ2求解無法解得特定垂度下索力與振動基頻的關(guān)系式，因此往往用能量法近似來實(shí)現(xiàn)這一目標(biāo)，得到的考慮垂度的索力-基頻的關(guān)系式如下：

(45)

式中，ω—拉索振動的基頻，Hz；H—索力的弦向分量，N；l—拉索計(jì)算長度，m；m—拉索單位長度質(zhì)量，kg/m。最終可以得出考慮索垂度和彈性影響的索力計(jì)算實(shí)用公式[29]：

(46)

上述公式在行業(yè)方面的認(rèn)可度高且較為基礎(chǔ)，但在拉索振型的考慮上存在局限，只是給出了基于基頻的索力公式，并未考慮拉索垂度對更高階振動頻率的影響。

李國強(qiáng)等[30]通過驗(yàn)證各因素對索力檢測的影響得出結(jié)論：垂度只影響索的奇數(shù)階振頻，尤其反映在一階頻率上，對高次奇數(shù)階頻率的影響較小，對偶數(shù)階頻率的影響可忽略不計(jì)。在涉及垂度作用下的拉索頻率選取中，該結(jié)論有著重要的參考價(jià)值。

同樣，在處理垂度對不同階頻率的影響時(shí)，冉志紅和李喬[31]通過采用奇異攝動解法這一種求解微分方程近似解析解的方法來求解拉索的動力平衡方程，給出了考慮拉索垂度影響的索力計(jì)算公式：

(47)

其中：

(48)

劉志軍[32]在研究用垂度修正索力的過程中，提出了采用傳遞矩陣法考慮垂度對拉索索力的影響，該方法基于拉索的微元體靜力平衡方程，將拉索分為了n+2個(gè)質(zhì)點(diǎn)和n+1個(gè)沒有質(zhì)量只有張力的弦，建立了求解索力的矩陣方程，方程的關(guān)鍵點(diǎn)在于修正垂度影響下的拉索自振頻率，進(jìn)而使得索力預(yù)測的結(jié)果更加準(zhǔn)確。傳遞矩陣法易于編程且精度可靠，在實(shí)際工程中具備快速測量的優(yōu)勢。

相較于用垂度修正索力，運(yùn)用拉索垂度特性反求索力也是一條新思路，葛俊穎等[33]通過有限元軟件分析證明了當(dāng)斜拉索的索力較大時(shí)，其拉索中部位置的索力和垂度存在著確定的函數(shù)關(guān)系，且不受拉索本身剛度和兩端支承條件的影響。當(dāng)選取一定長度的中間索段時(shí)索力的計(jì)算公式為：

(49)

式中，TM—該索段的平均索力，N；T0—最小索力，N；TL—最大索力，N，分別位于該段的最上端和最下端；L—索段的索長，m；α—對應(yīng)的傾角；fmax—最大垂度，后3項(xiàng)參數(shù)可由儀器測得，最后根據(jù)該段拉索的各項(xiàng)索力可確定整條拉索的索力情況，選取的索長越大結(jié)果越精確，且不受抗彎剛度和邊界條件的影響。

2.5 測量溫度影響

從拉索自身溫度效應(yīng)出發(fā)，拉索長度會隨材料的熱脹冷縮而產(chǎn)生細(xì)微的改變，從而造成索力和頻率的變化；整體上，橋梁拉索通常由鋼材制作而成，與之錨固的主梁則大多為混凝土構(gòu)件，材料間不同的膨脹系數(shù)決定了二者在同溫度變化下的變形存在差異，變形速率也不同，整體結(jié)構(gòu)的部分變形不一致也導(dǎo)致了索力隨溫度的變化，而實(shí)際的情況更為復(fù)雜，橋梁的溫度日變化幅度呈現(xiàn)出：氣溫>拉索溫度>梁、主塔溫度的梯度狀況，這使得索力關(guān)于溫度的變化值變得更難以捕捉。

一座橋梁在正常的大氣環(huán)境中會經(jīng)歷酷暑和嚴(yán)寒，年溫差通常會超過50℃。溫差會對頻率造成影響進(jìn)而影響索力，同時(shí)也會對垂度造成影響，所以拉索的索力、垂度、振動頻率、環(huán)境溫度之間的關(guān)系是復(fù)雜的，關(guān)鍵就在于能否將溫度的影響準(zhǔn)確量化，因此國外學(xué)者Kim[34]和Park[35]開始嘗試一種能夠結(jié)合更復(fù)雜因素的參數(shù)敏感性方法(PSM)來識別多因素影響下的索力。

國內(nèi)方面，候俊明[36]曾對某一單索面斜拉橋進(jìn)行索力監(jiān)測，實(shí)測結(jié)果顯示拉索索力會隨著溫度而異，在室外時(shí)隔4h的索力變化幅度可達(dá)到2.5%，由此可見溫度對于索力的影響不可忽略。但研究并沒有給出基于溫度變化的索力計(jì)算公式或修正方法，同時(shí)缺少了溫度對其他結(jié)構(gòu)參數(shù)影響的討論。任遠(yuǎn)[37]嘗試用限元模型對拉索進(jìn)行溫度效應(yīng)的監(jiān)測，雖然能夠有效反映索力的變化趨勢，但在施工控制中只能作為一種近似的手段。

在獲取了索力隨溫度變化趨勢的前提下，Chen[38]采用總體經(jīng)驗(yàn)?zāi)B(tài)分解法(EMD)對索力、環(huán)境溫度、索溫變化量的時(shí)域圖進(jìn)行時(shí)程分析，研究了基于環(huán)境氣溫和拉索索溫變化量的索力修正系統(tǒng)，并在EMD的基礎(chǔ)上更進(jìn)了一種基于噪聲輔助的數(shù)據(jù)分析方法(EEMD)，方法取得了顯著的效果，較大程度地消除了拉索的溫度效應(yīng)，但該技術(shù)的難度較大且局限于低塔斜拉橋的橋型，目前仍位于改進(jìn)階段。

考慮到需要建立一種可供現(xiàn)場快速使用的公式方法，Ma等[39]采用有限差分法來處理拉索動力方程中溫度的影響，并最終求得了索力關(guān)于溫度變化的數(shù)值解，但是公式仍較復(fù)雜；艾玉麒[40]提出了一種基于溫度的索力變化量公式：

ΔT1=KαEA(t1-t)

(50)

式中，K—溫度修正系數(shù)，可由相應(yīng)的有限元模型計(jì)算得到；α—線膨脹系數(shù)；t1—實(shí)時(shí)溫度，℃；t—設(shè)計(jì)溫度，℃。對比文獻(xiàn)[39]的公式，該方法原理簡明且相關(guān)參數(shù)易被采集，在涉及溫度的索力修正中其運(yùn)用更加便捷。

3 總結(jié)

本文通過總結(jié)影響頻率法索力測量的5項(xiàng)因素，得出以下結(jié)論：

(1)定義單一邊界條件的做法較為理想化，與實(shí)際不符。基于耦合邊界條件的剛度識別法能夠求得更貼近真實(shí)的拉索振型函數(shù)，使結(jié)構(gòu)分析更加合理。

(2)抗彎剛度對索力的影響主要針對短、粗拉索，可采用剛度修正公式、剛度系數(shù)折算法對抗彎剛度進(jìn)行求解。當(dāng)剛度難以測取時(shí)，可采用回避剛度的隱式索力計(jì)算法。

(3)計(jì)算長度對索力的影響主要針對端部約束大、長度短的拉索，根據(jù)經(jīng)驗(yàn)公式選取自由振動長度的方法雖快捷但不適用于短索，可采用基于固有振頻的有效振動長度計(jì)算法獲取計(jì)算長度。

(4)垂度的影響對象為長拉索，在面對具有一定垂度的拉索時(shí)應(yīng)盡量選取高階偶數(shù)階振頻計(jì)算索力。后續(xù)關(guān)于垂度的研究將逐步擺脫階數(shù)的限制，使得各階頻率計(jì)算下的索力結(jié)果都趨于穩(wěn)定。

(5)溫度影響的量化難度較大，不同構(gòu)件關(guān)于同一溫度變化所產(chǎn)生的溫度變形各不同。故測量索力時(shí)應(yīng)盡量避免在溫度較高的正午進(jìn)行?，F(xiàn)有頻率法中關(guān)于溫度影響的研究仍不成熟，未來需持續(xù)更進(jìn)。