• <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
    国产男女内射视频| 高清欧美精品videossex| 国产乱来视频区| 国产高清三级在线| 亚洲欧美精品自产自拍| 精品国产一区二区久久| 色哟哟·www| 成人无遮挡网站| 如何舔出高潮| 日韩电影二区| 久久久久久久久久久久大奶| 欧美国产精品一级二级三级| 亚洲欧美清纯卡通| 香蕉精品网在线| 久久亚洲国产成人精品v| 国产精品 国内视频| 在线观看三级黄色| av不卡在线播放| 亚洲国产精品国产精品| 亚洲经典国产精华液单| 中文字幕人妻熟女乱码| 九色成人免费人妻av| 国产成人a∨麻豆精品| 久久99热6这里只有精品| 久热这里只有精品99| 亚洲国产精品一区二区三区在线| 少妇的逼好多水| av不卡在线播放| 丝袜喷水一区| 日本-黄色视频高清免费观看| 亚洲综合色惰| 久久午夜福利片| 内地一区二区视频在线| 大陆偷拍与自拍| 亚洲美女视频黄频| 丰满少妇做爰视频| 成人综合一区亚洲| 午夜福利在线观看免费完整高清在| 国产成人91sexporn| 国产成人a∨麻豆精品| 久久久久久人妻| 成年人免费黄色播放视频| 亚洲一码二码三码区别大吗| 欧美97在线视频| av女优亚洲男人天堂| 9热在线视频观看99| 97精品久久久久久久久久精品| 国产亚洲一区二区精品| 久久人人97超碰香蕉20202| tube8黄色片| 18禁国产床啪视频网站| 亚洲,欧美,日韩| 女人久久www免费人成看片| 最近最新中文字幕大全免费视频 | 国产精品久久久久久久电影| 菩萨蛮人人尽说江南好唐韦庄| 啦啦啦视频在线资源免费观看| 国产精品久久久久久久久免| 少妇高潮的动态图| 亚洲国产av影院在线观看| 精品少妇内射三级| 国产精品久久久av美女十八| av网站免费在线观看视频| 国产成人精品婷婷| 丰满饥渴人妻一区二区三| 人妻一区二区av| 亚洲美女黄色视频免费看| 综合色丁香网| 天天躁夜夜躁狠狠躁躁| 国产av码专区亚洲av| 日日撸夜夜添| 9热在线视频观看99| 80岁老熟妇乱子伦牲交| 日本av免费视频播放| 校园人妻丝袜中文字幕| 国产爽快片一区二区三区| 国产在线免费精品| 少妇的丰满在线观看| 五月天丁香电影| 秋霞伦理黄片| 国产一区二区三区av在线| 久久久久久久精品精品| av黄色大香蕉| 久久久久久久久久久久大奶| 日本av手机在线免费观看| 国产精品蜜桃在线观看| 国产av码专区亚洲av| 国产一区亚洲一区在线观看| 秋霞在线观看毛片| 人人妻人人澡人人爽人人夜夜| 菩萨蛮人人尽说江南好唐韦庄| 国产在视频线精品| 午夜视频国产福利| 最黄视频免费看| 91精品伊人久久大香线蕉| 99久久人妻综合| 九九在线视频观看精品| 午夜激情久久久久久久| 丰满饥渴人妻一区二区三| 日本黄色日本黄色录像| 精品一区二区三区四区五区乱码 | 亚洲精品美女久久久久99蜜臀 | 免费少妇av软件| 欧美精品一区二区大全| 观看美女的网站| 中文字幕制服av| 久久国产精品大桥未久av| 色吧在线观看| 少妇人妻精品综合一区二区| 色5月婷婷丁香| 免费观看a级毛片全部| 黄色毛片三级朝国网站| 丝袜喷水一区| 久久99蜜桃精品久久| 最近中文字幕高清免费大全6| 婷婷色av中文字幕| 国产精品久久久久久久电影| 中国国产av一级| 蜜桃国产av成人99| 黑丝袜美女国产一区| 各种免费的搞黄视频| 丝袜脚勾引网站| 18禁国产床啪视频网站| 亚洲av男天堂| 麻豆乱淫一区二区| 夫妻性生交免费视频一级片| 有码 亚洲区| 国产精品一二三区在线看| 九草在线视频观看| 我的女老师完整版在线观看| 天堂俺去俺来也www色官网| 亚洲精品美女久久av网站| 久久这里只有精品19| 夜夜骑夜夜射夜夜干| 青春草视频在线免费观看| 欧美亚洲 丝袜 人妻 在线| 91在线精品国自产拍蜜月| 久久精品国产亚洲av天美| av.在线天堂| www.熟女人妻精品国产 | 99久久中文字幕三级久久日本| 老司机亚洲免费影院| 国产av一区二区精品久久| 欧美日韩成人在线一区二区| 午夜91福利影院| 精品一区二区三卡| 国产欧美日韩综合在线一区二区| 国产在视频线精品| 在线精品无人区一区二区三| 我的女老师完整版在线观看| 婷婷色av中文字幕| 捣出白浆h1v1| 爱豆传媒免费全集在线观看| 国产成人欧美| 国产高清不卡午夜福利| 高清av免费在线| 国产精品一二三区在线看| 90打野战视频偷拍视频| 成人影院久久| 纯流量卡能插随身wifi吗| av.在线天堂| 欧美精品av麻豆av| 极品少妇高潮喷水抽搐| 国产一区二区三区av在线| 大片免费播放器 马上看| 美女视频免费永久观看网站| 伊人久久国产一区二区| 啦啦啦视频在线资源免费观看| 久久久久久久精品精品| 韩国高清视频一区二区三区| 久热这里只有精品99| 久久久久久久久久久久大奶| www.色视频.com| 午夜福利视频精品| 热re99久久精品国产66热6| 久久精品熟女亚洲av麻豆精品| 亚洲精品美女久久av网站| 亚洲欧美成人综合另类久久久| 夜夜骑夜夜射夜夜干| 日本欧美视频一区| 亚洲成人av在线免费| 26uuu在线亚洲综合色| 黄色怎么调成土黄色| 高清在线视频一区二区三区| 中文欧美无线码| 99久久综合免费| 午夜免费观看性视频| 欧美人与性动交α欧美精品济南到 | 大话2 男鬼变身卡| 精品少妇久久久久久888优播| 老女人水多毛片| 亚洲精品日韩在线中文字幕| 咕卡用的链子| 考比视频在线观看| 日韩精品有码人妻一区| 人妻系列 视频| 亚洲经典国产精华液单| 免费av中文字幕在线| 男人舔女人的私密视频| 人成视频在线观看免费观看| 亚洲成人手机| 日韩欧美一区视频在线观看| 亚洲综合色惰| 久久久久久人人人人人| 1024视频免费在线观看| 又黄又粗又硬又大视频| 国产成人精品福利久久| 精品第一国产精品| 秋霞伦理黄片| 日韩三级伦理在线观看| 免费黄网站久久成人精品| 美女大奶头黄色视频| 国产成人免费无遮挡视频| 成年动漫av网址| 久久精品国产自在天天线| 欧美国产精品va在线观看不卡| 亚洲在久久综合| 日韩三级伦理在线观看| av线在线观看网站| 欧美人与性动交α欧美精品济南到 | 日日爽夜夜爽网站| 国产成人午夜福利电影在线观看| 一本—道久久a久久精品蜜桃钙片| 国产日韩欧美视频二区| 国产福利在线免费观看视频| 精品卡一卡二卡四卡免费| 免费看不卡的av| www.色视频.com| 建设人人有责人人尽责人人享有的| 最近中文字幕2019免费版| 亚洲av电影在线观看一区二区三区| 在线观看人妻少妇| 一区二区av电影网| 婷婷色av中文字幕| 18在线观看网站| 丰满少妇做爰视频| 性高湖久久久久久久久免费观看| 亚洲情色 制服丝袜| 一个人免费看片子| 草草在线视频免费看| 在线免费观看不下载黄p国产| 亚洲成色77777| 搡老乐熟女国产| 91成人精品电影| 成人手机av| 国产片特级美女逼逼视频| 国产在线免费精品| 校园人妻丝袜中文字幕| 日韩成人伦理影院| √禁漫天堂资源中文www| 丝瓜视频免费看黄片| 午夜日本视频在线| 美女国产高潮福利片在线看| 九九在线视频观看精品| 亚洲成人手机| 国产在线免费精品| 在线观看国产h片| 国产又爽黄色视频| 搡女人真爽免费视频火全软件| 国产黄色视频一区二区在线观看| 成年美女黄网站色视频大全免费| 日韩欧美精品免费久久| 18禁国产床啪视频网站| 精品久久久久久电影网| 王馨瑶露胸无遮挡在线观看| tube8黄色片| 亚洲精品乱码久久久久久按摩| 一级毛片我不卡| 亚洲久久久国产精品| 自拍欧美九色日韩亚洲蝌蚪91| 成年av动漫网址| 汤姆久久久久久久影院中文字幕| 久久精品夜色国产| 日韩人妻精品一区2区三区| videossex国产| 国产成人一区二区在线| 亚洲国产精品一区三区| 久久毛片免费看一区二区三区| 精品久久蜜臀av无| 免费黄频网站在线观看国产| 日本vs欧美在线观看视频| 999精品在线视频| 日韩一本色道免费dvd| av不卡在线播放| 看十八女毛片水多多多| 亚洲色图综合在线观看| 亚洲欧美成人精品一区二区| 国国产精品蜜臀av免费| 精品视频人人做人人爽| 一级片'在线观看视频| 国产在线一区二区三区精| 免费观看a级毛片全部| 亚洲精品第二区| 国产成人精品无人区| 性色avwww在线观看| 久久久久国产网址| av国产精品久久久久影院| 宅男免费午夜| av一本久久久久| 美女中出高潮动态图| www.色视频.com| 国产精品蜜桃在线观看| 欧美日韩精品成人综合77777| 亚洲精品第二区| 午夜福利网站1000一区二区三区| 日本欧美国产在线视频| 天堂中文最新版在线下载| 成人无遮挡网站| 少妇高潮的动态图| 在线观看免费高清a一片| 久久热在线av| kizo精华| 制服人妻中文乱码| 国产精品人妻久久久久久| 免费久久久久久久精品成人欧美视频 | 国产探花极品一区二区| 国产综合精华液| 午夜视频国产福利| 男女国产视频网站| 国产淫语在线视频| 一级a做视频免费观看| 国产日韩欧美在线精品| xxx大片免费视频| 高清在线视频一区二区三区| freevideosex欧美| 免费观看性生交大片5| 国产乱人偷精品视频| 22中文网久久字幕| 人人妻人人添人人爽欧美一区卜| 自线自在国产av| 狠狠婷婷综合久久久久久88av| 狂野欧美激情性xxxx在线观看| 在线免费观看不下载黄p国产| 日韩,欧美,国产一区二区三区| 我的女老师完整版在线观看| 韩国av在线不卡| 国产成人精品在线电影| 男女高潮啪啪啪动态图| 熟妇人妻不卡中文字幕| 成人毛片a级毛片在线播放| 黑人欧美特级aaaaaa片| 自线自在国产av| 精品熟女少妇av免费看| 美女xxoo啪啪120秒动态图| 高清毛片免费看| 亚洲av国产av综合av卡| 一区二区三区四区激情视频| 精品一区二区三卡| 成人黄色视频免费在线看| 国产免费一区二区三区四区乱码| 精品视频人人做人人爽| 婷婷色综合大香蕉| 国产亚洲精品久久久com| 伦精品一区二区三区| 亚洲精品色激情综合| 国产无遮挡羞羞视频在线观看| 美女脱内裤让男人舔精品视频| 久久久久人妻精品一区果冻| 国产乱人偷精品视频| 91国产中文字幕| 亚洲综合精品二区| 在线精品无人区一区二区三| 国产精品人妻久久久影院| 亚洲少妇的诱惑av| 一本—道久久a久久精品蜜桃钙片| av在线app专区| 久久这里只有精品19| 老熟女久久久| 国产精品国产三级国产av玫瑰| 日日啪夜夜爽| 亚洲av.av天堂| 97在线人人人人妻| 五月玫瑰六月丁香| 久久精品aⅴ一区二区三区四区 | 久久ye,这里只有精品| 少妇被粗大的猛进出69影院 | 少妇精品久久久久久久| 日韩av免费高清视频| 亚洲三级黄色毛片| 欧美少妇被猛烈插入视频| 少妇精品久久久久久久| 日本av免费视频播放| 夜夜骑夜夜射夜夜干| 久久亚洲国产成人精品v| 桃花免费在线播放| 丰满饥渴人妻一区二区三| 最近最新中文字幕大全免费视频 | 九九爱精品视频在线观看| 婷婷色麻豆天堂久久| 午夜福利视频精品| 精品一区二区三卡| 亚洲天堂av无毛| 亚洲av免费高清在线观看| 免费大片黄手机在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久久久国产电影| 欧美亚洲 丝袜 人妻 在线| 国产亚洲精品久久久com| 伦理电影大哥的女人| 久久影院123| 亚洲三级黄色毛片| 草草在线视频免费看| 丰满饥渴人妻一区二区三| 纵有疾风起免费观看全集完整版| 少妇猛男粗大的猛烈进出视频| 国产精品不卡视频一区二区| 亚洲av.av天堂| 纯流量卡能插随身wifi吗| 国产精品偷伦视频观看了| 女性被躁到高潮视频| 中文天堂在线官网| 韩国av在线不卡| 亚洲美女黄色视频免费看| 欧美 日韩 精品 国产| av福利片在线| 日本91视频免费播放| 黑人高潮一二区| 亚洲丝袜综合中文字幕| 久久久久人妻精品一区果冻| 黄色配什么色好看| 啦啦啦视频在线资源免费观看| 国产不卡av网站在线观看| 五月玫瑰六月丁香| 国产xxxxx性猛交| 高清av免费在线| 国产成人午夜福利电影在线观看| 久久99热这里只频精品6学生| 捣出白浆h1v1| 一级毛片我不卡| 久久99一区二区三区| 国产精品一区www在线观看| 欧美日韩成人在线一区二区| 狠狠精品人妻久久久久久综合| 一区二区av电影网| 国产xxxxx性猛交| 精品一区在线观看国产| 精品人妻熟女毛片av久久网站| 美女福利国产在线| 成年美女黄网站色视频大全免费| 少妇 在线观看| 久久99精品国语久久久| 久久久久久人人人人人| 亚洲成国产人片在线观看| h视频一区二区三区| 欧美精品av麻豆av| 大香蕉久久网| 国产精品熟女久久久久浪| 亚洲精品乱码久久久久久按摩| 亚洲av综合色区一区| 精品国产一区二区久久| 久久狼人影院| 久久热在线av| 亚洲av电影在线进入| 老司机亚洲免费影院| 日韩av不卡免费在线播放| 国产精品一二三区在线看| 精品国产乱码久久久久久小说| videos熟女内射| 色哟哟·www| 亚洲成av片中文字幕在线观看 | 91aial.com中文字幕在线观看| 中文字幕免费在线视频6| 人人妻人人澡人人看| 丁香六月天网| 国产片内射在线| 国产欧美亚洲国产| 久久国内精品自在自线图片| 最新的欧美精品一区二区| 国产一级毛片在线| 亚洲成国产人片在线观看| 少妇的逼水好多| 久久影院123| 国产不卡av网站在线观看| 国产亚洲午夜精品一区二区久久| 夜夜爽夜夜爽视频| 欧美精品国产亚洲| 伦精品一区二区三区| 亚洲高清免费不卡视频| 久久99精品国语久久久| 欧美国产精品va在线观看不卡| 亚洲精品自拍成人| 女的被弄到高潮叫床怎么办| 国产一区亚洲一区在线观看| 久久av网站| 免费观看无遮挡的男女| 国产福利在线免费观看视频| 成人二区视频| 女人精品久久久久毛片| 国产精品人妻久久久久久| 国产成人免费无遮挡视频| 国产免费福利视频在线观看| 超碰97精品在线观看| 国产一区二区在线观看av| 久久狼人影院| av国产精品久久久久影院| 免费观看在线日韩| 久久精品久久久久久久性| 国产成人91sexporn| 成年女人在线观看亚洲视频| 久久精品久久精品一区二区三区| 一级毛片黄色毛片免费观看视频| 在线免费观看不下载黄p国产| 久久精品久久久久久噜噜老黄| 国产日韩欧美在线精品| 综合色丁香网| 久久人妻熟女aⅴ| 建设人人有责人人尽责人人享有的| 日本wwww免费看| 精品视频人人做人人爽| xxxhd国产人妻xxx| 免费大片黄手机在线观看| 少妇猛男粗大的猛烈进出视频| 欧美最新免费一区二区三区| 毛片一级片免费看久久久久| 国产成人精品福利久久| 国精品久久久久久国模美| 日韩,欧美,国产一区二区三区| 综合色丁香网| 精品人妻熟女毛片av久久网站| 国产男女内射视频| 赤兔流量卡办理| 久久ye,这里只有精品| 国产伦理片在线播放av一区| 美女xxoo啪啪120秒动态图| 在线精品无人区一区二区三| 久久精品久久精品一区二区三区| 九草在线视频观看| 热re99久久国产66热| 波野结衣二区三区在线| 日韩欧美精品免费久久| 欧美日韩一区二区视频在线观看视频在线| 精品久久久精品久久久| 欧美xxⅹ黑人| 久久久久精品性色| 母亲3免费完整高清在线观看 | 久久97久久精品| 美女脱内裤让男人舔精品视频| 99热全是精品| 十分钟在线观看高清视频www| 免费看不卡的av| 99九九在线精品视频| 婷婷色麻豆天堂久久| 高清在线视频一区二区三区| 久久影院123| 你懂的网址亚洲精品在线观看| 人妻系列 视频| 成人国产av品久久久| 丝袜在线中文字幕| 午夜激情久久久久久久| av有码第一页| a 毛片基地| 亚洲欧美成人综合另类久久久| 人妻 亚洲 视频| 亚洲婷婷狠狠爱综合网| 欧美亚洲 丝袜 人妻 在线| 久久久国产欧美日韩av| 观看美女的网站| 满18在线观看网站| 18禁观看日本| 国产亚洲最大av| 人体艺术视频欧美日本| 另类精品久久| 国产亚洲欧美精品永久| 少妇高潮的动态图| 亚洲欧美日韩另类电影网站| www.色视频.com| 人妻一区二区av| 国产精品嫩草影院av在线观看| 波多野结衣一区麻豆| 日韩精品免费视频一区二区三区 | 亚洲国产精品999| 最后的刺客免费高清国语| 免费人妻精品一区二区三区视频| 国产精品麻豆人妻色哟哟久久| 国产成人av激情在线播放| 国产福利在线免费观看视频| 久久久久久人妻| 搡老乐熟女国产| 亚洲第一区二区三区不卡| 女的被弄到高潮叫床怎么办| 免费av中文字幕在线| 人体艺术视频欧美日本| 亚洲精品456在线播放app| 免费久久久久久久精品成人欧美视频 | 99香蕉大伊视频| 精品视频人人做人人爽| 18禁在线无遮挡免费观看视频| 在线观看免费日韩欧美大片| 国产一区有黄有色的免费视频| 校园人妻丝袜中文字幕| 精品人妻偷拍中文字幕| 精品亚洲成a人片在线观看| 久久久久久久精品精品| 天天躁夜夜躁狠狠久久av| 亚洲精品一区蜜桃| 成人漫画全彩无遮挡| 国产福利在线免费观看视频| 亚洲图色成人| 视频中文字幕在线观看| 国产永久视频网站| 国产亚洲av片在线观看秒播厂| 亚洲 欧美一区二区三区| 男女啪啪激烈高潮av片| 国产高清三级在线| 中文字幕av电影在线播放| 亚洲激情五月婷婷啪啪| 亚洲国产毛片av蜜桃av| 91精品国产国语对白视频| 亚洲精品自拍成人| 午夜影院在线不卡| 国产成人a∨麻豆精品| 成人手机av| 日韩欧美一区视频在线观看| 日韩免费高清中文字幕av| 国产熟女午夜一区二区三区|