• <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网站| 国产xxxxx性猛交| 日韩中文字幕欧美一区二区 | 精品一区在线观看国产| 免费在线观看黄色视频的| 交换朋友夫妻互换小说| 18禁国产床啪视频网站| 女的被弄到高潮叫床怎么办| 美女国产高潮福利片在线看| 两个人看的免费小视频| av福利片在线| 波野结衣二区三区在线| 欧美97在线视频| 久久久久久人人人人人| 欧美人与性动交α欧美软件| 看免费av毛片| 亚洲一级一片aⅴ在线观看| 亚洲av欧美aⅴ国产| 亚洲,欧美精品.| 99久久人妻综合| 黄色视频在线播放观看不卡| 99九九在线精品视频| 2018国产大陆天天弄谢| 波多野结衣一区麻豆| 黄频高清免费视频| 亚洲精华国产精华液的使用体验| 国产无遮挡羞羞视频在线观看| 9热在线视频观看99| 国产免费一区二区三区四区乱码| 亚洲精品国产一区二区精华液| 精品国产乱码久久久久久小说| 纵有疾风起免费观看全集完整版| 亚洲av电影在线进入| 国产成人精品福利久久| 亚洲人成网站在线观看播放| 亚洲图色成人| 久久人人爽人人片av| 欧美人与性动交α欧美软件| 一本久久精品| 久久久久久久久免费视频了| 色婷婷av一区二区三区视频| 亚洲精品久久久久久婷婷小说| 91精品三级在线观看| 欧美变态另类bdsm刘玥| 黄网站色视频无遮挡免费观看| 亚洲欧美清纯卡通| 激情五月婷婷亚洲| 十分钟在线观看高清视频www| 久久99蜜桃精品久久| 亚洲中文av在线| 久久久精品免费免费高清| 欧美 日韩 精品 国产| 老司机影院成人| 91在线精品国自产拍蜜月| 极品人妻少妇av视频| 国产成人精品无人区| av有码第一页| 欧美+日韩+精品| 精品一区二区免费观看| 亚洲欧美成人精品一区二区| 国产野战对白在线观看| 女的被弄到高潮叫床怎么办| 美女国产高潮福利片在线看| 丝袜美足系列| 午夜福利影视在线免费观看| 在线观看国产h片| 99久久中文字幕三级久久日本| 青草久久国产| 天天躁夜夜躁狠狠久久av| 久久久久久久久久久免费av| 久久久欧美国产精品| 欧美av亚洲av综合av国产av | 久久久久精品久久久久真实原创| 日韩av不卡免费在线播放| 亚洲人成77777在线视频| 国产成人一区二区在线| 国产成人精品一,二区| 丰满饥渴人妻一区二区三| 亚洲国产精品成人久久小说| 久久精品国产综合久久久| 91精品国产国语对白视频| 色视频在线一区二区三区| 亚洲av电影在线观看一区二区三区| 日韩av免费高清视频| 自拍欧美九色日韩亚洲蝌蚪91| 日本午夜av视频| 日韩人妻精品一区2区三区| 在线天堂最新版资源| 久久久久久久久久久免费av| 午夜福利,免费看| 少妇的逼水好多| 免费黄频网站在线观看国产| www.自偷自拍.com| 在线观看人妻少妇| 国产野战对白在线观看| 久久人人97超碰香蕉20202| 久久人人爽人人片av| 久久精品夜色国产| 女性被躁到高潮视频| 中文乱码字字幕精品一区二区三区| 春色校园在线视频观看| 免费观看性生交大片5| 成人漫画全彩无遮挡| 伦理电影免费视频| 一级片'在线观看视频| 男女边吃奶边做爰视频| 欧美日韩精品成人综合77777| 一区二区三区精品91| 亚洲精品久久久久久婷婷小说| 汤姆久久久久久久影院中文字幕| 亚洲激情五月婷婷啪啪| 亚洲av中文av极速乱| 精品福利永久在线观看| 精品国产乱码久久久久久男人| 精品国产露脸久久av麻豆| 亚洲精品美女久久av网站| 99国产精品免费福利视频| 少妇的逼水好多| 国产白丝娇喘喷水9色精品| 国产xxxxx性猛交| 久久精品aⅴ一区二区三区四区 | 肉色欧美久久久久久久蜜桃| 狠狠婷婷综合久久久久久88av| 久久99一区二区三区| freevideosex欧美| 午夜av观看不卡| 免费在线观看黄色视频的| 一区二区av电影网| 人人妻人人澡人人看| 国产欧美日韩综合在线一区二区| 中文字幕人妻熟女乱码| 性少妇av在线| 成人国产av品久久久| 午夜日韩欧美国产| 最近中文字幕高清免费大全6| 久久97久久精品| av在线播放精品| 国产1区2区3区精品| 免费av中文字幕在线| √禁漫天堂资源中文www| 欧美中文综合在线视频| 多毛熟女@视频| 国产日韩一区二区三区精品不卡| 我的亚洲天堂| 卡戴珊不雅视频在线播放| 久久青草综合色| 综合色丁香网| 久久久久久久精品精品| 男男h啪啪无遮挡| 精品少妇内射三级| 丰满少妇做爰视频| 下体分泌物呈黄色| 天堂俺去俺来也www色官网| 麻豆精品久久久久久蜜桃| 国产精品久久久久久精品电影小说| 国产精品99久久99久久久不卡 | 黄色视频在线播放观看不卡| 免费黄色在线免费观看| 欧美成人精品欧美一级黄| 在线观看人妻少妇| 国产精品久久久av美女十八| 国产精品麻豆人妻色哟哟久久| 熟妇人妻不卡中文字幕| 在线观看免费视频网站a站| kizo精华| 高清在线视频一区二区三区| 国产欧美亚洲国产| 国产精品偷伦视频观看了| 最新中文字幕久久久久| 国产极品粉嫩免费观看在线| 十分钟在线观看高清视频www| 国产激情久久老熟女| 亚洲精品国产色婷婷电影| 国产精品久久久久久av不卡| 美女福利国产在线| 日韩欧美精品免费久久| 老司机影院成人| 人妻少妇偷人精品九色| 日韩一区二区三区影片| 中文字幕亚洲精品专区| www.自偷自拍.com| 国产精品香港三级国产av潘金莲 | 亚洲熟女精品中文字幕| 亚洲欧美成人综合另类久久久| 一本—道久久a久久精品蜜桃钙片| 2018国产大陆天天弄谢| 国产av国产精品国产| 国产免费视频播放在线视频| 我要看黄色一级片免费的| 一级a爱视频在线免费观看| 精品国产超薄肉色丝袜足j| 欧美 日韩 精品 国产| 天天躁狠狠躁夜夜躁狠狠躁| www.自偷自拍.com| 国产av一区二区精品久久| 久久久久久久久久人人人人人人| 国产有黄有色有爽视频| 久久精品亚洲av国产电影网| 91aial.com中文字幕在线观看| 色婷婷av一区二区三区视频| 久久久精品94久久精品| 久久国产亚洲av麻豆专区| 在现免费观看毛片| 一级a爱视频在线免费观看| 韩国高清视频一区二区三区| 丝瓜视频免费看黄片| 乱人伦中国视频| 成人午夜精彩视频在线观看| 久久女婷五月综合色啪小说| 成人毛片60女人毛片免费| 欧美人与善性xxx| 亚洲少妇的诱惑av| 久久久久久人妻| 天天影视国产精品| 18禁国产床啪视频网站| 日本免费在线观看一区| 免费高清在线观看视频在线观看| 久久精品久久久久久噜噜老黄| 性少妇av在线| 国产激情久久老熟女| 黄色毛片三级朝国网站| 男人爽女人下面视频在线观看| 欧美日韩亚洲高清精品| 国产xxxxx性猛交| 亚洲久久久国产精品| 午夜福利网站1000一区二区三区| 国产国语露脸激情在线看| av网站在线播放免费| 观看美女的网站| 好男人视频免费观看在线| 亚洲国产精品成人久久小说| 97在线人人人人妻| kizo精华| 精品国产一区二区三区久久久樱花| 日产精品乱码卡一卡2卡三| 亚洲国产av新网站| 国产高清国产精品国产三级| 午夜日韩欧美国产| 老汉色∧v一级毛片| 亚洲第一av免费看| 欧美激情 高清一区二区三区| 叶爱在线成人免费视频播放| 国产精品无大码| 国产亚洲欧美精品永久| 国产1区2区3区精品| 中文乱码字字幕精品一区二区三区| 午夜精品国产一区二区电影| 丁香六月天网| 99国产精品免费福利视频| 亚洲精品aⅴ在线观看| 国产 一区精品| av网站在线播放免费| 天美传媒精品一区二区| 中文乱码字字幕精品一区二区三区| 纯流量卡能插随身wifi吗| 80岁老熟妇乱子伦牲交| 黄频高清免费视频| 婷婷色综合大香蕉| 国产亚洲欧美精品永久| 国产成人a∨麻豆精品| 97人妻天天添夜夜摸| 熟妇人妻不卡中文字幕| 精品人妻熟女毛片av久久网站| 亚洲精品aⅴ在线观看| 久久这里只有精品19| 中国国产av一级| 9191精品国产免费久久| 免费看av在线观看网站| 中文天堂在线官网| 亚洲国产av影院在线观看| 欧美另类一区| 丝袜美足系列| 国产爽快片一区二区三区| 国产伦理片在线播放av一区| 永久免费av网站大全| 免费观看a级毛片全部| 国产亚洲欧美精品永久| 午夜福利影视在线免费观看| 中文字幕亚洲精品专区| 蜜桃在线观看..| 亚洲欧美精品自产自拍| 女人精品久久久久毛片| av国产久精品久网站免费入址| 亚洲五月色婷婷综合| 天天影视国产精品| 大话2 男鬼变身卡| 日韩中字成人| 在线观看国产h片| 侵犯人妻中文字幕一二三四区| 成年女人在线观看亚洲视频| av网站免费在线观看视频| 日韩av免费高清视频| 水蜜桃什么品种好| 久久久精品免费免费高清| 在线看a的网站| 亚洲国产欧美在线一区| 亚洲精品国产av成人精品| 乱人伦中国视频| 你懂的网址亚洲精品在线观看| 亚洲av.av天堂| 91在线精品国自产拍蜜月| 成人二区视频| 熟女av电影| 免费人妻精品一区二区三区视频| 男女国产视频网站| 亚洲精品日本国产第一区| 涩涩av久久男人的天堂| 大陆偷拍与自拍| 日韩中文字幕视频在线看片| 又黄又粗又硬又大视频| 巨乳人妻的诱惑在线观看| 你懂的网址亚洲精品在线观看| xxx大片免费视频| 欧美精品亚洲一区二区| 国产精品不卡视频一区二区| 永久免费av网站大全| 国产一区亚洲一区在线观看| 久久 成人 亚洲| 亚洲国产精品国产精品| 国产男女超爽视频在线观看| 日韩不卡一区二区三区视频在线| 亚洲,一卡二卡三卡| 中文字幕制服av| 国产白丝娇喘喷水9色精品| 最新中文字幕久久久久| 好男人视频免费观看在线| 丁香六月天网| 精品国产一区二区三区久久久樱花| 七月丁香在线播放| 一本久久精品| 男女午夜视频在线观看| 精品国产乱码久久久久久男人| 亚洲av成人精品一二三区| 午夜福利视频精品| 韩国精品一区二区三区| 女性生殖器流出的白浆| 视频区图区小说| 精品少妇久久久久久888优播| 在线观看www视频免费| 久久久欧美国产精品| 街头女战士在线观看网站| 国精品久久久久久国模美| 国产爽快片一区二区三区| 亚洲经典国产精华液单| 又大又黄又爽视频免费| 咕卡用的链子| 国产色婷婷99| 欧美xxⅹ黑人| 午夜激情久久久久久久| 久久久久精品性色| 国产xxxxx性猛交| 91午夜精品亚洲一区二区三区| 男女边吃奶边做爰视频| 欧美日韩av久久| 久久毛片免费看一区二区三区| 婷婷色麻豆天堂久久| 校园人妻丝袜中文字幕| 999久久久国产精品视频| av有码第一页| 成人国产av品久久久| 黑丝袜美女国产一区| 久久久久国产一级毛片高清牌| 久久久久精品性色| 人人妻人人爽人人添夜夜欢视频| 国产欧美日韩综合在线一区二区| 日韩免费高清中文字幕av| 一区二区日韩欧美中文字幕| 亚洲精品,欧美精品| 久久国产亚洲av麻豆专区| 精品亚洲乱码少妇综合久久| 永久网站在线| 日韩av免费高清视频| 国产精品国产三级国产专区5o| 不卡视频在线观看欧美| 少妇被粗大猛烈的视频| 卡戴珊不雅视频在线播放| 欧美在线黄色| 大香蕉久久成人网| 久久国产亚洲av麻豆专区| 国产乱人偷精品视频| 99国产综合亚洲精品| 爱豆传媒免费全集在线观看| 国产片特级美女逼逼视频| 精品一区二区三卡| 最新的欧美精品一区二区| 嫩草影院入口| 精品一品国产午夜福利视频| 国产极品天堂在线| 成人18禁高潮啪啪吃奶动态图| 婷婷成人精品国产| 黄色一级大片看看| 桃花免费在线播放| 亚洲一区中文字幕在线| 亚洲精品一二三| 欧美日韩精品成人综合77777| 一区二区三区乱码不卡18| 亚洲第一青青草原| 国产一区二区三区av在线| 制服人妻中文乱码| 老熟女久久久| 寂寞人妻少妇视频99o| 国产精品一区二区在线不卡| 香蕉精品网在线| 美女脱内裤让男人舔精品视频| 中文精品一卡2卡3卡4更新| 97精品久久久久久久久久精品| 亚洲经典国产精华液单| 国产精品亚洲av一区麻豆 | 欧美日韩av久久| 老汉色av国产亚洲站长工具| 久久久精品国产亚洲av高清涩受| 久久精品国产综合久久久| 精品一区二区三卡| av在线观看视频网站免费| 观看美女的网站| 老女人水多毛片| 18在线观看网站| 亚洲精品国产av成人精品| 国产精品二区激情视频| 国产精品 欧美亚洲| 黄色 视频免费看| 看免费成人av毛片| 色吧在线观看| 精品少妇一区二区三区视频日本电影 | 精品人妻偷拍中文字幕| 黄片播放在线免费| 精品99又大又爽又粗少妇毛片| 欧美激情高清一区二区三区 | 欧美少妇被猛烈插入视频| 飞空精品影院首页| 国产成人精品久久久久久| 99久久中文字幕三级久久日本| 亚洲精品成人av观看孕妇| 欧美日韩综合久久久久久| 国产高清国产精品国产三级| 纯流量卡能插随身wifi吗| 永久免费av网站大全| 中文字幕av电影在线播放| 久久精品亚洲av国产电影网| 亚洲熟女精品中文字幕| 久久国产亚洲av麻豆专区| 在线观看一区二区三区激情| 午夜日本视频在线| 啦啦啦视频在线资源免费观看| 国产深夜福利视频在线观看| 纵有疾风起免费观看全集完整版| av国产精品久久久久影院| 新久久久久国产一级毛片| 日韩av不卡免费在线播放| 国精品久久久久久国模美| av在线老鸭窝| 1024视频免费在线观看| 国产成人a∨麻豆精品| 精品国产一区二区三区四区第35| 欧美日韩视频精品一区| 国产在视频线精品| 国产精品不卡视频一区二区| 亚洲欧美成人综合另类久久久| 成年人午夜在线观看视频| 99国产精品免费福利视频| 精品久久蜜臀av无| 婷婷成人精品国产| 久久精品国产亚洲av高清一级| 国产日韩一区二区三区精品不卡| 国产日韩欧美亚洲二区| 日产精品乱码卡一卡2卡三| 丝袜美足系列| 少妇的逼水好多| 亚洲伊人久久精品综合| 欧美日韩av久久| 国产成人精品福利久久| 午夜久久久在线观看| 久久久久久人人人人人| 永久免费av网站大全| 日日摸夜夜添夜夜爱| 亚洲第一av免费看| 欧美日韩精品成人综合77777| 精品少妇一区二区三区视频日本电影 | 欧美bdsm另类| 人妻少妇偷人精品九色| 性高湖久久久久久久久免费观看| 国产亚洲午夜精品一区二区久久| 免费黄频网站在线观看国产| 中文欧美无线码| av不卡在线播放| 老汉色av国产亚洲站长工具| 高清av免费在线| 在线亚洲精品国产二区图片欧美| 中文天堂在线官网| 男男h啪啪无遮挡| 熟女av电影| 男女边吃奶边做爰视频| 国产亚洲精品第一综合不卡| 日本色播在线视频| 亚洲国产精品成人久久小说| videosex国产| 色哟哟·www| av在线老鸭窝| 日韩中字成人| 少妇人妻精品综合一区二区| 国产av码专区亚洲av| 日本爱情动作片www.在线观看| 蜜桃在线观看..| 久久久国产一区二区| 新久久久久国产一级毛片| 国产欧美日韩综合在线一区二区| 一级黄片播放器| 制服丝袜香蕉在线| 精品国产国语对白av| 女性生殖器流出的白浆| 亚洲av电影在线观看一区二区三区| 午夜福利视频在线观看免费| 中文字幕人妻熟女乱码| 侵犯人妻中文字幕一二三四区| 国产成人精品久久二区二区91 | 久久精品国产自在天天线| 人妻系列 视频| 国产女主播在线喷水免费视频网站| 热re99久久国产66热| 热99久久久久精品小说推荐| 有码 亚洲区| 日韩电影二区| 各种免费的搞黄视频| 自线自在国产av| 成人黄色视频免费在线看| 亚洲精品第二区| 欧美精品人与动牲交sv欧美| 99久久精品国产国产毛片| 国产精品.久久久| 看免费av毛片| 黄片小视频在线播放| 精品一区在线观看国产| 免费av中文字幕在线| 18禁国产床啪视频网站| 亚洲国产精品一区三区| 精品少妇黑人巨大在线播放| 久久久久久久精品精品| 免费在线观看完整版高清| 在线观看三级黄色| 国产亚洲av片在线观看秒播厂| 欧美日韩一级在线毛片| 免费看av在线观看网站| 国产精品偷伦视频观看了| 久久精品夜色国产| av网站免费在线观看视频| www日本在线高清视频| 一级毛片黄色毛片免费观看视频| 久久 成人 亚洲| 十八禁网站网址无遮挡| 青青草视频在线视频观看| 18禁国产床啪视频网站| 国产片内射在线| 国产av一区二区精品久久| 男女高潮啪啪啪动态图| 午夜免费观看性视频| 免费久久久久久久精品成人欧美视频| 少妇精品久久久久久久| 视频区图区小说| 天天躁夜夜躁狠狠躁躁| 久久精品国产自在天天线| 免费看av在线观看网站| 制服诱惑二区| 麻豆精品久久久久久蜜桃| 超色免费av| 国产日韩欧美视频二区| 永久网站在线| 亚洲一区中文字幕在线| av网站免费在线观看视频| 亚洲欧美一区二区三区久久| 人人妻人人澡人人爽人人夜夜| 精品少妇内射三级| 看十八女毛片水多多多| 99九九在线精品视频| 国产有黄有色有爽视频| 性色avwww在线观看| 欧美激情极品国产一区二区三区| 国产av国产精品国产| 精品酒店卫生间| 精品第一国产精品| 18+在线观看网站| 国产精品欧美亚洲77777| 成人午夜精彩视频在线观看| 免费av中文字幕在线| 午夜福利一区二区在线看| 久久ye,这里只有精品| 久久久欧美国产精品| 亚洲内射少妇av| 蜜桃在线观看..| 人人澡人人妻人| 国产一区二区三区综合在线观看| 91精品国产国语对白视频| 国产精品一区二区在线观看99| 国精品久久久久久国模美| videosex国产| 精品久久久精品久久久| 肉色欧美久久久久久久蜜桃| 欧美 亚洲 国产 日韩一| 免费观看在线日韩| 中文乱码字字幕精品一区二区三区| 精品福利永久在线观看| 人妻系列 视频| 久久精品国产自在天天线| 亚洲图色成人| 日韩三级伦理在线观看| 人妻 亚洲 视频| 中国国产av一级|