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

    基于應變響應統(tǒng)計特征的海洋立管損傷診斷方法

    2017-02-14 09:26:31武文華歐進萍
    振動與沖擊 2017年1期
    關鍵詞:立管參量小波

    任 鵬, 周 智,, 武文華, 歐進萍

    (1.大連理工大學 土木工程學院,遼寧 大連 116024; 2.大連理工大學 海岸和近海工程國家重點實驗室,遼寧 大連 116023;3. 大連理工大學 工程力學系,工業(yè)裝備結構分析國家重點實驗室,遼寧 大連 116024)

    基于應變響應統(tǒng)計特征的海洋立管損傷診斷方法

    任 鵬1, 周 智1,2, 武文華3, 歐進萍2

    (1.大連理工大學 土木工程學院,遼寧 大連 116024; 2.大連理工大學 海岸和近海工程國家重點實驗室,遼寧 大連 116023;3. 大連理工大學 工程力學系,工業(yè)裝備結構分析國家重點實驗室,遼寧 大連 116024)

    考慮海洋立管的實際服役環(huán)境和現(xiàn)有監(jiān)測技術水平,發(fā)展一種基于應變響應統(tǒng)計特征的損傷診斷方法。推導以應變測量單元為基礎的細長梁式結構應變響應與結構模態(tài)參數(shù)的映射關系,提出連續(xù)測試數(shù)據(jù)的均方根應變作為海洋立管的損傷特征參量。針對損傷程度較小或量測噪聲影響較大可能引起的損傷誤判,根據(jù)小波多分辨率分析理論,利用小波去噪和分解處理該損傷特征參量的空間域數(shù)據(jù),進而實現(xiàn)準確地損傷定位。在此基礎上,進行了某張力腿式平臺頂部張緊式立管的數(shù)值模擬;依據(jù)損傷診斷策略對假定布設應變傳感器的上部500個應變測量單元分析得到上述損傷特征參量。結果表明:該方法對立管局部剛度下降進行了有效的識別,對單、多處損傷的損傷判別和定位的效果良好,損傷診斷效果受量測噪聲以及服役工況和頂張力等因素變化的影響較小,可為相關工程結構健康監(jiān)測提供理論和技術支持。

    海洋立管;損傷診斷;應變響應;統(tǒng)計特征;小波多分辨率分析

    海洋立管是海洋油氣鉆探和輸送的管線,是海洋平臺結構的關鍵功能設施。出于減少失效風險和增加生產(chǎn)效益的目的,國際各大油氣公司均非常重視對海洋立管、尤其是惡劣服役環(huán)境下的深海立管進行全壽命服役期間的完整性管理。結構健康監(jiān)測技術為立管完整性管理實時提供荷載、響應和損傷的原位測試數(shù)據(jù),為立管生產(chǎn)作業(yè)和風險分析提供及時反饋[1-2]。近年來,由于經(jīng)濟性與可靠性的提升,從英國健康與安全管理部門出臺柔性立管完整性監(jiān)測指導[3],到歐美各大石油公司在墨西哥灣、北海、西非以及巴西臨近海域進行的大量實踐,結構健康監(jiān)測技術已在立管系統(tǒng)設計驗證和安全預警方面取得了應用。其中,應變傳感技術的(準)分布式、各測點同步、在線連續(xù)測試以及直接關聯(lián)疲勞累積等優(yōu)勢,使其逐漸成為立管監(jiān)測的首選:圖1(a)為基于線性差動位移計(Linear Variable Differential Transformer,LVDT)的西非Akpo區(qū)塊立管應變監(jiān)測系統(tǒng)[4];圖1(b)為基于準分布式光纖光柵(Fiber Bragg Grating,FBG)的墨西哥灣Tahiti區(qū)塊應變/曲率監(jiān)測系統(tǒng)[5]。然而,作為健康監(jiān)測的目的之一,實現(xiàn)海洋立管與陸上結構目標一致的損傷判別與定位還比較困難,難以對其進行全壽命的安全評估和壽命預測[6]。缺乏面向實際服役環(huán)境以及匹配現(xiàn)有監(jiān)測技術水平的損傷診斷方法是重要的原因,且相關研究報道比較少見。

    圖1 海洋立管應變監(jiān)測的工程實例Fig.1 Engineering examples of strain monitoring for risers

    應變是對結構局部損傷演化十分敏感的物理量,結合應變傳感技術的損傷診斷方法是結構健康監(jiān)測的重點發(fā)展方向之一。GLISIC等[7]探討了基于長標距光纖光柵和布里淵分布式光纖傳感器獲取應變分布從而進行損傷探測的有效性并且實施了埋地管道監(jiān)測。作為結構振動系統(tǒng)的固有屬性而不依賴于荷載因素,與位移模態(tài)一樣,應變模態(tài)是同一變形能平衡狀態(tài)的另一種表現(xiàn)形式,同時還可作為高敏感性的損傷特征參量。李德葆等[8-9]系統(tǒng)地闡明了應變模態(tài)與位移模態(tài)的對應關系以及實驗應變模態(tài)分析方法。吳智深等[10-11]考慮長標距光纖光柵獲取應變響應并推導出長標距應變模態(tài)向量及其衍生損傷指標,發(fā)展了基于實測應變響應的細長梁式結構損傷診斷方法。在油氣管道方面,提出了相關的損傷診斷方法,主要包括:管道應變模態(tài)變化率的無損檢測法[12]、管道腐蝕檢測的應變敏度比法[13-14]與基于應變模態(tài)差和神經(jīng)網(wǎng)絡的管道損傷診斷[15];盡管上述方法均表現(xiàn)出對管道結構局部損傷的高敏感性,但尚未考慮如何以及能否通過實測響應獲取上述損傷特征參量。

    海洋立管的實際服役環(huán)境比起陸地管道更為復雜,海洋環(huán)境風、浪和海流等形成的外部激勵通常作為平穩(wěn)窄帶隨機過程,這與一般模態(tài)參數(shù)識別所需的白噪聲激勵條件相矛盾。另外,上部平臺運動、頂部張緊力和海床接觸等邊界條件的時變性以及立管作為細長柔性結構的幾何非線性造成模態(tài)類損傷特征參量獲取的困難。因此,本文從現(xiàn)有監(jiān)測技術能夠獲取的應變響應本身出發(fā),將連續(xù)測量數(shù)據(jù)的統(tǒng)計特征應用在損傷診斷方法中。建立均方根應變與結構模態(tài)參數(shù)的映射關系,將其作為海洋立管的損傷特征參量進行損傷的判別和定位。該方法亦能夠避免由于立管結構模型不確定因素給損傷診斷帶來的影響。

    小損傷與量測噪聲導致不顯著的損傷特征,文中根據(jù)小波多分辨率分析及其去噪原理,提取損傷特征參量的突變成分用以損傷定位。模擬某張力腿式平臺頂部張緊式生產(chǎn)立管的應變響應,分別針對兩類服役工況進行了單處和多處損傷的診斷,驗證了本文提出方法的有效性。

    1 基于應變響應統(tǒng)計特征的損傷診斷方法

    1.1 應變響應的統(tǒng)計特征

    對于細長梁式結構,離散單元的劃分一般與應變傳感器的測量標距相對應。文中稱此類依照幾何標距規(guī)定其尺度的單元為應變測量單元。

    根據(jù)模態(tài)疊加原理,n自由度彈性結構節(jié)點位移u

    (1)

    由應變與位移的幾何關系,得到

    (2)

    由模態(tài)解耦后的頻域分析可知

    (3)

    模態(tài)解耦后的時域分析可得杜哈梅爾卷積積分

    (4)

    (5)

    將式(4)代入式(5),得到

    (6)

    式(6)的具體推導過程參考文獻[16],式中RPrs(τ)為廣義力的自、互相關函數(shù),θ1和θ2為虛時間變量。

    上述平穩(wěn)隨機過程可實現(xiàn)兩對傅里葉變換,包括:功率譜密度Sε(ω)與自相關函數(shù)Rε(τ)、頻率響應函數(shù)Hr(ω)與脈沖響應函數(shù)hr(t)。將式(6)代入傅里葉變換對(詳見文獻[16]),得到

    (7)

    針對第j個應變測量單元的響應自譜密度,式(7)變?yōu)?/p>

    (8)

    利用上式可將各階模態(tài)的獨立貢獻和耦合影響分開

    (9)

    (10)

    1.2 損傷特征參量

    對于低阻尼體系,根據(jù)模態(tài)坐標系統(tǒng)振型互相獨立的性質,各階模態(tài)的響應是統(tǒng)計無關的,耦合影響可以忽略不計,此時式(9)和式(10)等式右邊交叉項為零[16]。

    如響應為平穩(wěn)隨機過程,如果滿足以下兩項識別條件,式(9)可用于峰值提取法識別模態(tài)參數(shù)[17]:

    二為隨機激勵為寬頻帶的近似白噪聲;對于激勵卓越周期處存在峰值的輸入功率譜,認為其能量峰值不能與結構共振范圍的峰值相比;

    驗證上述條件的方法:結構響應點自譜密度峰值曲線,與擬識別第r階振型形狀處于同一相位或相差180°,否則該峰值曲線對應激勵能量峰值。

    (11)

    式中:SPk(ω)為結構離散作用點的海洋環(huán)境荷載譜,如荷載譜已知,可用數(shù)值積分求解應變響應方差。另一方面,應變響應方差也可通過立管應變測量單元實測獲取[18];則上式可仿照隱式的靜態(tài)問題,即

    (12)

    基于上述推導過程,定義應變測量單元j無量綱均方根應變σεj(Root-Mean-Squared Strain, RMSS)為海洋立管的結構損傷特征參量。該損傷特征參量具有如下特性:

    1)當海洋環(huán)境荷載作為平穩(wěn)(或分段平穩(wěn))隨機過程處理時,立管振動響應中包含較多反映結構狀態(tài)的信息。式(11)建立的RMSS與結構模態(tài)參數(shù)的映射關系中,應變模態(tài)振型系數(shù)最為敏感。

    2)除海洋環(huán)境荷載外,RMSS僅與某一應變測量單元j的應變、位移模態(tài)振型系數(shù)有關;應變傳感器布置于易損區(qū)域即可,不用獲取全局性的模態(tài)振型。

    3)上部平臺運動、頂部張緊力和海床接觸等邊界條件改變導致立管軸向應變等發(fā)生變化,可通過剔除連續(xù)實測應變數(shù)據(jù)的統(tǒng)計平均趨勢,保留應變響應的脈動分量,從而忽略立管結構模型不確定因素給損傷診斷帶來的影響。

    4)根據(jù)式(11),RMSS的獲取與無損傷狀態(tài)的結構基準模型無關,且適合在線處理。

    1.3 基于小波多分辨率分析的改進

    1988年,MALLAT在構造正交小波基的工作中提出了多分辨率(尺度)分析的概念,其框架為基于同一空間不同尺度分析信號。對滿足一定條件的空間集合{Vl,l∈Z}存在尺度函數(shù)序列αl構成標準正交基,通過構建正交于尺度空間的小波空間{Wl,l∈Z}形成標準正交小波基函數(shù)βl。由于空間數(shù)目為無限可數(shù)的,因此多分辨率分析使得所感興趣的信號特征即使在某一分辨率檢測不到,在另一個分辨率卻很容易觀察處理,使得小波變換成為真正的時頻局部化方法。

    在MALLAT塔式算法中,使用分別等價于低通和高通濾波器的尺度和小波濾波系數(shù),從最高分辨率的尺度空間開始,根據(jù)雙尺度方程依次構造下一尺度的尺度函數(shù)和小波函數(shù)。該算法使得構造出的尺度基函數(shù)具有信號平均的成分,其對應的小波變換系數(shù)稱為近似系數(shù)Wla;小波基函數(shù)具有信號變化的成分,其對應的小波變換系數(shù)稱為細節(jié)系數(shù)Wld。從而原始信號s經(jīng)l尺度分解為小波近似成分a和小波細節(jié)成分d

    (13)

    本文涉及小波多分辨率分析以下三方面的應用

    1)信號平穩(wěn)處理:當獲取應變響應數(shù)據(jù)的統(tǒng)計平均變化較大時,采用小尺度范圍進行小波分解和重構,使得響應信號為脈動信號而趨于平穩(wěn)。

    2)信號去噪:實測應變響應不可避免地受到量測噪聲的影響,同時RMSS本身作為響應統(tǒng)計特征也具有一定的波動性。對于程度較小的損傷來說,可能不具有較為明顯的突變;此時認為損傷引起的空間域信號變化是振動系統(tǒng)物理參數(shù)變化的表征,屬于有用信號,可在某一尺度將其成分體現(xiàn)出來,所以采用分層閾值去噪方法,將無用的噪聲信號濾除,保留包含損傷的信號成分。

    (14)

    3)信號細節(jié)分析:對損傷特征參量RMSS處理后的空間域信號進行小波分解,提取其細節(jié)成分,從而放大現(xiàn)有信號的突變。

    1.4 損傷診斷策略

    1)依據(jù)應變傳感器布置獲取應變測量單元的應變響應數(shù)據(jù);由于海洋立管同時存在彎曲和拉伸變形,應變傳感器一般沿管壁環(huán)向布置,即采用90°或120°夾角方向的布設方案,如圖2所示;管壁應變還須考慮雙向應力狀態(tài)橫向效應的修正,具體換算方法詳見文獻[20-21],不重復贅述;

    圖2 應變傳感器沿管道環(huán)向布置示意圖Fig.2 Schematic diagram of hoop-deployed strain sensors

    2)提取連續(xù)測量的平穩(wěn)時間段應變響應;如響應的統(tǒng)計平均變化較大,采用小波分析進行處理,剔除統(tǒng)計平均趨勢,得到響應的脈動分量;

    3)直接計算各應變測量單元上述應變響應數(shù)據(jù)的功率譜密度Sεj(ω);

    5)通過RMSS離散點連接曲線的突變處判別損傷的出現(xiàn)和確定損傷的位置;

    6)當損傷程度較小或量測噪聲影響較大造成難以人為觀察突變時,采用小波去噪的方法,即通過小波分解、對分解后的細節(jié)系數(shù)進行閾值量化處理和重構等步驟保留包含損傷特征的空間域信號RMSSnew;此時式(14)中相關參數(shù)l,m和λ需重點考慮;

    7)對小波去噪處理后的RMSS空間域信號分解,利用其細節(jié)成分進行損傷定位。

    2 數(shù)值案例

    2.1 立管模型與模擬工況

    目標立管為南海海域某張力腿式平臺(TLP)頂部張緊式生產(chǎn)立管(Top Tensioned Riser,TTR)[22],立管模型參數(shù)如表1所示。根據(jù)典型TTR立管的設計分析結果,目標立管上部受波浪和上部浮體運動影響較大可視為易損區(qū)域。假設立管上部500 m布置標距為1 m的應變傳感器,并通過尺度為1 m的應變測量單元獲取得到應變響應。

    由于本文僅從損傷診斷方法的驗證方面考慮,假設波浪和海流為同向來流,應變傳感器布置于來流方向平面內,獲取應變響應ε(t),實際工況下依據(jù)圖2所示應變傳感器布置及其各應變換算關系,即可得到該主方向的ε(t)。TLP平臺與TTR立管的二維簡化模型如圖3所示。

    表1 立管模型參數(shù)

    圖3 TLP平臺及其TTR立管模型示意圖Fig.3 Schematic diagram of the TLP platform and its TTR riser

    考慮有限元數(shù)值模擬所需環(huán)境荷載資料有限,應盡量合理地設置TTR立管的服役工況。本文針對環(huán)境荷載和頂張力改變設置兩類工況驗證提出的損傷診斷策略。立管頂部張緊力分別設為400和1 000 kN。工況一中有義波高7.0 m,譜峰周期12.1 s;工況二有義波高13.8 m,譜峰周期16.1 s。拖拽力系數(shù)取1,慣性力系數(shù)取1.5。海流環(huán)境條件近似為表2。另外,假設上部平臺運動xTLP主要由二階慢漂力引起,并作用在平面內的縱蕩方向,慢漂周期為150 s。

    表2 兩類工況下的海流流速(m/s)

    隨機波浪被視為脈動隨機過程,具有各態(tài)歷經(jīng)性。不同海況的波浪頻譜(或稱能量譜)計算得到具有統(tǒng)計特性的有義波高和跨零周期(或譜峰周期);選擇合適的波浪理論,如線性Airy波或Stokes五階波理論等,通過求得的波浪運動的速度勢,計算流場中水質點的速度和加速度;在此基礎上,通過Morison公式可計算波浪對立管的環(huán)境荷載,如拖拽力和慣性力等。海流荷載計算方法同理。依據(jù)通用有限元軟件Abaqus提供的Aqua模塊進行定義加載,進而通過直接積分法模擬立管的動力響應。

    2.2 結果分析

    圖4 前兩階模態(tài)振型和功率譜峰值曲線的對比Fig.4 Comparisons of the first two mode shapes and peaks of PSDs

    如圖4為海洋立管服役工況一前兩階模態(tài)振型與應變測量單元(采用Abaqus的B21梁元模擬)響應功率譜密度前兩階峰值曲線的無量綱對比。可以看出,由于海洋環(huán)境荷載為窄帶隨機過程,立管結構響應受到激勵能量驅動,難以識別出模態(tài)參數(shù)。因此,本文采用基于響應統(tǒng)計特性的RMSS進行損傷判別與定位,并通過以下?lián)p傷狀態(tài)驗證其損傷診斷效果。

    圖5為根據(jù)第300號單元三種損傷狀態(tài)的100 s應變響應數(shù)據(jù)得到的損傷特征參量RMSS。其中,圖5(a)為服役工況一,單元剛度下降20%;圖5(b)為服役工況改變,單元剛度仍下降20%;圖5(c)為服役工況二,單元質量增加20%。可假定單元剛度下降為管壁裂紋造成,單元質量增加由海生物附著引起。結果說明:當服役環(huán)境發(fā)生改變時,只需采集一定時間段的平穩(wěn)隨機過程或脈動分量,即可使用RMSS進行損傷診斷;該損傷特征參量對剛度改變非常敏感并且不受到質量改變的影響,上述特性可以避免剛度和質量同時變化造成的損傷誤判。

    圖5 立管單處損傷狀態(tài)的損傷特征參量Fig.5 RMSS of the riser with a single damage

    圖6 立管單處損傷10%的RMSSFig.6RMSSoftheriserwithasingledamageof10%圖7 RMSS的小波去噪Fig.7Waveletde-noisingofRMSS圖8 RMSS小波去噪后的細節(jié)分解Fig.8DetaildecompositionofRMSSafterwaveletde-noising

    為進一步驗證提出方法的可行性,面向服役工況二,在第200、300和400號單元分別預設20%、10%和15%三處損傷;與此同時,對原應變響應數(shù)據(jù)加入10%高斯白噪聲,作為量測噪聲的模擬。圖9(a)和(b)分別為未加入和加入噪聲后損傷診斷效果,可見不同RMSS值與損傷程度變化趨勢相同,但量測噪聲對損傷特征參量影響較大。因此,根據(jù)上述損傷診斷策略中f和g兩步,進行小波去噪和分解,得到的RMSS改進如圖10所示。本文提出的方法能夠由量測噪聲影響下的應變測試數(shù)據(jù)對海洋立管的多處損傷進行定位。

    圖9 立管多處損傷的RMSSFig.9RMSSoftheriserwithmultipledamages圖10 量測噪聲影響下立管多處損傷的RMSSFig.10RMSSoftheriserwithmultipledamagesunderGausswhitemeasurementnoises

    3 結 論

    本文面向海洋立管實際服役環(huán)境和現(xiàn)有監(jiān)測技術水平,進行了細長梁式結構應變響應統(tǒng)計特征的理論推導,建立了應變響應與結構模態(tài)參數(shù)的映射關系,提出了應變測量單元均方根應變RMSS作為海洋立管的損傷特征參量;基于小波多分辨率分析對RMSS進行了改進,依據(jù)提出的損傷診斷策略實施例了目標立管的數(shù)值驗證。得到主要結論如下:

    (1)針對布設了應變傳感器的立管上部500個測量單元,單、多處損傷的損傷判別和定位的效果良好,RMSS的獲取無須考慮未損傷狀態(tài)的測試數(shù)據(jù)且適合在線處理。

    (2)依據(jù)損傷診斷策略對模擬的立管局部剛度下降進行了有效的識別,RMSS僅與應變測量單元處的應變響應數(shù)據(jù)相關。

    (3)對損傷程度較小或量測噪聲影響較大的情況,可應用小波多分辨率分析處理RMSS空間域數(shù)據(jù),進而實現(xiàn)準確地損傷定位,避免損傷誤判。

    本文方法或可推廣應用到承受脈動激勵的柔性結構,如索體、主桁架、高層建筑和輸電線塔等。

    [ 1 ] COOK H H,DOUGLAS E. Riser integrity management for deepwater developments [C]//Proceedings of Offshore Technology Conference,2006,Paper No. OTC-7891.

    [ 2 ] DNV. Riser Integrity Management [S]. Recommended Practice F206,2008.

    [ 3 ] HSE. Guidelines for Integrity Monitoring of Unbond Flexible Pipe [S]. Offshore Technology Report,OTO-98019,1998.

    [ 4 ] DOUARON S L,VUATTIER C,PERROMAT V. Akpo riser integrity monitoring system design, deployment, commissioning and start-up [C]//Proceedings of Deep Offshore Technology International Conference, New Orleans, Louisiana, 2009.

    [ 5 ] KARAYAKA M, CHEN J, RUF W. Tahiti online monitoring system for steel catenary risers and flowlines [C]//Proceedings of Offshore Technology Conference,2009,Paper No. OTC-19860.

    [ 6 ] GOLDSMITH B,F(xiàn)OYT E,HARIHARAN M. The role of offshore monitoring in an effective deepwater riser integrity management program [C]//Proceedings of 26th International Conference on offshore Mechanics and Arctic Engineering,San Diego,2007: 881-886.

    [ 7 ] GLISIC B, HUBBELL D, SIGURDARDOTTIR D H,et al. Damage detection and characterization using long-gauge and distributed fiber optic sensors [J]. Optical Engineering, 2013, 52(8): 1-12.

    [ 8 ] 李德葆,諸葛鴻程,王波. 實驗應變模態(tài)分析原理和方法 [J]. 清華大學學報(自然科學版),1990,30(2): 105-112. LI Debao,ZHUGE Hongcheng,WANG Bo. On the principle and technique of experimental strain modal analysis [J]. Journal of Tsinghua University(Science and Technology),1990,30(2): 105-112.

    [ 9 ] 李德葆. 實驗應變-應力模態(tài)分析若干問題的進展評述 [J]. 振動與沖擊,1996,15(1): 13-11. LI Debao. A general review on several fundamental points of experimental strain/stress modal analysis [J]. Journal of Vibration and Shock,1996,15(1): 13-11.

    [10] ADEWUYI A P,WU Z S,SERKER K. Strain measurements assessment of vibration-based damage identification methods using displacement and distributed [J]. Structural Health Monitoring-An International Journal, 2009, 8(6): 443-461.

    [11] 吳智深,李素貞,ADEWUYI A P. 基于應變分布響應的模態(tài)分析理論與應用 [J]. 科技導報,2010,28(8): 94-103. WU Zhishen,LI Suzhen,ADEWUYI A P. Modal analysis based on distributed strain measurements: theory and application [J]. Science and Technology Review,2010,28(8): 94-103.

    [12] 任權,李洪升,郭杏林. 基于應變模態(tài)變化率的壓力管道無損檢測[J]. 大連理工大學學報,2001,41(6): 648-652. REN Quan,LI Hongsheng,GUO Xinglin. Undamaged detection of pipe based on modal varied quotiety of strain [J]. Journal of Dalian University of Technology,2001,41(6): 648-652.

    [13] 李洪升,楊智榮,李亞民,等. 壓力管道腐蝕災變檢測的應變敏度比法 [J]. 工程力學,2007,24(3): 165-169. LI Hongsheng,YANG Zhirong,LI Yamin,et al. Detection of the corrosion damage in pressure pipeline using srain sensitivity ratio method [J]. Engineering Mechanics,2007,24(3): 165-169.

    [14] 許成祥,賈善坡,涂金釗,等. 基于工作應變模態(tài)的管道損傷識別試驗研究 [J]. 西南交通大學學報,2013,48(6): 1031-1037. XU Chengxiang, JIA Shanpo, TU Jinzhao, et al. Experimental study on damage identification of pipeline based on operational strain modal shape [J]. Journal of Southwest Jiaotong University,2013,48(6): 1031-1037.

    [15] 周邵萍,郝占峰,韓紅飛,等. 基于應變模態(tài)差和神經(jīng)網(wǎng)絡的管道損傷識別 [J]. 振動、測試與診斷,2015,35(2): 334-338. ZHOU Shaoping,HAO Zhanfeng,HAN Hongfei, et al. Damage identification in Straigh pipeline using strain modal difference and network [J]. Journal of Vibration, Measurement, and Diagnosis,2015,35(2): 334-338.

    [16] (美)克拉夫R,彭津J. 王光遠(譯). 結構動力學 [M]. 北京:高等教育出版社,2006: 339-345.

    [17] 寶志雯,陳志鷗. 從建筑物的脈動響應確定其動力特性 [J]. 深圳大學學報(理工版),1986,1(1): 36-49. BAO Zhiwen, CHEN Zhiou. Experimental determination of dynamic characteristics of a structure by means of its ambient response signals [J]. Journal of Shenzhen University(Science and Engineering), 1986,1(1): 36-49.

    [18] WANG D S, XIANG W, ZHU H P. Damage identification in beam type structures based on statistical moment using a two step method [J]. Journal of Sound and Vibration, 2014, 333: 745-760.

    [19] DONOHO D L,JOHNSTONE I M. Ideal spatial adaptation by wavelet shrinkage [J]. Biometrika,1994,81(3): 425-455.

    [20] 許學瑞,帥健,肖偉生. 滑坡多發(fā)區(qū)管道應變監(jiān)測應變計安裝方法 [J]. 油氣儲運,2010,29(10): 780-784. XU Xuerui, SHUAI Jian, XIAO Weisheng. The installation method of strain gage used to monitor pipeline strain in landslide-prone areas [J].Oil & Gas Storage and Transportation,2010,29(10): 780-784.

    [21] 李艷萍,冼國棟,王維斌,等. 管道應力應變監(jiān)測試驗規(guī)程研究 [J]. 油氣儲運,2008,27(10): 21-23. LI Yanping, XIAN Guodong, WANG Weibin, et al. Experimental specification research on pipeline stress and strain monitoring [J], Oil & Gas Storage and Transportation,2008,27(10): 21-23.

    [22] 王世圣,謝彬,李新仲. 在南海環(huán)境條件下深水典型TLP的運動響應分析 [J]. 中國造船,2011,52(增刊1): 94-101. WANG Shisheng, XIE Bin, LI Xinzhong. The motion response analysis of deep water typical TLP in environment conditions of South China Sea [J]. Ship Building of China,2011,52(Sup1): 94-101.

    Damagediagnosis strategy for offshore risers based on statistical features of strain responses

    REN Peng1, ZHOU Zhi1,2, WU Wenhua3, OU Jinping2

    (1. School of Civil Engineering, Dalian University of Technology, Dalian 116024, China;2. State Key Laboratory of Coastal and Offshore Engineering, Dalian University of Technology, Dalian 116023, China;3.State Key Laboratory of Structural Analysis for Industrial Equipment, Dalian University of Technology, Dalian 116024, China)

    Considering the service environment and the monitoring technology level of offshore risers, a damage diagnosis method was developed based on statistical characteristics of strain responses. A mapping relation between modal parameters and strain responses of a strain-measured element was derived for a beam-like structure. A root mean square strain based on a set of test data was proposed as the damage characteristic parameter of offshore risers. Aiming at a plamage misjudging due to a lower level of damage or a larger measurement noise, the data were processed based on the wavelet de-noising and decomposition. Furthermore, a TTR on a TLP was simulated to obtain the damage characteristic parameter through analyzing 500 strain-measured elements at the top of the TTR. The results showed that the method can be used to identify the local stiffness dropping of the TTR effectively althrough the damage diagnosis is interfered by the measurement noise and other uncertain factors; the damage diagnosis strategy can provide a support for structural health monitoring both theoretically and technically.

    offshore riser; damage diagnosis; strain response; statistical characteristics; wavelet multi-resolution

    國家重點基礎研究發(fā)展計劃(973計劃)(2011CB013705);國家自然科學基金自由申請項目(15572072)

    2015-09-07 修改稿收到日期:2016-12-31

    任鵬 男,博士生,1984年12月生

    周智 男,博士,教授,1973年7月生 E-mail:zhouzhi@dlut.edu.cn

    TU311.3;TE58

    A

    10.13465/j.cnki.jvs.2017.01.004

    猜你喜歡
    立管參量小波
    構造Daubechies小波的一些注記
    科技風(2021年19期)2021-09-07 14:04:29
    基于MATLAB的小波降噪研究
    電子制作(2019年13期)2020-01-14 03:15:32
    常見高層建筑物室內給水立管材質解析
    基于改進的G-SVS LMS 與冗余提升小波的滾動軸承故障診斷
    深水鋼懸鏈立管J型鋪設研究
    The Power of Integration
    Beijing Review(2015年43期)2015-11-25 03:12:04
    海洋立管濕模態(tài)振動分析
    環(huán)形光的形成與參量的依賴關系
    物理實驗(2015年9期)2015-02-28 17:36:51
    含雙參量的p-拉普拉斯邊值問題的多重解
    鎖定放大技術在參量接收陣中的應用
    聲學技術(2014年2期)2014-06-21 06:59:14
    国产熟女xx| 又黄又爽又免费观看的视频| 无限看片的www在线观看| 亚洲国产欧美网| 免费看光身美女| 久久久久国产精品人妻aⅴ院| 亚洲欧美精品综合一区二区三区| 亚洲成人免费电影在线观看| 我的老师免费观看完整版| www.自偷自拍.com| 欧美又色又爽又黄视频| 婷婷精品国产亚洲av在线| 此物有八面人人有两片| 亚洲片人在线观看| 91麻豆精品激情在线观看国产| 亚洲av免费在线观看| 国产亚洲欧美在线一区二区| 成年女人永久免费观看视频| 我要搜黄色片| a级毛片在线看网站| 两性夫妻黄色片| 禁无遮挡网站| 香蕉国产在线看| 欧美日韩福利视频一区二区| 女人被狂操c到高潮| 精品久久久久久成人av| 久久久国产成人精品二区| 亚洲精品乱码久久久v下载方式 | 黑人欧美特级aaaaaa片| 九色国产91popny在线| 美女被艹到高潮喷水动态| 国产亚洲精品一区二区www| 国产精品一区二区三区四区免费观看 | 亚洲欧美一区二区三区黑人| 最近最新免费中文字幕在线| 99riav亚洲国产免费| 在线观看免费午夜福利视频| 成人鲁丝片一二三区免费| 国产av一区在线观看免费| 亚洲在线自拍视频| 日韩欧美在线乱码| 老熟妇乱子伦视频在线观看| 999久久久精品免费观看国产| 成人国产一区最新在线观看| 91av网一区二区| 国产精品久久久久久精品电影| 天堂动漫精品| 俄罗斯特黄特色一大片| 俺也久久电影网| 精品无人区乱码1区二区| 中出人妻视频一区二区| 国产精品电影一区二区三区| 一级毛片高清免费大全| 欧美高清成人免费视频www| 亚洲成av人片在线播放无| 亚洲午夜精品一区,二区,三区| 免费在线观看影片大全网站| 午夜影院日韩av| 丰满的人妻完整版| 国产一区二区激情短视频| 亚洲av第一区精品v没综合| 国产精品 国内视频| 禁无遮挡网站| 国产精品美女特级片免费视频播放器 | 久久国产精品人妻蜜桃| 久久久久久久久中文| 老司机在亚洲福利影院| 一区二区三区高清视频在线| 一进一出抽搐gif免费好疼| 久久国产乱子伦精品免费另类| 他把我摸到了高潮在线观看| 国产一区二区在线av高清观看| 免费观看人在逋| 男人的好看免费观看在线视频| 99热6这里只有精品| 美女高潮的动态| 极品教师在线视频| 欧美成人免费av一区二区三区| 中文在线观看免费www的网站| 久久久久久久久中文| 你懂的网址亚洲精品在线观看 | 在现免费观看毛片| 特大巨黑吊av在线直播| 亚洲av中文字字幕乱码综合| 久久久久久久久久久免费av| 日本免费一区二区三区高清不卡| 欧美精品一区二区大全| 午夜福利高清视频| 国产高潮美女av| 欧美又色又爽又黄视频| 亚洲激情五月婷婷啪啪| 成人三级黄色视频| 久久婷婷人人爽人人干人人爱| 3wmmmm亚洲av在线观看| 99久久中文字幕三级久久日本| 菩萨蛮人人尽说江南好唐韦庄 | 最近手机中文字幕大全| 国产亚洲精品久久久com| 久久久久网色| 国产精品一二三区在线看| 久久精品国产亚洲av天美| 精品久久久久久久久久久久久| 联通29元200g的流量卡| 久久这里有精品视频免费| 国产精品99久久久久久久久| 久久精品国产亚洲av涩爱| 人妻制服诱惑在线中文字幕| 黄色一级大片看看| 亚洲经典国产精华液单| 成人亚洲欧美一区二区av| 精品人妻偷拍中文字幕| 欧美bdsm另类| 国产精品久久久久久精品电影| 亚洲18禁久久av| 噜噜噜噜噜久久久久久91| 成人午夜高清在线视频| 好男人视频免费观看在线| 国产在视频线精品| 中文字幕精品亚洲无线码一区| 最近2019中文字幕mv第一页| 午夜激情福利司机影院| 国产精品伦人一区二区| 男人狂女人下面高潮的视频| 亚洲美女搞黄在线观看| 国产av一区在线观看免费| 中文字幕av在线有码专区| 欧美最新免费一区二区三区| 寂寞人妻少妇视频99o| 中文字幕av成人在线电影| 中文欧美无线码| 精品久久久久久久久亚洲| 日韩三级伦理在线观看| 免费观看人在逋| 国产成人a∨麻豆精品| 国产精品嫩草影院av在线观看| 18禁裸乳无遮挡免费网站照片| 国产成人aa在线观看| 女的被弄到高潮叫床怎么办| 久久久久久久亚洲中文字幕| 免费无遮挡裸体视频| 成人毛片60女人毛片免费| 18+在线观看网站| 国产精品国产高清国产av| 全区人妻精品视频| 少妇人妻一区二区三区视频| 国内揄拍国产精品人妻在线| 国产一区二区在线av高清观看| 亚洲欧美清纯卡通| 成年av动漫网址| 黄色欧美视频在线观看| 亚洲经典国产精华液单| 免费看a级黄色片| 国产黄片视频在线免费观看| 成人二区视频| www.av在线官网国产| av播播在线观看一区| 国产国拍精品亚洲av在线观看| 18禁在线无遮挡免费观看视频| 非洲黑人性xxxx精品又粗又长| 国产精品不卡视频一区二区| 成人毛片a级毛片在线播放| 草草在线视频免费看| 狠狠狠狠99中文字幕| 啦啦啦啦在线视频资源| 特级一级黄色大片| 国产av在哪里看| 麻豆国产97在线/欧美| 日日干狠狠操夜夜爽| 日韩一区二区三区影片| 啦啦啦韩国在线观看视频| 三级经典国产精品| 欧美日韩精品成人综合77777| 国产精品爽爽va在线观看网站| 成人无遮挡网站| 中文字幕人妻熟人妻熟丝袜美| 永久网站在线| 美女大奶头视频| 校园人妻丝袜中文字幕| 建设人人有责人人尽责人人享有的 | 日本猛色少妇xxxxx猛交久久| 亚洲av免费高清在线观看| av在线老鸭窝| 久久久精品大字幕| 黄色配什么色好看| 国产三级在线视频| 天美传媒精品一区二区| 成人鲁丝片一二三区免费| 少妇熟女aⅴ在线视频| 国产不卡一卡二| 嫩草影院精品99| 大香蕉97超碰在线| 男人和女人高潮做爰伦理| 国产久久久一区二区三区| 大香蕉久久网| 91aial.com中文字幕在线观看| 日日干狠狠操夜夜爽| 中文字幕免费在线视频6| 联通29元200g的流量卡| 午夜福利网站1000一区二区三区| 免费电影在线观看免费观看| 欧美性感艳星| 亚洲最大成人手机在线| 日本一本二区三区精品| 国产女主播在线喷水免费视频网站 | av天堂中文字幕网| 日韩制服骚丝袜av| 深爱激情五月婷婷| 欧美成人精品欧美一级黄| 黄色配什么色好看| 国内揄拍国产精品人妻在线| 嫩草影院入口| videos熟女内射| 国产成人精品久久久久久| 国产三级中文精品| 亚州av有码| 国产黄色视频一区二区在线观看 | 淫秽高清视频在线观看| a级毛色黄片| 九草在线视频观看| 美女黄网站色视频| 一级毛片我不卡| 国产一区有黄有色的免费视频 | 天天一区二区日本电影三级| 国产精品综合久久久久久久免费| 中文字幕av成人在线电影| 亚洲精品国产成人久久av| 中文字幕精品亚洲无线码一区| 黄色欧美视频在线观看| 欧美+日韩+精品| 22中文网久久字幕| 亚洲真实伦在线观看| 亚洲欧洲日产国产| 亚洲欧美中文字幕日韩二区| 国产精品人妻久久久久久| 亚洲精品日韩av片在线观看| 亚洲aⅴ乱码一区二区在线播放| 极品教师在线视频| 在线a可以看的网站| 亚洲av男天堂| 搡老妇女老女人老熟妇| 国产亚洲一区二区精品| av在线亚洲专区| av在线观看视频网站免费| 成人漫画全彩无遮挡| 麻豆精品久久久久久蜜桃| 亚洲在久久综合| 色综合亚洲欧美另类图片| 日本欧美国产在线视频| 亚洲国产成人一精品久久久| 夫妻性生交免费视频一级片| 国产精品av视频在线免费观看| 亚洲av一区综合| 国产乱人偷精品视频| 纵有疾风起免费观看全集完整版 | 91av网一区二区| 69av精品久久久久久| 婷婷色av中文字幕| 永久网站在线| 亚洲自偷自拍三级| 国产成人aa在线观看| 久久久久免费精品人妻一区二区| 免费播放大片免费观看视频在线观看 | 亚洲国产欧美在线一区| 中文乱码字字幕精品一区二区三区 | 99热这里只有是精品在线观看| 人妻制服诱惑在线中文字幕| 99热网站在线观看| 亚洲国产色片| 欧美一级a爱片免费观看看| 日本免费一区二区三区高清不卡| 日本爱情动作片www.在线观看| 午夜福利成人在线免费观看| 国产熟女欧美一区二区| 亚洲欧洲日产国产| 青春草国产在线视频| 精品久久久久久久人妻蜜臀av| 亚洲精品成人久久久久久| 中文字幕熟女人妻在线| 久久精品国产99精品国产亚洲性色| 精品国产一区二区三区久久久樱花 | 波多野结衣巨乳人妻| 一区二区三区高清视频在线| 视频中文字幕在线观看| 国产精品伦人一区二区| av在线老鸭窝| 两个人的视频大全免费| 一级毛片电影观看 | 欧美成人午夜免费资源| 国内精品一区二区在线观看| 免费观看a级毛片全部| av卡一久久| 精品久久久久久久久久久久久| 免费av不卡在线播放| 色吧在线观看| 麻豆成人午夜福利视频| 麻豆成人av视频| 亚洲av.av天堂| 亚洲国产最新在线播放| 亚洲av中文字字幕乱码综合| 亚洲av成人精品一二三区| 久久久久九九精品影院| 最近最新中文字幕大全电影3| 国产成人精品久久久久久| 久久久久久久午夜电影| 色视频www国产| 白带黄色成豆腐渣| 亚州av有码| 干丝袜人妻中文字幕| 国产成人freesex在线| 中文精品一卡2卡3卡4更新| 在线播放国产精品三级| 日本wwww免费看| 老司机福利观看| 国产极品天堂在线| 成人性生交大片免费视频hd| 麻豆久久精品国产亚洲av| 国产在视频线精品| 你懂的网址亚洲精品在线观看 | 欧美精品一区二区大全| 亚洲av免费高清在线观看| 色网站视频免费| 最近中文字幕高清免费大全6| 久久亚洲国产成人精品v| 久久久久久国产a免费观看| 日本一本二区三区精品| 日韩三级伦理在线观看| 精品久久久噜噜| 色吧在线观看| 久久久久久久亚洲中文字幕| 国产一区亚洲一区在线观看| 国产精品麻豆人妻色哟哟久久 | 又黄又爽又刺激的免费视频.| 草草在线视频免费看| 午夜精品国产一区二区电影 | 国产探花极品一区二区| 免费观看人在逋| 国产亚洲精品av在线| 一级毛片电影观看 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | av国产久精品久网站免费入址| 亚洲欧美日韩东京热| 噜噜噜噜噜久久久久久91| 美女被艹到高潮喷水动态| av又黄又爽大尺度在线免费看 | 久久精品久久久久久噜噜老黄 | 亚洲国产欧美在线一区| 欧美97在线视频| 中文精品一卡2卡3卡4更新| 又黄又爽又刺激的免费视频.| 99热网站在线观看| 十八禁国产超污无遮挡网站| 久久亚洲精品不卡| 久久精品国产亚洲av天美| 最新中文字幕久久久久| 99热这里只有精品一区| 色综合亚洲欧美另类图片| 亚洲不卡免费看| 久久99精品国语久久久| 久久久国产成人精品二区| 亚洲欧美精品专区久久| 超碰97精品在线观看| 欧美日本亚洲视频在线播放| 嘟嘟电影网在线观看| 黄片无遮挡物在线观看| 特大巨黑吊av在线直播| 亚洲美女视频黄频| 中文在线观看免费www的网站| 亚洲伊人久久精品综合 | 国产精品无大码| 美女高潮的动态| 久久99蜜桃精品久久| 三级国产精品片| 2022亚洲国产成人精品| 久久精品久久久久久久性| 亚洲精品久久久久久婷婷小说 | 成人毛片60女人毛片免费| 1024手机看黄色片| 久久久亚洲精品成人影院| 三级毛片av免费| 男女啪啪激烈高潮av片| 99热全是精品| 国产精品无大码| 嘟嘟电影网在线观看| av免费观看日本| 18禁动态无遮挡网站| 精华霜和精华液先用哪个| 特大巨黑吊av在线直播| 日韩制服骚丝袜av| 日本一本二区三区精品| 亚洲国产欧美人成| 亚洲无线观看免费| 亚洲av熟女| 女人久久www免费人成看片 | 亚洲av中文av极速乱| 51国产日韩欧美| 亚洲精品日韩在线中文字幕| 婷婷色av中文字幕| 亚洲欧洲国产日韩| 十八禁国产超污无遮挡网站| 变态另类丝袜制服| 尾随美女入室| 丝袜美腿在线中文| 全区人妻精品视频| 乱系列少妇在线播放| 男女啪啪激烈高潮av片| 成人鲁丝片一二三区免费| 免费观看a级毛片全部| 国产v大片淫在线免费观看| 在线观看一区二区三区| 精品久久久久久久久久久久久| 国产精品,欧美在线| 久久久a久久爽久久v久久| 亚洲国产欧美在线一区| 久久久久久久久中文| 美女xxoo啪啪120秒动态图| 国产午夜精品论理片| 最近中文字幕高清免费大全6| 精品酒店卫生间| 日韩欧美精品免费久久| 亚洲av日韩在线播放| av女优亚洲男人天堂| 中文字幕av成人在线电影| 三级国产精品片| 少妇人妻一区二区三区视频| av线在线观看网站| 亚洲自偷自拍三级| 97超碰精品成人国产| 观看美女的网站| 99久久九九国产精品国产免费| 18禁在线无遮挡免费观看视频| 夜夜看夜夜爽夜夜摸| 麻豆成人午夜福利视频| 免费看a级黄色片| 18禁动态无遮挡网站| 久久99热这里只有精品18| 久久久a久久爽久久v久久| 亚洲国产欧美在线一区| 午夜视频国产福利| 一个人免费在线观看电影| 欧美激情久久久久久爽电影| 男人的好看免费观看在线视频| 欧美潮喷喷水| 国产成人福利小说| 91精品一卡2卡3卡4卡| 搡老妇女老女人老熟妇| 国产毛片a区久久久久| 六月丁香七月| 国产成人freesex在线| 日本爱情动作片www.在线观看| 波野结衣二区三区在线| 国产精华一区二区三区| 波多野结衣高清无吗| 极品教师在线视频| av免费在线看不卡| 国产激情偷乱视频一区二区| 日本爱情动作片www.在线观看| 国产精品1区2区在线观看.| 2021少妇久久久久久久久久久| 中文字幕av成人在线电影| 免费看av在线观看网站| 九九久久精品国产亚洲av麻豆| 国产成人精品一,二区| 天堂√8在线中文| 精品国内亚洲2022精品成人| 偷拍熟女少妇极品色| 白带黄色成豆腐渣| 99在线视频只有这里精品首页| 国产精品国产三级专区第一集| 女人被狂操c到高潮| 久久久成人免费电影| 国产成人精品一,二区| 综合色av麻豆| 亚洲国产精品合色在线| 亚洲av免费高清在线观看| 波多野结衣高清无吗| 女的被弄到高潮叫床怎么办| 国产一区亚洲一区在线观看| 久久亚洲国产成人精品v| 免费搜索国产男女视频| 成人毛片60女人毛片免费| 国产69精品久久久久777片| 亚洲美女视频黄频| 可以在线观看毛片的网站| 我要搜黄色片| 亚洲国产高清在线一区二区三| 乱系列少妇在线播放| 免费观看a级毛片全部| 亚洲天堂国产精品一区在线| 国产视频首页在线观看| 国产白丝娇喘喷水9色精品| 亚洲av免费高清在线观看| 亚洲真实伦在线观看| 色尼玛亚洲综合影院| АⅤ资源中文在线天堂| 两性午夜刺激爽爽歪歪视频在线观看| 美女cb高潮喷水在线观看| 亚洲精品亚洲一区二区| 成人亚洲精品av一区二区| 亚洲精品色激情综合| 国产亚洲精品av在线| 嫩草影院入口| 亚洲激情五月婷婷啪啪| 九九爱精品视频在线观看| 国产av不卡久久| 国产免费男女视频| 国产亚洲av嫩草精品影院| 国产av在哪里看| 亚洲国产精品国产精品| 成人特级av手机在线观看| 少妇高潮的动态图| 能在线免费观看的黄片| 国产不卡一卡二| 成人综合一区亚洲| 69人妻影院| 亚洲中文字幕一区二区三区有码在线看| a级毛色黄片| 久久久久久久久大av| 小蜜桃在线观看免费完整版高清| 久久韩国三级中文字幕| 精品久久久噜噜| 少妇裸体淫交视频免费看高清| 国产老妇女一区| 午夜激情欧美在线| 亚洲欧美成人综合另类久久久 | 国产黄色小视频在线观看| 九九在线视频观看精品| 亚洲欧美日韩高清专用| 国产免费一级a男人的天堂| 亚洲欧美中文字幕日韩二区| 亚洲欧美精品综合久久99| av.在线天堂| 亚洲国产精品成人久久小说| 国产女主播在线喷水免费视频网站 | 3wmmmm亚洲av在线观看| 99久国产av精品| 在线观看一区二区三区| 九九爱精品视频在线观看| 久久99热这里只有精品18| 尤物成人国产欧美一区二区三区| 亚洲av男天堂| 十八禁国产超污无遮挡网站| 精品午夜福利在线看| 国产探花在线观看一区二区| 九色成人免费人妻av| 18+在线观看网站| 免费大片18禁| 岛国毛片在线播放| 国模一区二区三区四区视频| 亚洲电影在线观看av| 毛片女人毛片| 国产免费福利视频在线观看| 国产黄色视频一区二区在线观看 | 人妻系列 视频| 欧美日韩在线观看h| 亚洲av不卡在线观看| 少妇丰满av| 亚洲av不卡在线观看| 九九久久精品国产亚洲av麻豆| 久久久a久久爽久久v久久| 国产黄a三级三级三级人| 伦理电影大哥的女人| 欧美xxxx性猛交bbbb| 美女国产视频在线观看| 色5月婷婷丁香| 免费观看的影片在线观看| 少妇人妻一区二区三区视频| 中文字幕久久专区| 国产免费视频播放在线视频 | 国产高清国产精品国产三级 | 麻豆精品久久久久久蜜桃| 日韩大片免费观看网站 | 国产视频首页在线观看| 欧美区成人在线视频| 麻豆av噜噜一区二区三区| 黑人高潮一二区| 免费看光身美女| 国产精品一区www在线观看| 草草在线视频免费看| 久久亚洲国产成人精品v| 2021少妇久久久久久久久久久| av国产免费在线观看| 麻豆一二三区av精品| 一边摸一边抽搐一进一小说| 91aial.com中文字幕在线观看| 欧美97在线视频| 欧美成人a在线观看| 2021天堂中文幕一二区在线观| 成人美女网站在线观看视频| 日日摸夜夜添夜夜爱| 男人舔奶头视频| 十八禁国产超污无遮挡网站| 99国产精品一区二区蜜桃av| 日韩av在线免费看完整版不卡| 亚洲天堂国产精品一区在线| 69人妻影院| 日韩av不卡免费在线播放| 老女人水多毛片| 91在线精品国自产拍蜜月| 欧美一区二区精品小视频在线| 天堂影院成人在线观看| 嫩草影院精品99| 国产免费福利视频在线观看| 亚洲不卡免费看| 国产精品女同一区二区软件| 边亲边吃奶的免费视频| 久久久久久久久久久丰满| 狂野欧美激情性xxxx在线观看| 日本免费在线观看一区| 欧美另类亚洲清纯唯美| 欧美性猛交黑人性爽| 亚洲电影在线观看av| 国产精品久久电影中文字幕| 亚洲av电影在线观看一区二区三区 | 欧美一级a爱片免费观看看| 一个人观看的视频www高清免费观看|