基于單目離面視覺的結(jié)構(gòu)振動檢測

喻其炳　尹愛軍　薛　磊　張　泉

1.重慶工商大學(xué)裝備系統(tǒng)服役健康保障國際聯(lián)合研究中心, 重慶，4000672.重慶大學(xué)機(jī)械傳動國家重點(diǎn)實(shí)驗(yàn)室，重慶，400044

?

基于單目離面視覺的結(jié)構(gòu)振動檢測

喻其炳1尹愛軍2薛磊2張泉2

1.重慶工商大學(xué)裝備系統(tǒng)服役健康保障國際聯(lián)合研究中心, 重慶，4000672.重慶大學(xué)機(jī)械傳動國家重點(diǎn)實(shí)驗(yàn)室，重慶，400044

針對常規(guī)視覺測量中光學(xué)系統(tǒng)復(fù)雜、計算量大等問題，提出了一種離面視覺測量方法。討論了離面振動成像的基本原理，結(jié)合模態(tài)分析理論，研究了常規(guī)成像條件下單相機(jī)結(jié)構(gòu)振動特征參數(shù)的識別及無激勵振型分析。懸臂梁的離面視覺振動測試實(shí)驗(yàn)、傳統(tǒng)傳感器測量實(shí)驗(yàn)及ANSYS仿真的對比結(jié)果表明，所提出方法可以準(zhǔn)確識別低階固有頻率、阻尼和振型等模態(tài)信息。最后分析了該方法的后續(xù)研究方向。

離面視覺；視覺測量；振動參數(shù)識別；振型

0　引言

結(jié)構(gòu)振動參數(shù)識別及振型檢測是結(jié)構(gòu)振動分析的重要內(nèi)容。傳統(tǒng)振動分析方法一般采用單個或多個振動傳感器提取結(jié)構(gòu)振動信息[1]，檢測系統(tǒng)復(fù)雜，空間分辨率低,且在一定程度上可能會改變結(jié)構(gòu)的原有動力學(xué)特性。視覺測量是一種非接觸式振動測量方法，它具有全場測量、不改變被測物的振動特性等優(yōu)點(diǎn)，在微電子器件自動檢測、工業(yè)產(chǎn)品質(zhì)量在線監(jiān)控以及結(jié)構(gòu)三維形狀測量等方面已得到廣泛應(yīng)用[2]。

視覺測量方法中，基于結(jié)構(gòu)光、單攝像機(jī)以及多目立體視覺的測量方法已得到廣泛的研究和應(yīng)用[3-9]。Jing等[5]利用結(jié)構(gòu)光測量系統(tǒng)，對復(fù)雜三維輪廓進(jìn)行了高精度的尺寸測量并進(jìn)行了像素級的結(jié)構(gòu)三維輪廓重建。Teyssieux等[6]利用單個高速CCD相機(jī)和顯微成像系統(tǒng)，實(shí)現(xiàn)了MEMS懸臂梁的面內(nèi)運(yùn)動位移和高頻振動的精確測量。Wang等[7]通過兩個高清晰度、高頻相機(jī)得到結(jié)構(gòu)三維振動散斑圖像，利用數(shù)字圖像相關(guān)法實(shí)現(xiàn)了結(jié)構(gòu)的全場振動模態(tài)、應(yīng)變等測量。Seungbae等[8]利用高速相機(jī)獲取PCB電路板的離面振動圖像，得到電路板的振動模態(tài)信息以及振動幅值隨時間變化的分布圖。

上述結(jié)構(gòu)振動視覺分析方法，有些需要額外的結(jié)構(gòu)光投射系統(tǒng)；有些需要復(fù)雜昂貴的光路系統(tǒng)與高速高精度成像系統(tǒng)；有些則需要進(jìn)行大量的匹配運(yùn)算等。本文針對這些問題，提出一種常規(guī)成像條件下的離面視覺測量方法，該方法利用單個相機(jī)獲取結(jié)構(gòu)離面振動信息，結(jié)合實(shí)驗(yàn)?zāi)B(tài)分析理論，識別結(jié)構(gòu)振動特征參數(shù)和振型。最后進(jìn)行了懸臂梁沖擊響應(yīng)實(shí)驗(yàn)。與ANSYS有限元分析的對比結(jié)果表明，所提出方法能夠準(zhǔn)確有效地實(shí)現(xiàn)結(jié)構(gòu)振動分析。

1　離面視覺測量原理

1.1離面振動成像原理

根據(jù)針孔成像模型[10-12]，給出離面視覺振動檢測原理，如圖1所示。相機(jī)垂直振動方向成像，離面位移造成物距大小的改變，使得成像大小及位置發(fā)生改變，通過分析這種成像變化即可獲得結(jié)構(gòu)振動特性。

圖1　攝像機(jī)成像等效針孔模型

(1)

由式(1)可知，隨著物距的增加，離面位移引起的測量變量會顯著減小。在本文的實(shí)際測量中，物距b一般遠(yuǎn)大于振動位移w，故式(1)可以近似簡化為

(2)

此時，離面位移量與相機(jī)成像的像素位移大小之間成近似的線性關(guān)系。

1.2振動信息提取

邊緣特征檢測法是一般圖像處理問題的關(guān)鍵[13]。一般情況下，邊緣信息在圖像上呈現(xiàn)出不連續(xù)性，如灰度突變、亮度變化等。經(jīng)典的圖像邊緣檢測方法有Roberts算子、Sobel算子、Laplace算子以及Canny算子等，這些方法主要是使用圖像梯度的極大值或一階、二階導(dǎo)數(shù)來檢測圖像邊緣，通過微分算子模板和圖像進(jìn)行卷積來完成。Canny算子在噪聲抑制和邊緣檢測之間取得較好的平衡，具有很好的邊緣檢測性能。

如圖2a所示，在被測結(jié)構(gòu)上設(shè)置具有反差顏色的線狀邊緣標(biāo)志，經(jīng)過圖像濾波和Canny邊緣提取算法處理得到包含邊緣位置信息的圖像(圖2b)。由圖2可知，相鄰兩幀圖像中的邊緣位置分別為H1和H2，通過計算每幅圖像中邊緣線位置的差值ΔH=H1-H2，即可得到線狀邊緣標(biāo)志物的運(yùn)動信息，并根據(jù)式(2)獲得此時結(jié)構(gòu)的振動信息。

(a)　　　　　　　　　　　(b)圖2　邊緣提取示意圖

1.3模態(tài)參數(shù)識別方法

在模態(tài)分析中，通常將系統(tǒng)離散化為N個自由度的模型進(jìn)行分析[14]。設(shè)X=(x1,x2,…,xN)T和F=(F1,F2,…,FN)T分別為系統(tǒng)各點(diǎn)的位移響應(yīng)向量和激振力向量，X(ω)和F(ω)分別為位移響應(yīng)與激振力的傅里葉變換，有

X(ω)=H(ω)F(ω)

(3)

即

其中，H(ω)稱為頻響函數(shù)矩陣，是一個對稱陣。

在理想脈沖激勵下，系統(tǒng)的輸入為F(ω)=[F1F2…FN]T，特別地，當(dāng)在p點(diǎn)進(jìn)行單點(diǎn)激勵時，有F(ω)= [0…0Fp0…0]T，則由式(3)得系統(tǒng)第r階脈沖響應(yīng)為

(4)

其中，rHij(ω)為該點(diǎn)第r階模態(tài)對頻響函數(shù)的貢獻(xiàn)，rX(ω)=[rx1j(ω)rx2j(ω)…rxNj(ω)]T，

rxij(ω)為該點(diǎn)第r階頻域脈沖響應(yīng)函數(shù)。由式(4)可以看出，在理想脈沖激勵下，由系統(tǒng)脈沖響應(yīng)rX(ω)得到的模態(tài)參數(shù)與由系統(tǒng)傳遞函數(shù)rH(ω)得到的模態(tài)參數(shù)之間相差一個常數(shù)。

由模態(tài)分析理論知[1]，系統(tǒng)歸一化的第r階模態(tài)振型向量為

(5)

式中，rRij為系統(tǒng)中對應(yīng)點(diǎn)的第r階留數(shù)值。

本文采用的歸一化方法為：以模態(tài)振型向量中最大元素取值為1，rRsj為留數(shù)向量中的最大值。

當(dāng)系統(tǒng)各階模態(tài)較為分散時，在實(shí)模態(tài)條件下，由單模態(tài)識別法得系統(tǒng)中各點(diǎn)的第r階留數(shù)值為

rRij=-2σrImrHij(ωr)

(6)

圖3　單模態(tài)識別法的幅頻曲線圖

系統(tǒng)的第r階阻尼比為

(7)

模態(tài)振型實(shí)質(zhì)上是指結(jié)構(gòu)上各測點(diǎn)之間的相對位移關(guān)系，是一組比值，非絕對量[15]。結(jié)合式(4)～式(6)可知，由系統(tǒng)脈沖響應(yīng)rX(ω)得到的模態(tài)振型與由系統(tǒng)傳遞函數(shù)rH(ω)得到的模態(tài)振型在進(jìn)行歸一化之后是相同的。

2　實(shí)驗(yàn)設(shè)計及數(shù)據(jù)分析

2.1實(shí)驗(yàn)系統(tǒng)

(a)實(shí)驗(yàn)示意圖

(b)實(shí)驗(yàn)裝置圖圖4　實(shí)驗(yàn)系統(tǒng)

實(shí)驗(yàn)系統(tǒng)如圖4所示，其中相機(jī)為德國Baumer公司的TXG03cCCD相機(jī)，最高分辨率為656pixel×490pixel，最高分辨率下的最高幀率為90幀/s。本研究分別采用150幀/s(200pixel×200pixel分辨率)和60幀/s兩種幀率模式進(jìn)行實(shí)驗(yàn)。實(shí)驗(yàn)中，由于是離面測量，相機(jī)應(yīng)正對xy面放置，相機(jī)中心軸與z軸平行。由于相機(jī)越靠近固定端，測得的像素位移越小，為了保證一定的分析精度，本實(shí)驗(yàn)懸臂梁振動檢測點(diǎn)為A點(diǎn)。此外，為了后續(xù)數(shù)據(jù)分析的方便，將懸臂梁A點(diǎn)分析區(qū)域固定在相機(jī)視場內(nèi)，實(shí)驗(yàn)中相機(jī)與懸臂梁之間的距離為35cm。選取拾振點(diǎn)時要盡量避免拾振點(diǎn)在模態(tài)振型的節(jié)點(diǎn)上，該實(shí)驗(yàn)取拾振點(diǎn)為B敲擊點(diǎn)處。懸臂梁為長(x向)550mm、寬(y向)30mm、高(z向)3mm的合金鋁板，其密度為2900kg/m3，彈性模量為63GPa。在懸臂梁xy面中間，為黑白交界的特征線。

實(shí)驗(yàn)中采用的激勵裝置為橡膠帽脈沖力錘，其頻率特性如圖5所示，力錘頻譜在截止頻率fc=75Hz以內(nèi)是平穩(wěn)的[1]。

(a)錘擊產(chǎn)生的　 　　(b)錘擊沖擊力頻譜圖　　沖擊力時域圖圖5　脈沖力錘頻率特性

2.2數(shù)據(jù)分析

2.2.1振動頻率測量

圖6a為140pixel×140pixel分辨率下，幀率為150幀/s時，在B點(diǎn)進(jìn)行錘擊激勵，懸臂梁A點(diǎn)的沖擊響應(yīng)時域波形；根據(jù)傅里葉變換得其頻率響應(yīng)曲線，如圖6b所示；圖6c為t=2s時的原始圖像。

在懸臂梁上B點(diǎn)位置處加裝ICP壓電式加速度傳感器，使用LabView軟件進(jìn)行懸臂梁在錘擊激勵下的振動信號測試。ICP壓電式加速度傳感器型號為PCB-333B45，頻率范圍為10～3000Hz，靈敏度為51.7mV/(mm·s-1)，測量時采樣頻率為1652Hz。測量數(shù)據(jù)如圖7a所示，根據(jù)傅里葉變換得到其頻率響應(yīng)，如圖7b。

(a)時域響應(yīng)信號

(b)頻率響應(yīng)

(c)原始圖像圖6　A點(diǎn)沖擊響應(yīng)頻率特性及原始圖像

(a)加速度傳感器時域響應(yīng)信號

(b)加速度傳感器頻率響應(yīng)圖7　加速度傳感器時域響應(yīng)信號和頻域響應(yīng)

由系統(tǒng)頻率響應(yīng)得到該試樣的一階、二階固有頻率。表1為該懸臂梁的一階、二階固有頻率與ANSYS有限元分析的結(jié)果對比。由表1可以看出，檢測結(jié)果與仿真結(jié)果基本一致。由于在懸臂梁結(jié)構(gòu)上加裝了加速度傳感器，故在一定程度上降低了結(jié)構(gòu)固有頻率，導(dǎo)致檢測頻率偏低。對比本文視覺測量方法的數(shù)據(jù)和傳統(tǒng)傳感器法的數(shù)據(jù)可知，視覺測量方法準(zhǔn)確有效地測量出了懸臂梁的一階、二階振動頻率。

表1　實(shí)驗(yàn)與仿真結(jié)果對比

2.2.2振型檢測

為了提高振動測量圖像分辨率，采用650pixel×200pixel分辨率，幀率為60幀/s進(jìn)行脈沖激勵實(shí)驗(yàn)。由于相機(jī)視場限制的原因，將懸臂梁在x方向分成3等份(OA、AB、BC)，如圖4所示，分三次用CCD相機(jī)采集每段梁的離面振動圖像信息。三次實(shí)驗(yàn)中，均在B點(diǎn)脈沖錘擊激勵。每段梁在特征線上等間距選取13個點(diǎn)作為響應(yīng)特征點(diǎn)，根據(jù)模態(tài)參數(shù)識別的理論方法，由其沖擊響應(yīng)得到三個一階振型向量X1(OA段)、X2(AB段)、X3(BC段)。圖8為AB段梁的一階振型向量X2(未歸一化)以及點(diǎn)10處的時域響應(yīng)圖和頻譜圖。圖9為三段梁分別在三次激勵下所得的一階振型向量。

(a)AB段的振型向量

(b)點(diǎn)10處時域響應(yīng)

(c)點(diǎn)10處幅值頻譜圖8　AB段梁的振型向量和點(diǎn)10處的時域響應(yīng)圖、頻譜圖

圖9　OA、AB、BC三段梁的振型向量

由式(5)及式(6)可知，在理想脈沖激勵下，系統(tǒng)第r階振型向量可由各點(diǎn)該階脈沖響應(yīng)函數(shù)值得到。由圖3可知，在60Hz采樣頻率范圍內(nèi)，脈沖激勵頻率特性為常數(shù)。由式(4)可知，在三次實(shí)驗(yàn)中，不同大小的脈沖激勵信號(不同大小的fp)對系統(tǒng)響應(yīng)的影響只是在響應(yīng)幅值上相差一個比例常數(shù)，且結(jié)構(gòu)的振型是連續(xù)的，故設(shè)懸臂梁的一階振型向量為

φ1=(X1，αX2，βX3)T

(8)

其中，α、β為比例常數(shù)。

圖10　懸臂梁的一階振型向量及其擬合曲線

圖11　有限元仿真分析與實(shí)測擬合曲線對比

3　結(jié)語

視覺測量方法具有非接觸、大范圍測量等優(yōu)點(diǎn)。本文提出一種常規(guī)成像條件下的離面視覺測量方法，實(shí)現(xiàn)單相機(jī)下的結(jié)構(gòu)離面振動分析。通過對懸臂梁實(shí)際測量分析，及與有限元仿真數(shù)據(jù)對比分析可知，該方法能夠以較高的精度進(jìn)行結(jié)構(gòu)固有頻率、振型等振動參數(shù)的檢測。

本文提出的離面視覺測量方法仍需在如下方面進(jìn)行研究：① 由于相機(jī)幀率的限制，僅進(jìn)行了前兩階固有頻率和一階振型的測量，對于高頻振動難以測量；② 利用邊緣提取獲取振動信息具有局限性，在復(fù)雜結(jié)構(gòu)上不易應(yīng)用，應(yīng)研究更一般的圖像特征提取方法；③ 結(jié)合平面和離面測量，建立三維成像模型，進(jìn)行三維振動的測量；④在本文實(shí)驗(yàn)中，假定力錘的輸入信號為理想脈沖信號，對測量結(jié)果有一定影響；⑤后續(xù)工作將對離面視覺振動檢測的應(yīng)用場合、檢測極限、相機(jī)分辨率等進(jìn)行深入分析研究。

[1]李德葆,陸秋海. 實(shí)驗(yàn)?zāi)B(tài)分析及其應(yīng)用[M]. 北京: 科學(xué)出版社, 2001.

[2]張廣軍. 視覺測量[M]. 北京: 科學(xué)出版社, 2008.

[3]王超峰,劉利，汪國保，等. 高頻振動方向和振幅測量的模糊圖像處理方法[J]. 振動與沖擊, 2006, 24(5): 59-62.

Wang Chaofeng, Liu Li,Wang Guobao,et al. Processing Method of High Frequency Vibration Amplitude and Direction Measurement of Fuzzy Image Processing Method[J].Journal of Vibration and Shock,2006, 24(5): 59-62.

[4]季云峰. 無目標(biāo)計算機(jī)視覺技術(shù)在斜拉索振動測試中的應(yīng)用研究[J]. 振動與沖擊, 2013, 32(20):184-202.

JI Yunfeng. Application of Non-target Computer-vision-based Technique in Stay-cable Vibration Measurement[J].Journal of Vibration and Shock,2013, 32(20):184-202.

[5]Jing X,Shao L,An W,et al.An Absolute Phase Technique for 3D Profile Measuremen Tusing Four-step Structured Light Pattern[J].Optics and Lasersin Engineering, 2012, 50: 1274-1280.

[6]Teyssieux D, Euphrasie S, Cretin B. MEMS In-plane Motion/Vibration Measurement System Based CCD Camera[J]. Measurement, 2011, 44(10): 2205-2216.

[7]Wang W, Mottershead J E. Frequency Response Functions of Shape Features from Full-field Vibration Measurements Using Digital Image Correlation[J]. Mechanical Systems and Signal Processing, 2012, 28: 333-347.

[8]Seungbae P, Chirag S, Jae K, et al.Transient Dynamic Simulation and Full-field Test Validation for a Slim-PCBOf Mobile Phone under Drop/Impact[C]//Electronic Components and Technology Conference. Reno, 2007: 914-923.

[9]王宏偉,趙國慶，王玉軍. 物體振動幅度的遠(yuǎn)距離非接觸測量技術(shù)[J]. 振動與沖擊, 2009,28(11):67-69.

Wang Hong wei, Zhao Guoqing, Wang Yujun. No-ncontact Measurement Technology for Vibration Detection of Distant Object[J]. Journal of Vibration and Shock, 2009,28(11):67-69.

[10]Réthoré J, Morestin F, Lafarge L, et al. 3D Displacement Measurements Using a Singlecamera[J]. Optics and Lasers in Engineering, 2014, 57: 20-27.

[11]徐飛鴻, 林皓. 單目CCD離面位移測量[J].長沙理工大學(xué)學(xué)報(自然科學(xué)版), 2013, 10(1):54-60.

Xu Feihong,Lin Hao.Out-of-plane Displacement Measurement Based on Single CCD[J]. Journal of Changsha University of Science and Technology(Natural Science), 2013, 10(1):54-60.

[12]孫偉,何小元，鄭翔,等. 基于數(shù)字圖像相關(guān)的三維剛體位移測量方法[J]. 光學(xué)學(xué)報, 2008,28(5): 894-901.

SunWei,HeXiaoyuan,ZhengXiang,et,al.ThreeDimensionalRigidBodyDisplacementMeasure-mentBasedonDigitalImageCorrelation[J].ActaOpticaSinica,2008,28(5): 894-901.

[13]張根輩, 臧朝平. 基于振動測試的非線性參數(shù)識別方法[J]. 振動與沖擊, 2013,32(1):83-88.

ZhangGenbei,ZangChaoping.ANovelMethodforNonlinearParametricIdentificationBasedonVibrationTests[J].JournalofVibrationandShock,2013,32(1):83-88.

[14]段虎明. 實(shí)驗(yàn)?zāi)B(tài)分析的前端信號精度研究及虛擬式模態(tài)分析儀的研制[D]. 重慶：重慶大學(xué), 2008.

(編輯盧湘帆)

AnalysisofStructureVibrationBasedonOut-of-planeMonocularVision

YuQibing1YinAijun2XueLei2ZhangQuan2

1.ResearchCenterofSystemHealthMaintenance,ChongqingTechnologyandBusinessUniversity,Chongqing,400067 2.StateKeyLaboratoryofMechanicalTransmissions,ChongqingUniversity,Chongqing,400044

Conventionalvisionmeasurementhadproblemsofcomplexsystemorlargeamountofcalculations.Themeasurementprinciplesofout-of-planevisionwereresearched.Andcombiningthetheoryofmodalanalysis,thestructuralvibrationparameterswereidentifiedbyaCCDcameraunderthenormallightconditions.AndcomparedtheexperimentalresultsandANSYSshowsthemethodcanmeasurethestructurevibrationparameters(naturalfrequency,dampingandvibrationmode)accurately.Thefutureresearchesofthismethodweregivenintheendofthepaper.

out-of-planevision;visualmeasurement;vibrationparameteridentification;vibrationmode

2015-10-26

國家自然科學(xué)基金資助項(xiàng)目(51105396)；裝備系統(tǒng)服役健康保障國際聯(lián)合研究中心開放基金資助項(xiàng)目(KFJJ201501002)；重慶市科委項(xiàng)目(cstc2015jcyjA70007)

TH113

10.3969/j.issn.1004-132X.2016.17.016

喻其炳，男，1972年生。重慶工商大學(xué)裝備系統(tǒng)服役健康保障國際聯(lián)合研究中心高級工程師。研究方向?yàn)闄C(jī)械設(shè)計與制造。尹愛軍(通信作者)，男，1978年生。重慶大學(xué)機(jī)械工程學(xué)院教授、博士研究生導(dǎo)師。薛磊，男，1990年生。重慶大學(xué)機(jī)械工程學(xué)院碩士研究生。張泉，男，1991年生。重慶大學(xué)機(jī)械工程學(xué)院碩士研究生。