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

    基于渦格法的螺旋槳剖面拱線設(shè)計

    2016-09-02 03:08:16饒志強(qiáng)楊晨俊
    中國艦船研究 2016年4期
    關(guān)鍵詞:方法設(shè)計

    饒志強(qiáng),楊晨俊

    1上海交通大學(xué)海洋工程國家重點實驗室,上海2002402高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心,上海200240

    基于渦格法的螺旋槳剖面拱線設(shè)計

    饒志強(qiáng)1,2,楊晨俊1,2

    1上海交通大學(xué)海洋工程國家重點實驗室,上海200240
    2高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心,上海200240

    船舶螺旋槳的葉剖面形狀與其負(fù)荷形式、效率、空泡等水動力性能密切相關(guān)。剖面拱線與剖面的負(fù)荷分布有直接關(guān)系。為了設(shè)計給定負(fù)荷分布的剖面拱線,提出用升力分布定義螺旋槳剖面拱線的負(fù)荷分布,在此基礎(chǔ)上對升力分布進(jìn)行參數(shù)化表達(dá),并進(jìn)行拱線設(shè)計。拱線的設(shè)計采用Newton-Raphson迭代方法使拱線的升力分布滿足給定值,拱線的性能計算采用基于薄翼理論的二維渦格法。計算表明:渦格法與Newton-Raphson迭代方法相結(jié)合能夠精確并快速設(shè)計滿足給定升力分布的拱線。最后,修改升力分布的部分參數(shù),并由此設(shè)計了SJ系列拱線供剖面設(shè)計參考。

    螺旋槳葉剖面;拱線設(shè)計;升力分布;渦格法;SJ系列拱線

    網(wǎng)絡(luò)出版地址:http://www.cnki.net/kcms/detail/42.1755.TJ.20160729.0945.012.html期刊網(wǎng)址:www.ship-research.com

    引用格式:饒志強(qiáng),楊晨俊.基于渦格法的螺旋槳剖面拱線設(shè)計[J].中國艦船研究,2016,11(4):38-43.

    RAO Zhiqiang,YANG Chenjun.Camber lines design approach of propeller blade sections based on the vortex lattice method[J].Chinese Journal of Ship Research,2016,11(4):38-43.

    0 引 言

    剖面形狀對螺旋槳的水動力性能有直接影響。在剖面設(shè)計中,NACA系列的剖面[1]因其優(yōu)良的抗空化性能而被設(shè)計者采用,比如NACA-66(mod)厚度分布和a=0.8拱線組合。但NACA剖面提出的時間很早,之后螺旋槳剖面(或稱翼剖面)設(shè)計也已經(jīng)發(fā)展了幾種方法,這些方法主要有單點設(shè)計法、多點設(shè)計法、面元法和優(yōu)化設(shè)計方法等[2]。單點法以Lighthill[3]方法為代表,做法是給定一個攻角和壓力分布,通過保角變換設(shè)計剖面。多點法以Eppler[4]方法為代表,它是在剖面表面分段給出不同攻角的速度分布、設(shè)計剖面。由于分段處理,它為改善剖面在設(shè)計和非設(shè)計工況下的空泡性能提供了潛力。一些學(xué)者對Eppler方法在剖面設(shè)計中的應(yīng)用進(jìn)行了研究。Eppler等[5]和Shen等[6]將該方法推廣到水翼設(shè)計中。王大政等[7-8]基于Eppler方法提出了可以控制最大厚度和拱度的剖面設(shè)計方法,并研究了剖面參數(shù)對空泡特性的影響。Zhou等[9]將Eppler方法設(shè)計的剖面直接應(yīng)用到螺旋槳上。面元設(shè)計方法是對初始剖面進(jìn)行擾動并做迭代計算,使其滿足給定壓力分布。處理方法主要有2種:一種是直接在剖面上布置控制點,并對其進(jìn)行攝動和迭代計算,這些工作以蘇玉民等[10]和蔡昊鵬等[11]的工作為代表;另一種是用B樣條擬合并擾動擬合曲線控制點,譚廷壽等[12]和李俊華等[13]都采用該方法設(shè)計剖面,所不同的是前者擬合拱線而厚度不變,后者對表面進(jìn)行擬合。前述設(shè)計方法的共同點都是需要按給定壓力或者速度分布進(jìn)行設(shè)計,但在設(shè)計初始階段合理的壓力和速度分布可能并不容易給定。給定值有缺陷可能對幾何或計算壓力分布的光順性造成影響,文獻(xiàn)[11]討論了這個問題。另外,給定壓力分布直接設(shè)計剖面可能會產(chǎn)生非常厚的剖面,文獻(xiàn)[13]中提出該問題并通過減小空泡裕度解決剖面過厚的問題。面對設(shè)計中的問題,優(yōu)化方法成為剖面設(shè)計的另一種選擇。黃勝等[14]和常欣等[15]在剖面上下表面疊加基函數(shù)并采用多目標(biāo)優(yōu)化算法設(shè)計剖面。曾志波[16]和龍文等[17]基于Eppler方法對螺旋槳剖面進(jìn)行了優(yōu)化設(shè)計。在優(yōu)化過程中,文獻(xiàn)[14-15]和文獻(xiàn)[16]分別對基函數(shù)的疊加量或者剖面最大厚度做了約束。

    分析部分針對空泡性能的研究成果,特別是Eppler優(yōu)化設(shè)計結(jié)果[16-17],可以看出滿足約束的最優(yōu)設(shè)計剖面有如下特征:最大厚度位置向?qū)н吳耙?,增加葉背和葉面片空泡裕度,最大拱度位置向隨邊后移將負(fù)荷向隨邊移動[16]。本文作者曾針對固定厚度分布不變,采用面元法設(shè)計了具有不同拱線的剖面,并考察載荷分布對空泡性能的影響[18]。計算表明:載荷的移動對剖面空泡斗的上下位置有明顯的影響。如此可以理解:載荷的移動在一定程度上能夠改變空泡性能。這樣問題的關(guān)鍵變成了如何設(shè)計具有給定載荷的拱線。如能實現(xiàn)給定載荷設(shè)計,拱線的幾何設(shè)計問題實際上也就變成了載荷設(shè)計。

    本文將針對上述問題對載荷分布用升力分布來定義,對其進(jìn)行參數(shù)化表達(dá),并采用Newton-Raphson迭代方法設(shè)計滿足給定升力分布要求的拱線,拱線性能的計算采用二維渦格法。最后,修改給定升力分布的部分參數(shù),設(shè)計SJ系列拱線供參考。

    1 拱線設(shè)計方法

    1.1升力分布的定義

    升力分布(CLx)定義為升力在拱線弦向的密度,具體定義為微段弦長上的升力系數(shù)與微段弦長的比,如式(1)所示。

    式中:ΔCL為對應(yīng)的升力系數(shù);Δc為微段弦長。由于拱線沒有厚度,因此可用薄翼理論來求解。

    在零厚度薄翼理論中,薄翼的升力系數(shù)(CL)為[19]

    式中:L為剖面總升力;ρ為流體密度;V為來流速度;c為弦長;xˉ為無因次的弦向坐標(biāo);γ()xˉ為弦向環(huán)量密度。

    根據(jù)式(1)升力分布的定義,容易得到零厚度拱線的升力分布為

    即在基于薄翼理論的拱線設(shè)計中,給定升力分布可以理解為2倍的給定環(huán)量密度在弦向的分布。采用二維渦格法可數(shù)值求解拱線的環(huán)量密度和升力分布。另外,通過升力分布的定義可知,剖面升力系數(shù)可由升力分布在弦向上的積分得到。

    1.2升力分布的參數(shù)化表達(dá)

    獲得升力分布的定義后,需要給定升力分布的形式。這里采用5段線來給定升力分布的形式。如圖1所示,升力分布由3條直線(L1~L3)和2條相切圓?。–1和C2)構(gòu)成,整條升力分布曲線各段光順連接。

    圖1 升力分布的參數(shù)化表達(dá)Fig.1 Parametric representation of lift distribution

    用5個參數(shù)可確定整個升力分布曲線,包括:

    1)第1個角點A的位置xa;

    2)第2個角點B的位置xb;

    3)圓弧C1的相對半徑,r1=Ar1r1 max;

    4)圓弧C2的相對半徑,r2=Ar2r2 max;

    5)第3段直線L3的斜率角θ3。

    式中:r1 max和r2 max分別為2段圓弧的最大半徑,此時,L1和L2這2條直線長度為0。參數(shù) Ar1和 Ar2為最大半徑的比率,取值在0~1之間,用于控制C1和C2的半徑。5個參數(shù)外加給定升力系數(shù)CL,即可通過迭代的方法獲得整條升力分布曲線。在這種參數(shù)化定義中,拱線導(dǎo)邊和隨邊的升力分布都為0,即載荷為0。第3段線的角度θ3可以控制升力分布的前移或者后移。如果調(diào)整A角和B角點位置,以及圓弧C1和C2的半徑大小還可以對最大升力分布的位置進(jìn)一步進(jìn)行移動。這是一種直接對載荷分布進(jìn)行控制的方法。盡管調(diào)整環(huán)量分布可以達(dá)到載荷調(diào)整的效果,但升力分布的定義更直觀。

    1.3拱線設(shè)計過程

    拱線的設(shè)計以升力分布為目標(biāo)。首先,選擇一個初始拱線,在初始拱線上設(shè)置M+1個控制點,如圖2所示。

    圖2 拱線上控制點布置Fig.2 Control nodes arrangement on camber line

    圖中,x/c和f/c分別為無因次弦長位置和拱度弦長比,導(dǎo)邊和隨邊為固定控制點,其他為可動控制點。對攻角和可動控制點進(jìn)行攝動,以獲得攝動量對弦向某些位置(稱為評價位置)升力分布的影響。本文中評價位置與可動控制點的位置重合。由于增加了攻角的攝動,因此在導(dǎo)邊和第一個可動控制點之間增加一個評估位置來構(gòu)建M階雅可比矩陣,如式(4)所示。

    對新的拱線和攻角再進(jìn)行攝動,直到評價位置上的升力分布計算值和給定值達(dá)到容差為止,本文容差為 1.0×10-4,這個設(shè)計過程即為Newton-Raphson迭代過程。在設(shè)計過程中,攝動初始拱線拱度為0,初始攻角為0.5°。

    拱線的計算方法采用二維渦格法,即將拱線分成若干單元(本文取120個單元)。在每個單元的1/4長度處布置一個點渦,并將3/4位置作為單元控制點的位置,物面不可穿透邊界條件在單元控制點上滿足。渦格法的具體實施可見文獻(xiàn)[19-21]的相關(guān)章節(jié)。

    這里舉一個例子詳細(xì)說明迭代結(jié)果。在這個例子中,θ3=0°,xa=0.1,xb=0.85,Ar1=0.2,Ar2=0.2,設(shè)計升力系數(shù)為CL=1.0。總控制點個數(shù)為31個(對應(yīng)雅可比矩陣階數(shù)M=30),其弦向位置由升力分布的5段線決定。圖3為設(shè)計過程中拱線隨迭代次數(shù)的變化??梢钥闯?,迭代3次后計算結(jié)果收斂。圖4為給定升力分布、設(shè)計拱線的升力分布以及評價位置處升力分布的比較??梢钥闯?,設(shè)計結(jié)果和給定的目標(biāo)十分吻合。這說明本方法迭代計算量少,設(shè)計結(jié)果精確。

    圖3 設(shè)計過程中拱線的變化Fig.3 Camber lines during iterations

    圖4 給定和設(shè)計拱線的升力分布比較Fig.4 Comparison of given and designed lift distributions

    在文獻(xiàn)[18]中也討論過剖面設(shè)計,但剖面的性能計算基于面元法。如果將剖面的厚度限定得很小,攝動剖面的拱線,也能通過Newton-Raphson方法設(shè)計出來。但由于厚度的影響不能完全忽略,實際計算時發(fā)現(xiàn)采用面元法的設(shè)計過程不太穩(wěn)定,設(shè)計拱線也可能會出現(xiàn)不光順的情況。這可能是由于文獻(xiàn)[18]中升力分布的定義和計算是在剖面的上下表面對應(yīng)位置插值造成的,而本文渦格法的計算中升力分布和環(huán)量密度直接相關(guān),無需插值計算。由于渦格法能夠獲得穩(wěn)定、光順的拱線,是一種比較理想的方法,這里僅討論渦格法的計算結(jié)果。

    2 SJ系列拱線

    針對不同的第3段直線斜率角θ3,給定不同的升力分布曲線,設(shè)計了系列拱線,其中,升力分布曲線的其他參數(shù)仍為:xa=0.1,xb=0.85,Ar1=0.2,Ar2=0.2,CL=1.0。這一系列的拱線命名為SJ系列拱線,其中SJ后第1位數(shù)字代表斜率角θ3的正負(fù)(0為負(fù),1為正),后2個數(shù)字代表角度。如SJ015,代表θ3=-15°時的設(shè)計拱線。圖5為部分SJ系列拱線的計算升力分布曲線。由于設(shè)計結(jié)果與給定升力分布十分吻合,故沒有提供給定升力分布的結(jié)果。

    圖5 SJ系列拱線的計算升力分布Fig.5 Calculated lift distributions of SJ series camber lines

    圖6為對應(yīng)的部分SJ系列的無因次拱線圖,其中f/fmax為弦向拱度與最大拱度比值。從圖中可以看出,拱線的最大拱度位置隨升力分布的前移或后移而分別前后移動。

    圖6 SJ系列無因次拱線圖Fig.6 Non-dimensional SJ series camber lines

    表1為SJ系列拱線的攻角α和拱度弦長比fmax/c。表2為SJ系列拱線的型值表,可供后續(xù)的設(shè)計參考。

    表1 SJ系列拱線攻角和拱度弦長比(CL=1.0)Tab.1 Angle of attack and camber ratio of SJ series camber lines(CL=1.0)

    表2 SJ系列拱線型值表Tab.2 Offset table of SJ series camber lines

    3 結(jié) 論

    本文對螺旋槳剖面拱線的升力分布進(jìn)行了定義,并提出了一種升力分布的參數(shù)化方法。采用Newton-Raphson方法設(shè)計了給定升力分布的系列拱線。主要結(jié)論和工作如下:

    1)對螺旋槳剖面拱線的載荷用升力分布來表達(dá),并提出了一種升力分布的參數(shù)化形式。

    2)結(jié)合Newton-Raphson迭代方法和渦格法可以快速、精確地設(shè)計給定升力分布的拱線。

    3)設(shè)計了SJ系列拱線并給出了幾何型值供設(shè)計參考。

    [1]ABBOTT I H,VON DOENHOFF A E.Theory of wingsections:including a summary of airfoil data[M].Dover,DE:Dover Publications Inc.,1959.

    [2]王獻(xiàn)孚.船用翼理論[M].北京:國防工業(yè)出版社,1998.

    [3]LIGHTHILL M J.A new method of two-dimensional aerodynamic design:M2112[R].London,UK:Aeronautical Research Council,1945.

    [4]EPPLER R,SOMERS D M.A computer program for thedesignandanalysisoflowspeedairfoils:NASA-TM-80210[R].Hampton,VA:NASA Langley Research Center,1980.

    [5]EPPLER R,SHEN Y T.Wing sections for hydrofoils-part 1:symmetrical profiles[J].Journal of Ship Research,1979,23(3):209-217.

    [6]SHEN Y T,EPPLER R.Wing sections for hydrofoils-part 2:nonsymmetrical profiles[J].Journal of Ship Research,1981,25(3):191-200.

    [7] 王大政,王言英,柴樹紅.新型葉剖面設(shè)計及葉剖面參數(shù)對空泡特性影響的研究——(1)新型葉剖面設(shè)計[J].水動力學(xué)研究與進(jìn)展,1999,14(2):189-200.

    WANG Dazheng,WANG Yanying,CHAI Shuhong. Study on a new approach to blade section design and effects of its parameters on cavitation characteristics partⅠ:a new approach to blade section design[J].Journal of Hydrodynamics,1999,14(2):189-200.

    [8]王大政,王言英.新型葉剖面設(shè)計及葉剖面參數(shù)對空泡特性影響的研究—(2)剖面參數(shù)對空泡特性的影響[J].水動力學(xué)研究與進(jìn)展,2000,15(3):320-328.

    WANG Dazheng,WANG Yanying.Study on a new approach to blade section design and effects of its parameters on cavitation characteristics partⅡ:effects of blade section parameters on cavitation characteristics [J].Journal of Hydrodynamics,2000,15(3):320-328.

    [9] ZHOU W X,WU Y H,DONG S T.A propeller design method with new blade section for improving cavitation inception under unsteady condition[J].Journal of Ship Mechanics,2001,5(6):18-26.

    [10]SU Y M,IKEHATA M,KAI H.A numerical method for designing three-dimensional wing based on surface panel method[J].Journal of the Society of Naval Architects of Japan,1997,182:39-47.

    [11]蔡昊鵬,蘇玉民,沈海龍.基于面元法的三維機(jī)翼數(shù)值設(shè)計[J].海軍工程大學(xué)學(xué)報,2009,21(6):30-35.

    CAI Haopeng,SU Yumin,SHEN Hailong.Numerical design of three-dimensional wing based on surface panel method[J].Journal of Naval University of Engineering,2009,21(6):30-35.

    [12]譚廷壽,劉飛.非均勻流場中槳葉剖面設(shè)計方法[J].武漢理工大學(xué)學(xué)報(交通科學(xué)與工程版),2006,30(3):437-440.

    TAN Tingshou,LIU Fei.Design method of blade section in non-uniform flow field[J].Journal of Wuhan University of Technology(Transportation Science& Engineering Edition),2006,30(3):437-440.

    [13]李俊華,唐登海,董世湯.指定壓力分布的螺旋槳設(shè)計[J].船舶力學(xué),2010,14(1/2):10-19.

    LI Junhua,TANG Denghai,DONG Shitang.Propeller design by prescribed pressure distribution[J]. Journal of Ship Mechanics,2010,14(1/2):10-19.

    [14]黃勝,任萬龍,王超.水翼剖面多目標(biāo)粒子群算法優(yōu)化[J].哈爾濱工程大學(xué)學(xué)報,2014,35(12):1451-1457.

    HUANGSheng,RENWanlong,WANGChao. Multi-objective particle swarm optimization of hydrofoil sections[J].Journal of Harbin Engineering University,2014,35(12):1451-1457.

    [15]常欣,郭春雨,孟祥印,等.基于粒子群算法的水翼剖面優(yōu)化設(shè)計[J].船舶工程,2010,32(5):1-3,26.

    CHANG Xin,GUO Chunyu,Meng Xiangyin,et al. Design of hydrofoil section based on particle swarm optimization[J].Ship Engineering,2010,32(5):1-3,26.

    [16] 曾志波.最大空泡起始航速及最小空泡脈動壓力螺旋槳設(shè)計[D].北京:中國艦船研究院,2014.

    [17]龍文,謝偉,楊向暉,等.基于新型剖面設(shè)計的螺旋槳多目標(biāo)性能權(quán)衡優(yōu)化[J].艦船科學(xué)技術(shù),2011,33(11):56-60.

    LONG Wen,XIE Wei,YANG Xianghui,et al. Multi-objective optimization of open water performance of propeller based on new anti-cavitation airfoil profile design[J].Ship Science and Technology,2011,33(11):56-60.

    [18] 饒志強(qiáng),楊晨俊,李巍.基于升力分布的二維翼型設(shè)計方法[C]//2015年船舶水動力學(xué)學(xué)術(shù)會議.哈爾濱:[s.n.],2015.

    [19] 王國強(qiáng),董世湯.船舶螺旋槳理論與應(yīng)用[M].哈爾濱:哈爾濱工程大學(xué)出版社,2007.

    [20] 徐華舫,張炳暄,朱自強(qiáng),等.亞、超音速定常位流的面元法[M].北京:國防工業(yè)出版社,1981.

    [21]KATZ J,PLOTKIN A.Low-speed aerodynamics [M].2nd ed.New York:Cambridge University Press,2001.

    Camber lines design approach of propeller blade sections based on the vortex lattice method

    RAO Zhiqiang1,2,YANG Chenjun1,2
    1 State Key Laboratory of Ocean Engineering,Shanghai Jiao Tong University,Shanghai 200240,China
    2 Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration,Shanghai 200240,China

    Blade sections correlates closely with marine propeller hydrodynamic performance,including loading status,efficiency,cavitation,and other factors,and there is a direct relationship between camber lines and the loading status of blade sections.For the sake of designing camber lines of blade sections with specified loading,a lift distribution is first defined as loading description.A parametric method is then proposed for lift distribution expression,after which the camber line is designed for the specified lift distribution through Newton-Raphson iteration.Finally,the flow characteristics of camber lines are calculated with 2D vortex lattice method.The results indicate that the lift distribution of camber lines designed by combining Newton-Raphson iteration and vortex lattice method matches the given one.Furthermore,certain parameters of the lift distribution are modified to generate camber lines called SJ series for the future design.

    propeller blade section;camber line design;lift distribution;vortex lattice method;SJ series camber lines

    U664.33

    A

    10.3969/j.issn.1673-3185.2016.04.006

    2015-09-30網(wǎng)絡(luò)出版時間:2016-7-29 9:45

    饒志強(qiáng)(通信作者),男,1986年生,博士生。研究方向:船舶推進(jìn)器水動力學(xué)。

    E-mail:hellostar@126.com

    楊晨俊,男,1964年生,博士,教授,博士生導(dǎo)師。研究方向:船舶推進(jìn)器水動力學(xué)。

    E-mail:cjyang@sjtu.edu.cn

    猜你喜歡
    方法設(shè)計
    何為設(shè)計的守護(hù)之道?
    《豐收的喜悅展示設(shè)計》
    流行色(2020年1期)2020-04-28 11:16:38
    學(xué)習(xí)方法
    瞞天過海——仿生設(shè)計萌到家
    設(shè)計秀
    海峽姐妹(2017年7期)2017-07-31 19:08:17
    有種設(shè)計叫而專
    Coco薇(2017年5期)2017-06-05 08:53:16
    可能是方法不對
    用對方法才能瘦
    Coco薇(2016年2期)2016-03-22 02:42:52
    四大方法 教你不再“坐以待病”!
    Coco薇(2015年1期)2015-08-13 02:47:34
    賺錢方法
    全区人妻精品视频| 男女那种视频在线观看| 久久久久久久精品吃奶| 神马国产精品三级电影在线观看| 亚洲在线观看片| 母亲3免费完整高清在线观看| 一a级毛片在线观看| 超碰成人久久| 国产成人欧美在线观看| 免费观看人在逋| 无遮挡黄片免费观看| 99国产精品99久久久久| 热99re8久久精品国产| 视频区欧美日本亚洲| 欧美+亚洲+日韩+国产| 天堂av国产一区二区熟女人妻| 女警被强在线播放| 久久香蕉国产精品| 国内少妇人妻偷人精品xxx网站 | 亚洲av成人不卡在线观看播放网| 怎么达到女性高潮| 成人av在线播放网站| 国产精品日韩av在线免费观看| 午夜福利在线在线| 久久天堂一区二区三区四区| 成年女人永久免费观看视频| 久久久国产精品麻豆| 90打野战视频偷拍视频| 国产私拍福利视频在线观看| 中文字幕高清在线视频| 成人av在线播放网站| 日韩 欧美 亚洲 中文字幕| 啦啦啦观看免费观看视频高清| 久久香蕉国产精品| 最新美女视频免费是黄的| 精品久久久久久久毛片微露脸| 亚洲最大成人中文| 亚洲人成电影免费在线| 久久久久久九九精品二区国产| 蜜桃久久精品国产亚洲av| 哪里可以看免费的av片| 超碰成人久久| 国产成人aa在线观看| 一进一出好大好爽视频| 欧美绝顶高潮抽搐喷水| 亚洲av电影在线进入| 精品国产亚洲在线| 91麻豆av在线| 亚洲欧美日韩无卡精品| 亚洲精品一区av在线观看| 亚洲国产高清在线一区二区三| 欧美极品一区二区三区四区| 亚洲熟妇熟女久久| 免费无遮挡裸体视频| 在线免费观看不下载黄p国产 | 五月伊人婷婷丁香| 丁香六月欧美| 在线观看午夜福利视频| 又紧又爽又黄一区二区| 特大巨黑吊av在线直播| 国产精品精品国产色婷婷| 国产亚洲av高清不卡| 色av中文字幕| 国产亚洲欧美在线一区二区| 国产成人精品久久二区二区免费| 一本一本综合久久| 国产免费av片在线观看野外av| 国产午夜精品论理片| 久久性视频一级片| 男女做爰动态图高潮gif福利片| 免费在线观看视频国产中文字幕亚洲| 一本一本综合久久| 日本精品一区二区三区蜜桃| 黑人巨大精品欧美一区二区mp4| 美女扒开内裤让男人捅视频| 精品国产乱码久久久久久男人| 一进一出抽搐gif免费好疼| 日本在线视频免费播放| 精品久久久久久久末码| 亚洲乱码一区二区免费版| 一级毛片高清免费大全| 好男人在线观看高清免费视频| 波多野结衣高清无吗| 一级毛片精品| 国产淫片久久久久久久久 | 黑人操中国人逼视频| 真人一进一出gif抽搐免费| 久久婷婷人人爽人人干人人爱| 搡老熟女国产l中国老女人| 操出白浆在线播放| 十八禁网站免费在线| 女同久久另类99精品国产91| 国产精品美女特级片免费视频播放器 | 亚洲一区高清亚洲精品| 综合色av麻豆| 我的老师免费观看完整版| 亚洲一区高清亚洲精品| av在线天堂中文字幕| 老汉色∧v一级毛片| 成年女人看的毛片在线观看| 久久精品亚洲精品国产色婷小说| 少妇的丰满在线观看| 91在线精品国自产拍蜜月 | 精品日产1卡2卡| 日本免费a在线| 久久久久亚洲av毛片大全| 国产成人精品久久二区二区免费| 亚洲人与动物交配视频| 免费高清视频大片| 999精品在线视频| 99热精品在线国产| 1024手机看黄色片| 国产成年人精品一区二区| 久久香蕉精品热| 久久国产精品人妻蜜桃| 给我免费播放毛片高清在线观看| 日本精品一区二区三区蜜桃| 亚洲精品在线美女| 蜜桃久久精品国产亚洲av| 天天躁日日操中文字幕| 午夜福利18| 国产午夜精品论理片| 首页视频小说图片口味搜索| 国产精品爽爽va在线观看网站| 日本五十路高清| 国产亚洲精品av在线| av片东京热男人的天堂| 黄频高清免费视频| 日韩欧美国产一区二区入口| 中文字幕高清在线视频| 国内精品久久久久久久电影| 日韩欧美一区二区三区在线观看| 在线观看66精品国产| 国产高清三级在线| 欧美成人性av电影在线观看| 国产精品综合久久久久久久免费| 国产1区2区3区精品| e午夜精品久久久久久久| 国内毛片毛片毛片毛片毛片| 99国产极品粉嫩在线观看| 国产熟女xx| 欧美另类亚洲清纯唯美| 欧美在线一区亚洲| 一区福利在线观看| 欧美乱色亚洲激情| 欧美极品一区二区三区四区| 一个人观看的视频www高清免费观看 | 亚洲午夜理论影院| 久久这里只有精品19| 最近视频中文字幕2019在线8| 久久久国产欧美日韩av| 亚洲国产欧美网| 国产亚洲欧美98| 亚洲精品中文字幕一二三四区| 身体一侧抽搐| 免费在线观看日本一区| 亚洲人成伊人成综合网2020| 一夜夜www| 国产精品久久久人人做人人爽| 午夜福利在线观看免费完整高清在 | 日本三级黄在线观看| 99riav亚洲国产免费| 深夜精品福利| 此物有八面人人有两片| 成年版毛片免费区| 亚洲天堂国产精品一区在线| 久9热在线精品视频| 成人特级av手机在线观看| 亚洲美女黄片视频| 免费无遮挡裸体视频| 亚洲欧美日韩高清在线视频| 久久久久久久精品吃奶| 国产精品久久久人人做人人爽| 成人av一区二区三区在线看| 在线观看舔阴道视频| 亚洲va日本ⅴa欧美va伊人久久| 亚洲va日本ⅴa欧美va伊人久久| 久久久久免费精品人妻一区二区| 日韩欧美在线二视频| 香蕉国产在线看| 欧美一区二区国产精品久久精品| 国产野战对白在线观看| 岛国在线观看网站| 国产精品 欧美亚洲| 亚洲色图 男人天堂 中文字幕| 亚洲av电影在线进入| 午夜影院日韩av| 亚洲熟女毛片儿| 在线观看舔阴道视频| 蜜桃久久精品国产亚洲av| 别揉我奶头~嗯~啊~动态视频| 18禁美女被吸乳视频| www.自偷自拍.com| 欧美不卡视频在线免费观看| 啦啦啦观看免费观看视频高清| 精品久久蜜臀av无| 毛片女人毛片| 叶爱在线成人免费视频播放| 国产91精品成人一区二区三区| 天堂影院成人在线观看| 老司机午夜十八禁免费视频| 色精品久久人妻99蜜桃| 久久久国产精品麻豆| 久久99热这里只有精品18| 亚洲欧美精品综合一区二区三区| 精品久久久久久久人妻蜜臀av| 精品久久蜜臀av无| 好男人在线观看高清免费视频| 欧美性猛交╳xxx乱大交人| 国产乱人视频| 精品福利观看| 9191精品国产免费久久| e午夜精品久久久久久久| 男人的好看免费观看在线视频| 久久香蕉国产精品| 亚洲 国产 在线| 激情在线观看视频在线高清| 成人一区二区视频在线观看| 成人亚洲精品av一区二区| 免费大片18禁| 中文字幕人妻丝袜一区二区| 国产一区二区在线观看日韩 | 亚洲av成人不卡在线观看播放网| 在线a可以看的网站| 黑人操中国人逼视频| 国产亚洲欧美98| 中国美女看黄片| 久久中文字幕人妻熟女| 国产91精品成人一区二区三区| av中文乱码字幕在线| 亚洲中文字幕一区二区三区有码在线看 | 法律面前人人平等表现在哪些方面| 热99re8久久精品国产| 性色avwww在线观看| 在线观看美女被高潮喷水网站 | 99热精品在线国产| 蜜桃久久精品国产亚洲av| 欧美丝袜亚洲另类 | 国产激情偷乱视频一区二区| 又爽又黄无遮挡网站| 成人av在线播放网站| 日本 欧美在线| 中文字幕熟女人妻在线| 精品乱码久久久久久99久播| 一进一出抽搐gif免费好疼| 国产精品1区2区在线观看.| 99久久99久久久精品蜜桃| 成人欧美大片| 中文在线观看免费www的网站| 久久久久国产一级毛片高清牌| 午夜激情欧美在线| 色播亚洲综合网| 黄色视频,在线免费观看| 国产高清视频在线观看网站| 美女午夜性视频免费| 神马国产精品三级电影在线观看| 国产精品久久久久久久电影 | 舔av片在线| 精品乱码久久久久久99久播| 成年版毛片免费区| 久久久精品大字幕| 亚洲国产精品久久男人天堂| xxxwww97欧美| 午夜激情欧美在线| 精品一区二区三区视频在线 | 国产亚洲精品一区二区www| 欧美一级a爱片免费观看看| 国产av一区在线观看免费| 日本一二三区视频观看| 18禁美女被吸乳视频| 精品国产乱码久久久久久男人| 午夜福利免费观看在线| 国产欧美日韩精品亚洲av| 婷婷六月久久综合丁香| 亚洲av成人不卡在线观看播放网| 我的老师免费观看完整版| 全区人妻精品视频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲欧美日韩东京热| 大型黄色视频在线免费观看| 精品久久久久久久久久免费视频| 超碰成人久久| 波多野结衣高清无吗| 18禁国产床啪视频网站| 午夜亚洲福利在线播放| 麻豆成人午夜福利视频| 法律面前人人平等表现在哪些方面| 免费在线观看视频国产中文字幕亚洲| 两性午夜刺激爽爽歪歪视频在线观看| 日本成人三级电影网站| 美女大奶头视频| 国产欧美日韩精品亚洲av| 母亲3免费完整高清在线观看| e午夜精品久久久久久久| 久久久久精品国产欧美久久久| 亚洲av片天天在线观看| 国产精品98久久久久久宅男小说| av在线蜜桃| 日韩欧美一区二区三区在线观看| cao死你这个sao货| 十八禁人妻一区二区| 麻豆成人午夜福利视频| 午夜激情欧美在线| 国产v大片淫在线免费观看| 国产 一区 欧美 日韩| 一级黄色大片毛片| 亚洲自拍偷在线| 天堂动漫精品| 日日干狠狠操夜夜爽| 欧美日韩福利视频一区二区| 欧美成人一区二区免费高清观看 | 亚洲性夜色夜夜综合| 国产精品女同一区二区软件 | 色综合亚洲欧美另类图片| 国产精品免费一区二区三区在线| 精品国产乱码久久久久久男人| 99精品欧美一区二区三区四区| 琪琪午夜伦伦电影理论片6080| 亚洲五月天丁香| 久久香蕉国产精品| 一夜夜www| 真人一进一出gif抽搐免费| 国产单亲对白刺激| 伦理电影免费视频| 亚洲在线自拍视频| 神马国产精品三级电影在线观看| 国产黄a三级三级三级人| 桃红色精品国产亚洲av| 丰满人妻熟妇乱又伦精品不卡| 日本黄大片高清| 岛国视频午夜一区免费看| 国模一区二区三区四区视频 | 一进一出抽搐gif免费好疼| 久久久成人免费电影| 九色成人免费人妻av| 男女之事视频高清在线观看| 青草久久国产| e午夜精品久久久久久久| 中文在线观看免费www的网站| 91老司机精品| 99热精品在线国产| 国产伦一二天堂av在线观看| 黑人操中国人逼视频| 国产黄片美女视频| 12—13女人毛片做爰片一| 国产又色又爽无遮挡免费看| 亚洲av成人av| 国产高清有码在线观看视频| 高清毛片免费观看视频网站| 久久草成人影院| 中文字幕久久专区| 国产午夜精品论理片| 国产视频一区二区在线看| 国产精品一区二区三区四区免费观看 | 日本与韩国留学比较| 少妇丰满av| 日韩三级视频一区二区三区| cao死你这个sao货| 麻豆国产av国片精品| 1024手机看黄色片| 日韩av在线大香蕉| 最近在线观看免费完整版| 日本黄色片子视频| 久久精品亚洲精品国产色婷小说| 国产激情偷乱视频一区二区| 岛国在线观看网站| 精品电影一区二区在线| 88av欧美| 国产精品久久久人人做人人爽| 极品教师在线免费播放| 午夜激情欧美在线| 久久伊人香网站| 韩国av一区二区三区四区| 国产精品久久久久久亚洲av鲁大| 国产又色又爽无遮挡免费看| 麻豆国产97在线/欧美| 男女之事视频高清在线观看| 日韩有码中文字幕| 国产亚洲欧美98| 日韩大尺度精品在线看网址| 国内精品久久久久精免费| 成人国产综合亚洲| 伊人久久大香线蕉亚洲五| 人妻夜夜爽99麻豆av| 亚洲美女视频黄频| 在线免费观看的www视频| 免费观看人在逋| 成人三级做爰电影| 99国产精品一区二区三区| 动漫黄色视频在线观看| 色综合亚洲欧美另类图片| or卡值多少钱| 成人无遮挡网站| 国产成人欧美在线观看| 欧美国产日韩亚洲一区| 日本熟妇午夜| 男女下面进入的视频免费午夜| 国产成人精品久久二区二区91| 一个人观看的视频www高清免费观看 | 国产伦人伦偷精品视频| 久久这里只有精品19| 18禁黄网站禁片午夜丰满| 熟女人妻精品中文字幕| 最近最新中文字幕大全电影3| 国产毛片a区久久久久| 国产免费av片在线观看野外av| 亚洲,欧美精品.| 性色av乱码一区二区三区2| 一边摸一边抽搐一进一小说| 久久久国产成人精品二区| 两个人看的免费小视频| 亚洲最大成人中文| 在线视频色国产色| 亚洲精品粉嫩美女一区| 欧美成人一区二区免费高清观看 | 欧美丝袜亚洲另类 | 成人鲁丝片一二三区免费| 看免费av毛片| 国产精品一区二区精品视频观看| 日韩人妻高清精品专区| 黄片小视频在线播放| 69av精品久久久久久| 美女cb高潮喷水在线观看 | 美女 人体艺术 gogo| 亚洲成a人片在线一区二区| 免费一级毛片在线播放高清视频| 99视频精品全部免费 在线 | 在线观看免费午夜福利视频| 国产69精品久久久久777片 | 国产视频内射| 欧美日韩福利视频一区二区| 18禁观看日本| 午夜免费激情av| 欧美乱码精品一区二区三区| 亚洲性夜色夜夜综合| 国产av麻豆久久久久久久| 日韩三级视频一区二区三区| 亚洲av片天天在线观看| 精品免费久久久久久久清纯| 老司机深夜福利视频在线观看| 黄片小视频在线播放| 亚洲人成电影免费在线| 曰老女人黄片| 国产黄a三级三级三级人| 国产97色在线日韩免费| 国产亚洲av高清不卡| 国产精品 国内视频| 精品日产1卡2卡| 国产真人三级小视频在线观看| 亚洲男人的天堂狠狠| 美女午夜性视频免费| 在线观看免费午夜福利视频| 国产毛片a区久久久久| 成年女人永久免费观看视频| 99久久成人亚洲精品观看| 色哟哟哟哟哟哟| 色播亚洲综合网| 亚洲中文日韩欧美视频| 国产精品,欧美在线| 色综合欧美亚洲国产小说| 99久久久亚洲精品蜜臀av| 日日干狠狠操夜夜爽| 小蜜桃在线观看免费完整版高清| 久久精品国产亚洲av香蕉五月| 精品99又大又爽又粗少妇毛片 | 黑人欧美特级aaaaaa片| 国产亚洲欧美在线一区二区| 超碰成人久久| 久久久久国内视频| ponron亚洲| 两人在一起打扑克的视频| 欧美在线黄色| 真实男女啪啪啪动态图| 久久精品综合一区二区三区| 免费在线观看视频国产中文字幕亚洲| 精品久久久久久久末码| 黄色 视频免费看| 欧美色欧美亚洲另类二区| 五月伊人婷婷丁香| 国产高清三级在线| 亚洲无线在线观看| www日本在线高清视频| 黄片大片在线免费观看| 日韩精品青青久久久久久| 亚洲欧美精品综合久久99| 一进一出好大好爽视频| 日日夜夜操网爽| 国产精品av久久久久免费| 激情在线观看视频在线高清| 欧美乱妇无乱码| 一区二区三区高清视频在线| 好男人在线观看高清免费视频| 亚洲成人免费电影在线观看| 国产成人aa在线观看| 18禁国产床啪视频网站| 亚洲国产中文字幕在线视频| 男女床上黄色一级片免费看| 亚洲国产日韩欧美精品在线观看 | 亚洲精品456在线播放app | 亚洲av成人一区二区三| 亚洲精品中文字幕一二三四区| 又黄又爽又免费观看的视频| 好男人在线观看高清免费视频| 村上凉子中文字幕在线| 高清在线国产一区| 国产三级黄色录像| 免费在线观看视频国产中文字幕亚洲| 中国美女看黄片| 国产精品综合久久久久久久免费| 国产成人欧美在线观看| 国产私拍福利视频在线观看| 日韩欧美在线乱码| 最近最新免费中文字幕在线| 伊人久久大香线蕉亚洲五| 一卡2卡三卡四卡精品乱码亚洲| 香蕉av资源在线| 看片在线看免费视频| 亚洲欧美一区二区三区黑人| 久久精品91蜜桃| 国产视频内射| 18禁国产床啪视频网站| 91九色精品人成在线观看| 国产精品 国内视频| 亚洲 欧美 日韩 在线 免费| 夜夜夜夜夜久久久久| 久久久成人免费电影| 亚洲中文日韩欧美视频| 中亚洲国语对白在线视频| 一区二区三区高清视频在线| 亚洲av电影不卡..在线观看| 哪里可以看免费的av片| 亚洲午夜精品一区,二区,三区| 亚洲色图 男人天堂 中文字幕| 国产野战对白在线观看| 日韩欧美国产在线观看| 精品久久久久久久毛片微露脸| 亚洲专区国产一区二区| 国产麻豆成人av免费视频| 男人舔女人下体高潮全视频| 亚洲 欧美一区二区三区| 岛国在线免费视频观看| 老熟妇仑乱视频hdxx| 欧美成人免费av一区二区三区| 久久精品影院6| 久久久久国产精品人妻aⅴ院| 真人一进一出gif抽搐免费| 无遮挡黄片免费观看| 男女那种视频在线观看| av在线天堂中文字幕| 日本黄色视频三级网站网址| 免费在线观看亚洲国产| 1024香蕉在线观看| 韩国av一区二区三区四区| 一边摸一边抽搐一进一小说| 免费看a级黄色片| av黄色大香蕉| 久久中文字幕一级| 日本黄大片高清| 久久草成人影院| 日本精品一区二区三区蜜桃| 免费看光身美女| 精品一区二区三区视频在线 | 国产精品久久久久久人妻精品电影| 人妻夜夜爽99麻豆av| 亚洲欧美精品综合久久99| 亚洲精品粉嫩美女一区| av福利片在线观看| 亚洲欧美激情综合另类| 国产精品一及| 特级一级黄色大片| 久久午夜亚洲精品久久| 精华霜和精华液先用哪个| 国产探花在线观看一区二区| 亚洲欧美日韩卡通动漫| 91字幕亚洲| 成人18禁在线播放| 嫩草影院精品99| 成人特级黄色片久久久久久久| 性欧美人与动物交配| 国产高清视频在线观看网站| 亚洲最大成人中文| 国产成人av教育| 18禁国产床啪视频网站| 免费看光身美女| 亚洲成人久久爱视频| 香蕉av资源在线| 老司机深夜福利视频在线观看| 别揉我奶头~嗯~啊~动态视频| avwww免费| 一级毛片高清免费大全| 18禁裸乳无遮挡免费网站照片| 国产野战对白在线观看| 啦啦啦韩国在线观看视频| 久久天躁狠狠躁夜夜2o2o| 亚洲自拍偷在线| 日韩大尺度精品在线看网址| 欧美激情在线99| 婷婷精品国产亚洲av在线| 最近最新中文字幕大全电影3| 91在线精品国自产拍蜜月 | 天天添夜夜摸| 动漫黄色视频在线观看| 成人亚洲精品av一区二区| 欧美性猛交黑人性爽| 九九在线视频观看精品| 亚洲一区高清亚洲精品| 波多野结衣高清作品| 欧美不卡视频在线免费观看| 老司机午夜十八禁免费视频| 天堂动漫精品| 免费无遮挡裸体视频| 国内精品久久久久久久电影| 亚洲av美国av| 手机成人av网站| 热99re8久久精品国产|