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

    半穿甲戰(zhàn)斗部對艦船目標(biāo)的侵徹毀傷數(shù)值仿真

    2024-01-03 16:01:40李冬琴劉家昊
    艦船科學(xué)技術(shù) 2023年22期
    關(guān)鍵詞:液艙安定性戰(zhàn)斗部

    李冬琴,張 宇,劉家昊

    (1. 江蘇科技大學(xué) 船舶與海洋工程學(xué)院,江蘇 鎮(zhèn)江 212100;2. 江南造船(集團(tuán))有限責(zé)任公司,上海 201913)

    0 引 言

    目前,大多數(shù)反艦導(dǎo)彈攜帶的是半穿甲戰(zhàn)斗部。半穿甲戰(zhàn)斗部在打擊艦船目標(biāo)時,彈體基于其自身強(qiáng)度和動能侵徹船體防護(hù)板,從而達(dá)到毀傷船體內(nèi)部設(shè)備的目標(biāo)。

    基于有限元仿真軟件Ansys 中的顯示動力學(xué)模塊,研究半穿甲戰(zhàn)斗部打擊艦船目標(biāo)的侵徹深度。在此之前對有限元軟件進(jìn)行仿真的有效性驗(yàn)證,然后導(dǎo)入相關(guān)具體模型進(jìn)行分析計算,研究彈體在不同初速度及不同攻角情況下對艦船的毀傷情況。在侵徹過程中,彈體會發(fā)生嚴(yán)重變形,出現(xiàn)早炸的情況導(dǎo)致武器失效,因此也同時研究分析了彈體侵徹過程中的安定性是否滿足要求。

    1 彈體侵徹液艙的有限元仿真有效性驗(yàn)證

    1.1 幾何模型

    基于Ansys 有限元軟件對艦船模型進(jìn)行分析。首先,依照參考文獻(xiàn)[1]中的試驗(yàn)數(shù)據(jù)從多個角度來驗(yàn)證有限元仿真的有效性。圖1 為該試驗(yàn)裝置示意圖,根據(jù)試驗(yàn)的具體結(jié)構(gòu),將幾何模型簡化為彈體、前靶板、流體區(qū)域和后靶板。幾何模型的具體尺寸按照文獻(xiàn)中試驗(yàn)?zāi)P偷某叽缭O(shè)置,前后靶板采用250 mm×180 mm×1.5 mm 的5A06 鋁合金材料;彈體采用平頭彈的形式,材料為直徑為16 mm 的45#鋼材;流體區(qū)域體積為300 mm×180 mm×250 mm。

    圖1 侵徹試驗(yàn)裝置示意圖Fig. 1 Schematic diagram of the apparatus for the penetration test

    網(wǎng)格劃分方面根據(jù)不同幾何體的體積尺寸,采用不同的網(wǎng)格尺寸。彈體采用3 mm 的四面體網(wǎng)格;由于要觀察前后靶板的破損情況,因此將前后靶板的尺寸設(shè)置為2 mm;流體區(qū)域體積較大,設(shè)置為8 mm 網(wǎng)格。彈體和前后靶板設(shè)置為拉格朗日網(wǎng)格。為保證流體區(qū)域的流動性,將流體區(qū)域設(shè)置為歐拉網(wǎng)格。具體模型如圖2 所示。

    圖2 液艙的有限元模型Fig. 2 Finite element modeling of liquid tank

    1.2 材料參數(shù)及本構(gòu)模型

    彈體和靶板采用Johnson-Cook 本構(gòu)模型。Johnson-Cook 本構(gòu)模型實(shí)質(zhì)上將溫度、應(yīng)變和應(yīng)變率這個3 變量進(jìn)行分開考慮,然后基于乘積關(guān)系將三者對動態(tài)屈服應(yīng)力的影響建立聯(lián)系。其本構(gòu)關(guān)系可表示為[2]:

    式中:A、B、n、C、m分別為初始屈服應(yīng)力常數(shù)、硬化模量、硬化指數(shù)、應(yīng)變率相關(guān)系數(shù)和溫度相關(guān)系數(shù);為材料的等效塑性應(yīng)變;為無量綱的塑性應(yīng)變率。T*=(T-Tr)(Tm-Tr)為相關(guān)的溫度項,Tr和Tm分別為室溫和2024 鋁合金的熔點(diǎn)。流域采用Polynomial 狀態(tài)方程。

    失效模型同樣采用Johnson-Cook 失效模型。當(dāng)損傷度達(dá)到臨界值時,應(yīng)力和壓力取為零值。單元的損傷度D定義為:

    式中,D的取值在0~1 之間,D=0 為初始未損傷,當(dāng)D=1 時,材料發(fā)生失效。Δ?P為一個時間步長的等效塑性應(yīng)變增量; ?f為當(dāng)前時刻的破壞應(yīng)變,其表達(dá)式為:

    表1 固體材料參數(shù)表Tab. 1 Parameter table for solid material

    表2 流體材料參數(shù)表Tab. 2 Parameter table for fluid materials

    1.3 侵徹液艙的仿真驗(yàn)證

    1)空泡效應(yīng)對比

    根據(jù)試驗(yàn)要求,將平頭彈的速度設(shè)置為178.61 m/s。隨著時間的增加,彈體的侵徹深度逐漸增加,同時空泡的直徑也逐漸增加,如圖3(b)所示。通過對比箭頭標(biāo)記位置的直徑發(fā)現(xiàn),試驗(yàn)和仿真產(chǎn)生的空泡大小基本吻合。

    圖3 空泡試驗(yàn)與仿真對比圖Fig. 3 Comparison plot for air bubble test and simulation

    2)靶板破壞形態(tài)對比

    以178.61 m/s 的平頭彈為例,通過對比前靶板的破壞形態(tài)來驗(yàn)證仿真結(jié)果的有效性。圖4 為試驗(yàn)和仿真所產(chǎn)生的彈孔對比圖,在破壞機(jī)理和直徑上基本吻合。

    圖4 靶板的破壞形態(tài)對比圖Fig. 4 Comparison of the diagram pattern of the target plate

    3)水中彈體速度對比

    根據(jù)試驗(yàn)要求分別將仿真計算的彈體速度設(shè)為178.61 m/s 和241.67 m/s,得到如圖5 所示的計算結(jié)果。通過圖像可看出,彈體的初速度越高,其在水中的速度衰減也越大。整體誤差在合理范圍內(nèi),也進(jìn)一步說明了顯示動力學(xué)仿真的可靠性和有效性。

    圖5 剩余速度仿真與試驗(yàn)對比圖Fig. 5 Comparison plot of simnlation and test for residual velocity

    2 彈體對液艙侵徹深度的影響分析

    2.1 典型半穿甲戰(zhàn)斗部及液艙的模型構(gòu)建

    1)半穿甲戰(zhàn)斗部

    以典型的半穿甲尖卵型戰(zhàn)斗部為研究對象。該型戰(zhàn)斗部的殼體材料為45#鋼,總長度1140 mm,殼體頭部長度560 mm,圓柱筒體580 mm,直徑380 mm。內(nèi)裝有105 kg 的TNT 炸藥,戰(zhàn)斗部總質(zhì)量為375 kg。圖6為該典型戰(zhàn)斗部的結(jié)構(gòu)示意圖。

    圖6 典型半穿甲戰(zhàn)斗部的結(jié)構(gòu)簡圖Fig. 6 Structural sketch of a typical semi-armor-piercing combat eleent

    2)液艙結(jié)構(gòu)

    艦船作為軍用船舶,在舷側(cè)設(shè)有一定厚度的防護(hù)液艙。一般防護(hù)液艙由4 層靶板構(gòu)成,分別為舷側(cè)外板、水艙外板、水艙內(nèi)板和防護(hù)縱壁。具體結(jié)構(gòu)如圖7所示[4]。

    圖7 艦船舷側(cè)防護(hù)液艙示意圖Fig. 7 Schematic diagram of the ship's side protective liquid tanks

    一般的船用鋼板按結(jié)構(gòu)可分為均質(zhì)板和加筋板,艦船的液艙靶板皆采用加強(qiáng)筋板架結(jié)構(gòu),其面板厚度25 mm,全長6000 mm,寬4000 mm,面板寬160 mm,厚16 mm,腹板高240 mm,小筋腹板高180 mm,厚8 mm。具體模型如圖8 所示。

    圖8 加強(qiáng)筋板架結(jié)構(gòu)模型圖Fig. 8 Structural model diagram of reinforced rib plate frame

    將采用厚度等效法將加筋板等效為均質(zhì)板,為后續(xù)的仿真實(shí)驗(yàn)做準(zhǔn)備。厚度等效法是將靶板上的T 型鋼和L 型鋼等結(jié)構(gòu),按照強(qiáng)度或質(zhì)量等效原則等效到靶板的厚度上,從而將加筋板轉(zhuǎn)化為均質(zhì)板來研究。

    設(shè)平面靶板與靶板上的加強(qiáng)筋為相同材質(zhì),因此它們的密度相等。加強(qiáng)筋的體積為:

    式中,SH、SG分別為桁材的面積和骨材的面積,hH、hG分別為桁材的高度和骨材的高度。

    加強(qiáng)筋的等效厚度為:

    式中,SB為平面靶板的面積。

    最終,加強(qiáng)筋板等效為均質(zhì)靶板的厚度為:

    式中,HB為平面靶板的厚度。

    根據(jù)上述理論,該典型艦船的防護(hù)液艙靶板的厚度約為30 mm,包括半穿甲戰(zhàn)斗部及靶板的其他尺寸均按原尺寸建立有限元模型。彈體和靶板分別設(shè)置為30 mm 的四面體網(wǎng)格和六面體網(wǎng)格。同時設(shè)置為拉格朗日網(wǎng)格;液體區(qū)域設(shè)置為60 mm 的四面體網(wǎng)格,同時設(shè)置為歐拉網(wǎng)格。每塊靶板的4 個側(cè)面設(shè)置固定約束,并對歐拉域進(jìn)行體積控制。具體有限元模型如圖9所示。

    圖9 船體舷側(cè)防護(hù)液艙有限元模型圖Fig. 9 Finite element model of hull's outboard protection tank

    彈體和靶板同樣采用J-C 本構(gòu)方程,水和空氣將分別采用Shock EOS Linear 和Ideal Gas EOS 狀態(tài)方程,相關(guān)材料參數(shù)如表3 和表4 所示[5-8]:

    表3 固體材料參數(shù)表Tab. 3 Parameter table of solid materials

    表4 流體材料參數(shù)表Tab. 4 Parameter table of fluid materials

    2.2 不同初速度下的侵徹深度分析

    圖10 為不同初速度下彈體在30 ms 后的侵徹深度。首先,由于戰(zhàn)斗部為尖卵型戰(zhàn)斗部,所以靶板的破口為花瓣型破口。其次,戰(zhàn)斗部在侵徹液艙時,產(chǎn)生了很明顯的空泡現(xiàn)象,在不同初速度下,戰(zhàn)斗部的侵徹深度存在明顯差異。當(dāng)初速度為300 m/s 時,戰(zhàn)斗部在30 ms 后只能侵徹到液艙內(nèi)部便不能繼續(xù)侵徹;當(dāng)初速度為400 m/s 時,戰(zhàn)斗部可擊穿液艙,但剩余動能無法侵徹防護(hù)縱壁;當(dāng)初速度為500 m/s 時,戰(zhàn)斗部剛好可侵徹整個舷側(cè)結(jié)構(gòu);當(dāng)初速度設(shè)為600 m/s 時,戰(zhàn)斗部可完全侵徹到舷側(cè)結(jié)構(gòu),并還有一定的剩余速度。因此可得出結(jié)論,當(dāng)彈體和目標(biāo)結(jié)構(gòu)確定的條件下,彈體初速度越高,侵徹深度和剩余速度越大。

    圖10 不同初速度下的侵徹深度Fig. 10 Depth of penetration at different initial velocities

    圖11 為不同初速度下彈體剩余速度的時歷曲線,可以看出,隨著計算時間的推移,彈體的剩余速度呈階梯式下降,這是由于彈體在侵徹4 層靶板時,每侵徹一層,彈體的剩余速度都會陡然下降。從不同的彈體初速度角度分析,初速度越高,彈體的剩余速度也越高。而初速度為300 m/s 和400 m/s 時,由于彈體沒有足夠多的動能,因此最終它們未能擊穿液艙靶板,其剩余速度皆小于0。當(dāng)彈體的初速度為300 m/s 時,由于30 ms 的時間彈體沒有穿透靶板,并且剩余速度仍大于0,因此將計算時間延長到40 ms,直到彈體的擊穿靶板或彈體速度小于0。

    圖11 不同初速度下彈體剩余速度的時歷曲線Fig. 11 Time history curves of the residual velocity of the projectile at different initial velocities

    從不同初速度的時歷曲線可以看出,彈體的初速度越高,最終的剩余速度也越高,當(dāng)彈體的初速度呈線性變化時,彈體的最終剩余速度并不呈現(xiàn)線性變化。若想獲得在一定區(qū)間內(nèi)任意初速度所對應(yīng)的剩余速度,只通過有限元仿真獲得,工作量巨大很難實(shí)現(xiàn)。因此,在計算每個初速度所對應(yīng)的剩余速度時,可采用近似擬合的方法計算。通過仿真實(shí)驗(yàn)獲得較多的數(shù)據(jù),然后將結(jié)果擬合成函數(shù)曲線求得在一定區(qū)間范圍內(nèi)任意初速度所對應(yīng)的剩余速度,從而根據(jù)剩余速度判斷出彈體能否擊穿舷側(cè)液艙靶板結(jié)構(gòu)。

    2.3 不同攻角下的侵徹深度分析

    彈體在打擊艦船目標(biāo)時,由于彈體飛行姿態(tài)的調(diào)整或者艦船目標(biāo)的運(yùn)動,彈體的速度方向不會總是與靶板的法線方向平行,速度方向和靶板的法線方向會形成一個夾角。討論在相同初速度條件下,不同攻角下彈體的侵徹深度及剩余速度。

    圖12 為初速度為500 m/s 時,30 ms 后彈體的侵徹深度??芍?,隨著攻角的增加,彈體的侵徹深度越來越小。為了進(jìn)一步確定侵徹深度,同樣以剩余速度作為彈體能否擊穿靶板的判斷依據(jù)。

    圖12 不同攻角下的侵徹深度Fig. 12 Depth of penetration at different argles of attack

    圖13 為相同初速度條件下,不同攻角的剩余速度時歷曲線??芍?,初速度相等的條件下,前15 ms的速度下降趨勢變化不大。在15~30 ms 時的速度變化存在明顯差異,隨著攻角的線性增加,彈體的剩余速度逐漸減小,剩余速度減小量呈非線性變化。

    圖13 不同攻角下的剩余速度時歷曲線Fig. 13 Residual velocity time history curves at different angles of attack

    3 彈體侵徹過程的安定性分析

    彈體的裝藥安定性是指彈體在受外界因素影響下,仍能保持其性質(zhì)的能力。國內(nèi)外學(xué)者在進(jìn)行彈體安定性分析中,大多采用有限元仿真的方法進(jìn)行彈體安定性評估。焦志剛等[9]基于Ls-dyna 進(jìn)行半穿甲戰(zhàn)斗部斜侵不同厚度鋼靶板的安定性分析,通過計算得到了殼體和裝藥的應(yīng)力曲線。王偉力等[10]同樣基于Lsdyna 分析了不同傾角下的彈體受力情況,得出了小傾角對安定性影響不大,大傾角影響較大的結(jié)論。本文主要研究的是彈體產(chǎn)生小傾角的情況,因此,將基于文獻(xiàn)[11]的安定性計算方法進(jìn)行彈體的安定性分析。

    3.1 安定性理論計算

    侵徹過程中,裝藥受到應(yīng)力的作用產(chǎn)生塑性應(yīng)變能。用應(yīng)力和應(yīng)變率表示單位體積的塑性應(yīng)變能變化率的能量守恒方程為:

    式中: ρ為材料密度;e為單位體積材料的應(yīng)變能;σii和 εii分別為各個方向的應(yīng)力分量和相應(yīng)的應(yīng)變率分量。

    將(7)兩邊積分得到:

    因?yàn)棣?σ(t),ε=ε(t),從而可得σ=σ(ε)。所以上式可寫為:

    式中: ε0為初始塑性應(yīng)變; εmax為最大應(yīng)變。取侵徹過程中,裝藥受到的最大應(yīng)力σmax和最大應(yīng)變值εmax,則有:

    在侵徹過程中,由于相互作用的時間非常短,可認(rèn)為裝藥的塑性應(yīng)變能轉(zhuǎn)化為內(nèi)能是一個絕熱升溫的過程,故有方程:

    式中:E為內(nèi)能;m為質(zhì)量;ρ為材料密度;V為體積;C''為比熱容;T0為初始溫度。

    由上兩式得:

    取T0=20℃,ρ=1.7 g/cm3,C''=1.372 J/(g·K)為TNT 的比熱容。在得知最大應(yīng)力σmax和最大應(yīng)變εmax情況下,可計算出裝藥的最大升溫。TNT 的爆點(diǎn)溫度為240℃,若計算溫度大于240℃時,彈體內(nèi)部裝藥由于壓縮變形過大,內(nèi)能升高導(dǎo)致彈體提前爆炸,安定性不符合要求;反之則滿足要求。

    3.2 彈體的最大應(yīng)力應(yīng)變分析

    圖14 為基于Ansys 數(shù)值仿真的彈體應(yīng)力應(yīng)變云圖??芍饴研蛻?zhàn)斗部的尖部在撞擊靶板時產(chǎn)生了較大的壓縮變形,但這對彈體內(nèi)部的裝藥影響程度較小。除此之外,隨著初速度的增加,彈體的應(yīng)變整體呈上升趨勢,應(yīng)變較大的區(qū)域面積也有所增大。

    圖14 不同初速度下的彈體應(yīng)變云圖Fig. 14 Strain map of the projectile at different initial velocities

    圖15 為彈體應(yīng)力應(yīng)變云圖。可知,隨著彈體攻角的變化,彈體的應(yīng)力分布也存在差異。應(yīng)力應(yīng)變最大的位置,在彈體的由圓柱收縮為尖卵的拐角處。一旦此處的變形過大,將會對內(nèi)部裝藥的安定性產(chǎn)生巨大影響。因此,可通過計算該位置的應(yīng)力應(yīng)變大小來評判彈體的安定性是否符合要求。

    圖15 不同攻角下的彈體應(yīng)變云圖Fig. 15 Strain maps of the projectlle at different angles of attack

    由于裝藥充滿整個彈體內(nèi)部,因此彈體在侵徹過程中所受的最大應(yīng)變可近似為內(nèi)部裝藥所受的最大應(yīng)變。得出裝藥所受最大應(yīng)變后,再求得最大應(yīng)力便可通過式(12)得到彈體在侵徹過程中,由于碰撞所產(chǎn)生的內(nèi)能,從而進(jìn)一步判斷彈體的安定性是否滿足要求。

    圖16 為TNT 的應(yīng)力應(yīng)變曲線,將仿真計算得到的最大應(yīng)變值代入到應(yīng)力應(yīng)變曲線中,即可得到相應(yīng)的最大應(yīng)力值。通過對比得知,該型半穿甲戰(zhàn)斗部在一定范圍內(nèi)改變其初速度和攻角其安定性均滿足要求[11]。

    圖16 TNT 的應(yīng)力應(yīng)變曲線Fig. 16 Stress-strain curve of TNT

    圖17 為彈體產(chǎn)生最大應(yīng)力應(yīng)變時所在的位置。可知,彈體在不同的攻角下,最大應(yīng)力應(yīng)變的位置存在差異。當(dāng)攻角為0°和5°時,彈體在侵徹第三層靶板時產(chǎn)生最大應(yīng)力;當(dāng)攻角為10°和15°時,彈體在侵徹第二層靶板產(chǎn)生最大應(yīng)力。當(dāng)彈體的安定性不符合要求時,便可得出引爆的相對位置,進(jìn)而判斷能否對艦船目標(biāo)的內(nèi)部設(shè)備造成毀傷。

    圖17 最大應(yīng)力應(yīng)變時的侵徹位置Fig. 17 Position of penetration at maximum stress-strain

    4 結(jié) 語

    本文基于Ansys 有限元仿真軟件進(jìn)行半穿甲戰(zhàn)斗部侵徹舷側(cè)液艙的仿真實(shí)驗(yàn),通過改變彈體不同初速度與不同攻角進(jìn)行仿真分析,得出以下結(jié)論:

    當(dāng)彈體和目標(biāo)結(jié)構(gòu)確定的條件下,初速度越高,侵徹深度和剩余速度越大。隨著計算時間的推移,彈體的剩余速度呈階梯式下降。當(dāng)彈體的初速度呈線性變化時,彈體的最終剩余速度并不呈現(xiàn)線性變化。隨著初速度的增加,彈體的應(yīng)變整體呈上升趨勢,應(yīng)變較大的區(qū)域面積也有所增大。

    隨著攻角的增加,彈體的剩余速度越來越小,而隨著攻角的線性增加,彈體剩余速度減小量呈非線性減小。隨著彈體攻角的變化,彈體的應(yīng)力分布也存在差異,應(yīng)力應(yīng)變最大位置,在彈體由圓柱收縮為尖卵的拐角處。最后根據(jù)彈體的最大應(yīng)力應(yīng)變仿真結(jié)果分析得到,裝藥安定性均滿足要求。

    猜你喜歡
    液艙安定性戰(zhàn)斗部
    B型LNG液艙支座縱骨趾端處表面裂紋擴(kuò)展計算
    基于CFD的大型船舶液艙晃蕩研究
    基于Android系統(tǒng)的汽油氧化安定性測定儀的設(shè)計
    電子制作(2018年11期)2018-08-04 03:25:50
    烴族組成對柴油儲存安定性影響研究
    考慮晃蕩效應(yīng)的獨(dú)立B型LNG液艙結(jié)構(gòu)多目標(biāo)優(yōu)化
    海洋工程(2016年2期)2016-10-12 05:08:07
    隨動定向戰(zhàn)斗部的抗振動沖擊控制研究
    FPSO與運(yùn)輸船旁靠時液艙晃蕩與船舶運(yùn)動耦合效應(yīng)分析
    幾種結(jié)構(gòu)的三段離散桿戰(zhàn)斗部數(shù)值模擬
    戰(zhàn)斗部遠(yuǎn)場水下爆炸對艦船沖擊損傷評估
    關(guān)于水泥檢驗(yàn)中安定性的時效問題
    国产成人影院久久av| 欧美激情在线99| 久久久久九九精品影院| 国产97色在线日韩免费| 9191精品国产免费久久| 亚洲电影在线观看av| 色噜噜av男人的天堂激情| 久久久久久久午夜电影| 午夜福利视频1000在线观看| 午夜福利成人在线免费观看| 淫妇啪啪啪对白视频| 国产欧美日韩精品一区二区| 国产激情偷乱视频一区二区| 国产aⅴ精品一区二区三区波| 日韩欧美在线二视频| 午夜久久久久精精品| 高潮久久久久久久久久久不卡| 国产精品一区二区三区四区久久| 欧美bdsm另类| av欧美777| 久久久久久久久大av| 日韩欧美一区二区三区在线观看| av中文乱码字幕在线| 在线观看66精品国产| 国产免费一级a男人的天堂| 好男人电影高清在线观看| 亚洲自拍偷在线| 色综合亚洲欧美另类图片| 少妇的逼好多水| 久久精品国产亚洲av香蕉五月| 他把我摸到了高潮在线观看| 亚洲18禁久久av| 国产精品一区二区三区四区久久| 亚洲欧美一区二区三区黑人| 成年人黄色毛片网站| 亚洲一区高清亚洲精品| 亚洲人成网站高清观看| 无限看片的www在线观看| 精品一区二区三区视频在线观看免费| 色在线成人网| 日本精品一区二区三区蜜桃| 亚洲欧美日韩卡通动漫| 午夜福利视频1000在线观看| 三级毛片av免费| 国产av一区在线观看免费| 精品久久久久久久人妻蜜臀av| 免费看光身美女| 免费无遮挡裸体视频| 精品国产亚洲在线| 人人妻人人澡欧美一区二区| 一进一出好大好爽视频| 久久精品综合一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 最近在线观看免费完整版| 亚洲avbb在线观看| 99热精品在线国产| 少妇高潮的动态图| 久久久久久人人人人人| 在线视频色国产色| 亚洲av免费在线观看| 淫秽高清视频在线观看| 69av精品久久久久久| 日韩有码中文字幕| 精品人妻一区二区三区麻豆 | 嫁个100分男人电影在线观看| tocl精华| 久久草成人影院| 中文字幕人妻熟人妻熟丝袜美 | 久久久久久久久久黄片| 亚洲七黄色美女视频| 男人和女人高潮做爰伦理| 亚洲第一欧美日韩一区二区三区| 免费在线观看日本一区| 欧美乱色亚洲激情| 国内精品美女久久久久久| 成人三级黄色视频| 97超视频在线观看视频| 搡老妇女老女人老熟妇| 一边摸一边抽搐一进一小说| 黄色片一级片一级黄色片| 级片在线观看| 亚洲欧美日韩东京热| 亚洲电影在线观看av| www日本黄色视频网| 国产一区二区亚洲精品在线观看| 免费人成视频x8x8入口观看| 18禁黄网站禁片午夜丰满| 中文字幕人妻熟人妻熟丝袜美 | 制服丝袜大香蕉在线| 男人和女人高潮做爰伦理| www.熟女人妻精品国产| 免费av观看视频| 又黄又爽又免费观看的视频| 精品一区二区三区人妻视频| eeuss影院久久| 少妇的逼好多水| 亚洲黑人精品在线| x7x7x7水蜜桃| 天堂av国产一区二区熟女人妻| 手机成人av网站| 国产爱豆传媒在线观看| 午夜福利欧美成人| 麻豆国产97在线/欧美| 日韩欧美免费精品| 久久久色成人| 免费在线观看亚洲国产| 国产精品 欧美亚洲| 天堂av国产一区二区熟女人妻| 少妇高潮的动态图| 国产精品综合久久久久久久免费| 97人妻精品一区二区三区麻豆| 长腿黑丝高跟| 亚洲五月天丁香| 又爽又黄无遮挡网站| 99热精品在线国产| 欧美高清成人免费视频www| 一个人免费在线观看的高清视频| 亚洲av成人av| 午夜亚洲福利在线播放| 日韩高清综合在线| 久久精品国产99精品国产亚洲性色| 欧美一级a爱片免费观看看| 色综合婷婷激情| 在线天堂最新版资源| 美女大奶头视频| 国产一区二区亚洲精品在线观看| 欧美成人一区二区免费高清观看| 亚洲人与动物交配视频| 观看免费一级毛片| 在线a可以看的网站| eeuss影院久久| 久久精品影院6| 午夜激情欧美在线| h日本视频在线播放| 18禁国产床啪视频网站| 999久久久精品免费观看国产| 免费av毛片视频| 一边摸一边抽搐一进一小说| 日日摸夜夜添夜夜添小说| 欧美一区二区精品小视频在线| 国产一区二区在线av高清观看| 夜夜夜夜夜久久久久| 在线观看舔阴道视频| www日本黄色视频网| 高清毛片免费观看视频网站| 欧美乱码精品一区二区三区| 国产精品香港三级国产av潘金莲| 51国产日韩欧美| 久久精品亚洲精品国产色婷小说| 日韩欧美三级三区| e午夜精品久久久久久久| 好看av亚洲va欧美ⅴa在| 搡老岳熟女国产| 欧美xxxx黑人xx丫x性爽| 观看免费一级毛片| 村上凉子中文字幕在线| 日本一二三区视频观看| 欧美日韩瑟瑟在线播放| 国产一区二区在线av高清观看| 国产单亲对白刺激| 日韩中文字幕欧美一区二区| 日日干狠狠操夜夜爽| 国产高清三级在线| 一a级毛片在线观看| 18禁裸乳无遮挡免费网站照片| 久久久久久久亚洲中文字幕 | 搡女人真爽免费视频火全软件 | 日本黄色视频三级网站网址| 观看美女的网站| 国产久久久一区二区三区| 三级国产精品欧美在线观看| 成人午夜高清在线视频| 真实男女啪啪啪动态图| 美女大奶头视频| 亚洲五月婷婷丁香| 一本综合久久免费| 亚洲av美国av| 亚洲avbb在线观看| 人人妻人人澡欧美一区二区| 99在线视频只有这里精品首页| 无限看片的www在线观看| 国产91精品成人一区二区三区| 日韩有码中文字幕| 国产毛片a区久久久久| 国产精品久久久人人做人人爽| 精品一区二区三区视频在线 | 女人十人毛片免费观看3o分钟| 级片在线观看| 亚洲国产精品sss在线观看| 日韩高清综合在线| 两性午夜刺激爽爽歪歪视频在线观看| 欧美一区二区国产精品久久精品| 观看免费一级毛片| 久久精品国产自在天天线| 亚洲成av人片免费观看| 国产毛片a区久久久久| 十八禁网站免费在线| 欧美丝袜亚洲另类 | 久久久久国产精品人妻aⅴ院| 美女黄网站色视频| 很黄的视频免费| 国产黄片美女视频| 国产高清激情床上av| 亚洲黑人精品在线| 欧美日韩亚洲国产一区二区在线观看| netflix在线观看网站| 在线天堂最新版资源| 欧美午夜高清在线| 国产精品久久视频播放| 亚洲avbb在线观看| 免费av不卡在线播放| 91麻豆精品激情在线观看国产| 亚洲精品影视一区二区三区av| 日本在线视频免费播放| 在线观看免费视频日本深夜| 精品国产三级普通话版| 亚洲aⅴ乱码一区二区在线播放| 一卡2卡三卡四卡精品乱码亚洲| 免费无遮挡裸体视频| www.色视频.com| 欧美成人一区二区免费高清观看| 九色国产91popny在线| 成人特级黄色片久久久久久久| 国产精品一及| 成熟少妇高潮喷水视频| 99国产极品粉嫩在线观看| 日本 欧美在线| 在线播放无遮挡| 热99在线观看视频| 国产97色在线日韩免费| 全区人妻精品视频| 欧美成人性av电影在线观看| 国产精品野战在线观看| 国产高清激情床上av| 午夜精品久久久久久毛片777| 757午夜福利合集在线观看| 国产一区二区在线av高清观看| 国产成人啪精品午夜网站| 国产高潮美女av| 69人妻影院| 国产探花在线观看一区二区| 日韩欧美三级三区| 成人18禁在线播放| 亚洲精品成人久久久久久| 欧美成狂野欧美在线观看| 欧美高清成人免费视频www| 日韩精品青青久久久久久| 免费av观看视频| 一级a爱片免费观看的视频| 狠狠狠狠99中文字幕| 久久国产精品人妻蜜桃| 精品电影一区二区在线| 琪琪午夜伦伦电影理论片6080| 成人av在线播放网站| 亚洲在线自拍视频| 午夜福利高清视频| 一个人免费在线观看的高清视频| 午夜激情欧美在线| 欧美+日韩+精品| 在线观看av片永久免费下载| 亚洲精品美女久久久久99蜜臀| 亚洲人成伊人成综合网2020| 色老头精品视频在线观看| 亚洲,欧美精品.| 好男人在线观看高清免费视频| 亚洲第一欧美日韩一区二区三区| 亚洲国产精品合色在线| 色综合亚洲欧美另类图片| 最近在线观看免费完整版| 9191精品国产免费久久| 久久精品夜夜夜夜夜久久蜜豆| 男人和女人高潮做爰伦理| 精品久久久久久久人妻蜜臀av| 亚洲国产精品sss在线观看| 极品教师在线免费播放| 成熟少妇高潮喷水视频| 精品国产亚洲在线| 深爱激情五月婷婷| 久久草成人影院| 国产精品影院久久| av在线天堂中文字幕| 色av中文字幕| 欧美日韩黄片免| 国产精品精品国产色婷婷| 最新在线观看一区二区三区| 激情在线观看视频在线高清| 村上凉子中文字幕在线| 亚洲五月天丁香| 级片在线观看| 身体一侧抽搐| 欧美黑人欧美精品刺激| 国产成人影院久久av| 精品久久久久久久人妻蜜臀av| 国产精品香港三级国产av潘金莲| 男人和女人高潮做爰伦理| 国产色爽女视频免费观看| 91久久精品电影网| 国产成人福利小说| 中亚洲国语对白在线视频| 久久久久久久亚洲中文字幕 | 精品人妻一区二区三区麻豆 | 亚洲无线观看免费| 亚洲成人中文字幕在线播放| 最后的刺客免费高清国语| 怎么达到女性高潮| 色在线成人网| 日日干狠狠操夜夜爽| 成熟少妇高潮喷水视频| 国产一区二区三区视频了| 亚洲精品色激情综合| 好男人在线观看高清免费视频| 久久久久久久久中文| 亚洲乱码一区二区免费版| 观看美女的网站| 成人18禁在线播放| 波多野结衣高清无吗| 一区二区三区免费毛片| 久久人妻av系列| 在线天堂最新版资源| 国产精品,欧美在线| 亚洲精华国产精华精| 亚洲国产精品久久男人天堂| 波野结衣二区三区在线 | 午夜免费激情av| 成人18禁在线播放| 在线天堂最新版资源| 欧美性猛交╳xxx乱大交人| 美女cb高潮喷水在线观看| 国产综合懂色| 三级国产精品欧美在线观看| 噜噜噜噜噜久久久久久91| 人妻久久中文字幕网| 白带黄色成豆腐渣| 哪里可以看免费的av片| 中出人妻视频一区二区| 亚洲专区中文字幕在线| 99久久无色码亚洲精品果冻| 两个人看的免费小视频| 亚洲专区中文字幕在线| 嫩草影院精品99| 亚洲国产色片| 人妻久久中文字幕网| 亚洲精品一卡2卡三卡4卡5卡| 一进一出好大好爽视频| 99riav亚洲国产免费| 亚洲不卡免费看| 国产一级毛片七仙女欲春2| 他把我摸到了高潮在线观看| 免费大片18禁| 99精品久久久久人妻精品| 国产探花在线观看一区二区| 色综合欧美亚洲国产小说| 观看免费一级毛片| 精品福利观看| 3wmmmm亚洲av在线观看| 亚洲av电影在线进入| 18美女黄网站色大片免费观看| 精品久久久久久成人av| 亚洲中文日韩欧美视频| 国产爱豆传媒在线观看| 少妇的逼好多水| 黄片大片在线免费观看| 国产综合懂色| 亚洲av二区三区四区| 搡老岳熟女国产| 18禁美女被吸乳视频| 国产精品久久电影中文字幕| 在线天堂最新版资源| 哪里可以看免费的av片| 欧美成狂野欧美在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 男女之事视频高清在线观看| 两人在一起打扑克的视频| 午夜日韩欧美国产| 老熟妇仑乱视频hdxx| 变态另类丝袜制服| 可以在线观看的亚洲视频| 国产精品野战在线观看| 亚洲专区中文字幕在线| 欧美日韩中文字幕国产精品一区二区三区| 日韩有码中文字幕| 中文字幕高清在线视频| 日韩免费av在线播放| 亚洲av日韩精品久久久久久密| 别揉我奶头~嗯~啊~动态视频| 天堂网av新在线| 99久久久亚洲精品蜜臀av| 久久精品影院6| 热99re8久久精品国产| 中文资源天堂在线| 国产精品久久久久久久电影 | 国产淫片久久久久久久久 | 别揉我奶头~嗯~啊~动态视频| 国产又黄又爽又无遮挡在线| 日韩人妻高清精品专区| 免费大片18禁| 女人被狂操c到高潮| 中文资源天堂在线| 2021天堂中文幕一二区在线观| 黄色视频,在线免费观看| 69人妻影院| 日本在线视频免费播放| 熟女电影av网| 18禁黄网站禁片免费观看直播| 免费观看人在逋| 特大巨黑吊av在线直播| 无限看片的www在线观看| 男女下面进入的视频免费午夜| tocl精华| 午夜精品一区二区三区免费看| 又粗又爽又猛毛片免费看| 哪里可以看免费的av片| 性欧美人与动物交配| 免费看光身美女| 中文字幕人妻丝袜一区二区| 在线观看日韩欧美| 午夜免费男女啪啪视频观看 | 国语自产精品视频在线第100页| 他把我摸到了高潮在线观看| 成年人黄色毛片网站| 欧美高清成人免费视频www| 国产主播在线观看一区二区| 桃红色精品国产亚洲av| 亚洲一区二区三区不卡视频| 亚洲av成人不卡在线观看播放网| 淫秽高清视频在线观看| 国产97色在线日韩免费| 国产探花极品一区二区| 亚洲精品日韩av片在线观看 | 欧美色视频一区免费| 波多野结衣巨乳人妻| 人人妻人人澡欧美一区二区| 我要搜黄色片| 欧洲精品卡2卡3卡4卡5卡区| 精品国产超薄肉色丝袜足j| 天堂动漫精品| 亚洲人成伊人成综合网2020| 特级一级黄色大片| 草草在线视频免费看| 国模一区二区三区四区视频| 欧洲精品卡2卡3卡4卡5卡区| 欧美激情久久久久久爽电影| 日本成人三级电影网站| 哪里可以看免费的av片| 老司机福利观看| 悠悠久久av| 51午夜福利影视在线观看| 蜜桃久久精品国产亚洲av| 最近最新免费中文字幕在线| 午夜两性在线视频| 床上黄色一级片| 国产成人欧美在线观看| 久久香蕉精品热| 极品教师在线免费播放| 一卡2卡三卡四卡精品乱码亚洲| 亚洲一区二区三区色噜噜| 免费看美女性在线毛片视频| 精品久久久久久成人av| 每晚都被弄得嗷嗷叫到高潮| а√天堂www在线а√下载| 国产成人a区在线观看| av国产免费在线观看| 毛片女人毛片| 熟妇人妻久久中文字幕3abv| 香蕉av资源在线| 国产精品美女特级片免费视频播放器| netflix在线观看网站| 热99在线观看视频| 在线观看免费午夜福利视频| 99国产精品一区二区三区| 亚洲av免费在线观看| 久久久久性生活片| 欧美最新免费一区二区三区 | 国产探花极品一区二区| 国产一区二区亚洲精品在线观看| www.熟女人妻精品国产| 最好的美女福利视频网| 观看美女的网站| 国产激情欧美一区二区| 欧美午夜高清在线| 色综合欧美亚洲国产小说| 亚洲精品亚洲一区二区| 国模一区二区三区四区视频| 韩国av一区二区三区四区| 中文字幕人妻丝袜一区二区| 亚洲不卡免费看| 欧美又色又爽又黄视频| 超碰av人人做人人爽久久 | a级一级毛片免费在线观看| av国产免费在线观看| 一进一出抽搐动态| 久久香蕉精品热| 亚洲av第一区精品v没综合| 国产欧美日韩精品一区二区| 神马国产精品三级电影在线观看| 最近最新免费中文字幕在线| 国产伦在线观看视频一区| 高清日韩中文字幕在线| 中文字幕精品亚洲无线码一区| 全区人妻精品视频| 制服丝袜大香蕉在线| 国产色婷婷99| 人人妻,人人澡人人爽秒播| 久久性视频一级片| 欧美不卡视频在线免费观看| 精品久久久久久久毛片微露脸| 淫妇啪啪啪对白视频| 国产免费av片在线观看野外av| 97超级碰碰碰精品色视频在线观看| 全区人妻精品视频| 免费观看人在逋| 中文字幕av在线有码专区| 美女黄网站色视频| 精品午夜福利视频在线观看一区| 国内少妇人妻偷人精品xxx网站| 亚洲国产精品成人综合色| 欧美+日韩+精品| 久久久久久大精品| 天堂网av新在线| 18禁国产床啪视频网站| 日本三级黄在线观看| 日韩欧美 国产精品| 18禁裸乳无遮挡免费网站照片| 久久精品国产亚洲av涩爱 | 久久国产精品人妻蜜桃| 亚洲成人中文字幕在线播放| 麻豆成人午夜福利视频| 亚洲av二区三区四区| 最近视频中文字幕2019在线8| 欧美一级a爱片免费观看看| 18禁黄网站禁片免费观看直播| 在线观看日韩欧美| 在线a可以看的网站| 少妇裸体淫交视频免费看高清| 99热精品在线国产| 黄色女人牲交| 特大巨黑吊av在线直播| 亚洲av成人不卡在线观看播放网| 十八禁人妻一区二区| 麻豆一二三区av精品| 他把我摸到了高潮在线观看| 久久草成人影院| 99国产极品粉嫩在线观看| 十八禁人妻一区二区| 午夜日韩欧美国产| 国产国拍精品亚洲av在线观看 | 特级一级黄色大片| 综合色av麻豆| 国内精品久久久久精免费| 美女高潮喷水抽搐中文字幕| 内地一区二区视频在线| 18禁黄网站禁片免费观看直播| 欧美色欧美亚洲另类二区| 国产精品综合久久久久久久免费| 国产亚洲精品久久久com| 中文亚洲av片在线观看爽| 国产高清有码在线观看视频| 日韩国内少妇激情av| 午夜福利免费观看在线| 国产伦精品一区二区三区四那| 床上黄色一级片| 窝窝影院91人妻| 国产一区二区亚洲精品在线观看| 脱女人内裤的视频| 国产一区二区在线av高清观看| 国内久久婷婷六月综合欲色啪| 午夜福利欧美成人| 久久久国产成人精品二区| 2021天堂中文幕一二区在线观| 国产视频内射| 伊人久久大香线蕉亚洲五| 国产精品美女特级片免费视频播放器| 亚洲一区二区三区色噜噜| 男人的好看免费观看在线视频| 波野结衣二区三区在线 | 精品熟女少妇八av免费久了| 欧美三级亚洲精品| 青草久久国产| 亚洲中文字幕一区二区三区有码在线看| 亚洲精华国产精华精| 又黄又爽又免费观看的视频| 免费观看人在逋| 欧美日本视频| 久久久精品欧美日韩精品| 91九色精品人成在线观看| 久久久国产精品麻豆| 日韩欧美在线乱码| 国产精品日韩av在线免费观看| a级一级毛片免费在线观看| 色综合站精品国产| 波多野结衣巨乳人妻| 成人特级av手机在线观看| 国产欧美日韩精品一区二区| 琪琪午夜伦伦电影理论片6080| 亚洲人成网站高清观看| 日韩成人在线观看一区二区三区| 久久久久九九精品影院| 成人特级黄色片久久久久久久| 国产精品 国内视频| 国产精品98久久久久久宅男小说| а√天堂www在线а√下载| 看免费av毛片| 美女高潮喷水抽搐中文字幕| 高潮久久久久久久久久久不卡| 3wmmmm亚洲av在线观看| 午夜精品久久久久久毛片777| 国产免费一级a男人的天堂| 亚洲国产色片| 亚洲 国产 在线| 老司机午夜十八禁免费视频| 人人妻人人澡欧美一区二区| 最新在线观看一区二区三区| 精品久久久久久,| 少妇的丰满在线观看| 啦啦啦免费观看视频1|