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

    湍流強度對風(fēng)電機組動力學(xué)特性及載荷的影響

    2020-08-12 15:01:52巫發(fā)明楊從新
    農(nóng)業(yè)工程學(xué)報 2020年13期
    關(guān)鍵詞:葉根湍流彎矩

    巫發(fā)明,楊從新,王 靛,楊 柳

    湍流強度對風(fēng)電機組動力學(xué)特性及載荷的影響

    巫發(fā)明1,2,楊從新1※,王 靛2,楊 柳2

    (1. 蘭州理工大學(xué)能源與動力工程學(xué)院,蘭州 730050;2. 中車株洲電力機車研究所有限公司風(fēng)電事業(yè)部,株洲 412000)

    為研究湍流強度對機組動力學(xué)特性及氣動載荷的影響,以3.3 MW三葉片水平軸風(fēng)電機組為研究對象,采用仿真及試驗相結(jié)合的方法,并對來流風(fēng)速和主導(dǎo)載荷進行功率譜分析。通過開展4種湍流強度0.10、0.12、0.14和0.16的計算仿真,結(jié)果表明隨著湍流強度的增加,風(fēng)電機組機艙振動加速度、載荷及等效疲勞載荷都有規(guī)律性變化。為驗證仿真結(jié)果的合理性,對某風(fēng)場的型式測試機組進行1a多的數(shù)據(jù)測試采集和分析。測試結(jié)果表明,機組運行在0.06、0.08、0.10和0.12這4種不同湍流強度下,其機組在不同風(fēng)速運行下的機組振動及載荷同樣出現(xiàn)有規(guī)律性的變化,仿真與實測結(jié)果的變化趨勢吻合度較高。該研究為風(fēng)電場風(fēng)電機組的微觀選址提供依據(jù),也對風(fēng)電機組設(shè)計有一定的指導(dǎo)意義。

    風(fēng);發(fā)電;動力學(xué);風(fēng)電機組;湍流強度;氣動載荷

    0 引 言

    湍流是普遍存在且被許多研究者持續(xù)探索的自然現(xiàn)象。大氣湍流對飛行器的飛行性能、結(jié)構(gòu)載荷、飛行安全的影響很大。因湍流引發(fā)的飛行事故時有發(fā)生,但通過現(xiàn)代技術(shù)可以有效避開強湍流或盡量降低危害程度。但風(fēng)電機組固定安裝于某個位置,當(dāng)天氣變化,導(dǎo)致機組可能受到不同程度湍流影響,為保證風(fēng)電機組的安全性和長期穩(wěn)定可靠運行,機組的設(shè)計需要考慮運行環(huán)境條件的湍流變化帶來的載荷問題[1-2]。在過去幾年里,越來越明顯的是風(fēng)電機組對風(fēng)能的轉(zhuǎn)化必須面對更高湍流的挑戰(zhàn),湍流強度對大型風(fēng)電機組的載荷的影響需要探討。在不同湍流強度下,其風(fēng)電機組載荷有相應(yīng)的變化[3]。

    近年來,風(fēng)電機組變得越來越大,陸地大功率機組直徑已達到160 m,海上機組直徑接近200 m,從流體力學(xué)的角度來看,風(fēng)電機組是在湍流大氣邊界層中運行的大型機械,這是非常大的挑戰(zhàn),因為風(fēng)場的多尺度特性,從剖面上的局部流動條件到風(fēng)電場尾流的相互作用[4-8]。復(fù)雜地形會使大氣邊界層流動產(chǎn)生擾動,湍流強度加局,從而影響風(fēng)電機組的效率,破壞機組的安全運行,縮短機組的壽命。因為湍流對風(fēng)電機組載荷、壽命和動力學(xué)特性有直接影響[9],當(dāng)湍流強度大時,增加風(fēng)電機組的疲勞載荷,加劇機組的振動,最終造成機組損壞。強烈的湍流擾動可能是導(dǎo)致機組斷裂損壞的主要原因[10-11]。

    湍流強度的增加使風(fēng)電機組的氣動力、重力和慣性載荷之間的耦合效應(yīng)更加復(fù)雜,國內(nèi)外一些研究工作者通過建立風(fēng)電機組模型了解風(fēng)電機組的動態(tài)總體性能、載荷情況,不少學(xué)者對風(fēng)電機組模擬的方法進行了研究,并成功地應(yīng)用于工程實踐。也有不少學(xué)者對如何減小風(fēng)電機組的載荷,提高風(fēng)電機組性能做了相關(guān)的研究,但針對湍流強度如何影響風(fēng)電機組載荷和動力學(xué)特性的研究還很少見。國際電工委員會發(fā)布的最新版IEC61400-1標(biāo)準[12],給出湍流強度指標(biāo)對風(fēng)電機組的載荷設(shè)計進行了分類,其湍流強度等級由第三版的3個等級增加到4個等級,正與國內(nèi)復(fù)雜地形相對應(yīng),特別容易出現(xiàn)個別機位的湍流強度特別大,出現(xiàn)超A類的現(xiàn)象。國際上Vera-Tudela等[13-14]通過仿真手段指出,當(dāng)湍流強度較大時,風(fēng)電機組輸出功率會相應(yīng)減少,風(fēng)電機組的疲勞載荷會增大,同時,還可能引起較大的瞬態(tài)載荷,最終造成風(fēng)電機組疲勞破壞,但是沒有實際測試數(shù)據(jù)支撐。葉杭冶等[15-18]對如何控制策略的方法,減小風(fēng)電機組的疲勞載荷,提高風(fēng)電機組性能及壽命進行了相關(guān)的仿真研究,但未見有關(guān)于湍流強度對風(fēng)電機組載荷隨風(fēng)速變化和應(yīng)用實測與仿真對比相結(jié)合的研究。Tony等[10]研究指出通常將作用在旋轉(zhuǎn)葉片上的穩(wěn)定風(fēng)速與波動風(fēng)速產(chǎn)生的載荷分開。將風(fēng)速在葉輪掃掠面積空間內(nèi)的變化,所產(chǎn)生的葉片上的周期性載荷稱為確定載荷部分,因為它是由一些參數(shù)唯一確定,如輪轂高度、轉(zhuǎn)速、風(fēng)剪切力等。另外,由風(fēng)速波動產(chǎn)生(即湍流)的隨機載荷只能通過概率預(yù)測來描述,因此稱其為隨機載荷部分。除了風(fēng)載,旋轉(zhuǎn)葉片還受重力和慣性載荷的作用。重力載荷只取決于葉片的角度和質(zhì)量分布,但慣性載荷可能受湍流的影響。從國內(nèi)外學(xué)者研究結(jié)果來看,有必要對湍流強度與風(fēng)電機組載荷、動力學(xué)特性之間的內(nèi)在聯(lián)系開展相應(yīng)的研究工作。本文以3.3MW變速變槳型風(fēng)電機組為例,主要研究湍流強度對大型風(fēng)電機組載荷、動力學(xué)特性的影響,并進行了仿真計算和測試驗證。

    1 研究方法及計算模型

    湍流強度反映了風(fēng)速的波動。湍流強度越大,氣流越不穩(wěn)定,風(fēng)速波動越大,使得風(fēng)電機組的氣動載荷變成了交變載荷。由于風(fēng)電機組承受的重力載荷、風(fēng)載、以及慣性載荷,最終都通過傳動鏈傳遞至塔筒。因此,對于風(fēng)電機組而言,葉根和塔底承受的彎矩載荷是最大的,值得重點關(guān)注的。依據(jù)文獻[18-19],塔底、葉根及傳動鏈的載荷都是考慮俯仰彎矩、傾覆彎矩及扭矩。另外,風(fēng)電機組是包含剛性體和柔性體的多體系統(tǒng),包含葉片、機艙和塔架等部件。各部件之間的相對運動以及各部件彈性變形相互影響和耦合,形成了風(fēng)電機組系統(tǒng)復(fù)雜的結(jié)構(gòu)動力學(xué)特性。湍流強度影響將直接帶來機組動力學(xué)特性問題。因此,本文以葉根,傳動鏈及塔底俯仰彎矩為例,對湍流強度與載荷之間的內(nèi)在聯(lián)系進行研究。同時,分析湍流強度變化對機艙振動加速度的影響。

    1.1 湍流強度

    湍流強度與離地高度,地表面粗糙度有關(guān)。然而,它也受地貌特征的影響,如高地和山脈,位于上風(fēng)向的樹和建筑物等。對于風(fēng)力機運行所面臨的大氣邊界層環(huán)境條件,主要功率譜有IEC Kaimal譜模型:

    中性大氣條件下,3個方向(=)的速度譜為

    式中hub為風(fēng)電機組的輪轂高度,m。

    式中ref為10 min平均風(fēng)速為15m/s時輪轂高度處湍流強度的期望值[12]。

    1.2 動力學(xué)模型

    對于個自由度具有阻尼的線性系統(tǒng)運動方程式為

    風(fēng)電機組機艙振動特性類似于阻尼強迫振動的特點[10],風(fēng)電機組動力學(xué)模型如下

    1.3 湍流的影響及功率譜

    另外,可采用功率頻譜描述在不同的頻率下能量的波動的分布。隨機性載荷能很方便地在頻域范圍內(nèi)進行分析,但是為了更加便于分析,相干函數(shù)通常采用指數(shù)相干模型,相干函數(shù)(Δ,)表達式如下:

    脈動風(fēng)速功率譜表現(xiàn)了脈動風(fēng)能量在整個頻率范圍內(nèi)的分布特征,是進行結(jié)構(gòu)風(fēng)荷載效應(yīng)隨機振動分析的前提之一。功率譜密度存在如下性質(zhì):輸出信號的功率譜密度等于輸入信號的功率譜密度與傳遞函數(shù)模值平方的乘積[20]。

    1.4 等效疲勞載荷

    等效疲勞載荷是對給定測量時間內(nèi)疲勞載荷效果的一種方便、簡明的描述。等效載荷是一個概括的信號載荷,在給定的時間歷程中以給定的頻率出現(xiàn)的總循環(huán)次數(shù)施加這個載荷時,它所造成的疲勞損傷與測得的載荷譜中所有雨流計算出的載荷幅值造成的損傷之各是相同的。

    雨流循環(huán)計算是用作結(jié)構(gòu)疲勞分析公認的方法。指定一個或多個材料S-N曲線的斜率與頻率,計算出的等效載荷作為恒頻正弦載荷的振幅,它將得到與原始信號一樣的疲勞損傷。載荷測試的10 min結(jié)果的1 Hz等效載荷由式(10)給出。

    2 仿真模型分析與驗證

    湍流對風(fēng)電機組壽命和可靠性有直接影響,影響風(fēng)電機組壽命的主要是機組承受的疲勞載荷和振動。因此本文主要分析了不同湍流度下機組的振動特性和疲勞載荷變化特性。

    2.1 計算模型

    為驗證湍流強度與機組動力學(xué)特性、載荷的相關(guān)性,本文運用風(fēng)電行業(yè)專業(yè)軟件GH-Bladed,對3.3MW風(fēng)電機組進行仿真計算。該機組的主要參數(shù)如表1所示。

    表1 3.3MW風(fēng)電機組參數(shù)

    計算條件為:平均風(fēng)速3~20 m/s,湍流強度ref為0.16(),0.14(),0.12()和0.10()。風(fēng)電機組振動特性及載荷分析中,需建立坐標(biāo)系才能對整機及部件進行建模和計算分析。國際電氣委員會(IEC)制定的相應(yīng)標(biāo)準中規(guī)定了坐標(biāo)系[12]。

    2.2 不同風(fēng)速下的機艙振動加速度對比分析

    為分析湍流強度變化對風(fēng)電機組動力學(xué)特性的影響。對正常運行風(fēng)速3~20 m/s下,不同湍流強度的風(fēng)電機組10 min內(nèi)最大值,最小值及平均值機艙振動加速度情況進行統(tǒng)計。

    從圖1a可知,機艙前后振動與湍流強度關(guān)聯(lián)較大,湍流強度上升0.02,其振動加速度值上升達到10%。從圖1b可知,其中機艙左右振動與湍流強度關(guān)聯(lián)也較大,且在額定風(fēng)速附近出現(xiàn)峰值點,原因是此處有轉(zhuǎn)矩控制到變槳控制切換的過程,存在一定的沖擊。由式(9)可知,湍流強度上升0.02,氣動力的波動變化可達到14%。再根據(jù)式(6)可知,機艙振動加速度與氣動力波動強度相關(guān)。

    圖1 不同湍流強度下機艙振動加速度隨風(fēng)速變化曲線(仿真)

    2.3 不同風(fēng)速下的運行載荷對比分析

    為分析湍流強度變化對風(fēng)電機組載荷的影響。對正常運行風(fēng)速3~20 m/s下,不同湍流強度的風(fēng)電機組10 min內(nèi)最大值,最小值及平均值彎矩載荷情況進行統(tǒng)計,同樣分析葉片、傳動鏈和塔筒的受載情況。葉片主要分析不同風(fēng)速下的葉根擺振彎矩(Mx)和葉根揮舞彎矩載荷(My);傳動鏈在不同風(fēng)速下的旋轉(zhuǎn)輪轂轉(zhuǎn)矩(Mx1)、俯仰彎矩(My1)和偏航彎矩(Mz1)的載荷;塔筒在不同風(fēng)速下塔頂左右彎矩(Mx2)、前后彎矩(My2)、扭矩(Mz2)和塔底左右彎矩(Mx3)、前后彎矩(My3)載荷。

    從圖2a和圖2b可知,葉根Mx與湍流強度關(guān)聯(lián)不大,但葉根揮舞彎矩My受湍流強度影響較為明顯,湍流強度上升0.02,其載荷上升達到10%。其Mx影響較小的原因是葉根擺振彎矩Mx源于葉片自身重力的影響,這是文獻[10]指出的確定載荷部分,與湍流強度關(guān)聯(lián)性不大。

    圖2 葉根彎矩隨風(fēng)速變化曲線(仿真)

    從圖3可知,旋轉(zhuǎn)輪轂Mx1與湍流強度關(guān)聯(lián)不大,但My1、Mz1受湍流強度影響較為明顯,湍流強度上升0.02,其載荷上升達到10%。其中旋轉(zhuǎn)輪轂Mx載荷影響較小的原因,轉(zhuǎn)矩是受轉(zhuǎn)速控制,能保持穩(wěn)定。

    從圖4可知,塔頂Mx2與湍流強度關(guān)聯(lián)較小,塔頂My2和扭矩Mz2受湍流強度影響較為明顯,湍流強度上升0.02,其載荷上升達到10%。

    從圖5可知,塔底Mx3與湍流強度關(guān)聯(lián)性隨著風(fēng)速的增加而增大,塔底My3受湍流強度影響較為明顯,湍流強度上升0.02,其載荷上升達到10%,特別是隨風(fēng)速的增大后更加明顯。

    2.4 功率譜對比分析

    針對風(fēng)速12m/s,湍流強度ref分別0.12和0.16。對風(fēng)速、葉片、輪轂和塔頂Mx、My載荷進行功率譜分析,具體如圖6,7所示。

    圖3 旋轉(zhuǎn)輪轂彎矩載荷隨風(fēng)速變化曲線(仿真)

    圖4 塔頂彎矩隨風(fēng)速變化曲線(仿真)

    從圖6、7的功率譜可明顯看出,風(fēng)速沒有明顯峰值,但葉片、輪轂和塔頂Mx,My主導(dǎo)載荷中出現(xiàn)明顯的尖峰,這主要來源于風(fēng)電機組轉(zhuǎn)頻和部件固有頻率,其中風(fēng)電機組轉(zhuǎn)頻1P約為0.175 Hz,塔筒1階約0.254 Hz,葉片揮舞1階約0.53 Hz,葉片擺振1階約0.92 Hz,傳動鏈1階或葉片揮舞2階約1.24 Hz,葉片擺振2階約2.5 Hz。其峰值都隨著湍流強度變大,主要原因是其波動能量隨著湍流強度變大而增加,但部件固有頻率值不發(fā)生改變。從功率譜變化情況說明湍流強度對風(fēng)電機組有明顯的影響,特別是葉片、輪轂和塔頂My載荷。

    圖5 塔底彎矩隨風(fēng)速變化曲線(仿真)

    2.5 等效疲勞載荷對比分析

    根據(jù)IEC61400-13標(biāo)準[21],在不同湍流強度下計算風(fēng)電機組的載荷,并進行等效疲勞載荷統(tǒng)計,表2給出葉根、傳動鏈和塔筒各關(guān)鍵載荷隨湍流強度變化的比值。隨著湍流強度上升,等效載荷有明顯的上升。湍流強度每上升0.02,傳動鏈(旋轉(zhuǎn)輪轂)和塔筒(塔頂和塔底)等效疲勞載荷基本有10%的上升。但葉根等效載荷Mx變化較小,其主要原因是葉根Mx載荷波動來源葉片質(zhì)量,與湍流強度相關(guān)性較小。

    圖7 風(fēng)速及主導(dǎo)載荷的功率譜(Iref-0.16)

    表2 不同湍流風(fēng)電機組等效疲勞載荷比值

    注:、、、為湍流強度0.16、0.14、0.12、0.10。

    Note:,,,are turbulence intensities 0.16, 0.14, 0.12, 0.10, respectively.

    3 試驗分析及驗證

    3.1 測試條件

    為了進一步研究和分析湍流強度引起的載荷問題,本文按照IEC61400-13的標(biāo)準要求[21],在青海都蘭風(fēng)電場的測試樣機上進行長達1年多的數(shù)據(jù)采集,并進行分析研究,測試樣機如圖8測試風(fēng)場的空氣密度0.9 kg/m3,年平均風(fēng)速7.5 m/s。測試機構(gòu)為國內(nèi)認證公司CGC。

    圖8 樣機

    為分析驗證湍流對正常運行下風(fēng)電機組的影響,主要通過測試機艙振動加速度,葉片和塔筒的受載情況來進行分析。本文測試分別在葉根、塔頂和塔底的安裝了載荷測試的應(yīng)變片傳感器,對載荷進行長期監(jiān)測。其測試的樣品量達到752個數(shù)據(jù)點,測試湍流區(qū)間主在0.06~0.12范圍內(nèi)653個數(shù)據(jù)點。圖9分別給出葉根Mx和My的測試散點圖,包括載荷的最大值,最小值,平均值和標(biāo)準偏差值。從圖9可以看出葉根Mx載荷集中性比較高,葉根My離散度相對高一些。

    圖9 葉根擺振和揮舞彎矩統(tǒng)計數(shù)據(jù)(測試)

    3.2 測試結(jié)果分析

    在不同風(fēng)速下,機艙振動前后加速度和左右加速度、葉根Mx和My載荷、塔筒頂部Mx2,My2,Mz2和塔筒底部Mx3,My3載荷。具體如圖10-圖13。

    圖10 機艙振動加速度隨風(fēng)速變化曲線(測試)

    從圖10a可知,機艙前后振動與湍流強度關(guān)聯(lián)較大,湍流強度上升0.02,其載荷上升達到10%。從圖10b可知,其中機艙左右振動與湍流強度關(guān)聯(lián)較大,且在額定風(fēng)速附近出現(xiàn)峰值點,仿真結(jié)果與試驗結(jié)果趨勢一致,其原因是此處有轉(zhuǎn)矩控制到變槳控制切換的過程,存在一定的沖擊。

    圖11 葉根彎矩隨風(fēng)速變化曲線(測試)

    從圖11測試結(jié)果可知,葉根載荷中Mx與湍流強度關(guān)聯(lián)不大,與仿真趨勢基本致。同樣My受湍流強度影響較為明顯,湍流強度上升0.02,其載荷上升達到10%。

    從圖12測試結(jié)果可知,塔頂My與Mz載荷與湍流強度相關(guān)性都較大,湍流強度上升0.02,其載荷都上升達到10%。但塔頂Mx的載荷在達到額定風(fēng)速后,其最大值基本不受湍流的影響,但最小值還隨湍流變化有一定的變化,其主要原因是達到額定后,最大轉(zhuǎn)矩已受到限定。與仿真結(jié)果趨一致。

    從圖13測試結(jié)果可知,塔底Mx、My載荷與湍流強度相關(guān)性較大,湍流強度上升0.02,其載荷都上升達到10%。

    從圖10-圖13測試結(jié)果來看,總體趨勢仿真基本一致,說明仿真的合理性,但由于測試的樣品量在高速和高湍流的數(shù)據(jù)量偏少,還存在一些奇異點。但總趨勢還是能反映出合理性。因此,風(fēng)電機組的優(yōu)化設(shè)需著重考慮風(fēng)場湍流強度影響。湍流強度是需要提前確定并在風(fēng)電機組載荷設(shè)計及風(fēng)電場選址過程中予以考慮的,對于湍流強度較大的地區(qū),風(fēng)電機組設(shè)計載荷應(yīng)考慮以適當(dāng)?shù)南禂?shù)進行修正,以保證風(fēng)電機組的疲勞壽命安全。

    圖12 塔頂彎矩隨風(fēng)速變化曲線(測試)

    圖13 塔底彎矩隨風(fēng)速變化曲線(測試)

    4 結(jié) 論

    本文研究風(fēng)電場湍流強度對風(fēng)電機組的影響,從分析情況可知,湍流強度的變化影響風(fēng)電機組的振動和疲勞載荷,并通過樣機測試進行驗證。主要結(jié)論如下:

    1)從仿真結(jié)果來看,隨著湍流強度的增加,風(fēng)電機組機艙振動加速度增加。湍流強度每上升0.02,機艙前后振動加速度值的增幅也達到10%。其中機艙左右振動與湍流強度關(guān)聯(lián)也較大,且在額定風(fēng)速附近出現(xiàn)峰值點,原因是此處有轉(zhuǎn)矩控制到變槳控制切換的過程,存在一定的沖擊。

    2)從仿真結(jié)果來看,隨著湍流強度的增加,風(fēng)電機組的重要部件的主導(dǎo)載荷和疲勞載荷增加。湍流強度每上升0.02,其葉根My、輪轂My1、塔頂和塔底My2、My3的等效疲勞載荷接近10%的上升;從功率譜分析來看,湍流強度對風(fēng)電機組有明顯的影響,特別是葉片、輪轂和塔頂My載荷。

    3)從測試結(jié)果來看,湍流強度與風(fēng)電機組機艙振動和疲勞等效載荷強相關(guān)。測試結(jié)果表明,機組在不同風(fēng)速運行下的機組振動及載荷同樣出現(xiàn)有規(guī)律性的變化,仿真與實測結(jié)果的變化趨勢吻合度較高。

    由以上分析可知,對于湍流強度較大的地區(qū),風(fēng)電機組設(shè)計載荷應(yīng)考慮以適當(dāng)?shù)南禂?shù)進行修正,以保證風(fēng)電機組的疲勞壽命安全。因此,對湍流高的地區(qū)的風(fēng)電機組部件的承受能力要適當(dāng)加強。這為風(fēng)電場風(fēng)電機組的微觀選址提供依據(jù),也對風(fēng)電機組設(shè)計有一定的指導(dǎo)意義。

    [1] 李鴻秀,朱瑞兆,王蕊,等. 不同地形風(fēng)電場湍流強度日變化和年變化分析[J]. 太陽能學(xué)報,2014,35(11):2327-2333.

    Li Hongxiu, Zhu Ruizhao, Wang Rui. Wind turbulence intensity of daily and yearly variation analysis in different kinds of terrain[J]. Acta Energiae Solaris Sinica, 2014, 35(11): 2327-2333. (in Chinese with English abstract)

    [2] 付德義,薛揚,焦渤,等. 湍流強度對風(fēng)電機組疲勞等效載荷的影響[J]. 華北電力大學(xué)學(xué)報,2015,42(1):45-50.

    Fu Deyi, Xue Yang, Jiao Bo, et al. Effect on the turbulence intensity to wind turbine fatigue equivalent load[J]. Journal of North China Electric Power University, 2015, 42(1): 45-50. (in Chinese with English abstract)

    [3] 鄧英,周峰,田德,等. 不同風(fēng)湍流模型的風(fēng)電機組載荷計算研究[J]. 太陽能學(xué)報,2014,35(12):2395-2400.

    Deng Ying, Zhou Feng, Tian De, et al. Research on load calculation of wind turbine for different wind turbulence model[J]. Acta Energiae Solaris Sinica, 2014, 35(12): 2395-2400. (in Chinese with English abstract)

    [4] 曹娜,于群,王偉勝,等. 風(fēng)電場尾流效應(yīng)模型研究[J]. 太陽能學(xué)報,2016,37(1):222-229.

    Cao Na, Yu Qun, Wang Weisheng, et al. Research on wakeeffect model on wind farm[J]. Acta Energiae Solaris Sinica, 2016, 37(1): 222-229. (in Chinese with English abstract)

    [5] 李少華,岳巍澎,匡青峰,等. 雙機組風(fēng)力機尾流互擾及陣列的數(shù)值模擬[J]. 中國電機工程學(xué)報,2011,31(5):101-107.

    Li Shaohua, Yue Weipeng, Kuang Qingfeng, et al. Numerical simulation of wake interaction and array of double wind turbine[J]. Proceedings of the CSEE, 2011, 31(5): 101-107. (in Chinese with English abstract)

    [6] 宋麗莉,吳戰(zhàn)平,秦鵬. 復(fù)雜山地近地層強風(fēng)特性分析[J]. 氣象學(xué)報,2009,67(3):452-460

    Song Lili, Wu Zhanping, Qin Peng. An analysis of the characteristics of strong winds in the surface layer over a complex terrain. Acta Meteorologica Sinica , 2009, 67(3): 452-460. (in Chinese with English abstract)

    [7] 吳豐林,方創(chuàng)琳. 中國風(fēng)能資源價值評估與開發(fā)階段劃分研究[J]. 自然資源學(xué)報,2009,24(8):1412-1421.

    Wu Fenglin, Fang Chuanglin. Wind power resource appraisa land development stage regional division of China[J]. Journal of Natural Resources, 2009, 24(8): 1412-1421. (in Chinese with English abstract)

    [8] 何鴻川,袁巍. 典型山地地形風(fēng)力機氣動性能評估及布局探索[J]. 太陽能學(xué)報,2017,38(10):2683-2692.

    He Hongchuan, Yuan Wei. Aerodynamic performance evaluation and layout exploration of wind turbine in typical mountainous terrain[J]. Acta Energiae Solaris Sinica, 2017, 38(10): 2683-2692. (in Chinese with English abstract)

    [9] 鄧英,李嘉楠,劉河生,等. 風(fēng)電機組尾流與疲勞載荷關(guān)系分析[J]. 農(nóng)業(yè)工程學(xué)報,2017,33(10):239-244.

    Deng Ying, Li Jianan, Liu Hesheng, et al. Analysis on relationship between wake and fatigue load of wind turbines[J]. Transactions of the Chinese Society of Agricultural Engineering, 2017, 33(10): 239-244. (in Chinese with English abstract)

    [10] Tony Burton, David Sharpe, Nick Jenkins. Wind Energy Handbook[M]. New York: A John Wiley and Sons, Ltd., Publication England, 2011

    [11] 廖明夫,宋文萍,王四季. 風(fēng)力機設(shè)計理論與結(jié)構(gòu)動力學(xué)[M]. 西安:西北工業(yè)大學(xué)出版社,2014

    [12] International Electro technical Commission. Wind turbines generator systems-part 1: Safety requirements, IEC61400-1.[S]. Geneva: International Electro technical Commission, 2019.

    [13] Vera-Tudela L, Kuhn M. Analysing wind turbine fatigue load prediction: The impact of wind farm flow conditions[J]. Renewable Energy, 2017, 107: 352-360.

    [14] Siyu Tao, Qingshan Xu. Integrated wind farm power curve and power curve distribution unction considering the wake Effect and Terrain Gradient[J]. Energies, 2019, 12: 2482.

    [15] 何偉. 湍流風(fēng)場模擬與風(fēng)力發(fā)電機組載荷特性研究[D]. 保定:華北電力大學(xué),2013.

    He Wei. Research on Turbulent Wind Field Simulation and Load Characteristics of Wind Turbines[D]. Baoding: North China Electric Power University, 2013. (in Chinese with English abstract)

    [16] 李晶,宋家驊,王偉勝. 大型變速恒頻風(fēng)力發(fā)電機組建模與仿真[J]. 中國電機工程學(xué)報,2004,24(6):100-105.

    Li Jing, Song Jiahua, Wang Weisheng. Modeling and dynamic simulation of variable speed wind turbine with large capacity[J]. Proceedings of the CSEE, 2004, 24(6): 100-105. (in Chinese with English abstract)

    [17] 林勇剛,李偉,葉杭冶,等. 變速恒頻風(fēng)力機組變槳距控制系統(tǒng)[J]. 農(nóng)業(yè)機械學(xué)報,2004,35(4):110-114.

    Lin Yonggang, Li Wei, Ye Hangye, et al. Pitchcontrolled system of a variable-speed and constant-frequency wind turbine[J]. Transactions of The Chinese Society for Agricultural Machinery, 2004, 35(4): 110-114. (in Chinese with English abstract)

    [18] 葉杭冶. 風(fēng)力發(fā)電機組的控制技術(shù)[M]. 北京:機械工業(yè)出版社,2002

    [19] 賀德馨. 風(fēng)工程與工業(yè)空氣動力學(xué)[M]. 北京:國防工業(yè)出版社,2006

    [20] 李尚遠,甘德強,欒某德. 一種基于功率譜密度的發(fā)電機機械功率擾動源定位方法[J]. 電網(wǎng)技術(shù),2019,43(1):236-243.

    Li Shangyuan, Gan Deqiang, Luan Moude. A novel method for locating forced oscillation source in mechanical power based on power spectral density[J]. Power System Technology, 2019, 43(1): 236-243. (in Chinese with English abstract)

    [21] International Electro technical Commission. Wind turbines generator systems-part 13: Conformity testing and certification, IEC61400-13[S]. Geneva: International Electro technical Commission, 2015.

    Effects of turbulence intensity on dynamic characteristics and load of wind turbine

    Wu Faming1,2, Yang Congxin1※,Wang Dian2, Yang Liu2

    (1.,,730050,;2..,.,412000,)

    In order to study the influence of turbulence intensity on the dynamic characteristics and aerodynamic load, simulation tests and experiments are carried out for a horizontal axis wind turbine with three blades whose rated power is 3.3MW. Furthermore, the power spectrum analysis of incoming wind speed and dominant load is carried out. Simulation tests are taken under the condition of four turbulence intensities of 0.10, 0.12, 0.14 and 0.16, whose results show that with the increase of turbulence intensity, the vibration acceleration, load and equivalent fatigue load of wind turbine nacelle change regularly. In order to verify the rationality of the simulation results, the data acquisition and analysis of a type-test wind turbine for more than one year were carried out. The main conclusions are as follows: 1) From the simulation results, with the increase of turbulence intensity, the vibration acceleration of wind turbine nacelle increases. When the turbulence intensity increases by 0.02, the acceleration increases by 10%. Besides, there is also a strong correlation between the side-to-side vibration and turbulence intensity where there is a peak point near the rated wind speed. The reason is that there is a switching process from torque control to pitch control, and there is a certain impact. 2) From the simulation results, with the increase of turbulence intensity, the key load and fatigue load of important components of wind turbine increase. When the turbulence intensity increases by 0.02, the equivalent fatigue loads of blade root MY, hub MY, tower top and tower bottom MY increase by nearly 10%. From the result of power spectrum analysis, the turbulence intensity has obvious influence on wind turbine, especially on MY load of blade, hub and tower top. 3) From the test results, the turbulence intensity is strongly related to the vibration and fatigue equivalent load of the wind turbine nacelle, and vibration and load of the unit also change regularly under the four different turbulence intensities of 0.06, 0.08, 0.10 and 0.12, which is in good agreement with the simulation results. From the above analysis, it can be seen that for areas with high turbulence intensity, the design load of wind turbine should be modified with appropriate coefficient in order to ensure the safety of fatigue life of wind turbine. Therefore, the bearing capacity of wind turbine components in high turbulence area should be strengthened. This provides a basis for the micro sitting of wind turbines in wind farms, and has a certain guiding significance for the design of wind turbines.

    wind;power generation; dynamics; wind turbine; turbulence; aerodynamic load

    巫發(fā)明,楊從新,王靛,等. 湍流強度對風(fēng)電機組動力學(xué)特性及載荷的影響[J]. 農(nóng)業(yè)工程學(xué)報,2020,36(13):48-55.doi:10.11975/j.issn.1002-6819.2020.13.006 http://www.tcsae.org

    Wu Faming, Yang Congxin, Wang Dian, et al. Effects of turbulence intensity on dynamic characteristics and load of wind turbine[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(13): 48-55. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2020.13.006 http://www.tcsae.org

    2020-01-06

    2020-05-10

    政府間國際科技創(chuàng)新合作重點專項(2017YFE0101900);國家自然科學(xué)基金:(11262011)

    巫發(fā)明,博士,高級工程師,主要從事空氣動力學(xué)、風(fēng)電機組設(shè)計與研究。Email:1wufaming@163.com

    楊從新,教授,博士導(dǎo)師,主要從事空氣動力學(xué)、流體機械的設(shè)計與研究。Email:ycxwind@163.com

    10.11975/j.issn.1002-6819.2020.13.006

    TK83

    A

    1002-6819(2020)-13-0048-08

    猜你喜歡
    葉根湍流彎矩
    戒賭迷局
    基于有限元模型仿真的風(fēng)電葉根T型螺母應(yīng)力計算方法研究
    機械工程師(2020年6期)2020-07-14 09:49:22
    三齒樅樹型葉根輪槽型線優(yōu)化設(shè)計
    熱力透平(2019年4期)2019-12-18 07:09:18
    零彎矩設(shè)計理論在連續(xù)梁橋中的應(yīng)用研究
    重氣瞬時泄漏擴散的湍流模型驗證
    精銑葉根的葉片測頻問題分析與對策
    東方汽輪機(2017年4期)2018-01-05 03:06:14
    CFRP-PCPs復(fù)合筋連續(xù)梁開裂截面彎矩計算方法研究
    鋼-混疊合連續(xù)梁負彎矩區(qū)計算分析
    板孔式有彎矩平衡梁應(yīng)用技術(shù)及研究
    “青春期”湍流中的智慧引渡(三)
    亚洲av成人av| 成在线人永久免费视频| 久9热在线精品视频| 校园春色视频在线观看| 欧美日韩瑟瑟在线播放| 一级a爱片免费观看的视频| 国产三级黄色录像| 中文字幕人妻熟女乱码| 国产精品电影一区二区三区 | 久久久国产精品麻豆| 久久性视频一级片| 免费久久久久久久精品成人欧美视频| av免费在线观看网站| 啦啦啦免费观看视频1| 极品教师在线免费播放| 国产区一区二久久| 不卡av一区二区三区| 国产男靠女视频免费网站| 中文字幕人妻丝袜一区二区| 亚洲熟女精品中文字幕| 91精品国产国语对白视频| 一进一出抽搐gif免费好疼 | 欧美精品高潮呻吟av久久| 国产一区二区三区视频了| av天堂久久9| 制服人妻中文乱码| 成人影院久久| www.熟女人妻精品国产| 黄色怎么调成土黄色| 欧美最黄视频在线播放免费 | bbb黄色大片| 黄色怎么调成土黄色| 妹子高潮喷水视频| 99热只有精品国产| 亚洲熟女精品中文字幕| 亚洲欧美日韩高清在线视频| 妹子高潮喷水视频| 久久久水蜜桃国产精品网| 男人的好看免费观看在线视频 | 亚洲精品在线美女| 亚洲熟妇熟女久久| 很黄的视频免费| 黄色毛片三级朝国网站| 男女午夜视频在线观看| 久久精品国产亚洲av高清一级| 中文字幕色久视频| 精品国内亚洲2022精品成人 | 成人免费观看视频高清| 日本黄色日本黄色录像| 高清毛片免费观看视频网站 | 久久精品人人爽人人爽视色| 亚洲欧美日韩另类电影网站| 岛国在线观看网站| 一级毛片女人18水好多| 国产日韩一区二区三区精品不卡| 老司机午夜福利在线观看视频| 人妻久久中文字幕网| 免费观看精品视频网站| 国产免费男女视频| 婷婷成人精品国产| 天天躁夜夜躁狠狠躁躁| 国产在线观看jvid| 建设人人有责人人尽责人人享有的| 精品电影一区二区在线| 美女视频免费永久观看网站| 韩国精品一区二区三区| 久久久国产欧美日韩av| 天堂动漫精品| 女人久久www免费人成看片| 午夜福利在线免费观看网站| 1024视频免费在线观看| 亚洲av欧美aⅴ国产| 欧美日本中文国产一区发布| 久久久久精品人妻al黑| 亚洲精品自拍成人| 亚洲成a人片在线一区二区| 国产欧美日韩综合在线一区二区| aaaaa片日本免费| 久久久久久久久免费视频了| 9热在线视频观看99| 日本五十路高清| 亚洲第一青青草原| 99热只有精品国产| 1024视频免费在线观看| 色综合欧美亚洲国产小说| 亚洲av片天天在线观看| 成年版毛片免费区| 精品国产超薄肉色丝袜足j| 亚洲精品在线观看二区| 国产免费av片在线观看野外av| 国产高清激情床上av| 久久狼人影院| 欧美午夜高清在线| 19禁男女啪啪无遮挡网站| av网站免费在线观看视频| 国产成人av教育| 91麻豆av在线| 国产精品一区二区免费欧美| 亚洲av成人一区二区三| 激情在线观看视频在线高清 | 久久久久久久久久久久大奶| tocl精华| 欧美黄色淫秽网站| 看片在线看免费视频| 亚洲三区欧美一区| 免费高清在线观看日韩| 多毛熟女@视频| 亚洲成人免费电影在线观看| 日韩三级视频一区二区三区| 久久久久国产一级毛片高清牌| 欧美激情久久久久久爽电影 | 欧美精品高潮呻吟av久久| cao死你这个sao货| 黑人欧美特级aaaaaa片| 婷婷精品国产亚洲av在线 | 一区二区三区激情视频| 人人妻人人澡人人爽人人夜夜| e午夜精品久久久久久久| 老司机亚洲免费影院| 久久亚洲精品不卡| 久99久视频精品免费| 亚洲欧美一区二区三区久久| 另类亚洲欧美激情| 精品国产一区二区久久| 老司机亚洲免费影院| 亚洲国产毛片av蜜桃av| 91国产中文字幕| 窝窝影院91人妻| 人人妻,人人澡人人爽秒播| 午夜亚洲福利在线播放| 精品一区二区三区av网在线观看| videos熟女内射| 亚洲黑人精品在线| 亚洲色图综合在线观看| 亚洲成人免费av在线播放| 国产免费现黄频在线看| 麻豆国产av国片精品| 久久久国产成人精品二区 | 欧美日韩黄片免| 亚洲国产精品sss在线观看 | 曰老女人黄片| 身体一侧抽搐| 99re在线观看精品视频| 1024香蕉在线观看| 国产成人精品久久二区二区91| 国产精品成人在线| 老司机午夜福利在线观看视频| 最新美女视频免费是黄的| 亚洲自偷自拍图片 自拍| 国产99白浆流出| 国产成人系列免费观看| 在线天堂中文资源库| 亚洲国产欧美一区二区综合| 久久这里只有精品19| 伦理电影免费视频| 亚洲国产精品合色在线| 亚洲国产毛片av蜜桃av| 欧美乱妇无乱码| 黑人操中国人逼视频| 久久亚洲精品不卡| 亚洲欧美激情综合另类| 亚洲成人免费电影在线观看| 国产有黄有色有爽视频| 国产在线精品亚洲第一网站| 国产不卡一卡二| 青草久久国产| 欧美国产精品va在线观看不卡| 国产欧美日韩一区二区三| 一进一出抽搐动态| 老司机在亚洲福利影院| 午夜福利在线免费观看网站| 黄片播放在线免费| 波多野结衣av一区二区av| 免费看十八禁软件| 亚洲一区二区三区不卡视频| 99精品欧美一区二区三区四区| x7x7x7水蜜桃| 久久香蕉精品热| 久久精品人人爽人人爽视色| 久久国产精品人妻蜜桃| 国产又爽黄色视频| 在线观看日韩欧美| 一级a爱视频在线免费观看| 在线天堂中文资源库| 久久草成人影院| 超碰97精品在线观看| 丝袜美腿诱惑在线| 国产成人av激情在线播放| 女人高潮潮喷娇喘18禁视频| 成人黄色视频免费在线看| 无限看片的www在线观看| 免费在线观看完整版高清| 99热网站在线观看| 悠悠久久av| 天堂中文最新版在线下载| 50天的宝宝边吃奶边哭怎么回事| 亚洲美女黄片视频| 国产色视频综合| 久久精品国产综合久久久| 80岁老熟妇乱子伦牲交| 国产成人欧美在线观看 | 色婷婷av一区二区三区视频| 身体一侧抽搐| 亚洲黑人精品在线| 久久国产亚洲av麻豆专区| 久久狼人影院| 亚洲免费av在线视频| 午夜两性在线视频| 交换朋友夫妻互换小说| 久久精品国产综合久久久| 国产不卡av网站在线观看| 国产精品秋霞免费鲁丝片| 又紧又爽又黄一区二区| 性色av乱码一区二区三区2| 淫妇啪啪啪对白视频| 色在线成人网| 亚洲人成电影免费在线| 国产精品亚洲一级av第二区| 最近最新中文字幕大全免费视频| 窝窝影院91人妻| 日韩欧美一区视频在线观看| e午夜精品久久久久久久| 丰满饥渴人妻一区二区三| 丝袜美腿诱惑在线| 一区二区三区国产精品乱码| 国产精品二区激情视频| 国产高清videossex| 每晚都被弄得嗷嗷叫到高潮| 欧美激情高清一区二区三区| 午夜免费鲁丝| 午夜免费观看网址| 亚洲avbb在线观看| 亚洲美女黄片视频| 99精品久久久久人妻精品| 9热在线视频观看99| 国产成人精品在线电影| 亚洲 欧美一区二区三区| 午夜福利欧美成人| 中文字幕色久视频| 男女午夜视频在线观看| 热99re8久久精品国产| 黄色成人免费大全| 久久人妻熟女aⅴ| 精品久久久久久久久久免费视频 | 亚洲av成人av| 久久国产精品人妻蜜桃| 精品一品国产午夜福利视频| 18禁裸乳无遮挡动漫免费视频| 久久婷婷成人综合色麻豆| 国产成人精品久久二区二区91| 少妇猛男粗大的猛烈进出视频| 他把我摸到了高潮在线观看| 国产精品免费大片| 久久国产精品男人的天堂亚洲| 18禁黄网站禁片午夜丰满| 一边摸一边做爽爽视频免费| 亚洲黑人精品在线| 久久久久国产精品人妻aⅴ院 | 女性被躁到高潮视频| 欧美日韩中文字幕国产精品一区二区三区 | 久久精品国产a三级三级三级| 又大又爽又粗| 成年版毛片免费区| 18禁美女被吸乳视频| 真人做人爱边吃奶动态| 亚洲av成人av| av天堂久久9| 巨乳人妻的诱惑在线观看| 免费在线观看黄色视频的| 又紧又爽又黄一区二区| 久久久久久久久久久久大奶| 免费在线观看视频国产中文字幕亚洲| 后天国语完整版免费观看| 岛国在线观看网站| 日本a在线网址| 丝袜美足系列| 老汉色av国产亚洲站长工具| 亚洲情色 制服丝袜| av在线播放免费不卡| 精品一区二区三区四区五区乱码| 久久久久久免费高清国产稀缺| 母亲3免费完整高清在线观看| netflix在线观看网站| 91麻豆精品激情在线观看国产 | 亚洲一区中文字幕在线| 欧美激情高清一区二区三区| 中文字幕高清在线视频| 久久久水蜜桃国产精品网| 国产高清国产精品国产三级| 国产又爽黄色视频| 欧美一级毛片孕妇| tube8黄色片| 国产极品粉嫩免费观看在线| 12—13女人毛片做爰片一| 国产精品久久视频播放| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av成人av| 怎么达到女性高潮| 精品人妻熟女毛片av久久网站| 欧洲精品卡2卡3卡4卡5卡区| 国产不卡一卡二| 在线av久久热| 欧美精品高潮呻吟av久久| 极品人妻少妇av视频| 久久婷婷成人综合色麻豆| 亚洲人成电影观看| 黑人巨大精品欧美一区二区mp4| 国产成人精品无人区| 在线观看舔阴道视频| 天堂√8在线中文| 91大片在线观看| 成年动漫av网址| 精品欧美一区二区三区在线| 欧美在线黄色| av网站在线播放免费| 午夜视频精品福利| 日韩熟女老妇一区二区性免费视频| 天堂俺去俺来也www色官网| 手机成人av网站| 天堂俺去俺来也www色官网| 久久久国产精品麻豆| 国产精品久久电影中文字幕 | 亚洲第一青青草原| 免费在线观看亚洲国产| 国产精品98久久久久久宅男小说| 免费在线观看完整版高清| 一进一出抽搐动态| 欧美精品啪啪一区二区三区| 啦啦啦免费观看视频1| 51午夜福利影视在线观看| 在线播放国产精品三级| 久久久久久久国产电影| 亚洲av日韩在线播放| 91老司机精品| 亚洲美女黄片视频| 十八禁人妻一区二区| 亚洲国产欧美网| av视频免费观看在线观看| 欧美日韩一级在线毛片| 国产精品乱码一区二三区的特点 | 亚洲一区中文字幕在线| 老汉色av国产亚洲站长工具| 中文字幕色久视频| 精品一区二区三区视频在线观看免费 | 又紧又爽又黄一区二区| 国产伦人伦偷精品视频| 免费日韩欧美在线观看| 亚洲欧美一区二区三区黑人| 色老头精品视频在线观看| 男女免费视频国产| 亚洲一区中文字幕在线| 久久久久久久精品吃奶| 欧美日韩精品网址| 亚洲一区二区三区欧美精品| 日韩中文字幕欧美一区二区| 制服诱惑二区| 91字幕亚洲| 视频区图区小说| 国产精品久久久久久精品古装| 黄频高清免费视频| 91精品三级在线观看| 久久中文字幕一级| 最新在线观看一区二区三区| 热99国产精品久久久久久7| 天天影视国产精品| 欧美人与性动交α欧美精品济南到| 日本黄色视频三级网站网址 | 色婷婷久久久亚洲欧美| 天天操日日干夜夜撸| 757午夜福利合集在线观看| 天堂中文最新版在线下载| 热99国产精品久久久久久7| 狠狠狠狠99中文字幕| 又黄又爽又免费观看的视频| 男女下面插进去视频免费观看| 热99国产精品久久久久久7| 亚洲国产精品一区二区三区在线| 国产精品一区二区精品视频观看| 十八禁人妻一区二区| 建设人人有责人人尽责人人享有的| 亚洲五月婷婷丁香| ponron亚洲| 成年人黄色毛片网站| 国产成人欧美| 一进一出抽搐动态| 午夜福利视频在线观看免费| 国产亚洲欧美精品永久| 午夜免费观看网址| 欧美最黄视频在线播放免费 | 亚洲aⅴ乱码一区二区在线播放 | 日本黄色日本黄色录像| 亚洲成人手机| 少妇裸体淫交视频免费看高清 | 老熟妇仑乱视频hdxx| 久99久视频精品免费| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品乱久久久久久| 亚洲国产欧美日韩在线播放| 精品国产一区二区三区久久久樱花| 久久久久久亚洲精品国产蜜桃av| 91字幕亚洲| 成人手机av| 又紧又爽又黄一区二区| 精品国产乱码久久久久久男人| 国产高清国产精品国产三级| 丝瓜视频免费看黄片| 91成年电影在线观看| 99久久综合精品五月天人人| 欧美日韩亚洲高清精品| 一边摸一边抽搐一进一小说 | 欧美日韩亚洲综合一区二区三区_| 一级毛片女人18水好多| 国产蜜桃级精品一区二区三区 | tube8黄色片| 国产精品美女特级片免费视频播放器 | 99国产极品粉嫩在线观看| 成年人黄色毛片网站| 久久国产精品影院| 两性夫妻黄色片| 每晚都被弄得嗷嗷叫到高潮| 操美女的视频在线观看| av欧美777| 成人18禁高潮啪啪吃奶动态图| www.999成人在线观看| av国产精品久久久久影院| 在线免费观看的www视频| 精品久久蜜臀av无| 黄色女人牲交| 精品少妇久久久久久888优播| 欧美激情极品国产一区二区三区| 久久精品国产综合久久久| 久久精品国产亚洲av高清一级| av免费在线观看网站| 国产精品偷伦视频观看了| 飞空精品影院首页| 十分钟在线观看高清视频www| 精品视频人人做人人爽| 精品久久久久久久久久免费视频 | 国产成人精品久久二区二区91| 一级黄色大片毛片| 久9热在线精品视频| 婷婷精品国产亚洲av在线 | 天堂动漫精品| 欧美乱码精品一区二区三区| 一级片免费观看大全| 999久久久精品免费观看国产| 日本一区二区免费在线视频| 巨乳人妻的诱惑在线观看| 国产片内射在线| 日本黄色日本黄色录像| 一本综合久久免费| 免费看十八禁软件| 天天躁夜夜躁狠狠躁躁| 色综合欧美亚洲国产小说| 亚洲精品一二三| 中文字幕av电影在线播放| 老司机影院毛片| 久久久久国产精品人妻aⅴ院 | 精品人妻1区二区| 成熟少妇高潮喷水视频| 欧美黑人欧美精品刺激| 91九色精品人成在线观看| 亚洲黑人精品在线| 悠悠久久av| 欧美丝袜亚洲另类 | 精品久久久久久,| 黑人猛操日本美女一级片| 亚洲av成人av| 国产精品久久久久成人av| 又大又爽又粗| 亚洲av成人不卡在线观看播放网| 亚洲综合色网址| 1024香蕉在线观看| 精品卡一卡二卡四卡免费| 人成视频在线观看免费观看| 黄色 视频免费看| x7x7x7水蜜桃| 黄色女人牲交| 一级毛片高清免费大全| 人成视频在线观看免费观看| 欧美 日韩 精品 国产| 18禁裸乳无遮挡动漫免费视频| 丰满的人妻完整版| 亚洲成人手机| 中亚洲国语对白在线视频| 悠悠久久av| 久久久久视频综合| 国产精品免费一区二区三区在线 | 亚洲久久久国产精品| 久久精品国产综合久久久| 电影成人av| 国产精品综合久久久久久久免费 | 精品福利永久在线观看| 12—13女人毛片做爰片一| 国产单亲对白刺激| 高清视频免费观看一区二区| 亚洲专区字幕在线| 欧美亚洲 丝袜 人妻 在线| av线在线观看网站| 老司机午夜福利在线观看视频| 人妻 亚洲 视频| 91老司机精品| 免费高清在线观看日韩| 色老头精品视频在线观看| a在线观看视频网站| 999精品在线视频| 精品久久久久久电影网| 一a级毛片在线观看| 久久人妻av系列| 国产欧美亚洲国产| 国产精品久久久久成人av| 久久精品91无色码中文字幕| 成人av一区二区三区在线看| 在线av久久热| 国产精品免费一区二区三区在线 | 亚洲欧美色中文字幕在线| 国产精品永久免费网站| 91精品三级在线观看| 精品国产国语对白av| 大型av网站在线播放| 天堂俺去俺来也www色官网| 一进一出抽搐动态| 欧美在线黄色| 精品久久久久久久久久免费视频 | 天堂√8在线中文| 欧美乱色亚洲激情| 色精品久久人妻99蜜桃| 一进一出抽搐gif免费好疼 | xxx96com| 久久国产精品大桥未久av| 亚洲午夜精品一区,二区,三区| 青草久久国产| 男女午夜视频在线观看| 高清在线国产一区| 国产精品.久久久| 老司机亚洲免费影院| 亚洲五月婷婷丁香| 91成人精品电影| 黄色丝袜av网址大全| 19禁男女啪啪无遮挡网站| 国产精品久久久久成人av| 又大又爽又粗| 亚洲av成人不卡在线观看播放网| 国产男女内射视频| 精品人妻熟女毛片av久久网站| 国产亚洲精品一区二区www | 久久久精品区二区三区| 极品人妻少妇av视频| 国产主播在线观看一区二区| 亚洲欧美日韩高清在线视频| 欧美激情高清一区二区三区| 一边摸一边做爽爽视频免费| 精品一区二区三区四区五区乱码| 亚洲精品久久成人aⅴ小说| 淫妇啪啪啪对白视频| 丝袜美腿诱惑在线| 久久香蕉精品热| 免费在线观看亚洲国产| av视频免费观看在线观看| 丁香六月欧美| 亚洲av日韩在线播放| 超碰97精品在线观看| 欧美黄色淫秽网站| 亚洲成人国产一区在线观看| 精品卡一卡二卡四卡免费| 国产精品 国内视频| 三上悠亚av全集在线观看| 成人三级做爰电影| 大香蕉久久网| 9热在线视频观看99| av福利片在线| 一a级毛片在线观看| 国产99久久九九免费精品| 亚洲欧美一区二区三区黑人| 国产一区二区三区视频了| 亚洲成国产人片在线观看| 99热只有精品国产| 国产在视频线精品| 99精品久久久久人妻精品| 国产免费av片在线观看野外av| 香蕉久久夜色| 后天国语完整版免费观看| 亚洲片人在线观看| 丝袜美足系列| 亚洲黑人精品在线| 热re99久久精品国产66热6| 国产一区二区激情短视频| 日本vs欧美在线观看视频| 国产精品影院久久| 亚洲第一欧美日韩一区二区三区| 国产精品香港三级国产av潘金莲| 国产日韩一区二区三区精品不卡| 一二三四社区在线视频社区8| 亚洲精品av麻豆狂野| 极品教师在线免费播放| 国产在视频线精品| 大型av网站在线播放| 水蜜桃什么品种好| 中文字幕人妻熟女乱码| 久久天堂一区二区三区四区| 亚洲va日本ⅴa欧美va伊人久久| 18禁黄网站禁片午夜丰满| 欧美激情 高清一区二区三区| 欧美激情久久久久久爽电影 | 99国产精品一区二区蜜桃av | 午夜久久久在线观看| 美女高潮喷水抽搐中文字幕| netflix在线观看网站| 久久精品成人免费网站| 一级a爱片免费观看的视频| 在线观看免费高清a一片| 久久99一区二区三区| 999精品在线视频| 狂野欧美激情性xxxx| 国产免费现黄频在线看|