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

    逆風(fēng)條件下柔性風(fēng)帆性能的流固耦合研究

    2022-07-05 03:41:22方澤江譚俊哲紀(jì)光英袁鵬孫澤王樹杰
    中國艦船研究 2022年3期
    關(guān)鍵詞:拱度背風(fēng)面風(fēng)帆

    方澤江,譚俊哲*,2,紀(jì)光英,袁鵬,2,孫澤,王樹杰,2

    1 中國海洋大學(xué) 工程學(xué)院,山東 青島 266100

    2 青島市海洋可再生能源重點(diǎn)實(shí)驗(yàn)室,山東 青島 266100

    0 引 言

    風(fēng)帆性能與帆船動力性能直接相關(guān)。柔性風(fēng)帆由薄織物制成,航行時呈現(xiàn)為帶有拱度的形狀,捕獲風(fēng)能,以產(chǎn)生帆船航行的驅(qū)動力。柔性風(fēng)帆在風(fēng)載荷下易發(fā)生非線性大變形,它的初始形狀與航行時的真實(shí)形狀存在很大差異。柔性風(fēng)帆航行時的真實(shí)形狀可簡稱為飛行形狀[1],風(fēng)帆的幾何形狀影響著周圍的空氣流動和產(chǎn)生的力。有必要研究柔性風(fēng)帆的變形并采集其飛行形狀[2],更準(zhǔn)確地預(yù)測柔性風(fēng)帆的空氣動力學(xué)性能。

    計算流體力學(xué)(CFD)能夠幫助工程師完成最佳風(fēng)帆設(shè)計[3]。用CFD 預(yù)測風(fēng)帆性能,絕大多數(shù)以風(fēng)帆初始形狀為形狀模型[4],風(fēng)帆形狀不隨風(fēng)況變化而改變[5-6]。柔性風(fēng)帆變形是非線性大變形問題,與風(fēng)載荷之間的相互作用屬于復(fù)雜的流固耦合現(xiàn)象。Lee 等[7]采用滑動網(wǎng)格方法對游艇風(fēng)帆的二維截面進(jìn)行流固耦合分析,研究風(fēng)帆變形對性能的影響。Bak 等[8]采用單向流固耦合方法研究了無拱度的織物材料的簡易風(fēng)帆變形的非線性行為,如褶皺現(xiàn)象。Bak 等[9-10]將風(fēng)帆視為殼體結(jié)構(gòu),基于集成的CFX和LS-DYNA對30 ft(9.144 m)游艇風(fēng)帆進(jìn)行部分雙向流固耦合分析,比較了初始風(fēng)帆和變形風(fēng)帆的空氣動力學(xué)性能;還研究了風(fēng)帆和桅桿變形對風(fēng)帆空氣動力學(xué)性能的影響。柔性風(fēng)帆變形的流固耦合試驗(yàn)研究也取得了許多進(jìn)展[11-14]。

    柔性風(fēng)帆性能試驗(yàn)成本高昂、試驗(yàn)條件受限;而CFD 數(shù)值模擬成本低、操作方便,但難以準(zhǔn)確解決柔性風(fēng)帆非線性大變形問題。為此,本文將對逆風(fēng)條件下30 ft 游艇主帆的變形進(jìn)行流固耦合研究。采用顯式有限元分析程序AUTODYN解決柔性風(fēng)帆非線性大變形問題,構(gòu)建變形后的柔性風(fēng)帆形狀。采用STAR-CCM+軟件預(yù)測柔性風(fēng)帆空氣動力學(xué)性能,比較柔性風(fēng)帆初始形狀和變形形狀表面的壓力分布、周圍的空氣流動情況、產(chǎn)生的升力和阻力,以及受力中心位置之間的差異,討論柔性風(fēng)帆變形對空氣動力學(xué)性能的影響。

    1 數(shù)值模擬方法

    1.1 流固耦合方法

    采用Euler-Lagrange 作用的流固耦合方法。在全耦合中,拉格朗日網(wǎng)格和歐拉網(wǎng)格進(jìn)行動態(tài)相互作用。流固耦合的結(jié)構(gòu)模型使用需定義人工厚度的殼單元,人工厚度至少是其周圍歐拉單元最小尺寸的2 倍。對于基于有限元方法(FEM)的柔性風(fēng)帆非線性大變形分析,在ANSYS Workbench 中通過集成的顯式動力學(xué)模塊Explicit Dynamics 和AUTODYN 進(jìn)行流固耦合研究。

    1.2 流體運(yùn)動控制

    風(fēng)帆周圍的湍流使用STAR-CCM+通過求解三維Navier-Stokes (N-S)方程組模擬。質(zhì)量守恒方程、動量守恒方程和能量守恒方程分別為[15]:

    對N-S 方程進(jìn)行時間積分,獲得雷諾平均Navier-Stokes (RANS)方程。湍流運(yùn)動計算時,還需附加湍流方程。本文研究采用SSTk-ω 湍流模型,因?yàn)镾STk-ω 考慮了湍流中剪切應(yīng)力的輸運(yùn),可以精確計算在逆壓梯度下從光滑表面流動分離的時刻和數(shù)量,能準(zhǔn)確模擬風(fēng)帆周圍的氣流分離,計算成本比更高階湍流模型低[16]。

    2 柔性風(fēng)帆模型

    2.1 形狀模型

    對Yoo 等[17]設(shè)計的30 ft 游艇主帆進(jìn)行數(shù)值模擬。風(fēng)帆的拱度形狀通過零厚度NACA 6 位數(shù)機(jī)翼截面a=0.8 的中弧線表示,中弧線通過載荷分布線得到,從x/c=0 至x/c=a載荷均勻分布,從x/c=a至x/c=1 載荷線性減小至0,函數(shù)表達(dá)式如式(4)所示[18],Lee 等[7],Bak 等[9],Yoo 等[19],Ciortan等[20], Bak 等[10,21]都采用此風(fēng)帆模型進(jìn)行了數(shù)值模擬或者比例模型試驗(yàn)研究。

    式中:y為中弧線的縱坐標(biāo);c為弦長;x為沿弦長方向的距離;a=0.8,為中弧線類型的標(biāo)記;cli=1,為設(shè)計升力系數(shù);g和l的表達(dá)式為

    NACAa=0.8 的中弧線如圖1 所示,采用5 個機(jī)翼截面的中弧線通過放樣構(gòu)建出風(fēng)帆形狀,5 個截面從風(fēng)帆底邊沿高度方向依次命名為foot,25%,50%,75%,top,風(fēng)帆幾何參數(shù)如表1 所示,幾何形狀如圖2 所示,帆高11.9 m,面積24.7 m2。

    圖1 NACA a=0.8 的中弧線Fig. 1 Mean camber line with NACA a=0.8

    表1 風(fēng)帆幾何參數(shù)Table 1 Geometrical parameters of the sail

    圖2 風(fēng)帆形狀模型Fig. 2 Shape model of the sail

    2.2 材料模型

    常用薄織物作為柔性風(fēng)帆材料,假設(shè)材料為各向同性(Isotropic)彈性材料,材料模型采用密度ρ、彈性模量E和泊松比μ來表征。Kevlar?49 織物的材料特性如表2 所示[9]。

    3 柔性風(fēng)帆的流固耦合分析

    3.1 流固耦合分析的條件設(shè)置

    基于AUTODYN 軟件采用Euler-Lagrange 作用的全耦合方法,設(shè)置以下條件進(jìn)行柔性風(fēng)帆非線性大變形分析。

    1) 以圖2 所示風(fēng)帆形狀作為柔性風(fēng)帆的初始形狀,表2 所示Kevlar?49 織物作為柔性風(fēng)帆材料。

    2) 為了保證較高的流固耦合精度,流體域的網(wǎng)格尺寸應(yīng)不大于結(jié)構(gòu)的網(wǎng)格尺寸。結(jié)構(gòu)的基本網(wǎng)格尺寸設(shè)置為60 mm,網(wǎng)格單元數(shù)量為6 988,計算域的網(wǎng)格數(shù)量為511 632。

    3) 相對風(fēng)向角為20°,風(fēng)速為10.3 m/s。

    4) 終止時間設(shè)置:流體域長度為4.3 m,風(fēng)速10.3 m/s,時間t=4.3/10.3=0.417 s,考慮風(fēng)和帆的耦合作用需要耗費(fèi)時間,所以取終止時間為0.5 s。

    5) 柔性風(fēng)帆的固定方案對變形具有很大影響,實(shí)際上主前緣固定在桅桿上,后角由橫桿末端的拉力支撐,因此在前緣和后角添加固定邊界條件約束所有移動,以近似代替桅桿和橫桿對主帆的固定和支撐,如圖3 所示。

    圖3 固定邊界條件Fig. 3 Fixed support conditions

    6) 風(fēng)帆被視為殼體結(jié)構(gòu),定義人工殼單元厚度為0.68 mm。

    7) 在風(fēng)帆的5 個截面上添加探測點(diǎn),如圖4所示,通過輸出探測點(diǎn)的空間坐標(biāo)數(shù)據(jù),近似構(gòu)建風(fēng)帆變形的形狀。

    圖4 帆面上的探測點(diǎn)Fig. 4 Gauges on sail surface

    3.2 流固耦合結(jié)果分析

    圖5 所示為風(fēng)帆變形的位移云圖,最大變形約為26 mm,在沒有固定邊界條件的風(fēng)帆邊緣變形較大。圖6 所示為風(fēng)帆變形前、后截面的拱度曲線,可見風(fēng)帆前緣附近拱度向背風(fēng)面變形更深,風(fēng)帆后緣附近帆高較低截面的拱度向背風(fēng)面靠近,帆高較高截面的拱度向迎風(fēng)面靠近,風(fēng)帆在風(fēng)載荷作用下發(fā)生明顯扭曲現(xiàn)象。風(fēng)帆初始形狀面積為24.7 m2,變形形狀面積為24.3 m2,面積減少約1.6%,說明風(fēng)帆面積變化很小,風(fēng)帆在風(fēng)載荷下的拉伸對風(fēng)帆變形的影響可以忽略不計。

    圖5 風(fēng)帆變形的位移云圖Fig. 5 Displacement contours of sail deformation

    圖6 初始風(fēng)帆和變形風(fēng)帆拱度曲線比較Fig. 6 A comparison on camber lines of the initial and deformed sails

    4 柔性風(fēng)帆的性能預(yù)測

    4.1 STAR-CCM+計算的條件設(shè)置

    為了說明柔性風(fēng)帆變形對其空氣動力學(xué)性能的影響,采用STAR-CCM+在相同條件下對初始風(fēng)帆和變形風(fēng)帆進(jìn)行性能預(yù)測。

    1) 風(fēng)帆被視為計算域內(nèi)零厚度的內(nèi)部面。

    2) 假設(shè)溫度為25℃,空氣密度為1.169 kg/m3,湍流邊界層中第1 層網(wǎng)格與帆面之間的距離y+=118。風(fēng)為縱向風(fēng),改變風(fēng)帆位置,使相對風(fēng)向角為20°,風(fēng)速設(shè)置為10.3 m/s。

    3) 采用多面體網(wǎng)格,加密風(fēng)帆附近的網(wǎng)格并設(shè)置邊界層網(wǎng)格(圖7),初始風(fēng)帆和變形風(fēng)帆的網(wǎng)格數(shù)量分別為36 893 和35 671,計算域的網(wǎng)格數(shù)量分別為3 557 850 和3 535 506。

    圖7 風(fēng)帆周圍網(wǎng)格加密和邊界層Fig. 7 Refinements of the meshes around the sail and prism layers

    4) 采用RANS 方程和k-ω SST 湍流模型。

    4.2 性能預(yù)測結(jié)果分析

    比較初始風(fēng)帆和變形風(fēng)帆表面的壓力分布、周圍的空氣流動情況、產(chǎn)生的升力和阻力、受力中心位置的變化。

    圖8 所示為初始風(fēng)帆和變形風(fēng)帆表面的壓力分布,變形風(fēng)帆迎風(fēng)面壓力值有所增大,前緣附近正壓值的范圍也有所擴(kuò)大,這解釋了如圖6 所示的風(fēng)帆前緣附近拱度向背風(fēng)面靠近。圖9 所示為初始風(fēng)帆在帆高25%截面處的壓力分布和速度流線,在背風(fēng)面后緣附近出現(xiàn)流動分離現(xiàn)象,流動分離導(dǎo)致失速使得背風(fēng)面后緣附近壓力值升高,約為?30 Pa,而迎風(fēng)面后緣附近壓力值約為5 Pa,因此帆高25%截面后緣附近的壓差為正,這解釋了如圖6 所示的帆高25% 截面后緣附近的拱度仍然向背風(fēng)面靠近。圖10 所示為初始風(fēng)帆和變形風(fēng)帆在帆高25%截面的速度分布,變形風(fēng)帆較初始風(fēng)帆在背風(fēng)面后緣附近產(chǎn)生的流動分離現(xiàn)象更嚴(yán)重。

    圖8 初始風(fēng)帆和變形風(fēng)帆的帆面壓力分布Fig. 8 Pressure distribution on the surfaces of the initial and deformed sails

    圖9 初始風(fēng)帆在帆高25%截面的壓力分布和速度流線Fig. 9 Pressure distribution and velocity lines around the initial sail on the section of 25% sail height

    圖10 初始風(fēng)帆和變形風(fēng)帆在帆高25%截面的速度分布Fig. 10 Velocity distribution around the initial and deformed sails on the section of 25% sail height

    通過比較初始風(fēng)帆和變形風(fēng)帆的升力系數(shù)(CL)和阻力系數(shù)(CD),分析風(fēng)帆變形對空氣動力學(xué)性能的影響。由表3 可知,與初始風(fēng)帆相比,變形風(fēng)帆的升力系數(shù)減小約0.98%,阻力系數(shù)增大約2.7%,表征空氣動力學(xué)性能的升阻比減小約3.6%。該結(jié)果與Bak 等[9]通過流固耦合研究的三維風(fēng)帆變形對性能的影響效果一致。分析表明,產(chǎn)生上述現(xiàn)象的主要原因在于:帆高較低的帆面所占面積較大,是決定風(fēng)帆性能的主要部分,在相對風(fēng)向角為20°的逆風(fēng)條件下,變形風(fēng)帆較低帆面的背風(fēng)面后緣附近產(chǎn)生了更嚴(yán)重的流動分離現(xiàn)象(圖10),較低帆面后緣附近的拱度向背風(fēng)側(cè)靠近(圖6),導(dǎo)致有效攻角減小。因此在該條件下拱度增大帶來的升力增大效應(yīng)相對較小,總體升力減小。

    表3 初始風(fēng)帆與變形風(fēng)帆的升、阻力系數(shù)Table 3 Lift and drag coefficients of the initial and deformed sails

    帆船在穩(wěn)定航行時,橫傾力矩等于扶正力矩,風(fēng)帆變形時其受力中心位置會發(fā)生變化,從而使帆船失去平衡。假設(shè)風(fēng)帆橫向投影面的質(zhì)心為風(fēng)帆受力中心[22],采用SolidWorks 評估風(fēng)帆初始形狀和變形形狀模型的質(zhì)心,變形風(fēng)帆的受力中心向右下方輕微移動。

    5 結(jié) 論

    本文采用非線性有限元方法計算柔性風(fēng)帆在風(fēng)載荷作用下的變形,以固定邊界條件代替桅桿和橫桿對風(fēng)帆的支撐和固定作用,在帆面上添加探測點(diǎn),通過輸出探測點(diǎn)的空間坐標(biāo)數(shù)據(jù)構(gòu)建變形風(fēng)帆的形狀,比較初始風(fēng)帆和變形風(fēng)帆的形狀特性;采用STAR-CCM+在相同條件下預(yù)測初始風(fēng)帆和變形風(fēng)帆的性能。主要結(jié)論如下:

    1) 變形風(fēng)帆較初始風(fēng)帆面積減小約1.6%,說明材料本身的拉伸對拱度變化的影響可忽略不計。

    2) 變形風(fēng)帆在帆面前緣附近拱度向背風(fēng)面靠近;在后緣附近帆高較低截面的拱度向背風(fēng)面靠近,帆高較高截面的拱度向迎風(fēng)面靠近,有明顯扭曲現(xiàn)象。

    3) 與初始風(fēng)帆相比,變形風(fēng)帆的迎風(fēng)面整體壓力值增大。

    4) 與初始風(fēng)帆相比,變形風(fēng)帆的升力系數(shù)減小約0.98%,阻力系數(shù)增大約2.7%,升阻比減小約3.6%。主要原因在于,變形風(fēng)帆較低帆面的背風(fēng)面后緣附近產(chǎn)生了更嚴(yán)重的流動分離現(xiàn)象,導(dǎo)致失速,使得分離區(qū)域的壓強(qiáng)增加,從而使風(fēng)帆兩側(cè)的壓強(qiáng)差降低,進(jìn)而引起升力減小;較低帆面后緣附近的拱度向背風(fēng)面靠近,造成有效攻角減小,因此風(fēng)帆的整體升力減小。

    5) 變形風(fēng)帆的受力中心向右下方輕微移動。

    6) 柔性風(fēng)帆變形的流固耦合研究適用于精確評估或者優(yōu)化風(fēng)帆性能。

    猜你喜歡
    拱度背風(fēng)面風(fēng)帆
    基于數(shù)值模擬研究預(yù)應(yīng)力T 梁橋上拱度的影響因素
    妹是風(fēng)帆哥是船
    心聲歌刊(2023年3期)2023-07-16 02:11:48
    起重機(jī)主梁上拱度的調(diào)整方法
    大型起重機(jī)主梁拱度研究
    潔白風(fēng)帆的集合體
    讓折桅的船兒重?fù)P風(fēng)帆
    非均勻等離子體Ka-Band傳輸性能中繼法優(yōu)化研究
    高超聲速風(fēng)洞子母彈大迎角拋殼投放試驗(yàn)
    高壓輸電鐵塔塔身背風(fēng)面風(fēng)荷載遮擋效應(yīng)研究
    靜靜的漁港 遠(yuǎn)航的風(fēng)帆
    中國火炬(2014年10期)2014-07-25 10:36:42
    国产男靠女视频免费网站| 精品久久久久久电影网| 最新的欧美精品一区二区| 波多野结衣一区麻豆| 国产熟女午夜一区二区三区| 欧美精品一区二区大全| 一级片免费观看大全| a在线观看视频网站| 国产黄频视频在线观看| 最黄视频免费看| 午夜精品久久久久久毛片777| 精品一区二区三区四区五区乱码| 国产日韩一区二区三区精品不卡| 欧美乱码精品一区二区三区| 少妇 在线观看| 成人国产一区最新在线观看| 十八禁高潮呻吟视频| 丁香欧美五月| 一进一出好大好爽视频| 亚洲熟妇熟女久久| 国产男女超爽视频在线观看| 91av网站免费观看| 这个男人来自地球电影免费观看| 午夜精品久久久久久毛片777| 桃花免费在线播放| 国产激情久久老熟女| 人人澡人人妻人| 99九九在线精品视频| 欧美精品亚洲一区二区| 国产高清视频在线播放一区| 人人妻人人澡人人爽人人夜夜| 美女主播在线视频| 亚洲av片天天在线观看| 亚洲欧美色中文字幕在线| 黄色视频,在线免费观看| 18禁黄网站禁片午夜丰满| 国产一区二区 视频在线| 亚洲精品成人av观看孕妇| 人妻 亚洲 视频| 大型av网站在线播放| 欧美一级毛片孕妇| 精品少妇一区二区三区视频日本电影| 18禁国产床啪视频网站| 国产欧美日韩一区二区精品| 精品少妇内射三级| 精品人妻1区二区| 成年人免费黄色播放视频| 香蕉国产在线看| 国产精品一区二区免费欧美| 中文字幕制服av| 日韩三级视频一区二区三区| 国产成人av教育| 啦啦啦视频在线资源免费观看| 亚洲国产欧美日韩在线播放| 人妻 亚洲 视频| 俄罗斯特黄特色一大片| 美女主播在线视频| 夫妻午夜视频| 男男h啪啪无遮挡| 亚洲色图av天堂| 超碰97精品在线观看| 伊人久久大香线蕉亚洲五| 他把我摸到了高潮在线观看 | 久久精品91无色码中文字幕| 国产欧美日韩一区二区精品| 大码成人一级视频| 亚洲情色 制服丝袜| 丰满少妇做爰视频| 亚洲五月色婷婷综合| 国产在线一区二区三区精| 国产av一区二区精品久久| 亚洲一码二码三码区别大吗| 菩萨蛮人人尽说江南好唐韦庄| 人人澡人人妻人| 超碰97精品在线观看| 亚洲av美国av| 热99re8久久精品国产| 色播在线永久视频| 国产精品香港三级国产av潘金莲| 黄片大片在线免费观看| 美女福利国产在线| 国产又爽黄色视频| 国产成人精品久久二区二区91| 人妻久久中文字幕网| 欧美日韩av久久| 人妻 亚洲 视频| 久久久久国内视频| 中文字幕av电影在线播放| 午夜视频精品福利| 国产精品av久久久久免费| 夜夜夜夜夜久久久久| 日韩免费高清中文字幕av| 热re99久久精品国产66热6| 国产精品亚洲av一区麻豆| 99久久国产精品久久久| av网站在线播放免费| 黑人欧美特级aaaaaa片| 欧美久久黑人一区二区| 老司机深夜福利视频在线观看| 亚洲欧美激情在线| 新久久久久国产一级毛片| 精品人妻在线不人妻| 午夜福利免费观看在线| 亚洲一区二区三区欧美精品| 久久这里只有精品19| 天堂8中文在线网| 热99国产精品久久久久久7| 久久久久久人人人人人| 欧美变态另类bdsm刘玥| 国产男女超爽视频在线观看| kizo精华| 国产精品av久久久久免费| 久久精品91无色码中文字幕| 无人区码免费观看不卡 | 在线观看www视频免费| 新久久久久国产一级毛片| 国产精品免费一区二区三区在线 | 精品乱码久久久久久99久播| 性色av乱码一区二区三区2| 男人操女人黄网站| av天堂久久9| 国产精品欧美亚洲77777| 亚洲七黄色美女视频| 欧美日韩一级在线毛片| 法律面前人人平等表现在哪些方面| 高清欧美精品videossex| 12—13女人毛片做爰片一| 高清视频免费观看一区二区| 久久久久久久精品吃奶| 亚洲成国产人片在线观看| 精品久久久精品久久久| 精品一区二区三区视频在线观看免费 | 久久精品亚洲熟妇少妇任你| 丰满少妇做爰视频| 纯流量卡能插随身wifi吗| videosex国产| 亚洲伊人久久精品综合| 黄色毛片三级朝国网站| 午夜精品国产一区二区电影| avwww免费| 天天添夜夜摸| 亚洲天堂av无毛| 国产一区二区 视频在线| 一级黄色大片毛片| 国产真人三级小视频在线观看| 亚洲第一青青草原| 久久久久久久精品吃奶| 欧美一级毛片孕妇| 精品久久久精品久久久| www.自偷自拍.com| 亚洲熟妇熟女久久| 久久久国产成人免费| 久久午夜综合久久蜜桃| 亚洲欧美日韩另类电影网站| 日韩 欧美 亚洲 中文字幕| 黑丝袜美女国产一区| 国产麻豆69| 黄色a级毛片大全视频| 热99国产精品久久久久久7| 中文字幕制服av| 国产伦理片在线播放av一区| 亚洲欧美一区二区三区黑人| 精品国产一区二区久久| 制服人妻中文乱码| 97人妻天天添夜夜摸| 一个人免费看片子| 亚洲av电影在线进入| 国产成人影院久久av| 日韩中文字幕视频在线看片| 成人永久免费在线观看视频 | 夫妻午夜视频| 日本五十路高清| 天堂中文最新版在线下载| 亚洲成国产人片在线观看| 亚洲人成电影观看| 免费看a级黄色片| 一进一出好大好爽视频| 久久精品国产亚洲av高清一级| 黄色片一级片一级黄色片| 亚洲中文日韩欧美视频| 老汉色∧v一级毛片| 美女午夜性视频免费| 亚洲精品一二三| 日韩视频一区二区在线观看| 99热国产这里只有精品6| 最新美女视频免费是黄的| 久久精品亚洲精品国产色婷小说| 好男人电影高清在线观看| 中文字幕人妻丝袜一区二区| 一级片免费观看大全| 777米奇影视久久| 少妇猛男粗大的猛烈进出视频| 免费在线观看影片大全网站| 日本精品一区二区三区蜜桃| 久久这里只有精品19| 777久久人妻少妇嫩草av网站| 国产在线免费精品| 久久影院123| 自线自在国产av| 老鸭窝网址在线观看| 国产精品电影一区二区三区 | 日本精品一区二区三区蜜桃| 午夜久久久在线观看| 飞空精品影院首页| 又大又爽又粗| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲精品成人av观看孕妇| 90打野战视频偷拍视频| 黄色怎么调成土黄色| 免费少妇av软件| 一区在线观看完整版| 99国产精品一区二区三区| 日本精品一区二区三区蜜桃| 成人18禁在线播放| 日本欧美视频一区| 69精品国产乱码久久久| 又紧又爽又黄一区二区| 久久中文字幕一级| 久久国产亚洲av麻豆专区| 国产高清激情床上av| 热99re8久久精品国产| 久久国产精品影院| 在线看a的网站| 免费在线观看日本一区| 蜜桃在线观看..| e午夜精品久久久久久久| 精品国产亚洲在线| 日韩人妻精品一区2区三区| 精品少妇内射三级| 捣出白浆h1v1| 国产一区二区在线观看av| 我要看黄色一级片免费的| 亚洲成av片中文字幕在线观看| 久久久久国产一级毛片高清牌| 12—13女人毛片做爰片一| 人成视频在线观看免费观看| 一区二区三区乱码不卡18| 亚洲av美国av| 黄色丝袜av网址大全| 亚洲少妇的诱惑av| 国产精品美女特级片免费视频播放器 | 成在线人永久免费视频| xxxhd国产人妻xxx| 国产精品免费一区二区三区在线 | 久久热在线av| 国产精品久久久人人做人人爽| 国产麻豆69| 激情视频va一区二区三区| 免费女性裸体啪啪无遮挡网站| 多毛熟女@视频| 9热在线视频观看99| 欧美黄色片欧美黄色片| 51午夜福利影视在线观看| 成人亚洲精品一区在线观看| 免费少妇av软件| 精品久久久久久久毛片微露脸| 91国产中文字幕| 熟女少妇亚洲综合色aaa.| 色老头精品视频在线观看| 搡老岳熟女国产| 满18在线观看网站| 国产97色在线日韩免费| 亚洲一区中文字幕在线| 一区福利在线观看| 男女边摸边吃奶| 国产一区二区 视频在线| 国产成人欧美| 视频区欧美日本亚洲| 别揉我奶头~嗯~啊~动态视频| 制服人妻中文乱码| 男人舔女人的私密视频| 人人妻人人爽人人添夜夜欢视频| 亚洲一卡2卡3卡4卡5卡精品中文| 91九色精品人成在线观看| 十八禁网站免费在线| √禁漫天堂资源中文www| 免费看a级黄色片| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩黄片免| 电影成人av| 久久精品亚洲精品国产色婷小说| 亚洲精品粉嫩美女一区| 91老司机精品| 久久精品国产亚洲av香蕉五月 | 国产精品一区二区在线不卡| 久久久精品免费免费高清| 狂野欧美激情性xxxx| 好男人电影高清在线观看| 欧美乱妇无乱码| 亚洲色图av天堂| 香蕉丝袜av| 91精品国产国语对白视频| 亚洲男人天堂网一区| 久久久国产一区二区| 极品少妇高潮喷水抽搐| 色婷婷av一区二区三区视频| 99久久人妻综合| 久久香蕉激情| 不卡一级毛片| 亚洲熟女精品中文字幕| 999久久久精品免费观看国产| 国产精品久久久久久人妻精品电影 | 宅男免费午夜| www.熟女人妻精品国产| 考比视频在线观看| 色尼玛亚洲综合影院| 久久国产精品男人的天堂亚洲| av片东京热男人的天堂| 国产野战对白在线观看| 国产97色在线日韩免费| 丝袜喷水一区| 成人国语在线视频| 无遮挡黄片免费观看| 一区二区av电影网| 亚洲伊人色综图| 日本一区二区免费在线视频| 18禁观看日本| 99热国产这里只有精品6| 国产精品久久久久久精品古装| 中国美女看黄片| 亚洲色图av天堂| 亚洲国产精品一区二区三区在线| 亚洲va日本ⅴa欧美va伊人久久| 黄色视频,在线免费观看| 一区二区三区国产精品乱码| 女性生殖器流出的白浆| 亚洲色图综合在线观看| 国产亚洲欧美在线一区二区| 啦啦啦视频在线资源免费观看| 十八禁网站网址无遮挡| 国产av国产精品国产| 99国产精品免费福利视频| 成人18禁高潮啪啪吃奶动态图| 精品人妻在线不人妻| 欧美另类亚洲清纯唯美| 国产xxxxx性猛交| 亚洲成国产人片在线观看| 国产野战对白在线观看| 可以免费在线观看a视频的电影网站| 男女午夜视频在线观看| 国产精品亚洲av一区麻豆| 91成人精品电影| 亚洲国产av影院在线观看| 999久久久国产精品视频| 欧美人与性动交α欧美精品济南到| 精品国产亚洲在线| 久9热在线精品视频| 婷婷成人精品国产| 美女主播在线视频| 黑人巨大精品欧美一区二区蜜桃| 色婷婷av一区二区三区视频| 丝袜喷水一区| av网站在线播放免费| 日韩视频在线欧美| 国产成人啪精品午夜网站| 久久精品熟女亚洲av麻豆精品| 91麻豆av在线| 国产高清激情床上av| 国产成人啪精品午夜网站| 日韩精品免费视频一区二区三区| www.熟女人妻精品国产| 亚洲av美国av| 日日爽夜夜爽网站| 色94色欧美一区二区| 色播在线永久视频| 中文字幕人妻丝袜制服| 成年人午夜在线观看视频| 99国产精品一区二区蜜桃av | 久久精品aⅴ一区二区三区四区| 美女扒开内裤让男人捅视频| 国产精品国产av在线观看| 亚洲欧洲日产国产| 国产不卡av网站在线观看| 多毛熟女@视频| 午夜福利欧美成人| 操出白浆在线播放| 一边摸一边做爽爽视频免费| 叶爱在线成人免费视频播放| 亚洲精品av麻豆狂野| 国产不卡av网站在线观看| 夜夜夜夜夜久久久久| 九色亚洲精品在线播放| 国产一区二区激情短视频| 人妻久久中文字幕网| 精品国产超薄肉色丝袜足j| 91av网站免费观看| 国产99久久九九免费精品| 巨乳人妻的诱惑在线观看| 别揉我奶头~嗯~啊~动态视频| 欧美 亚洲 国产 日韩一| 欧美黑人欧美精品刺激| 亚洲成国产人片在线观看| 建设人人有责人人尽责人人享有的| 俄罗斯特黄特色一大片| 最新美女视频免费是黄的| 亚洲精品一卡2卡三卡4卡5卡| 欧美日本中文国产一区发布| 国产激情久久老熟女| 成人国产一区最新在线观看| 精品熟女少妇八av免费久了| 久久久久久免费高清国产稀缺| 纵有疾风起免费观看全集完整版| 中文字幕高清在线视频| 黄片小视频在线播放| 亚洲国产av新网站| 精品国产亚洲在线| 交换朋友夫妻互换小说| 色播在线永久视频| 日本精品一区二区三区蜜桃| 美女扒开内裤让男人捅视频| videos熟女内射| av福利片在线| 少妇 在线观看| 男女之事视频高清在线观看| 国精品久久久久久国模美| 在线观看免费视频网站a站| 99久久精品国产亚洲精品| 婷婷成人精品国产| 18在线观看网站| 超碰成人久久| 午夜福利免费观看在线| 人妻久久中文字幕网| 999久久久精品免费观看国产| a级毛片在线看网站| 久久人人爽av亚洲精品天堂| 12—13女人毛片做爰片一| 这个男人来自地球电影免费观看| 99国产精品一区二区三区| videosex国产| 精品一区二区三区av网在线观看 | 久久精品91无色码中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 18在线观看网站| 国产精品国产高清国产av | 男女边摸边吃奶| 久久精品亚洲熟妇少妇任你| 丝袜美腿诱惑在线| 亚洲av片天天在线观看| 国产亚洲欧美在线一区二区| 亚洲情色 制服丝袜| 欧美黑人精品巨大| 精品久久蜜臀av无| 亚洲少妇的诱惑av| 女警被强在线播放| 欧美 日韩 精品 国产| 国产精品偷伦视频观看了| 波多野结衣av一区二区av| 色婷婷av一区二区三区视频| 在线观看免费视频网站a站| 欧美久久黑人一区二区| 亚洲av日韩精品久久久久久密| 一本久久精品| 亚洲精品自拍成人| 成人av一区二区三区在线看| 最新在线观看一区二区三区| 国产免费av片在线观看野外av| 青青草视频在线视频观看| 性高湖久久久久久久久免费观看| 99国产精品一区二区三区| 一边摸一边抽搐一进一小说 | 美女高潮到喷水免费观看| 国产在视频线精品| 电影成人av| 美女视频免费永久观看网站| 十分钟在线观看高清视频www| h视频一区二区三区| 天堂8中文在线网| 午夜福利在线免费观看网站| 成年版毛片免费区| 免费观看a级毛片全部| 国产男靠女视频免费网站| 亚洲国产中文字幕在线视频| 国产一区二区三区综合在线观看| 我要看黄色一级片免费的| 免费看十八禁软件| 嫩草影视91久久| 自拍欧美九色日韩亚洲蝌蚪91| 久久中文看片网| 国产成人精品久久二区二区免费| 精品福利永久在线观看| 欧美 日韩 精品 国产| 18禁裸乳无遮挡动漫免费视频| 另类亚洲欧美激情| 亚洲中文日韩欧美视频| 天天操日日干夜夜撸| 一区二区三区激情视频| e午夜精品久久久久久久| 两个人免费观看高清视频| 窝窝影院91人妻| 51午夜福利影视在线观看| 亚洲熟女毛片儿| 国产欧美日韩综合在线一区二区| 国产真人三级小视频在线观看| 黑人欧美特级aaaaaa片| 国产有黄有色有爽视频| 亚洲成人免费电影在线观看| 日本撒尿小便嘘嘘汇集6| 国产av一区二区精品久久| 亚洲精品av麻豆狂野| 久久人妻熟女aⅴ| 黄网站色视频无遮挡免费观看| 国产成人一区二区三区免费视频网站| 在线观看免费视频网站a站| 精品国内亚洲2022精品成人 | 欧美精品一区二区免费开放| 亚洲成人免费av在线播放| 久久午夜综合久久蜜桃| 国产av精品麻豆| 久热爱精品视频在线9| 日本一区二区免费在线视频| 精品一品国产午夜福利视频| 国产精品电影一区二区三区 | 国产精品久久电影中文字幕 | 国产福利在线免费观看视频| 老司机靠b影院| 精品久久蜜臀av无| 99精品欧美一区二区三区四区| 国产在视频线精品| 两性夫妻黄色片| 涩涩av久久男人的天堂| 91大片在线观看| 亚洲av电影在线进入| 国产精品免费视频内射| 国产一区二区在线观看av| 在线永久观看黄色视频| 王馨瑶露胸无遮挡在线观看| 又大又爽又粗| 精品亚洲成a人片在线观看| 国产精品一区二区在线不卡| 久久中文字幕一级| 亚洲精品在线美女| 中文字幕色久视频| 精品少妇黑人巨大在线播放| 亚洲va日本ⅴa欧美va伊人久久| 香蕉国产在线看| 国产精品麻豆人妻色哟哟久久| 久久人妻av系列| 免费观看人在逋| 亚洲精品国产区一区二| 色综合欧美亚洲国产小说| 亚洲熟妇熟女久久| 老司机在亚洲福利影院| 18禁黄网站禁片午夜丰满| 久久久精品免费免费高清| 亚洲av片天天在线观看| 国产成人影院久久av| 国产野战对白在线观看| 淫妇啪啪啪对白视频| 久久精品国产亚洲av高清一级| 飞空精品影院首页| 亚洲avbb在线观看| 午夜福利免费观看在线| 国产日韩欧美视频二区| 国产成人欧美在线观看 | 欧美日韩一级在线毛片| 久热这里只有精品99| 亚洲国产欧美网| 十八禁人妻一区二区| 亚洲美女黄片视频| 两个人看的免费小视频| 免费观看av网站的网址| 国产高清视频在线播放一区| 国产亚洲午夜精品一区二区久久| 免费高清在线观看日韩| 精品国产超薄肉色丝袜足j| 亚洲国产欧美在线一区| 五月开心婷婷网| 黑人巨大精品欧美一区二区mp4| 桃花免费在线播放| 精品少妇黑人巨大在线播放| 9热在线视频观看99| 亚洲国产av新网站| 亚洲男人天堂网一区| www日本在线高清视频| 69精品国产乱码久久久| 一区二区日韩欧美中文字幕| 在线观看舔阴道视频| 国产黄频视频在线观看| 精品国产一区二区三区四区第35| 一区二区三区国产精品乱码| 嫁个100分男人电影在线观看| 国产在线观看jvid| 真人做人爱边吃奶动态| 十八禁人妻一区二区| 天天影视国产精品| 午夜福利视频精品| 桃红色精品国产亚洲av| 日韩 欧美 亚洲 中文字幕| 热re99久久国产66热| 在线观看免费日韩欧美大片| 男女下面插进去视频免费观看| 18禁观看日本| 国产aⅴ精品一区二区三区波| 国产伦人伦偷精品视频| 纯流量卡能插随身wifi吗| 欧美精品av麻豆av| 亚洲欧美一区二区三区黑人| 久久人人97超碰香蕉20202| 国产日韩一区二区三区精品不卡| 满18在线观看网站| 亚洲伊人色综图| 人人妻人人添人人爽欧美一区卜| 夜夜爽天天搞| 亚洲成人免费av在线播放| 亚洲色图综合在线观看| 18在线观看网站| 国产真人三级小视频在线观看| 色综合婷婷激情| tocl精华| 国产精品久久久人人做人人爽| 国产精品久久久久成人av| 欧美日韩一级在线毛片| 国产精品免费一区二区三区在线 | 婷婷成人精品国产| 男女高潮啪啪啪动态图|