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

    浮筒式波浪發(fā)電平臺(tái)性能分析研究

    2016-10-11 07:25:17李德堂金豁然胡星辰趙春慧
    船舶與海洋工程 2016年4期
    關(guān)鍵詞:浮筒海況輸出功率

    呂 沁,李德堂,金豁然,胡星辰,魏 卓,趙春慧

    (1. 浙江海洋大學(xué) 船舶與海洋工程學(xué)院,浙江 舟山 316000;2. 浙江省近海海洋工程技術(shù)重點(diǎn)實(shí)驗(yàn)室,浙江 舟山 316000)

    浮筒式波浪發(fā)電平臺(tái)性能分析研究

    呂 沁1,2,李德堂1,2,金豁然1,2,胡星辰1,2,魏 卓1,2,趙春慧1,2

    (1. 浙江海洋大學(xué) 船舶與海洋工程學(xué)院,浙江 舟山 316000;2. 浙江省近海海洋工程技術(shù)重點(diǎn)實(shí)驗(yàn)室,浙江 舟山 316000)

    波浪能發(fā)電浮筒的研究關(guān)鍵是實(shí)現(xiàn)波浪能發(fā)電裝置的能量轉(zhuǎn)換最大化。根據(jù)“海院1號(hào)”浮筒式波浪發(fā)電平臺(tái)的運(yùn)動(dòng)特點(diǎn),提出了一種預(yù)測(cè)浮子垂蕩運(yùn)動(dòng)的時(shí)域方法。同時(shí),通過MATLAB模擬計(jì)算,獲得最佳系統(tǒng)壓力和平均傳輸功率。經(jīng)海測(cè)試驗(yàn)驗(yàn)證了數(shù)值模擬的正確性,為波浪發(fā)電平臺(tái)浮筒設(shè)計(jì)優(yōu)化提供依據(jù),具有工程應(yīng)用價(jià)值。

    波浪能;浮筒;MATLAB;平均輸出功率;取能效率

    0 引 言

    近年來,隨著社會(huì)的不斷發(fā)展,能源需求量快速增加,化石能源短缺、環(huán)境污染等問題逐漸成為世界性問題[1]。在當(dāng)下能源日趨緊張的情況下,尋找清潔、優(yōu)質(zhì)的新型能源成為人們關(guān)心的課題。占地球表面71%的海洋蘊(yùn)藏著巨大的可再生能源[2],主要包括海洋風(fēng)能、潮汐能、波浪能、海流能、溫差能和鹽差能等,這些新型海洋可再生能源的開發(fā)利用已成為研究熱點(diǎn)[3-5]。

    海洋波浪能是一種清潔的新型海洋能源,并且儲(chǔ)量十分豐富。波浪能主要由風(fēng)能傳遞形成。世界能源委員會(huì)的調(diào)查顯示,全球可利用的波浪能達(dá)到20億kW,相當(dāng)于目前世界發(fā)電能量的2倍[6]。我國擁有廣闊的海洋資源,波浪能的理論存儲(chǔ)量約為700萬kW,沿海波浪能能流密度為2~7kW/m[7-9]。有關(guān)數(shù)據(jù)顯示,我國浙江、福建、廣東和臺(tái)灣沿海均為波浪能豐富的地區(qū),其中浙江省所蘊(yùn)含的波浪能資源有巨大的開發(fā)潛力[10]。加快儲(chǔ)量豐富的波浪能資源的開發(fā)利用,將有效地緩解我國的能源短缺和環(huán)境污染等問題,具有重要的經(jīng)濟(jì)意義與社會(huì)意義。

    “海院1號(hào)”波浪發(fā)電平臺(tái)是浙江海洋大學(xué)船舶與海洋工程學(xué)院自主研發(fā)設(shè)計(jì)的自升式點(diǎn)吸收式波浪發(fā)電平臺(tái)(見圖1)。該平臺(tái)具有電力輸出穩(wěn)定、裝置可靠性好、發(fā)電效率高、管理和維護(hù)成本低等優(yōu)點(diǎn)。該波浪發(fā)電平臺(tái)主要由支撐平臺(tái)、3組波浪能收集模塊、液壓系統(tǒng)和電力輸出系統(tǒng)組成,其主要結(jié)構(gòu)見圖2和圖3。浮筒以樁柱為導(dǎo)向,上端裝有2塊波浪板,每塊波浪板通過帶有滾輪的群組油缸相配合。波浪能由浮筒獲取,因此浮筒的設(shè)計(jì)將直接影響到波浪發(fā)電平臺(tái)的發(fā)電效率,浮筒發(fā)電性能的研究對(duì)整個(gè)平臺(tái)至關(guān)重要。

    圖2 發(fā)電平臺(tái)整體結(jié)構(gòu)

    圖3 發(fā)電平臺(tái)能量轉(zhuǎn)換裝置

    采用MATLAB編程設(shè)計(jì),對(duì)“海院1號(hào)”波浪發(fā)電平臺(tái)浮筒在特定海況下進(jìn)行數(shù)值模擬,得到平均發(fā)電輸出功率與取能效率,分析發(fā)電效果。通過“海院1號(hào)”海上實(shí)測(cè)試驗(yàn),與該浮筒的數(shù)值模擬結(jié)果進(jìn)行比較,驗(yàn)證數(shù)值模擬的正確性。

    1 研究對(duì)象

    “海院1號(hào)”波浪發(fā)電平臺(tái)選取的發(fā)電浮筒是一個(gè)旋轉(zhuǎn)體(簡(jiǎn)稱浮筒1),其形狀見圖4。其結(jié)構(gòu)分為2部分:上半部分是一個(gè)直徑為3200mm,高度為625mm的圓柱體;下半部分是一個(gè)類似于半球的曲面,高度為 835mm。浮筒內(nèi)部設(shè)有加強(qiáng)結(jié)構(gòu)。

    圖4 發(fā)電平臺(tái)浮筒1尺度(mm)

    2 數(shù)學(xué)模型

    2.1 浮筒垂蕩運(yùn)動(dòng)方程

    規(guī)定浮筒在水面上達(dá)到靜平衡時(shí)的水線面中心為坐標(biāo)原點(diǎn),z軸豎直向上,x軸沿入射波方向?yàn)檎?/p>

    由于波浪發(fā)電裝置浮筒的運(yùn)動(dòng)具有復(fù)雜的非線性,因此采用一種時(shí)域法;規(guī)定入射波浪力為未受擾動(dòng)的入射波對(duì)浮筒產(chǎn)生的浮力,并且由于液壓系統(tǒng)阻尼力的作用,水動(dòng)力系數(shù)隨著濕表面變化而改變。在浮筒運(yùn)動(dòng)軌跡的計(jì)算中,既要考慮輻射系數(shù)、附加質(zhì)量系數(shù)及回復(fù)力系數(shù)與濕表面的變化關(guān)系,還要考慮濕表面變化所引起的浮力變化。根據(jù)上述方法,并結(jié)合船舶在波浪上運(yùn)動(dòng)的時(shí)域計(jì)算方法,可得到波浪發(fā)電裝置浮筒的垂蕩運(yùn)動(dòng)方程[11-12]為

    式(1)中:左邊第一項(xiàng)表示慣性力,其中,M為浮筒的質(zhì)量,z(t)為浮筒的垂直位移;左邊第二項(xiàng)表示附加質(zhì)量力,其相反數(shù)是與加速度成線性關(guān)系的輻射力,其中,ξ(x,y,t)為(x,y)處入射波引起的波面抬升,m∞為頻率趨向于無窮大時(shí)浮筒的附加質(zhì)量(由于浮筒濕表面隨著時(shí)間的變化而變化,故m∞為時(shí)間的函數(shù));左邊第三項(xiàng)表示阻尼力,其相反數(shù)是與速度成線性關(guān)系的輻射力,其中,h表示記憶效應(yīng)的有效時(shí)長,一重積分表示輻射的記憶效應(yīng),函數(shù)L(t)稱為記憶函數(shù)(也稱作時(shí)延函數(shù))。L(t)表示浮筒由于自由產(chǎn)生記憶效應(yīng)帶來的影響,其表達(dá)式[13]為

    式(2)中:B為輻射阻尼系數(shù),其值與入射波頻率以及相應(yīng)時(shí)刻浮筒的濕表面面積有關(guān),由于浮筒濕表面隨著時(shí)間的變化而變化,故B為時(shí)間的函數(shù);ε(ω) 為頻率ω的入射波初始相位;左邊第二項(xiàng)附加質(zhì)量力與第三項(xiàng)阻尼力共同組成輻射力;右邊第一項(xiàng)為瞬時(shí)浮力,其中S(t)為t時(shí)刻浮筒的濕表面面積;n為濕表面的外法向量;ρ為海水密度,取1025kg/m3;g為重力加速度,取9.81m/s2。而對(duì)濕表面的二重積分項(xiàng)表示在t時(shí)刻未受擾動(dòng)的入射波對(duì)浮筒產(chǎn)生的浮力;Mg為浮筒的重力;FPTO為液壓阻尼力。

    2.2 液壓阻尼力

    “海院1號(hào)”采用了雙向液壓系統(tǒng),并且使用了液壓蓄能器,其系統(tǒng)壓力穩(wěn)定,因此浮筒所受到的液壓阻尼力FPTO[14]可表示為

    式(3)中:Fr為廣義輻射力;v為浮筒的速度;k1為常數(shù),由系統(tǒng)壓力p和液壓活塞面積Ap決定,其中系統(tǒng)壓力p可以通過平臺(tái)液壓系統(tǒng)裝置人為控制,因此可以計(jì)算不同系統(tǒng)壓力下的液壓阻尼力。液壓活塞面積Ap為100cm2,即

    2.3 平均傳輸功率

    浮子瞬時(shí)傳輸功率 P (t)和一段時(shí)長 Δt浮筒的平均傳輸功率分別為

    3 數(shù)值模擬法

    3.1 參數(shù)選定

    在數(shù)值模擬中,海況數(shù)據(jù)采用浙江省舟山市普陀區(qū)朱家尖鎮(zhèn)東沙海域的數(shù)據(jù)。根據(jù)2014年4月20日的平臺(tái)實(shí)際測(cè)量,取浮筒吃水為0.835m,水深2.17m,波浪的最大波高為0.80m,有義波高≈0.65m,平均周期≈6.0s。

    數(shù)值模擬實(shí)際海況下的波浪情況有一定的難度,現(xiàn)認(rèn)定實(shí)際海況下的波浪由5個(gè)余弦波疊加而成[15]。為了使模擬數(shù)據(jù)與實(shí)際海況數(shù)據(jù)接近,選取5個(gè)余弦波的參數(shù) Ai,ωi和 εi分別為:A1= 0.05m,ω1= 2.0rad/s,ε1= 1.75;A2= 0.01m, ω2= 1.2rad/s,ε2= 2.25;A3= 0.05m, ω3= 1.4rad/s,ε3= 3.75;A4= 0.02m,ω4= 1.7rad/s,ε4= 4.25;A5= 0.33m, ω5= 1.05rad/s,ε5= 5.75。

    為描述不規(guī)則波的隨機(jī)性,初相位εi可選用均勻分布于[0,2π]中的隨機(jī)量[16]。其波浪的疊加公式為

    通過計(jì)算,模擬波浪的最大波高為0.805m,有義波高為0.64m,平均周期為5.98s,與實(shí)際海況波浪比較接近,選定以上數(shù)據(jù)為波浪參數(shù)較合理。

    3.2 數(shù)值求解方法

    由于浮筒受力和各個(gè)水動(dòng)力系數(shù)變化的復(fù)雜性,無法得到FPTO的解析解,因此采用數(shù)值方法求解。采用顯示改進(jìn)歐拉法對(duì)運(yùn)動(dòng)方程和記憶函數(shù)做離散化處理,根據(jù)所選參數(shù)與初始條件,通過MATLAB軟件編程,得到運(yùn)動(dòng)方程的解,最后根據(jù)式(6)得到浮筒在不同系統(tǒng)壓力下的平均輸出功率,其中規(guī)定初始條件為:

    數(shù)值模擬中,系統(tǒng)壓力選取為0~5MPa范圍內(nèi)以0.25MPa為間距的點(diǎn)。浮筒在每個(gè)選定的系統(tǒng)壓力上,計(jì)算一次平均輸出功率,其中計(jì)算時(shí)間間隔Δt設(shè)定為80s。最后,將浮筒不同系統(tǒng)壓力下的平均輸出功率繪制成曲線圖,并進(jìn)行浮筒取能效率的計(jì)算,定量與定性地分析浮筒的發(fā)電性能。單次數(shù)值模擬程序流程見圖5。

    圖5 MATLAB模擬程序流程

    4 數(shù)值模擬結(jié)果

    4.1 數(shù)值模擬結(jié)果

    通過MATLAB的數(shù)值模擬,可以得到浮筒系統(tǒng)壓力與平均輸出功率的關(guān)系曲線(見圖6);設(shè)定系統(tǒng)壓力為橫坐標(biāo),可以得到平臺(tái)發(fā)電效果最佳的系統(tǒng)壓力。

    圖6為浮筒在一段時(shí)間內(nèi)的平均輸出功率與系統(tǒng)壓力的關(guān)系圖,根據(jù) MATLAB數(shù)值模擬得到。當(dāng)系統(tǒng)壓力從0MPa增加時(shí),平均輸出功率也隨之增大,并且呈現(xiàn)一定的線性關(guān)系;當(dāng)系統(tǒng)壓力達(dá)到1.6MPa時(shí),平均輸出功率增加速度變慢,并且當(dāng)系統(tǒng)壓力達(dá)到2.0MPa時(shí),平均輸出功率達(dá)到最大值2425W;在平均輸出功率達(dá)到峰值之后,其值隨著系統(tǒng)壓力的增大開始減小;在系統(tǒng)壓力達(dá)到4MPa時(shí),平均輸出功率減慢幅度降低,逐漸達(dá)到穩(wěn)定。由圖6可知,浮筒的最大輸出功率為2425W。

    圖6 浮筒的系統(tǒng)壓力與平均輸出功率關(guān)系

    4.2 取能效率計(jì)算

    浮筒的取能效率是浮筒獲得的平均輸出功率與波浪總能量的比值,又稱為一級(jí)轉(zhuǎn)化效率,其中波浪總能量為

    式(9)中:ρ,g,T為海水密度、重力加速度及周期;H為有義波高; Cg為波速(淺水條件下為,d為水深);Aω為浮筒的水線面面積。

    通過計(jì)算可以得到模擬波浪能量為3807W。從而可以得到,當(dāng)系統(tǒng)壓力達(dá)到2.0MPa時(shí),浮筒的取能效率最大,為63.7%,可滿足發(fā)電平臺(tái)的需求。

    5 平臺(tái)海測(cè)試驗(yàn)驗(yàn)證

    實(shí)際海況下“海院1號(hào)”波浪發(fā)電平臺(tái)浮筒的取能效率可以通過波浪和浮筒運(yùn)動(dòng)數(shù)據(jù)得到。試驗(yàn)中,使用YF-YJ50激光位移傳感器記錄試驗(yàn)浮筒的運(yùn)動(dòng),得到浮筒運(yùn)動(dòng)曲線圖;同時(shí),通過YF-YJ50激光位移傳感器測(cè)試放置在海面上浮標(biāo)的運(yùn)動(dòng),得到波浪運(yùn)動(dòng)曲線圖,并且控制浮標(biāo)與浮筒的波浪入射方向、入射時(shí)刻一致?,F(xiàn)選取波浪較為穩(wěn)定的一個(gè)時(shí)間段內(nèi)的數(shù)據(jù)驗(yàn)證數(shù)值模擬結(jié)果的準(zhǔn)確性。

    圖7為選取的實(shí)際海況波浪運(yùn)動(dòng)曲線圖,測(cè)試時(shí)間是80s,水深2.17m,浮筒最大吃水0.835m,系統(tǒng)壓力1.8MPa。根據(jù)上跨零點(diǎn)法對(duì)圖7進(jìn)行分析,可以得到這個(gè)時(shí)間段內(nèi)波浪的最大波高為0.802m,有義波高為0.61m,波浪平均周期為6.15s。根據(jù)式(9)得到波浪的總能量為2820W。

    圖8為相對(duì)應(yīng)的浮筒運(yùn)動(dòng)曲線圖。由該圖可知,浮筒的運(yùn)動(dòng)總行程約為5.78m。由于“海院1號(hào)”波浪發(fā)電平臺(tái)采用雙向液壓系統(tǒng),即浮筒上下運(yùn)動(dòng)行程都是有效行程,根據(jù)式(3)和式(4)可以得到液壓阻尼力為18000N,根據(jù)式(10)得到浮筒的平均輸出功率為1301W。

    圖7 實(shí)際海況波浪

    圖8 實(shí)際海況浮筒運(yùn)動(dòng)

    根據(jù)上述計(jì)算結(jié)果,可以得到實(shí)際海況下的浮筒取能效率為46.1%。由圖6可知,浮筒在系統(tǒng)壓力為1.8MPa時(shí)的取能效率約為63.1%。數(shù)值模擬結(jié)果與實(shí)際海況下該浮筒取能效率存在一定的差距。但考慮到實(shí)際海況下具有較多不確定性的因素及液壓系統(tǒng)存在誤差,各個(gè)部件之間存在摩擦;而數(shù)值模擬較為理想化,模擬結(jié)果與實(shí)際測(cè)試存在17%的誤差合乎常理,數(shù)值模擬仍符合實(shí)際要求。由此,判定采用本文數(shù)值模擬分析浮筒的發(fā)電性能較合理。

    6 結(jié) 語

    通過MATLAB數(shù)值模擬“海院1號(hào)”浮筒式波浪發(fā)電平臺(tái)在某一海況條件下的輸出功率與取能效率,為波浪發(fā)電的浮筒設(shè)計(jì)優(yōu)化提供幫助。并且通過對(duì)平臺(tái)實(shí)際海況測(cè)試中負(fù)載變化的觀察分析和取能效率的計(jì)算,與模擬數(shù)據(jù)進(jìn)行比較,驗(yàn)證了數(shù)值模擬的正確性。從數(shù)值模擬與實(shí)際海況測(cè)試可以得到:

    1) 平臺(tái)浮筒具有良好的發(fā)電性能,可以滿足平臺(tái)的發(fā)電要求;

    2) 浮筒的發(fā)電功率隨液壓系統(tǒng)的系統(tǒng)壓力變化而變化,在系統(tǒng)壓力為2.0MPa時(shí),浮筒的輸出功率與取能效率最佳;

    3) 數(shù)值模擬方法可以較好地預(yù)測(cè)浮筒的輸出功率,具有工程實(shí)用價(jià)值。

    [1] 石陸豐,程建生,段金輝,等. 半潛式浮基風(fēng)電平臺(tái)設(shè)計(jì)及波浪動(dòng)力響應(yīng)分析[J]. 船舶與海洋工程,2015, 31 (1): 13-19.

    [2] 游亞戈,李偉,劉偉民,等. 海洋能發(fā)電技術(shù)的發(fā)展現(xiàn)狀與前景[J]. 電力系統(tǒng)自動(dòng)化,2010, 34 (14): 1-12.

    [3] ISAACS J D, SEYMOUR R J. The ocean as a power resource[J]. International Journal of Environmental Studies, 1973, 4 (3):201-205.

    [4] 王傳昆. 海洋能資源分析方法及儲(chǔ)量評(píng)估[M]. 北京:海洋出版社,2009.

    [5] AFDO F. Wave energy utilization:a review of the technologies[J]. Renewable and Sustainable Energy Reviews, 2010, 14 (3):899-918.

    [6] WICK G L, SCHMITT W R. Harvesting Ocean Energy[M]. Paris: UNESCO,1981.

    [7] 劉美琴,鄭源,趙振宙,等. 海洋能利用的發(fā)展與前景[J]. 海洋開發(fā)與管理,2010, 27(3): 80-82.

    [8] 游亞戈,盛松偉,吳必軍. 海洋波浪能發(fā)電技術(shù)現(xiàn)狀與前景[C]//第十五屆中國海洋(岸)工程學(xué)術(shù)討論會(huì)論文集(上),2011: 9-16.

    [9] 覃嶺,吳必軍. 雙圓柱形浮體波能裝置振動(dòng)特性分析[J]. 船舶力學(xué),2013, 17 (11): 1253-1261.

    [10] 葉欽,楊忠良,施偉勇. 浙江近海波浪能資源的初步研究[J]. 海洋學(xué)研究,2012, 30 (4): 13-19.

    [11] 劉應(yīng)中,繆國平. 船舶在波浪上的運(yùn)動(dòng)理論[M]. 上海:上海交通大學(xué)出版社,1987.

    [12] 胡可一,柳衛(wèi)東,陳兵,等. 垂直船艏在波浪中快速性綜合分析研究[J]. 船舶與海洋工程,2013, 29 (1): 1-5.

    [13] JOURNéE J M J, MASSIE W W. Offshore Hydromechanics[M]. Delft: Delft University of Technology, 2001.

    [14] LI D, LI D T, LI F, et al. Analysis of floating buoy of a wave power generating jack-up platform Haiyuan 1[J]. Advances in Mechanical Engineering, 2015, 5 (12): 105072.

    [15] DEAN R G, DALRYMPLE R A. Water wave mechanics for engineers and Scientists[M].Singapore: World Scientific,1991.

    [16] 張亮,李云波. 流體力學(xué)[M]. 哈爾濱:哈爾濱工程大學(xué)出版社,2006.

    Study on the Performance Analysis of Buoy-type Wave Power Generation Platform

    LYU Qin1,2, LI De-tang1,2, JIN Huo-ran1,2, HU Xing-chen1,2,WEI Zhuo1, 2, ZHAO Chun-hui1,2

    (1. School of Naval Architecture and Ocean Engineering, Zhejiang Ocean University, Zhoushan 316022, China;2. Key Laboratory of Offshore Engineering Technology of Zhejiang Province, Zhoushan 316022, China)

    The key technology of wave power generation buoy is to realize the maximum energy transformation of the wave power generation device. A time-domain method for predicting the vertical motion of the buoy is proposed based on the motion characteristics of the buoy-type wave power generation platform Haiyuan l. Meanwhile, the optimal system pressure and the average delivered power are obtained through MATLAB simulation, and the numerical simulation is verified correct through sea trial. This study provides the basis for the design optimization of wave power generation buoy,and has the practical engineering significance.

    wave energy; buoy; MATLAB; average delivered power; energy harvesting efficiency

    TK79

    A

    2095-4069 (2016) 04-0040-06

    10.14056/j.cnki.naoe.2016.04.009

    2015-08-24

    國家海洋局海洋可再生能源專項(xiàng)資金項(xiàng)目(ZJME2011BL04);上海交通大學(xué)海洋工程國家重點(diǎn)實(shí)驗(yàn)室研究基金資助項(xiàng)目(1205);舟山市科技局資助項(xiàng)目(2014C41013);浙江省自然科學(xué)基金(LY14E090003;Y5100180)2014年國家海洋經(jīng)濟(jì)創(chuàng)新發(fā)展區(qū)域示范項(xiàng)目

    呂沁,男,1990年生。浙江海洋大學(xué)船舶與海洋工程學(xué)院在讀研究生。

    猜你喜歡
    浮筒海況輸出功率
    浸沒式外輸浮筒的動(dòng)力響應(yīng)影響
    一種基于淺水靜水環(huán)境的浮筒構(gòu)架式水上鉆探平臺(tái)
    典型海況下艦載發(fā)射箱結(jié)構(gòu)強(qiáng)度仿真分析
    鋼質(zhì)浮筒在沉船人性化打撈中的應(yīng)用
    中國航海(2018年1期)2018-05-07 11:12:23
    鋼制浮筒在“世越號(hào)”打撈工程中的應(yīng)用
    惡劣海況下海洋石油116內(nèi)轉(zhuǎn)塔式FPSO裝配載優(yōu)化
    船海工程(2015年4期)2016-01-05 15:53:38
    極端海況下軟質(zhì)海崖大規(guī)模蝕退計(jì)算模型研究
    適用于智能電網(wǎng)的任意波形輸出功率源
    基于雙層BP神經(jīng)網(wǎng)絡(luò)的光伏電站輸出功率預(yù)測(cè)
    極限海況下單點(diǎn)系泊系統(tǒng)纜索動(dòng)張力研究
    国产精品欧美亚洲77777| 精品人妻一区二区三区麻豆| 国产免费又黄又爽又色| 妹子高潮喷水视频| 亚洲精品国产成人久久av| 国产av精品麻豆| 亚洲精品中文字幕在线视频 | 国产男人的电影天堂91| 成人免费观看视频高清| 色5月婷婷丁香| 色视频www国产| 久久精品国产亚洲网站| 日韩免费高清中文字幕av| 中文乱码字字幕精品一区二区三区| 国产精品一区www在线观看| 热re99久久精品国产66热6| 尤物成人国产欧美一区二区三区| 亚洲人成网站在线播| 老熟女久久久| 高清毛片免费看| 亚洲国产精品国产精品| 好男人视频免费观看在线| 亚洲伊人久久精品综合| 久久国产精品大桥未久av | 国产精品熟女久久久久浪| 大香蕉97超碰在线| 一本一本综合久久| www.av在线官网国产| 亚洲人与动物交配视频| 女人十人毛片免费观看3o分钟| 丰满乱子伦码专区| 国国产精品蜜臀av免费| 亚洲自偷自拍三级| 精品人妻视频免费看| 欧美三级亚洲精品| 国产乱人视频| 国产精品久久久久久久电影| av在线老鸭窝| 极品教师在线视频| 尤物成人国产欧美一区二区三区| 久久精品久久精品一区二区三区| 边亲边吃奶的免费视频| 内射极品少妇av片p| 又黄又爽又刺激的免费视频.| 国产日韩欧美亚洲二区| 午夜免费鲁丝| 免费人妻精品一区二区三区视频| 80岁老熟妇乱子伦牲交| 亚洲av成人精品一区久久| 全区人妻精品视频| 免费av不卡在线播放| 国产av精品麻豆| 纵有疾风起免费观看全集完整版| 午夜日本视频在线| 亚洲经典国产精华液单| 91精品一卡2卡3卡4卡| 夜夜爽夜夜爽视频| 久久久精品免费免费高清| 久久这里有精品视频免费| 成人综合一区亚洲| 在线观看一区二区三区激情| 特大巨黑吊av在线直播| 久久精品国产亚洲av天美| 久久久久久人妻| 国产高清国产精品国产三级 | 欧美少妇被猛烈插入视频| 久久精品久久久久久噜噜老黄| 两个人的视频大全免费| 久久毛片免费看一区二区三区| 国产高清三级在线| 日日啪夜夜爽| 少妇精品久久久久久久| 中国美白少妇内射xxxbb| 精品国产一区二区三区久久久樱花 | 国产一级毛片在线| 亚洲国产毛片av蜜桃av| 国产高清三级在线| 亚洲综合精品二区| 国产免费一区二区三区四区乱码| 精品久久久久久电影网| 亚洲精品亚洲一区二区| 国产精品偷伦视频观看了| 日本午夜av视频| 在线观看一区二区三区| av在线观看视频网站免费| 日韩,欧美,国产一区二区三区| 国产探花极品一区二区| av线在线观看网站| 黄色视频在线播放观看不卡| 日本一二三区视频观看| 亚洲色图综合在线观看| 噜噜噜噜噜久久久久久91| 91久久精品国产一区二区三区| 国产成人aa在线观看| 赤兔流量卡办理| 2022亚洲国产成人精品| 午夜免费鲁丝| 99久久人妻综合| 欧美最新免费一区二区三区| 国产精品成人在线| 九九久久精品国产亚洲av麻豆| 欧美精品国产亚洲| 日韩不卡一区二区三区视频在线| av免费观看日本| 国产一区二区在线观看日韩| 七月丁香在线播放| 亚州av有码| 久久久久久久精品精品| 嘟嘟电影网在线观看| 日韩视频在线欧美| 欧美三级亚洲精品| 亚洲精品国产色婷婷电影| 观看美女的网站| 国语对白做爰xxxⅹ性视频网站| 亚洲人与动物交配视频| 亚洲精品一区蜜桃| 在现免费观看毛片| 99久国产av精品国产电影| 国产美女午夜福利| 美女国产视频在线观看| 99久久中文字幕三级久久日本| 日韩强制内射视频| av一本久久久久| 视频中文字幕在线观看| 国产精品一区www在线观看| 国产色婷婷99| 国产免费又黄又爽又色| 免费看av在线观看网站| 一本久久精品| 国产精品欧美亚洲77777| 一边亲一边摸免费视频| 亚洲aⅴ乱码一区二区在线播放| 我要看日韩黄色一级片| 超碰av人人做人人爽久久| 日本黄色片子视频| 老司机影院毛片| 精华霜和精华液先用哪个| 久久久国产一区二区| 国产高清有码在线观看视频| 日韩中文字幕视频在线看片 | 亚洲高清免费不卡视频| 91精品伊人久久大香线蕉| 精品久久久久久久末码| 日韩av免费高清视频| 亚洲欧美一区二区三区黑人 | 国产亚洲5aaaaa淫片| 老司机影院毛片| 少妇丰满av| 精品久久国产蜜桃| 久久久久人妻精品一区果冻| 晚上一个人看的免费电影| 欧美成人午夜免费资源| 中文字幕久久专区| 国产精品99久久久久久久久| 亚洲精品aⅴ在线观看| 欧美日韩一区二区视频在线观看视频在线| 亚洲av中文字字幕乱码综合| 成人综合一区亚洲| 在线观看免费视频网站a站| 亚洲av中文字字幕乱码综合| 亚洲欧美成人精品一区二区| 亚洲欧美日韩东京热| 久久人妻熟女aⅴ| 国产精品秋霞免费鲁丝片| 精品午夜福利在线看| 18禁在线播放成人免费| 美女内射精品一级片tv| av女优亚洲男人天堂| 欧美一区二区亚洲| 亚洲av.av天堂| 777米奇影视久久| 国产精品爽爽va在线观看网站| 日韩强制内射视频| 亚洲精品色激情综合| 亚洲成人中文字幕在线播放| 精品国产乱码久久久久久小说| 乱码一卡2卡4卡精品| 男女无遮挡免费网站观看| 亚洲久久久国产精品| 成人美女网站在线观看视频| 欧美一级a爱片免费观看看| 午夜激情福利司机影院| 久久久久网色| 国产中年淑女户外野战色| 热re99久久精品国产66热6| 久久久久国产精品人妻一区二区| 在线免费十八禁| 大香蕉97超碰在线| 女性被躁到高潮视频| 免费大片18禁| 免费不卡的大黄色大毛片视频在线观看| 亚洲国产av新网站| 一级毛片电影观看| 七月丁香在线播放| av免费在线看不卡| 2022亚洲国产成人精品| 亚洲国产av新网站| 欧美另类一区| 成人18禁高潮啪啪吃奶动态图 | 欧美日韩精品成人综合77777| 插逼视频在线观看| 黄色一级大片看看| 日韩一本色道免费dvd| 一区二区三区四区激情视频| 国产伦精品一区二区三区四那| 亚洲成人中文字幕在线播放| 亚洲欧美中文字幕日韩二区| a级一级毛片免费在线观看| 熟女电影av网| 国产精品福利在线免费观看| 舔av片在线| 卡戴珊不雅视频在线播放| 国产精品免费大片| 成人漫画全彩无遮挡| 国产视频首页在线观看| 免费大片18禁| 乱系列少妇在线播放| 久久国产精品男人的天堂亚洲 | 久久午夜福利片| 亚洲精品成人av观看孕妇| 国产日韩欧美亚洲二区| 国产片特级美女逼逼视频| 91aial.com中文字幕在线观看| 一边亲一边摸免费视频| 亚洲va在线va天堂va国产| 婷婷色av中文字幕| 男女无遮挡免费网站观看| 亚洲经典国产精华液单| 七月丁香在线播放| 国产精品久久久久久av不卡| 欧美 日韩 精品 国产| 人妻一区二区av| 日本wwww免费看| 又大又黄又爽视频免费| xxx大片免费视频| 色视频在线一区二区三区| 多毛熟女@视频| 国产探花极品一区二区| 国产黄片美女视频| 久久热精品热| 亚洲欧美精品自产自拍| 精品熟女少妇av免费看| 亚洲激情五月婷婷啪啪| 一级毛片 在线播放| 久久精品久久久久久久性| 寂寞人妻少妇视频99o| 日日摸夜夜添夜夜爱| 一级毛片aaaaaa免费看小| 久久久久精品久久久久真实原创| 欧美丝袜亚洲另类| 亚洲av欧美aⅴ国产| 卡戴珊不雅视频在线播放| 国产午夜精品一二区理论片| 2018国产大陆天天弄谢| 亚洲av成人精品一区久久| 婷婷色综合大香蕉| 婷婷色av中文字幕| av在线观看视频网站免费| 99久国产av精品国产电影| 成人免费观看视频高清| 啦啦啦中文免费视频观看日本| 精品视频人人做人人爽| 最近最新中文字幕免费大全7| 国产精品伦人一区二区| 日本与韩国留学比较| 99久国产av精品国产电影| 国产成人精品福利久久| 超碰av人人做人人爽久久| 日韩电影二区| 自拍欧美九色日韩亚洲蝌蚪91 | 联通29元200g的流量卡| 26uuu在线亚洲综合色| 韩国av在线不卡| 亚洲欧美精品自产自拍| 夜夜骑夜夜射夜夜干| 国产在线视频一区二区| 在线观看美女被高潮喷水网站| 久久精品熟女亚洲av麻豆精品| 精品一区二区免费观看| 精品99又大又爽又粗少妇毛片| 久久久午夜欧美精品| 99久久精品国产国产毛片| 亚洲中文av在线| 欧美一区二区亚洲| 伊人久久国产一区二区| 亚洲一区二区三区欧美精品| 国产精品精品国产色婷婷| 亚洲高清免费不卡视频| 亚洲av成人精品一区久久| 亚洲成人av在线免费| 日产精品乱码卡一卡2卡三| 精品人妻一区二区三区麻豆| 午夜福利高清视频| 伦理电影免费视频| 久久久a久久爽久久v久久| 午夜日本视频在线| 99久久综合免费| 秋霞伦理黄片| 在线观看一区二区三区| 精品少妇久久久久久888优播| 精品亚洲成国产av| 国产成人精品福利久久| 欧美精品国产亚洲| 欧美 日韩 精品 国产| 国产高潮美女av| 日本猛色少妇xxxxx猛交久久| 欧美3d第一页| 欧美激情极品国产一区二区三区 | 激情五月婷婷亚洲| 香蕉精品网在线| 夫妻性生交免费视频一级片| 小蜜桃在线观看免费完整版高清| 国产高清不卡午夜福利| 免费观看a级毛片全部| 国产精品女同一区二区软件| 午夜激情福利司机影院| 一个人免费看片子| 国产有黄有色有爽视频| 国产精品免费大片| 美女内射精品一级片tv| 色婷婷久久久亚洲欧美| 国产精品麻豆人妻色哟哟久久| 观看免费一级毛片| 精品熟女少妇av免费看| 深爱激情五月婷婷| 国产视频内射| 久久99精品国语久久久| 最近的中文字幕免费完整| 永久免费av网站大全| 国内精品宾馆在线| 亚洲av二区三区四区| 午夜激情久久久久久久| 国产精品久久久久久久久免| 久久精品国产亚洲av涩爱| 日本wwww免费看| 老司机影院成人| 中文在线观看免费www的网站| 丝袜喷水一区| 欧美高清成人免费视频www| 日本黄大片高清| 亚洲熟女精品中文字幕| 中文字幕精品免费在线观看视频 | 一区二区av电影网| 国产精品久久久久久久久免| 看非洲黑人一级黄片| 婷婷色av中文字幕| 欧美日韩精品成人综合77777| 国产精品久久久久久精品古装| freevideosex欧美| 日韩免费高清中文字幕av| 日韩国内少妇激情av| 国产精品成人在线| 亚洲精品久久午夜乱码| 夜夜骑夜夜射夜夜干| 夜夜爽夜夜爽视频| 身体一侧抽搐| 久久人人爽人人爽人人片va| 久久韩国三级中文字幕| 最近2019中文字幕mv第一页| 啦啦啦啦在线视频资源| 黄色欧美视频在线观看| av网站免费在线观看视频| 日本猛色少妇xxxxx猛交久久| 中文字幕精品免费在线观看视频 | 五月天丁香电影| 国产久久久一区二区三区| 日日摸夜夜添夜夜添av毛片| 午夜福利在线在线| 久久精品国产亚洲网站| 深夜a级毛片| 男人狂女人下面高潮的视频| 少妇丰满av| 熟妇人妻不卡中文字幕| www.色视频.com| 日日啪夜夜爽| 三级国产精品欧美在线观看| 国产精品嫩草影院av在线观看| 国产免费福利视频在线观看| 在线免费十八禁| 国产爱豆传媒在线观看| 51国产日韩欧美| 在线天堂最新版资源| 舔av片在线| 成人美女网站在线观看视频| 亚洲av国产av综合av卡| 免费观看性生交大片5| 欧美xxxx性猛交bbbb| 亚洲av成人精品一区久久| 边亲边吃奶的免费视频| 久久6这里有精品| 午夜精品国产一区二区电影| 国产伦在线观看视频一区| 少妇人妻一区二区三区视频| 高清不卡的av网站| av专区在线播放| 26uuu在线亚洲综合色| 国产 精品1| 在线播放无遮挡| 国产精品av视频在线免费观看| 少妇人妻精品综合一区二区| 国产精品一区二区在线观看99| 国产男女内射视频| 亚洲精品乱码久久久v下载方式| 高清视频免费观看一区二区| 天天躁夜夜躁狠狠久久av| 国产男女内射视频| 国产在线一区二区三区精| 在线看a的网站| 日韩强制内射视频| 国产一级毛片在线| 2022亚洲国产成人精品| 女人十人毛片免费观看3o分钟| av视频免费观看在线观看| 国产免费一级a男人的天堂| 一级毛片电影观看| 少妇裸体淫交视频免费看高清| 欧美成人一区二区免费高清观看| 在线观看一区二区三区激情| 纵有疾风起免费观看全集完整版| 在线观看国产h片| 51国产日韩欧美| 久久亚洲国产成人精品v| 国产黄频视频在线观看| 国产av国产精品国产| 国产伦理片在线播放av一区| 久久 成人 亚洲| 人体艺术视频欧美日本| 国产淫片久久久久久久久| 中文字幕av成人在线电影| 麻豆精品久久久久久蜜桃| 亚洲熟女精品中文字幕| 免费看日本二区| 亚洲精品,欧美精品| 亚洲国产色片| 国产精品一区二区三区四区免费观看| 在线观看免费日韩欧美大片 | 国产精品久久久久久精品电影小说 | 精品人妻熟女av久视频| 极品少妇高潮喷水抽搐| 亚洲欧美成人综合另类久久久| 老女人水多毛片| 精品久久久久久久久av| 久久久a久久爽久久v久久| 男女下面进入的视频免费午夜| 插阴视频在线观看视频| 国产久久久一区二区三区| 亚洲国产日韩一区二区| 我的女老师完整版在线观看| 少妇熟女欧美另类| 免费看光身美女| 亚洲国产av新网站| 亚洲精品成人av观看孕妇| 日本黄色日本黄色录像| 国产精品欧美亚洲77777| 日韩电影二区| 七月丁香在线播放| 美女视频免费永久观看网站| 一区二区三区免费毛片| 精品一品国产午夜福利视频| 亚洲真实伦在线观看| 亚洲国产最新在线播放| 国产在线男女| 国产一区亚洲一区在线观看| 久久久久久久久久久丰满| 99热国产这里只有精品6| 啦啦啦啦在线视频资源| 日韩av在线免费看完整版不卡| 91精品国产国语对白视频| 小蜜桃在线观看免费完整版高清| 成人国产av品久久久| 亚洲伊人久久精品综合| 久久久久久人妻| 日韩中文字幕视频在线看片 | 毛片一级片免费看久久久久| 永久网站在线| 国产亚洲一区二区精品| 少妇被粗大猛烈的视频| 国产精品伦人一区二区| 亚洲欧美精品专区久久| 简卡轻食公司| 欧美日韩综合久久久久久| 我要看日韩黄色一级片| 国产日韩欧美在线精品| 亚洲美女视频黄频| 欧美3d第一页| 永久网站在线| 国产免费视频播放在线视频| av视频免费观看在线观看| 高清av免费在线| 一个人免费看片子| 日韩精品有码人妻一区| 最近2019中文字幕mv第一页| 成人毛片60女人毛片免费| 免费观看的影片在线观看| 免费看光身美女| 免费久久久久久久精品成人欧美视频 | 国产精品av视频在线免费观看| 久久久久久久大尺度免费视频| 中文天堂在线官网| 97热精品久久久久久| 青春草亚洲视频在线观看| 中文字幕制服av| 国产精品国产三级专区第一集| 日韩欧美精品免费久久| 日本黄色日本黄色录像| 伦理电影免费视频| 亚洲精华国产精华液的使用体验| 大码成人一级视频| 久久久a久久爽久久v久久| 少妇人妻久久综合中文| 男人添女人高潮全过程视频| 啦啦啦中文免费视频观看日本| 国产av精品麻豆| 高清欧美精品videossex| 欧美三级亚洲精品| 亚洲va在线va天堂va国产| 黄色一级大片看看| 久久久久久久精品精品| 蜜桃在线观看..| 中国三级夫妇交换| 一级av片app| 成人毛片a级毛片在线播放| 高清欧美精品videossex| 人人妻人人澡人人爽人人夜夜| 免费大片18禁| 欧美人与善性xxx| 久久久国产一区二区| 久久久久久久久久久丰满| 国产午夜精品久久久久久一区二区三区| 亚洲精品第二区| 日韩av在线免费看完整版不卡| 人妻少妇偷人精品九色| 亚洲av日韩在线播放| 五月天丁香电影| 国产极品天堂在线| 又爽又黄a免费视频| 激情 狠狠 欧美| 人体艺术视频欧美日本| 狂野欧美白嫩少妇大欣赏| 最新中文字幕久久久久| 最近最新中文字幕大全电影3| 一边亲一边摸免费视频| 久热这里只有精品99| 男女国产视频网站| 成年女人在线观看亚洲视频| 一本—道久久a久久精品蜜桃钙片| 久久久久久久大尺度免费视频| 日日啪夜夜撸| 国产精品一二三区在线看| 国产精品无大码| 亚洲欧美日韩另类电影网站 | 少妇被粗大猛烈的视频| 99热这里只有精品一区| 99久久精品热视频| tube8黄色片| 免费久久久久久久精品成人欧美视频 | 黑人高潮一二区| tube8黄色片| 国产男人的电影天堂91| 少妇精品久久久久久久| 中文精品一卡2卡3卡4更新| 老女人水多毛片| 亚洲av成人精品一二三区| 尤物成人国产欧美一区二区三区| 97超视频在线观看视频| 免费黄频网站在线观看国产| 七月丁香在线播放| 亚洲国产毛片av蜜桃av| 多毛熟女@视频| 亚洲久久久国产精品| 18禁在线无遮挡免费观看视频| 久久久久久久久久人人人人人人| 免费大片18禁| 久久精品国产亚洲av涩爱| 久久国产亚洲av麻豆专区| 99视频精品全部免费 在线| 亚洲av免费高清在线观看| 在线看a的网站| 国产 精品1| 亚洲av国产av综合av卡| 人妻系列 视频| 黄片wwwwww| 精品亚洲乱码少妇综合久久| 国产精品国产三级国产av玫瑰| 亚洲精品一区蜜桃| 大码成人一级视频| 80岁老熟妇乱子伦牲交| 偷拍熟女少妇极品色| 色综合色国产| 热99国产精品久久久久久7| 99热国产这里只有精品6| 蜜桃在线观看..| 国产精品熟女久久久久浪| 国产精品福利在线免费观看| 亚洲精品乱码久久久v下载方式| 国产欧美另类精品又又久久亚洲欧美| 一级毛片aaaaaa免费看小| 国内揄拍国产精品人妻在线| 亚洲图色成人| 97超视频在线观看视频| 美女中出高潮动态图| 国产精品欧美亚洲77777| 91精品一卡2卡3卡4卡| 日本爱情动作片www.在线观看| 国产高清不卡午夜福利| 亚洲国产色片| 久久99蜜桃精品久久| 蜜桃亚洲精品一区二区三区| 最黄视频免费看| 国产大屁股一区二区在线视频| 国产精品国产三级国产av玫瑰| 又黄又爽又刺激的免费视频.| av一本久久久久| 亚洲在久久综合| 国精品久久久久久国模美|