亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        新型組合式垂直軸風(fēng)機(jī)的優(yōu)化設(shè)計(jì)與性能研究

        2015-02-19 02:41:49霍玉雷楊洪偉談宏飛

        徐 璋,霍玉雷,陳 勇,楊洪偉,談宏飛

        (浙江工業(yè)大學(xué) 能源與動(dòng)力工程研究所,浙江 杭州 310014)

        新型組合式垂直軸風(fēng)機(jī)的優(yōu)化設(shè)計(jì)與性能研究

        徐璋,霍玉雷,陳勇,楊洪偉,談宏飛

        (浙江工業(yè)大學(xué) 能源與動(dòng)力工程研究所,浙江 杭州 310014)

        摘要:對(duì)于垂直軸風(fēng)機(jī),升力型風(fēng)能轉(zhuǎn)化效率高,但是其在低風(fēng)速下不能自啟動(dòng);阻力型可以在低風(fēng)速下自啟動(dòng),但其轉(zhuǎn)速低,而且其風(fēng)能利用率低.為了充分發(fā)揮二者性能,我們?cè)O(shè)計(jì)制作了一種新型組合式垂直軸風(fēng)力機(jī),通過采用超越離合器作為連接件將Darrieus風(fēng)力機(jī)和Savonius風(fēng)力機(jī)組合而成,并對(duì)其進(jìn)行了三維數(shù)值模擬和吹風(fēng)實(shí)驗(yàn).通過對(duì)比有無加裝S型風(fēng)機(jī)時(shí)風(fēng)機(jī)啟動(dòng)風(fēng)速和輸出功率的模擬結(jié)果與實(shí)驗(yàn)數(shù)據(jù),發(fā)現(xiàn)該新型組合式垂直軸風(fēng)機(jī)啟動(dòng)風(fēng)速明顯降低,風(fēng)能利用率僅略有下降,很好地實(shí)現(xiàn)了互補(bǔ).

        關(guān)鍵詞:CFD;垂直軸風(fēng)力機(jī);升力型;阻力型;超越離合器;吹風(fēng)實(shí)驗(yàn)

        Optimum design and study on the properties of a new combined

        type vertical axis wind turbine

        XU Zhang, HUO Yulei, CHEN Yong, YANG Hongwei, TAN Hongfei

        (Institute of Energy and Power Engineering, Zhejiang University of Technology, Hangzhou 310014, China)

        Abstract:For vertical axis wind turbines(VAWT),the lift type has a high conversion efficiency but it cannot start itself at low wind speed;the resistance type can start itself at low wind speed, but it also has a low fan speed and a low wind power utilization coefficient. In order to fully take advantage of their merits, a new-type VAWT, which joined a Savonius wind turbine and a Darrieus wind turbine by an overrunning clutch, was proposed. 3-d numerical simulation and experimental wind tunnel tests were carried out.By comparing the power output and start-up wind speed between the units with and without S-type wind turbine,it shows that the start-up wind speed of this new type VAWT decreases significantly, while the wind power utilization coefficient only reduces slightly, which combines the advantages of two types of wind turbines.

        Keywords:CFD; VAWT; lift type; resistance type; overrunning clutch; wind tunnel tests

        隨著能源與環(huán)境之間的矛盾日益加劇,作為可再生能源的風(fēng)能,已受到世界各國(guó)的歡迎和重視[1].但風(fēng)力發(fā)電在實(shí)際應(yīng)用推廣中遇到了入網(wǎng)難、棄風(fēng)率高、性能不穩(wěn)定等諸多問題[2],正因如此,分布式能源作為一種解決方案應(yīng)運(yùn)而生,各種新型的風(fēng)力發(fā)電機(jī),尤其是垂直軸風(fēng)機(jī)層出不窮.

        分別作為垂直軸風(fēng)力機(jī)升力型和阻力型最典型代表形式的薩窩紐斯型(Savonius type)和達(dá)里厄型(Darrieus type)風(fēng)力發(fā)電機(jī)組,有著各自的優(yōu)缺點(diǎn).為使Darrieus風(fēng)力機(jī)能在低風(fēng)速下啟動(dòng),采用加裝Savonius風(fēng)力機(jī)的方法[3].當(dāng)把S型風(fēng)力機(jī)串加在D型風(fēng)力機(jī)上后,D型風(fēng)力機(jī)由S型風(fēng)力機(jī)低速轉(zhuǎn)動(dòng)下產(chǎn)生的較大力矩帶動(dòng)從而進(jìn)行旋轉(zhuǎn),啟動(dòng)后,風(fēng)能利用率高的D型風(fēng)機(jī)起主要風(fēng)能轉(zhuǎn)化作用,保證其整體性能.根據(jù)以上原理,本研究設(shè)計(jì)出Darrieus-Savonius組合式垂直軸風(fēng)力機(jī)模型,以實(shí)現(xiàn)S型葉片和D型葉片的優(yōu)勢(shì)互補(bǔ).同時(shí),針對(duì)此類組合式垂直軸風(fēng)力機(jī)傳統(tǒng)設(shè)計(jì)中存在的高速運(yùn)行時(shí)S型與D型相互影響的問題,進(jìn)行了新型組合式垂直軸風(fēng)機(jī)的優(yōu)化設(shè)計(jì),并對(duì)其氣動(dòng)性能和風(fēng)能轉(zhuǎn)換效率進(jìn)行數(shù)值模擬與吹風(fēng)實(shí)驗(yàn),目的是研究及驗(yàn)證其可行性.

        1新型組合式垂直軸風(fēng)力機(jī)的優(yōu)化和結(jié)構(gòu)設(shè)計(jì)

        查閱了相關(guān)的一些文獻(xiàn)[4-8]后,通過進(jìn)行優(yōu)化和尺寸計(jì)算,確定了此組合式垂直軸風(fēng)力機(jī)中的的具體設(shè)計(jì)參數(shù),分別如表1,2所示.

        表1 S型風(fēng)機(jī)具體設(shè)計(jì)參數(shù)

        表2 D型風(fēng)機(jī)具體設(shè)計(jì)參數(shù)

        S型風(fēng)輪制作的主要材料為白鐵皮,為了最大化采集利用各個(gè)方向上的風(fēng)能,設(shè)計(jì)制作兩個(gè)相同規(guī)格的半圓形S型風(fēng)輪以上下90°方位差疊加布置.D型風(fēng)輪的葉片選擇玻璃鋼材質(zhì),形式為結(jié)構(gòu)簡(jiǎn)單,風(fēng)阻較小的NACA0015對(duì)稱翼型,葉片和支撐桿通過螺母螺栓與磁懸浮發(fā)電機(jī)頂蓋相連.S型風(fēng)機(jī)整體安裝于D型風(fēng)機(jī)之上,防止S型與D型風(fēng)力機(jī)之間迎風(fēng)面積相互影響.在此設(shè)計(jì)條件下,S型風(fēng)機(jī)在低風(fēng)速下便可以啟動(dòng),繼而帶動(dòng)D型風(fēng)力機(jī)工作.

        另外,組合式風(fēng)力機(jī)平穩(wěn)運(yùn)行之后,在氣動(dòng)性能方面,低速旋轉(zhuǎn)的S型會(huì)對(duì)高速旋轉(zhuǎn)的D型產(chǎn)生較大的消極影響.本設(shè)計(jì)擬采用帶軸承的CKZ型號(hào)的超越離合器,主動(dòng)軸與S型風(fēng)力機(jī)聯(lián)接,從動(dòng)軸與D型風(fēng)力機(jī)聯(lián)接.最終使得上下兩部分風(fēng)機(jī)在D型風(fēng)力機(jī)被S型風(fēng)力機(jī)順利帶動(dòng)后及時(shí)斷開,既保證了D型風(fēng)力機(jī)的啟動(dòng),又不使得D型的氣動(dòng)性能受S型的影響.安裝好的風(fēng)力機(jī)具體形式如圖1所示.

        圖1 組合式垂直軸風(fēng)力機(jī)Fig.1 Combined type vertical axis wind turbine

        2新型組合式垂直軸風(fēng)力機(jī)的數(shù)值模擬

        針對(duì)前述設(shè)計(jì)制作的風(fēng)力機(jī)模型,在有無S型風(fēng)機(jī)的情況下,分別建立模型,確定計(jì)算區(qū)域,在給定工況下,利用ANSYS Fluent軟件對(duì)模型進(jìn)行三維非定常數(shù)值模擬,得到風(fēng)機(jī)的模擬輸出轉(zhuǎn)矩和輸出功率,以及風(fēng)能利用系數(shù)等氣動(dòng)參數(shù),從而對(duì)比得出結(jié)論.

        2.1網(wǎng)格劃分與邊界條件的確定

        風(fēng)輪轉(zhuǎn)動(dòng)過程的非定常數(shù)值模擬均采用滑移網(wǎng)格模型,有S型風(fēng)機(jī)的總體網(wǎng)格數(shù)240萬,無S型風(fēng)機(jī)的三維網(wǎng)格數(shù)為110萬.出于忽略體積力和外部熱源的考慮,在模擬過程中均以非定常不可壓縮Navier-Stokes方程作為控制方程.然后依次設(shè)定如下參數(shù):

        均勻來流風(fēng)速為V0=8 m/s,D型、S型風(fēng)輪規(guī)格參考表1,2;設(shè)置速度入口邊界和壓力出口邊界;有S型風(fēng)機(jī)時(shí)選取壁面在兩部分風(fēng)機(jī)葉片以及封蓋上,無S型風(fēng)機(jī)時(shí)在D型風(fēng)機(jī)葉片上,且均設(shè)置無滑移壁面邊界條件;在外流場(chǎng)與兩類型風(fēng)力機(jī)各自的旋轉(zhuǎn)區(qū)域的交界面處,均設(shè)置Interface邊界條件;由于整機(jī)的三維數(shù)值模擬應(yīng)用滑移網(wǎng)格技術(shù),因此將葉片旋轉(zhuǎn)區(qū)域及其附近區(qū)域分別設(shè)定為不同的流體區(qū)域:D型葉片區(qū)域和S型葉片區(qū)域的旋轉(zhuǎn)速度分別為180 r/min和120 r/min;湍流模型采用標(biāo)準(zhǔn)k—ε湍流模型;計(jì)算方法采用SIMPLE算法,離散方法均采用二階迎風(fēng)格式[9].

        2.2三維CFD非定常數(shù)值模擬結(jié)果

        2.2.1風(fēng)力機(jī)的轉(zhuǎn)矩模擬

        根據(jù)模擬結(jié)果,分析其氣動(dòng)性能[10]可得:無S型風(fēng)機(jī)的情況下,平均轉(zhuǎn)矩系數(shù)大約為0.2;在180 r/min的轉(zhuǎn)速下,葉片受到周期性轉(zhuǎn)矩矩范圍內(nèi)葉片的平均轉(zhuǎn)矩為6.2 Nm.而有S型風(fēng)機(jī)的情況下,平均轉(zhuǎn)矩系數(shù)在0.3左右.同理可以得到,在D型180 r/min和S型120 r/min的條件下,葉片的平均轉(zhuǎn)矩為9.5 Nm.對(duì)比可以發(fā)現(xiàn),加裝S型以后風(fēng)力機(jī)輸出的轉(zhuǎn)矩系數(shù)均高于相同條件下的無S型,說明加裝S型的風(fēng)機(jī)在低風(fēng)速下比沒有加裝S型的更加容易啟動(dòng).

        2.2.2風(fēng)力機(jī)的風(fēng)能利用率模擬

        通過風(fēng)力機(jī)葉輪旋轉(zhuǎn)面的風(fēng)能,并不能全部被葉輪吸收利用,其利用量與總風(fēng)能量的比值,即風(fēng)能利用率,代表了風(fēng)力機(jī)將風(fēng)能轉(zhuǎn)化為機(jī)械能的能力,是在忽略轉(zhuǎn)換裝置效率的情況下,評(píng)價(jià)風(fēng)力機(jī)整個(gè)葉輪系統(tǒng)其氣動(dòng)性能最為直觀的指標(biāo),且一臺(tái)風(fēng)力機(jī)的風(fēng)能利用率隨尖速比變化,存在一個(gè)最佳尖速比.在不同的尖速比下,對(duì)有無S型風(fēng)機(jī)的情分別進(jìn)行三維非穩(wěn)態(tài)數(shù)值模擬,計(jì)算過程中對(duì)葉片所受到的轉(zhuǎn)矩系數(shù)進(jìn)行監(jiān)測(cè),依據(jù)風(fēng)能利用系數(shù)與轉(zhuǎn)矩系數(shù)的計(jì)算關(guān)系,對(duì)模擬結(jié)果進(jìn)行處理,可以得到有無S型風(fēng)機(jī)的情況下整機(jī)風(fēng)能利用系數(shù)變化關(guān)系曲線的對(duì)比結(jié)果,如圖2所示.

        圖2 有無S型風(fēng)機(jī)情況下整機(jī)風(fēng)能利用系數(shù)變化曲線Fig.2 Curve of wind power utilization coefficient of the wind turbine with and without S-type

        由圖2可見:有無S型情況下的風(fēng)能利用系數(shù)隨尖速比的變化趨勢(shì)基本一致,且均在尖速比3左右時(shí)風(fēng)能利用系數(shù)達(dá)到最大.加裝S型后整機(jī)的風(fēng)能利用系數(shù)均高于相應(yīng)尖速比下的無S型的情況,加裝S型的最高風(fēng)能利用率為0.27,平均風(fēng)能利用率為0.25;無S型的最佳風(fēng)能利用系數(shù)為0.29,平均風(fēng)能利用系數(shù)為0.25.從模擬的結(jié)果中可以看出:有S型的情況下風(fēng)能利用系數(shù)低于無S型的情況.其差異的主要原因是組合式垂直軸風(fēng)力機(jī)中加裝的兩個(gè)S型風(fēng)力機(jī),其屬于阻力型的風(fēng)機(jī),風(fēng)能利用系數(shù)較D型風(fēng)力機(jī)有一定程度的差距;且由于超越離合器本身主動(dòng)軸和從動(dòng)軸之間的相對(duì)轉(zhuǎn)動(dòng)存在較低程度的摩擦,風(fēng)力機(jī)的性能會(huì)受到一定影響,導(dǎo)致加裝S型后風(fēng)能利用系數(shù)有所降低.同樣是由于超越離合器的存在,在D型工作后使S型和D型斷開,各自旋轉(zhuǎn),D型起主要的風(fēng)能轉(zhuǎn)化作用,從而使得風(fēng)能利用系數(shù)只是略有降低,并未造成實(shí)質(zhì)性的影響.同時(shí),根據(jù)上節(jié)的結(jié)果來看,該組合式風(fēng)機(jī)能在低風(fēng)速下獲得更好的啟動(dòng)性能,所以總體來說利大于弊,符合實(shí)驗(yàn)?zāi)康?

        3吹風(fēng)實(shí)驗(yàn)及其結(jié)果分析

        在實(shí)際應(yīng)用中,由于風(fēng)力機(jī)的組裝、制造過程中的誤差以及工作時(shí)的各部件摩擦損耗,會(huì)使得風(fēng)力機(jī)的實(shí)際氣動(dòng)性能較模擬結(jié)果有一定程度的下降.為了對(duì)該組合式風(fēng)力機(jī)的實(shí)際氣動(dòng)性能進(jìn)行研究驗(yàn)證,制造1∶1模型并進(jìn)行了吹風(fēng)實(shí)驗(yàn):在不同風(fēng)速下,利用控制器遠(yuǎn)程管理軟件實(shí)時(shí)測(cè)量和記錄特定風(fēng)速下風(fēng)力機(jī)的轉(zhuǎn)速和輸出的電壓、電流、功率.為了能夠更加直觀地對(duì)比有無S型情況下風(fēng)機(jī)的性能,對(duì)測(cè)量數(shù)據(jù)進(jìn)行可視化處理,得到了有無S型情況下各參數(shù)之間的變化關(guān)系曲線圖.圖3為兩種情況下風(fēng)速與轉(zhuǎn)速變化關(guān)系曲線,圖4為兩種情況下轉(zhuǎn)速與輸出功率的變化關(guān)系曲線.

        圖3 有無S型風(fēng)機(jī)情況下風(fēng)速與轉(zhuǎn)速變化關(guān)系曲線圖Fig.3 Curve of wind speed and the fan speed of the wind turbine with and without S-type

        圖4 有無S型風(fēng)機(jī)情況下轉(zhuǎn)數(shù)與輸出功率變化關(guān)系曲線圖Fig.4 Curve of fan speed and Output power of the wind turbine with and without S-type

        對(duì)比圖3的兩條曲線并結(jié)合測(cè)量結(jié)果可以看出:有S型的情況下風(fēng)力機(jī)的啟動(dòng)風(fēng)速為0.782 m/s,遠(yuǎn)低于無S的情況下風(fēng)力機(jī)的啟動(dòng)速度1.768 m/s.在風(fēng)速達(dá)到2.5 m/s前,有S型的情況下轉(zhuǎn)速增加較為緩慢,這是因?yàn)镾型轉(zhuǎn)速大于D型轉(zhuǎn)速,兩者尚未脫離,S型風(fēng)機(jī)轉(zhuǎn)速的上升速度較D型慢,會(huì)對(duì)整機(jī)的造成一定的阻礙.但是在風(fēng)速達(dá)到2.5 m/s后,D型轉(zhuǎn)速開始大于S型,由于超越離合器的存在,S型脫離,變化關(guān)系曲線趨于一致.

        對(duì)比圖4的兩條曲線可以看出:兩種情況下風(fēng)機(jī)的輸出功率隨轉(zhuǎn)速具有相似的變化,只是在S型風(fēng)力機(jī)經(jīng)超越離合器脫離前,有S型的情況在數(shù)值上略微低于無S型的情況,這種差距在脫離后基本消失.

        根據(jù)以上的測(cè)量結(jié)果,可進(jìn)一步得到轉(zhuǎn)數(shù)與風(fēng)能利用系數(shù)的關(guān)系,如圖5所示.

        圖5 有無S型風(fēng)機(jī)情況下轉(zhuǎn)數(shù)與風(fēng)能利用率變化關(guān)系曲線圖Fig.5 Curve of fan speed and wind power utilization coefficient of the wind turbine with and without S-type

        由圖5可見:在200 r/min的轉(zhuǎn)速以內(nèi),加裝S型風(fēng)機(jī)前后,風(fēng)機(jī)的風(fēng)能利用系數(shù)伴隨著轉(zhuǎn)速具有相似的變化規(guī)律,臨界轉(zhuǎn)速都在160 r/min附近,超過160 r/min后,風(fēng)能利用系數(shù)逐漸降低.在數(shù)值上無S型的風(fēng)機(jī)整體偏低,有S型的情況下最高風(fēng)能利用率在0.25左右,高于無S型情況下的0.27,且都略低于模擬得到的最高風(fēng)能利用系數(shù).由于模擬時(shí)的三維模型經(jīng)過簡(jiǎn)化,忽略了中心轉(zhuǎn)軸及支撐桿的風(fēng)阻,因此其結(jié)果與預(yù)期是一致的.

        4結(jié)論

        對(duì)設(shè)計(jì)的新型組合式垂直軸風(fēng)力機(jī)進(jìn)行的CFD數(shù)值模擬實(shí)驗(yàn)表明了加裝Savonius風(fēng)力機(jī)后,組合式風(fēng)力機(jī)整體的啟動(dòng)風(fēng)速相對(duì)于單獨(dú)的Darrieus風(fēng)力機(jī)確有很大減小,雖然與此同時(shí),風(fēng)機(jī)的輸出功率,風(fēng)能利用率略有降低,但總體來看,利大于弊.吹風(fēng)實(shí)驗(yàn)由于考慮了各處的實(shí)際風(fēng)阻,同樣轉(zhuǎn)數(shù)下測(cè)得的風(fēng)能利用率比模擬的要小,啟動(dòng)風(fēng)速要比模擬的大,但其趨勢(shì)和模擬結(jié)果相同,很好地驗(yàn)證了模擬得出的結(jié)論.這一切都得益于創(chuàng)新性地提出使用超越離合器作為連接件,把影響降到了最低.但大多數(shù)部件經(jīng)由手工制作,還達(dá)不到實(shí)際應(yīng)用的水平,在性能上還有提升的空間,相信會(huì)得到更加廣闊的應(yīng)用和發(fā)展.

        參考文獻(xiàn):

        [1]葉瑞克,王學(xué)謙,鮑健強(qiáng).基于產(chǎn)業(yè)結(jié)構(gòu)轉(zhuǎn)型升級(jí)的浙江省節(jié)能減碳研究[J].浙江工業(yè)大學(xué)學(xué)報(bào):社會(huì)科學(xué)版,2011,10(3):283-289.

        [2]張運(yùn)洲,胡泊.“三北”地區(qū)風(fēng)電開發(fā)、輸送及消納研究[J].中國(guó)電力,2012,45(9):1-6.

        [3]何宗敬.Darrieus-Savonius組合風(fēng)輪氣動(dòng)性能的研究[J].太陽能學(xué)報(bào),1992,13(3):239-243.

        [4]趙振宙,鄭源,周大慶,等.基于數(shù)值模擬Savonius風(fēng)力機(jī)性能優(yōu)化研究[J].太陽能學(xué)報(bào),2010,31(7):907-911.

        [5]廖康平.垂直軸風(fēng)機(jī)葉輪空氣動(dòng)力學(xué)性能研究[D].哈爾濱:哈爾濱工程大學(xué)船舶工程學(xué)院,2006.

        [6]張國(guó)臣,徐志暉,張巖,等.一種復(fù)合型垂直軸風(fēng)力機(jī)實(shí)驗(yàn)研究[J].沈陽航空工業(yè)學(xué)院學(xué)報(bào),2010,27(5):29-33.

        [7]李健,孟憲舉.小型H型垂直軸風(fēng)機(jī)安裝角度的CFD數(shù)值模擬與分析[J].河北理工大學(xué)學(xué)報(bào),2011,33(4):56-59.

        [8]何宗敬.Darrieus-Savonius組合風(fēng)輪幾何尺寸確定方法探討[J].太陽能學(xué)報(bào),1993,14(4):325-329.

        [9]潘柏松,陳棟棟,韓斌,等.基于不同厚度兆瓦級(jí)風(fēng)力機(jī)葉片動(dòng)態(tài)特性研究[J].浙江工業(yè)大學(xué)學(xué)報(bào),2013,41(4):364-368.

        [10]張兆鑫,趙元虎,徐奔馳.基于CFD法的小型風(fēng)機(jī)非扭曲葉片氣動(dòng)性能分析[J].浙江工業(yè)大學(xué)學(xué)報(bào),2013,41(1):68-72.

        (責(zé)任編輯:劉巖)

        中圖分類號(hào):TK83

        文獻(xiàn)標(biāo)志碼:A

        文章編號(hào):1006-4303(2015)03-0261-04

        作者簡(jiǎn)介:徐璋(1973—),男,江蘇海門人,副教授,研究方向?yàn)槊喝紵廴疚锱欧趴刂萍靶履茉撮_發(fā),E-mail:xzzyn@zjut.edu.cn.

        收稿日期:2015-01-12

        中文字幕第一页人妻丝袜| 日日摸夜夜添狠狠添欧美| 国产精品成人午夜久久| 亚洲国产一区久久yourpan| 视频女同久久久一区二区| 天天做天天摸天天爽天天爱 | 亚洲精品无播放器在线播放| 国产91福利在线精品剧情尤物| 日韩精品一区二区av在线| 免费国产一区二区视频| 成人午夜福利视频后入| 在教室伦流澡到高潮h麻豆| 亚洲黄色性生活一级片| 亚洲国产大胸一区二区三区| 中文字幕亚洲精品一区二区三区| 痉挛高潮喷水av无码免费| 亚洲夜夜骑| 精品亚洲视频免费观看网站| 国产精品亚洲二区在线看| 久久人妻少妇嫩草av无码专区| 欧美性受xxxx黑人xyx性爽| 中文字幕亚洲精品码专区| 北条麻妃在线中文字幕| 国产高跟黑色丝袜在线| 国产在线丝袜精品一区免费| 免费人成黄页网站在线观看国产 | 嫩草伊人久久精品少妇av| 国产97色在线 | 日韩| 国产成人亚洲综合无码DVD| 国产成人久久精品二区三区| 日韩精品综合一本久道在线视频| 国产av人人夜夜澡人人爽麻豆| 91情侣视频| 国产一区二区三区在线av| 少妇被又大又粗又爽毛片久久黑人 | 婷婷成人亚洲综合国产| 精品少妇一区二区三区免费| 无人视频在线观看免费播放影院| 亚洲熟妇网| 日韩在线精品免费观看| 国产av一区二区三区天堂综合网|