視頻圖像測振技術對斜拉索索力的應用

徐旻杰， 袁向榮， 劉 輝， 陳澤賢， 任張晨

(廣州大學 土木工程學院, 廣州 510006)

視頻圖像測振技術對斜拉索索力的應用

徐旻杰， 袁向榮， 劉 輝， 陳澤賢， 任張晨

(廣州大學 土木工程學院, 廣州 510006)

以數碼攝像頭與PC機作為硬件設備,結合Matlab軟件設計視頻圖像測振系統，采用該系統識別了斜拉索模型的1～4階模態參數。采集索振動視頻，將其分解為圖像序列，經圖像處理，得索邊緣序列，取邊緣各點的振動時程，經振動模態分析，得到索的頻率、阻尼和振型，由頻率識別索力。以力傳感器所得索力數據作為參考，同通過視頻圖像測振技術得到的結果進行對比分析。結果表明，視頻圖像測振技術檢測索振動可行且具有一些優勢。隨著光學硬件的進一步發展和測振軟件系統的日益成熟，視頻圖像測振技術的工程應用前景光明。

斜拉索； 邊緣識別； 多項式擬合； 高階模態； 索力測量

0 引 言

在斜拉橋的施工過程中，監控斜拉索的索力對滿足施工過程中施工控制[1]的要求以及成橋后橋梁的結構內力分布起著至關重要的作用。在斜拉橋的使用過程中拉索由于振動、疲勞衰減和腐蝕等原因造成索力松弛[2]，拉索的損害會造成結構不可挽回的后果。索力的變化會影響結構的力學狀態和結構線型[3]，因此斜拉索的索力可以作為斜拉橋的結構健康狀態評估的重要指標。目前可供現場檢測索力的方法主要有：① 壓力表測定法；② 壓力傳感器測定法；③ 頻率法；④ 磁通量法；⑤ 光纖光柵振動測試法。本文探索索力檢測的視頻方法。

圖像測量技術具有表面全尺度、非接觸式、無負載效應、環境適應性強、重復可比性好、無設備耗損等優點，不僅適用于靜態測量，也適用于動態測量，還可進行全域高密度檢測。隨著數碼影像產品的快速進步和普及化、以及實際工程對圖像測量技術的應用與需求，數字測量技術迅速發展。例如有劉敏[4]提出的結構一維大變形識別的數字圖像邊緣檢測；袁向榮[ 5]提出的結構邊緣變形檢測一維DIC法，其邊緣識別精度可達到0.05像素；胡朝輝[6]提出的視頻圖像技術的結構振動測試；袁向榮[7]邊緣識別的正交多項式擬合及梁變形檢測。

本文基于數字視頻技術，以數碼攝像頭與PC機作為硬件設備結合Matlab軟件設計視頻圖像測振系統[8]，采用視頻圖像測振系統檢測索的振動，由圖像處理及模態分析方法識別索的頻率、阻尼及振型，由頻率識別索力。

1 視頻圖像檢測系統

視頻圖像檢測系統是基于圖像檢測技術的一種新的振動測試方法。其原理是通過拍攝設備的感光元件現場記錄被測對象的動態響應數據，后期通過已編程序對動態響應數據的提取和處理，得出所需要的物理信息。相較于傳統的振動測試，視頻圖像檢測法有著非接觸式的優點，因此免去了很多現場布線、傳感器布置的問題，能夠大大地減少現場的工作量，加快檢測速度，節省人力，對現場的測量人員的安全也有一定的提高。

1.1 硬件系統

視頻的拍攝可以采用單反相機、卡片相機，以及手機攝像頭。單反相機的優點是其感光元件大，精度高，適應相對惡劣的現場環境，提高工作范圍，但單反相機需要一定的學習成本，其便攜性也不佳；卡片相機的優點是可以拍攝更高幀數的視頻，目前民用卡片相機已經可以拍攝1 000幀的視頻，學習成本也較低；手機攝像頭的優點是便攜，操作方便，不需要學習成本，其精度也可以滿足大多數工程檢測。

1.2 軟件系統

軟件系統是基于Matlab編制的程序，包括視頻采集、視頻分解、信號分析4個模塊，如圖1所示。

2 動態邊緣的識別方法

經典邊緣檢測算子[10]只能識別整像素邊緣，要提高數據精度則需要采用亞像素邊緣定位方法。目前常用的亞像素邊緣識別方法有高斯曲線擬合法[11]、矩算子法[12]和多項式擬合法[13]。高斯曲線擬合法具有更高的定位精度，能更貼合實際的邊緣模型，但該方法過程繁瑣，計算量大，不易快速得出結果;矩算子法以理想的二級階躍型邊緣計算實際邊緣，因而存在原理上的誤差;多項式擬合法則集合了計算精度高，計算過程簡便的優點，易于理解，易于操作，易于得出結果。其基本思路是運用多項式函數擬合邊緣灰度的變化趨勢，根據情況提取多項式函數相應的數學特征視作邊緣點。多項式函數的數學表達式如下：

圖1 Matlab軟件流程圖

I(z)=c0+c1z+c2z2+…+cnzn

其中:z為邊緣振動方向上的像素點的位置；I為z點處的灰度值；c0,c1,…,cn為擬合多項式中的系數，通過最小二乘法來確定。

圖像邊緣[14]的灰度變化有階躍型、脈沖型[15]和屋脊型，如圖2所示。

(a) 階躍型 (b) 脈沖型 (c) 屋脊型

圖2 常見的3種邊緣類型

本實驗的對象是索的邊緣，是典型的脈沖型邊緣，因而采用6階多項式函數擬合邊緣[16]振動方向上的若干點的灰度，選取擬合曲線的最大值[17]視作為邊緣實際位置。

3 模型試驗及結果分析

本實驗采用一個斜拉橋模型，選取它的第3根斜拉索進行試驗。索的表面光滑完整，兩端固定，采用直徑為1 mm，長2 m的單根不銹鋼絲，剛度未知，均布質量為6.228 kg/m。拍攝機器使用蘋果6手機，用三腳架固定，機身角度與索的角度一致，拍攝240幀/s1 280×720像素的視頻。試驗模型如圖3所示。

在實驗開始前先對視頻設備的拍攝幀數進行標定。使用本次視頻設備拍攝一段電子秒表計時。電子秒表以0.044 s為單位跳動，由24.571 s開始25.627 s結束，時長1.056 s，共分解得254張圖像。經計算得該視頻幀數約為240幀/s。

視頻拍攝從激振前開始，一直拍攝至索基本穩定，采集激振前作為初始圖像，采集振動穩定后作為振動圖像，拍攝時長為30 s。視頻采集后利用程序將視頻分解為圖像序列。經剪輯的視頻時長21 s，錄制幀數為240幀/s，一共可以分解得5 040幅圖像，運用圖片剪切程序將圖片裁剪成合適的尺寸，留作分析用。

圖3 模型試驗場景圖

(a)第1張截圖

(b)第2張截圖

(c)第254張截圖

現選圖像序列101 079～101 083作為索振動圖像序列演示，如圖5所示。

(a) 101 079

(b) 101 080

(c) 101 081

(d) 101 082

(e) 101 083

圖5 索振動序列圖像

選取101 070～101 110亞像素索邊緣圖像序列作為演示，如圖6所示。

利用Matlab程序處理圖像序列，取6階多項式擬合邊緣灰度曲線，選取索1/2點、索1/4點的位移時程曲線與頻譜圖作為演示，如圖7～10所示。索的前4階模態的振型圖如圖11所示。 對索各點的振動信號進行模態分析得斜拉索的前4階振動頻率，數據如表1所示。從表1得出，該索的前4階模態的頻率，根據索力公式[18]N=4wL2(fk/k)2計算得出。其中：w是斜拉索的單位長度的質量；L是索的長度；k是階數。可以根據各階模態的頻率算出索力F1，再與力傳感器所測得數據F2作出對比，如表2所示。

(a) 101 070～101 079

(b) 101 080～101 089

(c) 101 090～101 090

(d) 101 091～101 100

圖7 1/2索處點位移時程圖

圖8 1/2索處點頻譜圖

圖9 1/4索處點位移時程圖

圖10 1/4索處點頻譜圖

(a) 1階振型圖

(b) 2階振型圖

(c) 3階振型圖

(d) 4階振型圖

本次實驗能測出斜拉索的前4階模態的頻率，由表1可知，斜拉索的前4階模態的實測頻率為：f1=14.656 Hz、f2=29.061 Hz、f3=44.531 Hz、f4=58.360 Hz。

表2 圖像法索力值與力傳感器檢測索力值比

前4階模態的頻率基本滿足線性關系。由前4階模態的頻率計算出的索力與力傳感器所采集的數據最大相差2.6%，誤差在合理可接受范圍之內。

4 結 論

本文基于數字圖像測量技術，用相關硬件與Matlab軟件完成了非接觸式圖像測振方法的系統設計，并通過實驗室模型實現了斜拉索模型的振動測試，為其進一步應用于實際結構提供了參考，現將試驗的主要結論及建議總結如下:

(1) 對于動態位移的測量，由本文方法測得的試驗數據均較合理，試驗測出對應索模型的1～4階模態的頻率。由該4階模態的頻率所計算出的索力與力傳感器所得的數據比值誤差較小。因此基于普通數碼攝像頭與PC機的測振技術可用于高頻振動系統的測試，且可以通過測量得到四階模態參數。

(2) 基于數字圖像測量技術在振動測試中測試數據應用于索力測量是可行的。

(3) 數字圖像測量技術有著非接觸式，無負載反應的先天優點。其測試數據在空間密度方面遠勝與傳統方法。

基于無人機技術的運用、高像素視頻設備的普及、計算機硬件設備的進步，以及工程上對數字圖像測量技術的應用與反饋，數字圖像測量技術的立體化、精準度、實用性和適用范圍將進一步發展。

[1] 劉金秋.斜拉橋施工過程中索力的確定及其影響因素[D].大連:大連理工大學，2004.

[2] 劉文會，李 強.拉索松弛產生的時變效應對斜拉橋的影響[J]吉林建筑工程學院學報,2009,26(2)：31-32.

[3] 趙 翔.拉索損傷對斜拉橋結構性能影響的研究[D].南京：東南大學，2005.

[4] 劉 敏.數字圖像處理技術在橋梁結構檢測中的應用研究[D].廣州：廣州大學，2009.

[5] 袁向榮.梁變形檢測的一維數字圖像相關法[J]廣州大學學報，2010,9(1)：54-56.

[6] 胡朝輝.數字圖像處理技術在橋梁結構振動測試中的應用研究[D].廣州：廣州大學，2010.

[7] 袁向榮.邊緣識別的正交多項式擬合及梁變形檢測[J].實驗室研究與探索，2013,32(10)：11-13.

[8] 胡朝輝，袁向榮，劉 敏.振動試驗視頻圖像測試技術[J].噪聲與振動控制，2011(3)：162-165.

[9] 劉晶波，杜修力.結構動力學[M].北京：機械工業出版社，2005.

[10] 任文杰.圖像邊緣檢測方法的研究[D].濟南：山東大學，2008.

[11] Fu G K， Mitchell O R. Edge location to subpixel values in digital imagery[J].IEEE-PAMI,1984(6)：188-201.

[12] Tabatabai A J, Moosa A G.An optical approach to structural displacement measurement ant its application[J]. Engineering Echanices,2002,128(5):511-520.

[13] Gongkang F U, Moosa A G.An optical approach to structural displacement measurement ant its application[J]. Eng Mech ASCE, 2002, 128(5):511-520.

[14] 孫即祥.圖像分析[M].北京：科學出版社，2004.

[15] 盛遵冰，崔賢玉，高國安.通用亞像素邊緣檢測算法[J].上海交通大學學報，2007,41(6)：911-915.

[16] 陳 靜，尚雅層，田軍委.快速多項式擬合亞像素邊緣檢測算法的研究[J].應用光學，2011,32(1)：91-95.

[17] 陳小衛，徐朝輝，郭海濤，等.利用極值梯度的通用亞像素邊緣檢測方法[J].測繪學報，2014,43(5)：500-507.

[18] 唐盛華，方 志，楊 索.考慮邊界條件的頻率法測索力實用公式[J].湖南大學學報，2012,39(8)：7-13.

ApplicationofVideoImageVibrationMeasurementTechniqueinAnalysisofCableForceofStayCables

XUMinjie,YUANXiangrong,LIUHui,CHENZexian,RENZhangchen

(School of Civil Engineering, Guangzhou University, Guangzhou 510006, China)

the paper introduces a new technology for vibration which uses digital camera and PC as the hardware, and Matlab as software. The system is used to identify the first-order, second-order, third-order and fourth-order modal parameters for stay cables. The cable vibration video is collected, and is decomposed into the image sequence.Through image processing, cable edge sequence can be obtained;the edge vibration at each point in time history and vibration mode are then analyzed.Then cable frequency, damping and vibration type can be gotten, thecable forcecan be analyzed by frequency identification. Comparing with the data obtained byforce sensor, the video image vibration technology is feasible to detect cable vibration and has some advantages.With the further development of optical hardware and software, video image vibration technology has a brightprospect inengineering application.

stay cable； edge recognition； polynomial fitting； higher-order modes； cable force measurement

U 441.3

A

1006-7167(2017)11-0007-04

2017-03-02

國家自然科學基金項目(51078093,51278137)；廣州市科技計劃項目(12C42011564)

徐旻杰(1993-)，男，廣東珠海人，碩士生，研究方向為橋梁工程、數字圖像處理。Tel.:15521376900；E-mail: 1326195334@qq.com

袁向榮(1957-)，男，河北故城人，博士，教授，碩士生導師，研究方向為橋梁工程分析計算檢測、結構動力學分析與測試，數字圖像處理。E-mail: rongxyuan@163.com