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

    基于CFD的實(shí)船阻力快速預(yù)估方法

    2014-10-30 08:13:08王志南沈興榮范佘明
    船舶與海洋工程 2014年3期
    關(guān)鍵詞:實(shí)船船模船型

    王志南,沈興榮,范佘明

    (中國(guó)船舶及海洋工程設(shè)計(jì)研究院,上海 200011)

    0 引 言

    船舶設(shè)計(jì)過程中,無論是初步設(shè)計(jì)階段的多方案選型,還是最終方案確定后的航速性能評(píng)估,阻力性能都是評(píng)判的主要依據(jù);不同方案快速性能的比較、主機(jī)的選型、航速的預(yù)估、能效指標(biāo)(諸如 EEDI、EEOI)的評(píng)判等都是建立在船舶阻力數(shù)據(jù)的基礎(chǔ)之上。艦船阻力性能的預(yù)報(bào)有多種方法,精度較高的有模型試驗(yàn)方法,相對(duì)粗略的有經(jīng)驗(yàn)公式回歸方法,還有當(dāng)下正迅速發(fā)展的基于 CFD的數(shù)值模擬方法,這些方法各有特點(diǎn),有的較為精確但費(fèi)時(shí)耗力,有的簡(jiǎn)單快捷卻精度稍差,在船舶設(shè)計(jì)的不同階段需要根據(jù)其各自特點(diǎn)選擇使用不同的預(yù)報(bào)方法。初步設(shè)計(jì)的多方案選型階段,船舶線型尚未最終確定,對(duì)每一方案都進(jìn)行水池試驗(yàn)顯然不可行,這就需要一種既快速方便、又相對(duì)準(zhǔn)確可靠的預(yù)報(bào)方法。

    1 船舶阻力預(yù)估方法

    船舶阻力的預(yù)估方法有多種,最為精確的是水池模型試驗(yàn)方法,除此之外有大量的根據(jù)船模系列試驗(yàn)結(jié)果或者是在總結(jié)、分析大量船模試驗(yàn)和實(shí)船試航數(shù)據(jù)基礎(chǔ)上得出的近似估算方法。模型試驗(yàn)是在滿足模型、實(shí)船弗勞德數(shù)相等、幾何相似并結(jié)合實(shí)船試航結(jié)果相關(guān)分析基礎(chǔ)上進(jìn)行的一種精確的阻力預(yù)報(bào)手段;其他近似的阻力估算方法諸如系列圖譜法、經(jīng)驗(yàn)公式法等,其預(yù)估的精度一方面取決于設(shè)計(jì)船與母型船或設(shè)計(jì)船與各圖譜所依據(jù)的船模系列之間的相似程度,另一方面取決于設(shè)計(jì)師的經(jīng)驗(yàn)。阻力近似估算法按計(jì)算內(nèi)容可分為兩類[1]:一類是直接近似估算總阻力或有效功率;另一類是估算剩余阻力,用相當(dāng)平板公式計(jì)算摩擦阻力,進(jìn)而得到總阻力。按照阻力近似估算方法的表達(dá)形式可分為圖譜法和回歸公式法兩種;而根據(jù)估算方法資料來源進(jìn)行分類,可分為母型船數(shù)據(jù)估算法、船模系列資料估算法和依據(jù)歸納船模和實(shí)船資料估算法等3類估算方法。

    長(zhǎng)期以來,船舶工程師根據(jù)大量系列船模的水池試驗(yàn)資料,建立了基于系列船模試驗(yàn)資料的船舶阻力預(yù)估方法,常用的比較有代表性的系列圖譜有泰勒(Taylor)系列、陶德(Todd)60系列、BSRA系列、SSPA系列、Delft系列、擴(kuò)展的泰勒系列等,國(guó)內(nèi)的也有CSSRC系列、淺吃水肥大型船系列、長(zhǎng)江船系列等[2]。系列圖譜法是上世紀(jì)80年代最為流行的阻力預(yù)估方法,其使用方便、直觀,可以直接生成船體線型方案。但這些早期的系列圖譜代表的線型技術(shù)已經(jīng)相對(duì)落后,其預(yù)報(bào)的船舶性能指標(biāo)不能真實(shí)反映現(xiàn)在的設(shè)計(jì)水平。

    隨著水池船模試驗(yàn)數(shù)據(jù)和實(shí)船試航數(shù)據(jù)的積累,在此基礎(chǔ)上進(jìn)行統(tǒng)計(jì)分析,回歸出阻力估算公式,這一類的阻力估算方法包括早期的艾亞(Ayre)法、蘭潑–凱勒(Lap-Keller)法,現(xiàn)在較為常用的有荷蘭MARIN水池的Holtop84方法、Hollenbach98方法、Davidson回歸法等[3],世界各國(guó)著名水池大都依據(jù)自身水池試驗(yàn)數(shù)據(jù)和實(shí)船試驗(yàn)結(jié)果建立了自己的阻力相關(guān)分析、預(yù)報(bào)方法?;谀P?實(shí)船試驗(yàn)數(shù)據(jù)統(tǒng)計(jì)回歸建立的阻力估算方法,其估算的精度一方面受到數(shù)據(jù)樣本積累數(shù)量的影響,另一方面僅僅代表了所收集樣本的平均水平,有可能偏于保守。

    船舶設(shè)計(jì)過程中,很多情況下是在母型船基礎(chǔ)上加以改動(dòng)從而得到新設(shè)計(jì)方案,新設(shè)計(jì)船型的阻力可以通過簡(jiǎn)單易行的母型船數(shù)據(jù)估算法進(jìn)行預(yù)估,如果新設(shè)計(jì)方案與母型船比較相似,則這類方法也可給出較為滿意的結(jié)果。常用的母型船數(shù)據(jù)估算法有海軍系數(shù)法、引伸比較定律法、基爾斯修正母型船剩余阻力法等等。目前很多船舶設(shè)計(jì)研究機(jī)構(gòu)都建立了船型數(shù)據(jù)庫,通過主要船型參數(shù)從數(shù)據(jù)庫中查詢與新設(shè)計(jì)方案相匹配的母型船,調(diào)取最匹配的母型船模型及實(shí)船試驗(yàn)資料,在此基礎(chǔ)上估算新設(shè)計(jì)方案的阻力性能,其估算的精度與數(shù)據(jù)庫中收集船型數(shù)量的多少及其覆蓋面有關(guān),這一方法目前被造船界廣泛認(rèn)同和接受。

    與模型試驗(yàn)中要求船體幾何完全相似不同,上述估算方法都是建立在主尺度或船型參數(shù)相近的基礎(chǔ)上的預(yù)估,因此這些方法使用起來較為靈活,但精度相對(duì)較差,僅可作粗略的評(píng)估比較之用。隨著計(jì)算機(jī)技術(shù)的發(fā)展和CFD應(yīng)用開發(fā)的不斷推進(jìn),基于數(shù)值模擬的船舶阻力預(yù)報(bào)方法研究取得了令人矚目的進(jìn)展。其中最為有名的是歐洲的虛擬水池研究項(xiàng)目,旨在通過研究開發(fā)一套綜合集成的船舶水動(dòng)力性能數(shù)值水池平臺(tái),提高歐洲先進(jìn)船舶設(shè)計(jì)競(jìng)爭(zhēng)力,增強(qiáng)歐洲水動(dòng)力學(xué)服務(wù)的范圍和質(zhì)量及其研發(fā)能力[4]。國(guó)內(nèi)的研究主要集中于模型尺度下船舶水動(dòng)力性能的數(shù)值模擬[5]。這一類方法考慮船舶真實(shí)幾何形狀,考慮流體介質(zhì)的特性,追求“真實(shí)水池”中進(jìn)行的物理模型試驗(yàn)功能,通過求解控制方程得到物體的受力情況及流場(chǎng)特性,其終極目標(biāo)是模擬真實(shí)外界環(huán)境下實(shí)尺度船舶水動(dòng)力性能的數(shù)值預(yù)報(bào),要達(dá)到這一目標(biāo)還有諸多困難需要解決,包括高雷諾數(shù)下流動(dòng)機(jī)理、湍流模型和網(wǎng)格技術(shù)等[6]。相比于水池試驗(yàn)的費(fèi)時(shí)耗力、長(zhǎng)周期而言,基于數(shù)值模擬的船舶水動(dòng)力性能預(yù)報(bào)顯然具有很大的優(yōu)勢(shì),但受限于目前的計(jì)算機(jī)硬件條件和CFD技術(shù)理論基礎(chǔ)的發(fā)展(包括湍流模型、邊界層處理和網(wǎng)格技術(shù)等),基于數(shù)值模擬的船舶水動(dòng)力性能預(yù)報(bào)技術(shù)尚不可完全取代水池試驗(yàn)。數(shù)值模擬方法考慮了船體真實(shí)的幾何形狀以及其周圍介質(zhì)的特性,因此較圖譜或回歸公式等近似的估算方法更為準(zhǔn)確。

    2 基于CFD的船舶阻力快速預(yù)報(bào)方法[7, 8]

    基于CFD的船舶阻力快速預(yù)報(bào)方法的思路是:通過CFD軟件快速求出船模各阻力分量的系數(shù)以及形狀因子,其中興波阻力系數(shù)Cw采用勢(shì)流方法計(jì)算得到;黏性阻力系數(shù)Cv在黏性流場(chǎng)中采用迭模方法計(jì)算得到,相比于水池試驗(yàn)法,認(rèn)為低速時(shí)興波可以忽略,從而測(cè)量得到黏性阻力;數(shù)值迭模方法避免了興波的影響,計(jì)算得到的黏性阻力數(shù)據(jù)更加穩(wěn)定,從而可以計(jì)算得到較為穩(wěn)定的形狀因子;摩擦阻力系數(shù)Cf采用1957ITTC公式計(jì)算,再結(jié)合以往船模–實(shí)船相關(guān)分析獲得的相應(yīng)船型的實(shí)船相關(guān)系數(shù)Ca,由下式得到實(shí)船的總阻力系數(shù)Cts,進(jìn)而可以估算實(shí)船有效功率 PE:

    式中:形狀因子(1+k)=Cv/Cf;Rs——實(shí)船阻力;Vs——實(shí)船航速。

    計(jì)算采用商業(yè)軟件SHIPFLOW,將船體周圍的流場(chǎng)劃分為如圖1所示的3個(gè)區(qū)域,以取得有效的近似流動(dòng)方程,提高計(jì)算效率。區(qū)域1采用基于Dawson方法的面元法計(jì)算勢(shì)流,通過迭代計(jì)算得到滿足邊界條件的解,算得興波阻力;在船體中前部貼近船體的薄邊界層區(qū)域2采用積分法沿邊界層厚度求解動(dòng)量的微分方程,采用勢(shì)流計(jì)算的壓力分布作為輸入?yún)?shù);在區(qū)域3采用RANS方程求解,將控制方程轉(zhuǎn)換到由泊松方程生成的貼體坐標(biāo)系中求解,RANS方程解的邊界條件由區(qū)域1、2的勢(shì)流和邊界層計(jì)算獲得,結(jié)合EASM、KWSST等湍流模型來封閉控制方程求解流場(chǎng)信息,通過區(qū)域3的計(jì)算可以求得黏性阻力,從而求出形狀因子。

    圖1 船體周圍求解區(qū)域

    基于CFD的實(shí)船阻力預(yù)估方法中的核心問題是興波阻力的計(jì)算和形狀因子的確定,對(duì)數(shù)值計(jì)算過程中影響興波阻力計(jì)算和形狀因子確定的因素需要進(jìn)行系統(tǒng)考察,研究中對(duì)網(wǎng)格的劃分、計(jì)算域的范圍、求解參數(shù)的設(shè)置等因素進(jìn)行了考察。圖2是興波阻力系數(shù)隨著網(wǎng)格密度變化和計(jì)算域上游邊界距離的變化趨勢(shì),可見網(wǎng)格質(zhì)量、計(jì)算域范圍對(duì)數(shù)值計(jì)算結(jié)果有較大的影響。要得到穩(wěn)定可靠的結(jié)果,需要對(duì)這些因素展開細(xì)致的研究。經(jīng)過系統(tǒng)研究,掌握了網(wǎng)格劃分技術(shù)、合適計(jì)算域范圍的選取、相關(guān)求解策略等關(guān)鍵技術(shù)。

    圖2 興波阻力的影響因素

    3 計(jì)算結(jié)果分析

    分別選取細(xì)長(zhǎng)型的集裝箱船和肥大型的散貨船預(yù)報(bào)阻力,船型的主要參數(shù)見表1。

    表1 對(duì)象船型的主要參數(shù)

    采用上述數(shù)值方法分別預(yù)估兩型船對(duì)應(yīng)試航工況附近速度點(diǎn)的阻力,對(duì)于集裝箱船和散貨船的相關(guān)系數(shù)Ca分別取–0.0001和–0.00012;根據(jù)模型水池試驗(yàn)結(jié)果得到的自航因子,由推進(jìn)效率和有效功率得到實(shí)船的收到功率。同時(shí)搜集了該兩型船實(shí)船試航數(shù)據(jù),以此來驗(yàn)證基于CFD的預(yù)報(bào)結(jié)果[9]。相關(guān)數(shù)據(jù)的比較見圖3、4。圖中實(shí)線是基于CFD的預(yù)報(bào)值,三角形點(diǎn)是實(shí)船試航結(jié)果,通過比較可見集裝箱船的阻力預(yù)報(bào)偏低,其中在75%主機(jī)負(fù)荷時(shí),預(yù)報(bào)值較實(shí)船試航結(jié)果偏低約3.8%,在90%主機(jī)負(fù)荷時(shí),預(yù)報(bào)值與試航結(jié)果吻合較好,預(yù)報(bào)值與試航結(jié)果的偏差在1%以內(nèi),滿負(fù)荷工況下,預(yù)報(bào)值較試航結(jié)果低約4.2%;散貨船預(yù)報(bào)結(jié)果偏高,75%和85%主機(jī)負(fù)荷的工況下,基于CFD的預(yù)報(bào)結(jié)果較實(shí)船試航結(jié)果偏大約2.5%,滿負(fù)荷工況下,預(yù)報(bào)值偏高約3.5%。兩型船基于CFD的航速預(yù)估在考察的工況下較試航結(jié)果偏差均小于0.3kn,該預(yù)報(bào)精度較以往經(jīng)驗(yàn)的方法有了較大提高,對(duì)于初步設(shè)計(jì)、多方案選型階段具有較強(qiáng)的工程實(shí)用性。

    圖3 2750TEU集裝箱船估算結(jié)果與試航結(jié)果比較

    圖4 92500dwt散貨船估算結(jié)果與試航結(jié)果比較

    總體來看,作為定量預(yù)報(bào)其精度尚有欠缺,但該方法較水池試驗(yàn)大大縮短了時(shí)間周期,且可提供諸如波形、船體表面壓力分布等流場(chǎng)信息,能夠?yàn)榇偷膬?yōu)化設(shè)計(jì)提供指導(dǎo)性的信息,同時(shí)又較傳統(tǒng)的圖譜法或回歸公式法預(yù)估的精度有進(jìn)一步的提高,具有較好的工程應(yīng)用前景。

    4 結(jié) 語

    對(duì)基于 CFD的實(shí)船阻力快速預(yù)估方法進(jìn)行了研究,應(yīng)用該方法對(duì)一型細(xì)長(zhǎng)型的集裝箱船以及一型肥大型的散貨船進(jìn)行了阻力估算,通過與實(shí)船試航結(jié)果的比較,得出如下結(jié)論:

    1) 基于CFD的實(shí)船阻力快速預(yù)估方法作為初步設(shè)計(jì)階段多方案選型的定性比較,以及作為航速、功率的初步估算是可行的;

    2) 估算中相關(guān)因子Ca的選擇需要進(jìn)一步研究分析,針對(duì)不同船型結(jié)合實(shí)船試航數(shù)據(jù),積累形成數(shù)據(jù)庫;

    3) 目前估算中自航因子是通過試驗(yàn)以及相關(guān)船型的資料分析得到,下一步需探討通過基于CFD的自航計(jì)算求取自航因子的可行性。

    該方法對(duì)兩型不同特點(diǎn)的船進(jìn)行了阻力預(yù)估,其結(jié)果滿足工程精度要求。為取得進(jìn)一步發(fā)展,基于CFD的實(shí)船阻力快速預(yù)估方法需不斷積累不同船型的算例,并與模型以及實(shí)船試驗(yàn)數(shù)據(jù)進(jìn)行相關(guān)分析,建立不同船型的相關(guān)系數(shù)數(shù)據(jù)庫,提高預(yù)報(bào)結(jié)果的可靠性和精度。

    [1] 姜次平,邵世明. 船舶阻力[M]. 上海:上海交通大學(xué)出版社,1986.

    [2] 陳昌運(yùn),董國(guó)祥,陳霞萍,等. 船舶綜合性能預(yù)報(bào)系統(tǒng)[J]. 上海船舶運(yùn)輸科學(xué)研究所學(xué)報(bào),2007, 30 (2): 81-87.

    [3] 中國(guó)船舶工業(yè)總公司. 船舶設(shè)計(jì)實(shí)用手冊(cè)(總體分冊(cè))[M]. 北京:國(guó)防工業(yè)出版社,1999.

    [4] 沈泓萃,趙 峰. 艦船綜合航行性能虛擬實(shí)驗(yàn)環(huán)境(數(shù)值水池)頂層研究[J]. 艦船科學(xué)技術(shù),2007, 29 (2): 17-21.

    [5] 朱德祥,沈泓萃,洪方文,等. 船模數(shù)值水池框架及其研究基礎(chǔ)[J].水動(dòng)力學(xué)研究與進(jìn)展,2008,A輯, 23 (1): 24-32.

    [6] 李勝忠,李 斌,趙峰,等. VIRTUE計(jì)劃研究進(jìn)展綜述[J]. 船舶力學(xué),2009, 13 (4): 662-675.

    [7] 倪崇本,朱仁傳,繆國(guó)平,等. 一種基于CFD的船舶總阻力預(yù)報(bào)方法[J]. 水動(dòng)力學(xué)研究與進(jìn)展,2010,A輯,25 (5): 579-586.

    [8] 余建偉,朱仁傳,繆國(guó)平,等. 漁政船的阻力計(jì)算與預(yù)報(bào)研究[J]. 水動(dòng)力學(xué)研究與進(jìn)展,2011,A輯,26 (2): 133-139.

    [9] 王 健,李海濤. 計(jì)算流體力學(xué)方法在船舶領(lǐng)域的實(shí)用性研究[J]. 船舶與海洋工程,2012, (4): 6-11.

    猜你喜歡
    實(shí)船船模船型
    船舶阻力試驗(yàn)實(shí)船換算的蒙特卡羅仿真
    基于模糊PID的船模航向控制研究
    長(zhǎng)江1 000車滾裝船EEDI驗(yàn)證分析
    船海工程(2019年6期)2019-12-25 06:12:50
    5 668 TEU集裝箱船球艏改型與浮態(tài)節(jié)能的實(shí)船驗(yàn)證
    船模靜水橫搖試驗(yàn)的不確定度分析
    中國(guó)古代著名的海船船型
    ——福船
    中國(guó)古代著名的海船船型—廣船
    Anti-ageing effects of a new Dimethylaminoethanol-based formulation on DGalactose induced skin ageing model of rat
    西洋船模王——童鑑良
    航海(2014年6期)2014-12-12 10:36:03
    內(nèi)傾船型的非線性橫搖仿真分析
    动漫黄色视频在线观看| 精品久久久精品久久久| 欧美久久黑人一区二区| 曰老女人黄片| 人人妻,人人澡人人爽秒播| 日韩欧美一区二区三区在线观看 | 成人av一区二区三区在线看 | 亚洲精品av麻豆狂野| 老司机深夜福利视频在线观看 | 国产成人欧美在线观看 | 高潮久久久久久久久久久不卡| 97人妻天天添夜夜摸| 国产精品麻豆人妻色哟哟久久| 免费高清在线观看日韩| 18禁黄网站禁片午夜丰满| 精品国产乱子伦一区二区三区 | 最黄视频免费看| 欧美大码av| 一级a爱视频在线免费观看| 日韩 欧美 亚洲 中文字幕| 日韩欧美免费精品| 黄片大片在线免费观看| 啦啦啦中文免费视频观看日本| 久久国产精品影院| 欧美在线黄色| 伊人亚洲综合成人网| 操美女的视频在线观看| 人妻人人澡人人爽人人| 黑人巨大精品欧美一区二区蜜桃| 国产一区二区三区在线臀色熟女 | 国产真人三级小视频在线观看| 国产亚洲精品久久久久5区| 国产欧美日韩综合在线一区二区| 91精品国产国语对白视频| 欧美 日韩 精品 国产| av欧美777| 日韩三级视频一区二区三区| 巨乳人妻的诱惑在线观看| 婷婷色av中文字幕| 制服人妻中文乱码| av一本久久久久| 久久久久精品国产欧美久久久 | 久久精品国产综合久久久| 国产免费视频播放在线视频| 搡老熟女国产l中国老女人| 亚洲精品国产av成人精品| 老司机在亚洲福利影院| 在线 av 中文字幕| 飞空精品影院首页| 精品人妻1区二区| 成人亚洲精品一区在线观看| 国产麻豆69| 亚洲中文日韩欧美视频| 亚洲av成人不卡在线观看播放网 | 热99国产精品久久久久久7| 亚洲国产毛片av蜜桃av| 久久天躁狠狠躁夜夜2o2o| 国产亚洲午夜精品一区二区久久| 亚洲精品成人av观看孕妇| 久久久欧美国产精品| 亚洲精品日韩在线中文字幕| 黄色a级毛片大全视频| 亚洲精品在线美女| 黄色视频不卡| 国产国语露脸激情在线看| 啦啦啦视频在线资源免费观看| 亚洲五月色婷婷综合| 蜜桃在线观看..| 18在线观看网站| 国产一区二区 视频在线| 性色av一级| 日本精品一区二区三区蜜桃| 国产熟女午夜一区二区三区| 丁香六月欧美| 欧美av亚洲av综合av国产av| 中国美女看黄片| 欧美精品av麻豆av| 国产一级毛片在线| 叶爱在线成人免费视频播放| 高清av免费在线| 黄色视频在线播放观看不卡| 一区二区日韩欧美中文字幕| 麻豆国产av国片精品| 日日夜夜操网爽| 国产成人欧美| 午夜福利免费观看在线| 国产极品粉嫩免费观看在线| 欧美人与性动交α欧美软件| 亚洲视频免费观看视频| 久久香蕉激情| 丰满少妇做爰视频| 精品人妻在线不人妻| 又大又爽又粗| 国产成人a∨麻豆精品| 乱人伦中国视频| 狠狠婷婷综合久久久久久88av| 可以免费在线观看a视频的电影网站| 亚洲专区字幕在线| 久久天堂一区二区三区四区| 中文字幕人妻熟女乱码| av不卡在线播放| 日韩精品免费视频一区二区三区| 亚洲精品日韩在线中文字幕| 两个人看的免费小视频| 精品人妻1区二区| 精品免费久久久久久久清纯 | 日韩有码中文字幕| 亚洲欧美一区二区三区久久| 日韩制服丝袜自拍偷拍| 少妇被粗大的猛进出69影院| 日韩欧美一区二区三区在线观看 | 免费在线观看视频国产中文字幕亚洲 | 少妇精品久久久久久久| 深夜精品福利| av网站免费在线观看视频| 日韩中文字幕视频在线看片| 最新在线观看一区二区三区| 青春草亚洲视频在线观看| 欧美国产精品va在线观看不卡| 黑人巨大精品欧美一区二区mp4| 99久久人妻综合| 女人高潮潮喷娇喘18禁视频| 91成人精品电影| 欧美黑人精品巨大| 男女之事视频高清在线观看| 国产精品成人在线| 脱女人内裤的视频| 欧美 亚洲 国产 日韩一| 国产日韩欧美视频二区| h视频一区二区三区| cao死你这个sao货| 亚洲精品国产一区二区精华液| 久久精品国产综合久久久| 丝瓜视频免费看黄片| 精品人妻熟女毛片av久久网站| 亚洲中文av在线| 国产国语露脸激情在线看| 午夜成年电影在线免费观看| 热re99久久国产66热| 嫩草影视91久久| 国产欧美亚洲国产| 久久精品亚洲av国产电影网| 亚洲伊人久久精品综合| 久久精品国产综合久久久| 天天躁夜夜躁狠狠躁躁| 亚洲精品美女久久av网站| 在线观看免费午夜福利视频| 王馨瑶露胸无遮挡在线观看| 精品亚洲成a人片在线观看| 丝袜美腿诱惑在线| 两性夫妻黄色片| av网站在线播放免费| 手机成人av网站| 丁香六月天网| 深夜精品福利| 国产极品粉嫩免费观看在线| 国产成人精品久久二区二区免费| 午夜免费成人在线视频| 亚洲中文av在线| 日韩大码丰满熟妇| 夫妻午夜视频| 欧美日韩黄片免| 国产伦理片在线播放av一区| 色精品久久人妻99蜜桃| 在线av久久热| 伦理电影免费视频| av不卡在线播放| 美女中出高潮动态图| 黄色视频不卡| 制服诱惑二区| 色视频在线一区二区三区| 老熟妇乱子伦视频在线观看 | 日韩三级视频一区二区三区| 免费看十八禁软件| 9色porny在线观看| 免费在线观看影片大全网站| 欧美激情极品国产一区二区三区| 18禁黄网站禁片午夜丰满| 大片免费播放器 马上看| av天堂在线播放| 99精品欧美一区二区三区四区| 麻豆av在线久日| 91老司机精品| 飞空精品影院首页| 亚洲av日韩精品久久久久久密| 最新的欧美精品一区二区| av一本久久久久| 婷婷色av中文字幕| 国产欧美亚洲国产| 日韩有码中文字幕| 超碰成人久久| 国产av一区二区精品久久| 亚洲精品久久久久久婷婷小说| 丰满迷人的少妇在线观看| 老司机在亚洲福利影院| 在线观看免费视频网站a站| 久久精品国产亚洲av高清一级| 一区在线观看完整版| 久久ye,这里只有精品| 丁香六月欧美| 久久久久久人人人人人| 亚洲av电影在线进入| 操美女的视频在线观看| 91九色精品人成在线观看| 亚洲av成人不卡在线观看播放网 | 桃花免费在线播放| 久久国产亚洲av麻豆专区| 每晚都被弄得嗷嗷叫到高潮| 高清视频免费观看一区二区| 一区在线观看完整版| 亚洲精品av麻豆狂野| 老司机靠b影院| 考比视频在线观看| 五月开心婷婷网| 亚洲av国产av综合av卡| 高清在线国产一区| 侵犯人妻中文字幕一二三四区| 国产av一区二区精品久久| 青草久久国产| 淫妇啪啪啪对白视频 | 丝瓜视频免费看黄片| 国产在线一区二区三区精| 国产亚洲欧美在线一区二区| 久久ye,这里只有精品| 97在线人人人人妻| 后天国语完整版免费观看| 无限看片的www在线观看| 99精品欧美一区二区三区四区| 久久久久视频综合| 久久人人97超碰香蕉20202| 啪啪无遮挡十八禁网站| 99久久99久久久精品蜜桃| 他把我摸到了高潮在线观看 | 国产在视频线精品| 19禁男女啪啪无遮挡网站| 国产亚洲欧美精品永久| 国产野战对白在线观看| 视频区图区小说| 99热全是精品| 一级a爱视频在线免费观看| 久久久国产精品麻豆| 欧美精品av麻豆av| 午夜免费观看性视频| 国产成人av教育| 精品少妇一区二区三区视频日本电影| 少妇裸体淫交视频免费看高清 | 国产av精品麻豆| 18禁黄网站禁片午夜丰满| 欧美久久黑人一区二区| 下体分泌物呈黄色| kizo精华| 亚洲色图 男人天堂 中文字幕| 99国产精品99久久久久| 精品人妻一区二区三区麻豆| 国产男女超爽视频在线观看| 国产97色在线日韩免费| 曰老女人黄片| 日本精品一区二区三区蜜桃| 人人妻人人澡人人看| 高清欧美精品videossex| 国产精品香港三级国产av潘金莲| a级毛片在线看网站| 淫妇啪啪啪对白视频 | 老司机影院毛片| 狂野欧美激情性xxxx| 老熟妇乱子伦视频在线观看 | 在线观看免费午夜福利视频| 久久精品亚洲熟妇少妇任你| 国产精品av久久久久免费| 精品国内亚洲2022精品成人 | 亚洲欧美日韩另类电影网站| 69av精品久久久久久 | 免费黄频网站在线观看国产| 国产精品秋霞免费鲁丝片| 欧美成人午夜精品| 中国美女看黄片| 免费日韩欧美在线观看| 免费观看av网站的网址| av电影中文网址| 捣出白浆h1v1| 青青草视频在线视频观看| 久热这里只有精品99| 成人国语在线视频| 精品人妻一区二区三区麻豆| 亚洲国产av新网站| 美女视频免费永久观看网站| 成人手机av| 亚洲专区国产一区二区| a级毛片在线看网站| 超碰成人久久| 成年美女黄网站色视频大全免费| av电影中文网址| 性高湖久久久久久久久免费观看| 久久青草综合色| 亚洲精品一二三| 亚洲精品国产一区二区精华液| 国产亚洲av片在线观看秒播厂| 免费在线观看黄色视频的| 777久久人妻少妇嫩草av网站| 国产高清视频在线播放一区 | 精品熟女少妇八av免费久了| 女性被躁到高潮视频| 色婷婷久久久亚洲欧美| 岛国毛片在线播放| 日韩中文字幕视频在线看片| 国产精品影院久久| 亚洲,欧美精品.| a级毛片在线看网站| 宅男免费午夜| 久久av网站| 天天操日日干夜夜撸| 国产精品一区二区在线不卡| 亚洲欧美清纯卡通| 国产又色又爽无遮挡免| 日韩大码丰满熟妇| 久久ye,这里只有精品| 我要看黄色一级片免费的| 在线亚洲精品国产二区图片欧美| a级片在线免费高清观看视频| 中国美女看黄片| 国产成人欧美在线观看 | 日韩三级视频一区二区三区| 中亚洲国语对白在线视频| 亚洲精品国产色婷婷电影| 一二三四在线观看免费中文在| 少妇精品久久久久久久| 久久人妻熟女aⅴ| 人妻久久中文字幕网| 九色亚洲精品在线播放| 91麻豆精品激情在线观看国产 | 欧美一级毛片孕妇| 精品人妻1区二区| av福利片在线| 欧美 日韩 精品 国产| 免费高清在线观看视频在线观看| 侵犯人妻中文字幕一二三四区| 在线看a的网站| 欧美日韩国产mv在线观看视频| 成年人免费黄色播放视频| 国产成人av激情在线播放| av视频免费观看在线观看| 色综合欧美亚洲国产小说| 一个人免费看片子| 亚洲精品乱久久久久久| 无遮挡黄片免费观看| 午夜老司机福利片| 丝瓜视频免费看黄片| 国产欧美日韩一区二区精品| 久久久久久免费高清国产稀缺| 日韩人妻精品一区2区三区| 热99国产精品久久久久久7| 国产精品自产拍在线观看55亚洲 | 在线观看舔阴道视频| 少妇被粗大的猛进出69影院| 18禁黄网站禁片午夜丰满| 麻豆乱淫一区二区| 久久天躁狠狠躁夜夜2o2o| 岛国在线观看网站| 高清黄色对白视频在线免费看| 一本综合久久免费| 91av网站免费观看| 超碰成人久久| 国产免费福利视频在线观看| 老司机亚洲免费影院| 久久人妻福利社区极品人妻图片| 午夜福利影视在线免费观看| 欧美激情高清一区二区三区| 亚洲自偷自拍图片 自拍| 成在线人永久免费视频| 精品久久久精品久久久| 国产有黄有色有爽视频| 又黄又粗又硬又大视频| 天堂中文最新版在线下载| 午夜福利影视在线免费观看| 俄罗斯特黄特色一大片| 亚洲欧美清纯卡通| 69av精品久久久久久 | 国产精品久久久av美女十八| 老司机午夜福利在线观看视频 | 久久亚洲国产成人精品v| 欧美少妇被猛烈插入视频| av欧美777| √禁漫天堂资源中文www| 国产精品香港三级国产av潘金莲| 日本猛色少妇xxxxx猛交久久| 欧美97在线视频| 久久久久国产一级毛片高清牌| 日本91视频免费播放| e午夜精品久久久久久久| 人人妻,人人澡人人爽秒播| 精品国产一区二区久久| 精品人妻一区二区三区麻豆| 国产97色在线日韩免费| 国产成人系列免费观看| 老司机福利观看| 精品国产一区二区三区四区第35| 国产激情久久老熟女| 成年动漫av网址| 国产男女内射视频| 日韩欧美一区二区三区在线观看 | 中国美女看黄片| 亚洲欧美激情在线| 国产av又大| 国产精品久久久久久精品古装| 嫁个100分男人电影在线观看| 国产视频一区二区在线看| 又大又爽又粗| 丝袜在线中文字幕| bbb黄色大片| 自线自在国产av| 一区二区av电影网| videos熟女内射| 国产区一区二久久| 欧美在线黄色| 男女高潮啪啪啪动态图| 99国产精品免费福利视频| 日本av手机在线免费观看| 十八禁网站网址无遮挡| 国产免费一区二区三区四区乱码| 中文字幕人妻丝袜制服| 曰老女人黄片| 91国产中文字幕| 亚洲精品国产区一区二| 国产欧美日韩综合在线一区二区| 国产精品影院久久| av网站在线播放免费| 一区福利在线观看| 中文字幕最新亚洲高清| 亚洲精品美女久久久久99蜜臀| 久久国产亚洲av麻豆专区| videosex国产| 中文字幕精品免费在线观看视频| 亚洲国产欧美一区二区综合| 一边摸一边做爽爽视频免费| 亚洲天堂av无毛| 亚洲精品国产av蜜桃| 久久这里只有精品19| 精品亚洲成a人片在线观看| 爱豆传媒免费全集在线观看| 午夜福利乱码中文字幕| 久久热在线av| 久久女婷五月综合色啪小说| 我要看黄色一级片免费的| 伊人亚洲综合成人网| 久热爱精品视频在线9| 香蕉国产在线看| 在线亚洲精品国产二区图片欧美| 国产精品亚洲av一区麻豆| 日韩视频在线欧美| 国产人伦9x9x在线观看| 老汉色av国产亚洲站长工具| 亚洲av成人一区二区三| 交换朋友夫妻互换小说| 国产亚洲一区二区精品| 国产欧美日韩综合在线一区二区| 亚洲欧美色中文字幕在线| 熟女少妇亚洲综合色aaa.| 日韩欧美一区二区三区在线观看 | 中文字幕人妻熟女乱码| 极品少妇高潮喷水抽搐| 午夜影院在线不卡| 久久久精品94久久精品| 一二三四社区在线视频社区8| 少妇猛男粗大的猛烈进出视频| 大陆偷拍与自拍| 色婷婷久久久亚洲欧美| 青草久久国产| 99久久人妻综合| 大码成人一级视频| 亚洲精品一二三| 亚洲 欧美一区二区三区| 国产有黄有色有爽视频| 国产亚洲av高清不卡| 免费在线观看黄色视频的| tube8黄色片| 国产精品欧美亚洲77777| 亚洲男人天堂网一区| 真人做人爱边吃奶动态| 中文字幕另类日韩欧美亚洲嫩草| 九色亚洲精品在线播放| 精品国产乱码久久久久久小说| 国产亚洲精品第一综合不卡| 正在播放国产对白刺激| 午夜福利免费观看在线| 岛国在线观看网站| 一区在线观看完整版| 日本五十路高清| 熟女少妇亚洲综合色aaa.| 中文字幕人妻丝袜一区二区| netflix在线观看网站| 国产成人精品久久二区二区91| 激情视频va一区二区三区| 人人妻人人澡人人看| 国产在线观看jvid| 久久影院123| 欧美日韩中文字幕国产精品一区二区三区 | 一级,二级,三级黄色视频| 婷婷成人精品国产| 成人国产av品久久久| 中文字幕人妻丝袜一区二区| 午夜福利影视在线免费观看| 欧美亚洲 丝袜 人妻 在线| 狠狠精品人妻久久久久久综合| 亚洲成国产人片在线观看| 黄片小视频在线播放| 美女福利国产在线| a级毛片在线看网站| 三上悠亚av全集在线观看| 亚洲精品乱久久久久久| 深夜精品福利| 国产亚洲精品一区二区www | 成年美女黄网站色视频大全免费| 麻豆乱淫一区二区| av片东京热男人的天堂| 国产精品久久久人人做人人爽| 亚洲色图综合在线观看| 性色av乱码一区二区三区2| 久久人人爽人人片av| 日韩有码中文字幕| 精品国产一区二区三区四区第35| 精品国内亚洲2022精品成人 | 精品视频人人做人人爽| 一级毛片精品| 天堂中文最新版在线下载| 韩国高清视频一区二区三区| 日本撒尿小便嘘嘘汇集6| 国产又爽黄色视频| 天堂8中文在线网| 黑丝袜美女国产一区| 秋霞在线观看毛片| 日韩制服骚丝袜av| av视频免费观看在线观看| 中文字幕高清在线视频| 一级片免费观看大全| 女人被躁到高潮嗷嗷叫费观| 婷婷丁香在线五月| 91精品三级在线观看| 国产精品熟女久久久久浪| 国产色视频综合| videos熟女内射| 国产一区二区 视频在线| 黑丝袜美女国产一区| 午夜两性在线视频| 两人在一起打扑克的视频| 久久久欧美国产精品| 亚洲成人手机| 麻豆av在线久日| 日韩制服丝袜自拍偷拍| 国产精品欧美亚洲77777| 久久女婷五月综合色啪小说| 国产视频一区二区在线看| 桃花免费在线播放| 日本精品一区二区三区蜜桃| 亚洲自偷自拍图片 自拍| 后天国语完整版免费观看| 欧美av亚洲av综合av国产av| 国产99久久九九免费精品| 亚洲七黄色美女视频| 欧美日韩亚洲高清精品| 亚洲精品国产一区二区精华液| 日韩三级视频一区二区三区| 欧美精品亚洲一区二区| 精品国产一区二区久久| 亚洲成人国产一区在线观看| 欧美 亚洲 国产 日韩一| 亚洲一卡2卡3卡4卡5卡精品中文| 国产主播在线观看一区二区| 9191精品国产免费久久| 十八禁人妻一区二区| 纯流量卡能插随身wifi吗| 国产男女内射视频| 亚洲人成77777在线视频| 秋霞在线观看毛片| 狂野欧美激情性xxxx| 欧美激情极品国产一区二区三区| 9热在线视频观看99| 女警被强在线播放| 考比视频在线观看| 日韩免费高清中文字幕av| 欧美精品一区二区免费开放| 在线观看一区二区三区激情| 免费一级毛片在线播放高清视频 | www.熟女人妻精品国产| 老司机影院毛片| 妹子高潮喷水视频| 日本猛色少妇xxxxx猛交久久| 国产一区二区三区在线臀色熟女 | 亚洲欧美成人综合另类久久久| 蜜桃在线观看..| 亚洲va日本ⅴa欧美va伊人久久 | 久久久久久免费高清国产稀缺| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av日韩精品久久久久久密| 亚洲精品第二区| 可以免费在线观看a视频的电影网站| 久久久久精品人妻al黑| 丝袜脚勾引网站| 亚洲国产av新网站| 黄片小视频在线播放| a级毛片黄视频| 热99re8久久精品国产| 国产一卡二卡三卡精品| 欧美人与性动交α欧美软件| 国产极品粉嫩免费观看在线| 少妇被粗大的猛进出69影院| 亚洲自偷自拍图片 自拍| 日本一区二区免费在线视频| 黄片小视频在线播放| a级毛片黄视频| 操美女的视频在线观看| 亚洲成av片中文字幕在线观看| 97人妻天天添夜夜摸| 国产精品亚洲av一区麻豆| 国产精品.久久久| 精品国产超薄肉色丝袜足j|