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

    水下發(fā)射氣幕提前噴射時間對載荷影響仿真研究

    2017-01-16 02:02:52張曉樂楊興林
    艦船科學(xué)技術(shù) 2016年12期
    關(guān)鍵詞:氣液航行流場

    張曉樂,程 棟,楊興林

    (1. 中國船舶重工集團(tuán)公司 第七一三研究所,河南 鄭州 450015;2. 江蘇科技大學(xué),江蘇 鎮(zhèn)江 212003)

    水下發(fā)射氣幕提前噴射時間對載荷影響仿真研究

    張曉樂1,程 棟1,楊興林2

    (1. 中國船舶重工集團(tuán)公司 第七一三研究所,河南 鄭州 450015;2. 江蘇科技大學(xué),江蘇 鎮(zhèn)江 212003)

    建立結(jié)構(gòu)網(wǎng)格和流場網(wǎng)格耦合的航行體出筒過程氣幕降載計算模型。通過有/無氣幕工況的計算對比,發(fā)現(xiàn)氣幕提前噴射時機(jī)對降載效果影響較大。通過對氣幕提前航行體發(fā)射 0.25 s、0.5 s 和 0.72 s 工況的計算,發(fā)現(xiàn)航行體橫向載荷隨著提前噴射時間延長,呈現(xiàn)先降低后升高的變化規(guī)律,選擇最佳噴射時機(jī)可提高降載效果。該規(guī)律可用于指導(dǎo)氣幕降載工程設(shè)計。

    水下發(fā)射;氣幕;耦合仿真

    0 引 言

    水下航行體垂直出筒過程中,一方面受到發(fā)射筒及適配裝置的約束進(jìn)行軸向運(yùn)動[1],另一方面受到海水橫向流的作用[2],橫向與軸向相對來流使得航行體受到較強(qiáng)非定常流體動力作用[3–5],航行體出筒階段載荷水平較高。為改善航行體出筒受力環(huán)境,將由發(fā)射筒口周圍布置的若干個燃?xì)獍l(fā)生器產(chǎn)生的燃?xì)鈬娚淙牒K校谕部诟浇纬蓺饽?,在航行體外部附著氣體,降低航行體出筒過程受到的海水橫向載荷,改善發(fā)射環(huán)境。

    程棟等[6]結(jié)合理論與噴射試驗總結(jié)出了計算氣幕保護(hù)上下邊界的計算方法,用以指導(dǎo)氣幕發(fā)生裝置設(shè)計。尚書聰?shù)萚7–9]基于動量和動量矩定理建立了導(dǎo)彈水下發(fā)射出筒過程動力學(xué)模型,計算了有無氣幕情況下力學(xué)環(huán)境的改善情況,并對導(dǎo)彈出筒姿態(tài)進(jìn)行數(shù)值仿真。上述研究都是基于理論推導(dǎo)的數(shù)值計算研究,未考慮高能燃?xì)庠谒械拿}動對航行體的影響。

    本文采取剛體結(jié)構(gòu)網(wǎng)格和流場網(wǎng)格耦合的方法[10–11],結(jié)合一個算例開展航行體出筒過程氣幕降載效果的三維數(shù)值仿真計算。計算有/無氣幕和多個氣幕噴射時間提前航行體發(fā)射的工況,研究氣幕噴射時機(jī)對航行體橫向載荷的影響。

    1 計算模型

    1.1 流場模型

    由于算例中流體域相對于過軸線的來流速度方向?qū)ΨQ,為減小計算量建立對稱計算模型。由于筒內(nèi)燃?xì)鈱Τ鐾策^程沒有影響,流場計算模型只建立了海水域,不考慮發(fā)射筒及筒內(nèi)流場,海水域流場區(qū)域 9 m × 4 m × 16 m,見圖 1。為保證計算精度,在航行體運(yùn)行軌跡及氣幕噴射區(qū)域進(jìn)行 Euler 網(wǎng)格加密,網(wǎng)格密度小于 1.7% 航行體直徑。

    圖 1 流場計算網(wǎng)格Fig. 1 Mesh of fluid field

    計算區(qū)域及邊界條件如圖 2 所示。流場區(qū)域上下為壓力邊界,左右為海水速度入口和出口邊界,海水橫向流動速度為 X 正向,速度為航行體出筒速度的6%。最底部為發(fā)射平臺,設(shè)置為 wall 邊界。對稱面設(shè)置為對稱邊界條件。

    1.2 航行體及噴口模型

    航行體沿軸線建立一半對稱模型,如圖 3 所示。航行體采用 Lagrange 網(wǎng)格離散,并設(shè)置為剛體,除了垂直向上的運(yùn)動外,其余 5 個自由度固定。航行體筒內(nèi)勻加速運(yùn)動,航行體出筒速度為 VO。航行體圓柱段直徑為 Dm,總長為 7 Dm。

    在航行體迎流側(cè),以航行體軸線為圓心布置 5 個氣幕噴管,如圖 4 所示。噴口同樣采用 Lagrange 網(wǎng)格離散且 6 個自由度固定。噴口頂部網(wǎng)格設(shè)置為氣幕燃?xì)獾鸟詈厦娉隹?。氣幕燃?xì)庥晒腆w發(fā)動機(jī)產(chǎn)生,速度入口穩(wěn)定段流速為 2.3 VO,穩(wěn)定段壓力 1.9 MPa。

    圖 2 計算域及邊界條件Fig. 2 An outline of the computational domain with boundary conditions

    圖 3 航行體模型Fig. 3 Model of vehicle

    圖 4 氣幕噴管位置Fig. 4 Location of gas screen nozzle

    圖 5 耦合計算固體與流體重疊網(wǎng)格Fig. 5 Overlap mesh of solid and fluid field in coupling simulation

    算例初始時刻航行體和氣幕噴口的 Lagrange 計算網(wǎng)格與流場區(qū)域計算 Euler 網(wǎng)格部分重疊(見圖 5)。固體網(wǎng)格表面為流體-固體耦合面,計算中耦合面內(nèi)重疊區(qū)域 Euler 網(wǎng)格材料自動為空。圖 6 為采用本文計算方法與某航行體水下出筒試驗中,航行體背流面某測點壓力結(jié)果對比。計算模型不考慮航行體筒內(nèi)運(yùn)行過程的壓力,故初始段計算壓力為 0。由圖 6 可知,出筒后計算壓力曲線與試驗曲線吻合較好,說明采用的計算方法合理可信。

    3 計算結(jié)果分析

    3.1 有無氣幕計算結(jié)果對比

    航行體按照設(shè)定軌跡出筒,整個時間設(shè)置為 0.8 s。圖 7 為無氣幕防護(hù)狀態(tài)航行體出筒過程圖(出筒 0.46 s時狀態(tài),為便于對比以下各圖同)。由于出筒速度較高,航行體肩部會產(chǎn)生少量空泡。圖 8 為計算的無氣幕狀態(tài)下航行體出筒過程受力曲線。隨著航行體入水長度的增加,航行體的 X 方向受力在不斷變化。為對比各種工況氣幕對航行體的載荷,取航行體出筒過程X 方向平均受力作為比較依據(jù)。根據(jù)計算結(jié)果無氣幕狀態(tài),航行體出筒過程 X 向平均受力 30 140 N。

    圖 9 為氣幕-發(fā)射同步啟動工況航行體出筒過程氣液界面圖。航行體發(fā)射筒內(nèi)運(yùn)動的零時刻,氣幕同步開始噴射。由圖可以看出,產(chǎn)生的氣幕主要分布在航行體 X 向迎流面,部分將航行體與來流海水隔離。圖 10為計算得到的氣幕-發(fā)射同步啟動工況下航行體出筒過程受力曲線。根據(jù)計算,氣幕-發(fā)射同步啟動航行體出筒過程 X 向平均受力 52 374 N,高于無氣幕狀態(tài)航行體出筒過程平均受力 73.8%。說明此工況下氣幕噴射沒有起到降低橫向載荷的效果。

    圖 6 計算與試驗航行體背流面測點壓力對比Fig. 6 Pressure comparison between simulation and test point which located in the back of vehicle

    圖 7 無氣幕航行體出筒過程速度矢量圖Fig. 7 Velocity vector of the vehicle launch with no gas screen

    圖 8 無氣幕航行體出筒過程 X 方向受力曲線Fig. 8 Loads curve in X direction of the vehicle launch with no gas screen

    圖 9 氣幕噴射-發(fā)射同步啟動航行體出筒氣液界面圖Fig. 9 Fluid phases of vehicle launch with gas screen jet simultaneity

    圖 10 氣幕噴射-發(fā)射同步啟動航行體 X 方向受力曲線Fig. 10 Loads curve in X direction of the vehicle launch with gas screen jet simultaneity

    氣幕噴口需要將高溫燃?xì)庖暂^高的射流速度注入海水。在形成穩(wěn)定的氣幕前,高溫燃?xì)膺M(jìn)入海水膨脹,排開海水做功。氣幕噴口位于航行體的迎流面,當(dāng)氣幕噴射與航行體發(fā)射同步啟動,氣幕膨脹做功推動海水的同時也給航行體 X 正向作用力,疊加海水 X方向作用力后,加劇了航行體的橫向載荷。為了解決此問題,采取提前航行體發(fā)射時間零點進(jìn)行噴射氣幕,使高溫氣幕在海水中充分釋放能量,形成較為穩(wěn)定的氣幕形態(tài),以實現(xiàn)航行體出筒過程降載。

    3.2 氣幕提前噴射計算

    為研究氣幕提前時間對降載效果的影響,分別進(jìn)行了氣幕提前 0.25 s、提前 0.5 s 和提前 0.72 s 工況的仿真計算。圖 11 給出了氣幕提前 0.25 s 噴射航行體出筒過程氣液界面圖。相對于無提前噴射工況,氣幕體積明顯增大。從圖 12 氣幕提前 0.25 s 噴射航行體 X 方向受力曲線可以計算出,該工況下行形體出筒過程 X 向平均受力 22 294 N,相對于無氣幕工況和氣幕-發(fā)射同步啟動工況均有大幅減小。

    圖 13 給出了氣幕提前 0.5 s 噴射航行體出筒過程氣液界面圖。由氣液界面分布可知,此工況下航行體在出筒過程中已經(jīng)部分穿過氣幕。從圖 14 氣幕提前0.5 s 噴射航行體 X 方向受力曲線可以計算出,該工況下航行形體出筒過程 X 向平均受力 27 285 N,相對于氣幕提前 0.25 s 噴射工況有所增加。

    圖 15 給出了氣幕提前 0.72 s 噴射航行體出筒過程氣液界面圖。由氣液界面分布可知,在航行體接觸氣幕時,大部分提前噴射的氣幕已經(jīng)運(yùn)動到航行體橫向背流面。從圖 16 氣幕提前 0.72 s 噴射航行體 X 方向受力曲線可以計算出,該工況下行形體出筒過程 X 向平均受力 32 388 N,此值已經(jīng)和無氣幕狀態(tài)航行體橫向受力狀態(tài)接近,并略有增加。

    圖 11 氣幕提前 0.25 s 噴射航行體出筒氣液界面圖Fig. 11 Fluid phases of gas screen jet 0.25 s before vehicle launch

    圖 12 氣幕提前 0.25 s 噴射航行體 X 方向受力曲Fig. 12 Vehicle loads curve in X direction of gas screen jet0.25 s before launch

    圖 13 氣幕提前 0.5 s 噴射航行體出筒氣液界面圖Fig. 13 Fluid phases of gas screen jet 0.5 s before vehicle launch

    圖 14 氣幕提前 0.5 s 噴射航行體 X 方向受力曲線Fig. 14 Vehicle loads curve in X direction of gas screen jet 0.5 s before launch

    圖 15 氣幕提前 0.72 s 噴射航行體出筒氣液界面圖Fig. 15 Fluid phases of gas screen jet 0.72 s before vehicle launch

    圖 16 氣幕提前 0.72 s 噴射航行體 X 方向受力曲線Fig. 16 Vehicle loads curve in X direction of gas screen jet 0.72 s before launch

    3.3 仿真結(jié)果統(tǒng)計

    表 1 給出了本文仿真計算各個工況結(jié)果的統(tǒng)計,并將統(tǒng)計結(jié)果繪制在圖 17 上。由圖可知,氣幕能否達(dá)到降載效果與噴射時機(jī)有很大關(guān)系,隨著提前噴射時間延長,航行體橫向載荷呈現(xiàn)先降低后升高的變化趨勢。氣幕噴射-發(fā)射同步啟動工況會加劇航行體的受力,氣幕噴射時間提前發(fā)射零點過早也可能達(dá)不到降載效果。因此,為實現(xiàn)最佳氣幕降載效果,需要選擇恰當(dāng)?shù)膰娚鋾r機(jī)。既要使氣幕燃?xì)饽芰吭诤K谐浞轴尫?,形成穩(wěn)定的氣幕,又要避免過早噴射導(dǎo)致氣幕提前越過了航行體。對于本文的算例來說,提前 0.25 s噴射位于最佳噴射時機(jī)附近,選擇此提前時間可以有效降低航行體出筒載荷。

    表 1 各工況降載效果統(tǒng)計Tab. 1 Statistics of loads in simulation result

    圖 17 各工況提前噴射時間與降載效果曲線Fig. 17 Curve of vehicle loads with different jet time before launch

    4 結(jié) 語

    本文采用 Lagrange 結(jié)構(gòu)網(wǎng)格和 Euler 流場網(wǎng)格耦合的方法,對潛射航行器出筒過程氣幕降載效果進(jìn)行了仿真計算,并通過多工況對比分析了降載效果隨提前噴射時間的變化關(guān)系。得到如下結(jié)論:

    1)采用高能燃?xì)庾鳛闅饽粚叫畜w進(jìn)行降載,從開始噴射到形成穩(wěn)定氣幕需要一個過程,在此過程可能會加劇航行體載荷。

    2)氣幕能否達(dá)到降載效果與氣幕噴射提前航行體發(fā)射零點的時間有很大關(guān)系,航行體橫向載荷隨著提前噴射時間延長,呈現(xiàn)先降低后升高的變化趨勢。

    3)選擇合適的噴射時機(jī)(如本算例提前 0.25 s 附近)可以有效降低航行體出筒橫向載荷。

    [1]倪火才. 潛地彈道導(dǎo)彈發(fā)射裝置構(gòu)造[M]. 哈爾濱: 哈爾濱工程大學(xué)出版社, 1998: 40–50.

    [2]趙世平, 蔡體敏. 橫向流對潛艇垂直發(fā)射導(dǎo)彈的影響[J]. 船舶力學(xué), 2006, 10(4): 33–37. ZHAO Shi-ping, CAI Ti-min. Effects of lateral flow to the missile vertical launched from a submarine[J]. Journal of Ship Mechanics, 2006, 10(4): 33–37.

    [3]BURGDORF O. Hydrodynamics of unsteady underwaterlaunched missiles with trailing cavities and crossflow drag[C]//Proceedings of the 19th AIAA, Fluid Dynamics, Plasma Dynamics, and Lasers Conference. Honolulu, HI, USA: AIAA, 1987: 3–8.

    [4]尚書聰, 張宇文, 袁緒龍. 導(dǎo)彈水下動機(jī)座垂直發(fā)射過程的水動力計算[J]. 艦船科學(xué)技術(shù), 2009, 31(9): 56–60. SHANG Shu-cong, ZHANG Yu-wen, YUAN Xu-long. The hydrodynamic computation on moving base vertical launching of underwater missile[J]. Ship Science and Technology, 2009, 31(9): 56–60.

    [5]黃建春, 葉取源, 朱世權(quán). 不同發(fā)射深度下導(dǎo)彈水下點火氣水流體動力計算[J]. 應(yīng)用力學(xué)學(xué)報, 1994, 11(3): 19–24. HUANG Jian-chun, YE Qu-yuan, ZHU Shi-quan. Gas-water dynamic calculation for the underwater ignition of a missile at different depths[J]. Chinese Journal of Applied Mechanics, 1994, 11(3): 19–24.

    [6]程棟, 何國強(qiáng), 邢軍, 等. 水下發(fā)射筒口氣幕上下邊界計算方法研究[J]. 彈箭與制導(dǎo)學(xué)報, 2010, 30(6): 165–167. CHENG Dong, HE Guo-qiang, XING Jun, et al. Underwater launch tube outlet on the upper and lower gas screen border calculation method[J]. Journal of Projectiles, Rockets, Missiles and Guidance, 2010, 30(6): 165–167.

    [7]尚書聰, 孫建中, 秦麗萍. 潛載導(dǎo)彈水下發(fā)射出筒橫向動力學(xué)特性研究[J]. 振動與沖擊, 2012, 31(23): 82–86.

    [8]尚書聰, 孫建中, 程棟, 等. 筒口氣幕環(huán)境的導(dǎo)彈出筒過程受力影響[J]. 哈爾濱工程大學(xué)學(xué)報, 2012, 33(11): 1423–1427, 1434.

    [9]尚書聰, 孫建中, 程棟. 筒口氣幕環(huán)境下的導(dǎo)彈出筒姿態(tài)數(shù)值仿真[J]. 應(yīng)用力學(xué)學(xué)報, 2013, 30(3): 428–433. SHANG Shu-cong, SUN Jian-zhong, CHENG Dong. The numerical simulation for the outlet attitude of the missile based on the gas screen environment[J]. Chinese Journal of Applied Mechanics, 2013, 30(3): 428–433.

    [10]蔡斯淵, 侯海量, 吳林杰. 隔層設(shè)置對防雷艙液艙防護(hù)能力的影響[J]. 哈爾濱工程大學(xué)學(xué)報, 2016, 37(4): 527–532.

    [11]吳國民, 周心桃, 李俊. 水下爆炸作用下固支平板動態(tài)響應(yīng)分析[J]. 艦船科學(xué)技術(shù), 2013, 35(4): 25–28, 39. WU Guo-min, ZHOU Xin-tao, LI Jun. Dynamical response of flat plates subjected to underwater explosion[J]. Ship Science and Technology, 2013, 35(4): 25–28, 39.

    [12]司海龍, 虞昊, 李政杰, 等. 船體尾壓浪板砰擊載荷分析[J].艦船科學(xué)技術(shù), 2015, 37(2): 19–23, 29. SI Hai-long, YU Hao, LI Zheng-jie, et al. The research of slamming pressure impacted on wave suppression board[J]. Ship Science and Technology, 2015, 37(2): 19–23, 29.

    The numerical simulation of gas screen jet time ahead of underwater-launch which affect vehicle loads

    ZHANG Xiao-le1, CHENG Dong1, YANG Xing-lin2
    (1. The 713 Research Institute of CSIC, Zhengzhou 450015, China; 2. Jiangsu University of Science and Technology, Zhenjiang 212003, China)

    The gas screen jet in the underwater-launch was simulated , using a 3D symmetric model based on the coupling of Lagrange structure mesh and Euler fluid mesh. It was found the jet time ahead of underwater-launch is important to the vehicle horizontal loads, through the simulations with /or without the gas screen. Ahead time of 0.25 s、0.5 s and 0.72 s for gas screen jet were simulated, and vehicle horizontal loads change tendency with growing after reducing was found. It shows that the vehicle loads can be reduced by choosing the proper time. The result can be used in the design of gas screen.

    underwater-launch;gas screen;coupling simulation

    TJ399

    A

    1672–7619(2016)12–0178–05

    10.3404/j.issn.1672–7619.2016.12.037

    2016–10–10

    “十二五”預(yù)研資助項目(1010403050102)

    張曉樂(1981–),男,博士,高級工程師,主要研究方向為水下發(fā)射技術(shù)。

    猜你喜歡
    氣液航行流場
    到慧骃國的航行
    大型空冷汽輪發(fā)電機(jī)轉(zhuǎn)子三維流場計算
    微重力下兩相控溫型儲液器內(nèi)氣液界面仿真分析
    氣液分離罐液位計接管泄漏分析
    轉(zhuǎn)杯紡排雜區(qū)流場與排雜性能
    小舟在河上航行
    基于HYCOM的斯里蘭卡南部海域溫、鹽、流場統(tǒng)計分析
    航行
    青年歌聲(2017年6期)2017-03-13 00:57:56
    CO2 驅(qū)低液量高氣液比井下氣錨模擬與優(yōu)化
    基于瞬態(tài)流場計算的滑動軸承靜平衡位置求解
    国产av麻豆久久久久久久| 天美传媒精品一区二区| 亚洲片人在线观看| 啦啦啦免费观看视频1| 精品电影一区二区在线| 亚洲中文字幕一区二区三区有码在线看| 欧美乱码精品一区二区三区| 国产伦在线观看视频一区| 日本一本二区三区精品| 丁香六月欧美| 在线播放无遮挡| 亚洲最大成人中文| 在线观看免费视频日本深夜| 在线天堂最新版资源| 每晚都被弄得嗷嗷叫到高潮| 午夜免费成人在线视频| 啦啦啦免费观看视频1| www.999成人在线观看| 亚洲精品粉嫩美女一区| 亚洲欧美一区二区三区黑人| 嫁个100分男人电影在线观看| 日本免费一区二区三区高清不卡| 亚洲美女黄片视频| 日韩欧美国产一区二区入口| 91av网一区二区| 免费av毛片视频| 国产精品日韩av在线免费观看| 51午夜福利影视在线观看| 禁无遮挡网站| 最新在线观看一区二区三区| 高清在线国产一区| 午夜日韩欧美国产| 色在线成人网| 亚洲av第一区精品v没综合| 很黄的视频免费| 国产av麻豆久久久久久久| 国产免费一级a男人的天堂| 岛国视频午夜一区免费看| 国模一区二区三区四区视频| 亚洲人成网站高清观看| 日本免费一区二区三区高清不卡| 天天添夜夜摸| 女人被狂操c到高潮| 亚洲午夜理论影院| 精品久久久久久,| 免费看美女性在线毛片视频| 欧美中文日本在线观看视频| 日韩精品中文字幕看吧| 国产欧美日韩一区二区三| 色精品久久人妻99蜜桃| 黑人欧美特级aaaaaa片| 久久久久久九九精品二区国产| 久久精品影院6| 久9热在线精品视频| 人妻丰满熟妇av一区二区三区| 亚洲精品在线美女| 欧美一区二区国产精品久久精品| 国产69精品久久久久777片| 韩国av一区二区三区四区| 婷婷精品国产亚洲av在线| 亚洲无线在线观看| 国产成人欧美在线观看| 亚洲色图av天堂| 国产主播在线观看一区二区| 久久午夜亚洲精品久久| 人人妻人人看人人澡| 亚洲精品亚洲一区二区| 搡老妇女老女人老熟妇| 黄色丝袜av网址大全| av天堂中文字幕网| 天堂√8在线中文| 一级作爱视频免费观看| 亚洲无线观看免费| av欧美777| 天天躁日日操中文字幕| 亚洲av成人精品一区久久| 亚洲 欧美 日韩 在线 免费| 男女床上黄色一级片免费看| 国内精品久久久久久久电影| 国产毛片a区久久久久| 国产精品香港三级国产av潘金莲| 亚洲aⅴ乱码一区二区在线播放| 人妻久久中文字幕网| 婷婷亚洲欧美| 五月玫瑰六月丁香| 麻豆国产av国片精品| 免费一级毛片在线播放高清视频| 在线视频色国产色| 免费人成视频x8x8入口观看| 国产精品久久电影中文字幕| 国产精品98久久久久久宅男小说| 他把我摸到了高潮在线观看| 桃红色精品国产亚洲av| 午夜免费男女啪啪视频观看 | 天堂√8在线中文| 可以在线观看的亚洲视频| 看片在线看免费视频| 91字幕亚洲| 国产精品自产拍在线观看55亚洲| aaaaa片日本免费| 国产精品亚洲美女久久久| 美女免费视频网站| 99久久无色码亚洲精品果冻| 久久人妻av系列| 欧美日韩黄片免| 国产欧美日韩精品亚洲av| 狂野欧美白嫩少妇大欣赏| netflix在线观看网站| 亚洲欧美日韩无卡精品| 国产成人影院久久av| 99久久综合精品五月天人人| 国产成人啪精品午夜网站| 俄罗斯特黄特色一大片| 欧美成人a在线观看| 亚洲成av人片免费观看| 老司机福利观看| 99精品在免费线老司机午夜| 婷婷亚洲欧美| 欧美日本亚洲视频在线播放| www.999成人在线观看| 亚洲精华国产精华精| 久久久国产成人精品二区| 啦啦啦韩国在线观看视频| 网址你懂的国产日韩在线| 国产av在哪里看| 国产一区二区在线av高清观看| 十八禁人妻一区二区| 午夜精品一区二区三区免费看| 久久精品人妻少妇| 国产精品 国内视频| 观看免费一级毛片| 亚洲乱码一区二区免费版| 丝袜美腿在线中文| 国产激情欧美一区二区| 久久人人精品亚洲av| 亚洲av熟女| 国产在视频线在精品| 日本黄大片高清| 欧美日韩中文字幕国产精品一区二区三区| 校园春色视频在线观看| 亚洲成人免费电影在线观看| 有码 亚洲区| 两人在一起打扑克的视频| 国内精品一区二区在线观看| 亚洲av成人av| 中出人妻视频一区二区| 亚洲 欧美 日韩 在线 免费| 一夜夜www| av黄色大香蕉| 色在线成人网| 1024手机看黄色片| 国产成+人综合+亚洲专区| 国产不卡一卡二| 国内精品美女久久久久久| 日本精品一区二区三区蜜桃| 国产精品久久久久久久久免 | 一个人看视频在线观看www免费 | 亚洲一区二区三区不卡视频| 久久精品夜夜夜夜夜久久蜜豆| 亚洲五月婷婷丁香| 午夜免费男女啪啪视频观看 | 免费观看精品视频网站| 亚洲专区中文字幕在线| 18禁美女被吸乳视频| 天堂影院成人在线观看| 日本成人三级电影网站| 少妇的逼好多水| 97超级碰碰碰精品色视频在线观看| 老熟妇乱子伦视频在线观看| or卡值多少钱| 淫妇啪啪啪对白视频| av天堂在线播放| 日韩av在线大香蕉| 成人精品一区二区免费| 高清日韩中文字幕在线| 最新在线观看一区二区三区| 精品一区二区三区视频在线观看免费| 最近最新中文字幕大全免费视频| 丁香六月欧美| 男女下面进入的视频免费午夜| 久久久色成人| 国产毛片a区久久久久| 国产精品一区二区三区四区久久| 亚洲成a人片在线一区二区| 婷婷精品国产亚洲av| 国产一区二区在线av高清观看| www国产在线视频色| 欧美三级亚洲精品| 91久久精品国产一区二区成人 | 婷婷六月久久综合丁香| 国产精品嫩草影院av在线观看 | 蜜桃亚洲精品一区二区三区| 国产高清视频在线播放一区| 午夜老司机福利剧场| 小蜜桃在线观看免费完整版高清| 香蕉av资源在线| 午夜激情福利司机影院| 国产精品亚洲一级av第二区| 两个人视频免费观看高清| 怎么达到女性高潮| 他把我摸到了高潮在线观看| 欧美一级毛片孕妇| 男人舔奶头视频| 一区二区三区免费毛片| 久久精品国产亚洲av涩爱 | 国产精品嫩草影院av在线观看 | 法律面前人人平等表现在哪些方面| 每晚都被弄得嗷嗷叫到高潮| 久久久久国内视频| 好看av亚洲va欧美ⅴa在| 国产精品久久久久久精品电影| 99在线视频只有这里精品首页| 成人高潮视频无遮挡免费网站| 深夜精品福利| 国内久久婷婷六月综合欲色啪| 精品电影一区二区在线| 久久性视频一级片| 757午夜福利合集在线观看| 十八禁网站免费在线| 成人鲁丝片一二三区免费| 精品一区二区三区视频在线观看免费| 日本 欧美在线| 午夜日韩欧美国产| 少妇的逼好多水| 亚洲久久久久久中文字幕| 亚洲 国产 在线| 免费看十八禁软件| 在线十欧美十亚洲十日本专区| www.999成人在线观看| 色av中文字幕| 人人妻人人澡欧美一区二区| 桃红色精品国产亚洲av| 黄色日韩在线| 一个人免费在线观看的高清视频| 中文字幕久久专区| 18禁国产床啪视频网站| 操出白浆在线播放| 91在线观看av| 2021天堂中文幕一二区在线观| 国产av不卡久久| 桃红色精品国产亚洲av| 免费在线观看日本一区| 免费看a级黄色片| 午夜福利视频1000在线观看| 国产乱人伦免费视频| 国产精品日韩av在线免费观看| 成人18禁在线播放| 亚洲自拍偷在线| 国产欧美日韩一区二区精品| 久久香蕉精品热| 在线观看美女被高潮喷水网站 | 两性午夜刺激爽爽歪歪视频在线观看| 九九热线精品视视频播放| 精品99又大又爽又粗少妇毛片 | 精品一区二区三区人妻视频| 俺也久久电影网| 日韩国内少妇激情av| 国产国拍精品亚洲av在线观看 | 色在线成人网| 很黄的视频免费| 久久精品91蜜桃| 国产成年人精品一区二区| 免费观看的影片在线观看| 国产av麻豆久久久久久久| 人人妻人人澡欧美一区二区| 51午夜福利影视在线观看| 国产精品永久免费网站| 亚洲国产高清在线一区二区三| 有码 亚洲区| 麻豆国产av国片精品| 国产探花在线观看一区二区| 成人国产一区最新在线观看| 欧美性感艳星| 国产精品久久电影中文字幕| 午夜两性在线视频| 国产一区二区在线av高清观看| 免费在线观看成人毛片| 国产午夜精品论理片| 一a级毛片在线观看| 国产成人影院久久av| 精品国内亚洲2022精品成人| 国产精品影院久久| 国产美女午夜福利| 欧美日韩国产亚洲二区| 在线观看一区二区三区| 国产精华一区二区三区| 亚洲人成网站在线播放欧美日韩| 一个人观看的视频www高清免费观看| 久久精品国产综合久久久| 男女床上黄色一级片免费看| 国产在线精品亚洲第一网站| 99热只有精品国产| 久久久久精品国产欧美久久久| 午夜福利18| av女优亚洲男人天堂| 国产老妇女一区| 欧美zozozo另类| 国产亚洲精品一区二区www| 国产av不卡久久| 国产又黄又爽又无遮挡在线| 国产精品久久久人人做人人爽| 久久精品国产99精品国产亚洲性色| 午夜福利高清视频| 美女cb高潮喷水在线观看| 国产黄片美女视频| 亚洲av熟女| 欧美+日韩+精品| 久久久久久大精品| 精品久久久久久久人妻蜜臀av| 深夜精品福利| 久久久精品大字幕| 五月玫瑰六月丁香| 亚洲欧美日韩高清在线视频| 国产欧美日韩精品亚洲av| 久久人妻av系列| 国产精品久久久人人做人人爽| 午夜福利成人在线免费观看| 日韩 欧美 亚洲 中文字幕| 首页视频小说图片口味搜索| 国产真实伦视频高清在线观看 | 一本久久中文字幕| 熟女电影av网| 日韩高清综合在线| 国产老妇女一区| 成人午夜高清在线视频| 窝窝影院91人妻| 久久婷婷人人爽人人干人人爱| 51国产日韩欧美| 少妇高潮的动态图| 少妇的逼好多水| 国产精品久久久久久久久免 | 小蜜桃在线观看免费完整版高清| 国产精品久久久久久亚洲av鲁大| 成人午夜高清在线视频| 亚洲在线观看片| 久久人人精品亚洲av| 99riav亚洲国产免费| 中文字幕高清在线视频| 亚洲精品成人久久久久久| 俺也久久电影网| 香蕉久久夜色| 亚洲av免费在线观看| 国产单亲对白刺激| 午夜两性在线视频| 乱人视频在线观看| 精品人妻一区二区三区麻豆 | 91字幕亚洲| 毛片女人毛片| 日韩精品中文字幕看吧| 99在线人妻在线中文字幕| 国产美女午夜福利| 成人鲁丝片一二三区免费| 欧美区成人在线视频| 天天添夜夜摸| 看免费av毛片| 亚洲片人在线观看| 欧美黄色淫秽网站| 欧美+亚洲+日韩+国产| 99久久无色码亚洲精品果冻| 亚洲第一电影网av| 欧美又色又爽又黄视频| 俺也久久电影网| 岛国在线观看网站| 少妇的逼好多水| 欧美一级毛片孕妇| 一区二区三区高清视频在线| 国内毛片毛片毛片毛片毛片| 高清在线国产一区| 免费看a级黄色片| 国产伦人伦偷精品视频| 真人做人爱边吃奶动态| 神马国产精品三级电影在线观看| 精品一区二区三区视频在线 | 丰满乱子伦码专区| 亚洲第一电影网av| 午夜影院日韩av| 欧美成人一区二区免费高清观看| 成人精品一区二区免费| 日本 av在线| 亚洲人与动物交配视频| 亚洲精品一卡2卡三卡4卡5卡| 19禁男女啪啪无遮挡网站| 成熟少妇高潮喷水视频| ponron亚洲| 国产精品久久视频播放| 国产熟女xx| 宅男免费午夜| 成人鲁丝片一二三区免费| 亚洲中文字幕一区二区三区有码在线看| 观看美女的网站| 欧美最黄视频在线播放免费| 国产不卡一卡二| 18禁黄网站禁片午夜丰满| 日本成人三级电影网站| 一个人免费在线观看电影| 脱女人内裤的视频| www.色视频.com| 欧美中文综合在线视频| 国产私拍福利视频在线观看| xxx96com| 国产国拍精品亚洲av在线观看 | 丰满人妻一区二区三区视频av | 国产主播在线观看一区二区| 中文字幕人妻熟人妻熟丝袜美 | 午夜免费男女啪啪视频观看 | 国产精品一区二区三区四区免费观看 | 91在线精品国自产拍蜜月 | 男人的好看免费观看在线视频| 在线a可以看的网站| 女人十人毛片免费观看3o分钟| 日本熟妇午夜| 亚洲aⅴ乱码一区二区在线播放| www日本在线高清视频| 啪啪无遮挡十八禁网站| 网址你懂的国产日韩在线| 中亚洲国语对白在线视频| 久久香蕉精品热| 性色avwww在线观看| 亚洲国产高清在线一区二区三| 国产精品嫩草影院av在线观看 | 国产高清三级在线| 床上黄色一级片| 啦啦啦韩国在线观看视频| 欧美性猛交黑人性爽| 精品久久久久久久人妻蜜臀av| 亚洲国产欧美人成| 内射极品少妇av片p| 制服丝袜大香蕉在线| 熟妇人妻久久中文字幕3abv| 中亚洲国语对白在线视频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品影视一区二区三区av| 欧美区成人在线视频| 国产亚洲精品久久久com| 欧美三级亚洲精品| 午夜激情福利司机影院| 色在线成人网| 老司机福利观看| 香蕉av资源在线| 最新美女视频免费是黄的| xxxwww97欧美| 18禁黄网站禁片免费观看直播| 精品久久久久久成人av| 99精品欧美一区二区三区四区| 精品国内亚洲2022精品成人| 叶爱在线成人免费视频播放| 白带黄色成豆腐渣| 嫁个100分男人电影在线观看| 亚洲国产高清在线一区二区三| 好男人在线观看高清免费视频| 两个人的视频大全免费| 黄色片一级片一级黄色片| 少妇人妻一区二区三区视频| 欧美+日韩+精品| 国产精品精品国产色婷婷| 亚洲无线观看免费| 夜夜爽天天搞| 人妻夜夜爽99麻豆av| 日韩有码中文字幕| av国产免费在线观看| 99久久综合精品五月天人人| 男女做爰动态图高潮gif福利片| 男女那种视频在线观看| 天堂√8在线中文| 9191精品国产免费久久| 亚洲黑人精品在线| 一级黄色大片毛片| 亚洲人与动物交配视频| 亚洲成人中文字幕在线播放| 亚洲天堂国产精品一区在线| 久久久久久久午夜电影| 内射极品少妇av片p| 欧美成人一区二区免费高清观看| 午夜福利在线观看免费完整高清在 | 午夜视频国产福利| 久久精品人妻少妇| 欧美乱妇无乱码| 国产黄片美女视频| 我要搜黄色片| 18+在线观看网站| 成人特级av手机在线观看| 天堂影院成人在线观看| 日韩欧美精品v在线| 国产欧美日韩一区二区精品| 白带黄色成豆腐渣| 久久6这里有精品| 91九色精品人成在线观看| 老汉色av国产亚洲站长工具| 亚洲天堂国产精品一区在线| 免费av毛片视频| 亚洲男人的天堂狠狠| 国产午夜精品久久久久久一区二区三区 | 少妇人妻一区二区三区视频| 欧美色欧美亚洲另类二区| 在线观看av片永久免费下载| 成人性生交大片免费视频hd| 老汉色∧v一级毛片| 十八禁人妻一区二区| 熟女电影av网| 国模一区二区三区四区视频| 老司机在亚洲福利影院| 色精品久久人妻99蜜桃| 嫩草影视91久久| 欧美一级毛片孕妇| 亚洲av成人av| 欧美精品啪啪一区二区三区| 亚洲最大成人手机在线| 在线视频色国产色| 制服人妻中文乱码| 叶爱在线成人免费视频播放| 在线a可以看的网站| 国产成人av教育| 在线观看一区二区三区| 看免费av毛片| 亚洲欧美日韩卡通动漫| 国产99白浆流出| 欧美中文日本在线观看视频| 久久中文看片网| 国产亚洲精品综合一区在线观看| 欧美zozozo另类| 一区二区三区高清视频在线| 一卡2卡三卡四卡精品乱码亚洲| 男女床上黄色一级片免费看| 99国产精品一区二区三区| a在线观看视频网站| 日本五十路高清| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品国产精品久久久不卡| 日韩精品中文字幕看吧| 99在线视频只有这里精品首页| 日韩欧美精品v在线| 亚洲人成网站高清观看| 香蕉丝袜av| 国产午夜福利久久久久久| 成年免费大片在线观看| 在线免费观看的www视频| 看片在线看免费视频| 欧美黑人巨大hd| 久久久久久九九精品二区国产| 亚洲欧美精品综合久久99| 波多野结衣高清作品| 国产亚洲精品综合一区在线观看| 午夜老司机福利剧场| 精品一区二区三区视频在线观看免费| 国产69精品久久久久777片| 欧美黄色片欧美黄色片| 国产成+人综合+亚洲专区| 毛片女人毛片| 欧美av亚洲av综合av国产av| 啦啦啦观看免费观看视频高清| 亚洲七黄色美女视频| 色哟哟哟哟哟哟| 午夜福利成人在线免费观看| 亚洲内射少妇av| 老熟妇乱子伦视频在线观看| 亚洲人成电影免费在线| 国内揄拍国产精品人妻在线| a级一级毛片免费在线观看| 老鸭窝网址在线观看| 大型黄色视频在线免费观看| 尤物成人国产欧美一区二区三区| 欧美性感艳星| 成年女人看的毛片在线观看| 蜜桃亚洲精品一区二区三区| 免费高清视频大片| 少妇熟女aⅴ在线视频| 欧美日韩福利视频一区二区| 久久亚洲精品不卡| 91在线观看av| 男人和女人高潮做爰伦理| 亚洲成人中文字幕在线播放| 国产精品野战在线观看| 国模一区二区三区四区视频| 日韩欧美国产一区二区入口| 九色国产91popny在线| 免费观看人在逋| 悠悠久久av| 国产精品永久免费网站| 丰满乱子伦码专区| 夜夜看夜夜爽夜夜摸| 特大巨黑吊av在线直播| 天堂动漫精品| 网址你懂的国产日韩在线| 国产日本99.免费观看| 美女黄网站色视频| 国内精品一区二区在线观看| 三级毛片av免费| 国产精品久久久人人做人人爽| 香蕉久久夜色| 一级作爱视频免费观看| 草草在线视频免费看| 国产亚洲精品一区二区www| 欧美乱妇无乱码| 国内揄拍国产精品人妻在线| 韩国av一区二区三区四区| 男插女下体视频免费在线播放| 美女免费视频网站| 99久久精品一区二区三区| 成人无遮挡网站| 少妇人妻一区二区三区视频| 免费人成在线观看视频色| 一边摸一边抽搐一进一小说| 欧美又色又爽又黄视频| 亚洲在线观看片| 色播亚洲综合网| 国产精品久久久久久亚洲av鲁大| 国产综合懂色| 亚洲aⅴ乱码一区二区在线播放| 婷婷精品国产亚洲av在线| 制服人妻中文乱码| 一夜夜www| 久久欧美精品欧美久久欧美| 国产精品久久电影中文字幕| 亚洲av美国av| 91在线观看av|