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

    船舶艉流場受預旋定子影響試驗研究

    2021-11-05 05:39:42郭春雨徐鵬韓陽王超郭欣雨
    哈爾濱工程大學學報 2021年10期
    關鍵詞:盤面螺旋槳流場

    郭春雨, 徐鵬, 韓陽, 王超, 郭欣雨

    (1.哈爾濱工程大學 船舶工程學院, 黑龍江 哈爾濱 150001; 2.水下測控技術重點實驗室,遼寧 大連 116013)

    在經(jīng)濟全球化的今天,航運作為貨物運輸?shù)囊粋€重要渠道,對拉動經(jīng)濟發(fā)展具有重要作用。然而,燃油成本的日趨提高,增加了企業(yè)負擔。鑒于此,航運企業(yè)必須大力倡導節(jié)能技術,以達到降本增效的目的。另一方面,全球變暖趨勢日益增加,船舶能效指數(shù)實施勢在必行,船舶減排的重要性也受到了廣泛關注。水動力節(jié)能附體效率高、成本低、易于實施等優(yōu)勢逐漸凸顯[1]。預旋定子是水動力節(jié)能附體中的一種,其由多個固定不動且沿周向分布的機翼型葉片構成,它能夠在槳前產(chǎn)生預旋流,降低螺旋槳能量損失,從而提高船舶的推進性能[2]。Celik等[3]以升力線理論為基礎設計了一個螺旋槳后置定子程序并對螺旋槳的水動力性能進行分析。黃樹權等[4]通過結合數(shù)值模擬和模型試驗2種方法研究了預旋定子設計參數(shù)對水動力性能的影響。苗飛等[5]利用數(shù)值模擬的方法研究了定子周向布置形式、單個定子周向影響范圍以及來流攻角對定子預旋作用的影響,并指出2個定子最佳布置角度為45°。苗飛等[6]通過分析前置預旋定子葉剖面阻力特性提出一種新型的葉剖面設計方法,最終船模試驗證實了方案可行性。楊帆等[7]采用數(shù)值模擬和模型試驗的方法評估了安裝有預旋定子的某大型散貨船節(jié)能效果,并通過速度矢量圖和伴流等值圖闡述節(jié)能機理。楊帆等[8]在先前研究基礎上設計出一種預旋定子,并通過CFD方法對影響定子性能的3個主要參數(shù)(襟翼角度、襟翼寬度、定子直徑)進行分析。凌乃俊等[9]以升力線理論為基礎,結合CFD數(shù)值方法,通過引進加權因子,預報了設計環(huán)量對船舶節(jié)能效果的影響,并得出最佳設計環(huán)量,并與模型試驗對比驗證了設計方法的可行性。Gao等[10]對預旋定子的相關參數(shù)(周向布置形式、葉片數(shù)量、葉片攻角)對節(jié)能效果影響進行研究,并設計出一種預旋定子。

    進入21世紀,國內外多家艦船研究機構都對前置預旋定子進行了較深入的研究,如中國船舶科學研究中心、大宇造船海洋株式會社、三菱重工業(yè)有限公司、荷蘭海事研究所等。德國漢堡水池、瑞典船舶研究中心以及韓國船舶海洋工程研究所在船模水池進行過一系列模型試驗,表明預旋定子的節(jié)能效果可達3%~5%,而且目前已在實船上獲得應用[9],并取得了較好的節(jié)能收益。在國內,對預旋定子的研究起步較晚,整體研究水平較國外有一定差距,且研究大多集中在數(shù)值模擬,對模型試驗驗證較少。

    通過模型試驗獲得準確的船舶艉流場信息一直是船舶與海洋工程實驗流體力學領域研究的熱點[11]。由于船艉幾何形狀的特點,艉流場流動復雜,具有較強的三維流動分離現(xiàn)象,這對艉流場信息的獲取帶來了很大的挑戰(zhàn)。如今,粒子圖像測速技術(particle image velocimetry,PIV)作為一種瞬態(tài)、全局、無接觸的測量方法備受流體力學屆關注。PIV技術從20世紀80年代開始發(fā)展至今,相關技術越來越成熟、越來越完善,同時準確性得到了眾多學者的驗證,試驗結果也得到了廣泛認可。因此,PIV技術迅速成為試驗流體力學精細流場測量的首要方法。為此,本研究基于哈爾濱工程大學船模拖曳水池水下體視粒子圖像測速系統(tǒng)(stereoscopic particle image velocimetry,SPIV)系統(tǒng),對船舶艉流場進行測量,通過對有無預旋定子的伴流場中三向速度分布、渦量場分布以及旋渦強度進行艉流場流動特性分析。

    1 試驗設備、試驗模型與測量工況

    1.1 基本試驗設備與條件

    本文依托哈爾濱工程大學船模拖曳水池開展某散貨船艉流場PIV試驗,試驗時遵循國際拖曳水池會議(international towing tank conference,ITTC)阻力測量規(guī)則,滿足Fr相似,雷諾數(shù)Re>2×106?;緦嶒炘O備如下:

    拖曳水池:長度108 m;寬度7 m;水深3.5 m;

    拖車:車速V≤6.5 m/s;穩(wěn)速范圍:0.1~6.5 m/s;精度:0.1%;平均加速度:a+>0.09g;平均減速度:a->0.15g。

    通過4自由度適航儀連接船模與拖車,4自由度適航儀型號為GEL-421-1,精度為1%。

    1.2 試驗模型與測量工況

    本文進行PIV流場測量試驗的模型為某型號單槳散貨船,實船總長度為324.9 m,船??傞LLOA=5.415 m,縮尺比λ=60,主要尺度如表1所示。

    表1 某型號單槳散貨船模型主要參數(shù)Table 1 The main parameters of ship model

    散貨船模型船艉部以及預旋定子模型均進行噴涂處理,防止試驗時激光強反光對試驗結果影響。試驗模型及船體坐標系如圖1所示,現(xiàn)規(guī)定原點位于螺旋槳盤面中心位置,X軸沿船模中縱剖面指向船艉,Z軸垂直向上,Y軸遵循右手法則指向右舷。預旋定子由上海船舶運輸科學研究院提供,為不對稱型葉片,葉片左右分布不對稱,左側有3個葉片,每個葉片剖面不對稱,根據(jù)圖紙采用 PLA材料經(jīng)3D打印制作而成,并對預旋定子模型表面進行光滑處理。預旋定子葉片相對位置以及實物如圖2所示。試驗工況如表2所示,測量截面為螺旋槳盤面處(x=0 cm)及槳盤面后7.5 cm處(x=7.5 cm)。

    表2 實驗工況Table 2 Test conditions

    圖1 試驗模型及船體坐標系Fig.1 Test model and coordinate system

    圖2 預旋定子葉片相對位置以及實物Fig.2 Pre-swirl stator and relative position of the blade of the stator

    2 船舶艉流場PIV試驗測量

    2.1 DANTEC隨車式水下SPIV系統(tǒng)

    本實驗所用船舶艉流場測量設備為DANTEC公司為哈爾濱工程大學船模拖曳水池定制的隨車式水下SPIV系統(tǒng),布置于拖車右側,SPIV系統(tǒng)規(guī)格參數(shù)如下:

    CCD相機分辨率:2 048 pixel×2 048 pixel;

    雙脈沖Nd:YAG激光器:最大脈沖頻率15 Hz;最大脈沖激光能量:1 200 mJ;常規(guī)工作能量:200 mJ;

    激光束持續(xù)時長:4 ns;激光波長:532 nm;脈沖激光片厚度:0.6 mm;

    最大測量范圍:400 mm×400 mm;

    示蹤粒子:聚酰胺顆粒。

    該套水下SPIV系統(tǒng)包含兩臺CCD相機、激光器、同步器和控制器。CCD相機以及激光光學元件布置于水下雷體中,激光導光臂和相機線纜鋪設在支柱1和支柱2中,激光器、同步器以及控制器安裝在拖車上。SPIV測試系統(tǒng)及拆解示意圖如圖3所示。

    圖3 SPIV測試系統(tǒng)Fig.3 The SPIV system of the probe in the submersible PIV system

    2.2 PIV流場數(shù)據(jù)采集、測量與分析

    采用50 μm的聚酰胺顆粒作為示蹤粒子,每次試驗開始前通過示蹤粒子播撒裝置將示蹤粒子撒入待測水體,使每次測量過程中平均粒子濃度約為150 g/m3,即每個查問區(qū)域具有15~20個粒子,這樣能夠保證互相關分析的精確度。

    系統(tǒng)工作時,由同步器控制激光器和2臺CCD相機,在發(fā)出激光的瞬間2臺相機從不同的角度拍攝測量區(qū)域中被激光照亮的示蹤粒子,得到2組不同角度的平面二維流場矢量結果,通過標定信息最終重構出該平面測量區(qū)域的三維速度矢量。SPIV艉流場測量圖與測量示意圖如圖4所示。

    圖4 SPIV艉流場測量實物與測量建模Fig.4 Measurement diagram and schematic of the SPIV measurement of the wake field

    相機與脈沖激光器通過同步器控制,將SPIV系統(tǒng)的采集速率調至最大即7.5 Hz,保證進行測量時達到最大的樣本采集量,每次測量過程中一共記錄250組圖像,拍攝總時長為33.33 s,數(shù)據(jù)分析時依據(jù)此250組瞬時照片進行時間平均得到最終結果。

    CCD相機為雙幀拍攝模式,2幀之間間隔為400 μs,這樣能夠保證示蹤粒子在單位時間間隔內運動的最大位移不超過半個查詢區(qū)域。采用Dantec Studio 6.2軟件對試驗圖片進行批量處理與后續(xù)分析。自適應互相關算法獲得流場矢量信息,查詢區(qū)域大小設置為128 pixel×128 pixel,64 pixel×64 pixel,32 pixel×32 pixel,重疊率為50%。

    吳鐵成[12]應用該SPIV系統(tǒng)在拖曳水池中進行了大量艉流場試驗,對該系統(tǒng)不確定度進行分析,并對KCS標模多工況下伴流場進行測量,與公開數(shù)據(jù)進行對比,驗證了該系統(tǒng)的準確性。

    3 船舶伴流場PIV試驗結果及對比分析

    3.1 船舶艉部伴流場速度分布

    圖5為某散貨船縮比模型無預旋定子工況槳盤面處伴流場測試結果,利用拖曳航速U對結果進行無量綱化處理。由于激光反光及槳軸遮擋,部分流場結果無法獲取,依據(jù)船舶的對稱性,僅給出一半的流場測試結果,從圖5中可以看出,散貨船伴流場由于受到船舶艉部形狀的影響,會產(chǎn)生舭渦,受舭渦影響,速度云圖等值線具有明顯的“鉤狀”結構,對應0.25≤u/U≤0.35,軸向速度等值線之所以會產(chǎn)生“鉤狀”畸變,是由于舭渦的存在將船體伴流中動能較低的流體傳遞至船體中心附近,弱化了該處原有速度,形成局部“鉤狀”結構。與此同時還觀察到,一個與舭渦旋向相反、位于螺旋槳槳軸下方的假轂轂帽渦,上述現(xiàn)象與文獻[13]中描述類似。

    圖5 無預旋定子工況槳盤面處艉流場Fig.5 The wake field without pre-swirl stator at propeller plane

    圖6為存在預旋定子工況下槳盤面處伴流場測試結果,利用拖曳航速U對結果進行無量綱化處理。槳盤面上方存在大量脈動速度是由于粒子濃度不均造成的。通過圖6可以看出,預旋定子對艉流具有干擾作用,主要影響區(qū)域位于有葉片存在的左半部分,上方“鉤狀”結構被破壞,表現(xiàn)出速度等值線向槳盤面外部偏移。對于無葉片存在的右半部分,流場結果大致與圖5相似,速度等值線具有明顯的“鉤狀”結構,流場結構具有一致性。但兩者仍有差別,表現(xiàn)出無葉片一側的軸向速度增加,這是由于安裝預旋定子以后,在槳盤面處,定子對水流有匯聚作用[14],使更多的水流匯聚到槳盤面,間接使得無定子一側區(qū)域的平均軸向速度增加??梢缘贸?,預旋定子不僅僅影響葉片存在的一側,對沒有葉片的一側也會產(chǎn)生影響。同樣,預旋定子對展向、垂向速度都有較大影響,使得流場等值線向槳盤面中心收縮,增加了槳盤面左上方流場不均勻程度。通過流線圖可以看出,預旋定子產(chǎn)生的預旋流改變了槳軸下方的流線,降低了假轂轂帽渦強度。

    圖6 有預旋定子工況槳盤面處艉流場Fig.6 The wake field with pre-swirl stator at propeller plane

    通過矢量坐標轉換,得到了槳盤面處艉流場周向速度,如圖7所示。利用拖曳航速U對結果進行無量綱化處理,并給出速度矢量大小。通過圖7可以看出,預旋定子能夠顯著增加槳盤面外圍的周向速度,對靠近槳盤面中心區(qū)域影響較小。對比有無預旋定子周向速度發(fā)現(xiàn),有預旋定子工況無葉片一側周向速度結構與無預旋定子工況基本一致,在存在葉片一側,預旋定子使得周向速度顯著增加,槳盤面外圍速度等值線變得密集,切向速度梯度增加,產(chǎn)生與螺旋槳旋向相反的預旋流,這能夠增加槳盤面進流,改善螺旋槳工作條件。

    圖7 槳盤面處艉流場周向速度對比Fig.7 Comparison diagram of tangential velocity of wake field at propeller plane

    3.2 船舶艉部伴流場漩渦分布

    圖8為槳盤面處無預旋定子船舶艉流場的漩渦強度(swirling strength)與渦量(vorticity(X))測量結果,依據(jù)船舶的對稱性,僅給出一半流場結果。其中,漩渦強度被定義為速度梯度張量J的復雜特征值虛部部分:

    圖8 無預旋定子工況槳盤面處漩渦分布Fig.8 Vortex distribution of wake field without pre-swirl stator

    (1)

    由于測的的結果為平面內數(shù)據(jù),沿船長方向的平面數(shù)據(jù)梯度不能被計算,把它們設為零可以簡化特征值計算,所以虛部的平方值可以計算為[16]:

    漩渦強度的局部最小負值可以用來識別渦核,而正值表示流場,剪切力可能會被顯示但沒有旋轉運動。旋渦強度的單位為1/s2。

    圍繞X軸的漩渦的渦量:

    (3)

    同樣由于X方向速度梯度無法計算,所以僅圍繞X軸的渦量可以顯示為:

    (4)

    渦量的單位為1/s。

    從圖8可以看出,舭渦與螺旋槳假轂轂帽渦可以被很好地識別,圖中渦量正值表示漩渦為逆時針旋轉,負值表示旋渦為順時針旋轉??梢钥闯龃S蚁咸幍聂皽u為逆時針旋轉,左舷處的舭渦為順時針旋轉,且是由船體中部舭部產(chǎn)生并傳遞至螺旋槳盤面。螺旋槳槳軸下方的假轂轂帽渦在右舷處為順時針旋轉,左舷為逆時針旋轉,與舭渦旋向方向相反,螺旋槳假轂轂帽渦由船體艉軸處產(chǎn)生且隨流體傳遞到螺旋槳盤面處,整體漩渦位置與范圍與圖5中的流線圖相互對應。

    圖9(a)為加裝預旋定子工況螺旋槳盤面處漩渦強度分布,圖9(b)為螺旋槳盤面處渦量分布結果。

    從圖8和9(a)中可以看出,無預旋定子的船艉伴流場體現(xiàn)出了較明顯的2組渦(舭渦及假轂轂帽渦),而有預旋定子的船艉伴流場整體上雖然也體現(xiàn)出此2組渦,但其周圍還混雜著一些強度較低的渦;同時,從圖9(b)中可以看出,已出現(xiàn)正反渦相互摻雜的現(xiàn)象,尤其是下方的一對假轂轂帽渦,考慮到螺旋槳旋向為右旋,而此紫色正值渦明顯有擴散趨勢,其旋向為左旋,以此可以看出,預旋定子能夠在槳前產(chǎn)生與槳旋向相反的預旋流,從而改善螺旋槳工作條件,提高推進效率。

    3.3 不同半徑處的軸向速度分布

    受槳軸遮擋作用影響,槳盤面位置處的流場存在數(shù)據(jù)缺失,因此取槳盤面后7.5 cm處流場作為有無預旋定子存在的軸向速度差異對比,從而分析出每一角度上有無預旋定子存在的軸向速度差異。

    根據(jù)軸向速度云圖,截取0.2R、0.4R、0.6R、0.8R、1.0R圓周范圍內若干點上的軸向速度,每一圓周與預旋定子的相對位置關系如圖10所示。將圓周視為鐘表,規(guī)定0°位置為12點鐘方向(指向船艉正上方),之后每間隔10°取一次點,從船艉前視,順時針方向為正方向。最終,將軸向速度u使用拖曳航速U作無量綱化處理。軸向速度對比結果如圖11所示。

    圖10 0.2R~1.0R圓周與預旋定子的相對位置關系Fig.10 Relative position relationship between 0.2R~1.0R circumference and pre rotating stator

    圖11 x=7.5 cm處不同半徑下軸向速度對比Fig.11 Comparison of axial velocity under different radius at x=7.5 cm

    從圖11中可以看出,隨著船艉伴流向更大半徑處發(fā)展,流場受船艉阻礙作用減小,對應船底部即180°速度提升最為明顯,其逐漸恢復至外流場流速。另外,受“鉤狀”速度分布影響,在0.6R~1.0R內0°~45°、315°~360°內總會存在先增加后降低的速度分布規(guī)律。相比較而言,0.2R~0.4R內的相同位置處速度分布較為平緩。

    對比有無預旋定子帶來的速度分布差異,在0.8R與1.0R的數(shù)據(jù)中可以看出,2條曲線基本重合,差異主要來自0.6R及以內的部分。裸船體的艉流場速度分布大致關于180°對稱;受非對稱定子阻礙作用影響,存在預旋定子工況的速度分布并不明顯對稱,0°~180°內,2條曲線比較吻合。然而在另半圓周內,180°對應的速度峰值偏移10°~20°,此角度對應于下方葉片位置角度,且存在預旋定子時此偏移的速度峰值較裸船體180°的速度峰值提高10%~20%。觀察上方2葉片附近位置的流場速度,與裸船體相同位置處的速度相比均有減少,但幅度較小,整體上下方葉片附近的流場差異更大一些。這與船舶艉部的舭渦、假轂榖帽渦的位置與旋向息息相關。不同葉片位置處相反的渦旋向導致了流場速度分布的不同。以上現(xiàn)象說明此預旋定子葉片布置形式,對于葉片存在的半圓區(qū)域(180°~360°)影響較大,下方葉片的存在對流場分布的影響更大,它導致了在0.6R內,下方葉片附近位置的流場加速現(xiàn)象,槳盤面流體進流速度增加,改善了伴流條件,對提升螺旋槳效率起到了積極作用。但同時,加裝預旋定子對于槳盤面的流場均勻性產(chǎn)生了不利影響,在設計時應充分考慮。

    4 結論

    1)PIV技術是船舶與海洋工程領域重要的測量技術手段之一。本文利用PIV技術精確捕捉到某散貨船艉流場中舭渦以及“鉤狀”速度輪廓等流場特征,“鉤狀”速度云圖等值線對應的速度值約為u/U=0.30。充分顯現(xiàn)了其具有瞬態(tài)、全局、無接觸以及高精度等測量優(yōu)點。

    2)在槳盤面位置處,預旋定子對艉流具有干擾作用,主要影響區(qū)域位于有葉片存在的左半部分,上方“鉤狀”結構被破壞,表現(xiàn)出速度等值線向槳盤面外部偏移。對于無葉片存在的右半部分,速度等值線具有明顯的“鉤狀”結構。預旋定子不僅僅影響葉片存在的一側,對沒有葉片的一側也會產(chǎn)生影響。預旋定子使得周向速度顯著增加,產(chǎn)生與螺旋槳旋向相反的預旋流,能夠增加槳盤面進流,改善螺旋槳工作條件。在渦量圖中,有預旋定子存在時已出現(xiàn)正反渦相互摻雜的現(xiàn)象,尤其是下方的一對假轂轂帽渦,由槳旋向以及渦的擴散趨勢,可以看出預旋定子能夠在槳前產(chǎn)生與槳旋向相反的預旋流。

    3)對比有無預旋定子存在的艉流場軸向速度分布,可以發(fā)現(xiàn)此種布置形式的預旋定子在0.6R范圍內的影響較大。受船艉伴流及艉渦的影響,中下方葉片對流場分布的影響更大,且造成附近速度峰值角度偏移、峰值增加,槳盤面流體進流速度增加,改善了伴流條件。

    猜你喜歡
    盤面螺旋槳流場
    大型空冷汽輪發(fā)電機轉子三維流場計算
    大電機技術(2021年2期)2021-07-21 07:28:24
    淺談制動盤面端跳動測量
    基于CFD的螺旋槳拉力確定方法
    轉杯紡排雜區(qū)流場與排雜性能
    基于HYCOM的斯里蘭卡南部海域溫、鹽、流場統(tǒng)計分析
    數(shù)獨競猜
    數(shù)獨競猜
    基于瞬態(tài)流場計算的滑動軸承靜平衡位置求解
    3800DWT加油船螺旋槳諧鳴分析及消除方法
    廣東造船(2015年6期)2015-02-27 10:52:46
    試論棋例裁決難點——無關聯(lián)①
    棋藝(2014年7期)2014-09-09 09:41:50
    亚洲成a人片在线一区二区| 日日干狠狠操夜夜爽| 欧美中文日本在线观看视频| 偷拍熟女少妇极品色| 精品一区二区三区视频在线观看免费| 午夜久久久久精精品| 女人十人毛片免费观看3o分钟| 全区人妻精品视频| 国产免费男女视频| 美女被艹到高潮喷水动态| 亚洲美女搞黄在线观看 | 99在线人妻在线中文字幕| 日本 av在线| 国产黄片美女视频| 日韩 亚洲 欧美在线| 国产亚洲91精品色在线| 波多野结衣高清无吗| 91狼人影院| a级一级毛片免费在线观看| 中亚洲国语对白在线视频| 一级毛片久久久久久久久女| 久久久久久大精品| 美女高潮喷水抽搐中文字幕| 精品久久久久久,| 日本五十路高清| av在线亚洲专区| 久久久精品大字幕| 欧美极品一区二区三区四区| 校园春色视频在线观看| 国产日本99.免费观看| 久99久视频精品免费| 丝袜美腿在线中文| 免费无遮挡裸体视频| 国内精品久久久久久久电影| 久久人人爽人人爽人人片va| 亚洲av二区三区四区| 午夜爱爱视频在线播放| 岛国在线免费视频观看| 久久久久久久久久成人| 欧美色欧美亚洲另类二区| av视频在线观看入口| 俄罗斯特黄特色一大片| 我要看日韩黄色一级片| 丰满乱子伦码专区| 有码 亚洲区| 国产精品一及| 97人妻精品一区二区三区麻豆| 亚洲专区中文字幕在线| 亚洲成人免费电影在线观看| 日韩在线高清观看一区二区三区 | 久久精品久久久久久噜噜老黄 | 性插视频无遮挡在线免费观看| 久久久久久久久大av| 欧美三级亚洲精品| 亚洲三级黄色毛片| 免费在线观看成人毛片| 舔av片在线| 丰满的人妻完整版| 又粗又爽又猛毛片免费看| 白带黄色成豆腐渣| 国产av麻豆久久久久久久| 成人亚洲精品av一区二区| 久久久久九九精品影院| 久久久久久久久大av| 亚洲成人久久性| 在线天堂最新版资源| 午夜免费男女啪啪视频观看 | 国产一区二区亚洲精品在线观看| 日韩欧美三级三区| 国产免费av片在线观看野外av| 久久婷婷人人爽人人干人人爱| 国产蜜桃级精品一区二区三区| ponron亚洲| 97碰自拍视频| 夜夜看夜夜爽夜夜摸| 亚洲久久久久久中文字幕| 久久久久九九精品影院| 久久99热这里只有精品18| 日韩中字成人| 国产一区二区三区av在线 | 国产毛片a区久久久久| 欧美日韩瑟瑟在线播放| 国产免费男女视频| 国产熟女欧美一区二区| 成人特级黄色片久久久久久久| 黄色一级大片看看| 赤兔流量卡办理| 少妇的逼水好多| 小蜜桃在线观看免费完整版高清| 精品久久久噜噜| 久久久精品大字幕| 国产亚洲精品综合一区在线观看| 日日夜夜操网爽| 免费搜索国产男女视频| 国产精品综合久久久久久久免费| 又爽又黄a免费视频| 一本一本综合久久| 亚洲一级一片aⅴ在线观看| 国产成人一区二区在线| 内地一区二区视频在线| 97热精品久久久久久| 亚洲黑人精品在线| 精品免费久久久久久久清纯| 日本 欧美在线| 女人十人毛片免费观看3o分钟| 亚洲国产精品成人综合色| 欧美日韩中文字幕国产精品一区二区三区| 国产一区二区激情短视频| 久久久久久国产a免费观看| 欧美日韩中文字幕国产精品一区二区三区| 淫秽高清视频在线观看| 在线观看美女被高潮喷水网站| 九色成人免费人妻av| 日韩一本色道免费dvd| 又黄又爽又刺激的免费视频.| 久久久久九九精品影院| 99久久久亚洲精品蜜臀av| 日日夜夜操网爽| 午夜福利视频1000在线观看| 在线免费观看的www视频| 日韩精品中文字幕看吧| 久久久久久久亚洲中文字幕| 亚洲avbb在线观看| 国产视频内射| 日本 欧美在线| 日韩欧美三级三区| 免费电影在线观看免费观看| 日韩人妻高清精品专区| av天堂在线播放| 亚洲专区中文字幕在线| 中文亚洲av片在线观看爽| 又黄又爽又免费观看的视频| 国产视频一区二区在线看| 欧美区成人在线视频| 亚洲 国产 在线| 婷婷色综合大香蕉| 国产免费av片在线观看野外av| 国产精品,欧美在线| АⅤ资源中文在线天堂| 国产黄a三级三级三级人| 中文字幕熟女人妻在线| 色综合站精品国产| www.色视频.com| 国产精品久久久久久亚洲av鲁大| 在线a可以看的网站| 国产熟女欧美一区二区| 91午夜精品亚洲一区二区三区 | 亚洲精品在线观看二区| 嫩草影视91久久| 亚洲成人中文字幕在线播放| 97人妻精品一区二区三区麻豆| 亚洲av中文av极速乱 | 久久久色成人| 久久99热6这里只有精品| 91麻豆av在线| 国产又黄又爽又无遮挡在线| 大又大粗又爽又黄少妇毛片口| 可以在线观看的亚洲视频| 在线观看美女被高潮喷水网站| 午夜福利在线观看吧| 在现免费观看毛片| 美女免费视频网站| 三级国产精品欧美在线观看| 日日夜夜操网爽| 女人被狂操c到高潮| 欧美另类亚洲清纯唯美| 国产精品久久视频播放| 国产高清视频在线观看网站| 村上凉子中文字幕在线| 91av网一区二区| 午夜a级毛片| 香蕉av资源在线| 国产亚洲91精品色在线| 嫩草影院精品99| 久久6这里有精品| 国产亚洲av嫩草精品影院| 午夜精品在线福利| 很黄的视频免费| 老司机福利观看| 久久久久久久亚洲中文字幕| 国产探花在线观看一区二区| 最近最新中文字幕大全电影3| 99热6这里只有精品| 一区二区三区免费毛片| 两个人视频免费观看高清| 久久久久久久久大av| av国产免费在线观看| 日韩欧美精品v在线| av中文乱码字幕在线| 国产精品福利在线免费观看| 日本与韩国留学比较| 亚洲自拍偷在线| 日本色播在线视频| 国产一区二区在线观看日韩| 久久这里只有精品中国| 我的老师免费观看完整版| 国产一级毛片七仙女欲春2| 人妻夜夜爽99麻豆av| av天堂在线播放| 国产 一区精品| 国产真实伦视频高清在线观看 | 狠狠狠狠99中文字幕| 夜夜爽天天搞| 国产视频内射| 亚洲人成网站高清观看| 欧美精品啪啪一区二区三区| 成人国产综合亚洲| 日韩人妻高清精品专区| 国产伦精品一区二区三区视频9| eeuss影院久久| 我的女老师完整版在线观看| 精品国产三级普通话版| 色噜噜av男人的天堂激情| 久久久久久伊人网av| 国产精品精品国产色婷婷| 免费人成视频x8x8入口观看| 亚洲电影在线观看av| 国产精品乱码一区二三区的特点| 日本成人三级电影网站| 不卡一级毛片| 一a级毛片在线观看| 久久精品影院6| 最近中文字幕高清免费大全6 | 干丝袜人妻中文字幕| 尤物成人国产欧美一区二区三区| 午夜视频国产福利| 成年版毛片免费区| 一进一出抽搐gif免费好疼| avwww免费| 午夜福利在线在线| 免费在线观看日本一区| 日本三级黄在线观看| 国产视频内射| 又爽又黄无遮挡网站| 久久国内精品自在自线图片| 久久国内精品自在自线图片| 精品久久久久久成人av| 国产高清三级在线| 国产高清三级在线| 在线播放无遮挡| 搡老熟女国产l中国老女人| 夜夜看夜夜爽夜夜摸| www日本黄色视频网| 国产一区二区亚洲精品在线观看| 日韩高清综合在线| 欧美丝袜亚洲另类 | 亚洲精品亚洲一区二区| 日本 av在线| 国产一区二区亚洲精品在线观看| 嫩草影院入口| 一级黄片播放器| 中文字幕人妻熟人妻熟丝袜美| 老师上课跳d突然被开到最大视频| 国产成人a区在线观看| 又黄又爽又刺激的免费视频.| 国产色婷婷99| 亚洲久久久久久中文字幕| 亚洲在线自拍视频| 免费一级毛片在线播放高清视频| 国产亚洲精品综合一区在线观看| 天天一区二区日本电影三级| h日本视频在线播放| 麻豆久久精品国产亚洲av| 亚洲av中文av极速乱 | 国内精品宾馆在线| 美女高潮的动态| 精品99又大又爽又粗少妇毛片 | 无遮挡黄片免费观看| 18禁黄网站禁片免费观看直播| 国产男人的电影天堂91| 欧美性猛交黑人性爽| 69人妻影院| 国产精品一区www在线观看 | 成人国产一区最新在线观看| 男人舔女人下体高潮全视频| 一级毛片久久久久久久久女| 少妇人妻精品综合一区二区 | av天堂中文字幕网| 午夜福利欧美成人| 91在线精品国自产拍蜜月| 欧美色视频一区免费| 人妻久久中文字幕网| 欧美日韩亚洲国产一区二区在线观看| 国产成人a区在线观看| 亚洲,欧美,日韩| 久久欧美精品欧美久久欧美| 88av欧美| 亚洲va在线va天堂va国产| 麻豆精品久久久久久蜜桃| 中文在线观看免费www的网站| 在线观看66精品国产| 3wmmmm亚洲av在线观看| 搡老妇女老女人老熟妇| 最新在线观看一区二区三区| 午夜精品一区二区三区免费看| 国产精品亚洲一级av第二区| 久久精品夜夜夜夜夜久久蜜豆| 一夜夜www| 97超级碰碰碰精品色视频在线观看| 国产伦精品一区二区三区四那| 亚洲中文字幕一区二区三区有码在线看| 亚洲美女视频黄频| 69av精品久久久久久| 久久精品国产亚洲网站| 啦啦啦韩国在线观看视频| 国产亚洲精品久久久久久毛片| 日韩中字成人| 在线观看美女被高潮喷水网站| 乱人视频在线观看| 在线观看午夜福利视频| www.www免费av| 麻豆av噜噜一区二区三区| 欧美3d第一页| 日本一本二区三区精品| 亚洲在线观看片| 免费av观看视频| 欧美激情国产日韩精品一区| 白带黄色成豆腐渣| 亚洲专区中文字幕在线| 国产私拍福利视频在线观看| 91在线精品国自产拍蜜月| 如何舔出高潮| 久久久久国产精品人妻aⅴ院| 国内久久婷婷六月综合欲色啪| 人人妻人人看人人澡| 97碰自拍视频| 99九九线精品视频在线观看视频| 国产69精品久久久久777片| 在线国产一区二区在线| 成人国产一区最新在线观看| 日韩强制内射视频| 99热这里只有是精品在线观看| 欧美高清成人免费视频www| 国产三级中文精品| 少妇的逼好多水| 日本 欧美在线| 久久久国产成人免费| 亚洲人成网站在线播| 成人欧美大片| 在线播放无遮挡| 国产一区二区三区视频了| av在线观看视频网站免费| 变态另类丝袜制服| 最后的刺客免费高清国语| 婷婷丁香在线五月| 中文资源天堂在线| 久久人人精品亚洲av| 岛国在线免费视频观看| 亚洲精品乱码久久久v下载方式| 内射极品少妇av片p| 他把我摸到了高潮在线观看| 欧美成人一区二区免费高清观看| 亚洲男人的天堂狠狠| 如何舔出高潮| 免费在线观看影片大全网站| 午夜亚洲福利在线播放| 亚洲图色成人| 99国产极品粉嫩在线观看| 精品久久久久久久末码| 亚洲中文日韩欧美视频| 午夜a级毛片| 日本一二三区视频观看| 中文资源天堂在线| 草草在线视频免费看| 精品人妻熟女av久视频| 国产午夜精品论理片| 搡女人真爽免费视频火全软件 | 国产毛片a区久久久久| 亚洲国产欧洲综合997久久,| 99久国产av精品| 国内精品一区二区在线观看| 久久这里只有精品中国| 国产亚洲精品久久久com| 窝窝影院91人妻| 国产乱人视频| 毛片一级片免费看久久久久 | 日韩欧美国产在线观看| 亚洲国产日韩欧美精品在线观看| 亚洲精品日韩av片在线观看| 亚洲四区av| netflix在线观看网站| 99久久精品热视频| 久久久色成人| 国产黄片美女视频| 我要看日韩黄色一级片| 欧美国产日韩亚洲一区| 我的老师免费观看完整版| 国产一区二区亚洲精品在线观看| 亚洲一区高清亚洲精品| 免费观看在线日韩| 九九热线精品视视频播放| 亚洲人成网站在线播| 国产高潮美女av| 国产高清有码在线观看视频| 国产成人a区在线观看| 亚洲电影在线观看av| 亚洲最大成人av| 亚洲国产色片| 国产精品国产高清国产av| 变态另类丝袜制服| 久久精品国产亚洲av天美| 我的老师免费观看完整版| 日本 av在线| 简卡轻食公司| 国产精品99久久久久久久久| 自拍偷自拍亚洲精品老妇| 亚洲av一区综合| 中文在线观看免费www的网站| 日韩精品青青久久久久久| 久久6这里有精品| 免费在线观看日本一区| 亚洲精品日韩av片在线观看| 99久久无色码亚洲精品果冻| 天堂网av新在线| 欧美日本亚洲视频在线播放| 春色校园在线视频观看| 日韩精品中文字幕看吧| 精品久久国产蜜桃| 午夜爱爱视频在线播放| 久久久久久久久久成人| 成人国产麻豆网| 在线观看午夜福利视频| 99热这里只有是精品50| av在线蜜桃| 91久久精品国产一区二区三区| 一级黄片播放器| 内射极品少妇av片p| 一a级毛片在线观看| 久久午夜福利片| 国产精品一区二区免费欧美| 亚洲自拍偷在线| 欧美色视频一区免费| 国产精品乱码一区二三区的特点| 我的女老师完整版在线观看| 国产精品一区www在线观看 | 免费看a级黄色片| 99久久无色码亚洲精品果冻| 国产人妻一区二区三区在| 亚洲av成人精品一区久久| 自拍偷自拍亚洲精品老妇| 日本爱情动作片www.在线观看 | 久久午夜福利片| 国产精品一区二区免费欧美| 亚洲欧美日韩无卡精品| 最近在线观看免费完整版| 成人精品一区二区免费| 精品国产三级普通话版| 熟女人妻精品中文字幕| 国产高清不卡午夜福利| 国产一区二区激情短视频| 午夜福利成人在线免费观看| 在线观看66精品国产| 九色成人免费人妻av| 亚洲欧美精品综合久久99| 日韩中文字幕欧美一区二区| 亚洲,欧美,日韩| 亚洲va在线va天堂va国产| 中文字幕av在线有码专区| 国产一区二区三区av在线 | a级一级毛片免费在线观看| 在线观看一区二区三区| 亚洲av免费高清在线观看| 欧美人与善性xxx| 久久久久久久久久黄片| 免费看日本二区| 成年美女黄网站色视频大全免费 | 成年女人在线观看亚洲视频| 亚洲熟女精品中文字幕| 亚洲国产色片| 色综合色国产| 亚洲内射少妇av| 国产人妻一区二区三区在| 亚洲丝袜综合中文字幕| 亚洲av欧美aⅴ国产| 国产亚洲午夜精品一区二区久久| 狂野欧美白嫩少妇大欣赏| 国产伦精品一区二区三区四那| 中文字幕av成人在线电影| 有码 亚洲区| 超碰97精品在线观看| 亚洲久久久国产精品| av国产免费在线观看| 国产成人免费无遮挡视频| 97在线人人人人妻| 国产精品福利在线免费观看| 春色校园在线视频观看| 国产伦在线观看视频一区| 啦啦啦在线观看免费高清www| 美女视频免费永久观看网站| av在线蜜桃| 国产免费视频播放在线视频| 80岁老熟妇乱子伦牲交| 国产精品成人在线| 日本一二三区视频观看| 国产探花极品一区二区| 亚洲精品色激情综合| 欧美xxxx性猛交bbbb| 久久国产精品大桥未久av | 自拍欧美九色日韩亚洲蝌蚪91 | 麻豆精品久久久久久蜜桃| 人体艺术视频欧美日本| 偷拍熟女少妇极品色| 新久久久久国产一级毛片| 22中文网久久字幕| 一二三四中文在线观看免费高清| 日日啪夜夜爽| 婷婷色麻豆天堂久久| 亚洲国产最新在线播放| 夜夜爽夜夜爽视频| 中文欧美无线码| 18禁在线无遮挡免费观看视频| 久久精品熟女亚洲av麻豆精品| 黄片无遮挡物在线观看| 成人无遮挡网站| 国产又色又爽无遮挡免| 亚洲美女视频黄频| 人妻少妇偷人精品九色| 熟女av电影| 最近最新中文字幕免费大全7| 国产精品一二三区在线看| 你懂的网址亚洲精品在线观看| 美女xxoo啪啪120秒动态图| 极品少妇高潮喷水抽搐| 男的添女的下面高潮视频| 国产精品精品国产色婷婷| 亚洲第一av免费看| 26uuu在线亚洲综合色| 亚洲精品国产av成人精品| 亚洲欧美成人精品一区二区| 一级毛片 在线播放| 国产精品久久久久久精品古装| 永久免费av网站大全| 啦啦啦视频在线资源免费观看| 成人午夜精彩视频在线观看| 国产精品国产三级国产av玫瑰| 91精品国产九色| 亚洲精品亚洲一区二区| 97热精品久久久久久| 精品一品国产午夜福利视频| www.色视频.com| 99热这里只有是精品在线观看| 纵有疾风起免费观看全集完整版| 丰满少妇做爰视频| 女人久久www免费人成看片| 观看美女的网站| 人人妻人人看人人澡| 国产国拍精品亚洲av在线观看| 老师上课跳d突然被开到最大视频| 欧美日韩精品成人综合77777| 秋霞在线观看毛片| 日韩一本色道免费dvd| 亚洲国产色片| 成年美女黄网站色视频大全免费 | 韩国av在线不卡| 高清午夜精品一区二区三区| 舔av片在线| 久久久久久久亚洲中文字幕| 97超视频在线观看视频| 精品人妻一区二区三区麻豆| 成人国产麻豆网| 亚洲av不卡在线观看| 精品人妻视频免费看| 国产亚洲一区二区精品| 成人毛片a级毛片在线播放| 国产免费一区二区三区四区乱码| 日本黄色日本黄色录像| 久久99热这里只频精品6学生| 尾随美女入室| 成人影院久久| 97超碰精品成人国产| 国产高清国产精品国产三级 | 亚洲国产av新网站| 国产精品久久久久久久久免| 少妇人妻一区二区三区视频| 一级黄片播放器| 国产亚洲午夜精品一区二区久久| 欧美成人一区二区免费高清观看| 精品酒店卫生间| 女人十人毛片免费观看3o分钟| 人妻一区二区av| 丰满迷人的少妇在线观看| 亚洲精华国产精华液的使用体验| 最黄视频免费看| 少妇 在线观看| 日日啪夜夜撸| 午夜激情福利司机影院| 国产永久视频网站| 国产黄色视频一区二区在线观看| 精品熟女少妇av免费看| 91精品国产国语对白视频| 赤兔流量卡办理| 国产国拍精品亚洲av在线观看| 国产在线一区二区三区精| 亚洲欧美精品专区久久| 最近手机中文字幕大全| 人妻制服诱惑在线中文字幕| 搡老乐熟女国产| 在线观看免费视频网站a站| av免费在线看不卡| 欧美高清性xxxxhd video| 欧美日韩在线观看h| 国产人妻一区二区三区在| 亚洲电影在线观看av| 免费观看无遮挡的男女| 日日摸夜夜添夜夜添av毛片| 亚洲精品国产色婷婷电影| 青春草视频在线免费观看| 一级a做视频免费观看| 人人妻人人看人人澡| 伦理电影免费视频| 尾随美女入室| 嫩草影院新地址| 黄色一级大片看看| 少妇猛男粗大的猛烈进出视频| 色综合色国产|