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

    基于Tikhonov正則化迭代求解的結(jié)構(gòu)損傷識別方法

    2019-09-17 06:46:56夏志鵬王樹青徐明強王皓宇
    振動與沖擊 2019年17期
    關(guān)鍵詞:正則模態(tài)程度

    夏志鵬, 王樹青, 徐明強, 王皓宇

    (中國海洋大學 海洋工程系, 山東 青島 266100)

    海洋平臺結(jié)構(gòu)長期服役在惡劣的海洋環(huán)境中,容易產(chǎn)生各種形式的損傷,使結(jié)構(gòu)的承載能力下降,甚至導致平臺失效,造成巨大經(jīng)濟損失及人員傷亡[1]。因此,針對海洋平臺結(jié)構(gòu)的健康監(jiān)測與損傷識別非常重要。目前結(jié)構(gòu)損傷檢測的方法眾多[2-3],基于振動測試的結(jié)構(gòu)健康監(jiān)測技術(shù)相對較為簡單且成本較低,是非常具有發(fā)展前景的損傷識別技術(shù)[4]。其中,基于模態(tài)參數(shù)的損傷檢測是近年來新興且有效的檢測手段。

    在某些情況下,基于模態(tài)參數(shù)的結(jié)構(gòu)損傷識別過程可以簡化為線性方程組Cα=b的求解問題。當不考慮測量噪聲或噪聲干擾較小時,該系統(tǒng)的求解往往可以得到滿意的結(jié)果;而當結(jié)構(gòu)測量模態(tài)信息受噪聲影響較為嚴重時,系統(tǒng)的求解結(jié)果往往會振蕩發(fā)散,導致檢測方法失效。因此,噪聲魯棒性是此類損傷識別技術(shù)在發(fā)展過程中必須考慮的問題。

    從數(shù)學的角度看,利用結(jié)構(gòu)的振動測試數(shù)據(jù)識別其損傷是求解反問題的過程,其不適定性體現(xiàn)在所構(gòu)建系統(tǒng)的病態(tài)上,即微小的測量誤差都可能導致解的振蕩發(fā)散。為解決這一問題,學者們做了大量研究[5-11]。其中,基于Tikhonov正則化[12]的方法應(yīng)用較為廣泛,其基本思想是:用一族與原問題相“鄰近”的適定問題的解去逼近原問題的真實解。王藝霖等將Tikhonov 正則化方法用于結(jié)構(gòu)的損傷識別,一定程度上改善了系統(tǒng)的不適定性,提高了損傷識別精度,并指出正則化方法的作用可通過剛度矩陣條件數(shù)的減小來明確衡量。Hua等將Tikhonov正則化與基于靈敏度分析的有限元模型修正方法相結(jié)合,進行了簡單框架結(jié)構(gòu)的損傷識別,該方法體現(xiàn)出較好的抗噪性。應(yīng)用Tikhonov正則化雖然能在一定程度上抑制噪聲,改善損傷識別結(jié)果的穩(wěn)定性,但它的解是過度光滑的,不具有稀疏性[13]。換言之,應(yīng)用該方法雖然能篩選出真實的損傷信息,但也引入了過多的虛假損傷信息,給損傷識別帶來困難。為改善這一問題,張純等[14]在Tikhonov罰函數(shù)項中引入光滑函數(shù),并結(jié)合基于靈敏度分析的模型修正方法進行損傷識別研究,明顯改善了損傷識別效果;張純等[15-16]先后將L1和L1/2范數(shù)正則化模型修正方法應(yīng)用于結(jié)構(gòu)損傷識別,有效地改善了基于Tikhonov正則化損傷識別結(jié)果過度光滑的缺陷。

    利用正則化方法求解線性不適定系統(tǒng),關(guān)鍵在于正則化參數(shù)的選擇。上述正則化及其改進方法大都采用L曲線法選取正則化參數(shù),而L曲線的繪制過程往往需要進行大量的試算,對于大型的結(jié)構(gòu)系統(tǒng),試算過程更加復(fù)雜。對于Tikhonov正則化,其參數(shù)的合理取值范圍一般較小,系統(tǒng)求解效果對正則化參數(shù)的依賴性很高,并且當噪聲水平較高時,L曲線趨于平直,很難通過該方法選取合適的正則化參數(shù)[17]。

    解的過度光滑及正則化參數(shù)選擇時試算量過大是制約Tikhonov正則化應(yīng)用于損傷檢測的兩大難題。對此,本文綜合上述正則化方法的優(yōu)勢及缺陷,提出了一種基于Tikhonov正則化的迭代求解方法,用于測量噪聲影響下?lián)p傷識別線性系統(tǒng)的求解。該方法基于Tikhonov正則化理論,能夠在迭代過程中利用二分法自適應(yīng)調(diào)整正則化參數(shù),并通過迭代的方式重構(gòu)正則化權(quán)重矩陣,能夠充分抑制噪聲,保留實際損傷信息。通過一個數(shù)值算例,驗證了本文方法的有效性。

    1 基本理論

    1.1 交叉模態(tài)應(yīng)變能方法

    交叉模態(tài)應(yīng)變能(CMSE)方法[18-21]是一種典型的以線性系統(tǒng)求解結(jié)構(gòu)損傷信息的識別方法。該方法計算簡單,并且可以同時識別結(jié)構(gòu)的損傷位置和損傷程度。與其他的檢測方法相比,該方法有如下優(yōu)點:不需要質(zhì)量歸一化的振型;僅利用結(jié)構(gòu)損傷前后的較少幾階模態(tài)即可判斷結(jié)構(gòu)的健康狀況。因而,應(yīng)用CMSE方法進行損傷檢測具有很好的發(fā)展前景。

    設(shè)K和M分別為結(jié)構(gòu)剛度和質(zhì)量矩陣,λi和Φi分別表示第i階特征值和特征向量。結(jié)構(gòu)損傷前后的特征方程可表示為

    KΦi=λiMΦi

    (1)

    (2)

    本文中,上標“*”用來表示損傷后的情況。

    (3)

    (4)

    對式(3)兩邊取轉(zhuǎn)置,并考慮矩陣K和M的對稱性得

    (5)

    聯(lián)立式(4)、(5),并考慮到M*=M,可以得到:

    (6)

    本文中,單元損傷通過桿件模量的等效折減進行模擬,即對第n個結(jié)構(gòu)單元,其損傷后的模量值

    (7)

    式中,En為損傷前第n個單元的模量,αn為相應(yīng)的損傷程度。則損傷后結(jié)構(gòu)的剛度矩陣可寫為

    (8)

    式中:Ne為結(jié)構(gòu)單元的總數(shù);Kn為第n個單元損傷前的單元剛度矩陣。

    將式(8)代入式(6)并整理后得

    (9)

    用Ni和Nj分別表示結(jié)構(gòu)損傷前后的模態(tài)階數(shù),可以構(gòu)建Nq=Ni×Nj個線性方程,將其寫成矩陣形式

    Cα=b

    (10)

    (11)

    可采用奇異值分解法(SVD)尋求最小二乘解。對系數(shù)矩陣C作奇異值分解

    (12)

    (13)

    當測量模態(tài)空間不完備時,采用Guyan擴階法[22]對不完備振型進行如下擴階處理

    (14)

    1.2 Tikhonov正則化方法

    Tikhonov方法的思想是將最小二乘解和最小范數(shù)解綜合考慮。對于式(9)所示的線性系統(tǒng),可建立目標函數(shù)如下

    (15)

    其中L為Ne階正則化矩陣,一般作變換將L轉(zhuǎn)化為單位矩陣I[23]。ξ為正則化參數(shù),用于控制最小二乘解和最小范數(shù)解之間的平衡。

    當L為單位矩陣時,式(14)等價于求解

    (CTC+ξ2I)α=CTb

    (16)

    該方程的解為

    (17)

    對系數(shù)矩陣C作式(11)所示的奇異值分解,則式(17)可表示為

    (18)

    其中

    (19)

    稱為Tikhonov過濾因子,顯然,f(σi)的大小依賴于奇異值σi和正則化參數(shù)ξ。而當ξ?σi時,式(18)轉(zhuǎn)化式(13)所示的最小二乘解。

    2 基于Tikhonov正則化的迭代求解方法

    Tikhonov方法的正則化矩陣一般為單位矩陣,這表明Tikhonov正則化方法對各單元損傷參數(shù)進行了相同的加權(quán)平滑處理,其解不具有稀疏性。而對于損傷結(jié)構(gòu)而言,其損傷參數(shù)是稀疏分布的,即除個別損傷單元外,解向量的大部分元素均為零。因此,傳統(tǒng)的Tikhonov正則化應(yīng)用于損傷檢測時存在一定的缺陷。適合損傷識別的理想正則化方法應(yīng)具有如下特性:對于真實的損傷信息,正則化的平滑作用較小以便保留結(jié)構(gòu)的實際損傷信息,即正則化矩陣中對應(yīng)的權(quán)重較小;當出現(xiàn)虛假損傷時,增強平滑作用以消除噪聲導致的虛假損傷信息,即正則化矩陣中對應(yīng)的權(quán)重較大。為了達到這一效果,本文考慮通過具有稀疏性的損傷參數(shù)來重新構(gòu)造正則化矩陣,構(gòu)造形式如下

    L(α)=diag{(|α1|+ε)-1,…,(|αn|+ε)-1}

    (20)

    其中αi為第i個單元對應(yīng)的損傷參數(shù),ε為一正常數(shù),避免因αi→0而導致矩陣對角元素趨于∞。

    當L為一般的對角矩陣時,式(15)等價于求解

    (CTC+ξ2LTL)α=CTb

    (21)

    該方程的解為

    (22)

    由于正則化矩陣L(α)需用損傷參數(shù)α來構(gòu)造,而α為待求的未知量,因此,需要采用迭代的方式進行求解。

    需要注意的是,正則化參數(shù)ξ的選取是利用正則化方法求解線性不適定問題的關(guān)鍵。前述研究中,學者們一般采用L曲線法選取正則化參數(shù),而L曲線的繪制通常需要進行大量的試算,尤其是對于復(fù)雜的大型結(jié)構(gòu)系統(tǒng)。此外,當噪聲水平較高時,L曲線趨于平直,很難通過該方法選取合適的正則化參數(shù)。為解決正則化參數(shù)選取時需要大量試算的問題,提出了一種奇異值二分法用于本文新方法中正則化參數(shù)的選取。

    基于Tikhonov正則化迭代求解(Tikhonov Regularization Iterative Method,TRIM)的損傷識別過程如下:

    步驟1構(gòu)建損傷檢測線性系統(tǒng),求得系統(tǒng)的初始解

    根據(jù)結(jié)構(gòu)損傷前后的模態(tài)信息利用損傷檢測方法(如CMSE)構(gòu)建方程組Cα=b。對系數(shù)矩陣C進行奇異值分解,獲得系統(tǒng)的r個奇異值σ1≥σ2≥…≥σr≥0。

    (23)

    步驟2調(diào)整正則化參數(shù),重構(gòu)正則化矩陣,進行第k次迭代求解

    (1) 基于奇異值二分法調(diào)整正則化參數(shù)ξ

    (24)

    相應(yīng)的正則化參數(shù)調(diào)整為

    (25)

    上述選擇和調(diào)整正則化參數(shù)的方法,我們稱之為奇異值二分法。

    (26)

    (27)

    此處ε統(tǒng)一取為0.01%,以減小對損傷檢測結(jié)果的影響。

    (3) 進行第k次求解

    應(yīng)用新的迭代參數(shù)ξk和矩陣Lk,計算新的解向量

    (28)

    步驟3重復(fù)步驟2所述過程,直到前后兩次迭代求解的向量差值滿足如下條件

    (29)

    其中δ為收斂閾值,本文建議δ=0.1%,在保證求解精度的同時,控制所需要的迭代步數(shù)?;谛路椒ǖ淖罱K損傷識別結(jié)果為

    (30)

    基于Tikhonov正則化迭代求解的CMSE損傷檢測(簡稱CMSE-TRIM)流程如圖1所示。

    圖1 CMSE-TRIM損傷檢測流程圖

    3 數(shù)值算例

    3.1 海洋平臺結(jié)構(gòu)

    本算例的研究對象為一個近海導管架平臺結(jié)構(gòu)[24],如圖2所示。該結(jié)構(gòu)由36個外徑為17.8 cm,壁厚為0.89 cm的均勻鋼管構(gòu)件組成,每根鋼管劃分為一個單元。結(jié)構(gòu)三層標高分別為9.14 m,18.29 m和27.43 m。底部和頂部的邊長分別為10.97 m×10.97 m和3.66 m×3.66 m。平臺的四條腿固定于地面。選用材料為鋼材,彈性模量2.1×1011Pa,密度7 850 kg/m3,截面面積2.825×10-3m2,截面慣性矩2.89×10-6m4。按照有限元方法組建整體質(zhì)量矩陣(采用集中質(zhì)量法)和整體剛度矩陣,進行特征值分析,得到有限元模型的前三階頻率分別為7.51 Hz、8.80 Hz和8.88 Hz。平臺前三階振型如圖3所示。

    圖2 海洋平臺結(jié)構(gòu)示意圖

    圖3 有限元模型前三階振型

    3.2 損傷工況設(shè)置

    由于實際測試中結(jié)構(gòu)的高階模態(tài)一般難以激勵,轉(zhuǎn)動自由度信息難以測量,因此本文算例中假設(shè)僅能測得損傷結(jié)構(gòu)1~3階頻率以及振型中的平動自由度信息,并且所測振型信息含有一定水平的噪聲。

    為了研究本文提出的方法在不同因素影響下的損傷識別效果,算例中設(shè)置了四類影響因素:① 損傷位置敏感性;② 損傷程度敏感性;③ 噪聲水平敏感性;④ 模態(tài)階數(shù)敏感性。分別探討不同影響因素下的損傷識別效果。結(jié)構(gòu)損傷后的振型信息添加一定水平噪聲,噪聲添加方式如下

    φi,j=φi,j(1+μGi,j)

    (31)

    式中,μ為噪聲水平;Gi,j為均值為0、方差為1的高斯隨機數(shù);φi,j為對應(yīng)于第i個自由度的第j階振型值。實測頻率信息一般較為準確,因此算例中不考慮噪聲對頻率的影響。

    4 損傷識別與結(jié)果分析

    本文采用Guyan擴階方法對仿真模擬得到的損傷后模態(tài)振型進行擴階處理,利用損傷前后的模態(tài)信息構(gòu)建CMSE方程,求解36個單元(n=36)的損傷參數(shù)。并對比CMSE-Tikhonov(基于Tikhonov正則化求解的CMSE損傷檢測)和CMSE-TRIM的損傷識別效果。

    4.1 損傷位置敏感性

    設(shè)置噪聲水平1%,損傷程度25%,選用結(jié)構(gòu)損傷后1~3階模態(tài),研究損傷方法對不同類型損傷構(gòu)件(工況S1:13,工況S2:20,工況S3:23)的損傷識別效果。

    工況S1設(shè)置13號單元發(fā)生25%損傷,在測量模態(tài)空間不完備及1%噪聲影響下利用CMSE-Tikhonov方法進行損傷識別,具體求解過程如下:首先根據(jù)1.1節(jié)內(nèi)容構(gòu)建CMSE方程組Cα=b,對系數(shù)矩陣C進行奇異值分解,并將得到的所有奇異值分別假設(shè)為正則化參數(shù),根據(jù)式(16)進行試算并繪制L曲線圖(如圖4);找出L曲線最大曲率點所對應(yīng)的奇異值,作為最終選用的正則化參數(shù),進而求解損傷向量。由L曲線法選取的正則化參數(shù)為ξ=7.20×107,根據(jù)式(16)進行求解,最終識別結(jié)果如圖5(a)所示,其中,深灰色表示干擾項,淺灰色表示損傷單元的損傷指標。此外,圖中DS(damage severity),表示實際損傷單元的損傷程度;RE(relative error),表示評估的損傷程度與真實值的相對誤差;c表示迭代的總次數(shù)。

    應(yīng)用CMSE-TRIM對工況A進行損傷識別的過程如下:構(gòu)建CMSE方程組Cα=b,對矩陣C進行奇異值分解,并直接選取矩陣C的中間奇異值σ18=1.06×108作為初始的正則化參數(shù);按照圖1所示流程進行迭代求解,迭代9次后計算結(jié)果達到收斂閾值,最終求解結(jié)果如圖5(b)所示(圖中c表示迭代次數(shù))。工況B和工況C損傷識別過程與工況A類似,具體流程不再贅述,最終求解結(jié)果分別見圖6和圖7。應(yīng)用TRIM方法所得損傷單元各迭代步下的計算結(jié)果如圖8所示。

    圖4 工況S1 “L曲線”圖

    對比三種不同損傷位置下CMSE-Tikhonov求解結(jié)果,可以看出,該方法對立柱單元損傷檢測效果最好,斜撐次之,橫撐檢測效果最差。相比之下,CMSE-TRIM方法對上述三工況的損傷識別結(jié)果都非常好,無論是在損傷定位還是損傷程度識別方面都明顯優(yōu)于CMSE-Tikhonov方法。從圖8中損傷單元各迭代步計算結(jié)果可以看出,13號橫撐損傷時迭代次數(shù)最多(9次),20號斜撐次之(5次),識別23號立柱損傷所需迭代次數(shù)最少(4次)。迭代次數(shù)的多少側(cè)面反映出不同類型損傷單元的損傷識別難易程度,即相同條件下,立柱損傷最容易識別,斜撐次之,橫撐損傷最不易識別。此外,從圖8中13號單元迭代結(jié)果變化可以看出,該計算結(jié)果經(jīng)歷了由大變小,由小變大,最后逐漸趨于穩(wěn)定這幾個階段,體現(xiàn)出TRIM方法能夠自適應(yīng)調(diào)節(jié)計算結(jié)果的特性。

    (a) CMSE-Tikhonov

    (b) CMSE-TRIM

    (a) CMSE-Tikhonov

    (b) CMSE-TRIM

    (a) CMSE-Tikhonov

    (b) CMSE-TRIM

    圖8 各損傷單元不同迭代步的損傷程度計算結(jié)果

    4.2 損傷程度敏感性

    設(shè)置噪聲水平1%,損傷單元23,選用結(jié)構(gòu)損傷后1~3階模態(tài),研究損傷方法對不同損傷程度(工況D1:5%,工況D2:15%,工況D3:25%)的損傷識別效果。兩種方法識別結(jié)果如圖9~圖11所示。

    由圖9~圖11可以看出,當單元的損傷程度較小時,使用CMSE-Tikhonov方法會出現(xiàn)大量的虛假損傷掩蓋真實的損傷信息,從而無法完全損傷檢測;隨著損傷程度的增加,求解結(jié)果中的干擾項逐漸減少,損傷程度識別誤差也逐漸減小。相比之下,CMSE-TRIM對三種不同損傷程度的工況均能正確識別出損傷單元的位置,而且沒有出現(xiàn)虛假損傷。在損傷程度識別方面,CMSE-TRIM的識別誤差明顯小于CMSE-Tikhonov方法。此外,從計算所需的迭代次數(shù)上可以看出,在其他影響因素相同的條件下,TRIM方法需要較多的迭代次數(shù)來計算小損傷工況。

    4.3 噪聲水平敏感性

    設(shè)置損傷單元23,損傷程度25%,選用結(jié)構(gòu)損傷后1~3階模態(tài),研究損傷方法在不同噪聲水平(工況N1:1%,工況N2:3%,工況N3:5%)下的損傷識別效果。識別結(jié)果如圖12~圖14所示。

    分析圖12~圖14可以看出,當噪聲水平相對較低時,使用CMSE-Tikhonov方法能夠較好地識別出損傷位置及損傷程度,干擾項較少。隨著噪聲水平的提高,該方法損傷識別效果逐漸變差:在損傷定位方面,虛假損傷干擾項逐漸增多,掩蓋了真實的損傷信息;損傷程度識別誤差也逐漸增大。與之相比,CMSE-TRIM在不同噪聲水平影響下均能正確判定出損傷位置,不摻雜任何虛假損傷,而且損傷程度識別精度明顯高于CMSE-Tikhonov。此外,從計算所需的迭代次數(shù)上可以看出,在其他影響因素相同的條件下,當影響實測模態(tài)的噪聲水平提高時,TRIM方法計算所需的迭代次數(shù)相應(yīng)增加。

    (a) CMSE-Tikhonov

    (b) CMSE-TRIM

    (a) CMSE-Tikhonov

    (b) CMSE-TRIM

    (a) CMSE-Tikhonov

    (b) CMSE-TRIM

    (a) CMSE-Tikhonov

    (b) CMSE-TRIM

    (b) CMSE-TRIM

    (a) CMSE-Tikhonov

    (b) CMSE-TRIM

    4.4 模態(tài)階數(shù)敏感性

    設(shè)置噪聲水平1%,損傷單元23,損傷程度25%,研究使用結(jié)構(gòu)損傷后不同階數(shù)模態(tài)(工況M1:1,工況M2:1~2,工況M3:1~3)的損傷識別效果。兩種方法識別結(jié)果如圖15~圖17所示。

    從圖15~圖17可以看出,當取用的損傷后模態(tài)只有第1階時,使用CMSE-Tikhonov和CMSE-TRIM都無法完成損傷識別;當取用的損傷后模態(tài)為前2階時,CMSE-Tikhonov方法雖然能夠識別出23號單元發(fā)生損傷,但是求解結(jié)果中摻雜多處虛假損傷干擾項,損傷程度求解誤差為14.20%;當取用的損傷后模態(tài)為前3階時,CMSE-Tikhonov求解結(jié)果中干擾項明顯減少,損傷程度識別精度有所提高。相比之下,當取用的損傷后模態(tài)為前2階或前3階時,CMSE-TRIM方法均能正確判定出損傷位置,不摻雜任何虛假損傷,而且損傷程度識別精度明顯高于CMSE-Tikhonov。此外,從計算所需的迭代次數(shù)上可以看出,在其他影響因素相同的條件下,取用的損傷后模態(tài)階數(shù)越多,TRIM方法計算所需的迭代次數(shù)越少。

    (a) CMSE-Tikhonov

    (b) CMSE-TRIM

    (a) CMSE-Tikhonov

    (b) CMSE-TRIM

    (a) CMSE-Tikhonov

    (b) CMSE-TRIM

    5 結(jié) 論

    針對測量噪聲影響下?lián)p傷檢測“病態(tài)”方程組的求解問題進行了研究,提出了一種基于Tikhonov正則化的迭代求解方法來抑制測量噪聲對損傷檢測的影響。選取一個海洋平臺結(jié)構(gòu)開展數(shù)值模擬,探討了新方法在不同因素(損傷位置,損傷程度,噪聲水平,取用模態(tài)階數(shù))影響下的損傷識別效果,并同傳統(tǒng)Tikhonov正則化求解下的損傷檢測進行了對比,結(jié)論如下:

    傳統(tǒng)的Tikhonov正則化方法和本文提出的基于Tikhonov正則化的迭代求解方法均能在一定程度上解決噪聲影響下“病態(tài)”方程組的求解問題。相對于傳統(tǒng)的Tikhonov正則化求解方法,本文方法能夠有效地抑制虛假損傷、不遺漏真實損傷,具有更好的噪聲魯棒性。此外,本文方法直接選取方程組系數(shù)矩陣的中間奇異值作為初始正則化參數(shù),通過迭代來自適應(yīng)調(diào)整正則化參數(shù)及矩陣,無需通過大量試算來繪制L曲線圖,減少了計算工作量。數(shù)值分析表明,CMSE-TRIM方法對海洋平臺結(jié)構(gòu)不同類型及不同損程度的損傷工況均能正確地識別,具有較好的噪聲魯棒性。損傷位置不敏感、損傷程度較小、噪聲水平較高以及取用的模態(tài)階數(shù)較少都會增加損傷檢測的難度,相應(yīng)地會增加TRIM方法的迭代計算次數(shù)。

    猜你喜歡
    正則模態(tài)程度
    男女身高受歡迎程度表
    意林(2021年2期)2021-02-08 08:32:47
    剩余有限Minimax可解群的4階正則自同構(gòu)
    類似于VNL環(huán)的環(huán)
    國內(nèi)多模態(tài)教學研究回顧與展望
    基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識別
    有限秩的可解群的正則自同構(gòu)
    由單個模態(tài)構(gòu)造對稱簡支梁的抗彎剛度
    計算物理(2014年2期)2014-03-11 17:01:39
    斷裂對油氣富集程度的控制作用
    斷塊油氣田(2014年6期)2014-03-11 15:33:53
    幸福的程度
    多模態(tài)話語模態(tài)的協(xié)同及在外語教學中的體現(xiàn)
    外語學刊(2010年2期)2010-01-22 03:31:03
    在线观看av片永久免费下载| 中文在线观看免费www的网站| 超碰97精品在线观看| 一级av片app| 成人亚洲精品av一区二区| 久久精品国产鲁丝片午夜精品| 欧美 日韩 精品 国产| 激情五月婷婷亚洲| 亚洲真实伦在线观看| 波野结衣二区三区在线| 久久久久久久久久人人人人人人| 极品少妇高潮喷水抽搐| 亚洲av成人精品一二三区| 婷婷色麻豆天堂久久| 深夜a级毛片| 汤姆久久久久久久影院中文字幕| 中文天堂在线官网| 欧美亚洲 丝袜 人妻 在线| 深夜a级毛片| 99久久中文字幕三级久久日本| 亚洲av成人精品一二三区| 国产成年人精品一区二区| 女人被狂操c到高潮| 亚洲久久久久久中文字幕| 亚洲国产色片| av女优亚洲男人天堂| 国产大屁股一区二区在线视频| 午夜福利视频精品| 丝袜脚勾引网站| 嫩草影院入口| 男人和女人高潮做爰伦理| 免费在线观看成人毛片| 丝瓜视频免费看黄片| tube8黄色片| 国产在线男女| 成人高潮视频无遮挡免费网站| 亚洲精品乱码久久久久久按摩| 爱豆传媒免费全集在线观看| 精品一区在线观看国产| 国产成人91sexporn| 国产成人91sexporn| 日本猛色少妇xxxxx猛交久久| 一级二级三级毛片免费看| 丰满少妇做爰视频| 国产乱来视频区| 干丝袜人妻中文字幕| 蜜桃亚洲精品一区二区三区| 新久久久久国产一级毛片| 涩涩av久久男人的天堂| 国产在线男女| 黄色配什么色好看| 身体一侧抽搐| 精品熟女少妇av免费看| 成人毛片60女人毛片免费| 交换朋友夫妻互换小说| 欧美日韩视频精品一区| 丰满人妻一区二区三区视频av| 久久久精品欧美日韩精品| 大片免费播放器 马上看| 久久久久久久精品精品| 永久免费av网站大全| 久久99热这里只有精品18| 亚洲av免费在线观看| 啦啦啦中文免费视频观看日本| 好男人在线观看高清免费视频| 麻豆成人午夜福利视频| 国产一级毛片在线| 三级国产精品片| 极品少妇高潮喷水抽搐| 国产成人一区二区在线| 国产有黄有色有爽视频| 97人妻精品一区二区三区麻豆| 各种免费的搞黄视频| 日韩成人av中文字幕在线观看| 日本一二三区视频观看| 国产精品蜜桃在线观看| 中文字幕久久专区| 亚洲电影在线观看av| 久久精品国产自在天天线| 丰满人妻一区二区三区视频av| 女人十人毛片免费观看3o分钟| 久久99热6这里只有精品| 夜夜爽夜夜爽视频| 亚洲成人av在线免费| 成人毛片a级毛片在线播放| 亚洲欧美日韩无卡精品| av又黄又爽大尺度在线免费看| 亚洲精品久久午夜乱码| 人妻少妇偷人精品九色| 久热久热在线精品观看| 成人美女网站在线观看视频| 国产熟女欧美一区二区| www.色视频.com| 舔av片在线| 中文乱码字字幕精品一区二区三区| 高清在线视频一区二区三区| 欧美激情国产日韩精品一区| 两个人的视频大全免费| 亚洲国产精品专区欧美| 97在线视频观看| 国语对白做爰xxxⅹ性视频网站| 成人鲁丝片一二三区免费| 可以在线观看毛片的网站| 成人高潮视频无遮挡免费网站| 国产午夜精品久久久久久一区二区三区| 99久久中文字幕三级久久日本| 超碰av人人做人人爽久久| 亚洲最大成人手机在线| 国产欧美日韩一区二区三区在线 | 18禁裸乳无遮挡免费网站照片| 又爽又黄无遮挡网站| 尤物成人国产欧美一区二区三区| 欧美激情久久久久久爽电影| 日本熟妇午夜| 精品国产露脸久久av麻豆| 久久99蜜桃精品久久| 成人特级av手机在线观看| 国产淫片久久久久久久久| 性插视频无遮挡在线免费观看| 人妻系列 视频| 亚洲国产精品国产精品| 国产精品国产av在线观看| 成人免费观看视频高清| 亚洲av欧美aⅴ国产| 亚洲欧美精品专区久久| 最近最新中文字幕免费大全7| a级毛片免费高清观看在线播放| 嫩草影院新地址| 观看免费一级毛片| 日韩国内少妇激情av| 亚洲精品日韩av片在线观看| 亚洲精品国产色婷婷电影| 久久久久精品性色| 中文在线观看免费www的网站| 国产一区有黄有色的免费视频| 日韩在线高清观看一区二区三区| 亚洲av欧美aⅴ国产| 黑人高潮一二区| 久久综合国产亚洲精品| 欧美成人一区二区免费高清观看| 秋霞在线观看毛片| 日韩电影二区| 日本-黄色视频高清免费观看| 婷婷色麻豆天堂久久| 亚洲国产精品专区欧美| 精品久久国产蜜桃| 久久韩国三级中文字幕| 搡女人真爽免费视频火全软件| 亚洲av免费在线观看| 男插女下体视频免费在线播放| 性色av一级| 亚洲人成网站在线观看播放| 国产午夜福利久久久久久| 99热网站在线观看| 午夜福利在线在线| 久久久久久久午夜电影| 狂野欧美激情性bbbbbb| 亚洲不卡免费看| 国产精品人妻久久久久久| 在线看a的网站| 精品人妻视频免费看| 久久精品综合一区二区三区| 亚洲四区av| 涩涩av久久男人的天堂| 99久久九九国产精品国产免费| 汤姆久久久久久久影院中文字幕| 激情五月婷婷亚洲| 亚洲国产精品国产精品| 男女啪啪激烈高潮av片| 国产精品三级大全| 天天躁日日操中文字幕| 日本一二三区视频观看| 国产精品蜜桃在线观看| 男人添女人高潮全过程视频| 欧美区成人在线视频| 最新中文字幕久久久久| 久久久久久久精品精品| 如何舔出高潮| 制服丝袜香蕉在线| 在线观看美女被高潮喷水网站| 精品99又大又爽又粗少妇毛片| 三级国产精品欧美在线观看| 日日摸夜夜添夜夜爱| 日本一本二区三区精品| 日本午夜av视频| 国产在视频线精品| 国产在线男女| 国产亚洲5aaaaa淫片| 午夜激情久久久久久久| 午夜福利视频1000在线观看| 蜜桃久久精品国产亚洲av| 在线a可以看的网站| 大码成人一级视频| 欧美高清性xxxxhd video| 久久久久久久久久成人| 日韩不卡一区二区三区视频在线| 久久久久国产精品人妻一区二区| 日韩视频在线欧美| 久久热精品热| 国产精品爽爽va在线观看网站| 国产一区有黄有色的免费视频| 国产日韩欧美亚洲二区| 人人妻人人澡人人爽人人夜夜| 天天躁夜夜躁狠狠久久av| 欧美少妇被猛烈插入视频| 男人和女人高潮做爰伦理| 男女无遮挡免费网站观看| 欧美最新免费一区二区三区| 久久国产乱子免费精品| 欧美一区二区亚洲| 亚洲精品第二区| 国产又色又爽无遮挡免| 欧美亚洲 丝袜 人妻 在线| 国产精品爽爽va在线观看网站| 一本久久精品| 97在线人人人人妻| 身体一侧抽搐| 日韩强制内射视频| 亚洲精品中文字幕在线视频 | 夜夜看夜夜爽夜夜摸| 日日撸夜夜添| 91精品伊人久久大香线蕉| 超碰97精品在线观看| 免费少妇av软件| 日韩欧美 国产精品| av一本久久久久| 最后的刺客免费高清国语| 香蕉精品网在线| 国产伦在线观看视频一区| 午夜福利高清视频| 亚洲内射少妇av| 欧美xxxx性猛交bbbb| 国产成人一区二区在线| 永久网站在线| 欧美少妇被猛烈插入视频| 久久99热6这里只有精品| 99久久九九国产精品国产免费| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 高清欧美精品videossex| 观看免费一级毛片| 国产成人精品福利久久| 久久精品久久久久久噜噜老黄| 亚洲精品久久午夜乱码| 欧美极品一区二区三区四区| 少妇的逼水好多| 精品人妻视频免费看| 日韩欧美精品免费久久| 久久精品熟女亚洲av麻豆精品| 亚洲精品久久久久久婷婷小说| 亚洲精品成人av观看孕妇| 中文字幕久久专区| 国产黄色免费在线视频| 激情五月婷婷亚洲| 成年免费大片在线观看| 人人妻人人看人人澡| 国产精品国产三级专区第一集| 美女xxoo啪啪120秒动态图| 91精品国产九色| 97超碰精品成人国产| 两个人的视频大全免费| 久久久国产一区二区| 美女视频免费永久观看网站| 日产精品乱码卡一卡2卡三| 亚洲四区av| 国产精品一区二区三区四区免费观看| av一本久久久久| 日本一本二区三区精品| 日韩伦理黄色片| 色婷婷久久久亚洲欧美| 91狼人影院| av免费在线看不卡| 水蜜桃什么品种好| 日日摸夜夜添夜夜添av毛片| eeuss影院久久| 午夜福利视频1000在线观看| 国产探花在线观看一区二区| 日日啪夜夜撸| 国产精品久久久久久av不卡| 国产探花在线观看一区二区| 亚洲精品,欧美精品| 啦啦啦中文免费视频观看日本| 99热这里只有是精品50| 午夜福利在线在线| 亚洲高清免费不卡视频| 免费黄频网站在线观看国产| 久久久午夜欧美精品| 性插视频无遮挡在线免费观看| 久久久a久久爽久久v久久| 成人毛片60女人毛片免费| 成人特级av手机在线观看| 亚洲成人久久爱视频| 另类亚洲欧美激情| 国产精品国产三级专区第一集| 欧美zozozo另类| 欧美老熟妇乱子伦牲交| 综合色av麻豆| 久久这里有精品视频免费| 亚洲天堂av无毛| xxx大片免费视频| 欧美激情久久久久久爽电影| av免费观看日本| 国产成人a∨麻豆精品| 亚洲精品一区蜜桃| 精品国产一区二区三区久久久樱花 | 国精品久久久久久国模美| 欧美日本视频| 久久久精品欧美日韩精品| 久久精品熟女亚洲av麻豆精品| av免费在线看不卡| 高清毛片免费看| 少妇人妻一区二区三区视频| 国产精品一二三区在线看| 成年女人看的毛片在线观看| 亚洲av电影在线观看一区二区三区 | 成人免费观看视频高清| 美女xxoo啪啪120秒动态图| 久久久精品免费免费高清| 丰满人妻一区二区三区视频av| 欧美成人午夜免费资源| 欧美xxⅹ黑人| 亚洲精品成人久久久久久| 中国国产av一级| 欧美3d第一页| 大陆偷拍与自拍| 国产一区二区三区av在线| 国产成人免费观看mmmm| 欧美激情在线99| 成年人午夜在线观看视频| 小蜜桃在线观看免费完整版高清| 只有这里有精品99| 王馨瑶露胸无遮挡在线观看| 在线观看一区二区三区激情| 汤姆久久久久久久影院中文字幕| 一二三四中文在线观看免费高清| 免费黄网站久久成人精品| 亚洲成人av在线免费| 各种免费的搞黄视频| 成年免费大片在线观看| 国产精品偷伦视频观看了| 97精品久久久久久久久久精品| 九九在线视频观看精品| 欧美国产精品一级二级三级 | 免费大片18禁| av在线亚洲专区| 噜噜噜噜噜久久久久久91| 精品视频人人做人人爽| 又粗又硬又长又爽又黄的视频| 涩涩av久久男人的天堂| 七月丁香在线播放| 极品少妇高潮喷水抽搐| 日韩大片免费观看网站| 国产高清三级在线| 亚洲国产欧美在线一区| 久久精品熟女亚洲av麻豆精品| 日韩大片免费观看网站| 97在线人人人人妻| 18禁裸乳无遮挡动漫免费视频 | 日本猛色少妇xxxxx猛交久久| 日本一本二区三区精品| 午夜激情久久久久久久| 99久久精品一区二区三区| 国产亚洲91精品色在线| 中文字幕久久专区| 男女边吃奶边做爰视频| 免费黄频网站在线观看国产| 亚洲精品久久久久久婷婷小说| 午夜老司机福利剧场| 日韩国内少妇激情av| 成人特级av手机在线观看| 久久久精品免费免费高清| 有码 亚洲区| 秋霞伦理黄片| 久久久精品欧美日韩精品| 成人午夜精彩视频在线观看| 久久久精品94久久精品| 免费av观看视频| 男女边吃奶边做爰视频| 欧美老熟妇乱子伦牲交| 激情 狠狠 欧美| 午夜免费男女啪啪视频观看| av一本久久久久| 午夜免费男女啪啪视频观看| 一级毛片久久久久久久久女| 精品少妇黑人巨大在线播放| 91在线精品国自产拍蜜月| av专区在线播放| 欧美另类一区| 色5月婷婷丁香| 搡老乐熟女国产| 熟女av电影| 久久久精品欧美日韩精品| 中国美白少妇内射xxxbb| 舔av片在线| 人妻一区二区av| 日韩成人av中文字幕在线观看| 日韩国内少妇激情av| 久久6这里有精品| 国产白丝娇喘喷水9色精品| 久久韩国三级中文字幕| 一个人看的www免费观看视频| 91久久精品国产一区二区成人| 欧美激情久久久久久爽电影| 97热精品久久久久久| 91精品伊人久久大香线蕉| 久久女婷五月综合色啪小说 | 亚洲av中文字字幕乱码综合| 91久久精品国产一区二区三区| 91精品一卡2卡3卡4卡| 一级毛片久久久久久久久女| 国产精品一区二区三区四区免费观看| 亚洲婷婷狠狠爱综合网| 国产av码专区亚洲av| 99久久精品国产国产毛片| 日本免费在线观看一区| 精品人妻熟女av久视频| av播播在线观看一区| 热re99久久精品国产66热6| 亚洲成人中文字幕在线播放| 亚洲精品乱码久久久v下载方式| 亚洲国产精品专区欧美| 久久久久久久亚洲中文字幕| a级毛片免费高清观看在线播放| 国产成人午夜福利电影在线观看| 久久久久久久精品精品| 黄色日韩在线| 99久久精品国产国产毛片| 两个人的视频大全免费| 日韩一区二区三区影片| 日韩伦理黄色片| 日日啪夜夜爽| 日韩av免费高清视频| 波多野结衣巨乳人妻| 嫩草影院精品99| 欧美人与善性xxx| 丝袜美腿在线中文| 国产熟女欧美一区二区| 国产黄色视频一区二区在线观看| 国产真实伦视频高清在线观看| 成人亚洲精品av一区二区| 亚洲欧洲国产日韩| 真实男女啪啪啪动态图| 在线免费十八禁| 国产免费一级a男人的天堂| 久久久久国产精品人妻一区二区| 久久99精品国语久久久| 女的被弄到高潮叫床怎么办| 少妇人妻久久综合中文| 欧美潮喷喷水| 99九九线精品视频在线观看视频| 看十八女毛片水多多多| 一区二区三区免费毛片| 两个人的视频大全免费| 性色av一级| 国内揄拍国产精品人妻在线| 精品久久久噜噜| 国产视频首页在线观看| 不卡视频在线观看欧美| 久久久久久久精品精品| 在线播放无遮挡| 国国产精品蜜臀av免费| 国产亚洲av片在线观看秒播厂| 亚洲欧美日韩无卡精品| 97热精品久久久久久| 亚洲无线观看免费| 有码 亚洲区| 久久久久久久亚洲中文字幕| 免费av毛片视频| 亚洲人成网站在线播| 精品人妻视频免费看| 国产精品国产三级国产专区5o| 欧美精品人与动牲交sv欧美| 亚洲va在线va天堂va国产| 大香蕉久久网| 嫩草影院入口| 国产黄片美女视频| 国产精品久久久久久av不卡| 1000部很黄的大片| 日本黄色片子视频| 婷婷色综合大香蕉| 日本一本二区三区精品| 亚洲av电影在线观看一区二区三区 | 国产精品秋霞免费鲁丝片| 国产日韩欧美在线精品| 99九九线精品视频在线观看视频| 波野结衣二区三区在线| 熟女人妻精品中文字幕| 97人妻精品一区二区三区麻豆| 精品人妻一区二区三区麻豆| 欧美少妇被猛烈插入视频| 久久鲁丝午夜福利片| 一级毛片我不卡| 国产永久视频网站| 美女视频免费永久观看网站| 免费不卡的大黄色大毛片视频在线观看| 久久这里有精品视频免费| 亚洲精品久久午夜乱码| 免费观看在线日韩| 性色av一级| 99热网站在线观看| 免费看不卡的av| 最近中文字幕高清免费大全6| 国产老妇女一区| 97精品久久久久久久久久精品| 国产一区二区在线观看日韩| 成人综合一区亚洲| 午夜亚洲福利在线播放| 肉色欧美久久久久久久蜜桃 | 搡女人真爽免费视频火全软件| 亚洲三级黄色毛片| 久久久久久久精品精品| 国产午夜福利久久久久久| 国产精品一二三区在线看| 菩萨蛮人人尽说江南好唐韦庄| 男人添女人高潮全过程视频| 久久99热这里只频精品6学生| 国产有黄有色有爽视频| 国产视频内射| 日韩伦理黄色片| 男人添女人高潮全过程视频| 91在线精品国自产拍蜜月| 亚洲第一区二区三区不卡| 青春草国产在线视频| 国产精品秋霞免费鲁丝片| 久久久精品免费免费高清| 亚洲aⅴ乱码一区二区在线播放| 国内精品美女久久久久久| 亚洲精品中文字幕在线视频 | 制服丝袜香蕉在线| 成人鲁丝片一二三区免费| 亚洲一级一片aⅴ在线观看| 中国美白少妇内射xxxbb| 亚洲国产最新在线播放| 欧美亚洲 丝袜 人妻 在线| 97在线人人人人妻| 久久影院123| .国产精品久久| 在线免费观看不下载黄p国产| 久久精品久久久久久噜噜老黄| 久久精品熟女亚洲av麻豆精品| 久久精品久久精品一区二区三区| 国产精品一区二区性色av| 啦啦啦啦在线视频资源| 亚洲av免费在线观看| 国产精品99久久久久久久久| 日本欧美国产在线视频| 亚洲欧美日韩无卡精品| 精品久久久噜噜| 久久精品久久精品一区二区三区| 我的老师免费观看完整版| 男人狂女人下面高潮的视频| 在线观看三级黄色| 日韩国内少妇激情av| 中文乱码字字幕精品一区二区三区| 亚洲成人中文字幕在线播放| 日韩电影二区| 91精品伊人久久大香线蕉| 天天躁夜夜躁狠狠久久av| 丝瓜视频免费看黄片| 一级毛片电影观看| 99精国产麻豆久久婷婷| 久久精品久久久久久久性| 成人国产av品久久久| 校园人妻丝袜中文字幕| 午夜福利在线观看免费完整高清在| 婷婷色av中文字幕| 久久久国产一区二区| 女人久久www免费人成看片| 国产黄色视频一区二区在线观看| 久久久成人免费电影| 中国美白少妇内射xxxbb| 国产亚洲精品久久久com| 哪个播放器可以免费观看大片| 国产探花在线观看一区二区| 一个人观看的视频www高清免费观看| 男女边摸边吃奶| 涩涩av久久男人的天堂| 午夜亚洲福利在线播放| 国产高潮美女av| 国产欧美日韩一区二区三区在线 | 天天一区二区日本电影三级| 久久久久精品性色| 日韩制服骚丝袜av| 亚洲欧美成人精品一区二区| 春色校园在线视频观看| 精华霜和精华液先用哪个| 可以在线观看毛片的网站| 好男人视频免费观看在线| 国产黄频视频在线观看| 国产日韩欧美亚洲二区| 天美传媒精品一区二区| 交换朋友夫妻互换小说| 亚洲精品一区蜜桃| 亚洲av一区综合| 边亲边吃奶的免费视频| 久久精品国产鲁丝片午夜精品| 偷拍熟女少妇极品色| 三级国产精品欧美在线观看| 亚洲国产成人一精品久久久| 91久久精品国产一区二区三区| 夜夜看夜夜爽夜夜摸| 亚洲国产成人一精品久久久| 亚洲精品国产成人久久av| 别揉我奶头 嗯啊视频| 在线观看国产h片| 高清午夜精品一区二区三区| 晚上一个人看的免费电影| 精品一区在线观看国产| 肉色欧美久久久久久久蜜桃 | 国产男人的电影天堂91| 亚洲在久久综合| 亚洲人与动物交配视频| 18禁裸乳无遮挡动漫免费视频 | 亚洲精品日本国产第一区| 少妇人妻 视频| 高清视频免费观看一区二区|