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

    海上風(fēng)機(jī)結(jié)構(gòu)振動監(jiān)測試驗與特性分析

    2019-01-14 07:06:44董霄峰練繼建王海軍
    關(guān)鍵詞:停機(jī)葉輪臺風(fēng)

    董霄峰,練繼建,王海軍

    ?

    海上風(fēng)機(jī)結(jié)構(gòu)振動監(jiān)測試驗與特性分析

    董霄峰1, 2,練繼建1,王海軍1

    (1. 天津大學(xué)水利工程仿真與安全國家重點(diǎn)實驗室,天津 300072;2. 天津大學(xué)前沿技術(shù)研究院有限公司,天津 301700)

    以某海上風(fēng)電試驗樣機(jī)為研究對象,基于現(xiàn)場原型觀測獲取整體風(fēng)機(jī)結(jié)構(gòu)在停機(jī)、正常運(yùn)行、開(停)機(jī)及臺風(fēng)工況下的振動響應(yīng)數(shù)據(jù),系統(tǒng)地分析海上風(fēng)機(jī)結(jié)構(gòu)在不同運(yùn)行狀態(tài)下的振動響應(yīng)時頻域特性與變化規(guī)律,探討影響海上風(fēng)機(jī)結(jié)構(gòu)振動的關(guān)鍵因素及其對風(fēng)機(jī)振動安全性的影響規(guī)律.結(jié)果表明:在風(fēng)機(jī)處于停機(jī)狀態(tài)下,結(jié)構(gòu)振動隨環(huán)境風(fēng)速的增加而顯著增大;而運(yùn)行狀態(tài)時,葉輪轉(zhuǎn)速對結(jié)構(gòu)振動影響效應(yīng)明顯,對振動變化起到主導(dǎo)作用.海上風(fēng)機(jī)在開(停)機(jī)過渡過程中振動顯著大于正常運(yùn)行狀態(tài),臺風(fēng)期風(fēng)機(jī)在額定轉(zhuǎn)速運(yùn)行時的振動超過實測最大風(fēng)速時停機(jī)工況的響應(yīng).塔筒頂部與基礎(chǔ)頂部振動響應(yīng)同步性說明整體風(fēng)機(jī)具有很好的變形協(xié)調(diào)性.

    海上風(fēng)電;現(xiàn)場實測;振動特性;運(yùn)行因素;臺風(fēng)工況

    海上風(fēng)電工程具有大風(fēng)速、高產(chǎn)出及長壽命等優(yōu)勢,近年來得到了快速發(fā)展[1].然而,隨著風(fēng)機(jī)單機(jī)裝機(jī)容量逐漸增加,相應(yīng)的配套葉片長度與塔筒高度也隨之增大,例如我國福建平海灣海上風(fēng)電場[2](容量5MW、葉片64.0m、塔筒90.0m)與英國Westermost海上風(fēng)電場[3](容量6MW、葉片77.0m、塔筒102.0m)等.海上風(fēng)機(jī)逐漸呈現(xiàn)大容量、高塔筒、長葉片的發(fā)展趨勢,由于風(fēng)機(jī)結(jié)構(gòu)與葉輪葉片均屬于柔性結(jié)構(gòu),因此,隨這種變化而來的一個重要問題就是風(fēng)機(jī)在運(yùn)行過程中的結(jié)構(gòu)振動安全問題.

    近年來,隨著測試方法與監(jiān)測技術(shù)逐漸成熟,學(xué)術(shù)界基于實測數(shù)據(jù)來獲悉風(fēng)機(jī)結(jié)構(gòu)振動特性方面的研究迅速發(fā)展.Kusiak等[4]通過對兩臺風(fēng)機(jī)實測數(shù)據(jù)進(jìn)行時、頻域分析,以說明結(jié)構(gòu)振動與識別反饋參數(shù)間的關(guān)系.Bassett等[5]針對一臺2.3MW風(fēng)機(jī),實測開機(jī)與穩(wěn)定運(yùn)行狀態(tài)下的加速度數(shù)據(jù),通過小波分解獲得結(jié)構(gòu)振動特征.馬人樂等[6]對3臺風(fēng)機(jī)塔筒進(jìn)行現(xiàn)場測試,提出“槳葉-輪轂-機(jī)艙-塔筒”整體耦合建模方法,指導(dǎo)風(fēng)機(jī)的風(fēng)致動力響應(yīng)分析與振動控制.Yang等[7]考慮數(shù)據(jù)間相關(guān)性,提出一種有效的結(jié)構(gòu)監(jiān)測技術(shù)并實現(xiàn)了風(fēng)機(jī)的安全評估.Schlechtingen等[8-9]基于風(fēng)機(jī)監(jiān)測與數(shù)據(jù)采集系統(tǒng)(SCADA)提出一種新的自適應(yīng)神經(jīng)模糊推理系統(tǒng)模型,以測試不同種類監(jiān)測信號.Kilic等[10]利用無線傳感器建立一套新型的風(fēng)電測試系統(tǒng)SESHMS,并通過現(xiàn)場測試進(jìn)行驗證.Hu等[11-12]對5MW風(fēng)機(jī)結(jié)構(gòu)體系進(jìn)行長達(dá)兩年的現(xiàn)場測試,通過獲得的動態(tài)信號研究結(jié)構(gòu)共振、環(huán)境與運(yùn)行因素對風(fēng)機(jī)動力特性的影響.Currie等[13]將結(jié)構(gòu)整體無線測試技術(shù)引入到原型觀測中,通過位移測試評估風(fēng)機(jī)基礎(chǔ)穩(wěn)定性,并對嵌入式連接環(huán)的變形進(jìn)行預(yù)警.姚興佳[14]通過對某風(fēng)電場風(fēng)機(jī)塔筒進(jìn)行長期監(jiān)測,發(fā)現(xiàn)結(jié)構(gòu)動態(tài)特征與振動隨著機(jī)組葉輪轉(zhuǎn)速的變化而發(fā)生顯著改變.

    綜上所述,目前研究對整體海上風(fēng)機(jī)及支撐結(jié)構(gòu)的振動監(jiān)測較少,持續(xù)的現(xiàn)場觀測數(shù)據(jù)相對匱乏,海上風(fēng)機(jī)結(jié)構(gòu)的基本振動特性還缺乏全面深入的研究.在前期研究中,雖然對塔筒頂部振動特性開展探討并實現(xiàn)風(fēng)機(jī)工作模態(tài)變化規(guī)律的研究[15],但從結(jié)構(gòu)整體安全角度出發(fā),還缺少針對風(fēng)機(jī)底部振動特性及響應(yīng)異同性的研究.因此,本文針對海上風(fēng)機(jī)結(jié)構(gòu)的停機(jī)狀態(tài)、正常運(yùn)行、開停機(jī)及臺風(fēng)等不同工況,開展風(fēng)機(jī)時頻振動特征及影響規(guī)律研究,對比分析影響海上風(fēng)機(jī)頂部與底部振動的關(guān)鍵因素,并探討整體風(fēng)機(jī)結(jié)構(gòu)在不同運(yùn)行狀態(tài)下的振動安全性.本文第1部分介紹研究海上風(fēng)電試驗樣機(jī)的工程概況,第2部分則通過數(shù)值模擬方法掌握海上風(fēng)機(jī)自振特性,第3、4、5部分基于實測響應(yīng)數(shù)據(jù)詳細(xì)分析風(fēng)機(jī)在不同工況下的時、頻振動特性,第6部分為文章結(jié)論.

    1?工程概況

    測試所用海上風(fēng)電試驗樣機(jī)位于我國黃海海域內(nèi),額定功率2.5MW,額定轉(zhuǎn)速18r/min,基礎(chǔ)采用復(fù)合式筒型基礎(chǔ)形式[16-17].風(fēng)機(jī)葉輪直徑93.4m,輪轂高度80.0m,塔筒為3段式安裝,機(jī)艙與塔筒及相鄰塔筒之間設(shè)有工作平臺.本次測試在塔筒內(nèi)部由上至下布置5個測點(diǎn),位置為距離工作平臺高約1.5m處的塔筒壁上,具體布置位置如圖1所示.傳感器選擇低頻振動位移傳感器,頻響范圍最低可達(dá)到0.1Hz,為三向測振,滿足采集結(jié)構(gòu)低頻多向振動信號的試驗?zāi)康模?/p>

    圖1?工程位置與原型觀測測點(diǎn)布置

    2?海上風(fēng)電試驗樣機(jī)自振特性

    為了更好地說明風(fēng)機(jī)結(jié)構(gòu)振動特性,先按風(fēng)機(jī)實際尺寸在有限元軟件ADINA中構(gòu)建數(shù)值模型.基礎(chǔ)與塔筒分別賦予鋼筋混凝土與鋼材屬性,機(jī)艙、輪轂、葉片采用模型等效處理的方式建模,即保證結(jié)構(gòu)整體剛度與質(zhì)量不變的前提下采用較為簡單的體型來替代復(fù)雜的風(fēng)機(jī)機(jī)艙、輪轂、葉輪等結(jié)構(gòu),各材料密度按每部分結(jié)構(gòu)實際質(zhì)量與模型中體積計算.風(fēng)機(jī)葉片材料參考GFRP材料的各項參數(shù),采用正交各向異性材料形式模擬,選定葉片展向模量40.0GPa,剪切模量5.5GPa,泊松比取為0.22,常阻尼比0.02,密度經(jīng)簡化計算為146.08kg/m3.表1列出了數(shù)值模型主要結(jié)構(gòu)的材料參數(shù).

    表1?模型各部分結(jié)構(gòu)的等效材料參數(shù)

    Tab.1?Equivalent material parameters of each model part

    考慮到現(xiàn)場筒型基礎(chǔ)大部分埋于地基中,并經(jīng)過長時間作用與周圍土體形成了一個穩(wěn)定的整體,故為簡化計算將忽略基礎(chǔ)與土質(zhì)地基間的相互作用,同時將模型基礎(chǔ)下部與土體接觸部分的邊界條件按照全約束處理.在整體坐標(biāo)系下建立的風(fēng)機(jī)結(jié)構(gòu)數(shù)值模型如圖2所示.通過模態(tài)分析獲得海上風(fēng)機(jī)結(jié)構(gòu)的前3階模態(tài)頻率與振型如圖3所示.由圖3可知,實測海上風(fēng)機(jī)結(jié)構(gòu)前3階振動頻率分別為0.36Hz、1.46Hz與2.56Hz,與基于實測數(shù)據(jù)的模態(tài)識別結(jié)果基本一致[18],主要體現(xiàn)風(fēng)機(jī)支撐結(jié)構(gòu)的第1、2階擺動振動形式.

    圖2?實際海上風(fēng)機(jī)結(jié)構(gòu)與數(shù)值模型示意

    圖3?海上風(fēng)機(jī)模態(tài)特性

    3?正常工況風(fēng)機(jī)振動特性

    3.1?停機(jī)工況

    選取風(fēng)機(jī)頂部1#測點(diǎn)與底部5#測點(diǎn)作為參考,以實測風(fēng)速在2.0~10.0m/s范圍內(nèi)結(jié)構(gòu)停機(jī)狀態(tài)下20組振動響應(yīng)為研究對象,采樣時長為60s,以水平、向信號的均方根值作為研究標(biāo)準(zhǔn)[19].圖4給出了海上風(fēng)機(jī)振動隨環(huán)境風(fēng)速的變化規(guī)律,由于風(fēng)機(jī)在停機(jī)狀態(tài)工況下僅受到以風(fēng)荷載為主的環(huán)境激勵影響,因此其振幅與風(fēng)速間呈現(xiàn)出明顯的正相關(guān)關(guān)系,且在頻域上主要體現(xiàn)風(fēng)機(jī)1階固有模態(tài)頻率0.324~0.342Hz(如圖5所示).對于1#測點(diǎn)水平向振動,當(dāng)風(fēng)速低于5.0m/s時振動幅度變化較小,低于2.0mm,當(dāng)風(fēng)速超過5.0m/s時向振動幅度迅速增加,并在風(fēng)速為10.0m/s時達(dá)到最大均方根值11.04mm.水平向振動則呈現(xiàn)出逐漸增加的趨勢,最大均方根值為5.95mm.塔筒底部振動與頂部具有明顯的同步性,說明整體風(fēng)機(jī)結(jié)構(gòu)變形協(xié)調(diào)性較好,、向振動最大均方根值分別為0.061mm與0.042mm.

    圖4?停機(jī)工況下水平振動隨風(fēng)速變化趨勢

    圖5?停機(jī)工況下主頻率隨風(fēng)速變化趨勢

    3.2?運(yùn)行工況

    為證明所選運(yùn)行工況可以真實有效地表達(dá)各運(yùn)行因素與結(jié)構(gòu)振動間的關(guān)系,繪制現(xiàn)場實測機(jī)組負(fù)荷與理論機(jī)組負(fù)荷比對如圖6所示,再繪制葉輪轉(zhuǎn)速隨實測風(fēng)速的變化如圖7所示.圖6說明隨風(fēng)速變化現(xiàn)場實測機(jī)組負(fù)荷分布趨勢與理論機(jī)組負(fù)荷分布基本保持一致,且滿足風(fēng)機(jī)機(jī)組理論功率曲線[20].圖7說明葉輪轉(zhuǎn)速隨著運(yùn)行風(fēng)速的增加而逐漸增大,直至風(fēng)速達(dá)到額定風(fēng)速12.0m/s時穩(wěn)定在18r/min運(yùn)行.由此可知,測試數(shù)據(jù)滿足直驅(qū)式風(fēng)力發(fā)電機(jī)的運(yùn)行風(fēng)速與機(jī)組負(fù)荷、葉輪轉(zhuǎn)速間變化關(guān)系,實測數(shù)據(jù)真實可靠.

    圖6?實測負(fù)荷與理論負(fù)荷對比

    圖7?葉輪轉(zhuǎn)速隨風(fēng)速變化趨勢

    圖8給出了兩測點(diǎn)水平、向振動均方根值隨風(fēng)速的變化趨勢.由圖可知,風(fēng)機(jī)在正常運(yùn)行范圍內(nèi),風(fēng)速與水平向振動幅度總體上呈現(xiàn)正相關(guān)性,當(dāng)風(fēng)速達(dá)到風(fēng)機(jī)切入風(fēng)速3.0m/s后,風(fēng)機(jī)開始運(yùn)行并產(chǎn)生明顯振動;隨著風(fēng)速的不斷增大,風(fēng)機(jī)葉輪轉(zhuǎn)速同時不斷增加,結(jié)構(gòu)振動愈發(fā)劇烈;當(dāng)風(fēng)速達(dá)到額定風(fēng)速12.0m/s后,葉輪轉(zhuǎn)速基本保持不變,但由于運(yùn)行風(fēng)速的繼續(xù)增大導(dǎo)致作用在風(fēng)機(jī)上的荷載增加,故結(jié)構(gòu)振動幅度會不斷增加.由圖可知,在運(yùn)行風(fēng)速均超過10m/s時,海上風(fēng)機(jī)1#與5#測點(diǎn)的水平向、向振動最大,其均方根值分別為41.61mm、57.38mm與0.259mm、0.235mm.

    圖8?水平向振動隨風(fēng)速變化趨勢

    由圖9可知,海上風(fēng)機(jī)結(jié)構(gòu)在運(yùn)行狀態(tài)下振動的主要頻域信息具有顯著的變化規(guī)律.當(dāng)風(fēng)機(jī)葉輪轉(zhuǎn)速未超過8r/min時,結(jié)構(gòu)振動主要體現(xiàn)環(huán)境激勵下的1階固有模態(tài)頻率(1#測點(diǎn)為0.319~0.358Hz,5#測點(diǎn)為0.324~0.374Hz);而當(dāng)轉(zhuǎn)速繼續(xù)增加后,風(fēng)機(jī)振動在體現(xiàn)自振特性的同時,也會同時體現(xiàn)與實時葉輪轉(zhuǎn)頻相關(guān)的諧波信息;而當(dāng)轉(zhuǎn)速超過14r/min后,葉輪轉(zhuǎn)動引起諧波激勵,對風(fēng)機(jī)振動起到主導(dǎo)作用,實測振動信號中主要包含的是轉(zhuǎn)頻信息.而相比于1#測點(diǎn),位于風(fēng)機(jī)底部的5#測點(diǎn)同時呈現(xiàn)出其他頻域規(guī)律.一方面是在整個轉(zhuǎn)速范圍內(nèi)其振動主頻率呈現(xiàn)固有模態(tài)的次數(shù)明顯多于1#測點(diǎn);另一方面是兩種荷載聯(lián)合激勵時,5#測點(diǎn)主頻率出現(xiàn)了0.42~0.60Hz的情況.前者主要是由于5#測點(diǎn)距離風(fēng)機(jī)葉輪較遠(yuǎn),上部諧波激勵對其影響程度有限;后者則可能是諧波頻率經(jīng)傳播在風(fēng)機(jī)底部形成疊加而產(chǎn)生了高能量的倍頻振動.結(jié)果可以說明:隨著運(yùn)行風(fēng)速與葉輪轉(zhuǎn)速的增加,誘發(fā)海上風(fēng)機(jī)結(jié)構(gòu)振動的激勵變化遵循由環(huán)境荷載激勵轉(zhuǎn)為環(huán)境荷載激勵和葉輪轉(zhuǎn)動聯(lián)合作用,再到完全由葉輪轉(zhuǎn)動產(chǎn)生諧波激勵影響的規(guī)律.

    圖9?運(yùn)行工況下水平振動主頻率隨風(fēng)速變化趨勢

    4?開停機(jī)工況風(fēng)機(jī)振動特性

    4.1?開機(jī)過程工況

    針對風(fēng)機(jī)開機(jī)過程,分別進(jìn)行了9組不同風(fēng)速工況的監(jiān)測試驗,每次測試時間為150s,并于第60s時啟動機(jī)組開機(jī)運(yùn)行,圖10給出了兩測點(diǎn)在不同風(fēng)速條件下開機(jī)過程中水平向振動位移.

    圖10?不同風(fēng)速下水平z向開機(jī)過程中的振動位移

    由圖10可知,所有測試工況的開機(jī)過程均可分為停機(jī)狀態(tài)、過渡過程、開機(jī)過程及運(yùn)行狀態(tài)4個典型階段.在停機(jī)狀態(tài)時海上風(fēng)機(jī)僅受到環(huán)境荷載作用,由于葉輪處于靜止?fàn)顟B(tài),因此整體結(jié)構(gòu)振動幅度較?。陲L(fēng)機(jī)啟動進(jìn)入過渡過程后,機(jī)組轉(zhuǎn)矩會隨著葉輪的轉(zhuǎn)動而緩慢增加,由于此時機(jī)組并未發(fā)電且轉(zhuǎn)速較低,結(jié)構(gòu)振動幅度相對停機(jī)狀態(tài)會略有增加但依然較?。?dāng)進(jìn)入開機(jī)過程后,葉輪轉(zhuǎn)速繼續(xù)增大,作用在葉輪上的荷載持續(xù)增加,從而外部能量集中轉(zhuǎn)化為結(jié)構(gòu)振動能量而使得振幅會迅速增加.直至達(dá)到穩(wěn)定轉(zhuǎn)速而進(jìn)入穩(wěn)定的運(yùn)行狀態(tài)后,風(fēng)機(jī)電力系統(tǒng)完全接入工作,外部能量將部分轉(zhuǎn)化為電能,結(jié)構(gòu)振動趨于穩(wěn)定.由圖10可以看出,風(fēng)機(jī)在開機(jī)過程中的振動相比運(yùn)行狀態(tài)下均明顯增大,但相比于1#測點(diǎn),位于風(fēng)機(jī)底部的5#測點(diǎn)在開機(jī)過程中振動雖會增大,但突變性并不明顯,這也是由其距離風(fēng)機(jī)制動系統(tǒng)較遠(yuǎn)造成的,結(jié)構(gòu)頂部突變荷載對其影響程度有限.9種工況下兩測點(diǎn)水平方向振動位移均方根值增加百分比見表2,1#與5#測點(diǎn)最大增比分別為42.00%與48.54%,均出現(xiàn)在風(fēng)速為5.5m/s時.

    表2 開機(jī)過程中各工況水平向振動位移均方根值增加百分比

    Tab.2 RMS value increase percentage of horizontal vi-bration displacement under different startup con-ditions

    4.2?停機(jī)過程工況

    針對風(fēng)機(jī)停機(jī)過程,同樣分別進(jìn)行了9組不同風(fēng)速工況的監(jiān)測試驗,測試時間為150s,并選取合理的停機(jī)時間,圖11給出了兩個測點(diǎn)在不同風(fēng)速條件下停機(jī)過程中水平向振動位移.

    由圖11可知,所有測試工況的停機(jī)過程均可分為運(yùn)行狀態(tài)、停機(jī)過程、慣性緩沖過程及停機(jī)狀態(tài)4個典型階段.風(fēng)機(jī)停機(jī)前處于運(yùn)行狀態(tài),結(jié)構(gòu)振動較為平穩(wěn),振動信號呈現(xiàn)出明顯的大幅值、多周期疊加的特征.而風(fēng)機(jī)開始停機(jī)后,由于不再有電能輸出,外部能量會集中轉(zhuǎn)化為動能使得振動幅值在短時間內(nèi)會明顯增加,并因受到機(jī)組剎車制動的影響而快速進(jìn)入慣性緩沖過程,此時風(fēng)機(jī)不受機(jī)組動力影響,而僅受到自身慣性作用逐漸趨于在外部環(huán)境荷載影響下平穩(wěn)響應(yīng)的停機(jī)狀態(tài),停機(jī)后的停機(jī)狀態(tài)下信號具有單一周期、小幅值的振動特性.由圖11可知,風(fēng)機(jī)在停機(jī)過程中振動相比于運(yùn)行狀態(tài)下明顯增大,但5#測點(diǎn)仍由于距離風(fēng)機(jī)制動系統(tǒng)較遠(yuǎn)而振動突變性較弱.9種工況下兩測點(diǎn)水平方向振動位移均方根值增加百分比見表3,1#與5#測點(diǎn)最大增比分別為256.72%與163.02%,均出現(xiàn)在風(fēng)速為7.6m/s時.

    圖11 測點(diǎn)不同風(fēng)速下水平z向停機(jī)過程中的振動位移

    表3 停機(jī)過程中各工況水平向振動位移均方根值增加百分比

    Tab.3 q RMS value increase percentage of horizontal vi-bration displacement under different shutdown conditions

    5?臺風(fēng)工況風(fēng)機(jī)振動特性

    本次試驗觀測到3組臺風(fēng)影響下海上風(fēng)機(jī)振動響應(yīng),各組臺風(fēng)期間實測獲得10min平均風(fēng)速分別為11.0m/s、12.0m/s、13.0m/s,最大瞬時風(fēng)速分別為17.0m/s、26.0m/s、20.0m/s,各測點(diǎn)水平、方向最大振動幅值見表4,向振動最大幅值206.00mm與向最大幅值274.15mm分別出現(xiàn)在臺風(fēng)蘇拉&達(dá)維與布拉萬影響期間.塔筒底部的4#測點(diǎn)和法蘭環(huán)下部靠近基礎(chǔ)的5#測點(diǎn)最大水平位移分別為1.32mm和0.99mm,說明在較強(qiáng)臺風(fēng)作用下風(fēng)機(jī)結(jié)構(gòu)各部位振動滿足安全標(biāo)準(zhǔn).圖12給出3次臺風(fēng)期間1#與5#測點(diǎn)水平、向5min最大絕對振動位移及運(yùn)行風(fēng)速與葉輪轉(zhuǎn)速的變化曲線,羅馬數(shù)字Ⅰ~Ⅴ分別表示單個工況中風(fēng)機(jī)典型的運(yùn)行階段.由圖可見,1#測點(diǎn)與5#測點(diǎn)水平向振動具有明顯的同步性,說明海上風(fēng)機(jī)結(jié)構(gòu)振動的整體性較好,同時進(jìn)一步論證了臺風(fēng)期海上風(fēng)機(jī)結(jié)構(gòu)振動特性與變化規(guī)律[15].

    本次測試臺風(fēng)期運(yùn)行階段主要?dú)w納為正常停機(jī)階段、臺風(fēng)停機(jī)階段、正常運(yùn)行階段與限制運(yùn)行階段.以圖12(a)中的蘇拉&達(dá)維雙臺風(fēng)為例,風(fēng)機(jī)上、下兩個測點(diǎn)振動位移的區(qū)域Ⅰ均體現(xiàn)測試臺風(fēng)過程中正常停機(jī)階段的振動特性,此時由于風(fēng)速未達(dá)到機(jī)組運(yùn)行的切入風(fēng)速,風(fēng)機(jī)無法啟動運(yùn)行而處于轉(zhuǎn)速為零的停機(jī)狀態(tài),結(jié)構(gòu)振動幅度僅與風(fēng)速變化相關(guān),類似狀態(tài)也體現(xiàn)于圖12(b)中的區(qū)域Ⅳ.而圖12(b)的區(qū)域Ⅳ平均風(fēng)速接近7.0m/s,風(fēng)機(jī)仍處于未運(yùn)行狀態(tài),可能是事故檢修造成的正常停機(jī)狀態(tài).

    表4?臺風(fēng)工況下各測點(diǎn)水平向最大振動幅值

    Tab.4 The maximum horizontal vibration amplitude of different measured points under typhoon condi-tion

    注:3#測點(diǎn)、向數(shù)據(jù)超出傳感器量程,故在此不予考慮。

    臺風(fēng)停機(jī)階段體現(xiàn)在如圖12(a)中的區(qū)域Ⅲ,此時實測臺風(fēng)期風(fēng)速較大,但轉(zhuǎn)速與兩測點(diǎn)振動幅度較小,其原因為風(fēng)機(jī)受控制進(jìn)入順槳停機(jī)狀態(tài),結(jié)構(gòu)振動僅受環(huán)境荷載作用,以避免在大風(fēng)速運(yùn)行下產(chǎn)生較大的瞬時荷載導(dǎo)致葉片或結(jié)構(gòu)破壞,類似狀態(tài)也體現(xiàn)于圖12(b)中的區(qū)域Ⅱ與圖12(c)中的區(qū)域Ⅱ.

    圖12(a)中的區(qū)域Ⅱ?qū)儆谂_風(fēng)期的正常運(yùn)行階段,此階段葉輪轉(zhuǎn)速與運(yùn)行風(fēng)速具有明顯的同步性,風(fēng)速范圍往往低于風(fēng)機(jī)額定風(fēng)速,結(jié)構(gòu)振動同時受到葉輪轉(zhuǎn)動產(chǎn)生的諧波激勵與環(huán)境荷載的聯(lián)合影響,類似狀態(tài)也體現(xiàn)于圖12(b)中的區(qū)域Ⅰ與圖12(c)的區(qū)域Ⅲ和Ⅴ.而臺風(fēng)??陂g振動最大值即出現(xiàn)在該階段中圖12(b)中的區(qū)域Ⅰ.

    圖12?臺風(fēng)影響期間風(fēng)機(jī)振動位移時程

    臺風(fēng)期的限制運(yùn)行階段主要體現(xiàn)在圖12(a)的區(qū)域Ⅳ和Ⅴ,此階段運(yùn)行風(fēng)速大多高于風(fēng)機(jī)額定風(fēng)速,而轉(zhuǎn)速與風(fēng)速卻呈現(xiàn)非同步的現(xiàn)象.雖然兩區(qū)域風(fēng)速呈現(xiàn)下降趨勢,但隨著區(qū)域Ⅴ葉輪轉(zhuǎn)速增加,風(fēng)機(jī)振動相比于區(qū)域Ⅳ卻明顯增大,為典型的通過固定機(jī)組負(fù)荷控制葉輪轉(zhuǎn)速的情況,臺風(fēng)達(dá)維&蘇拉期間風(fēng)機(jī)振動的最大值即出現(xiàn)在這樣的階段.類似的現(xiàn)象也出現(xiàn)在圖12(b)的區(qū)域Ⅲ和Ⅴ與圖12(c)的區(qū)域Ⅰ中,前者體現(xiàn)的是在不同固定負(fù)荷條件下運(yùn)行時結(jié)構(gòu)振幅隨葉輪轉(zhuǎn)速的變化而變化.后者則出現(xiàn)了臺風(fēng)布拉萬期間的振動最大值,其開始階段機(jī)組以額定狀態(tài)運(yùn)行,隨后由于機(jī)組負(fù)荷值的限定,葉輪轉(zhuǎn)速依然保持穩(wěn)定的變化而不受風(fēng)速增加的影響,風(fēng)機(jī)振動也隨著轉(zhuǎn)速減小而減弱.

    由以上研究發(fā)現(xiàn),葉輪轉(zhuǎn)速變化對海上風(fēng)機(jī)結(jié)構(gòu)振動的影響要比運(yùn)行風(fēng)速帶來的影響更加明顯,這一點(diǎn)從12(a)的區(qū)域Ⅲ中也可獲得驗證,期間很小的轉(zhuǎn)速增加即可引起結(jié)構(gòu)振動的明顯變化.而在轉(zhuǎn)速穩(wěn)定的情況下,風(fēng)速變化才會對風(fēng)機(jī)振動產(chǎn)生影響,此規(guī)律在12(a)的區(qū)域Ⅳ與12(c)的區(qū)域Ⅰ中有所體現(xiàn),特別是對于位置較低的5#測點(diǎn),振動尤為明顯.

    6?結(jié)?論

    本文以某海上風(fēng)力發(fā)電試驗樣機(jī)為研究對象,以原型觀測的振動位移為研究數(shù)據(jù),全面、深入、系統(tǒng)地分析海上風(fēng)機(jī)整體結(jié)構(gòu)在不同運(yùn)行狀態(tài)下的振動響應(yīng)特性及變化規(guī)律,主要獲得以下結(jié)論.

    (1) 不同工況下海上風(fēng)機(jī)結(jié)構(gòu)振動響應(yīng)特性研究中,風(fēng)機(jī)結(jié)構(gòu)振動幅度均在限定范圍以內(nèi),滿足安全運(yùn)行標(biāo)準(zhǔn).風(fēng)機(jī)塔筒頂部與基礎(chǔ)頂部測點(diǎn)振動響應(yīng)具有很好的同步性,頻域上也呈現(xiàn)出相似的變化特征,說明整體風(fēng)機(jī)結(jié)構(gòu)變形協(xié)調(diào)性較好.

    (2) 海上風(fēng)機(jī)結(jié)構(gòu)在停機(jī)工況下振動僅由環(huán)境荷載決定,水平、方向最大振動均方根值分別為11.04mm與5.95mm,在頻域上體現(xiàn)風(fēng)機(jī)自身模態(tài)屬性,1階模態(tài)頻率為0.324~0.342Hz.

    (3) 在正常運(yùn)行工況下,運(yùn)行風(fēng)速、葉輪轉(zhuǎn)速與機(jī)組負(fù)荷間存在明顯同步性,隨著運(yùn)行因素的增加,風(fēng)機(jī)振動呈現(xiàn)明顯增大趨勢,水平、方向最大振動位移均方根值分別為41.61mm與57.38mm.頻域信息則因外部荷載的改變體現(xiàn)從固有模態(tài)頻率向諧波激勵頻率的過渡過程,底部結(jié)構(gòu)由于遠(yuǎn)離發(fā)電機(jī)而受到更多環(huán)境荷載與諧波倍頻的影響.

    (4) 在整個開、停機(jī)工況下,停機(jī)狀態(tài)與運(yùn)行狀態(tài)均存在一個過渡階段與開(停)機(jī)過程階段,且開(停)機(jī)過程階段相比于運(yùn)行狀態(tài)的結(jié)構(gòu)振動幅度多數(shù)呈現(xiàn)增大趨勢,水平向的振動均方根值分別最大增加42.00%與256.72%,而底部測點(diǎn)振動突變性由于激勵能量傳遞損耗而相對較弱.

    (5) 在臺風(fēng)達(dá)維&蘇拉與布拉萬期間,風(fēng)機(jī)結(jié)構(gòu)振動的最大值出現(xiàn)在額定工作狀態(tài)下,最大振幅分別為206.00mm與274.15mm,而臺風(fēng)??r下振動最大值出現(xiàn)在風(fēng)機(jī)無干擾正常運(yùn)行階段,最大振幅為125.80mm,此結(jié)果可以說明在運(yùn)行狀態(tài)下風(fēng)機(jī)結(jié)構(gòu)振動主要受到葉輪轉(zhuǎn)速的作用,而在轉(zhuǎn)速穩(wěn)定的情況下,風(fēng)速變化才會對風(fēng)機(jī)振動產(chǎn)生影響,對基礎(chǔ)底部測點(diǎn)尤為明顯.

    [1] Ahmed N A,Cameron M. The challenges and possible solutions of horizontal axis wind turbines as a clean energy solution for the future[J]. Renewable & Sustainable Energy Reviews,2014,38(5):439-460.

    [2] Rodrigues S,Restrepo C,Kontos E,et al. Trends of offshore wind projects[J]. Renewable & Sustainable Energy Reviews,2015,49:1114-1135.

    [3] Colmenar-Santos A,Perera-Perez J,Borge-Diez D,?et al. Offshore wind energy:A review of the current status,challenges and future development in Spain [J]. Renewable & Sustainable Energy Reviews,2016,64:1-18.

    [4] Kusiak A,Zhang Zijun. Analysis of wind turbine vibrations based on SCADA data[J]. Journal of Solar Energy Engineering,2010,132(3):1-12.

    [5] Basset K,Rupp C,Ting D S K. Vibration analysis of 2.3 MW wind turbine operation using the discrete wavelet transform[J]. Wind Engineering,2010,34(4):375-388.

    [6] 馬人樂,馬躍強(qiáng),劉慧群. 風(fēng)電機(jī)組塔筒模態(tài)的環(huán)境脈動實測與數(shù)值模擬研究[J]. 振動與沖擊,2011,30(5):152-155.

    Ma Renle,Ma Yueqiang,Liu Huiqun. Ambient vibration test and numerical simulation for modes of wind turbine towers[J]. Journal of Vibration and Shock,2011,30(5):152-155(in Chinese).

    [7] Yang W X,Court R,Jiang J K. Wind turbine condition monitoring by the approach of SCADA data analysis[J]. Renewable Energy,2013,53(9):365-376.

    [8] Schlechtingen M,Santos I F,Achiche S. Wind turbine condition monitoring based on SCADA data using normal behavior models. Part 1:System description[J]. Applied Soft Computing Journal,2013,13(1):259-270.

    [9] Schlechtingen M,Santos I F,Achiche S. Wind turbine condition monitoring based on SCADA data using normal behavior models. Part 2:Application examples[J]. Applied Soft Computing Journal,2014,14(1):447-460.

    [10] Kilic G,Unluturk M S. Testing of wind turbine towers using wireless sensor network and accelerometer[J]. Renewable Energy,2015,75:318-325.

    [11] Hu W H,Th?ns S,Rohrmann R G,et al. Vibration-based structural health monitoring of a wind turbine system. Part I:Resonance phenomenon[J]. Engineering Structures,2015,89:260-272.

    [12] Hu W H,Th?ns S,Rohrmann R G,et al. Vibration-based structural health monitoring of a wind turbine system Part II:Environmental/operational effects on dynamic properties[J]. Engineering Structures,2015,89:273-290.

    [13] Currie M,Saafi M,Tachtatzis C,et al. Structural integrity monitoring of onshore wind turbine concrete foundations[J]. Renewable Energy,2015,83(4):1131-1138.

    [14] 姚興佳. 風(fēng)力發(fā)電測試技術(shù)[M]. 北京:電子工業(yè)出版社,2011.

    Yao Xingjia. Measurement Technology in Wind Power [M]. Beijing:Electronic Industry Press,2011(in Chinese).

    [15] Dong X,Lian J,Wang H,et al. Structural vibration monitoring and operational modal analysis of offshore wind turbine structure[J]. Ocean Engineering,2018,150:280-297.

    [16] Lian Jijian,Ding Hongyan,Zhang Puyang. Design of large-scale prestressing bucket foundation for offshore wind turbines[J]. Transactions of Tianjin University,2012,18(2):79-84.

    [17] Lian Jijian,Sun Liqiang,Zhang Jinfeng,et al. Bearing capacity and technical advantages of composite bucket foundation of offshore wind turbines[J]. Transactions of Tianjin University,2011,17(2):132-137.

    [18] 董霄峰,練繼建,楊?敏,等. 海上風(fēng)機(jī)結(jié)構(gòu)工作模態(tài)識別的組合降噪方法[J]. 天津大學(xué)學(xué)報:自然科學(xué)與工程技術(shù)版,2015,48(3):203-208.

    Dong Xiaofeng,Lian Jijian Yang Min,et al. A compound de-noising method for operational modal identification of offshore wind turbine structure[J]. Journal of Tianjin University:Science and Technology,2015,48(3):203-208(in Chinese).

    [19] 練繼建,王海軍,秦?亮. 水電站廠房結(jié)構(gòu)研究[M]. 北京:中國水利水電出版社,2007.

    Lian Jijian,Wang Haijun,Qin Liang. Research on Hydropower House Structure[M]. Beijing:China Water & Power Press,2007(in Chinese).

    [20] 董霄峰,練繼建,王海軍. 運(yùn)行狀態(tài)下海上風(fēng)機(jī)結(jié)構(gòu)振源特性研究[J]. 振動與沖擊,2017,36(17):21-28.

    Dong Xiaofeng,Lian Jijian,Wang Haijun. Study on vibration source features of offshore wind power structure under the operational conditions[J]. Journal of Vibration and Shock,2017,36(17):21-28(in Chinese).

    (責(zé)任編輯:樊素英)

    Monitoring Experiment and Characteristic Analysis of Structural Vibration of Offshore Wind Turbine

    Dong Xiaofeng1, 2,Lian Jijian1,Wang Haijun1

    (1. State Key Laboratory of Hydraulic Engineering Simulation and Safety,Tianjin University,Tianjin 300072,China;2. Frontier Technology Research Institute Limited Company,Tianjin University,Tianjin 301700,China)

    This study considered one offshore wind turbine prototype as the research object,and its structural vibration responses under standstill,normal operational,startup,shutdown,and typhoon conditions were obtained through field measurement. First,the variation characteristics of the vibration responses under different operational conditions are examined systematically through time and frequency domain analyses. Furthermore,the key factors influencing structural vibration and their influence rules on the vibration safety of offshore wind turbine are explained clearly. Results show that the structural vibration increases with the increase in wind speed when the wind turbine operates under the standstill condition. The research results also show that the rotation speed of the blades has an obvious influence on structural vibration and plays a dominant role during vibration under the normal operational condition. Thereafter,the structural vibration of offshore wind turbine under the transition process of the startup(shutdown) condition is significantly greater than that under the normal operational condition with the same variables. Simultaneously,the observed vibration under the typhoon condition with rated rotation speed exceeds the corresponding responses under the standstill condition with maximum measured wind speed. Finally,the synchronization of the vibration responses measured from the tower and foundation indicates that the overall structure has good deformation coordination.

    offshore wind power;field measurement;vibration characteristic;operational factor;typhoon condition

    the Innovation Method Fund of China(No. 2016IM030100),the National Natural Science Foundation of China (No. 51709202),the Program of Introducing Talents of Discipline to Universities(No. B14012),Tianjin Science and Technology Program(No. 16PTGCCX00160).

    10.11784/tdxbz201708007

    TK8

    A

    0493-2137(2019)02-0191-09

    2017-08-05;

    2018-08-28.

    董霄峰(1986— ),男,博士,講師.

    董霄峰,tju_dongxf@126.com.

    國家科技計劃項目創(chuàng)新方法工作專項資助(2016IM030100);國家自然科學(xué)基金資助項目(51709202);高等學(xué)校學(xué)科創(chuàng)新引智計劃資助項目(B14012);天津市科技計劃資助項目(16PTGCCX00160).

    猜你喜歡
    停機(jī)葉輪臺風(fēng)
    臺風(fēng)過韓
    質(zhì)量管理工具在減少CT停機(jī)天數(shù)中的應(yīng)用
    1.4317 QT2鋼在高能泵葉輪上的應(yīng)用
    臺風(fēng)來了
    小讀者(2020年4期)2020-06-16 03:33:46
    臺風(fēng)愛搗亂
    應(yīng)用石膏型快速精密鑄造技術(shù)制造葉輪
    臺風(fēng)來時怎樣應(yīng)對
    離心泵葉輪切割方法
    雷克薩斯NX200t車停機(jī)和起動系統(tǒng)解析
    基于CFD/CSD耦合的葉輪機(jī)葉片失速顫振計算
    日本av免费视频播放| 午夜老司机福利剧场| 搡女人真爽免费视频火全软件| 亚洲经典国产精华液单| 99热这里只有是精品50| 特大巨黑吊av在线直播| 汤姆久久久久久久影院中文字幕| 国产精品蜜桃在线观看| 免费av中文字幕在线| 亚洲精品日韩在线中文字幕| 午夜影院在线不卡| 三级国产精品片| 人妻少妇偷人精品九色| 午夜福利影视在线免费观看| 一二三四中文在线观看免费高清| 免费av中文字幕在线| 黑人猛操日本美女一级片| 3wmmmm亚洲av在线观看| 精品国产露脸久久av麻豆| 日日摸夜夜添夜夜添av毛片| 久久这里有精品视频免费| 免费人成在线观看视频色| 精品亚洲成国产av| 免费看av在线观看网站| 天堂俺去俺来也www色官网| 久久久久网色| av国产久精品久网站免费入址| 久久久精品94久久精品| 美女国产视频在线观看| 亚洲精品日韩av片在线观看| 一级,二级,三级黄色视频| 黑丝袜美女国产一区| 久久久欧美国产精品| 如何舔出高潮| 国产淫片久久久久久久久| 人人妻人人添人人爽欧美一区卜| 色吧在线观看| 欧美日韩在线观看h| 一级毛片黄色毛片免费观看视频| 国产又色又爽无遮挡免| 久久久a久久爽久久v久久| 狠狠精品人妻久久久久久综合| 精品人妻熟女av久视频| 亚洲成人手机| 观看美女的网站| 97在线视频观看| 少妇人妻 视频| 69精品国产乱码久久久| 深夜a级毛片| 男女啪啪激烈高潮av片| 久久99一区二区三区| 三级经典国产精品| 99国产精品免费福利视频| 亚洲精华国产精华液的使用体验| 成年女人在线观看亚洲视频| 新久久久久国产一级毛片| 99九九在线精品视频 | 亚州av有码| 久久久久久久久大av| 国产极品粉嫩免费观看在线 | 97超碰精品成人国产| 97超碰精品成人国产| 在线精品无人区一区二区三| 男人添女人高潮全过程视频| 黄色欧美视频在线观看| 十八禁高潮呻吟视频 | 国产黄色免费在线视频| 美女大奶头黄色视频| 免费高清在线观看视频在线观看| 亚洲av二区三区四区| 亚洲,欧美,日韩| 99国产精品免费福利视频| 欧美+日韩+精品| 一个人看视频在线观看www免费| 中国美白少妇内射xxxbb| 国产在线男女| 中国美白少妇内射xxxbb| 国产探花极品一区二区| 麻豆成人av视频| 大香蕉久久网| 大香蕉久久网| 国产成人精品福利久久| 老司机影院毛片| 男女免费视频国产| 午夜福利在线观看免费完整高清在| 久久久国产精品麻豆| 日韩成人伦理影院| 国产在线男女| 亚洲,欧美,日韩| 欧美+日韩+精品| 中国美白少妇内射xxxbb| 黄色一级大片看看| 麻豆成人av视频| 99re6热这里在线精品视频| 国产欧美亚洲国产| 亚洲成人av在线免费| 看免费成人av毛片| 美女脱内裤让男人舔精品视频| 久久人妻熟女aⅴ| 亚洲激情五月婷婷啪啪| freevideosex欧美| 国产一区二区在线观看av| 99久久中文字幕三级久久日本| av又黄又爽大尺度在线免费看| 国产精品人妻久久久影院| 国产成人91sexporn| 高清av免费在线| 久久午夜综合久久蜜桃| 欧美性感艳星| 亚洲色图综合在线观看| 午夜福利影视在线免费观看| 高清欧美精品videossex| 久久精品国产亚洲网站| 久久精品久久久久久久性| 中文字幕亚洲精品专区| 极品少妇高潮喷水抽搐| 中文字幕久久专区| 蜜桃在线观看..| 欧美人与善性xxx| 岛国毛片在线播放| 制服丝袜香蕉在线| freevideosex欧美| 国产一区二区在线观看日韩| 精品国产一区二区久久| 国产伦精品一区二区三区视频9| 桃花免费在线播放| 97在线视频观看| 视频中文字幕在线观看| 久久午夜福利片| 久久久久网色| 嘟嘟电影网在线观看| 99久久人妻综合| 免费人成在线观看视频色| 中国美白少妇内射xxxbb| 久久精品国产亚洲av涩爱| a级片在线免费高清观看视频| 精品亚洲成国产av| 婷婷色av中文字幕| 热re99久久精品国产66热6| 久久精品国产a三级三级三级| 秋霞伦理黄片| 亚洲性久久影院| 亚洲av二区三区四区| 亚洲欧美一区二区三区黑人 | 交换朋友夫妻互换小说| 免费高清在线观看视频在线观看| 国产伦精品一区二区三区视频9| 免费久久久久久久精品成人欧美视频 | 国产成人精品无人区| 大话2 男鬼变身卡| 亚洲真实伦在线观看| 菩萨蛮人人尽说江南好唐韦庄| av在线老鸭窝| 日韩一区二区视频免费看| 伦精品一区二区三区| av天堂中文字幕网| 国产女主播在线喷水免费视频网站| 亚洲精品日韩在线中文字幕| 26uuu在线亚洲综合色| 五月玫瑰六月丁香| 国产亚洲最大av| 热re99久久国产66热| 99久国产av精品国产电影| 亚洲成人一二三区av| 亚洲国产精品国产精品| 两个人免费观看高清视频 | 亚洲欧洲精品一区二区精品久久久 | 亚洲成人手机| 最近手机中文字幕大全| 日本黄色片子视频| 涩涩av久久男人的天堂| 午夜免费观看性视频| 国产成人精品福利久久| 中文天堂在线官网| 美女主播在线视频| 亚洲国产欧美在线一区| 精品午夜福利在线看| 日本爱情动作片www.在线观看| 婷婷色综合www| 青春草国产在线视频| 高清毛片免费看| 熟女av电影| 天堂中文最新版在线下载| 精华霜和精华液先用哪个| 亚洲国产欧美在线一区| 久久99热这里只频精品6学生| 我要看黄色一级片免费的| 超碰97精品在线观看| 成人美女网站在线观看视频| 九草在线视频观看| 欧美精品一区二区大全| 搡女人真爽免费视频火全软件| 国产精品一区二区在线不卡| 人妻人人澡人人爽人人| 男人狂女人下面高潮的视频| 夜夜爽夜夜爽视频| 午夜福利,免费看| 日本-黄色视频高清免费观看| 日本色播在线视频| 国国产精品蜜臀av免费| 成人无遮挡网站| 亚洲精品成人av观看孕妇| 日本91视频免费播放| 高清视频免费观看一区二区| 国产成人精品久久久久久| 午夜av观看不卡| 十分钟在线观看高清视频www | 啦啦啦视频在线资源免费观看| 热99国产精品久久久久久7| 中文字幕亚洲精品专区| 国产av一区二区精品久久| 黄色配什么色好看| 中文天堂在线官网| 日韩精品有码人妻一区| 成人18禁高潮啪啪吃奶动态图 | www.色视频.com| 青春草亚洲视频在线观看| 国产精品嫩草影院av在线观看| 嘟嘟电影网在线观看| 国产日韩欧美亚洲二区| 国产精品国产av在线观看| 亚洲国产精品一区二区三区在线| 婷婷色综合大香蕉| 97在线视频观看| 国产精品久久久久久av不卡| a级毛色黄片| 国产日韩欧美在线精品| 国产精品.久久久| 久久99蜜桃精品久久| 又粗又硬又长又爽又黄的视频| 日本av手机在线免费观看| 97精品久久久久久久久久精品| 精品亚洲乱码少妇综合久久| 久久av网站| 91午夜精品亚洲一区二区三区| 一本—道久久a久久精品蜜桃钙片| 在线观看美女被高潮喷水网站| 亚洲四区av| 伊人久久国产一区二区| 少妇高潮的动态图| 免费高清在线观看视频在线观看| 午夜久久久在线观看| 在线观看美女被高潮喷水网站| 国产 精品1| 国产乱人偷精品视频| 亚洲国产日韩一区二区| 亚洲精品久久久久久婷婷小说| 在线观看三级黄色| 国产在视频线精品| 亚洲av在线观看美女高潮| 精品人妻熟女毛片av久久网站| 男人舔奶头视频| av免费在线看不卡| 91精品国产国语对白视频| 特大巨黑吊av在线直播| 亚洲国产日韩一区二区| 亚洲一区二区三区欧美精品| 丰满人妻一区二区三区视频av| 国精品久久久久久国模美| 免费在线观看成人毛片| 成年女人在线观看亚洲视频| 夜夜看夜夜爽夜夜摸| 黑人高潮一二区| 午夜激情久久久久久久| 国产免费视频播放在线视频| 亚洲精品456在线播放app| 一个人免费看片子| 黄色日韩在线| 日本av免费视频播放| 国产综合精华液| 久久久久久久大尺度免费视频| 尾随美女入室| 国产黄色免费在线视频| 精品午夜福利在线看| 王馨瑶露胸无遮挡在线观看| 国产精品久久久久久久久免| 国产在线男女| 91午夜精品亚洲一区二区三区| 国产成人免费观看mmmm| 在线观看国产h片| 少妇熟女欧美另类| 蜜桃在线观看..| 在线天堂最新版资源| 国内揄拍国产精品人妻在线| av天堂久久9| 在线免费观看不下载黄p国产| 少妇被粗大的猛进出69影院 | 国产视频首页在线观看| 国产 精品1| 亚洲丝袜综合中文字幕| 亚洲精品aⅴ在线观看| 精品午夜福利在线看| 99久久精品国产国产毛片| 精品亚洲成a人片在线观看| 人体艺术视频欧美日本| 国产亚洲91精品色在线| 国产精品一区www在线观看| 国产欧美日韩精品一区二区| 久久久久久久久久久久大奶| 夫妻午夜视频| 国产日韩欧美视频二区| 欧美+日韩+精品| 精品久久久久久久久亚洲| 久久97久久精品| 成年av动漫网址| 亚洲经典国产精华液单| 少妇人妻 视频| 欧美日韩国产mv在线观看视频| 国产成人精品一,二区| 国产精品一区二区性色av| 久久久久人妻精品一区果冻| 久久99热6这里只有精品| 亚洲丝袜综合中文字幕| 久久免费观看电影| 久久午夜综合久久蜜桃| 男人狂女人下面高潮的视频| 伦理电影大哥的女人| 成人国产麻豆网| 婷婷色综合www| 免费人成在线观看视频色| 中文天堂在线官网| 日本猛色少妇xxxxx猛交久久| 日韩 亚洲 欧美在线| 日韩在线高清观看一区二区三区| 国产黄色免费在线视频| 国产黄片视频在线免费观看| av线在线观看网站| 欧美三级亚洲精品| 欧美另类一区| 欧美3d第一页| 看十八女毛片水多多多| 亚洲国产av新网站| 亚洲婷婷狠狠爱综合网| 丝瓜视频免费看黄片| 99精国产麻豆久久婷婷| 亚洲欧洲日产国产| 少妇 在线观看| 午夜福利,免费看| 九九久久精品国产亚洲av麻豆| 一本久久精品| av专区在线播放| 国产黄频视频在线观看| 91久久精品国产一区二区成人| tube8黄色片| 欧美精品人与动牲交sv欧美| 国产精品女同一区二区软件| 乱人伦中国视频| 欧美精品一区二区大全| 成人漫画全彩无遮挡| 中文字幕制服av| 国产av码专区亚洲av| 免费av中文字幕在线| 久热久热在线精品观看| 夜夜看夜夜爽夜夜摸| 国产亚洲最大av| 香蕉精品网在线| 在线天堂最新版资源| 亚洲精品国产色婷婷电影| 亚洲av综合色区一区| 亚洲国产日韩一区二区| 亚洲国产欧美在线一区| 久久99热6这里只有精品| 少妇高潮的动态图| 爱豆传媒免费全集在线观看| 少妇 在线观看| 女性生殖器流出的白浆| 国产精品人妻久久久影院| 国产中年淑女户外野战色| 色网站视频免费| 成人18禁高潮啪啪吃奶动态图 | 国产免费一区二区三区四区乱码| 日韩视频在线欧美| 欧美日韩在线观看h| 日韩免费高清中文字幕av| 97超碰精品成人国产| 搡女人真爽免费视频火全软件| 五月天丁香电影| 久久久久久久久久成人| av免费在线看不卡| 国产精品国产三级国产专区5o| 黄色一级大片看看| 久久久久视频综合| 国产成人aa在线观看| 黄色一级大片看看| 免费久久久久久久精品成人欧美视频 | 午夜福利影视在线免费观看| 赤兔流量卡办理| 九九久久精品国产亚洲av麻豆| 免费黄色在线免费观看| www.av在线官网国产| 国产亚洲91精品色在线| 美女视频免费永久观看网站| 中文天堂在线官网| 91久久精品国产一区二区成人| 国产av一区二区精品久久| 亚洲av成人精品一二三区| 少妇 在线观看| 国产一级毛片在线| 最近中文字幕2019免费版| 成人18禁高潮啪啪吃奶动态图 | 亚洲欧美日韩另类电影网站| 久久女婷五月综合色啪小说| a级毛片免费高清观看在线播放| 国产欧美日韩一区二区三区在线 | 国产精品成人在线| 在现免费观看毛片| 内射极品少妇av片p| 一本—道久久a久久精品蜜桃钙片| 久久99热6这里只有精品| 午夜福利在线观看免费完整高清在| 精品少妇久久久久久888优播| 日本欧美视频一区| 欧美日韩精品成人综合77777| 熟女人妻精品中文字幕| 色94色欧美一区二区| 精品久久久久久久久av| 乱人伦中国视频| 国产午夜精品久久久久久一区二区三区| 丰满乱子伦码专区| 国产乱人偷精品视频| 22中文网久久字幕| 六月丁香七月| 欧美老熟妇乱子伦牲交| 日韩制服骚丝袜av| 99热6这里只有精品| 欧美日韩综合久久久久久| 一级毛片久久久久久久久女| 久久精品熟女亚洲av麻豆精品| 精品久久久久久久久亚洲| 亚洲欧美中文字幕日韩二区| 国产黄片美女视频| 在线观看国产h片| av在线老鸭窝| 少妇被粗大猛烈的视频| 2021少妇久久久久久久久久久| 久久亚洲国产成人精品v| 久久久久视频综合| 纯流量卡能插随身wifi吗| 日韩av不卡免费在线播放| 国语对白做爰xxxⅹ性视频网站| 久久av网站| 欧美精品国产亚洲| 久久精品国产a三级三级三级| 一区二区三区精品91| 亚洲精品日本国产第一区| 亚洲成人av在线免费| 黄色怎么调成土黄色| 亚洲人成网站在线观看播放| 久久久久久伊人网av| 国产精品蜜桃在线观看| 日韩一本色道免费dvd| 99久久综合免费| 亚洲久久久国产精品| 色哟哟·www| 一级爰片在线观看| 日韩不卡一区二区三区视频在线| 亚洲精品成人av观看孕妇| 乱人伦中国视频| 一本—道久久a久久精品蜜桃钙片| 国产日韩欧美视频二区| 日本-黄色视频高清免费观看| av卡一久久| 国产成人一区二区在线| 日韩精品有码人妻一区| 精品久久久精品久久久| 老女人水多毛片| 老司机影院毛片| 成人毛片60女人毛片免费| 日本av免费视频播放| 一本色道久久久久久精品综合| 亚洲av电影在线观看一区二区三区| 亚洲高清免费不卡视频| 欧美精品亚洲一区二区| 黄色配什么色好看| 欧美性感艳星| 男人添女人高潮全过程视频| 国产欧美日韩综合在线一区二区 | 久热这里只有精品99| 18禁裸乳无遮挡动漫免费视频| 18禁动态无遮挡网站| 亚洲精品日韩在线中文字幕| 国产淫片久久久久久久久| 国产黄片美女视频| 日本爱情动作片www.在线观看| 久热这里只有精品99| 国产精品国产三级国产av玫瑰| 成年av动漫网址| 日本色播在线视频| 久久精品国产亚洲网站| 男的添女的下面高潮视频| 日本欧美国产在线视频| 人人妻人人澡人人爽人人夜夜| 在线观看美女被高潮喷水网站| av视频免费观看在线观看| 国产亚洲91精品色在线| 久久久久久久久久久丰满| 一区在线观看完整版| 久久精品久久久久久久性| 国内少妇人妻偷人精品xxx网站| 色婷婷久久久亚洲欧美| 国产有黄有色有爽视频| 好男人视频免费观看在线| 国产精品.久久久| 久久99热这里只频精品6学生| 中国美白少妇内射xxxbb| 大码成人一级视频| 久久久久久久久久久丰满| 日韩中字成人| 久久久久国产精品人妻一区二区| 大陆偷拍与自拍| 男男h啪啪无遮挡| 成年人免费黄色播放视频 | 午夜激情福利司机影院| 亚洲精品aⅴ在线观看| 国模一区二区三区四区视频| 女人久久www免费人成看片| 人人妻人人添人人爽欧美一区卜| 国产精品不卡视频一区二区| 激情五月婷婷亚洲| 国产精品不卡视频一区二区| 黄色欧美视频在线观看| 五月玫瑰六月丁香| 国产日韩欧美亚洲二区| 日韩精品有码人妻一区| 中文字幕亚洲精品专区| 久久国产乱子免费精品| 建设人人有责人人尽责人人享有的| 亚洲av欧美aⅴ国产| 国产精品嫩草影院av在线观看| 欧美日韩一区二区视频在线观看视频在线| 18禁在线无遮挡免费观看视频| 亚洲欧美日韩东京热| 久久综合国产亚洲精品| 91精品一卡2卡3卡4卡| 国产高清国产精品国产三级| 美女中出高潮动态图| 日韩三级伦理在线观看| 日日啪夜夜爽| 一区二区av电影网| 亚洲精品国产av蜜桃| 看非洲黑人一级黄片| 亚洲国产成人一精品久久久| 69精品国产乱码久久久| 啦啦啦啦在线视频资源| 自拍偷自拍亚洲精品老妇| 国产乱来视频区| 亚洲精品自拍成人| 国产视频内射| 日本午夜av视频| 少妇熟女欧美另类| freevideosex欧美| 人妻系列 视频| 国产黄色免费在线视频| a 毛片基地| 久久久久久久久久久久大奶| 下体分泌物呈黄色| av网站免费在线观看视频| 简卡轻食公司| 国产女主播在线喷水免费视频网站| 亚洲欧美中文字幕日韩二区| 国产女主播在线喷水免费视频网站| 亚洲国产精品成人久久小说| 成人毛片a级毛片在线播放| 欧美日韩av久久| 亚洲精品视频女| 性色av一级| 成人黄色视频免费在线看| 香蕉精品网在线| 免费播放大片免费观看视频在线观看| 免费看av在线观看网站| 午夜av观看不卡| 高清毛片免费看| 我的女老师完整版在线观看| 日本欧美视频一区| 国产国拍精品亚洲av在线观看| 国产在线男女| 美女xxoo啪啪120秒动态图| 久久亚洲国产成人精品v| 日韩成人av中文字幕在线观看| 一个人看视频在线观看www免费| 久久97久久精品| 亚洲经典国产精华液单| 嫩草影院入口| 九九在线视频观看精品| 亚洲av免费高清在线观看| 天天躁夜夜躁狠狠久久av| 午夜激情福利司机影院| 日本av免费视频播放| 熟妇人妻不卡中文字幕| 99久久中文字幕三级久久日本| 夜夜爽夜夜爽视频| 我要看日韩黄色一级片| 精品人妻熟女毛片av久久网站| 国产免费视频播放在线视频| 美女脱内裤让男人舔精品视频| 女性被躁到高潮视频| 日韩av不卡免费在线播放| 婷婷色av中文字幕| 搡女人真爽免费视频火全软件| 乱码一卡2卡4卡精品| 国产成人精品婷婷| 亚洲国产精品999| 少妇猛男粗大的猛烈进出视频| 日本免费在线观看一区| 国产淫片久久久久久久久| 久久99精品国语久久久| 欧美日韩亚洲高清精品| 日本vs欧美在线观看视频 | 丰满人妻一区二区三区视频av| 成年女人在线观看亚洲视频| av天堂久久9| 五月开心婷婷网| 国产精品久久久久成人av| 久久免费观看电影| 男人狂女人下面高潮的视频|