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

    基于泰勒展開邊界元法的近水面潛艇垂向二階波浪力(矩)計算

    2017-02-17 09:36:32段文洋王隸加陳紀康趙彬彬
    關(guān)鍵詞:橢球元法泰勒

    段文洋, 王隸加, 陳紀康, 趙彬彬

    (哈爾濱工程大學(xué) 船舶工程學(xué)院,黑龍江 哈爾濱 150001)

    基于泰勒展開邊界元法的近水面潛艇垂向二階波浪力(矩)計算

    段文洋, 王隸加, 陳紀康, 趙彬彬

    (哈爾濱工程大學(xué) 船舶工程學(xué)院,黑龍江 哈爾濱 150001)

    針對潛艇在近水面航行難以定深控制問題,利用一階泰勒展開邊界元方法求解切向誘導(dǎo)速度精度高的優(yōu)勢對其垂向二階波浪力(矩)研究。該方法對邊界積分方程中的偶極強度進行泰勒展開并保留至一階導(dǎo)數(shù)項,同時在格林第三公式中關(guān)于場點沿邊界取切向?qū)?shù)封閉方程組,直接求解出速度勢及其沿物面的切向速度。計算迎浪狀態(tài)下橢球體的垂向二階波浪力(矩),與現(xiàn)有成果吻合度較高;進而計算迎浪狀態(tài)下潛艇模型的垂向二階波浪力(矩)。數(shù)值計算結(jié)果表明:該方法有較高的計算精度且收斂速度快。

    泰勒展開邊界元法;垂向二階波浪載荷;潛艇;頻域;切向誘導(dǎo)速度

    潛艇在近水面航行,有時需要保持定深以完成指定任務(wù)。文獻[1]表明:垂向平面內(nèi)的二階波浪力(矩)會使?jié)撏е饾u向波面上浮,不利于定深控制和操縱。特別是在迎浪狀態(tài)下,二階波浪力表現(xiàn)為升力,將潛艇推向水面,甚至出現(xiàn)“拋甩”現(xiàn)象;此外還使?jié)撏Мa(chǎn)生縱傾,甚至艉部出水,對螺旋槳的性能造成影響。因此,有必要對近水面潛艇的垂向二階波浪力(矩)展開研究。

    繆國平等[2]用切片法將二階波浪力從二維拓展到三維,計算潛體在不同潛深和浪向角時二階垂向波浪力的變化趨勢。林青山等[3]利用二階波浪力切片理論對潛艇二階縱傾力矩等進行了考察,并分析了指揮臺、艉舵等對整條潛艇的影響。馮學(xué)知等[4]利用STF流體動力切片理論和Frank源相結(jié)合的方法對近水面潛體在波浪力作用下的運動響應(yīng)建模并計算,表明二階波浪力對潛體在垂向平面內(nèi)的上浮起主要作用。馮學(xué)知等[5]基于細長體假設(shè),結(jié)合三維修正系數(shù),計算不同潛深和浪向下的運動響應(yīng)和二階波浪力(矩)。雖然切片法廣泛應(yīng)用于潛艇受力和運動預(yù)報,但本質(zhì)上仍基于二維理論,對艇體首尾帶來的影響值得商榷;而通過三維修正系數(shù)雖然能夠滿足工程要求,但該方法對不同模型的適用性需要驗證。

    Lee等[6]利用三維面元法對迎浪狀態(tài)下橢球體的一階波浪力(矩)和二階波浪力(矩)進行數(shù)值計算,并與細長體理論比較。結(jié)果表明:隨著面元數(shù)量增加,二階垂向波浪力(矩)計算收斂極慢。雖然利用三維分布源頻域方法可以求解二階波浪力,但基于低階分布源面元法求解二階波浪力收斂較差,尤其是對非光滑邊界物體,切向誘導(dǎo)速度誤差極大[7]。針對該困難,段文洋[8]提出了泰勒展開邊界元方法(taylor expansion boundary element method,TEBEM)的思想,即對源-偶混合分布方法中偶極強度進行一階泰勒展開。其中,偶極切向變化恰好是切向速度,該方法對尖角處誘導(dǎo)速度有很好的模擬。Duan等[9-11]又將該方法展開至二階,處理繞射輻射等問題,證明了泰勒展開邊界元法能夠有效地提高非光滑邊界處切向誘導(dǎo)速度的計算精度。

    本文基于一階泰勒展開邊界元方法,對近水面潛艇二階定常波浪力(矩)進行計算,考慮該方法對潛艇計算的適應(yīng)性。

    1 近水面潛體二階波浪力的定解問題

    通常將空間速度勢分解為入射勢和擾動勢兩部分,即φ=φ0+φP。其中,入射速度勢的表達式為

    (1)

    式中:A、ω、k0分別表示入射波浪的波幅、圓頻率、波數(shù),β表示入射波的傳播方向與x軸正向之間的夾角(浪向角),g表示重力加速度。

    根據(jù)疊加原理,擾動勢可分解為7個組成部分:

    (2)

    式中:vj為潛體6個自由度的運動速度,φj表示由于物體單位運動引起的輻射勢,φ7為繞射勢。

    擾動速度勢φj(j=1,2,…,7)的定解條件為

    (3)

    式中:D表示流域;n表示物面法向,由流體域指向潛體內(nèi)部;SH表示潛體平均濕表面。

    根據(jù)所求得的速度勢和速度,利用近場公式沿潛體表面進行積分,可以得到作用在潛體上的一階和二階平均波浪力(矩):

    (4)

    (5)

    (6)

    (7)

    2 一階泰勒展開邊界元方法

    傳統(tǒng)基于源-偶混合分布的低階面元法以各網(wǎng)格面元中點速度勢作為整個面元速度勢。利用無限水深三維頻域無航速格林函數(shù)可以得到潛體平均濕表面上的速度勢邊界積分方程[1],即:

    (8)

    (9)

    為了構(gòu)建封閉方程組,對一階TEBEM方法進行補充。在面元所在的局部坐標系下作垂直于法向的平面,并在該平面內(nèi)取兩個相互正交方向?qū)鳇c求一階導(dǎo)數(shù)。得到以下方程組:

    (10)

    將式(9)代入到式(8)、(10)中聯(lián)立,并將積分邊界SH離散為N個面元,可以得到一階TEBEM方法的離散格式。對于任意面元i可以得到如下離散方程組,i=1,2,…,N:

    (11)

    式中:i和j表示面元編號;k= 1, 2, 3;m=1,2, 3。詳細推導(dǎo)見文獻[9-11]。

    3 水下橢球體的結(jié)果與比較

    為了驗證TEBEM方法對于求解二階波浪力(矩)的準確性,本文與Lee等[6]利用常數(shù)元方法計算的水下橢球體數(shù)值結(jié)果進行對比。

    該算例假定長軸L與短軸D之比L/D=10的零航速橢球回轉(zhuǎn)體在迎浪狀態(tài)下重心位于靜水面(0,0,-D)處,計算不同波浪頻率下受到的波浪力(矩)。

    圖1~4分別表示利用一階TEBEM方法在不同頻率下計算的一階垂蕩波浪力、一階縱搖波浪力矩、二階垂向波浪力和二階縱傾波浪力矩。

    圖1 橢球體一階垂蕩波浪力Fig.1 First-order heave force of submerged spheroid

    圖2 橢球體一階縱搖波浪力矩Fig.2 First-order pitch moment of submerged spheroid

    圖3 橢球體二階垂向波浪力Fig.3 Vertical drift force of submerged spheroid

    圖4 橢球體二階縱傾波浪力矩Fig.4 Pitch drift moment of submerged spheroid

    圖1和圖2分別給出了不同面元數(shù)下得到的一階波浪力(矩)的計算結(jié)果,為了便于比較,利用ρgA(D/2)2和ρgA(D/2)3作無因次化處理,將結(jié)果與Newman利用三維面元法(面元數(shù)為1 024)得到的結(jié)果對比。從圖中可以看出,一階TEBEM方法得到的計算結(jié)果具有較好的收斂性;與Newman的結(jié)果比較,兩者的主要差別集中在長波階段,但誤差均控制在3%以內(nèi)。

    圖3和圖4分別給出了不同面元數(shù)下的二階垂向波浪力和二階縱傾波浪力矩,并利用ρgA2(D/2)和ρgA2(D/2)2作無因次化處理。由于Newman利用三維面元法計算的二階波浪力(矩)并沒有達到收斂,故而將其所有面元數(shù)的計算結(jié)果全部比較。容易得出,利用一階TEBEM方法求得的二階波浪力(矩)具有較好的收斂性;并且除了中長波段(0.4<λ/L<1.2)以外,利用一階TEBEM方法的結(jié)果與Newman 1 024網(wǎng)格數(shù)的結(jié)果基本一致。

    4 潛艇波浪力計算實例

    由上述第3節(jié)計算結(jié)果可知,基于一階TEBEM方法能夠得到令人滿意的計算結(jié)果。本節(jié)以國外某一公開潛艇模型[12]為例,與上一節(jié)進行類似內(nèi)容的計算。

    由于并未查閱到有關(guān)該模型重心位置、慣性系數(shù)等潛艇模型參數(shù),本節(jié)主要參考實際潛艇,估計相關(guān)參數(shù)。由于缺少該潛艇模型在無拘束狀態(tài)下的波浪力計算相關(guān)文獻,本節(jié)計算迎浪狀態(tài)下,不同網(wǎng)格數(shù)的一階波浪力(矩)和二階波浪力(矩)。表1給出了潛艇模型的主要參數(shù),圖5給出了潛艇模型表面網(wǎng)格劃分示意圖。

    表1 潛艇模型主要參數(shù)

    圖6和圖7給出了不同面元數(shù)下一階垂蕩力和一階縱搖力矩關(guān)于波浪頻率的計算結(jié)果,并分別利用ρgA(D/2)2和ρgA(D/2)3作無因次化處理。

    圖5 潛艇模型表面網(wǎng)格劃分Fig.5 The mesh on the submarine model surface

    圖6 Fr=0時,潛艇一階垂蕩波浪力Fig.6 First-order heave force of submarine as Fr=0

    圖7 Fr=0時,潛艇一階縱搖波浪力矩Fig.7 First-order pitch moment of submarine as Fr=0

    從圖6和圖7可以看出,一階垂蕩波浪力和一階縱搖波浪力矩利用1 000左右網(wǎng)格數(shù)就能得到收斂結(jié)果。

    圖8 Fr=0時,潛艇二階垂向波浪力Fig.8 Vertical drift force of submarine as Fr=0

    圖9 Fr=0時,潛艇二階縱傾波浪力矩Fig.9 Pitch drift moment of submarine as Fr=0

    圖8和圖9給出了在相同工況下,二階垂向波浪力和二階縱傾波浪力矩的計算結(jié)果。從圖中可以看出,二階垂向波浪力的計算結(jié)果很快得到收斂,而二階縱傾波浪力矩的計算結(jié)果收斂較慢。同時,二階垂向波浪力和二階縱傾波浪力矩在長波段的數(shù)值結(jié)果也能在一定程度上反映潛艇在垂直面內(nèi)的上浮甚至“拋甩”的現(xiàn)象。

    5 結(jié)論

    本文在頻域范疇內(nèi)利用泰勒展開邊界元方法計算潛艇一階和二階垂向波浪力和波浪力矩。以無限水深下近水面的橢球體和潛艇模型來驗證一階TEBEM方法的計算精度和收斂性。結(jié)論如下:

    1)通過與橢球體的計算結(jié)果對比,證明一階TEBEM方法具有較高的計算精度和收斂性。

    2)采用一階TEBEM方法對潛艇的波浪力(矩)計算:面元數(shù)1 000左右可以得到一階垂蕩波浪力和一階縱搖波浪力的收斂結(jié)果,面元數(shù)2 000左右可以得到二階垂向波浪力的收斂結(jié)果,二階縱傾波浪力矩則需要3 000左右的面元數(shù)。

    [1]戴遺山, 段文洋. 船舶在波浪中運動的勢流理論[M]. 北京: 國防工業(yè)出版社, 2008: 78-84, 107-113. DAI Yishan, DUAN Wenyang. Potential flow theory of ship motions in waves[M]. Beijing: National Defense Industry Press, 2008: 78-84, 107-113.

    [2]繆國平, 劉應(yīng)中, 糜振星. 二階波浪力的切片理論與潛體上的垂向定常力[J]. 水動力學(xué)研究與進展, 1993, 8(4): 435-447. MIAO Guoping, LIU Yingzhong, MI Zhenxing. A strip theory for the second-order wave forces and the steady vertical forces on submerged slender bodies[J]. Journal of hydrodynamics, 1993, 8(4): 435-447.

    [3]林青山, 繆國平, 李誼樂, 等. 潛艇上二階波浪力的考察[J]. 船舶力學(xué), 1998, 2(6): 15-23. LIN Qingshan, MIAO Guoping, LI Yile, et al. On the second-order wave forces on submarines[J]. Journal of ship mechanics, 1998, 2(6): 15-23.

    [4]馮學(xué)知, 繆泉明, 蔣強強. 近水面波浪力作用下潛體的非線性運動響應(yīng)[J]. 中國造船, 1996(3): 29-35. FENG Xuezhi, MIAO Quanming, JIANG Qiangqiang. Prediction of nonlinear motion responses of a submerged slender body running near free-surface due to wave-exciting forces[J]. Shipbuilding of China, 1996(3): 29-35.

    [5]馮學(xué)知, 蔣強強, 繆泉明, 等. 潛體波浪中近水面不同潛深和航向時運動和波浪力計算[J]. 船舶力學(xué), 2002, 6(2): 1-14. FENG Xuezhi, JIANG Qiangqiang, MIAO Quanming, et al. Computation of motion and wave forces for a submarine running near free surface in different depth of immersion and direction[J]. Journal of ship mechanics, 2002, 6(2): 1-14.

    [6]LEE C H, NEWMAN J N. First and second order wave effects on a submerged spheroid[J]. Journal of ship research, 1991, 35(3): 183-190.

    [7]徐剛, 段文洋. 常數(shù)分布Rankine源法與二階繞射問題精度研究[J]. 哈爾濱工程大學(xué)學(xué)報, 2010, 31(9): 1144-1152. XU Gang, DUAN Wenyang. Numerical investigation of second-order wave diffraction based on the Rankine source method[J]. Journal of Harbin Engineering University, 2010, 31(9): 1144-1152.

    [8]DUAN Wenyang. Taylor expansion boundary element method for floating body hydrodynamics[C]//Proceedings of the 27th International Workshop on Water Waves and Floating Bodies. Copenhagen, Denmark, 2012.

    [9]段文洋, 陳紀康, 趙彬彬. 基于泰勒展開邊界元法的深水浮體二階平均漂移力計算[J]. 哈爾濱工程大學(xué)學(xué)報, 2015, 36(3): 302-306. DUAN Wenyang, CHEN Jikang, ZHAO Binbin. Calculation of second-order mean drift loads for the deepwater floating body based on the Taylor expansion boundary element method[J]. Journal of Harbin Engineering University, 2015, 36(3): 302-306.

    [10]DUAN Wenyang, CHEN Jikang, ZHAO Binbin. Second-order Taylor expansion boundary element method for the second-order wave diffraction problem[J]. Engineering Analysis with Boundary Elements, 2015, 58: 140-150.

    [11]DUAN Wenyang, CHEN Jikang, ZHAO Binbin. Second-order Taylor expansion boundary element method for the second-order wave radiation problem[J]. Applied Ocean Research, 2015, 52: 12-26.

    [12]GROVES N C, HUANG T T, CHANG M S. Geometric characteristics of DARPA SUBOFF models (DTRC model Nos. 5470 and 5471). DTRC/SHD-1298-01[R]. Bethesda, Maryland: David Taylor Research Centre, 1989.

    Calculation of vertical second-order drift loads on a submarine floating near the free water surface based on Taylor expansion boundary element method

    DUAN Wenyang,WANG Lijia,CHEN Jikang,ZHAO Binbin

    (College of Shipbuilding Engineering, Harbin Engineering University, Harbin 150001, China)

    A submarine traveling near the free water surface experiences difficulty in depth control. Thus, the vertical second-order mean force and moment acting on the submarine were studied. One-order Taylor expansion boundary element method (TEBEM) is used to improve the solution accuracy of the tangential induced velocity.This method keeps the first-order derivative for the dipole strength of the border integral equation and then takes the two tangential derivatives with respect to the field points on the boundary to form the closed equations in the third formula of the Green function. The velocity potential and the tangential induced velocity along object surface can be solved directly. This method was applied to calculate the vertical second-order drift loads of a submerged spheroid in head seas. The results show good agreement with previous numerical solutions. The vertical second-order drift loads of a submarine model were also calculated. TEBEM was found to achieve higher accuracy and quick convergence.

    Taylor expansion boundary element method; second-order vertical drift force; second-order pitch drift moment; submarine; frequency-domain; tangential induced velocity

    2016-06-11.

    時間:2016-12-12.

    國家自然科學(xué)基金資助項目(11272097).

    段文洋(1967-),男,教授,博士生導(dǎo)師; 陳紀康(1986-),男,博士后.

    陳紀康,E-mail: cjkhrb@sina.com.

    10.11990/jheu.201606032

    O352

    A

    1006-7043(2017)01-0008-05

    段文洋, 王隸加, 陳紀康,等. 基于泰勒展開邊界元法的近水面潛艇垂向二階波浪力(矩)計算[J]. 哈爾濱工程大學(xué)學(xué)報, 2017, 38(1): 8-12. DUAN Wenyang,WANG Lijia,CHEN Jikang,et al. Calculation of vertical second-order drift loads on a submarine floating near the free water surface based on Taylor expansion boundary element method[J]. Journal of Harbin Engineering University, 2017, 38(1): 8-12.

    網(wǎng)絡(luò)出版地址:http://www.cnki.net/kcms/detail/23.1390.u.20161212.0920.008.html

    猜你喜歡
    橢球元法泰勒
    獨立坐標系橢球變換與坐標換算
    橢球槽宏程序編制及其Vericut仿真
    智能制造(2021年4期)2021-11-04 08:54:44
    換元法在解題中的運用
    基于離散元法的礦石對溜槽沖擊力的模擬研究
    重型機械(2019年3期)2019-08-27 00:58:46
    橢球精加工軌跡及程序設(shè)計
    基于外定界橢球集員估計的純方位目標跟蹤
    一起綿羊泰勒焦蟲病的診斷治療經(jīng)過
    換元法在解題中的應(yīng)用
    “微元法”在含電容器電路中的應(yīng)用
    泰勒公式的簡單應(yīng)用
    河南科技(2014年14期)2014-02-27 14:12:08
    交换朋友夫妻互换小说| 国产老妇女一区| 午夜福利在线在线| 欧美国产精品一级二级三级 | 纵有疾风起免费观看全集完整版| 日韩三级伦理在线观看| 69av精品久久久久久| 性色av一级| 国产探花在线观看一区二区| 交换朋友夫妻互换小说| 男女国产视频网站| 精品一区二区三区视频在线| 九九爱精品视频在线观看| 国产探花极品一区二区| 精品国产乱码久久久久久小说| 成人漫画全彩无遮挡| 最近的中文字幕免费完整| 青春草视频在线免费观看| 夜夜看夜夜爽夜夜摸| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美日韩卡通动漫| 成人鲁丝片一二三区免费| 国产午夜精品一二区理论片| 亚洲av国产av综合av卡| 深爱激情五月婷婷| 18禁动态无遮挡网站| 内射极品少妇av片p| 久久99精品国语久久久| 久久久久久久精品精品| 18禁动态无遮挡网站| 久久精品国产亚洲av涩爱| 欧美亚洲 丝袜 人妻 在线| 久久久精品94久久精品| 成年免费大片在线观看| 亚洲国产日韩一区二区| 久久久久精品久久久久真实原创| 欧美激情在线99| 如何舔出高潮| 网址你懂的国产日韩在线| 最新中文字幕久久久久| 午夜老司机福利剧场| 亚洲久久久久久中文字幕| 亚洲丝袜综合中文字幕| 中文天堂在线官网| 欧美性猛交╳xxx乱大交人| 国产精品国产三级国产专区5o| 99久久精品一区二区三区| 1000部很黄的大片| 丰满乱子伦码专区| 嘟嘟电影网在线观看| 人人妻人人看人人澡| 中文字幕av成人在线电影| 国产黄a三级三级三级人| 色婷婷久久久亚洲欧美| 国产精品国产av在线观看| 亚洲精品视频女| 精品一区在线观看国产| 干丝袜人妻中文字幕| 日韩成人伦理影院| 欧美bdsm另类| 婷婷色综合www| 国产色爽女视频免费观看| 在线天堂最新版资源| 七月丁香在线播放| 亚洲最大成人手机在线| 99热这里只有是精品50| 午夜福利在线观看免费完整高清在| 成人国产av品久久久| 国产一区二区三区综合在线观看 | 久久精品久久精品一区二区三区| 中文字幕制服av| 99久久九九国产精品国产免费| 看非洲黑人一级黄片| 色播亚洲综合网| 欧美高清性xxxxhd video| 白带黄色成豆腐渣| 2018国产大陆天天弄谢| 国产中年淑女户外野战色| kizo精华| 在线观看一区二区三区| 婷婷色综合www| 精品久久久噜噜| 午夜福利网站1000一区二区三区| 亚洲最大成人av| 青春草国产在线视频| 国产伦理片在线播放av一区| 九九久久精品国产亚洲av麻豆| 久久午夜福利片| av一本久久久久| 在线免费观看不下载黄p国产| 免费av观看视频| 亚洲真实伦在线观看| 久久综合国产亚洲精品| 国产精品蜜桃在线观看| 欧美区成人在线视频| 亚洲精品色激情综合| 亚洲av免费高清在线观看| 成人综合一区亚洲| 亚洲欧美精品专区久久| 国产黄片视频在线免费观看| 最近最新中文字幕大全电影3| 永久网站在线| 激情五月婷婷亚洲| 免费黄网站久久成人精品| av天堂中文字幕网| 精品久久久久久久久av| 水蜜桃什么品种好| 99久久九九国产精品国产免费| 精品久久久久久电影网| 最近2019中文字幕mv第一页| 亚洲国产日韩一区二区| 中文精品一卡2卡3卡4更新| 日韩 亚洲 欧美在线| 欧美高清成人免费视频www| av在线观看视频网站免费| 在线观看美女被高潮喷水网站| 女的被弄到高潮叫床怎么办| 麻豆久久精品国产亚洲av| 精品久久国产蜜桃| 亚洲久久久久久中文字幕| 国产有黄有色有爽视频| 欧美精品国产亚洲| 美女视频免费永久观看网站| 身体一侧抽搐| 亚洲精品国产色婷婷电影| 欧美丝袜亚洲另类| 成人亚洲精品一区在线观看 | 欧美性感艳星| tube8黄色片| 22中文网久久字幕| 美女xxoo啪啪120秒动态图| 国产精品蜜桃在线观看| 成人毛片a级毛片在线播放| 亚洲av电影在线观看一区二区三区 | 亚洲av成人精品一区久久| 永久网站在线| 亚洲色图综合在线观看| 成人综合一区亚洲| 精品国产乱码久久久久久小说| 插阴视频在线观看视频| 中文资源天堂在线| 免费电影在线观看免费观看| 日本黄色片子视频| 国产av码专区亚洲av| 国产黄a三级三级三级人| 欧美丝袜亚洲另类| 午夜福利在线观看免费完整高清在| 国产一区二区亚洲精品在线观看| 狠狠精品人妻久久久久久综合| 插逼视频在线观看| 免费大片18禁| 97在线人人人人妻| 纵有疾风起免费观看全集完整版| 午夜福利视频精品| 国产精品国产av在线观看| av国产免费在线观看| 性色avwww在线观看| 日本一本二区三区精品| 日本午夜av视频| 在线观看免费高清a一片| 天堂网av新在线| freevideosex欧美| av国产久精品久网站免费入址| 成人综合一区亚洲| 精品国产三级普通话版| 久久久久国产精品人妻一区二区| 人人妻人人澡人人爽人人夜夜| 国产av不卡久久| 久久精品久久久久久久性| 亚洲成色77777| 中文精品一卡2卡3卡4更新| 国产亚洲一区二区精品| 亚洲第一区二区三区不卡| 国产成人精品婷婷| 亚洲美女搞黄在线观看| 男人舔奶头视频| 69av精品久久久久久| 中文字幕久久专区| 亚洲自拍偷在线| 日韩欧美精品免费久久| 中国美白少妇内射xxxbb| 亚洲国产av新网站| 青春草国产在线视频| 在线a可以看的网站| 成人特级av手机在线观看| 亚洲精品国产av蜜桃| 制服丝袜香蕉在线| 色播亚洲综合网| 亚洲丝袜综合中文字幕| 美女国产视频在线观看| 欧美丝袜亚洲另类| 麻豆久久精品国产亚洲av| 午夜亚洲福利在线播放| 99久久精品一区二区三区| 少妇丰满av| av女优亚洲男人天堂| 伊人久久精品亚洲午夜| 有码 亚洲区| 亚洲精品国产av蜜桃| 日韩欧美 国产精品| 亚洲,一卡二卡三卡| 国产色婷婷99| 菩萨蛮人人尽说江南好唐韦庄| 97精品久久久久久久久久精品| 岛国毛片在线播放| 国产精品爽爽va在线观看网站| 亚洲精品国产av蜜桃| 男的添女的下面高潮视频| 亚洲欧美日韩卡通动漫| 免费看不卡的av| 亚洲美女视频黄频| 国产视频内射| 麻豆久久精品国产亚洲av| 久久精品久久精品一区二区三区| 久久久欧美国产精品| 亚洲欧美成人综合另类久久久| 下体分泌物呈黄色| 欧美区成人在线视频| av在线天堂中文字幕| 最近最新中文字幕大全电影3| 国产精品偷伦视频观看了| 精品久久久精品久久久| 日本熟妇午夜| 九色成人免费人妻av| 欧美日韩亚洲高清精品| 亚洲精品,欧美精品| 日韩欧美精品免费久久| 五月开心婷婷网| 男女边吃奶边做爰视频| 高清毛片免费看| 久久久亚洲精品成人影院| 天天躁日日操中文字幕| 内射极品少妇av片p| 久久人人爽人人爽人人片va| 蜜桃久久精品国产亚洲av| 国产爱豆传媒在线观看| 禁无遮挡网站| 久久鲁丝午夜福利片| 九色成人免费人妻av| 亚洲精品乱码久久久v下载方式| 日韩一本色道免费dvd| 成人无遮挡网站| 免费在线观看成人毛片| 国产午夜福利久久久久久| av国产精品久久久久影院| 精品一区二区三卡| 我的女老师完整版在线观看| 丝瓜视频免费看黄片| 亚洲精品久久久久久婷婷小说| 亚洲综合精品二区| www.av在线官网国产| 有码 亚洲区| 97在线人人人人妻| 精品国产三级普通话版| 国产精品蜜桃在线观看| 亚洲成人精品中文字幕电影| 少妇熟女欧美另类| 久久久久久九九精品二区国产| 嫩草影院入口| 亚洲欧美日韩东京热| 精品久久国产蜜桃| 欧美97在线视频| 九九久久精品国产亚洲av麻豆| 日本黄色片子视频| 亚洲激情五月婷婷啪啪| 国产人妻一区二区三区在| 亚洲经典国产精华液单| 国产高清有码在线观看视频| 尾随美女入室| 97超视频在线观看视频| 亚洲精品乱码久久久v下载方式| av在线观看视频网站免费| 日韩亚洲欧美综合| 亚洲经典国产精华液单| 亚洲国产精品成人综合色| 亚洲精品一二三| 中文精品一卡2卡3卡4更新| av线在线观看网站| 亚洲美女视频黄频| 欧美高清性xxxxhd video| av播播在线观看一区| 在线免费十八禁| 人体艺术视频欧美日本| 久久午夜福利片| 十八禁网站网址无遮挡 | 亚洲精品色激情综合| 亚洲人成网站在线播| 国产精品国产三级专区第一集| 国产老妇伦熟女老妇高清| 国产精品一区二区在线观看99| 亚洲av一区综合| 成人午夜精彩视频在线观看| 有码 亚洲区| av在线老鸭窝| 久久99热这里只有精品18| 嘟嘟电影网在线观看| 婷婷色麻豆天堂久久| 日本欧美国产在线视频| 肉色欧美久久久久久久蜜桃 | 亚洲av日韩在线播放| 午夜福利视频精品| 国产毛片a区久久久久| 嫩草影院精品99| 夜夜爽夜夜爽视频| xxx大片免费视频| 蜜桃亚洲精品一区二区三区| 国产精品无大码| 亚洲欧美日韩另类电影网站 | 亚洲熟女精品中文字幕| 韩国高清视频一区二区三区| 一区二区三区精品91| 久久久久精品性色| 日韩欧美 国产精品| 久久国内精品自在自线图片| 国产午夜福利久久久久久| 欧美日韩国产mv在线观看视频 | 国产一区二区三区av在线| 最近手机中文字幕大全| 尾随美女入室| 99久久中文字幕三级久久日本| 视频中文字幕在线观看| 亚洲怡红院男人天堂| 禁无遮挡网站| 尾随美女入室| 亚洲国产av新网站| 岛国毛片在线播放| 亚洲av免费高清在线观看| 国产欧美日韩精品一区二区| 亚洲av成人精品一二三区| 黄色视频在线播放观看不卡| 少妇的逼水好多| 各种免费的搞黄视频| 国产亚洲精品久久久com| 国产精品99久久久久久久久| 各种免费的搞黄视频| 国产欧美日韩一区二区三区在线 | 亚洲在久久综合| 国产成人精品婷婷| 又爽又黄a免费视频| 日本午夜av视频| 亚洲国产精品成人久久小说| 99re6热这里在线精品视频| 777米奇影视久久| 午夜精品国产一区二区电影 | 卡戴珊不雅视频在线播放| 中文资源天堂在线| 一区二区av电影网| 在线 av 中文字幕| 国产白丝娇喘喷水9色精品| 国产一区二区亚洲精品在线观看| 亚洲图色成人| 欧美日韩一区二区视频在线观看视频在线 | 久久影院123| 成年av动漫网址| 午夜福利网站1000一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产精品999| 美女国产视频在线观看| 水蜜桃什么品种好| 国产精品一区二区性色av| 精品久久久久久久久av| 老师上课跳d突然被开到最大视频| 免费观看a级毛片全部| 乱系列少妇在线播放| 国产高清不卡午夜福利| 高清av免费在线| 欧美高清性xxxxhd video| 亚洲精品久久久久久婷婷小说| 国产片特级美女逼逼视频| 免费av不卡在线播放| 校园人妻丝袜中文字幕| 91精品一卡2卡3卡4卡| 又粗又硬又长又爽又黄的视频| 国产日韩欧美亚洲二区| 精品国产三级普通话版| 欧美3d第一页| 国产极品天堂在线| 久久久久精品性色| 插逼视频在线观看| 亚洲av一区综合| 中文字幕久久专区| 菩萨蛮人人尽说江南好唐韦庄| 在线 av 中文字幕| 久久久久国产网址| 国产亚洲91精品色在线| 国产精品蜜桃在线观看| 欧美精品人与动牲交sv欧美| 亚洲精品久久午夜乱码| 久久久久久久亚洲中文字幕| 欧美日韩在线观看h| 国产亚洲精品久久久com| 成人高潮视频无遮挡免费网站| 自拍欧美九色日韩亚洲蝌蚪91 | 高清av免费在线| 如何舔出高潮| 国内精品美女久久久久久| 国产av国产精品国产| 日韩视频在线欧美| 免费观看a级毛片全部| 乱系列少妇在线播放| 高清午夜精品一区二区三区| 毛片女人毛片| 亚洲精品自拍成人| 国产精品久久久久久精品古装| 亚洲欧美日韩无卡精品| 丝瓜视频免费看黄片| 亚洲综合色惰| 联通29元200g的流量卡| 狂野欧美激情性xxxx在线观看| 国产成人午夜福利电影在线观看| 久久人人爽av亚洲精品天堂 | 99热国产这里只有精品6| 九九爱精品视频在线观看| 99热网站在线观看| 男女无遮挡免费网站观看| 国产一区二区亚洲精品在线观看| 熟女av电影| 日本爱情动作片www.在线观看| 伦理电影大哥的女人| av国产久精品久网站免费入址| 免费观看av网站的网址| 日韩 亚洲 欧美在线| 只有这里有精品99| 特大巨黑吊av在线直播| 我的女老师完整版在线观看| 一级毛片aaaaaa免费看小| 国国产精品蜜臀av免费| 欧美人与善性xxx| 69人妻影院| 三级经典国产精品| 久久久色成人| 免费看a级黄色片| 免费看日本二区| 国产精品秋霞免费鲁丝片| 国产 一区 欧美 日韩| 高清av免费在线| 亚洲精品一二三| 日本爱情动作片www.在线观看| 国产白丝娇喘喷水9色精品| 插阴视频在线观看视频| 中文字幕av成人在线电影| 建设人人有责人人尽责人人享有的 | 国产免费一级a男人的天堂| 熟女人妻精品中文字幕| 亚洲精品日韩av片在线观看| 极品少妇高潮喷水抽搐| 精品久久国产蜜桃| 日韩国内少妇激情av| 久久午夜福利片| 三级国产精品欧美在线观看| 麻豆精品久久久久久蜜桃| 3wmmmm亚洲av在线观看| av在线老鸭窝| 激情 狠狠 欧美| 伊人久久精品亚洲午夜| 国产在视频线精品| 婷婷色麻豆天堂久久| 特大巨黑吊av在线直播| 爱豆传媒免费全集在线观看| 波多野结衣巨乳人妻| 久久久久久久精品精品| 国产熟女欧美一区二区| 51国产日韩欧美| av播播在线观看一区| 亚洲美女搞黄在线观看| 久久人人爽人人爽人人片va| 伦精品一区二区三区| 性插视频无遮挡在线免费观看| 亚洲欧美一区二区三区黑人 | 简卡轻食公司| 日本色播在线视频| 国产日韩欧美亚洲二区| 男女国产视频网站| 成年版毛片免费区| 亚洲天堂av无毛| 亚洲国产成人一精品久久久| 国产成年人精品一区二区| 国产成人午夜福利电影在线观看| 卡戴珊不雅视频在线播放| 永久网站在线| 欧美3d第一页| 最近最新中文字幕免费大全7| a级一级毛片免费在线观看| 久久99热这里只频精品6学生| 三级男女做爰猛烈吃奶摸视频| 国产精品熟女久久久久浪| 日产精品乱码卡一卡2卡三| 蜜桃久久精品国产亚洲av| 国产亚洲91精品色在线| 99re6热这里在线精品视频| 亚洲国产最新在线播放| 国产乱人视频| 久久久久久久大尺度免费视频| 我的女老师完整版在线观看| 国产精品人妻久久久久久| 少妇猛男粗大的猛烈进出视频 | 亚洲人与动物交配视频| 日本午夜av视频| 99热这里只有是精品50| 国产精品女同一区二区软件| 三级男女做爰猛烈吃奶摸视频| 午夜爱爱视频在线播放| eeuss影院久久| 免费少妇av软件| 超碰av人人做人人爽久久| 人体艺术视频欧美日本| 久久精品国产亚洲网站| 久久久久久伊人网av| 赤兔流量卡办理| 色播亚洲综合网| 日韩欧美 国产精品| 如何舔出高潮| 大片电影免费在线观看免费| 伊人久久精品亚洲午夜| 18禁在线无遮挡免费观看视频| 两个人的视频大全免费| 国产精品伦人一区二区| 国产精品国产三级国产av玫瑰| 久久久欧美国产精品| 六月丁香七月| 天天躁夜夜躁狠狠久久av| 国产爽快片一区二区三区| 国产黄a三级三级三级人| 免费大片黄手机在线观看| 在现免费观看毛片| 亚洲精品日韩av片在线观看| 黄色欧美视频在线观看| 久久久久久久精品精品| 在线看a的网站| 国产v大片淫在线免费观看| 五月玫瑰六月丁香| 在线精品无人区一区二区三 | 在线观看一区二区三区激情| 欧美性感艳星| 欧美xxⅹ黑人| 晚上一个人看的免费电影| 大码成人一级视频| 看免费成人av毛片| 在线 av 中文字幕| 少妇被粗大猛烈的视频| 男女无遮挡免费网站观看| 国精品久久久久久国模美| 亚洲av日韩在线播放| 欧美 日韩 精品 国产| 五月玫瑰六月丁香| 人妻制服诱惑在线中文字幕| 国产人妻一区二区三区在| 又大又黄又爽视频免费| 99久久精品一区二区三区| 毛片一级片免费看久久久久| 午夜日本视频在线| 国产亚洲av片在线观看秒播厂| 99九九线精品视频在线观看视频| 肉色欧美久久久久久久蜜桃 | 免费不卡的大黄色大毛片视频在线观看| 亚洲精品一二三| 日韩免费高清中文字幕av| 中文字幕av成人在线电影| 直男gayav资源| 久久影院123| 免费观看在线日韩| 91精品国产九色| 少妇的逼好多水| 成人亚洲精品一区在线观看 | 国产在线一区二区三区精| 我要看日韩黄色一级片| 欧美日韩一区二区视频在线观看视频在线 | 中文资源天堂在线| 国产亚洲av嫩草精品影院| 国产 一区精品| 免费人成在线观看视频色| 久久韩国三级中文字幕| 亚洲av国产av综合av卡| 真实男女啪啪啪动态图| av在线播放精品| 中文天堂在线官网| 精品亚洲乱码少妇综合久久| 青青草视频在线视频观看| 欧美xxⅹ黑人| 国产爽快片一区二区三区| 国产高清国产精品国产三级 | 日韩三级伦理在线观看| 国产欧美亚洲国产| 97精品久久久久久久久久精品| 国产精品一及| 欧美人与善性xxx| 下体分泌物呈黄色| 国产成人91sexporn| 国产伦理片在线播放av一区| 在线 av 中文字幕| 国产精品秋霞免费鲁丝片| 国产毛片在线视频| 亚洲国产色片| www.色视频.com| 内射极品少妇av片p| 精品人妻熟女av久视频| 七月丁香在线播放| 国产成人freesex在线| 最近2019中文字幕mv第一页| 成人国产av品久久久| 欧美zozozo另类| 久久精品人妻少妇| 男女那种视频在线观看| 一二三四中文在线观看免费高清| 亚洲精品亚洲一区二区| 又黄又爽又刺激的免费视频.| 国精品久久久久久国模美| 午夜日本视频在线| 在线观看一区二区三区激情| 91精品伊人久久大香线蕉| 国产黄频视频在线观看| 成人一区二区视频在线观看| 午夜老司机福利剧场| 日韩中字成人| 97超碰精品成人国产|