• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    利用改進交叉模型交叉模態(tài)的隨機模型修正方法

    2023-07-10 18:17:59王炎陳輝黃斌柴滿
    振動工程學報 2023年2期
    關鍵詞:交叉修正模態(tài)

    王炎 陳輝 黃斌 柴滿

    摘要 將混合攝動?伽遼金方法和改進的交叉模型交叉模態(tài)技術相結合,提出了一種隨機模型修正方法。該方法有效緩解了模型修正過程中測量數據有限和測量誤差不確定的影響??紤]到實測模態(tài)數據具有不確定性,基于改進的交叉模型交叉模態(tài)方法,建立了一個新的描述結構隨機參數和隨機響應關系的模型修正方程。利用混合攝動?伽遼金方法求解該隨機修正方程,進而得到結構隨機修正參數的統(tǒng)計特征。簡支梁的數值結果表明,該方法在測量數據不確定性較大時仍能保持很高的修正精度,同時計算效率比蒙特卡羅模擬方法高出一個數量級。在測量模態(tài)數據較少的情況下,該方法比單獨的混合攝動?伽遼金修正方法修正效果好,且比交叉模型交叉模態(tài)法的修正精度更高。框架試驗的結果表明,該方法可以同時修正結構的剛度和質量,修正后的結構參數與預設工況基本吻合,同時能復現結構的測量模態(tài),從而驗證了所提方法的有效性。

    關鍵詞 隨機模型修正; 隨機混合?攝動伽遼金方法; 改進的交叉模型交叉模態(tài)方法

    引 言

    近幾十年來,基于動力測量數據的有限元模型修正方法越來越受到關注。許多研究人員在這一領域進行了廣泛的研究,并取得了大量研究成果[1?4]。

    在動力有限元模型修正中,修正參數的選擇對修正結果有很大影響。如果修正參數過多,在修正過程中往往會出現病態(tài)問題,所以在修正模型之前首先要排除不敏感的修正參數[5]。關于動力有限元模型修正方法,Hu等[6]提出了一種基于交叉模型交叉模態(tài)(CMCM)方法的模型修正技術。與傳統(tǒng)的模型修正方法不同,該方法可以同時修正結構的剛度矩陣、質量矩陣和阻尼矩陣。此外,該方法不用迭代計算,計算效率較高。在CMCM方法中,通過將結構的實測模態(tài)和計算模態(tài)相乘,就可以僅用少量的測量模態(tài)構建多個模型修正方程。李世龍等[7]利用CMCM方法,提出了一種有效識別子結構邊界約束狀態(tài)的模型修正方法。Wang等[8]使用了CMCM方法對海上平臺進行了試驗研究,證明了當結構的實際測量模態(tài)不完整且只有低階測量模態(tài)可用時CMCM方法的有效性。在已有的CMCM方法的基礎上,Liu等[9]提出了一種基于改進的交叉模型交叉模態(tài)(ICMCM)的模型修正方法,該方法充分利用實測數據,進一步增加了修正方程的個數。然而,這些方法僅涉及確定性有限元模型修正,當結構參數的不確定性或者測量噪聲無法避免時,現有的CMCM方法將不適用。因此,充分利用CMCM方法的優(yōu)點,并將它融入隨機模型修正中,是一項非常有意義的工作。

    在隨機模型修正領域中,蒙特卡羅方法、攝動法以及貝葉斯方法被廣泛使用。Schu?ller等[10]使用了具有大樣本的蒙特卡羅模擬來計算模型修正的統(tǒng)計特性。宗周紅等[11]在對下白石連續(xù)剛構橋進行模型修正的過程中,利用蒙特卡羅模擬方法和有限元方法進行不確定性量化分析,并評價模型的預測精度,實現對于連續(xù)剛構橋的有限元模型確認。但是對于大型結構而言,這種方法的計算效率太低,耗時過長。與蒙特卡羅方法不同的是,攝動法具有推導簡單、計算效率高的特點。Hua等[12]使用一種改進的攝動法,利用隨機實測模態(tài)數據對桁架橋有限元模型進行修正,并估計了結構參數的均值和均方差。盡管攝動法的計算效率比較高,但其對測量誤差的變異性要求比較苛刻。隨機模型修正方法中,另一個具有代表性的方法是基于馬爾可夫鏈與蒙特卡羅抽樣的貝葉斯方法[13?15],但是基于此種抽樣的貝葉斯方法會面臨較大的挑戰(zhàn),即需要非常耗時的重復有限元計算。為了提高計算效率,Wan等[16]和Fang等[17]分別采用了高斯代理模型和隨機響應面模型對原有的貝葉斯方法進行了改進。與上述方法不同,Huang等[18]提出了一種基于混合攝動?伽遼金方法(HPG)的隨機模型修正方法(HPG?SMUM),該方法在測量變異性較大情況下具有比較高的修正精度和效率,此方法也為確定性模型修正方法擴展到隨機領域提出了一個新的思路和完整的框架。

    本文將隨機攝動?伽遼金方法與改進的交叉模型交叉模態(tài)方法結合,提出一種隨機模型修正方法。該方法可利用含測量誤差的少量模態(tài)測量數據實現結構有限元模型的有效修正。文中用一個簡支梁的數值算例來驗證該方法的有效性和不同模態(tài)組合的穩(wěn)定性,并利用七層框架的模態(tài)試驗來驗證所提方法在較少測量模態(tài)情況下仍能同時有效地修正結構剛度和質量。

    1 基于ICMCM的隨機模型修正方程

    考慮具有N個自由度的無阻尼結構,該結構初始模型滿足以下特征值方程:

    式中 Ka和Ma分別為初始結構模型的整體剛度矩陣和質量矩陣;λi和?i分別為初始模型的第i階特征值和特征向量;nc為初始模型的計算模態(tài)的個數。

    類似地,實際結構的特征值方程可以表示為:

    式中 Kd和Md分別為實際結構模型的整體剛度矩陣和質量矩陣;λ?j和??j分別為實際模型的第j階特征值和特征向量;nm為實際模型的計算模態(tài)的個數。

    初始結構與實際結構的質量矩陣、剛度矩陣存在以下關系:

    式中 Ne為結構的單元個數;Kn和Mn分別為結構第n個單元的N×N單元組裝矩陣;αn和βn分別為結構第n個單元的剛度和質量的修正系數,表示實際結構的單元剛度和質量相對于初始矩陣的變化率。

    通過文獻[6]可以得到確定性的基于交叉模型交叉模態(tài)的模型修正方程為:

    對式(5)進行因式變換可以得到:

    通過求解式(6)所示的方程組可以得到剛度和質量的修正系數αn和βn。但是由于在實際的模態(tài)測量中只能精確測量出前幾階模態(tài),使得修正系數方程組的方程個數比較少,導致求解結果不正確且不穩(wěn)定。因此Liu等[9]對傳統(tǒng)的CMCM方法進行改進,充分利用測量模態(tài)數據,在式(2)方程兩邊同時乘??Tj,得到如下所示的基于ICMCM的模型修正方程:

    對式(7)進行因式變換,可以得到:

    顯然,式(8)也含有與方程(6)相同的修正系數,結合式(6)和(8),就可以得到更多的修正方程,確保修正方程的適定性。

    在實際結構的模態(tài)試驗過程中,不可避免地會遇到測量誤差。假定第j階的特征值和特征向量可以表示為:

    式中 λ?0j和??0j分別為測量的第j階特征值和特征向量均值;λ?1j和??1j分別為第j階測量誤差的確定性部分;ξj為與測量誤差相關的隨機變量,且隨機變量的分布類型由實測數據的統(tǒng)計特征或者是工程經驗決定。

    假設所有隨機變量ξj完全相關,并且表示為隨機變量ξ,則第n個剛度單元和質量單元的修正系數可以分別用下式表示:

    上述方法就是本文所提出的結合HPG和ICMCM的隨機模型修正方法(HPG?ICMCM)。假設用[C(0)E(0)]和[C(1)E(1)]代替式(19),(22)和(25)中的[C(0)IE(0)I]和[C(1)IE(1)I],相應的向量γ(0),γ(1),γ(2)也可以通過上述方式遞推得到。此時,HPG?ICMCM方法退化為HPG?CMCM方法。

    需要注意的是,實際結構的轉角模態(tài)往往難以測量。此外,由于測量條件的限制,僅能測量包括部分測點的振型。因此,本文使用文獻[19]的模態(tài)擴階方法得到完備振型的均值和標準差。同時,在求解方程組(19),(22)和(25)的過程中,采用截斷奇異值分解或者L1正則化技術[20]避免方程病態(tài)的問題。

    3 數值算例

    考慮一個簡支梁,如圖2所示。簡支梁的跨度為6 m,截面為0.2 m×0.25 m,彈性模量為2.8×1010 Pa,密度為2.5×103 kg/m3。將該Euler?Bernoulli梁的有限元模型沿梁長度方向劃分為15個相同的單元,每個節(jié)點包含豎向位移和轉角兩個自由度。

    根據測量經驗,可以假設實測模態(tài)數據服從某種概率分布,如正態(tài)分布或者β分布。由于實測數據是有界的,本文假設實測模態(tài)數據服從β分布。依據工程經驗,動力特性測試數據的變異系數一般在0.01~0.02之間,這里假設變異系數為0.02。

    首先,考慮結構質量不發(fā)生變化,單元1,3,7,9和11的剛度分別減小30%,15%,20%,20%和30%,其余單元的剛度和初始模型相同。取前六階初始模型的計算模態(tài)和前六階測量模態(tài),使用HPG?ICMCM方法和HPG?CMCM方法對模型進行修正,同時利用與本文所提出的方法對應的蒙特卡羅模擬方法(MC?ICMCM)和Huang等[18]的HPG?SMUM方法求解上述方法中的修正系數的統(tǒng)計特性,在求解過程中使用奇異值分解正則化技術降低矩陣求逆的不適定性,以提高計算精度。修正結果如圖3和4所示。

    觀察圖3和4,不難看出,當結構的自由度比較多但測量模態(tài)有限時,通過HPG?SMUM方法得到的修正系數均值和預設的真實值差別比較大。例如,HPG?SMUM方法得到的單元1剛度修正系數均值為-0.02,和MC?ICMCM方法結果相比,相對誤差接近90%。同時,單元2,15的剛度修正系數均小于-0.1,出現了明顯誤判。對于單元8,10,13和15,HPG?SMUM方法的修正系數標準差結果和MC?ICMCM方法結果最大相對誤差達到400%,說明在這種情況下HPG?SMUM方法修正效果不能令人滿意。而通過HPG?ICMCM方法和HPG?CMCM方法得到的各單元修正系數與MC?ICMCM相比較,均值的絕對誤差均未超過0.03,標準差的相對誤差基本小于20%??梢哉f明統(tǒng)計結果和仿真試驗預設的結果基本吻合,并且HPG?ICMCM方法的均值結果吻合更好。

    為了驗證同時修正質量和剛度時本文方法的有效性,假設單元3,5,6,8,9,11和13的實際質量分別增加10%,20%,20%,20%,20%,20%和10%,同時,單元1,3,5,7,9,11和13的彈性模量分別降低30%,15%,20%,20%,20%,30%和30%,其余單元的質量和剛度和初始模型相同。選擇這15個單元的質量和彈性模量作為待修正的參數。首先假設測量得到了被測結構的前六階模態(tài)的頻率和豎向位移振型,再通過模態(tài)擴階方法得到被測模態(tài)的完整形式。之后對于初始模型,通過計算得到其前七階模態(tài)數據。這里分別使用MC?ICMCM,HPG?ICMCM和HPG?CMCM三種方法進行模型修正。修正系數的統(tǒng)計特性如圖5~8所示。

    從圖5~8中可以看出,一方面,在剛度和質量出現變化的單元里,由HPG?ICMCM方法得到的剛度和質量修正系數均值與MC?ICMCM方法得到結果的相對誤差均小于10%。并且,除了單元3之外,各單元剛度與質量修正系數標準差與MC?ICMCM方法相比均小于30%,這個現象說明所提出的HPG?ICMCM方法的修正精度和效果是令人滿意的。另一方面,HPG?CMCM方法的修正系數均值和預設的值相差較大,特別是在修正剛度時,除單元11和13以外,均出現了明顯的誤判。從而可以說明HPG?ICMCM方法得到的修正系數的統(tǒng)計特性比HPG?CMCM方法更加接近假定的真值,并且與MC?ICMCM方法得到的結果非常接近。除此之外,為了分析測量誤差變異系數對HPG?ICMCM方法和HPG?CMCM方法的影響,圖9給出了變異系數為0.02時,修正后結構的前五階頻率的概率密度函數。

    從圖9中可以看出,HPG?ICMCM方法與蒙特卡羅模擬方法的結果吻合,而通過HPG?CMCM方法得到的修正頻率不符合仿真預設的實測頻率。這進一步說明了HPG?ICMCM方法的優(yōu)越性。此外,基于三萬個樣本的MC?ICMCM方法在CPU為i5?10400、運行內存16 GB的個人計算機上計算時間超過了1800 s,而本文提出的HPG?ICMCM方法僅用時120 s,二者對比說明了此方法的高效率。

    對于不同的模態(tài)組合,文獻[9]指出當測量模態(tài)數據較少時,確定性CMCM方法可能會導致修正結果不穩(wěn)定。接下來,將驗證在有限實測數據的情況下,HPG?ICMCM方法的穩(wěn)定性。假設質量和剛度的折減量與之前簡支梁仿真算例的預設值完全相同,不進行模態(tài)擴階,考慮計算和測量模態(tài)的不同組合工況。不同工況下,模態(tài)組合如表1所示。

    在這四種工況下,采用HPG?ICMCM方法對簡支梁進行模型修正,得到單元修正系數的統(tǒng)計特性,如圖10~13所示。圖10~13結果顯示在四種不同的模態(tài)組合中,獲得的修正系數統(tǒng)計特性非常接近。

    從圖14中可以看到,選取不同的模態(tài)組合都可以得到較準確的修正結果。由于采用了ICMCM方法增加了修正方程數量,盡管測量模態(tài)的數量逐漸減少,修正后的頻率仍然能很好地與測量結果吻合,說明了本文提出方法的穩(wěn)定性。

    4 七層框架試驗

    為了驗證HPG?ICMCM方法的有效性,制作了一個七層框架,如圖15(a)所示。該七層框架層高為150 mm,框架動力模型采用葫蘆串模型,如圖15(b)所示,各單元質量為每層鋁合金質量塊及低頻傳感器和夾具組成。層間剛度由兩側的鋼板提供,兩側側板均采用1 mm厚的304不銹鋼板切割成型制作??蚣艿牡撞渴褂寐萁z緊固在試驗臺上。層間鋼板材料的彈性模量為194 GPa、密度為7.93 g/cm3,泊松比為0.3。每層側板的寬度為100 mm,在框架的模態(tài)試驗中,將第2,4和6層間兩側的鋼板分別切除30%,10%和20%,用來模擬剛度退化??蚣艿母鲉卧馁|量如表2所示,在單元2,4和5處附加質量塊模擬質量變化。在模態(tài)測試中,使用5個加速度傳感器分兩批測量。由于傳感器的重量不能忽略,為了使測試過程中每層質量相同,因此在沒有布置傳感器的層中布置與傳感器等重的配重塊。試驗中,采取了6種不同的傳感器布置方式進行了6組測量。

    每一組測量均采用不測力法對框架結構進行模態(tài)測試。在采集了7個測量點的加速度數據之后,使用增強型頻域分解方法[21]識別該框架的模態(tài),并采用測量軟件內5種不同的分析點數(512,1024,2048,4096,8192)進行模態(tài)分析。對30組樣本進行統(tǒng)計分析之后,得到前三階測量模態(tài)的均值,并且得出測量頻率的變異系數為0.01。預計在實際工程測量中變異系數會更大。

    選擇前三階實測模態(tài)和初始模型的前四階計算模態(tài)用于模型修正,將7個單元的彈性模量和質量作為修正系數,總共14個修正參數。其中,七層框架的剛度修正系數從下到上編號為1~7,每層對應的質量修正系數編號為8~14。采用HPG?ICMCM方法進行計算,并使用L1正則化技術降低求解過程中矩陣求解的不適定性,得到修正系數的統(tǒng)計特性如圖16和17所示。

    從圖16中可以看出,修正參數的均值與預設工況基本吻合。由于測量誤差的隨機性,修正后的參數也具有隨機性,修正系數的標準差如圖17所示。從圖17中可以看出,修正系數的標準差最大值為0.03,最小值為0.005。用修正后的參數計算結構頻率的概率密度函數,如圖18所示。從圖18中可以看出,本文方法修正的結構頻率概率密度與測量結果基本一致。這說明了本文方法對于試驗框架結構是有效的。

    5 結 論

    本文提出了一種交叉模型交叉模態(tài)隨機有限元模型修正方法。該方法成功地將確定性的改進交叉模型交叉模態(tài)模型修正方法拓展到隨機領域。建立了基于ICMCM方法的隨機模型修正方程,并對方程進行了求解。該方法同時具備了ICMCM方法僅用少量模態(tài)即可構造大量修正方程的優(yōu)點,以及能夠考慮測量誤差的隨機性,并能用混合攝動?伽遼金方法高效求解隨機模型修正方程的優(yōu)點。

    簡支梁算例的結果表明,本文方法可以有效處理測量數據中較大的不確定性,并且計算效率要比直接采用蒙特卡羅模擬方法高出1個數量級。當測量數據較少時,新的方法比已有的混合攝動?伽遼金修正方法修正效果好,且比交叉模型交叉模態(tài)法的修正精度更高。七層框架結構試驗表明了本文方法對實際結構模型修正的有效性。

    參考文獻

    1Mottershead J E, Friswell M I. Model updating in structural dynamics: a survey[J]. Journal of Sound & Vibration, 1993, 167(2): 347-375.

    2李輝, 丁樺. 結構動力模型修正方法研究進展[J]. 力學進展, 2005, 35(2): 170-180.

    LI Hui, DING Hua. Progress in model updating for structural dynamics[J]. Advances in Mechanics, 2005, 35(2): 170-180.

    3姜東, 費慶國, 吳邵慶. 基于區(qū)間分析的不確定性結構動力學模型修正方法[J]. 振動工程學報, 2015, 28(3): 352-358.

    JIANG Dong, FEI Qingguo, WU Shaoqing. Updating of structural dynamics model with uncertainty based on interval analysis[J]. Journal of Vibration Engineering, 2015, 28(3): 352-358.

    4Wan H P, Ren W X. A residual-based Gaussian process model framework for finite element model updating[J]. Computers and Structures, 2015, 156: 149-159.

    5Wan H P, Ren W X. Parameter selection in finite element model updating by global sensitivity analysis using Gaussian process metamodel[J]. Journal of Structural Engineering, 2015, 141(6): 04014164.

    6Hu S L J, Li Huajun, Wang Shuqing. Cross-model cross-mode method for model updating[J]. Mechanical Systems & Signal Processing, 2007, 21(4): 1690-1703.

    7李世龍, 馬立元, 李永軍, 等.一種新的子結構邊界約束模型修正方法及其應用[J]. 振動工程學報, 2015, 28(5): 730-740.

    LI Shilong, MA Liyuan, LI Yongjun, et al. A method for model updating of substructure boundary constraints and its application[J]. Journal of Vibration Engineering, 2015, 28(5): 730-740.

    8Wang Shuqing, Li Yingchao, Li Huajun. Structural model updating of an offshore platform using the cross model cross mode method: an experimental study[J]. Ocean Engineering, 2015, 97: 57-64.

    9Liu K, Yan R J, Guedes Soares C. An improved model updating technique based on modal data[J]. Ocean Engineering, 2018, 154: 277-287.

    10Schu?ller G I, Calvi A, Pellissetti M F, et al. Uncertainty analysis of a large-scale satellite finite element model[J]. Journal of Spacecraft and Rockets, 2009, 46(1): 191-202.

    11宗周紅, 高銘霖, 夏樟華. 基于健康監(jiān)測的連續(xù)剛構橋有限元模型確認(Ⅱ)——不確定性分析與模型精度評價[J]. 土木工程學報, 2011, 44(3): 85-92.

    ZONG Zhouhong, GAO Minglin, XIA Zhanghua. Finite element model validation of the continuous rigid frame bridge based on structural health monitoring—part Ⅱ: uncertainty analysis and evaluation of model accuracy[J]. China Civil Engineering Journal, 2011, 44(3): 85-92.

    12Hua X G, Ni Y Q, Chen Z Q, et al. An improved perturbation method for stochastic finite element model updating[J]. International Journal for Numerical Methods in Engineering, 2008, 73(13): 1845-1864.

    13Beck J L, Katafygiotis L S. Updating models and their uncertainties. I: Bayesian statistical framework[J]. Journal of Engineering Mechanics, 1998, 124(4): 455-461.

    14Katafygiotis L S, Beck J L. Updating models and their uncertainties. II: model identifiability[J]. Journal of Engineering Mechanics, 1998, 124(4): 463-467.

    15Lam H F, Yang J H, Au S K. Markov chain Monte Carlo-based Bayesian method for structural model updating and damage detection[J]. Structural Control and Health Monitoring, 2018, 25(4): e2140.

    16Wan H P, Ren W X. Stochastic model updating utilizing Bayesian approach and Gaussian process model[J]. Mechanical Systems and Signal Processing, 2016,70-71: 245-268.

    17Fang S E, Chen S, Lin Y Q, et al. Probabilistic damage identification incorporating approximate Bayesian computation with stochastic response surface[J]. Mechanical Systems and Signal Process, 2019, 128: 229-243.

    18Huang B, Chen H. A new approach for stochastic model updating using the hybrid perturbation-Galerkin method[J]. Mechanical Systems and Signal Process, 2019, 129: 1-19.

    19張德文. 改進Guyan~遞推減縮技術[J]. 計算結構力學及其應用, 1996, 13(1): 90-94.

    ZHANG Dewen. An improved Guyan reduction and successive reduction procedure of dynamic model[J]. Chinese Journal of Computational Mechanics, 1996, 13(1): 90-94.

    20Zhou X Q, Xia Y, Weng S. L1 regularization approach to structural damage detection using frequency data[J]. Structural Health Monitoring, 2015, 14(6): 571-582.

    21Brincker R, Zhang L, Andersen P. Modal identification from ambient responses using frequency domain decomposition[A]. Proceedings of IMAC?ⅩⅧ: A Conference on Structural Dynamics[C]. USA: Society for Experimental Mechanics, 2000: 625-630.

    Stochastic model updating method using the improved cross-model cross-mode technique

    WANG Yan 1 ?CHEN Hui 1,2HUANG Bin 1 ?CHAI Man 1

    1. School of Civil Engineering and Architecture, Wuhan University of Technology, Wuhan 430070, China;

    2. College of Post and Telecommunication, Wuhan Institute of Technology, Wuhan 430073, China

    Abstract In this paper, a new stochastic model updating method is proposed, which combines the random hybrid perturbation-Galerkin method with the improved cross-model cross-mode technique. This method effectively alleviates the impaction of limited measurement data and uncertain measurement errors on model updating. Considering the uncertainty of the measured modal data, a new stochastic updating equation with update coefficient vector is established based on the improved cross-model cross-mode method. Using the hybrid perturbation-Galerkin method to solve the stochastic updated equation, the update coefficient vector is obtained. The statistical characteristics of the update coefficients can then be determined. The numerical results of the simply supported beam show that the proposed method can effectively deal with the relatively large uncertainty in the actual measurement data, and shows relatively strong stability in the case of different modal combinations, and has a higher computational efficiency than the Monte Carlo method. Considering the rank deficit, the improved cross-model cross-mode method proposed in this paper can get better updating results than the cross-model cross-mode method. The experimental results of the frame show that the new method can simultaneously modify the stiffness and the quality of the structure, and the updated model can be used to obtain modal data consistent with the measured results, thus verifying the effectiveness of the proposed method.

    Keywords stochastic model updating; hybrid perturbation-Galerkin method; improved cross-model cross-mode technique

    猜你喜歡
    交叉修正模態(tài)
    Some new thoughts of definitions of terms of sedimentary facies: Based on Miall's paper(1985)
    修正這一天
    快樂語文(2021年35期)2022-01-18 06:05:30
    合同解釋、合同補充與合同修正
    法律方法(2019年4期)2019-11-16 01:07:28
    “六法”巧解分式方程
    軟件修正
    連一連
    國內多模態(tài)教學研究回顧與展望
    基于Fast-ICA的Wigner-Ville分布交叉項消除方法
    計算機工程(2015年8期)2015-07-03 12:19:54
    基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識別
    由單個模態(tài)構造對稱簡支梁的抗彎剛度
    計算物理(2014年2期)2014-03-11 17:01:39
    欧美精品一区二区免费开放| 十八禁高潮呻吟视频| 热re99久久精品国产66热6| 老司机在亚洲福利影院| 色老头精品视频在线观看| 国产精品1区2区在线观看. | 亚洲专区国产一区二区| 水蜜桃什么品种好| 12—13女人毛片做爰片一| 国产精品国产三级国产专区5o| 亚洲伊人色综图| 亚洲国产日韩一区二区| 桃红色精品国产亚洲av| 国产精品 欧美亚洲| 热re99久久国产66热| 亚洲中文字幕日韩| 又紧又爽又黄一区二区| 欧美精品人与动牲交sv欧美| 日韩视频一区二区在线观看| 老汉色av国产亚洲站长工具| 国产精品九九99| www.av在线官网国产| 国产在线免费精品| 国产av精品麻豆| 免费观看av网站的网址| 亚洲一区二区三区欧美精品| 亚洲天堂av无毛| 日本撒尿小便嘘嘘汇集6| 亚洲欧洲日产国产| 亚洲情色 制服丝袜| 国产成人免费无遮挡视频| 我要看黄色一级片免费的| 久久人人爽人人片av| 日韩制服丝袜自拍偷拍| 久久精品熟女亚洲av麻豆精品| 黄色视频不卡| 又大又爽又粗| 十分钟在线观看高清视频www| 亚洲欧洲日产国产| 久久人人爽人人片av| 他把我摸到了高潮在线观看 | 欧美日韩福利视频一区二区| 中文字幕高清在线视频| 老熟妇仑乱视频hdxx| 91精品三级在线观看| 91国产中文字幕| 国产日韩一区二区三区精品不卡| 美女国产高潮福利片在线看| 亚洲专区字幕在线| 色精品久久人妻99蜜桃| 国产福利在线免费观看视频| 欧美精品一区二区大全| 青草久久国产| 国产麻豆69| 亚洲成国产人片在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 久久九九热精品免费| 亚洲七黄色美女视频| 日韩人妻精品一区2区三区| 欧美日韩福利视频一区二区| 美女扒开内裤让男人捅视频| 看免费av毛片| 亚洲欧美精品综合一区二区三区| 中文字幕制服av| 在线观看免费视频网站a站| 99国产精品免费福利视频| 日本a在线网址| 日本wwww免费看| 日韩制服骚丝袜av| 国产一区二区在线观看av| 亚洲精品乱久久久久久| 午夜免费鲁丝| 又紧又爽又黄一区二区| 久久精品熟女亚洲av麻豆精品| 国产一区二区激情短视频 | 美女福利国产在线| 国产精品99久久99久久久不卡| 热re99久久精品国产66热6| 国产主播在线观看一区二区| 99国产综合亚洲精品| 久久九九热精品免费| 久久久国产一区二区| 丰满少妇做爰视频| av免费在线观看网站| 极品少妇高潮喷水抽搐| 亚洲视频免费观看视频| 成年美女黄网站色视频大全免费| 精品熟女少妇八av免费久了| 久久热在线av| 搡老熟女国产l中国老女人| 日韩中文字幕欧美一区二区| 亚洲国产成人一精品久久久| h视频一区二区三区| 色94色欧美一区二区| 电影成人av| 一本色道久久久久久精品综合| 91麻豆av在线| 看免费av毛片| 欧美日韩av久久| 亚洲三区欧美一区| 亚洲 国产 在线| 日韩精品免费视频一区二区三区| 十八禁网站网址无遮挡| 国产成人精品无人区| 日本av免费视频播放| 老司机福利观看| av国产精品久久久久影院| 手机成人av网站| 男人舔女人的私密视频| 久久狼人影院| av欧美777| 国产免费现黄频在线看| 国产精品久久久久久人妻精品电影 | 50天的宝宝边吃奶边哭怎么回事| 亚洲欧美清纯卡通| 亚洲欧美激情在线| 波多野结衣av一区二区av| 在线精品无人区一区二区三| 男女免费视频国产| 亚洲色图 男人天堂 中文字幕| 亚洲综合色网址| 国产成人免费无遮挡视频| 国产日韩欧美视频二区| 国产亚洲av片在线观看秒播厂| 国产91精品成人一区二区三区 | 国产精品 国内视频| 亚洲五月婷婷丁香| 狂野欧美激情性xxxx| xxxhd国产人妻xxx| 在线十欧美十亚洲十日本专区| 老汉色∧v一级毛片| 久久亚洲精品不卡| 捣出白浆h1v1| 电影成人av| 美女主播在线视频| 夫妻午夜视频| 亚洲色图综合在线观看| 国产精品99久久99久久久不卡| 久久久国产欧美日韩av| 国产精品二区激情视频| 中文字幕人妻丝袜一区二区| 亚洲av日韩精品久久久久久密| 捣出白浆h1v1| 91av网站免费观看| 久久久精品国产亚洲av高清涩受| 国产成人a∨麻豆精品| 12—13女人毛片做爰片一| 中文字幕人妻丝袜一区二区| 日韩人妻精品一区2区三区| 亚洲欧洲日产国产| 水蜜桃什么品种好| 80岁老熟妇乱子伦牲交| 亚洲av欧美aⅴ国产| 日韩一卡2卡3卡4卡2021年| 久久人人爽av亚洲精品天堂| 欧美精品亚洲一区二区| 国产成人系列免费观看| av欧美777| 日韩,欧美,国产一区二区三区| 亚洲精华国产精华精| 一区二区三区精品91| 女人精品久久久久毛片| 精品少妇久久久久久888优播| 777久久人妻少妇嫩草av网站| 精品人妻一区二区三区麻豆| 欧美黑人欧美精品刺激| 欧美日韩成人在线一区二区| 国产一区二区 视频在线| 黄色a级毛片大全视频| 亚洲中文日韩欧美视频| 国产精品久久久久成人av| 天天躁夜夜躁狠狠躁躁| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人系列免费观看| 80岁老熟妇乱子伦牲交| 精品人妻熟女毛片av久久网站| 亚洲国产日韩一区二区| 国产免费福利视频在线观看| tocl精华| 十分钟在线观看高清视频www| 国产精品欧美亚洲77777| 日日摸夜夜添夜夜添小说| 一本大道久久a久久精品| 99热全是精品| 无遮挡黄片免费观看| 午夜免费观看性视频| 大香蕉久久网| 国内毛片毛片毛片毛片毛片| 狠狠狠狠99中文字幕| 大香蕉久久成人网| 中国国产av一级| 午夜福利影视在线免费观看| 免费日韩欧美在线观看| 97精品久久久久久久久久精品| 自线自在国产av| 免费黄频网站在线观看国产| 国产一级毛片在线| 国产在线一区二区三区精| 日韩欧美一区二区三区在线观看 | 日本猛色少妇xxxxx猛交久久| 日本猛色少妇xxxxx猛交久久| 久久热在线av| 欧美人与性动交α欧美软件| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲男人天堂网一区| 少妇裸体淫交视频免费看高清 | 国产日韩欧美视频二区| 亚洲成人国产一区在线观看| 日本a在线网址| 欧美中文综合在线视频| 亚洲欧美激情在线| 国产淫语在线视频| 国产伦人伦偷精品视频| 爱豆传媒免费全集在线观看| 日日夜夜操网爽| 激情视频va一区二区三区| 欧美精品高潮呻吟av久久| 亚洲精品美女久久久久99蜜臀| 又大又爽又粗| 中国国产av一级| 九色亚洲精品在线播放| 日韩精品免费视频一区二区三区| 免费高清在线观看视频在线观看| 中文字幕最新亚洲高清| 天堂俺去俺来也www色官网| 亚洲av美国av| 久久性视频一级片| 日本猛色少妇xxxxx猛交久久| 中文字幕精品免费在线观看视频| 久久久久久久大尺度免费视频| 国产伦理片在线播放av一区| 亚洲av国产av综合av卡| 天堂8中文在线网| 久久久久久人人人人人| 91av网站免费观看| 中国国产av一级| 精品人妻1区二区| 一级,二级,三级黄色视频| 国精品久久久久久国模美| 蜜桃国产av成人99| 国产欧美日韩一区二区三 | 日本一区二区免费在线视频| 久久久久久久久免费视频了| 伊人久久大香线蕉亚洲五| 精品国产一区二区三区久久久樱花| 菩萨蛮人人尽说江南好唐韦庄| 91精品伊人久久大香线蕉| 纯流量卡能插随身wifi吗| www.999成人在线观看| 欧美日韩一级在线毛片| 国产一区二区 视频在线| 在线观看www视频免费| 欧美激情 高清一区二区三区| 18在线观看网站| 天堂8中文在线网| 欧美日韩亚洲高清精品| 亚洲av日韩精品久久久久久密| 国产麻豆69| 亚洲欧美清纯卡通| 搡老岳熟女国产| 男女下面插进去视频免费观看| www.精华液| 天堂俺去俺来也www色官网| 美女高潮喷水抽搐中文字幕| 国产国语露脸激情在线看| 一区二区三区乱码不卡18| 日本五十路高清| 免费女性裸体啪啪无遮挡网站| 亚洲免费av在线视频| 99久久综合免费| 亚洲少妇的诱惑av| 欧美xxⅹ黑人| 午夜激情久久久久久久| 狂野欧美激情性bbbbbb| 日韩欧美一区二区三区在线观看 | 欧美乱码精品一区二区三区| 国产日韩欧美在线精品| 精品福利永久在线观看| 国产日韩一区二区三区精品不卡| 另类精品久久| videosex国产| 国产成人一区二区三区免费视频网站| 久久精品国产综合久久久| 在线观看一区二区三区激情| 操美女的视频在线观看| 91字幕亚洲| 亚洲人成电影观看| 黑人巨大精品欧美一区二区mp4| 免费观看av网站的网址| 99久久综合免费| 十八禁网站网址无遮挡| 青青草视频在线视频观看| 午夜免费观看性视频| 国产精品亚洲av一区麻豆| 欧美 日韩 精品 国产| 亚洲国产欧美一区二区综合| 丝袜脚勾引网站| 久久久国产精品麻豆| av一本久久久久| 午夜影院在线不卡| 亚洲精品乱久久久久久| 日本五十路高清| 最黄视频免费看| 亚洲欧美日韩另类电影网站| 亚洲激情五月婷婷啪啪| 水蜜桃什么品种好| 伊人亚洲综合成人网| 亚洲自偷自拍图片 自拍| 人人妻人人澡人人爽人人夜夜| 精品少妇久久久久久888优播| 国产精品久久久人人做人人爽| av又黄又爽大尺度在线免费看| 久久精品国产a三级三级三级| 婷婷色av中文字幕| 啪啪无遮挡十八禁网站| 成年人午夜在线观看视频| 香蕉丝袜av| 久久国产亚洲av麻豆专区| 欧美激情 高清一区二区三区| 亚洲国产中文字幕在线视频| 成年动漫av网址| 日韩大片免费观看网站| 日本91视频免费播放| 国产精品秋霞免费鲁丝片| 婷婷色av中文字幕| 亚洲天堂av无毛| 韩国高清视频一区二区三区| 免费人妻精品一区二区三区视频| 亚洲熟女毛片儿| 一级,二级,三级黄色视频| 岛国在线观看网站| cao死你这个sao货| 在线 av 中文字幕| 交换朋友夫妻互换小说| a 毛片基地| 久久99一区二区三区| 一本大道久久a久久精品| 久久国产精品大桥未久av| 久久国产亚洲av麻豆专区| 日本黄色日本黄色录像| 看免费av毛片| 大香蕉久久网| 丝袜美足系列| 黄色视频,在线免费观看| 黄色怎么调成土黄色| 人人妻,人人澡人人爽秒播| 午夜激情av网站| 精品人妻在线不人妻| 久久九九热精品免费| 中文字幕高清在线视频| 国产精品香港三级国产av潘金莲| 青草久久国产| 丝袜美足系列| 成年动漫av网址| 午夜精品久久久久久毛片777| 亚洲va日本ⅴa欧美va伊人久久 | 精品国产一区二区三区久久久樱花| 欧美另类一区| 成人国产一区最新在线观看| 啪啪无遮挡十八禁网站| 亚洲精品av麻豆狂野| 日本欧美视频一区| 亚洲色图 男人天堂 中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 久久午夜综合久久蜜桃| 国产av一区二区精品久久| 欧美国产精品一级二级三级| 91麻豆精品激情在线观看国产 | 国产欧美日韩精品亚洲av| 国产精品亚洲av一区麻豆| 欧美久久黑人一区二区| 黄色毛片三级朝国网站| 极品少妇高潮喷水抽搐| 18在线观看网站| 精品久久久精品久久久| 国产伦理片在线播放av一区| 老司机午夜十八禁免费视频| 国产精品免费视频内射| 久久久国产精品麻豆| 国产高清videossex| 日本五十路高清| 亚洲人成电影观看| 夜夜骑夜夜射夜夜干| 国产一区二区三区综合在线观看| 99久久精品国产亚洲精品| 亚洲五月色婷婷综合| 黑丝袜美女国产一区| 久久热在线av| 国产真人三级小视频在线观看| 在线观看免费午夜福利视频| 欧美日韩av久久| 黑丝袜美女国产一区| 精品国产国语对白av| 国产精品自产拍在线观看55亚洲 | 大片电影免费在线观看免费| 另类亚洲欧美激情| 9色porny在线观看| 99国产精品免费福利视频| av在线播放精品| 久久精品人人爽人人爽视色| 一区二区三区乱码不卡18| 免费观看av网站的网址| 日韩 亚洲 欧美在线| a在线观看视频网站| 另类精品久久| 日本wwww免费看| 91成人精品电影| 自拍欧美九色日韩亚洲蝌蚪91| 午夜福利乱码中文字幕| 久久午夜综合久久蜜桃| 成年人免费黄色播放视频| 国产人伦9x9x在线观看| 在线观看人妻少妇| 免费观看av网站的网址| 亚洲黑人精品在线| av在线播放精品| 精品视频人人做人人爽| 亚洲五月婷婷丁香| 九色亚洲精品在线播放| 亚洲av日韩在线播放| 美女扒开内裤让男人捅视频| 老熟女久久久| 日韩大码丰满熟妇| 不卡av一区二区三区| 亚洲人成77777在线视频| 精品欧美一区二区三区在线| 久久久久久人人人人人| 国产黄频视频在线观看| 国产精品99久久99久久久不卡| 精品国产乱码久久久久久男人| 脱女人内裤的视频| 国产精品国产av在线观看| 人人妻人人澡人人看| 亚洲精品第二区| 精品福利永久在线观看| 我要看黄色一级片免费的| 成人亚洲精品一区在线观看| 久久99热这里只频精品6学生| 久久精品aⅴ一区二区三区四区| 亚洲第一欧美日韩一区二区三区 | 中国国产av一级| 日韩欧美一区视频在线观看| 成人国语在线视频| 成年av动漫网址| 亚洲一卡2卡3卡4卡5卡精品中文| 国产亚洲av高清不卡| 久久久久久久久久久久大奶| 一本大道久久a久久精品| 亚洲av电影在线进入| 中文字幕av电影在线播放| 高清视频免费观看一区二区| 午夜久久久在线观看| 国产成人系列免费观看| 在线av久久热| 一区福利在线观看| 亚洲情色 制服丝袜| 女性生殖器流出的白浆| 不卡av一区二区三区| 18在线观看网站| 丰满饥渴人妻一区二区三| 三级毛片av免费| 动漫黄色视频在线观看| 三级毛片av免费| 日本欧美视频一区| 亚洲情色 制服丝袜| 国产免费现黄频在线看| 亚洲午夜精品一区,二区,三区| 十八禁人妻一区二区| 午夜精品国产一区二区电影| 一本—道久久a久久精品蜜桃钙片| 一区二区日韩欧美中文字幕| 黄色a级毛片大全视频| 90打野战视频偷拍视频| 9热在线视频观看99| 成人手机av| 自线自在国产av| 中文精品一卡2卡3卡4更新| 他把我摸到了高潮在线观看 | 欧美 日韩 精品 国产| 最新在线观看一区二区三区| 久久毛片免费看一区二区三区| 色视频在线一区二区三区| 日韩熟女老妇一区二区性免费视频| 99国产精品一区二区三区| 日韩大码丰满熟妇| 99久久人妻综合| av网站在线播放免费| 人妻一区二区av| 国产免费一区二区三区四区乱码| 电影成人av| 欧美日韩黄片免| 国产成人免费观看mmmm| 男人添女人高潮全过程视频| 美女扒开内裤让男人捅视频| 欧美黑人精品巨大| 国产高清国产精品国产三级| 蜜桃在线观看..| 十八禁人妻一区二区| e午夜精品久久久久久久| 狠狠精品人妻久久久久久综合| 在线观看免费日韩欧美大片| 亚洲va日本ⅴa欧美va伊人久久 | 丝袜人妻中文字幕| 精品少妇一区二区三区视频日本电影| a级毛片在线看网站| 国产精品99久久99久久久不卡| 精品福利永久在线观看| 老熟妇乱子伦视频在线观看 | 午夜成年电影在线免费观看| 日韩视频在线欧美| 成人三级做爰电影| 国产一卡二卡三卡精品| 色综合欧美亚洲国产小说| 免费少妇av软件| 91麻豆av在线| 男女无遮挡免费网站观看| 中文精品一卡2卡3卡4更新| 久久国产精品男人的天堂亚洲| 午夜久久久在线观看| 日韩三级视频一区二区三区| 色老头精品视频在线观看| 制服诱惑二区| 自线自在国产av| 国产激情久久老熟女| 黄色视频不卡| 欧美日韩av久久| 国产欧美日韩综合在线一区二区| 在线观看免费视频网站a站| 亚洲精品av麻豆狂野| 69av精品久久久久久 | 亚洲精品国产精品久久久不卡| 老司机影院成人| 69精品国产乱码久久久| 国产精品久久久av美女十八| 日本av手机在线免费观看| 亚洲精品中文字幕一二三四区 | 97精品久久久久久久久久精品| 国产欧美日韩综合在线一区二区| 一区二区日韩欧美中文字幕| 大片免费播放器 马上看| 一级,二级,三级黄色视频| 丝瓜视频免费看黄片| www.熟女人妻精品国产| 黑丝袜美女国产一区| 国产精品国产三级国产专区5o| 桃花免费在线播放| 午夜福利,免费看| 免费在线观看视频国产中文字幕亚洲 | 久久精品国产综合久久久| 国产黄色免费在线视频| 女人高潮潮喷娇喘18禁视频| 中文字幕最新亚洲高清| 91av网站免费观看| 美女视频免费永久观看网站| 久久久久国产一级毛片高清牌| 欧美黑人欧美精品刺激| 欧美精品啪啪一区二区三区 | 少妇 在线观看| 日本撒尿小便嘘嘘汇集6| 久久亚洲国产成人精品v| 又黄又粗又硬又大视频| 国产熟女午夜一区二区三区| 精品视频人人做人人爽| 亚洲中文日韩欧美视频| 午夜福利在线观看吧| 国产在线一区二区三区精| 久久午夜综合久久蜜桃| 性色av乱码一区二区三区2| 激情视频va一区二区三区| 精品国产乱码久久久久久小说| 一个人免费看片子| 法律面前人人平等表现在哪些方面 | 精品第一国产精品| 国产不卡av网站在线观看| 色视频在线一区二区三区| 精品一品国产午夜福利视频| 国产欧美日韩综合在线一区二区| 国产亚洲一区二区精品| 亚洲天堂av无毛| 精品亚洲乱码少妇综合久久| 国产精品亚洲av一区麻豆| 女人高潮潮喷娇喘18禁视频| 欧美日韩亚洲国产一区二区在线观看 | 国产伦理片在线播放av一区| 国精品久久久久久国模美| 俄罗斯特黄特色一大片| 久久久久网色| 国产精品偷伦视频观看了| 少妇粗大呻吟视频| 国产区一区二久久| 在线观看免费视频网站a站| 美女扒开内裤让男人捅视频| 免费在线观看完整版高清| 日日摸夜夜添夜夜添小说| 91av网站免费观看| 狂野欧美激情性bbbbbb| 后天国语完整版免费观看| 自线自在国产av| 久久久久视频综合| 欧美黑人欧美精品刺激| 9热在线视频观看99| 亚洲av男天堂| 在线天堂中文资源库| 欧美精品亚洲一区二区| 十分钟在线观看高清视频www| 在线天堂中文资源库| 亚洲五月色婷婷综合| av又黄又爽大尺度在线免费看| 女性生殖器流出的白浆| 亚洲中文字幕日韩| 2018国产大陆天天弄谢| 人人妻人人澡人人爽人人夜夜| 啦啦啦啦在线视频资源| 亚洲avbb在线观看| 亚洲九九香蕉|