基于陀螺儀轉(zhuǎn)角傳感器的動(dòng)態(tài)信號(hào)測量及物理參數(shù)時(shí)域識(shí)別

周云 曾雅麗思 趙瑜 易偉建

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

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

中圖分類號(hào)：TU317? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：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)識(shí)別方法始于20世紀(jì)70年代，引入土木工程領(lǐng)域已有近40余年的發(fā)展歷史[1].該方法是工程結(jié)構(gòu)損傷診斷和狀態(tài)評(píng)估的重要手段，它可以從建筑物動(dòng)力響應(yīng)中識(shí)別出結(jié)構(gòu)的當(dāng)前狀態(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]. 進(jìn)行結(jié)構(gòu)物理參數(shù)識(shí)別的方法，主要分為頻域識(shí)別和時(shí)域識(shí)別兩種，與頻域識(shí)別方法相比，時(shí)域識(shí)別更加直接，實(shí)測的時(shí)程數(shù)據(jù)中的結(jié)構(gòu)信息更加豐富，且增加時(shí)程響應(yīng)采樣點(diǎn)數(shù)可以有效增加時(shí)域識(shí)別方程的個(gè)數(shù)[3]. 目前制約時(shí)域識(shí)別應(yīng)用于實(shí)際工程中的一個(gè)主要原因是結(jié)構(gòu)響應(yīng)信息不完備，通常時(shí)域方程中主要使用位移類信息，轉(zhuǎn)角類信息很難測試得到，往往只能測量結(jié)構(gòu)的一個(gè)狀態(tài)量，加速度、速度或位移響應(yīng)的一種，而且轉(zhuǎn)角信息難以準(zhǔn)確測量.如何在測量信息有限的情況下識(shí)別參數(shù)以及如何進(jìn)行信息不完備情況下的識(shí)別問題是研究重點(diǎn)之一[4].

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

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

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

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

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

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

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

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

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

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

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

3? ?基于最小二乘遞推算法的結(jié)構(gòu)物理參數(shù)時(shí)域識(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ù)識(shí)別已經(jīng)比較成熟，本文研究對象為一般框架結(jié)構(gòu)，在水平荷載作用下，其框架節(jié)點(diǎn)將發(fā)生轉(zhuǎn)動(dòng)，因此需要考慮框架節(jié)點(diǎn)轉(zhuǎn)動(dòng)的影響.

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

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

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

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

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

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

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

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

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

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

4? ?算例分析

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

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

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

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

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

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

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

5? ?試驗(yàn)研究

5.1? ?MEMS陀螺儀傳感器的角度測量精度驗(yàn)證

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

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

5.2? ?MEMS陀螺儀傳感器的時(shí)域識(shí)別

5.2.1? ?試驗(yàn)?zāi)Ｐ?/p>

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

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

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

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

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

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

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

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

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

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

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

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

6? ?結(jié)? ?論

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

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

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

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

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

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

參考文獻(xià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]? ? 李國強(qiáng)，李杰. 工程結(jié)構(gòu)動(dòng)力檢測理論與應(yīng)用[M]. 北京：科學(xué)出版社，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]? ? 謝獻(xiàn)忠. 結(jié)構(gòu)動(dòng)力學(xué)系統(tǒng)時(shí)域辨識(shí)理論與試驗(yàn)研究[D]. 長沙：湖南大學(xué)土木工程學(xué)院，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)時(shí)域系統(tǒng)識(shí)別研究[D]. 上海：同濟(jì)大學(xué)土木工程學(xué)院，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ù)識(shí)別中的應(yīng)用[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），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]? ? 陳雋. 高層建筑損傷檢測中的復(fù)合反演理論與試驗(yàn)研究[D]. 上海：同濟(jì)大學(xué)建筑工程系，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ù)識(shí)別方法研究[J]. 工程力學(xué)，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)定平臺(tái)中陀螺技術(shù)的發(fā)展現(xiàn)狀和應(yīng)用研究[J]. 飛航導(dǎo)彈，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]? 陳華，來強(qiáng)濤，郭桂良. 基于正交誤差對準(zhǔn)的MEMS陀螺測控電路設(shè)計(jì)[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版），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]. 傳感技術(shù)學(xué)報(bào)，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)抗震設(shè)計(jì)原理[M]. 上海：上?？茖W(xué)技術(shù)出版社，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]? 陳為真，汪秉文，胡曉婭. 基于時(shí)域積分的加速度信號(hào)處理[J].華中科技大學(xué)學(xué)報(bào)（自然科學(xué)版），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）