亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于經(jīng)驗模態(tài)分解和模型縮聚的動力響應重構(gòu)方法研究

        2022-02-11 10:44:36鄒云峰付正億何旭輝盧玄東陽勁松
        工程力學 2022年2期
        關(guān)鍵詞:子結(jié)構(gòu)均方測點

        鄒云峰,付正億,何旭輝,盧玄東,陽勁松,周 帥

        (1.中南大學土木工程學院,湖南,長沙 410075;2.中南大學交通運輸工程學院,湖南,長沙 410075;3.中國建筑第五工程局有限公司,湖南,長沙 410004)

        隨著新技術(shù)、新材料在土木工程結(jié)構(gòu)中的應用,人們對土木工程結(jié)構(gòu)的安全性及可靠性提出更高要求。結(jié)合安裝在工程結(jié)構(gòu)上的傳感器采集所得數(shù)據(jù),可對工程結(jié)構(gòu)進行健康監(jiān)測及使用壽命評估[1-3],但結(jié)構(gòu)健康監(jiān)測及使用壽命評估需要解決的關(guān)鍵問題是合理地監(jiān)控關(guān)鍵位置的健康狀態(tài)[4-6],因此需要依賴于傳感器系統(tǒng)。而大多數(shù)工程結(jié)構(gòu)有限元模型自由度數(shù)目過多[7-10],此外,鑒于工程結(jié)構(gòu)的幾何復雜性及部件多樣性,某些位置(例如結(jié)構(gòu)交界面、狹縫等)不易安裝傳感器,而這些位置往往是結(jié)構(gòu)健康監(jiān)測及使用壽命評估的關(guān)鍵部位[11-13]。

        基于上述問題,近年來通過已安裝傳感器位置的響應對未安裝傳感器位置響應進行重構(gòu)的方法得到廣泛重視。Ribeiro 等[14]針對多自由度系統(tǒng)提出一種基于廣義傳遞矩陣的響應重構(gòu)法,在頻域內(nèi)建立兩組響應的傳遞矩陣,進而實現(xiàn)響應待測點的重構(gòu)工作。Law 等[15-17]基于頻響函數(shù)將上述位移傳遞率推廣至加速度傳遞率,完成了全模型或模型子結(jié)構(gòu)上的響應重構(gòu)。在進行子結(jié)構(gòu)基礎(chǔ)上的響應重構(gòu)時,將子結(jié)構(gòu)間的界面力視為外力并建立傳遞函數(shù),從而完成響應重構(gòu)。以上工作為頻域上的響應重構(gòu),而完成頻域基礎(chǔ)上的響應重構(gòu)工作時,需要多次運用傅里葉變換及逆變換,就計算效率而言,時域上的重構(gòu)方法更具優(yōu)越性。為此,Zhang 等[18-19]通過模態(tài)分析法,在時域基礎(chǔ)上建立兩組響應的模態(tài)轉(zhuǎn)換矩陣,并結(jié)合歸一化方法,完成了結(jié)構(gòu)應變、位移響應重構(gòu)。He 等[20-21]拓展了該重構(gòu)方法,提出一種基于經(jīng)驗模態(tài)分解(Empirical Mode Decomposition,EMD)的響應重構(gòu)方法,該方法通過間歇準則與經(jīng)驗模態(tài)分解法,得到測得響應的各單頻時域響應,結(jié)合模態(tài)傳遞矩陣完成各單頻響應的重構(gòu),再通過模態(tài)疊加法得到未安裝傳感器位置的動力響應。由于該方法不能完成應力及應變的重構(gòu),He 等[22]又基于歐拉伯努利梁理論,結(jié)合懸臂梁模擬及試驗,進一步實現(xiàn)了應力及應變的重構(gòu)。但現(xiàn)有的基于EMD 的時域響應重構(gòu)方法,需要獲取整個結(jié)構(gòu)的全局剛度矩陣與質(zhì)量矩陣。土木工程復雜結(jié)構(gòu)的自由度數(shù)目較大,在提取剛度矩陣、質(zhì)量矩陣時需要大量的計算資源與計算時間。同時,由于矩陣維數(shù)過大,也會導致計算傳遞矩陣時間過長,計算效率過低的問題,難以推廣應用于工程實際中的實時監(jiān)測[23-26]。

        針對以上問題,本文提出基于經(jīng)驗模態(tài)分解和模型縮聚的動力響應重構(gòu)法。不同于傳統(tǒng)的基于帶有間歇性準則的經(jīng)驗模態(tài)分解法,該響應重構(gòu)在超單元的基礎(chǔ)上完成。根據(jù)響應采集點及待測點所在位置劃分子結(jié)構(gòu),建立各子結(jié)構(gòu)的有限元模型,對于工程結(jié)構(gòu)而言,一般選取低階模態(tài)作為子結(jié)構(gòu)的主模態(tài),通過坐標變換,得到各子結(jié)構(gòu)的超單元模型[27-28]。然后,進行二次坐標變換,將各子結(jié)構(gòu)耦合為全模型,提取該模型的模態(tài)轉(zhuǎn)換矩陣。結(jié)合帶有間歇性準則的經(jīng)驗模態(tài)分解法,將采集點的響應轉(zhuǎn)換成各個單頻響應,通過模態(tài)轉(zhuǎn)換矩陣即可得到待重構(gòu)響應的各單頻響應,并通過模態(tài)疊加法即可完成超單元基礎(chǔ)上的重構(gòu),此時所得的重構(gòu)響應即為待測點的響應。

        1 基于經(jīng)驗模態(tài)分解和模型縮聚的動力響應重構(gòu)理論推導

        1.1 子結(jié)構(gòu)劃分

        一般情況下,根據(jù)響應采集位置和待測點所處位置劃分子結(jié)構(gòu),將這兩個測點所處位置分別作為各自子結(jié)構(gòu)的界面,由此形成子結(jié)構(gòu)。通過有限元法提取子結(jié)構(gòu)的剛度矩陣和質(zhì)量矩陣,便于后期計算。

        對于一般的子結(jié)構(gòu),其有限元基礎(chǔ)上的動力學運動方程可表示為:

        根據(jù)Craig-Bampton 模態(tài)綜合法,提取子結(jié)構(gòu)s的模態(tài)轉(zhuǎn)換矩陣 Φs,結(jié)合布爾矩陣L,通過兩次坐標變換,可得到該模型的超單元基礎(chǔ)上的動力學運動方程式(2)。各子結(jié)構(gòu)耦合為超單元模型后,各子結(jié)構(gòu)間的界面力gs(t)耦合為0。

        式中,各向量主要由兩部分構(gòu)成,各子結(jié)構(gòu)界面自由度集及子結(jié)構(gòu)內(nèi)部單元各自由度主模態(tài)集。皆為超單元模型的質(zhì)量矩陣、阻尼矩陣及剛度矩陣。模型的超單元位移響應可表示為:

        式中:p為超單元模型的位移響應集;pin為第n個子結(jié)構(gòu)的內(nèi)部單元自由度的主模態(tài)位移響應集;pjn為第n個子結(jié)構(gòu)的界面自由度的位移響應集。

        1.2 子結(jié)構(gòu)的動力響應重構(gòu)理論

        假設子結(jié)構(gòu)界面自由度m的響應Xm(t)已知,需要重構(gòu)另一子結(jié)構(gòu)邊界自由度r的響應Xr(t)。首先,結(jié)合坐標轉(zhuǎn)換將已知響應Xm(t)轉(zhuǎn)換為pm(t),由于邊界自由度為約束模態(tài),則Xm(t)等于pm(t)?;趲в虚g歇準則的EMD 方法,將pm(t)分解為各單頻的時域響應。其中,帶有間歇準則的經(jīng)驗模態(tài)分解法(EMD)通過以下步驟實現(xiàn)。

        首先,確定輸入信號pm(t)的極大值及極小值點;然后,由各極值點擬合上包絡線及下包絡線,計算上下包絡線的平均值m1(t);再由式h1(t)計算第一個分量,h1(t)可表示為:

        h1(t)=pm(t)-m1(t)(4)

        判斷h1(t)是否滿足IMF 條件,若否,則將h1(t)視為新的信號重新篩;建立信號h1(t)的包絡線,計算上下包絡線的平均值m11,得到分量h11(t),h11(t)可表示為:

        重復k次以上篩選過程,直到分量h1k為本征模態(tài)函數(shù)(IMF),h1k可表示為:

        設數(shù)據(jù)中篩選得到的第一個IMF 分量f1(t)=h1k,把剩余信號pm(t)-f1(t)重復以上篩選步驟獲得第二個IMFf2(t)。再將剩余信號不斷循環(huán),依次得到其他剩余的IMF,直到最后剩余信號r(t)為單調(diào)函數(shù),則停止篩選,剩余信號r(t)可表示為:

        原始信號pm(t)可表示為:

        通過EMD 分解得到的IMF 可能含有多個頻率成分,為了得到單頻的IMF,需要設置帶通濾波器。對于大型復雜結(jié)構(gòu),結(jié)構(gòu)動力響應重構(gòu)的階數(shù)一般較多,通過傅里葉變換,可得到每階模態(tài)的頻率ωi,由此確定帶通濾波器的濾波區(qū)間[ωiLωiH],其中ωiL<ωi<ωiH。將pm(t)通過濾波器后,再由EMD 分解得到各個單頻的IMF,則時域信號可表示為:

        式中:di(t) 為第i階模態(tài)響應;si(t)為其余非模態(tài)響應的IMF。

        界面自由度m的響應pm(t)及自由度r的響應pr(t)可表示為:

        超單元的r自由度響應與其對應的有限元響應關(guān)系可表示為:

        式中,Xr表示超單元的r自由度響應對應的有限元響應。

        具體重構(gòu)過程如圖1 所示。

        圖1 重構(gòu)流程圖Fig.1 Flowchart of the response reconstruction method

        2 數(shù)值模擬分析

        2.1 位移響應重構(gòu)

        通過懸臂梁仿真算例驗證本文提出的基于經(jīng)驗模態(tài)分解和模型縮聚的動力響應重構(gòu)法的有效性。如圖2 所示,該懸臂梁由鋁材制成,鋁材彈性模量為69 600 MPa,密度2730 kg/m3,長寬高分別為5 m、0.5 m 和0.05 m。采用ANSYS 建立有限元模型,單元類型為beam3 單元,共有40 個單元,40 個節(jié)點,120 個自由度。以自由端為位移響應采集點,跨中節(jié)點為位移響應待測點,如圖3 所示,以1 單元~20 單元為子結(jié)構(gòu)1,21 單元~40 單元為子結(jié)構(gòu)2。

        圖3 子結(jié)構(gòu)劃分示意圖Fig.3 The schematic diagram of the substructures

        該重構(gòu)方法與外界激勵無關(guān),將豎向隨機荷載均布于懸臂梁上,該隨機荷載由6 階低通巴特沃斯濾波器過濾的白噪聲模擬。每階模態(tài)阻尼設為1%,采樣頻率為1 kHz。為評估該重構(gòu)方法的有效性,將均方誤差(MSE)及相關(guān)系數(shù)(correlation coefficient)設為評估重構(gòu)效果的指標,其中MSE可表示為:

        式中,z(t)、k和分別表示待測點的重構(gòu)響應、矩陣z(t)的列數(shù)及待測點響應的理論值。

        如圖4(a)所示為模擬生成的60 s 響應采集點的位移響應,圖4(b)所示為該位移響應的頻譜圖,由該頻譜圖可知已識別得到的懸臂梁固有頻率為1.63 Hz 和10.38 Hz。根據(jù)該結(jié)構(gòu)的頻率可得表1所示的帶通濾波器。如圖4(c)和圖4(d)所示,結(jié)合帶有間歇準則的經(jīng)驗模態(tài)分解法(EMD),可將響應采集點的位移響應分解為多個單頻時域信號。

        圖4 響應采集點的位移Fig.4 The displacement measurement data of the response gauge location

        表1 帶通濾波器的頻率范圍Table 1 the frequency range of the band-pass filter

        由式(13)可得重構(gòu)所得的響應待測點的位移響應,如圖5 所示為響應理論值與重構(gòu)值的對比圖,其中圖5(a)表示60 s 時域的對比圖,為更清晰地反映重構(gòu)效果,如圖5(b)所示,任取其中10 s的對比結(jié)果,由圖可知,重構(gòu)所得的響應與理論值吻合較好。且響應采集點的理論值與重構(gòu)值的相關(guān)系數(shù)為0.986(接近1),均方誤差(MSE)為1.28×10-5遠小于3%),以上表明了該重構(gòu)方法的有效性及精確性。就計算效率而言,結(jié)合模型縮聚,懸臂梁有限元模型轉(zhuǎn)換為超單元模型,剛度矩陣和質(zhì)量矩陣大小由原來的120×120 降低至10×10,由此,計算量減小,計算效率得到提高。基于有限元基礎(chǔ)上的經(jīng)驗模態(tài)分解法完成該懸臂梁響應重構(gòu)所需時間為80.42 s,而基于經(jīng)驗模態(tài)分解和模型縮聚的動力響應重構(gòu)方法完成該懸臂梁響應重構(gòu)所需時間為6.89 s,由此可知,計算效率得以提高。

        圖5 響應采集點的理論值與重構(gòu)值對比圖Fig.5 The comparison of the theoretical results and the reconstructed results

        2.2 噪聲對重構(gòu)精度的影響

        結(jié)合圖2 所示的懸臂梁數(shù)值模擬算例分析四個不同噪聲等級(噪聲等級分別為0%、5%、10%、20%)對懸臂梁響應重構(gòu)精度的影響。為模擬傳感器的測量噪聲,生成位移響應后,將噪聲單元加入各位移響應。其中,噪聲單元指高斯脈沖過程中的加速度的均方根占最大均方根百分比。通過噪聲等級來表征噪聲單元,如:10%的噪聲等級表示該噪聲為高斯脈沖過程中的加速度的均方根占最大均方根百分比為10%?;谠撝貥?gòu)方法,通過懸臂梁自由端的響應重構(gòu)懸臂梁跨中節(jié)點的響應。如圖6 所示,為四個不同噪聲等級下待測點的響應重構(gòu)值與理論值對比圖。定義兩個指標以衡量噪聲等級對重構(gòu)精度的影響,即均方誤差和相關(guān)系數(shù)。表2、表3 分別表示不同噪聲等級下的響應重構(gòu)均方誤差和重構(gòu)響應的理論值與重構(gòu)值的相關(guān)系數(shù)。圖7、圖8 分別為重構(gòu)響應的理論值與重構(gòu)值的均方誤差、相關(guān)系數(shù)隨噪聲等級的變化曲線圖。

        表2 不同噪聲等級下的響應重構(gòu)均方誤差Table 2 Mean square errors of the four cases

        表3 不同噪聲等級下的響應重構(gòu)相關(guān)系數(shù)Table 3 Correlation coefficients the four cases

        圖6 四個不同噪聲等級下待測點的響應重構(gòu)結(jié)果Fig.6 The reconstructed results at four different noise levels

        圖7 不同噪聲等級下的均方誤差Fig.7 Reconstruction performance measured in MSE under different levels

        圖8 不同噪聲等級下的相關(guān)系數(shù)Fig.8 Reconstruction performance measured in correlation coefficient under different levels

        由不同噪聲等級下的重構(gòu)對比圖、MSE(皆小于3%)及相關(guān)系數(shù)(皆接近于1)可知,不同噪聲等級下,響應重構(gòu)的精度皆能得到保證。但隨著噪聲等級增加,MSE 隨之增大,相關(guān)系數(shù)隨之減小,EMD 的邊界效應越來越顯著,由此可知,噪聲等級對重構(gòu)精確度存在一定的影響。

        2.3 主模態(tài)數(shù)量對重構(gòu)精度的影響

        結(jié)合圖2 所示的懸臂梁數(shù)值模擬算例,分析不同主模態(tài)數(shù)量對懸臂梁響應重構(gòu)精度的影響,態(tài))下的重構(gòu)效果,以研究主模態(tài)數(shù)量對重構(gòu)精度的影響。結(jié)合該重構(gòu)方法,通過懸臂梁自由端的響應重構(gòu)懸臂梁跨中節(jié)點的響應。如圖9 所示,為4 個工況下待測點的響應重構(gòu)值與理論值對比圖,仍以均方誤差(MSE)及相關(guān)系數(shù)來衡量主模態(tài)數(shù)量對重構(gòu)精度的影響。如表4、表5 所示,分別為不同主模態(tài)數(shù)量下響應重構(gòu)的均方誤差(MSE)和相關(guān)系數(shù)。由圖9、表4 及表5 可知,4 種工況的均方誤差(MSE)均較小(小于3%)及相關(guān)系數(shù)討論4 種工況(工況1:1 階主模態(tài);工況2:2 階主模態(tài);工況3:3 階主模態(tài);工況4:4 階主模均接近1。由此可得,主模態(tài)的數(shù)量對該重構(gòu)方法的精度影響較小,主要是由于子結(jié)構(gòu)界面的模態(tài)轉(zhuǎn)換矩陣由約束模態(tài)計算所得,與主模態(tài)的相關(guān)參數(shù)并無太大關(guān)系。

        表4 各工況下的響應重構(gòu)均方誤差Table 4 Mean square errors of the four cases

        表5 各工況下的響應重構(gòu)相關(guān)系數(shù)Table 5 Correlation coefficient of the four cases

        圖9 各工況下待測點的響應重構(gòu)Fig.9 The comparison of the theoretical and reconstructed results in the different conditions

        3 結(jié)論

        本文的研究工作將有限元基礎(chǔ)上的響應重構(gòu)拓展到超單元,將模態(tài)綜合法和帶有間歇準則的經(jīng)驗模態(tài)分解法結(jié)合以解決大型復雜工程結(jié)構(gòu)的重構(gòu)問題。得到的主要結(jié)論如下:

        (1)該重構(gòu)方法無需考慮各子結(jié)構(gòu)邊界條件,有效地降低了有限元模型相關(guān)參數(shù)數(shù)學矩陣的維度,可以更好地適應大型復雜工程結(jié)構(gòu)的動力響應重構(gòu),且較大地提高了響應重構(gòu)的計算效率。

        (2)基于懸臂梁仿真算例,通過該響應重構(gòu)方法得到的待測點響應重構(gòu)值與理論值較接近,其相關(guān)系數(shù)均接近于1 且均方誤差(MSE)小于3%,由此驗證了該方法的有效性。

        (3)噪聲及主模態(tài)數(shù)量影響研究表明,隨著噪聲等級的提高,該重構(gòu)方法的精度能得到保證,但仍對重構(gòu)精度有一定影響;就主模態(tài)數(shù)量對重構(gòu)精度影響而言,改變主模態(tài)數(shù)量對重構(gòu)精度影響不大,由于子結(jié)構(gòu)界面的模態(tài)轉(zhuǎn)換矩陣由約束模態(tài)計算得到,與主模態(tài)相關(guān)參數(shù)并無太大關(guān)系。

        猜你喜歡
        子結(jié)構(gòu)均方測點
        液壓支架整機靜強度試驗及等效應力分析
        一類隨機積分微分方程的均方漸近概周期解
        完全對換網(wǎng)絡的結(jié)構(gòu)連通度和子結(jié)構(gòu)連通度
        基于CATIA的汽車測點批量開發(fā)的研究與應用
        Beidou, le système de navigation par satellite compatible et interopérable
        鋼框架腹板雙角鋼連接梁柱子結(jié)構(gòu)抗倒塌性能分析
        基于子結(jié)構(gòu)的柴油機曲軸有限元建模方法研究
        基于抗差最小均方估計的輸電線路參數(shù)辨識
        電力建設(2015年2期)2015-07-12 14:15:59
        基于隨機牽制控制的復雜網(wǎng)絡均方簇同步
        拱壩結(jié)構(gòu)損傷的多測點R/S分析
        久久精品中文字幕第23页| 国产极品粉嫩福利姬萌白酱| 亚洲一区二区懂色av| 久久久久人妻精品一区二区三区| 欧美黑人群一交| 香蕉视频一级| 视频一区中文字幕亚洲| 在线观看亚洲av每日更新影片| 国产精品兄妹在线观看麻豆| 亚洲av熟妇高潮30p| 精品久久久无码不卡| 人妖国产视频一区二区| 色欲aⅴ亚洲情无码av| 亚洲av鲁丝一区二区三区| 2020最新国产激情| 日韩免费精品在线观看| 久久精品国产久精国产果冻传媒| 亚洲精品国产美女久久久| 久久综合视频网站| 精品一区二区三区老熟女少妇| 爽爽影院免费观看| 国产又滑又嫩又白| 久久亚洲国产成人精品v| 男男做h嗯啊高潮涩涩| 国产不卡视频一区二区三区| 亚洲暴爽av天天爽日日碰| 精品国产福利一区二区三区| 日本免费一区二区三区在线播放| 免费看黑人男阳茎进女阳道视频| 五月激情婷婷丁香| 蜜桃av福利精品小视频| 亚洲熟妇av一区二区三区| 看国产黄大片在线观看| 大屁股少妇一区二区无码| 久久热免费最新精品视频网站| 久久婷婷人人澡人人喊人人爽| 欧美精品一级| 国产精品中文字幕日韩精品| 亚洲av不卡免费在线| 55夜色66夜色国产精品视频| 2022国内精品免费福利视频|