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

    降落傘穩(wěn)降階段的SPH方法數(shù)值模擬

    2016-10-14 08:25:50姚向茹余莉吳瓊
    航天返回與遙感 2016年3期
    關(guān)鍵詞:傘衣氣動力降落傘

    姚向茹余莉吳瓊

    ?

    降落傘穩(wěn)降階段的SPH方法數(shù)值模擬

    姚向茹余莉吳瓊

    (南京航空航天大學航空宇航學院,南京 210016)

    隨著計算機技術(shù)的不斷發(fā)展以及數(shù)值模擬技術(shù)的興起,對降落傘力學特性的模擬日益精確。為了消除傳統(tǒng)降落傘模擬中因網(wǎng)格大變形而產(chǎn)生負體積的影響,文章采用光滑粒子流體動力學方法,對選定降落傘模型開展了穩(wěn)降階段的流固耦合數(shù)值模擬,計算得到的氣動力系數(shù)與風洞試驗結(jié)果較為一致。研究了空氣粒子密度對數(shù)值結(jié)果的影響,結(jié)果表明,粒子密度過少,降落傘結(jié)構(gòu)單元難以和流場實現(xiàn)有效耦合,計算精度降低;但粒子密度超過一定值,不僅會增加計算消耗,也會增加舍入誤差而影響計算精度,確定了最佳粒子數(shù)。在此基礎(chǔ)上,模擬了不同初速度下降落傘的下降過程,獲得了降落傘外形、投影直徑、下降速度、降落傘氣動力及傘衣擺角的動態(tài)變化情況。數(shù)值結(jié)果表明傘衣的呼吸頻率、穩(wěn)定后的平衡速度和氣動力受速度影響較小,但速度增加,傘衣呼吸強度及擺角增加,且擺動角度的衰減性要低于投影直徑的衰減能力。SPH模擬方法有效地回避了降落傘流固耦合的網(wǎng)格大變形問題,在柔性大變形流固耦合問題上有廣闊的應(yīng)用前景。

    下降過程 流固耦合 光滑粒子流體動力學方法 降落傘

    0 引言

    降落傘系統(tǒng)作為一種高效、可靠的回收裝置,在載人飛船、返回式衛(wèi)星的回收著陸以及火星探測器的著陸等過程中,都得到了廣泛的應(yīng)用。降落傘工作過程是流場結(jié)構(gòu)劇烈作用的高度非線性過程,長期以來大量依賴于工程試驗。隨著數(shù)值算法的不斷進步和計算機硬件水平的提高,采用流場結(jié)構(gòu)耦合分析方法進行降落傘工作過程研究受到國內(nèi)外學者的日益關(guān)注[1-3]。

    早期基于時空離散的降落傘流固耦合方法主要是CFD/MSD模型,將降落傘離散成一系列由彈性力、阻尼力聯(lián)系的質(zhì)點,在流場力的作用下,質(zhì)點發(fā)生運動[4-5]。之后,基于DSD/SST界面跟蹤技術(shù)進行降落傘流固耦合問題的研究成為比較通用的一種方法,在翼傘、十字形傘和環(huán)帆傘等的三維流固耦合計算中得到了應(yīng)用[6-9]。上述方法,流場計算均是采用基于網(wǎng)格的控制體積法或有限元法進行,對于大變形開傘過程,易產(chǎn)生負體積;而對于傘載系統(tǒng)發(fā)生運動時,更需要巨大的計算消耗。

    光滑粒子流體動力學(SPH)方法不需要生成網(wǎng)格,完全消除自由邊界上的數(shù)值發(fā)散,更避免了網(wǎng)格扭曲與網(wǎng)格重構(gòu)等問題,能夠十分方便地模擬網(wǎng)格大變形問題[10],在高速碰撞數(shù)值模擬[11-12]、高能炸藥的爆炸及水下爆炸沖擊等問題[13-15]均得到了很好的應(yīng)用。由于SPH方法粒子的運動與流體或氣體運動相似,SPH粒子方法也應(yīng)用在粘性流體模擬[16]、風沙運動模擬[17]以及液固撞擊數(shù)值模擬[18]方面,文獻[19]中采用SPH方法模擬了跨膜流動,但該方法在降落傘的流固耦合模擬上還鮮有報道。

    本文采用SPH方法成功實現(xiàn)了降落傘穩(wěn)降階段流固耦合數(shù)值模擬,得到了傘衣外形、應(yīng)力分布、速度、氣動力以及擺動角度的變化情況,數(shù)值結(jié)果符合工程實際情況;分析了空氣粒子密度、降落傘速度對穩(wěn)降過程的影響。SPH模擬方法為降落傘流固耦合模擬提供了一個新的研究手段。

    1 數(shù)學模型

    1.1 流場SPH控制方程

    連續(xù)流體的基本守恒方程組:

    SPH離散可以寫成:

    式中分別為粒子和粒子的三維坐標向量;為狄拉克函數(shù);為定義光滑函數(shù)的影響區(qū)域的光滑長度。

    1.2 FEM-SPH耦合算法實現(xiàn)

    對于降落傘,采用有限元方法(FEM)求解,其控制方程為:

    耦合計算通過罰函數(shù)接觸算法實現(xiàn),在流場SPH粒子與有限元結(jié)構(gòu)單元表面之間引入一個界面接觸力,限制SPH粒子對結(jié)構(gòu)單元的穿透,從而模擬固體和流體之間的相互作用。

    界面接觸力(法向F,切向F)的計算公式為:

    2 數(shù)值模型

    本文計算以風洞試驗的傘衣幅數(shù)為8幅的平面圓形傘作為數(shù)值模擬對象[20],載荷質(zhì)量為1.2kg,空氣密度為1.2kg/m3,傘頂孔直徑0.1m,具體計算參數(shù)如表1所示。為了簡化問題,傘衣初始外形假設(shè)為半球形(圖1)。傘衣選用薄殼三角形單元,網(wǎng)格數(shù)為1 674,傘繩采用一維梁單元,網(wǎng)格數(shù)為400。流場計算采用長方體計算區(qū)域,長、寬、高方向尺寸為10m×10m×10m。流場模型如圖2所示。

    表1 降落傘計算參數(shù)

    Tab.1 Calculating parameters of parachute

    圖1 降落傘模型

    圖2 流場模型

    3 數(shù)值結(jié)果及分析

    3.1 粒子密度對數(shù)值結(jié)果的影響

    為了研究粒子密度對數(shù)值結(jié)果的影響,在其他參數(shù)保持不變的情況下,對流場域進行了三種粒子密度的計算,粒子數(shù)目分別為125 000、300 763和512 000,初始速度為40m/s。

    圖3為三種粒子密度下降落傘的穩(wěn)降過程。圖4和圖5分別表示降落傘在不同粒子密度下運動的速度和受力變化情況,穩(wěn)定后的數(shù)值結(jié)果平均值見表2所示。從中可以看出,假設(shè)條件以及計算軟件精度會對結(jié)果產(chǎn)生一定影響,造成計算誤差。此外,粒子數(shù)過稀,計算初期傘繩和流場無法發(fā)生有效的相互作用,傘衣前方受擾流場小,流場與傘衣的相互作用減弱,穩(wěn)態(tài)時氣動力和氣動力系數(shù)均較小,和文獻[20]風洞試驗數(shù)據(jù)(氣動力系數(shù)為0.738)相比,誤差接近19%。當粒子數(shù)達到51萬時,阻力系數(shù)計算誤差高于30萬粒子數(shù)的2.42%,同時計算機時增加了0.6倍。

    圖3 降落傘穩(wěn)降過程

    圖4 降落傘的速度(不同粒子密度)

    圖5 降落傘的受力情況(不同粒子密度)

    表2 穩(wěn)態(tài)數(shù)值結(jié)果

    Tab.2 Numerical results

    可見,粒子數(shù)過稀,傘衣和流場無法發(fā)生有效的流固耦合作用,導致數(shù)值計算結(jié)果誤差較大;但當粒子數(shù)過高時,不僅會大大增加計算消耗,同時也會因為離散點較多,求解計算次數(shù)增加而導致舍入誤差過大影響計算精度,對于本文研究對象,30萬粒子數(shù)較為合理。

    3.2 不同速度穩(wěn)降階段工作特性分析

    在30萬流場粒子數(shù)下,本文模擬了30m/s、40m/s和50m/s初始速度下降落傘的下降情況,計算結(jié)果如圖6~圖9所示。

    圖6 降落傘投影直徑變化

    圖7 降落傘速度

    圖8 降落傘受力情況

    圖9 降落傘軸線擺動角度

    圖6表示三種工況下降落傘投影直徑的變化。降落傘在下降過程中傘衣底邊會收縮,達到最小,然后又慢慢張開達到最大,傘衣的呼吸頻率受速度影響很小,但速度增加,呼吸強度加強。

    圖7~圖8分別為降落傘下降速度及軸向氣動力隨時間的變化曲線。速度越大,初始氣動阻力越大,減速度越快,但達到0.2s時,三者差別不大,其穩(wěn)態(tài)時的平均氣動力和下降速度幾乎趨于一致,和工程實際一致。

    圖9表示下降過程中降落傘對稱軸與下降速度方向之間的夾角。降落傘在下降過程中,傘衣前后左右擺動,初始速度越大,擺角越大。由于運動慣性的影響,擺角的衰減能力遠遠低于氣動力及投影直徑的衰減能力,且擺動周期高于投影直徑的呼吸周期。

    4 結(jié)束語

    本文以文獻[20]風洞試驗用傘為研究對象,采用SPH方法展開了不同粒子數(shù)目、不同初始速度下降落傘穩(wěn)降階段的流固耦合數(shù)值模擬,氣動力系數(shù)和風洞試驗結(jié)果較為一致,得出如下結(jié)論:

    1)粒子密度對計算結(jié)果有一定影響,數(shù)目過少將不能充分地實現(xiàn)結(jié)構(gòu)單元和流場之間的耦合,數(shù)目過多不僅會過大增加計算消耗,也會增加數(shù)值計算誤差。

    2)得到了降落傘投影直徑的變化情況,表明呼吸頻率受速度影響很小,但速度增加,呼吸強度增加。

    3)速度越大,初始氣動阻力越大,減速度越快,但達到平衡時,降落傘的氣動阻力和平衡速度趨于一致。

    4)速度越大,降落傘擺動角度越大,穩(wěn)定性降低。由于慣性影響,擺角的衰減能力遠遠低于氣動力及投影直徑的衰減能力。

    [1] STEIN K R, BENNEY R J, STEEVES E C. A Computational Model that Couples Aerodynamic Structural Dynamic Behavior of Parachutes during the Opening Process[R].ADA264115, NASA, 1993.

    [2] 郭叔偉, 王海濤, 董楊彪, 等. 降落傘呼吸現(xiàn)象研究[J]. 航天返回與遙感, 2010, 31(1): 18-23. GUO Shuwei, WANG Haitao, DONG Yangbiao, et al. Research on Parachute Breath Behavior [J]. Spacecraft Recovery & Remote Sensing, 2010, 31(1): 18-23. (in Chinese)

    [3] 連亮, 王中陽, 張紅英, 等. 基于ALE方法的群傘穩(wěn)降階段的數(shù)值模擬[J]. 航天返回與遙感, 2014, 35(1):21-28. LIAN Liang, WANG Zhongyang, ZHANG Hongying, et al. Numerical Simulation of Cluster Parachute System During Steady-state Descent Phase Based on ALE Method[J]. Spacecraft Recovery & Remote Sensing, 2014, 35(1): 21-28. (in Chinese)

    [4] YU L, MING X. Study on Transient Aerodynamic Characteristics of Parachute Opening Process[J]. Acta Mechanica Sinica, 2007, 23(6): 627-633.

    [5] 余莉, 史獻林, 明曉. 降落傘充氣過程的數(shù)值模擬[J]. 航空學報, 2007, 28(1): 52-57. YU Li, SHI Xianlin, MING Xiao. Numerical Simulation of Parachute during Opening Process[J]. Acta Aeronautica et Astronautica Sinica, 2007, 28(1): 52-57. (in Chinese)

    [6] KALRO V, TEZDUYAR T E. A Parallel 3D Computational Method for Fluid-structure Interactions in Parachute Systems[J]. Computer Methods in Applied Mechanics and Engineering, 2000, 190(3): 321-332.

    [7] STEIN K, BENNEY R, KALRO V, et al. Parachute Fluid-structure Interactions 3-D Computation[J]. Computer Methods in Applied Mechanics and Engineering, 2000, 190(3): 373-386.

    [8] STEIN K, BENNEY R, TEZDUYAR T, et al. Fluid-structure Interaction of a Cross Parachute Numerical Simulation[J]. Computer Methods in Applied Mechanics and Engineering, 2001, 191(6): 673-687.

    [9] TEZDUYAR TE, SATHE S, SCHWAAB M, et al. Fluid-structure Interaction Modeling of Ringsail Parachutes[J]. Computational Mechanics, 2008, 43(1): 133-142.

    [10] FOSTER N, FEDKIW R. Practical Animation of Liquids[C]. Proceedings of International Conference on Computer Graphics and Interactive Techniques. Standford: Stanford University, 2001: 23-30.

    [11] LIBERSKY L D, RANDLES P W, CAMEY TC, et al. Recent Improvements in SPH Modeling of Hypervelocity impact[J]. International Journal of Impact Engineering, 1997, 20(6): 525-532.

    [12] LIBERSKY L D, PETSCHEK A G, CAMEY T C, et al. High Strain Lagrangian Hydrodynamics: A Three-dimensional SPH Code for Dynamic Material Response[J]. Journal of Computational Physics, 1993, 109(1): 67-75.

    [13] LIU M B, LIU G R, ZONG Z, LAM KY. Computer Simulation of High Explosive Explosion Using Smoothed Particle Hydrodynamics Methodology[J]. Computers and Fluids, 2003, 32(3): 305-322.

    [14] LIU M B, LIU G R, ZONG Z, et al. Smoothed Particle Hydrodynamics for Numerical Simulation of Underwater Explosions[J]. Computational Mechanics, 2003, 30(2): 106-118.

    [15] LIU M B, LIU G R, LAM K Y. Investigations into Water Mitigations using a Meshless Particle Method[J]. Shock Waves, 2002, 12(3):181-195.

    [16] 金善勤, 鄭興, 段文洋. 改進 SPH 方法在GPU并行中對粘性流場的模擬[J]. 哈爾濱工程大學學報, 2015, 36(8): 1-7. JIN Shanqin, ZHENG Xing, DUAN Wenyang. Viscosity Flow Simulation by Using Improved SPH Method Based on GPU Parallel Calculation[J]. Journal of Harbin Engineering University, 2015, 36(8): 1-7. (in Chinese)

    [17] 陳福振, 強洪夫, 高巍然. 風沙運動問題的SPH-FVM耦合方法數(shù)值模擬研究[J]. 物理學報, 2014, 63(13): 14-26. CHEN Fuzhen, QIANG Hongfu, GAO Weiran. Smoothed Particle Hydrodynamics Algorithm and Finite Element Method Applied in Numerical Simulation of the Sand Movement Problem[J]. Acta Phys. Sin., 2014, 63(13): 14-26. (in Chinese)

    [18] 王安文, 徐緋, 張岳青. SPH方法在液固撞擊數(shù)值模擬中的應(yīng)用[J]. 計算物理, 2012, 29(4): 525-533. WANG Anwen, XU Fei, ZHANG Yueqing. Smoothed Particle Hydrodynamics Algorithm Applied in Numerical Simulation of Fluid-solid Impact[J]. Chinese Journal of Computational Physics, 2012, 29(4): 525-533. (in Chinese)

    [19] 初陽, 盧文強. 多孔介質(zhì)中流體輸運的孔尺度SPH模擬[J]. 工程熱物理學報, 2009, 30(3): 437-440. CHU Yang, LU Wenqiang. Porescale SPH Simulation of Flow in Porous Media[J]. Journal of Engineering Thermophysics, 2009, 30(3): 437-440. (in Chinese)

    [20] 余莉, 明曉, 陳麗君. 不同透氣情況降落傘的流場試驗研究[J]. 空氣動力學學報, 2008, 26(1):19-25. YU Li, MING Xiao, CHEN Lijun. Experimental Investigation of Parachute in Different Breathability[J]. Acta Aerodynamica Sinica, 2008, 26(1):19-25. (in Chinese)

    Numerical Simulation of Steady-state Descent Phase of Parachute Based on SPH Method

    YAO Xiangru YU Li WU Qiong

    (College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China)

    With the rise of computer technology and numerical simulation, it is possible to describe accurately the dynamic chatacteristics of parachute. In order to eliminate the influence of negative volume due to large mesh deformation in the traditional parachute simulation, this article uses Smoothed Particle Hydrodynamics (SPH) method. The steady-state descent phase of parachute is numerically simulated, and the aerodynamic coefficients are consistent with wind tunnel test results. The influence of air particle density on the numerical results is investigated. It turns out that an effective coupling of the parachute structure with the flow field is difficult to be achieved in low particle density condition, which reduces the calculation accuracy. However, when the particle density exceeds a certain value, the computational consumption will increase and the calculation accuracy will be affected because of the rounding error. Therefore, the optimal particle number is determined. On this basis, the descent processes of parachute with different initial velocities are simulated. The dynamic change of the canopy shape, the projected diameter, the drop speed, the aerodynamic force and the swinging angle are obtained. It is found that the velocity has less impact on the parachute breath rate, the balancing speed as well as the aerodynamic force. The canopy breath amplitude and the swinging angle increase as the velocity increase, and the attenuation property of swinging angle is lower than the project diameter. SPH method can effectively avoid large grid deformation problems in the aspect of parachute fluid-structure interaction and has wide application prospects in the field of large deformation on flexible fabric with fluid-structure interaction method.

    descent process; fluid-structure interaction; smoothed particle hydrodynamics method; parachute

    (編輯:劉穎)

    V244.21

    A

    1009-8518(2016)03-0048-07

    10.3969/j.issn.1009-8518.2016.03.006

    姚向茹,女,1992年生,南京航空航天大學航空宇航學院碩士生。研究方向為飛行器救生及生命保障,回收系統(tǒng)動力學。E-mail:15150678158@163.com。

    2015-11-17

    國家自然科學基金(11172137)資助項目;江蘇高校優(yōu)勢學科建設(shè)工程資助項目(PAPD);南京航空航天大學研究生創(chuàng)新基地(實驗室)開放基金資助(kfjj20150104-01);中央高校基本科研業(yè)務(wù)費專項資金資助

    猜你喜歡
    傘衣氣動力降落傘
    小型降落傘充氣和傘衣塌陷過程
    飛行載荷外部氣動力的二次規(guī)劃等效映射方法
    設(shè)計參數(shù)及大氣參數(shù)對降落傘充氣性能的影響
    降落傘
    側(cè)風對拍動翅氣動力的影響
    傘衣透氣性對翼傘氣動特性的影響
    降落傘
    誰為你折疊降落傘
    NASA"s Mars Landing Experiment Success
    科技傳播(2014年15期)2014-08-30 17:00:55
    高速鐵路接觸線覆冰后氣動力特性的風洞試驗研究
    成年女人在线观看亚洲视频| 免费大片黄手机在线观看| av.在线天堂| 蜜臀久久99精品久久宅男| 久久久久久伊人网av| 国精品久久久久久国模美| 亚洲国产精品成人久久小说| 精品久久久噜噜| 久久久久国产网址| 国产欧美日韩综合在线一区二区| 久久97久久精品| 午夜激情福利司机影院| h视频一区二区三区| 国产精品久久久久久久电影| 97超视频在线观看视频| 美女福利国产在线| 青春草国产在线视频| 欧美一级a爱片免费观看看| 久久久久久久久久人人人人人人| 日本wwww免费看| av国产精品久久久久影院| 麻豆精品久久久久久蜜桃| 99久久精品一区二区三区| 亚洲成人av在线免费| 精品国产乱码久久久久久小说| 99九九线精品视频在线观看视频| 超碰97精品在线观看| 午夜免费男女啪啪视频观看| 国产精品国产三级专区第一集| 免费久久久久久久精品成人欧美视频 | 哪个播放器可以免费观看大片| 久久久精品94久久精品| 最近的中文字幕免费完整| 国产免费又黄又爽又色| 国产精品久久久久久久电影| 精品久久久久久久久av| 老司机影院成人| 午夜福利视频在线观看免费| 亚洲av日韩在线播放| 99热网站在线观看| 精品少妇久久久久久888优播| 日本欧美国产在线视频| 午夜免费观看性视频| 满18在线观看网站| 九色成人免费人妻av| 亚洲四区av| 丝袜喷水一区| 在线天堂最新版资源| 26uuu在线亚洲综合色| 另类精品久久| 免费高清在线观看视频在线观看| 综合色丁香网| 亚洲人成网站在线播| 久久久久国产网址| 精品国产一区二区久久| 热re99久久国产66热| 一级爰片在线观看| 成人毛片60女人毛片免费| 日韩中字成人| 亚洲欧美成人综合另类久久久| 亚洲熟女精品中文字幕| 黄片无遮挡物在线观看| 精品国产一区二区久久| 一级毛片电影观看| 热re99久久国产66热| 亚洲精品国产av成人精品| 嫩草影院入口| 国产精品嫩草影院av在线观看| a 毛片基地| 婷婷成人精品国产| 亚洲精品国产av蜜桃| 精品亚洲成a人片在线观看| 国产精品麻豆人妻色哟哟久久| 国产色婷婷99| 日韩在线高清观看一区二区三区| 久久毛片免费看一区二区三区| 国产免费现黄频在线看| 亚洲怡红院男人天堂| 一二三四中文在线观看免费高清| 精品亚洲成国产av| 三级国产精品欧美在线观看| 亚洲欧洲日产国产| 亚洲人与动物交配视频| 亚洲,欧美,日韩| av播播在线观看一区| 黄色怎么调成土黄色| 久久精品久久久久久久性| 色5月婷婷丁香| 街头女战士在线观看网站| 亚洲国产精品一区三区| 蜜桃久久精品国产亚洲av| 国产精品.久久久| 国产伦精品一区二区三区视频9| 久久综合国产亚洲精品| 精品一区二区三区视频在线| 亚洲av中文av极速乱| 三级国产精品欧美在线观看| 欧美精品一区二区免费开放| 久久人人爽人人爽人人片va| 热re99久久国产66热| 老司机影院成人| 国产爽快片一区二区三区| 国产精品不卡视频一区二区| 伊人久久精品亚洲午夜| 成人亚洲欧美一区二区av| 国产色婷婷99| 啦啦啦中文免费视频观看日本| 欧美精品国产亚洲| 中国三级夫妇交换| 少妇人妻 视频| 精品一区二区免费观看| 亚洲av男天堂| 18禁裸乳无遮挡动漫免费视频| 国产男女内射视频| 人人妻人人添人人爽欧美一区卜| 天美传媒精品一区二区| xxxhd国产人妻xxx| 99久久综合免费| 高清欧美精品videossex| 午夜激情福利司机影院| 亚洲精品一区蜜桃| 国产高清三级在线| 国产成人精品福利久久| 一区二区三区精品91| 午夜激情久久久久久久| 国产精品女同一区二区软件| 亚洲人成网站在线观看播放| 夜夜看夜夜爽夜夜摸| 免费人妻精品一区二区三区视频| 黑人猛操日本美女一级片| 国产成人aa在线观看| 熟妇人妻不卡中文字幕| 国产毛片在线视频| 成人亚洲欧美一区二区av| av免费观看日本| 搡老乐熟女国产| 国产成人aa在线观看| 日韩一区二区视频免费看| 久久亚洲国产成人精品v| 国产男女超爽视频在线观看| 国产永久视频网站| 只有这里有精品99| 99热网站在线观看| 一级毛片电影观看| 一区二区三区免费毛片| 好男人视频免费观看在线| 亚洲四区av| 中文字幕亚洲精品专区| 午夜影院在线不卡| 人妻 亚洲 视频| 波野结衣二区三区在线| 2022亚洲国产成人精品| 亚洲av.av天堂| 天堂8中文在线网| 国模一区二区三区四区视频| 免费黄色在线免费观看| 少妇人妻久久综合中文| 一级片'在线观看视频| 一二三四中文在线观看免费高清| 下体分泌物呈黄色| 久久狼人影院| 国产av码专区亚洲av| 少妇被粗大猛烈的视频| 在线观看美女被高潮喷水网站| 免费看av在线观看网站| 久久久久久久亚洲中文字幕| 国产免费福利视频在线观看| 哪个播放器可以免费观看大片| 欧美少妇被猛烈插入视频| 嫩草影院入口| 三级国产精品片| 男女无遮挡免费网站观看| 97超碰精品成人国产| 国产男女超爽视频在线观看| 两个人的视频大全免费| 亚洲欧美成人综合另类久久久| 男女高潮啪啪啪动态图| 亚洲欧美清纯卡通| 晚上一个人看的免费电影| 国产高清三级在线| 成年人免费黄色播放视频| 一区二区三区精品91| √禁漫天堂资源中文www| 搡老乐熟女国产| 寂寞人妻少妇视频99o| 国产伦精品一区二区三区视频9| 国产精品欧美亚洲77777| 91精品伊人久久大香线蕉| 高清欧美精品videossex| 欧美 日韩 精品 国产| 亚洲欧美色中文字幕在线| 久久精品人人爽人人爽视色| 日本色播在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 免费高清在线观看视频在线观看| 一级a做视频免费观看| 搡老乐熟女国产| 久久精品国产亚洲网站| 伦理电影大哥的女人| 久久久久精品久久久久真实原创| 欧美性感艳星| 国产精品蜜桃在线观看| 亚洲人成77777在线视频| 国产国语露脸激情在线看| 男男h啪啪无遮挡| 精品99又大又爽又粗少妇毛片| 不卡视频在线观看欧美| 日韩人妻高清精品专区| 国产国语露脸激情在线看| 黑人欧美特级aaaaaa片| 国内精品宾馆在线| av在线老鸭窝| 大香蕉久久网| 日本午夜av视频| 亚洲国产精品一区三区| 韩国av在线不卡| 久久久欧美国产精品| 亚洲一区二区三区欧美精品| 九色成人免费人妻av| a 毛片基地| 精品人妻在线不人妻| 91午夜精品亚洲一区二区三区| 日韩av免费高清视频| 下体分泌物呈黄色| 97超碰精品成人国产| 国产 一区精品| 美女内射精品一级片tv| 亚洲精品美女久久av网站| 国产精品人妻久久久影院| 欧美3d第一页| 午夜日本视频在线| 亚洲国产精品成人久久小说| 欧美日韩一区二区视频在线观看视频在线| 久久ye,这里只有精品| 韩国av在线不卡| 国产探花极品一区二区| 日韩中字成人| 王馨瑶露胸无遮挡在线观看| 在线观看人妻少妇| 91午夜精品亚洲一区二区三区| 亚洲精品久久午夜乱码| 亚洲国产av影院在线观看| 午夜91福利影院| 国产色婷婷99| 欧美亚洲 丝袜 人妻 在线| 黄片播放在线免费| 免费黄网站久久成人精品| 大陆偷拍与自拍| 精品久久国产蜜桃| a级毛片黄视频| 久久精品熟女亚洲av麻豆精品| 日本黄大片高清| 97超视频在线观看视频| 亚洲国产最新在线播放| 999精品在线视频| 蜜桃在线观看..| 成人毛片60女人毛片免费| 国产成人免费观看mmmm| 嫩草影院入口| av线在线观看网站| 黄色一级大片看看| 少妇人妻 视频| 免费观看无遮挡的男女| 久久精品国产亚洲av天美| 欧美日韩综合久久久久久| 一级a做视频免费观看| 婷婷成人精品国产| 中文字幕人妻丝袜制服| av卡一久久| 人人妻人人爽人人添夜夜欢视频| 亚洲精品色激情综合| 国产精品久久久久久精品古装| 久久人妻熟女aⅴ| av在线播放精品| 国产乱人偷精品视频| 中文欧美无线码| 黄色视频在线播放观看不卡| 欧美xxⅹ黑人| 我要看黄色一级片免费的| 久久av网站| 3wmmmm亚洲av在线观看| 久久精品人人爽人人爽视色| 久久久亚洲精品成人影院| 精品久久久久久久久av| 成年人午夜在线观看视频| 国产精品一二三区在线看| 国产探花极品一区二区| 国产熟女午夜一区二区三区 | 少妇人妻精品综合一区二区| 另类亚洲欧美激情| 91国产中文字幕| 亚洲av电影在线观看一区二区三区| 青春草国产在线视频| 免费看光身美女| 视频区图区小说| 高清欧美精品videossex| 精品视频人人做人人爽| 亚洲av中文av极速乱| 日日爽夜夜爽网站| 亚洲图色成人| 欧美日韩综合久久久久久| 欧美日韩av久久| 国产av一区二区精品久久| 亚洲综合色惰| 十八禁高潮呻吟视频| 亚洲天堂av无毛| 日韩 亚洲 欧美在线| 国产精品久久久久久久电影| 亚洲综合精品二区| 在线天堂最新版资源| 国产精品久久久久久av不卡| 亚洲经典国产精华液单| 精品少妇内射三级| 成人免费观看视频高清| 日日摸夜夜添夜夜爱| 欧美+日韩+精品| 日韩强制内射视频| 免费看不卡的av| 国产亚洲午夜精品一区二区久久| 精品一区二区三卡| 18禁裸乳无遮挡动漫免费视频| 亚洲美女搞黄在线观看| 亚洲欧美成人综合另类久久久| 久热久热在线精品观看| 国产 一区精品| 青春草国产在线视频| 91午夜精品亚洲一区二区三区| 五月天丁香电影| 日本-黄色视频高清免费观看| videossex国产| 久久久久国产网址| 国产免费视频播放在线视频| 久久久久久久久久久免费av| 国产成人免费无遮挡视频| 天天影视国产精品| 男人操女人黄网站| 欧美三级亚洲精品| 日韩制服骚丝袜av| 国产亚洲精品久久久com| 日韩精品免费视频一区二区三区 | 老女人水多毛片| 天天操日日干夜夜撸| 色视频在线一区二区三区| 丝袜喷水一区| 亚洲精品乱码久久久久久按摩| 精品少妇久久久久久888优播| 少妇精品久久久久久久| 国产黄色免费在线视频| 高清不卡的av网站| 一个人免费看片子| 婷婷色麻豆天堂久久| 欧美日韩国产mv在线观看视频| 亚洲精华国产精华液的使用体验| 极品人妻少妇av视频| 精品人妻熟女av久视频| 国产探花极品一区二区| 久久久久人妻精品一区果冻| 午夜视频国产福利| 成人18禁高潮啪啪吃奶动态图 | 亚洲综合色网址| 免费少妇av软件| 大码成人一级视频| 日韩视频在线欧美| 欧美另类一区| 人妻 亚洲 视频| √禁漫天堂资源中文www| 精品一品国产午夜福利视频| 91久久精品国产一区二区三区| av视频免费观看在线观看| 亚洲内射少妇av| 国产淫语在线视频| 午夜久久久在线观看| 欧美97在线视频| 少妇的逼好多水| 大又大粗又爽又黄少妇毛片口| 国产一区亚洲一区在线观看| 国产色婷婷99| 男女高潮啪啪啪动态图| 91精品一卡2卡3卡4卡| 中文字幕最新亚洲高清| 亚洲在久久综合| 校园人妻丝袜中文字幕| 久久精品人人爽人人爽视色| 丝袜美足系列| av.在线天堂| 国产精品熟女久久久久浪| 日本-黄色视频高清免费观看| 一本一本综合久久| 交换朋友夫妻互换小说| 日韩一区二区三区影片| 69精品国产乱码久久久| 九九爱精品视频在线观看| 国产国语露脸激情在线看| 丝瓜视频免费看黄片| 国产成人91sexporn| 高清在线视频一区二区三区| 在线观看三级黄色| 秋霞伦理黄片| 久久人人爽av亚洲精品天堂| 久久鲁丝午夜福利片| 卡戴珊不雅视频在线播放| 亚洲人成77777在线视频| 国语对白做爰xxxⅹ性视频网站| 99视频精品全部免费 在线| 中文字幕av电影在线播放| 国产伦理片在线播放av一区| 啦啦啦在线观看免费高清www| 精品亚洲成国产av| 亚洲国产精品专区欧美| 欧美精品亚洲一区二区| 中文字幕人妻丝袜制服| 精品少妇久久久久久888优播| 在线观看免费视频网站a站| 亚洲国产av影院在线观看| 999精品在线视频| 久久免费观看电影| 国产免费一区二区三区四区乱码| 国产成人一区二区在线| 男女边吃奶边做爰视频| 国产乱人偷精品视频| 免费观看a级毛片全部| 狂野欧美白嫩少妇大欣赏| 午夜福利在线观看免费完整高清在| 黄色毛片三级朝国网站| 欧美激情 高清一区二区三区| 大片电影免费在线观看免费| 成年av动漫网址| 黄片无遮挡物在线观看| 九色亚洲精品在线播放| 99精国产麻豆久久婷婷| 一本久久精品| 美女大奶头黄色视频| 青春草国产在线视频| 日本午夜av视频| 久久久久久久久久久丰满| av免费在线看不卡| 久久狼人影院| 人人澡人人妻人| 五月开心婷婷网| 性色avwww在线观看| 精品国产露脸久久av麻豆| 国产黄片视频在线免费观看| 精品亚洲乱码少妇综合久久| av线在线观看网站| 男女啪啪激烈高潮av片| 一本久久精品| 国产欧美亚洲国产| 国产男人的电影天堂91| 久久精品夜色国产| 永久免费av网站大全| √禁漫天堂资源中文www| 99久久人妻综合| 男人操女人黄网站| 久久人人爽人人片av| 人体艺术视频欧美日本| 国产黄片视频在线免费观看| 午夜视频国产福利| 少妇的逼好多水| 久久久久久久亚洲中文字幕| 九九久久精品国产亚洲av麻豆| 飞空精品影院首页| 免费人妻精品一区二区三区视频| 丝袜脚勾引网站| 在线看a的网站| 国产在视频线精品| 国产男女超爽视频在线观看| 91久久精品国产一区二区成人| 在线精品无人区一区二区三| 久久午夜综合久久蜜桃| 美女中出高潮动态图| 精品一区二区免费观看| 蜜桃久久精品国产亚洲av| 国产黄色免费在线视频| 51国产日韩欧美| 亚洲欧美一区二区三区黑人 | 久久这里有精品视频免费| 18在线观看网站| 国产精品久久久久久精品古装| 精品国产乱码久久久久久小说| 久久久a久久爽久久v久久| 久久久精品94久久精品| 国产一区二区三区综合在线观看 | 中文字幕人妻熟人妻熟丝袜美| 精品久久国产蜜桃| 黄色配什么色好看| 老司机影院成人| 街头女战士在线观看网站| 久久精品国产鲁丝片午夜精品| 亚洲内射少妇av| 十八禁网站网址无遮挡| 国产亚洲精品久久久com| 汤姆久久久久久久影院中文字幕| 特大巨黑吊av在线直播| 在线观看三级黄色| 男人添女人高潮全过程视频| 一本—道久久a久久精品蜜桃钙片| 纯流量卡能插随身wifi吗| 欧美精品一区二区免费开放| 国产精品国产三级国产专区5o| 2022亚洲国产成人精品| 亚洲精品第二区| 亚洲精品乱码久久久久久按摩| 日韩 亚洲 欧美在线| 精品人妻偷拍中文字幕| 亚洲精品日本国产第一区| 天堂中文最新版在线下载| 亚洲高清免费不卡视频| videosex国产| 啦啦啦视频在线资源免费观看| 九草在线视频观看| 国产成人av激情在线播放 | 亚洲色图 男人天堂 中文字幕 | 一本久久精品| 免费观看在线日韩| 夜夜骑夜夜射夜夜干| 久久久久久久久久人人人人人人| 9色porny在线观看| 一本大道久久a久久精品| 日韩成人伦理影院| av又黄又爽大尺度在线免费看| 成年女人在线观看亚洲视频| 黑人猛操日本美女一级片| 成人午夜精彩视频在线观看| 久久久精品94久久精品| 久久毛片免费看一区二区三区| 中国国产av一级| 大香蕉久久网| 丰满少妇做爰视频| 国产免费一级a男人的天堂| 国产 精品1| 国产极品天堂在线| 久久99热6这里只有精品| 日本av免费视频播放| 国产一级毛片在线| 五月伊人婷婷丁香| 久久精品国产a三级三级三级| 女人精品久久久久毛片| 久久久国产欧美日韩av| 久久人人爽人人爽人人片va| 久热这里只有精品99| 久久午夜综合久久蜜桃| 99国产精品免费福利视频| 亚洲精品久久午夜乱码| 女性生殖器流出的白浆| 亚洲国产最新在线播放| 欧美xxⅹ黑人| 国产成人aa在线观看| 免费人妻精品一区二区三区视频| 亚洲精品456在线播放app| 精品久久久久久久久亚洲| 99久久人妻综合| 各种免费的搞黄视频| av专区在线播放| 欧美性感艳星| 大香蕉久久网| 一本久久精品| 欧美精品人与动牲交sv欧美| 久久精品国产亚洲网站| 日本色播在线视频| 亚洲av在线观看美女高潮| 99久久中文字幕三级久久日本| 久久国产精品男人的天堂亚洲 | 九色成人免费人妻av| 伦精品一区二区三区| 亚洲婷婷狠狠爱综合网| 99国产综合亚洲精品| 老司机影院成人| 亚洲精品久久久久久婷婷小说| 美女视频免费永久观看网站| 欧美日韩精品成人综合77777| 亚洲国产精品国产精品| 久久久久视频综合| 80岁老熟妇乱子伦牲交| 国产精品女同一区二区软件| 国产欧美亚洲国产| 亚洲欧美清纯卡通| 一级,二级,三级黄色视频| 蜜桃国产av成人99| www.色视频.com| 日韩精品免费视频一区二区三区 | 自拍欧美九色日韩亚洲蝌蚪91| 热99久久久久精品小说推荐| 久久99热这里只频精品6学生| 久久精品人人爽人人爽视色| 99久国产av精品国产电影| freevideosex欧美| 亚洲成人一二三区av| av福利片在线| 又黄又爽又刺激的免费视频.| 丰满乱子伦码专区| 在线精品无人区一区二区三| 好男人视频免费观看在线| 天天操日日干夜夜撸| 在线观看www视频免费| 极品人妻少妇av视频| 日韩精品免费视频一区二区三区 | 午夜激情久久久久久久| 色5月婷婷丁香| 午夜福利视频精品| 成人毛片60女人毛片免费| 热re99久久精品国产66热6| 久久久久精品久久久久真实原创| 国产一区有黄有色的免费视频| 久久ye,这里只有精品| 女人久久www免费人成看片| 免费日韩欧美在线观看| 视频在线观看一区二区三区| 国产欧美日韩一区二区三区在线 | 夜夜爽夜夜爽视频| 黄片播放在线免费| 国产成人精品在线电影| 国产精品嫩草影院av在线观看| 免费人妻精品一区二区三区视频| 国产深夜福利视频在线观看| 国产国拍精品亚洲av在线观看|