周云 曾雅麗思 趙瑜 易偉建

摘? ?要：針對結(jié)構(gòu)識別算法應用于實際工程時，結(jié)構(gòu)的轉(zhuǎn)角信息難于準確測量及轉(zhuǎn)角自由度通常容易被忽略的問題，本文研究了使用陀螺儀轉(zhuǎn)角傳感器測量動態(tài)信號的方法及響應信息不完備條件下的結(jié)構(gòu)物理參數(shù)識別.首先，針對結(jié)構(gòu)轉(zhuǎn)動響應信息測量困難這一問題，提出采用商業(yè)級的微機電系統(tǒng)（MEMS）陀螺儀傳感器測量角度和角速度響應，并基于最小二乘遞推算法對結(jié)構(gòu)物理參數(shù)識別方法進行了理論公式推導.然后以一座4層框架結(jié)構(gòu)為算例進行分析，設置由廣義逆方法重構(gòu)轉(zhuǎn)角和采用轉(zhuǎn)角真實值兩種工況，并對結(jié)構(gòu)物理參數(shù)進行識別，從而驗證了理論推導的正確性.同時，對兩種工況下所識別的物理參數(shù)進行比較，結(jié)果表明重構(gòu)轉(zhuǎn)動響應時物理參數(shù)識別的效果不夠理想，故考慮測量轉(zhuǎn)動響應.先對MEMS陀螺儀傳感器在受到?jīng)_擊振動下的動態(tài)精度進行了試驗驗證，在結(jié)構(gòu)的初位移小于10 mm時，動態(tài)角度測量的精度為0.1°.在此基礎上，通過一個3層2跨的鋼框架模型的動力試驗實測數(shù)據(jù)和分析結(jié)果，驗證了使用MEMS陀螺儀傳感器直接測量轉(zhuǎn)動響應相比于重構(gòu)轉(zhuǎn)動響應對彎剪型結(jié)構(gòu)進行剛度參數(shù)時域識別的效果更好.

關鍵詞：陀螺儀;動轉(zhuǎn)角;參數(shù)識別;時域;框架結(jié)構(gòu)

中圖分類號：TU317? ? ? ? ? ? ? ? ? ? ? ? 文獻標志碼：A

文章編號：1674—2974（2020）09—0010—13

Abstract：When the structure identification algorithm was applied to practical engineering， the structure angle information was difficult to accurately measure and the degree of freedom of rotation was often neglected. The method of measuring dynamic signals with gyroscope sensor and the identification of structural physical parameters under the condition of incomplete output information were studied in this paper. Firstly， a commercial Micro Electro Mechanical Systems （MEMS） gyroscope sensor was proposed to measure the angular velocity and angle responses， and the theoretical formulas of structural physical parameters identification method based on the least square recursive method were deduced. Then， taking a four-story frame structure as an example for analysis， two working conditions that the rotation angle reconstructed by the generalized inverse method and the adoption of the true value of rotation angle were set up， and the physical parameters was identified. The correctness of the theoretical derivation was verified. The results of physical parameters identified under two working conditions were compared， which reflected that the effect of physical parameter identification was not very ideal when the angular response was reconstructed， so we can consider measuring the rotational responses. The dynamic accuracy of MEMS gyroscope sensor under impact vibration was verified by experiments firstly. The results demonstrated that the accuracy of dynamic angle measurement was 0.1° when the initial displacement of the structure was less than 10 mm. On this basis， the measured dynamic test data and analysis results of a three-story and two-span steel frame model verified that using MEMS gyroscope sensor to directly measure the rotational response was better than using reconstructed rotational response for the time domain identification of stiffness parameters for the bending-shear structure.

Key words：gyroscope;angle of rotation;parameter identification;timedomain;frame structure

系統(tǒng)識別方法始于20世紀70年代，引入土木工程領域已有近40余年的發(fā)展歷史[1].該方法是工程結(jié)構(gòu)損傷診斷和狀態(tài)評估的重要手段，它可以從建筑物動力響應中識別出結(jié)構(gòu)的當前狀態(tài).結(jié)構(gòu)的性態(tài)通常用結(jié)構(gòu)模態(tài)參數(shù)（主要為自振頻率和振型）和結(jié)構(gòu)物理參數(shù)（主要為剛度參數(shù)）來描述. 結(jié)構(gòu)物理參數(shù)的變化相比模態(tài)參數(shù)的變化更為直觀地反映了結(jié)構(gòu)狀態(tài)的改變[2]. 進行結(jié)構(gòu)物理參數(shù)識別的方法，主要分為頻域識別和時域識別兩種，與頻域識別方法相比，時域識別更加直接，實測的時程數(shù)據(jù)中的結(jié)構(gòu)信息更加豐富，且增加時程響應采樣點數(shù)可以有效增加時域識別方程的個數(shù)[3]. 目前制約時域識別應用于實際工程中的一個主要原因是結(jié)構(gòu)響應信息不完備，通常時域方程中主要使用位移類信息，轉(zhuǎn)角類信息很難測試得到，往往只能測量結(jié)構(gòu)的一個狀態(tài)量，加速度、速度或位移響應的一種，而且轉(zhuǎn)角信息難以準確測量.如何在測量信息有限的情況下識別參數(shù)以及如何進行信息不完備情況下的識別問題是研究重點之一[4].

為了解決這些問題，眾多學者展開了相關研究.擴展卡爾曼濾波[5]方法可以求解響應信息不完備下的結(jié)構(gòu)參數(shù)識別問題，但結(jié)構(gòu)參數(shù)的識別精度低，甚至會出現(xiàn)參數(shù)發(fā)散現(xiàn)象.傳統(tǒng)方法主要通過重構(gòu)技術解決時域識別時測量信息的不完備問題.平動信息不完備的處理方法如下：Li等[6]在獲得系統(tǒng)的位移響應后，通過采用微分算子變換方法求得速度與加速度響應，完成系統(tǒng)的參數(shù)識別. 李國強等[2]采用積分算子變換方法，實現(xiàn)了在無需初速度信息的情況下，由加速度響應重建位移、速度響應.轉(zhuǎn)動信息的計算問題主要采用靜力凝聚法[2]和廣義逆方法[7]通過水平位移測量值重構(gòu)轉(zhuǎn)動信息，前者需要預先知道結(jié)構(gòu)的參數(shù)信息，而后者不需要，但是兩種方法計算效果近似，因此廣義逆方法更具適用性. 此外，趙昕等[8]采用一組正交基函數(shù)對梁式構(gòu)件的平動反應曲線進行擬合，從而獲得其轉(zhuǎn)角響應信息.

關于測量動態(tài)旋轉(zhuǎn)角度的研究不多，主要是因為轉(zhuǎn)動響應的最大振幅通常在毫弧度范圍以內(nèi)，所以要求儀器必須非常靈敏，而且需對測量的運動以外的其他運動非常不敏感[9]. 通過文獻調(diào)研發(fā)現(xiàn)，Ewins等[10]在文獻中對轉(zhuǎn)角自由度間接測量方法進行了詳細的分類介紹.包括采用附加T模塊法[11]、激光多普勒測振儀[12]、激光位移計[13]等測量轉(zhuǎn)角，但這些測量方法和設備過于復雜，難以運用于自由度較多的情況. Lee等[14]采用商用相機和低成本的幀捕獲器，利用事先標定參數(shù)的圖像處理技術，測量了五層試驗塔在激振器下的角度響應.然而，攝像機的方法需要額外的計算圖像處理，如像素掃描、目標識別和輪廓定義，實際應用并不方便. 也有學者嘗試采用應變推算轉(zhuǎn)角，Xia等[15]通過在廣州塔上安裝振動應變儀，利用實時監(jiān)測的應變數(shù)據(jù)計算彎曲梁型結(jié)構(gòu)頂部的傾斜角度，Dimsdale[9]采用在桿件兩端粘貼應變片的方式間接計算桿件兩端的彎曲轉(zhuǎn)角，該方法的應用前提是與其他桿件的連接方式不影響桿件的相對轉(zhuǎn)動.測量轉(zhuǎn)角自由度的傳感器技術遠不如測量位移自由度的傳感器成熟，目前采用較多的是傾角儀和陀螺儀，Hou等[16]利用MEMS傾角儀監(jiān)測鋼筋混凝土柱在地震荷載作用下的結(jié)構(gòu)旋轉(zhuǎn)角，該傳感器對被測物體的振動頻率非常敏感，僅適用于1 Hz以內(nèi)的外激勵頻率.傾角儀適合低頻的振動，而陀螺儀在高頻振動下的性能更為穩(wěn)定，是測量動態(tài)角度的理想傳感器.

由于轉(zhuǎn)角響應難以獲得，目前結(jié)構(gòu)識別算法的研究主要集中于剪切型結(jié)構(gòu)和桁架結(jié)構(gòu)[17-18]，而以彎剪變形為主的結(jié)構(gòu)研究較少.彎剪型結(jié)構(gòu)在正問題分析中通常利用靜力凝聚的方法縮聚掉轉(zhuǎn)角自由度，在反問題求解過程中，主要利用平動響應重構(gòu)轉(zhuǎn)動響應.目前關于如何通過直接測量得到動態(tài)角度信息的研究十分有限，測量轉(zhuǎn)角自由度的傳感器價格也普遍較高，不適合大量應用于工程結(jié)構(gòu)中.在時域識別方面，傳統(tǒng)算法在重構(gòu)狀態(tài)信息以及廣義逆求轉(zhuǎn)角時，計算模型所使用的方程式存在對沖效應，需要應用超定方程組修正結(jié)果，方法過于復雜且精度不能完全滿足要求.因此本文以一座4層框架結(jié)構(gòu)為例，在對該結(jié)構(gòu)進行轉(zhuǎn)動信息重構(gòu)的前提下，基于最小二乘遞推算法由結(jié)構(gòu)的平動響應數(shù)據(jù)在時域內(nèi)識別了結(jié)構(gòu)的物理參數(shù).對比運用轉(zhuǎn)動響應真實值的情形下所識別的物理參數(shù)，研究使用轉(zhuǎn)角重構(gòu)值對結(jié)構(gòu)物理參數(shù)時域識別的影響.為獲得準確的轉(zhuǎn)角測量值，本文利用MEMS陀螺儀傳感器，依據(jù)由平動響應信息擬合轉(zhuǎn)動響應的原理，對其在受到?jīng)_擊振動下的動態(tài)精度進行了試驗驗證. 在此基礎上，通過一個3層2跨的鋼框架模型的動力試驗，比較了采用廣義逆方法重構(gòu)轉(zhuǎn)角和直接測量轉(zhuǎn)角兩種情況下，阻尼參數(shù)和剛度參數(shù)的識別效果.

1? ?MEMS陀螺儀傳感器簡介

陀螺儀是一種可以用來測量被測物體相對于慣性空間的角速度和角位移的慣性器件.其基本的工作原理為，高速旋轉(zhuǎn)剛體的旋轉(zhuǎn)軸在不受外力影響時，其所指方向不改變. 陀螺儀在載體姿態(tài)測量方面具有穩(wěn)定性好、精度高、對環(huán)境的適應能力極強等諸多優(yōu)點，被廣泛應用于航空飛行器姿態(tài)測量、軍事軌跡測量、車輛駕駛控制、數(shù)碼產(chǎn)品如相機防抖、手機體感操作功能等領域.

MEMS陀螺儀使用現(xiàn)代的微機電技術，具有體積小、質(zhì)量輕、成本低、集成度高和精度較高的特點，具有良好的應用前景[19].本文的研究選用的是商業(yè)級的MEMS陀螺儀傳感器，深圳維特智能科技的HWT901B姿態(tài)測量傳感器，它集成了高精度的陀螺儀、加速度計、地磁場傳感器等模塊. MEMS陀螺儀傳感器通過采集加速度、角速度以及磁場信號的原始數(shù)據(jù)，使用高性能的微處理器與卡爾曼動態(tài)濾波算法，通過四元數(shù)方法解算得到圍繞X、Y、Z軸旋轉(zhuǎn)的俯仰角、偏航角和翻滾角. 因此，陀螺儀傳感器實際采集的是角速度響應，通過解算才能獲得角位移.HWT901B陀螺儀傳感器的性能參數(shù)如表1所示，其采集數(shù)據(jù)后，通過TTL電平串口進行數(shù)據(jù)傳輸，支持串口數(shù)字接口. 考慮到實際測量時距離較遠，研究團隊對其進行改造加入了HC-12無線收發(fā)傳輸串口通訊模塊，整套動態(tài)角度測量系統(tǒng)的連接方法如圖1所示.

2? ?基于轉(zhuǎn)角擬合的角度測量精度校驗

MEMS陀螺儀傳感器體積微小，內(nèi)部含有多個運動單元，因此其性能對工作環(huán)境變化較為敏感.測量過程中的外界沖擊振動，環(huán)境干擾和自身漂移率都會影響陀螺儀的精度.目前，可通過提高微加工工藝消除尺寸偏差[20]、合理設計優(yōu)化陀螺儀結(jié)構(gòu)降低耦合誤差[21]、建立相應誤差模型進行對準[22]、采用數(shù)字信號處理技術減小隨機噪聲的影響[23]等方法來有效提高其精度，HWT901B正是采用卡爾曼動態(tài)濾波算法減少了其漂移誤差.考慮到動力試驗時，結(jié)構(gòu)將產(chǎn)生高頻的往復振動，有必要對MEMS陀螺儀傳感器在受到外界沖擊振動下的動態(tài)精度進行校驗，一般采用陀螺測試轉(zhuǎn)臺來檢驗和評定陀螺儀性能，轉(zhuǎn)臺為精密測試設備，一般土木工程應用很少，考慮購買轉(zhuǎn)臺或?qū)ν勇輧x的精度進行檢測的送檢價格都過高，受實驗條件所限制，故采用了一種角度擬合方法[8]對MEMS陀螺儀傳感器的精度進行了驗證，通過對激光位移計實測的結(jié)構(gòu)平動響應曲線擬合得到結(jié)構(gòu)的轉(zhuǎn)動響應曲線.

如圖2所示的等截面梁單元，它聯(lián)結(jié)著兩個節(jié)點，一般不考慮軸向變形的影響，則該單元在平面內(nèi)有4個自由度，包含2個平動自由度和2個轉(zhuǎn)動自由度.設有如圖3所示的受到脈沖激勵F，跨度為L的歐拉-伯努利梁，其撓度曲線上橫坐標為x的任一點的縱坐標為y，已知t時刻該梁上n個測點的平動反應測量值為

在有限單元法中，梁單元變形曲線為三次多項式，對于荷載作用下的歐拉梁，應當取不低于三次的多項式擬合梁的撓曲線.為求解三次多項式中的4個參數(shù)，應至少在梁上布置4個測點.

根據(jù)以上方程可知，只要通過激光位移計測試得到t時刻4個測點的位移y1、y2、y3、y4，基于梁單元變形曲線為三次多項式的假定，且已知4個測點的位置x1、x2、x3、x4，可以擬合得到梁在t時刻的撓曲線方程（1），求解出a1、a2、a3，將其代入式（2）可求得梁在t時刻的轉(zhuǎn)角反應，即可計算出t時刻橫坐標為xk處的轉(zhuǎn)角θ.

3? ?基于最小二乘遞推算法的結(jié)構(gòu)物理參數(shù)時域識別

框架結(jié)構(gòu)是由梁、柱構(gòu)件組成的結(jié)構(gòu)，可根據(jù)橫梁與柱的線剛度比值劃分為剪切型框架和一般框架結(jié)構(gòu).若線剛度比值大于5，則在結(jié)構(gòu)分析中，可假定橫梁剛度無窮大，忽略梁的變形，視為剪切型結(jié)構(gòu)[24]. 若框架梁柱的線剛度比值較小，則為一般框架結(jié)構(gòu).剪切型框架的參數(shù)識別已經(jīng)比較成熟，本文研究對象為一般框架結(jié)構(gòu)，在水平荷載作用下，其框架節(jié)點將發(fā)生轉(zhuǎn)動，因此需要考慮框架節(jié)點轉(zhuǎn)動的影響.

如圖4（a）所示的完整自由度模型一，對于鋼框架結(jié)構(gòu)的梁柱連接節(jié)點，節(jié)點的軸向變形和剪切變形與它的轉(zhuǎn)動變形相比，可以忽略不計，因此框架結(jié)構(gòu)計算模型采用梁單元的剛度矩陣為4 × 4階矩陣.對于單跨框架結(jié)構(gòu)，由于結(jié)構(gòu)對稱，左右兩側(cè)的轉(zhuǎn)動響應一致，可將模型一簡化為如圖4（b）所示的含N個平動自由度和N個轉(zhuǎn)動自由度的模型二. 對于多跨框架結(jié)構(gòu)，同一層梁柱節(jié)點的轉(zhuǎn)角響應并不完全相同，根據(jù)實際情況，可將模型一簡化為如圖4（c）所示的含N個平動自由度和2 N個轉(zhuǎn)動自由度的模型三.

在外力f（t）作用下，對于需要考慮節(jié)點轉(zhuǎn)動影響的框架結(jié)構(gòu)，動力平衡方程為

上述論述表明：利用系統(tǒng)的輸入、輸出響應信息可以得到結(jié)構(gòu)的物理參數(shù)估計值.

3.1? ?無轉(zhuǎn)角測量值的求解——結(jié)構(gòu)的轉(zhuǎn)角信息重構(gòu)

以相對簡單的單跨框架結(jié)構(gòu)為例說明，在沒有任何實測的轉(zhuǎn)動信息的基礎上，為了進行式（11）的計算，可采用廣義逆理論由水平位移測量值計算轉(zhuǎn)角，計算模型選取模型二.

由于該式中的Kφx 和Kφφ包含有待求參數(shù)，因此無法直接用來計算轉(zhuǎn)角響應. 陳雋[7]在文獻中給出了廣義逆方法計算轉(zhuǎn)角的公式：

由平動信息重構(gòu)角度、角速度、角加速度信息后，可進行結(jié)構(gòu)的參數(shù)識別.

3.2? ?有轉(zhuǎn)角測量值的求解

對于可以直接測量轉(zhuǎn)角值的框架結(jié)構(gòu)，同樣以模型二為例說明，選取的待識別參數(shù)為各層的層抗彎剛度（EI）i和各層的阻尼參數(shù)b（EI）i. 根據(jù)前述理論介紹，可構(gòu)造如式（10）所示的方程，第i層和第i+1層的參數(shù)識別方程如下式：

由式（19）可知，已知框架結(jié)構(gòu)的各層層高、質(zhì)量矩陣、所受外激勵以及各節(jié)點的平動、轉(zhuǎn)動時程響應，進而建立參數(shù)識別標準方程進行最小二乘遞推法參數(shù)識別，即可求得結(jié)構(gòu)的物理參數(shù).有轉(zhuǎn)角測量值即可省去采用廣義逆方法重構(gòu)轉(zhuǎn)角信息的過程，更為方便.

4? ?算例分析

以一個4層彎剪型框架結(jié)構(gòu)為例，在結(jié)構(gòu)左側(cè)頂點水平施加一個大小為100 N的脈沖力F，采樣頻率設為200 Hz，共采集12 000個數(shù)據(jù)點. 計算模型如圖5所示，假設參數(shù)如表2所示.采用Newmark-β時程分析方法進行正問題計算得到結(jié)構(gòu)的動力響應.動力方程中采用與剛度相關的比例阻尼矩陣和協(xié)調(diào)質(zhì)量矩陣.采用忽略軸向變形和剪切變形的4×4階的桿系單元矩陣，組裝成16×16階的總剛度矩陣.模型共16個自由度，由于結(jié)構(gòu)完全對稱，可簡化為含4個平動自由度和4個轉(zhuǎn)動自由度的框架模型，則總剛度矩陣變?yōu)?×8階的矩陣. 為了驗證算法的抗噪性能，在結(jié)構(gòu)響應分析中分別加入1%、5%的白噪聲，加入白噪聲的水平由結(jié)構(gòu)響應的均方根值和噪聲均方根值的百分比確定.接著使用最小二乘遞推算法對1～4層的層間抗彎剛度K1、K2、K3、K4和各層的剛度阻尼參數(shù)C1、C2、C3、C4進行識別.

根據(jù)實際觀測輸出信息，設置通過廣義逆方法重構(gòu)轉(zhuǎn)角的工況1和采用轉(zhuǎn)角真實值的工況2. 其中輸出的加速度響應和角度、角速度響應都采用真實值.

工況1：采用Newmark-β積分法計算結(jié)構(gòu)上4個平動自由度的加速度響應，對加速度響應進行傅里葉變換，將信號從時域轉(zhuǎn)換到頻域，再通過將時域內(nèi)的積分轉(zhuǎn)換為頻域計算，在頻域內(nèi)進行運算過后再通過傅里葉逆變換將信號從頻域還原到時域，并取其實部就可以得到信號在時域的速度和位移響應[25]，采用3.1節(jié)所述廣義逆方法，根據(jù)式（15）（17）和（18）由位移響應重構(gòu)角度、角速度、角加速度響應.

工況2：基于圖5的計算模型，采用Newmark-β積分法計算結(jié)構(gòu)上4個平動自由度的加速度響應，4個轉(zhuǎn)動自由度的角度、角速度響應.然后，同工況1的方法由頻域積分法獲得速度和位移響應.對角速度響應利用向前差分的數(shù)值微分法重構(gòu)角加速度的響應，初始角速度為0.

獲得完備的輸入和輸出響應信息后，假設結(jié)構(gòu)的待識別參數(shù)初值，將以上數(shù)據(jù)和結(jié)構(gòu)的已知參數(shù)代入公式（19）構(gòu)造結(jié)構(gòu)參數(shù)識別方程，根據(jù)式（11）運用最小二乘遞推算法可求出結(jié)構(gòu)物理參數(shù)的識別值，參數(shù)識別結(jié)果列入表3，圖6為不考慮噪聲下，工況1和工況2識別結(jié)構(gòu)物理參數(shù)的收斂過程.

比較工況1和工況2在無噪聲時的結(jié)構(gòu)參數(shù)識別結(jié)果，可以發(fā)現(xiàn)工況1重構(gòu)轉(zhuǎn)角響應時剛度參數(shù)識別的最大誤差為13.04%;而工況2采用角度和角速度響應真值，識別得到的剛度參數(shù)最大誤差為0.06%.說明通過廣義逆方法由平動響應重構(gòu)出的角度，角速度和角加速度響應準確度不高，基于最小二乘遞推算法進一步進行識別的剛度參數(shù)結(jié)果也不夠理想.從識別的物理量來看，阻尼的識別效果整體要差于剛度參數(shù)，考慮噪聲后誤差更大，阻尼真值很容易被測試噪聲所湮沒.

考慮5%的噪聲存在時，對轉(zhuǎn)角進行重構(gòu)的工況1的剛度參數(shù)識別最大誤差達到21.13%，分析其原因，主要是因為加速度積分成速度和位移產(chǎn)生的誤差和廣義逆方法本身的誤差累積，使得誤差進一步放大.考慮5%噪聲時，采用角度、角速度真實值的工況2剛度參數(shù)識別的最大誤差為2.94%，說明了最小二乘遞推算法具有一定的抗噪能力.整體來說，采用重構(gòu)轉(zhuǎn)角值的工況1所識別的參數(shù)誤差較大，因此可以考慮通過傳感器直接測量轉(zhuǎn)動響應.

5? ?試驗研究

5.1? ?MEMS陀螺儀傳感器的角度測量精度驗證

有效識別物理參數(shù)的前提是保證測量的響應信息準確，因此需對動力試驗條件下MEMS陀螺儀傳感器的角度測量精度進行驗證. 試驗結(jié)構(gòu)模型為如圖7所示的門式鋼框架結(jié)構(gòu)，選取其右側(cè)柱構(gòu)件為試驗對象. 柱尺寸為B × W × L = 5 0 mm × 3 mm ×1 060 mm，角度信號由維特智能HWT901B陀螺儀傳感器采集，共布置2個MEMS陀螺儀傳感器，位于門式框架右側(cè)柱構(gòu)件中間位置的傳感器用來采集柱的轉(zhuǎn)動響應，左側(cè)柱構(gòu)件上對稱布置的傳感器用來保持框架結(jié)構(gòu)質(zhì)量對稱.位移信號由基恩士IL300+IL1000激光位移計采集，其量程為280 mm，精度為0.7 mm. 試驗共布置5個激光位移計，頂端布置1個用于控制初位移，后4個以23 cm為間距等距離布置，用于測量柱的平動響應.試驗時，通過左側(cè)絞盤牽引施加2 ～ 10 mm不等的初位移.然后剪斷繩子瞬時釋放，結(jié)構(gòu)恢復自由振動狀態(tài). 每隔2 mm設置一次工況，共5組試驗. 由第2節(jié)所述的轉(zhuǎn)動響應曲線擬合方法可以得到激光位移計測得的轉(zhuǎn)動時程響應，與MEMS陀螺儀傳感器測得的轉(zhuǎn)動時程響應結(jié)果對比如圖8和圖9所示. （篇幅所限，僅展示初位移分別為4 mm和10 mm時的測量結(jié)果.）

通過上述圖片中MEMS陀螺儀傳感器與激光位移計輸出的角度的對比，可以明顯觀察出：隨著初位移的增大，MEMS陀螺儀傳感器與激光位移計輸出的角度差值也隨之增大.分析誤差產(chǎn)生的原因主要為以下兩點：1）根據(jù)動轉(zhuǎn)角擬合原理，應保證激光位移計發(fā)射的激光，門式框架受激勵后振動的二維平面，絞盤上繩子施加力的方向均在同一個平面內(nèi).由于現(xiàn)場試驗條件限制，很容易產(chǎn)生誤差. 2）激光位移計所測量構(gòu)件平面的表面不夠平整. 綜合來看，在結(jié)構(gòu)初位移為10 mm以內(nèi)，可以基本保證MEMS陀螺儀傳感器進行動態(tài)角度測量的精度在0.1°.

5.2? ?MEMS陀螺儀傳感器的時域識別

5.2.1? ?試驗模型

試驗結(jié)構(gòu)模型為如圖10所示的一個3層2跨的鋼框架結(jié)構(gòu)，框架模型整體尺寸為1 m × 1.6 m，首層層高為0.6 m，二、三層層高為0.5 m，框架柱的間距0.5 m，梁、柱均采用50 mm × 3 mm的截面尺寸的鋼薄板. 梁柱節(jié)點處采用尺寸為∟18 mm × 18 mm × 3 mm的角鋼加M6（螺紋直徑6 mm）×16（螺紋長度16 mm）的螺栓連接，如圖11所示. 框架模型構(gòu)件采用Q235鋼，取0.5 m長的鋼梁構(gòu)件做測試，實測質(zhì)量密度為7 832.60 kg/m3，用脈沖錘擊法識別梁的固有頻率，然后利用固有頻率識別梁的彈性模量，經(jīng)過10次測試平均得到鋼梁的彈性模量為2.072 × 1011 N/m2，基于彎剪型結(jié)構(gòu)計算模型，可得鋼框架結(jié)構(gòu)的層間抗彎剛度參數(shù)EI的計算值為22.622 1 N/m.

5.2.2? ?測試方法及結(jié)果

用脈沖錘擊法對框架結(jié)構(gòu)進行動力測試，由于結(jié)構(gòu)完全對稱，可將框架左側(cè)與右側(cè)梁柱連接節(jié)點的轉(zhuǎn)動響應視為相等，僅在一側(cè)布置MEMS陀螺儀傳感器.沖擊激勵使用美國PCB力錘，加速度信號由PCB系列壓電式加速度傳感器333B40采集，角度及角速度信號由維特智能HWT901B陀螺儀傳感器采集，使用LMS Cadax-8動態(tài)信號分析儀進行試驗數(shù)據(jù)結(jié)果的采集. 為保證框架結(jié)構(gòu)質(zhì)量對稱，在每根梁上未布置MEMS陀螺儀傳感器的一側(cè)，均放置了100 g的配重磁鐵用于模擬MEMS陀螺儀傳感器的質(zhì)量.

試驗中采用力錘在鋼框架的3#節(jié)點位置水平從左至右激勵一次，力錘的錘擊力測量結(jié)果經(jīng)放大后如圖12（a）所示，采樣頻率為100 Hz，采樣點數(shù)為12 000個.

傳感器測點布置如圖11所示，加速度傳感器布置于7#、8#、9#測點，采集得到的加速度響應信號如圖12（b）所示. MEMS陀螺儀傳感器布置在4#至9#測點，采集的框架中間4#、5#、6#測點的角度和角速度響應信號分別如圖12（c）和圖12（d）所示. （篇幅所限，未展示框架右邊7#、8#、9#測點的角度和角速度信號）. 對測試得到的各層加速度響應分別進行傅里葉變換，對相應的各階頻率取平均值得到結(jié)構(gòu)的前三階模態(tài)頻率，分別為2.073 Hz、6.704 Hz、11.773 Hz.利用采集得到的結(jié)構(gòu)自由振動時的加速度時程曲線的峰值的對數(shù)衰減率求得阻尼比，并根據(jù)第1和第2階振型的頻率，可以得到剛度阻尼系數(shù)b，乘以剛度參數(shù)EI可以得到阻尼參數(shù)的計算值為0.024 4 N·s/m.

5.2.3? ?結(jié)構(gòu)參數(shù)時域識別

框架模型的梁與柱構(gòu)件等剛度，可視為一般框架結(jié)構(gòu).模型中所有梁柱構(gòu)件具有相同的截面尺寸和材料特性，梁的跨度相等，整體結(jié)構(gòu)對稱，基于以下假定對其進行參數(shù)識別：1）框架模型的材料為線彈性;2）梁與柱連接節(jié)點處采用角鋼加螺栓連接，假設梁與柱的連接有足夠的剛度，如圖13所示，當框架發(fā)生彈性變形時，能保持相交的梁柱桿件之間原有的角度不變[26]. 3）框架底層用膨脹螺絲固定于地面，在結(jié)構(gòu)分析中忽略框架的平面外變形.4）梁與柱的軸向變形和剪切變形較之轉(zhuǎn)動變形可以忽略，因此假設框架同一層內(nèi)的各梁柱連接節(jié)點的水平位移均相等.

試驗結(jié)構(gòu)模型共18個自由度，包含9個平動自由度和9個轉(zhuǎn)動自由度.基于同一層內(nèi)節(jié)點的水平位移均相等的假設，可將上部結(jié)構(gòu)的平動自由度簡化成3個. 考慮轉(zhuǎn)角重構(gòu)和轉(zhuǎn)角實測兩種工況，由于前者將由平動響應重構(gòu)轉(zhuǎn)動響應，所以將9個轉(zhuǎn)動自由度簡化為3個，后者則考慮所有的轉(zhuǎn)動自由度.最終兩個工況的計算模型如圖14所示.

已知MEMS陀螺儀傳感器和配重磁鐵的質(zhì)量均為100 g，根據(jù)鋼梁構(gòu)件的尺寸和質(zhì)量密度，可求出梁單元的質(zhì)量，同理可求出柱單元的質(zhì)量，進而可求得4×4階的單元質(zhì)量矩陣. 將單元質(zhì)量矩陣組裝后可獲得框架結(jié)構(gòu)最終協(xié)調(diào)質(zhì)量矩陣，轉(zhuǎn)角重構(gòu)和轉(zhuǎn)角實測工況分別采用6×6階和12×12階的協(xié)調(diào)質(zhì)量矩陣.根據(jù)第3節(jié)所述方法，可以構(gòu)造參數(shù)識別的方程式（19），轉(zhuǎn)角重構(gòu)工況下直接采用3.2節(jié)的公式，轉(zhuǎn)角實測工況下，當識別第1、2層時，式（19）中：

將響應信息、外激勵、質(zhì)量矩陣和層高代入上式，由最小二乘遞推算法識別出結(jié)構(gòu)的剛度和阻尼參數(shù)如表4所示，剛度參數(shù)與阻尼參數(shù)的收斂曲線如圖15所示.

從識別結(jié)果可看出，轉(zhuǎn)角實測工況下的剛度參數(shù)識別效果明顯優(yōu)于轉(zhuǎn)角重構(gòu)工況，轉(zhuǎn)角重構(gòu)工況下的最大和最小誤差分別為20.36%和14.38%，轉(zhuǎn)角實測工況下的最大和最小誤差分別為13.92%和7.84%.阻尼參數(shù)的識別結(jié)果受噪聲影響較大，失去參考價值.

分析轉(zhuǎn)角實測工況下剛度參數(shù)識別的誤差產(chǎn)生的原因有以下可能;其一，由加速度響應重構(gòu)速度、位移響應的過程中，由于存在趨勢項和噪聲，會使得積分結(jié)果產(chǎn)生誤差.其二，由于該框架結(jié)構(gòu)的轉(zhuǎn)動變形較小，而MEMS陀螺儀傳感器的動態(tài)測量精度僅為0.1°，因此識別效果不夠理想.其三，受實驗條件限制，由于MEMS陀螺儀傳感器為接觸式傳感器，不能完全測得節(jié)點連接處的實際轉(zhuǎn)動變形，只能以節(jié)點附近的轉(zhuǎn)動變形近似替代.但使用MEMS陀螺儀傳感器測量的轉(zhuǎn)角精度仍然高于用廣義逆方法重構(gòu)出的轉(zhuǎn)角，這為轉(zhuǎn)動信息的準確測量提供了新的思路.

6? ?結(jié)? ?論

基于廣義逆方法重構(gòu)轉(zhuǎn)角信息和采用撓曲線擬合的方式重構(gòu)轉(zhuǎn)角信息，提出采用商業(yè)級的MEMS陀螺儀傳感器測量角度和角速度響應，改變了傳統(tǒng)的僅測試位移類響應的方法，對響應信息進行了直觀補充.然后通過最小二乘遞推算法識別彎剪型結(jié)構(gòu)的剛度和阻尼參數(shù)，得出的結(jié)論如下：

1）基于撓曲線擬合方法對MEMS陀螺儀傳感器的角度測量精度校驗，主要用來驗證動力試驗條件下測量的角度是否準確.本文采用的陀螺儀傳感器HWT901B在初位移為10 mm內(nèi)可保證動態(tài)角度精度為0.1°，試驗時，測得鋼框架結(jié)構(gòu)的最大位移為4 mm左右，受到力錘激勵后的最大角位移為0.3°左右，其動態(tài)角度精度約為0.02°，精度誤差約為6.67%，該傳感器基本可以滿足實驗要求.

2）實際工程中的結(jié)構(gòu)形式更復雜，剛度也更大，相應的轉(zhuǎn)動響應的振幅更小，因此對使用的傳感器靈敏度要求更高，花費也相對較高.部分位置如節(jié)點連接處的轉(zhuǎn)角如何準確測量、如何采用信號處理算法及誤差補償技術提高MEMS陀螺儀的性能，從而提高其精度仍是值得進一步研究的課題.

3）根據(jù)算例分析和試驗研究發(fā)現(xiàn)，基于廣義逆理論重構(gòu)轉(zhuǎn)動信息的方法需要使用位移、速度、加速度響應分別來重構(gòu)角位移、角速度、角加速度響應. 試驗時，通常只測量位移類響應中的加速度響應，加速度響應積分成速度、位移響應的過程中會產(chǎn)生誤差. 由位移類響應重構(gòu)轉(zhuǎn)角類響應的過程中，會將在位移類響應中累積的誤差帶入轉(zhuǎn)角類響應. 同時，廣義逆方法本身也存在誤差. 以上原因，導致重構(gòu)轉(zhuǎn)角下的剛度參數(shù)識別效果不如直接測量轉(zhuǎn)角理想.

4）采用MEMS陀螺儀傳感器直接測量轉(zhuǎn)動響應在實測框架中的平均誤差為11.13%.相比于由廣義逆方法重構(gòu)轉(zhuǎn)動響應的平均誤差17.07%，誤差減少了5.94%.說明使用最小二乘遞推算法對彎剪型結(jié)構(gòu)進行剛度參數(shù)時域識別時，使用轉(zhuǎn)角響應測量值的方法更好，而且與廣義逆重構(gòu)轉(zhuǎn)動響應的方法相比，無需進行復雜計算也不受平動響應的影響.

本文在常規(guī)的僅測量平動信息的基礎上補充了轉(zhuǎn)動信息的實測數(shù)據(jù)，使時域識別的結(jié)構(gòu)物理參數(shù)更為準確，對結(jié)構(gòu)的健康監(jiān)測具有一定參考價值.

參考文獻

[1]? ? HART G C，YAO J T P. System identification in structural dynamics[J]. Journal of the Engineering Mechanics Division，1977，103（6）：1089—1104.

[2]? ? 李國強，李杰. 工程結(jié)構(gòu)動力檢測理論與應用[M]. 北京：科學出版社，2002：185—191.

LI G Q，LI J. Theory and application of dynamic tests in engineering[M]. Beijing： Science Press，2002：185—191. （In Chinese）

[3]? ? 謝獻忠. 結(jié)構(gòu)動力學系統(tǒng)時域辨識理論與試驗研究[D]. 長沙：湖南大學土木工程學院，2005：1—2.

XIE X Z. Theoretical and experimental researeh ondynamic systemidetification in time domain[D]. Changsha： College of Civil Engineering，Hunan University，2005：1—2. （In Chinese）

[4]? ? 趙昕. 工程結(jié)構(gòu)時域系統(tǒng)識別研究[D]. 上海：同濟大學土木工程學院，2003：7.

ZHAO X. Research on time domain system identification of engineering structures[D]. Shanghai：College of Civil Engineering，Tongji University，2003：7.（In Chinese）

[5]? ? 潘芹，易偉建. 卡爾曼濾波在結(jié)構(gòu)物理參數(shù)識別中的應用[J]. 湖南大學學報（自然科學版），2002，29（3）：137—142.

PAN Q，YI W J. Application of extended Kalman filter toidentification of structural parameters[J]. Journal of Hunan University （Natural Sciences），2002，29（3）：137—142. （In Chinese）

[6]? ? LI J，ROBERTS J B. Stochastic structural system identification，Part 1：mean parameter estimation[J]. Computational Mechanics，1999，24（3）：211—215.

[7]? ? 陳雋. 高層建筑損傷檢測中的復合反演理論與試驗研究[D]. 上海：同濟大學建筑工程系，1999：63—66.

CHEN J. Compound inversion theory and experimental study on damage detection of tall buildings[D]. Shanghai：College of Architectural Engineering，Tongji University，1999：63—66. （In Chinese）

[8]? ? 趙昕，李杰. 轉(zhuǎn)角信息未知條件下的結(jié)構(gòu)參數(shù)識別方法研究[J]. 工程力學，2003，20（4）：55—59.

ZHAO X，LI J. Structural parameter estimation with unknown rotational response[J]. Engineering Mechanics，2003，20（4）：55—59.（In Chinese）

[9]? ? DIMSDALE J S. System identification of structures with joint rotation[M]. Berkeley：University of California，1983：34—35.

[10]? EWINS D J，SAINSBURY M G. Mobility measurements for the vibration analysis of connected structures[J]. Shock and Vibration Bulletin，1972，44：105—122.

[11]? MOTTERSHEAD J E，KYPRIANOU A，OUYANG H. Estimation of rotational frequency responses[J]. Key Engineering Materials，2003，245/246：157—166.

[12]? RATCLIFFE M J，LIEVEN N A J. Measuring rotational degreesof freedom using a laser Doppler vibrometer[J]. Journal of Vibration and Acoustics，2000，122：12—20.

[13]? SUN J H，ZHANG J，LIU Z，et al. High-dynamic angle measurement based on laser displacement sensors[J]. Applied Optics，2013，52（23）：5676—5685.

[14]? LEE J J，HO H N，LEE J H. A vision-based dynamic rotational angle measurement system for large civil structures[J]. Sensors，2012，12（12）：7326—7336.

[15]? XIA Y，ZHANG P，NI Y Q，et al. Deformation monitoring of a super-tall structure using real-time strain data[J]. Engineering Structures，2014，67：29—38.

[16]? HOU S，ZENG C S，ZHANG H B，et al.Monitoring interstory drift in buildings under seismic loading using MEMS inclinometers[J]. Construction and Building Materials，2018，185：453—467.

[17]? XU B，CHEN G，WU Z S. Parametric identification for a truss structure using axial strain[J]. Computer-Aided Civil and Infrastructure Engineering，2007，22（3）：210—222.

[18]? SERDAR S，MARIA Q F. Instantaneous damagedetection of bridge structures and experimentalverification[J]. Journal of Structural Control and Health Monitoring，2008，15（7）：958—973.

[19]? 楊業(yè)飛，申文濤. 慣性穩(wěn)定平臺中陀螺技術的發(fā)展現(xiàn)狀和應用研究[J]. 飛航導彈，2011（2）：72—79.

YANG Y F，SHENG W T. Development and application of gyroscope technology in inertial stabilization platform[J].Aerodynamic Missile Journal，2011（2）：72—79.（In Chinese）

[20]? VISWANATH A，LI T，GIANCHANDANI Y B. High resolution micro ultrasonic machining （HR-USM） for post-fabrication trimming of fused silica 3-D microstructures[C]// IEEE International Conference on MICRO Electro Mechanical Systems. San Francisco，CA，USA：MEMS 2014，2014：494—497.

[21]? SINGH T P，TSUCHIYA T，TABATA O. Frequency response of in-plane coupled resonators for investigating the acceleration sensitivity of MEMS tuning fork gyroscopes[J]. Microsystem Technologies，2012，18（6）：797—803.

[22]? 陳華，來強濤，郭桂良. 基于正交誤差對準的MEMS陀螺測控電路設計[J]. 湖南大學學報（自然科學版），2017，44（10）：114—121.

CHEN H，LAI Q T，GUO G L. Design of MEMS gyroscope interface circuit based on quadrature error alignment[J]. Journal of Hunan University （Natural Sciences），2017，44（10）：114—121.（In Chinese）

[23]? 趙鶴，王喆垚. 基于UKF的MEMS傳感器姿態(tài)測量系統(tǒng)[J]. 傳感技術學報，2011，24（5）：642—646.

ZHAO H，WANG Z Y. MEMS sensors based attitude measurement system using UKF[J]. Chinese Journal of Sensors and Actuators，2011，24（5）：642—646.（In Chinese）

[24]? 朱伯龍，屠成松，許哲明. 工程結(jié)構(gòu)抗震設計原理[M]. 上海：上海科學技術出版社，1982：173—184.

ZHU B L，TU C S，XU Z M. Seismic design principles of engineering structures[M]. Shanghai：Shanghai Science and Technology Press，1982：173—184.（In Chinese）

[25]? AISC-LRFD99 Load and resistance factor designspecification for structural steel buildings[S]. Chicago：American Institute of Steel Construction，1999：1—5.

[26]? 陳為真，汪秉文，胡曉婭. 基于時域積分的加速度信號處理[J].華中科技大學學報（自然科學版），2010，38（1）：1—4.

CHEN W Z，WANG B W，HU X Y. Acceleration signal processing by aumerical integration [J].Journal of Huazhong University of Science and Technology（Natural Sciences），2010，38（1）：1—4.（In Chinese）