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

    彈體正侵徹混凝土靶動態(tài)開坑作用過程*

    2023-09-23 09:15:52李鵬程張先鋒王桂吉劉均偉鄧宇軒
    爆炸與沖擊 2023年9期
    關鍵詞:彈頭彈體流線

    李鵬程,張先鋒,王桂吉,劉 闖,劉均偉,鄧宇軒,盛 強

    (1.南京理工大學機械工程學院,江蘇 南京 210094;2.中國工程物理研究院流體物理研究所,四川 綿陽 621999)

    隨著防御系統(tǒng)的不斷發(fā)展,混凝土作為防護結構材料廣泛應用于軍事及民用領域,其中地面多層建筑物、地下軍事設施(如指揮所,導彈發(fā)射井等)依靠復雜的混凝土復合結構而具有優(yōu)良的防護性能。鉆地武器侵徹該類混凝土目標的終點效應一直是新概念武器研發(fā)和工程防護領域關注的熱點問題[1]。無論是混凝土厚靶還是薄靶,當受到彈體撞擊時,迎彈面均會產生明顯的開坑破壞效應。迎彈面開坑階段是完整侵徹過程的初始階段,對于混凝土薄靶的貫穿,迎彈面開坑階段更是整個貫穿過程中的重要組成部分[2],此時開坑破壞效應將直接反應彈體的侵徹毀傷能力。因此,研究彈體侵徹混凝土的開坑作用過程,對于全面掌握彈體侵徹作用過程以及定量分析開坑破壞效應均具有十分重要的意義,可為彈體非正侵徹混凝土靶的開坑作用過程研究提供參考。

    針對彈體侵徹混凝土介質的開坑問題,學者們主要對彈體侵徹混凝土介質的開坑破壞最終形貌進行了深入研究,獲得了豐富的試驗數(shù)據(jù)[3-12],建立了可計算最終開坑直徑、深度、體積等參數(shù)的經驗公式[13-20]、理論模型[21]和數(shù)值模擬方法[22-25],能夠分析彈靶初始條件對開坑破壞最終形貌的影響規(guī)律。而在彈體侵徹混凝土靶的開坑作用過程方面,Liu 等[26]、Yu 等[27]、Feng 等[28]開展了彈體侵徹混凝土靶試驗,獲得了連續(xù)的彈體侵徹混凝土靶的開坑作用過程;Forrestal 等[14]認為彈體在侵徹混凝土靶的開坑作用過程中,靶體阻力線性增長,基于開坑深度等于兩倍彈體直徑假設,建立了開坑阻力線性增長模型,該模型只能得到彈體開坑過程的合阻力,不能給出彈體表面的應力形式;Warren[29]為了能夠在彈靶分離計算方法中反映開坑作用過程,基于開坑深度等于兩倍彈體直徑假設,結合空腔膨脹理論應力表達式和應力逐層增大法建立了開坑阻力數(shù)值計算方法,該方法能夠對開坑阻力進行數(shù)值計算分析;Rosenberg 等[30]以二維數(shù)值模擬計算結果為基礎,建立了開坑階段靶體阻力解析模型,但該模型不能反映彈體的開坑作用過程;Yankelevsky 等[31]認為在彈體侵徹混凝土靶的開坑作用過程中可以忽略材料破碎行為、裂紋擴展和自由面效應的影響,開坑深度等于彈體頭部長度,并結合空腔膨脹理論建立了開坑阻力計算模型,該模型能夠計算開坑過程的靶體阻力;柴傳國等[32]基于高速錄像記錄的彈體速度隨開坑深度的變化關系,采用最小二乘法擬合得到了開坑過程中彈頭表面的應力形式;薛建鋒等[2]根據(jù)開坑過程中阻力做功的等效關系得到了新阻力模型,并分析了彈體侵徹半無限混凝土靶的侵徹性能,但該方法不能反映彈體侵徹開坑的作用過程;王麗梅等[33]忽略開坑破壞過程的影響,基于混凝土靶的最終破壞形貌,結合滑移線理論和沖量定理推導了開坑阻力計算公式,計算結果與試驗在變化趨勢和持續(xù)時間上一致。綜上所述,學者們針對混凝土侵徹初始階段的研究多關注開坑破壞的最終形貌,較少關注動態(tài)開坑的作用過程,在理論分析和工程模型中多采用簡化等效處理方法,缺乏彈體侵徹混凝土動態(tài)開坑作用過程的分析和討論。

    為進一步分析彈體侵徹混凝土靶的動態(tài)開坑作用過程,本文中將基于開坑破壞過程分析對開坑階段進行劃分,分析各階段對彈體侵徹的影響;結合彈體頭部形狀函數(shù)、Maxwell 的Z 模型[34]流線場分布以及法向膨脹理論[35],建立考慮混凝土飛濺過程影響的開坑阻力分析計算模型,研究典型彈靶參數(shù)對動態(tài)開坑作用過程的影響規(guī)律。

    1 考慮混凝土飛濺過程影響的開坑阻力計算模型

    彈體侵徹混凝土靶開坑階段是整個侵徹過程中的初始階段,也是整個侵徹過程中力學分析最困難和被研究最少的階段。本文中基于試驗[26]得到的開坑破壞過程對彈體侵徹混凝土靶開坑階段進行劃分,分析各階段對彈體侵徹的影響;結合彈體頭部形狀函數(shù)、Z 模型[34]流線場分布、法向膨脹理論[35]等,建立考慮混凝土飛濺過程影響的開坑阻力分析計算模型。

    1.1 開坑階段的劃分與基本假設

    Liu 等[26]、Yu 等[27]、Feng 等[28]開展了大量的彈體侵徹混凝土靶試驗研究,并利用高速錄像記錄了彈體侵徹混凝土靶的開坑作用過程。圖1 為Liu 等[26]拍攝的彈體侵徹混凝土靶的開坑作用過程,從圖1(a)~(b)中均能看出,彈體侵入混凝土靶的開坑過程中,彈著點附近存在明顯的混凝土介質粉碎飛濺現(xiàn)象,從100 μs 到300 μs,彈體完全進入靶內,粉碎介質按一定角度從靶面飛濺,粉碎飛濺區(qū)域增大,靶面隆起較小,未見明顯靶面裂紋擴展;300 μs 以后,混凝土靶迎彈面介質破壞形式發(fā)生改變,靶面隆起和裂紋擴展現(xiàn)象逐漸明顯,伴隨有混凝土塊剝落飛出現(xiàn)象。通過高速錄像可以發(fā)現(xiàn),發(fā)生剝落現(xiàn)象的區(qū)域大于粉碎飛濺區(qū)域,剝落飛出的混凝土塊尺寸明顯大于粉碎飛濺的混凝土介質尺寸;混凝土介質粉碎飛濺現(xiàn)象和大塊剝落現(xiàn)象存在明顯的時間先后順序,粉碎飛濺現(xiàn)象發(fā)生在彈體入靶過程中,而混凝土大塊剝落現(xiàn)象發(fā)生在彈體完全入靶后。結合彈體侵徹混凝土靶開坑作用過程中的破壞現(xiàn)象,依據(jù)混凝土介質破壞形式的不同,可將開坑過程分為兩個不同的階段:粉碎飛濺成坑階段(particle ejection)和大塊層裂成坑階段(spallation),如圖2 所示。

    圖1 彈體侵徹混凝土靶動態(tài)開坑過程[26]Fig.1 The dynamic process during the cratering stage of a projectile penetrating into a concrete target[26]

    圖2 彈體侵徹混凝土靶開坑階段劃分Fig.2 Two phases in the cratering stage during projectile penetration into concrete target

    粉碎飛濺成坑階段發(fā)生在彈頭侵入混凝土靶過程中,由混凝土介質破碎流動形成[36],如圖2(a)所示。在由于彈頭高速擠入混凝土介質產生的剪切和高壓作用下,彈著點附近混凝土發(fā)生粉碎性破碎[16],并在彈頭擠壓以及迎彈面拉伸波作用下,粉碎介質從靶體內部向靶面流動,在靶面發(fā)生飛濺,飛濺方向與彈體侵徹方向相反并成一定夾角;隨著彈頭不斷侵入混凝土靶,粉碎性飛濺成坑區(qū)域逐漸增大,當彈頭侵入混凝土靶一定深度后,粉碎飛濺成坑區(qū)域范圍趨于穩(wěn)定,隨后進入大塊層裂階段。

    大塊層裂成坑階段由混凝土靶內部裂紋擴展以及迎彈面附近殘余拉伸波作用形成,如圖2(b)所示;混凝土介質是典型的脆性材料,受彈體高速沖擊的影響,彈靶接觸面附近粉碎區(qū)以外一定范圍內,混凝土介質會經歷較大的剪切變形,形成大量微裂紋,隨著時間的推移,微裂紋逐漸擴展并相互貫通形成宏觀裂紋[16];受自由面效應的影響,迎彈面附近靶體內部宏觀裂紋逐漸向更容易擴展的迎彈面運動,當裂紋擴展到迎彈面時,受靶體內部殘余拉伸波的影響,混凝土介質發(fā)生隆起,伴隨大塊的層裂剝落現(xiàn)象,形成最終的彈坑。此外,對于開坑階段裂紋擴展,Xu 等[25]開展了彈體侵徹混凝土的數(shù)值模擬,分析指出,當彈體已經侵徹大約3 倍彈丸直徑的深度(約135 μs)時,才可以觀察到早期的裂紋擴展,這些裂紋后來發(fā)展成完整彈坑。Yankelevsky 等[31]分析認為,開坑區(qū)裂紋擴展是一個耗時的過程,取決于混凝土材料的低裂縫擴展速度、沿裂縫表面發(fā)展的抗剪能力。結合圖1 所拍攝的彈體入靶過程也可以看出,當彈體完全進入混凝土靶后才出現(xiàn)明顯的層裂剝落現(xiàn)象。因此,與彈體在大塊層裂區(qū)的運動時間尺度相比,大塊層裂區(qū)混凝土介質發(fā)生宏觀裂紋擴展并形成層裂是一個相對緩慢且滯后的過程,在出現(xiàn)宏觀拉伸裂紋擴展時,彈體已經通過該區(qū)域,進入隧道區(qū)[31]。

    結合混凝土靶開坑階段分析以及彈體運動相對位置可知,混凝土粉碎飛濺發(fā)生在彈頭侵入混凝土靶的過程中,這將直接影響彈靶實際接觸面積的大小變化,從而影響彈體侵徹的阻力;而混凝土大塊層裂主要發(fā)生在彈頭通過開坑區(qū)以后,該階段存在明顯的時間滯后效應,對彈體侵徹影響小[25,31]。因此,假設彈頭侵入混凝土靶的開坑過程僅受飛濺效應的影響,彈頭侵入開坑階段靶體的響應區(qū)分別為飛濺區(qū)、壓實區(qū)[35]和未擾動區(qū),飛濺區(qū)與壓實區(qū)均為粉碎區(qū),如圖3 所示。飛濺區(qū)對彈體的侵徹阻力為零,開坑過程中彈體所受阻力由壓實區(qū)提供,忽略自由面效應對壓實區(qū)的影響,壓實區(qū)介質對彈體的阻力可采用描述隧道區(qū)阻力形式的法向膨脹理論對彈靶壓實接觸面積積分近似求解;當彈頭侵入混凝土靶的開坑過程不受飛濺效應的影響時,飛濺區(qū)達到最大,隨后彈頭侵徹的開坑過程僅為簡單的彈頭嵌入壓實區(qū)過程;當彈頭表面與混凝土靶壓實區(qū)完全接觸時,彈體開坑階段結束[31],隨后進入穩(wěn)定侵徹階段。對于開坑飛濺現(xiàn)象,Maxwell[34]于1977 年建立了描述該現(xiàn)象的Z 模型。該模型認為,撞擊發(fā)生后,靶體介質內瞬時形成一簇以撞擊點為中心的流線。在撞擊過程中,靶體內的飛濺介質依次沿固定流線運動,且不同流線上的材料顆粒相互之間不存在干擾,并在靶面沿固定角度飛濺[37]。因此,為了能夠評估開坑過程中飛濺區(qū)的大小,結合彈體頭部形狀函數(shù)、Z 模型流線場分布,對開坑飛濺效應進行分析。為了簡化問題,做出如下幾點假設:

    圖3 彈體侵徹混凝土靶的開坑過程Fig.3 The process of projectile penetration into concrete target

    (1) 彈體開坑作用過程中,彈頭不發(fā)生變形或侵蝕,即為剛性彈體;

    (2) 彈體開坑作用過程中,粉碎區(qū)介質在靶體內部運動時可視為密度為混凝土極限壓縮密度[35]的不可壓縮流體;

    (3) 彈體開坑作用過程中,粉碎區(qū)介質不發(fā)生飛濺現(xiàn)象時,介質沿彈體表面法線方向運動;當粉碎區(qū)介質發(fā)生飛濺現(xiàn)象時,介質改變運動方向,沿著Z 模型的流線運動,彈頭的侵入不改變飛濺介質的運動方向;

    (4) 飛濺效應主要由粉碎介質所受的拉伸波和壓縮波共同決定,拉伸波僅改變粉碎介質的運動方向,不改變介質的動能,且作用時間短,作用速度快。

    1.2 彈體頭部形狀函數(shù)

    彈體正侵徹混凝土靶時,所受的阻力主要來自彈體頭部與混凝土接觸的表面。因此,忽略彈身的影響,基于局部相互作用理論,引入局部笛卡爾坐標系描述尖卵形彈體的頭部形狀,如圖4 所示。圖中,r為彈體半徑,Lm為彈體頭部長度,R為頭部曲率半徑,ψ 為彈頭表面法向角,J為彈體頭部曲率半徑R與彈體直徑2r的比值,即彈體頭部形狀系數(shù)。因此,尖卵形彈體頭部形狀函數(shù)在局部笛卡爾坐標系中可表示為:

    圖4 彈體頭部形狀的結構參數(shù)示意圖Fig.4 Structural parameters of projectile

    考慮到尖卵形彈體為旋轉對稱彈體,為便于分析計算,將彈體頭部形狀函數(shù)轉換至柱坐標系下。借助柱坐標系與笛卡爾坐標系間的轉換關系:

    式(1)可以轉換為:

    式中:ρ 為極半徑, φ 為極角。

    1.3 開坑飛濺流線場基本方程

    基于以上假設,結合Z 模型[34]流線場分布以及彈體頭部形狀函數(shù),建立彈體開坑過程中彈頭表面微元對應的飛濺流線極坐標方程[34]:

    飛濺流線極坐標方程與彈體結構的關系如圖5 所示。其中,h為彈體侵徹深度;L0為撞擊點O(極坐標原點)到彈頭表面微元的直線長度;L為撞擊點O到流線任意位置的極半徑;?0為撞擊點O到彈頭表面微元直線與彈體軸線的夾角;?為極半徑L與彈體軸線的夾角;Lf為彈頭表面微元對應的流線靶面飛濺位置到撞擊點O的直線距離;K為沖擊波衰減相關系數(shù),與靶面飛濺角γe有關,γe為介質飛濺方向與迎彈面的夾角。為了簡化問題分析,Z 模型中假設任意流線介質靶面飛濺角均相同。靶面飛濺角為0°時,K=2;靶面飛濺角為45°時,K=3;靶面飛濺角為65°時,K=4。L0、?0、Lf與彈體結構特性的關系如下:

    圖5 基于極坐標的飛濺流線方程Fig.5 Ejection streamline equation in polar coordinates

    結合流線極坐標方程(4)~(7),得到球坐標系(球坐標中心與極坐標中心O重合,如圖5 所示)中彈體表面微元對應的流管體積方程和質量方程:

    式中:Q1為粉碎區(qū)混凝土的極限壓縮密度, ? φ 為彈體表面微元對應的圓心角。如圖5 所示,彈體表面微元弧長為Δb,彈體表面微元對應的流管積分區(qū)間為[Lf,Lf+ΔLf],ΔLf可近似表示為流管截面寬度[38],Kurosawa 等[38]認為ΔLf=10-4r,由于Δb和ΔLf均為極小的微元量,本文中假設ΔLf=Δb,則ΔLf可以表示為:

    1.4 開坑飛濺流線場能量方程

    建立開坑飛濺流線場是確定飛濺區(qū)大小的基礎,而判斷流管是否發(fā)生飛濺的能量關系是確定飛濺區(qū)大小的關鍵。Kurosaua 等[38]認為,沖擊開坑過程中介質發(fā)生飛濺主要由流管本身的強度能、重力勢能以及輸入流管的動能決定,當輸入流管的動能大于流線強度能和重力勢能時,介質沿流管發(fā)生飛濺。參考Kurosaua 等[38]的工作,忽略重力勢能的影響,結合流管體積方程(8)以及混凝土材料的強度特性,建立彈體表面微元對應流管的強度能方程:

    式中:ε 為混凝土介質完全壓實對應的體應變,σcf為混凝土抗壓強度。

    沖擊開坑過程中,彈靶接觸面附近粉碎區(qū)介質飛濺所需的能量主要由該區(qū)域產生的壓縮波和拉伸波做功提供,壓縮波和拉伸波的作用過程存在明顯的先后順序[39]。彈頭侵徹靶體介質時,彈靶接觸面附近介質發(fā)生粉碎性破壞,壓縮波和粉碎介質均沿著彈體表面微元法線方向運動[31,35];隨后撞擊點附近自由面產生的拉伸波對粉碎介質做功,改變介質運動方向,介質在拉伸波和壓縮波的共同作用下發(fā)生飛濺運動[39]。結合基本假設(4),則輸入流管的動能等于彈體表面微元對應的法向運動介質壓縮動能。法向膨脹理論[35]充分考慮了混凝土的材料行為、空腔膨脹理論和壓縮介質波傳播的影響,能夠計算壓縮波作用后彈體表面微元對應法向運動介質的動能。因此,忽略自由面效應的影響,彈體表面微元對應的法向運動介質動能采用法向膨脹理論進行計算,則輸入彈體表面微元對應流管的動能方程為[35]:

    式中:l為響應介質波相對于彈體表面的傳播距離,w為靶體介質到彈體表面微元的法向距離,w的取值范圍為[0,l],vw為距離彈頭表面微元的法向距離為w處的介質的法向運動速度,其與彈體表面法向速度的關系為[35]:

    式中:vn為彈體表面的法向速度。

    聯(lián)立式(11)~(12),即可得到輸入流線的動能克服強度能后的剩余動能為:

    聯(lián)立式(4)~(14),結合Kurosaua 等[38]建立的飛濺能量關系可知,彈體開坑過程中,當彈頭表面微元對應流線剩余動能小于零時,介質不發(fā)生飛濺;當彈頭表面微元對應流線剩余動能大于零時,介質沿對應流線發(fā)生飛濺;當彈體表面微元對應流線剩余動能等于零時,則該微元對應流線為飛濺區(qū)邊界,結合式(4)可得,彈體開坑過程中飛濺區(qū)的直徑、深度分別為:

    式中:Le為撞擊點到飛濺區(qū)邊界彈體表面微元的直線長度,?e為撞擊點到飛濺區(qū)邊界彈體表面微元的直線與彈體軸線的夾角,De為飛濺區(qū)直徑,He為飛濺區(qū)深度。式(15)~(16)中相關參數(shù)定義如圖6 所示。

    圖6 開坑飛濺區(qū)示意圖Fig.6 Schematic diagram of ejection region

    1.5 開坑阻力計算方法

    開坑飛濺區(qū)的大小直接影響開坑阻力的求解,由式(15)~(16)即可得到彈體開坑過程中飛濺區(qū)范圍變化關系,則彈體開坑過程中彈頭表面微元所受的應力表達式為[35]:

    式中:σd為混凝土的動態(tài)壓縮極限應力,Q0為混凝土材料的初始密度,λ(l)為波陣面相關系數(shù),plock為混凝土材料的鎖定壓力,v0為彈體初始速度。

    彈體開坑過程中,彈頭表面微元應力可分解為垂直于微元法線的法向應力σn和沿微元切線方向的切向應力στ。根據(jù)庫倫摩擦定律,開坑過程中彈頭表面微元摩擦應力為:

    式中:μ為彈靶間的滑動摩擦因數(shù),彈體侵徹混凝土靶時的滑動摩擦因數(shù)具有一定的速度依賴性,為了盡量減小開坑過程求解的復雜性,滑動摩擦因數(shù)可取為常數(shù)[40],本文中取μ=0.01[41]。

    侵徹過程中,彈體表面微元在侵徹方向上的合應力為:

    結合式(21),將彈體表面微元在侵徹方向上的合應力σz對開坑過程中壓實區(qū)的彈靶接觸面積積分,即可得到彈體開坑過程中所受的軸向合阻力:

    式中:?為開坑階段壓實區(qū)彈靶接觸面積,ds為彈體表面微元面積。有:

    基于剛性彈假設,結合式(21),開坑階段彈體在局部坐標系下的運動方程為:

    式中:M為彈體總質量,a為彈體運動加速度。由于式(22)、(24)為非線性微分方程,對于開坑過載、彈體速度以及侵徹深度的計算沒有簡單的解析解,因此開坑過程采用數(shù)值方法迭代求解。彈體開坑過程計算流程如圖7 所示,當彈靶接觸面積剛好等于彈頭表面積時,跳出循環(huán),開坑階段結束。此時彈靶剛好完全接觸,接觸深度為:

    圖7 彈體侵徹混凝土靶體的開坑效應計算流程Fig.7 Calculation process of crater effect of projectile penetrating concrete target

    式中:Hc為彈靶完全接觸深度,Hemax為飛濺區(qū)最大深度。

    開坑結束后,彈體進入彈靶完全接觸的穩(wěn)定侵徹階段,為了保證開坑階段和穩(wěn)定侵徹階段阻力的連續(xù)性,穩(wěn)定侵徹階段阻力可直接采用文獻[35]中的方法計算。

    基于圖7 的流程計算典型尖卵形彈體(r=38.1 mm,J=3,M=13 kg,v0=500 m/s)侵徹半無限混凝土靶(σcf=39 MPa)的動態(tài)開坑過程。如圖8 所示,為了簡化飛濺區(qū)的形狀繪制,圖8 中粉碎飛濺區(qū)形狀曲線是由式(15)~(16)得到的飛濺區(qū)邊界流線對應的起點(彈體表面與飛濺區(qū)邊界相交點)、終點(靶面臨界飛濺位置)以及兩點之間的中間點所建立的樣條曲線,能夠直觀反映開坑過程中飛濺區(qū)的直徑和深度的變化規(guī)律,不代表真實的粉碎飛濺區(qū)形狀。從圖8 可以看出,隨著彈頭的侵入,粉碎飛濺區(qū)逐漸增大,這與Liu 等[26]、Yu 等[27]、Feng 等[28]在試驗中觀察到的現(xiàn)象相同;此外,如圖8 所示,與Forrestal 等[14]建立的開坑阻力線性增長模型、Yankelevsky 等[31]建立的彈頭簡單嵌入模型相比,本文中建立的考慮混凝土飛濺過程的開坑阻力分析模型能夠計算開坑阻力的同時,還能進一步描述彈體侵徹混凝土靶動態(tài)開坑過程中飛濺區(qū)的大小變化過程。

    圖8 彈體侵徹混凝土靶開坑過程的計算結果對比Fig.8 Comparison of calculation results of penetration process during the cratering stage

    2 考慮混凝土飛濺過程影響的開坑阻力計算模型可靠性驗證

    2.1 開坑阻力計算可靠性驗證

    為了驗證考慮混凝土飛濺過程影響的開坑阻力計算模型的正確性,采用文獻中的過載試驗數(shù)據(jù)[14]和彈體入靶過程中速度隨侵徹深度的變化試驗數(shù)據(jù)[32]對理論模型進行驗證。文獻中的彈靶主要參數(shù)如表1 所示。

    表1 彈靶主要參數(shù)Table 1 Parameters of projectiles and targets

    由于上述文獻中均未給出混凝土的極限密度Q1以及混凝土材料的鎖定壓力plock,參考文獻[35],取混凝土極限密度Q1=2 640 kg/m3,混凝土材料的鎖定壓力plock與混凝土強度σcf的比值為16.7,用于上文理論模型的計算;此外,上述文獻也并未給出開坑粉碎飛濺的飛濺角度且不具有可便于測量飛濺角度的高速錄像照片,結合Liu 等[26]、Yu 等[27]、Feng 等[28]的試驗對開坑飛濺過程中的飛濺角進行測量,得到平均飛濺角分別為37.50°、39.12°、47.78°。因此,可以看出,彈體結構對開坑飛濺角度的影響較小,不同試驗飛濺角均接近45°,式(4)中沖擊波衰減相關系數(shù)K可取為3。隨后,將表1 中Forrestal 等[14]的試驗彈靶參數(shù)分別代入本文開坑阻力計算模型、阻力線性增長模型[14]以及忽略開坑飛濺影響的阻力計算模型,得到彈體開坑過程中過載的計算結果。結合文獻中的試驗數(shù)據(jù)對考慮飛濺過程影響的開坑阻力計算模型的可靠性進行驗證,并與其他模型進行對比分析。其中,阻力線性增長模型中取混凝土靶的靜阻力為360 MPa[14,31]。

    彈體侵徹混凝土靶開坑過載計算結果與試驗結果的對比如圖9~10 所示。從圖9~10 中可以看出,考慮飛濺效應影響的開坑阻力模型的過載計算結果與試驗結果均吻合較好;若開坑階段采用阻力線性增長模型[14]計算開坑阻力,彈體頭部形狀系數(shù)為3 時,計算結果與試驗結果相比明顯偏低,而彈體頭部形狀系數(shù)為6 時,計算結果與試驗結果相比存在偏高的現(xiàn)象;此外,從圖9(c)~(d)以及圖10(b)可以看出,若開坑階段忽略飛濺效應的影響,計算結果將過高的估計開坑階段的靶體阻力。

    圖9 彈體頭部形狀系數(shù)為3 時開坑過載計算結果與試驗對比Fig.9 Comparison of acceleration between calculated and test results during the cratering stage of J=3

    圖10 彈體頭部形狀系數(shù)為6 時開坑過載計算結果與試驗對比Fig.10 Comparison of acceleration between calculated and test results during the cratering stage of J=6

    此外,將表1 中柴傳國等[32]的試驗彈靶參數(shù)代入考慮開坑飛濺過程影響的阻力模型中,得到速度隨侵徹深度變化的計算結果。圖11 為彈體入靶過程中速度隨侵徹深度變化的計算結果與試驗結果的對比情況。從圖11 可以看出,模型計算的開坑階段速度隨侵徹深度變化的計算結果與試驗結果吻合較好,進一步說明了模型能夠可靠地預測開坑過程中的靶體阻力。

    圖11 彈體入靶過程中速度衰減的計算結果與試驗結果的對比Fig.11 Comparison of velocity change between calculation results and test results

    2.2 開坑飛濺區(qū)參數(shù)計算可靠性驗證

    為了進一步驗證考慮混凝土飛濺過程影響的開坑阻力計算模型的可靠性,對文獻[26-27]中高速錄像拍攝的開坑過程照片進行處理,由于文獻中高攝錄像照片數(shù)量有限且不能觀測到飛濺區(qū)深度的變化過程,因此主要通過高速錄像獲得了未發(fā)生明顯大塊層裂和靶面隆起前的開坑直徑作為飛濺區(qū)的最大直徑,用于驗證模型計算的有效性,飛濺區(qū)最大直徑的測量方法如圖12 所示。

    忽略文獻[26-27]中混凝土所含纖維等復合材料對開坑過程的影響,將文獻中彈靶參數(shù)分別代入考慮飛濺過程影響的開坑阻力計算模型,得到彈體開坑過程中最大飛濺區(qū)直徑,結合文獻中的試驗數(shù)據(jù)對模型的可靠性進行驗證。計算結果與試驗數(shù)據(jù)對比如表2 所示。從表2 中可以看出,理論計算最大飛濺區(qū)直徑與試驗數(shù)據(jù)吻合較好,考慮飛濺過程影響的開坑阻力計算模型能夠表征彈體侵徹混凝土靶動態(tài)開坑過程中的飛濺區(qū)大小。

    3 動態(tài)開坑作用過程分析

    彈體侵徹混凝土靶動態(tài)開坑過程受多種因素的影響,包括彈體初始撞擊速度、彈體頭部形狀系數(shù)、彈體直徑、混凝土強度等。為了進一步研究彈體侵徹混凝土靶的動態(tài)開坑作用過程,下面討論初始撞擊速度、彈體頭部形狀系數(shù)、混凝土強度對動態(tài)開坑作用過程的影響。

    3.1 初始撞擊速度的影響

    基于已驗證過的理論模型和2.1 節(jié)中Forrestal 等[14]的彈靶參數(shù)(r=38.1 mm,J=3,M=13 kg,σcf=39 MPa)開展撞擊速度分別為200、400、600 和800 m/s 的理論計算,獲得了不同撞擊速度下彈體侵徹混凝土靶的動態(tài)開坑作用過程。撞擊速度為400、800 m/s 的動態(tài)開坑作用過程如圖13(a)所示;圖13(b)~(c)為不同撞擊速度下飛濺區(qū)直徑和深度隨侵徹深度的變化過程,圖13 中t1為飛濺區(qū)直徑和深度最大時刻,t2為彈體侵徹開坑結束時刻,。

    圖13 不同撞擊速度下混凝土靶的動態(tài)開坑作用過程Fig.13 Cratering process of concrete under different impact velocities

    由圖13 可知,初始撞擊速度范圍為200~800 m/s 時,飛濺區(qū)最大直徑約為2r~7r,最大深度約為0.2r~0.7r;隨著撞擊速度的增大,飛濺區(qū)增大,動態(tài)開坑作用時間以及飛濺區(qū)范圍達到穩(wěn)定的時間均縮短,而飛濺區(qū)范圍達到穩(wěn)定的時間占動態(tài)開坑作用時間的比例增大;初始撞擊速度從200 m/s 增大到800 m/s 時,飛濺區(qū)直徑增大了264%,飛濺區(qū)深度增大了238%,而飛濺區(qū)范圍達到穩(wěn)定的時間縮短了54%,動態(tài)開坑作用時間縮短了73%。撞擊速度對彈體侵徹混凝土靶的開坑作用過程影響顯著。

    結合式(1)~(16)可知,彈體頭部結構和混凝土強度相同時,開坑過程中彈頭表面對應的流線強度能相同;而輸入流線的動能主要由彈靶相互作用過程中彈靶接觸面上的壓力做功提供,初始撞擊速度越大,彈靶作用壓力越大,彈體輸入混凝土介質的流線動能越大。因此,隨著初始撞擊速度的增大,流線動能克服強度能后的剩余動能越大,彈體侵徹混凝土靶動態(tài)開坑作用過程中的飛濺區(qū)越大。

    3.2 彈體頭部形狀系數(shù)的影響

    以Forrestal 等[14]的試驗彈體結構(r=38.1 mm,M=13 kg)為基礎,保持彈體直徑、質量相同,設計3 種頭部形狀系數(shù)分別為2、4、6 的尖卵形彈體。隨后基于理論模型,獲得了撞擊速度為600 m/s 時,3 種彈體侵徹39 MPa 混凝土靶的動態(tài)開坑作用過程。圖14(a)為彈體頭部形狀系數(shù)分別為2 和6 時彈體動態(tài)開坑作用過程;圖14(b)~(c)為不同頭部形狀系數(shù)彈體動態(tài)開坑過程中飛濺區(qū)直徑和深度隨侵徹深度的變化過程。

    圖14 不同頭部形狀彈體撞擊混凝土靶的動態(tài)開坑作用過程Fig.14 Cratering process of concrete under impact of projectile with different ogive nose shapes

    改變彈體頭部形狀系數(shù)對彈體頭部結構影響較大。由式(1)~(10)可知,混凝土強度相同時,增大彈體頭部形狀系數(shù),彈頭長細比增大,開坑過程中彈頭表面對應的流線體積增大,從而導致流線強度能增大;另一方面,開坑飛濺區(qū)的擴大主要發(fā)生在彈頭侵入初期,時間小于118 μs 時,如圖14(b)~(c)所示,此時彈頭均未完全進入混凝土靶;彈頭侵入深度相同時,彈體頭部形狀系數(shù)越大,彈頭進入混凝土靶部分的表面積越??;此外,由式(12)~(13)、(17)~(18)可知,撞擊速度相同時,增大彈體頭部形狀系數(shù),彈頭表面法向角增大,彈靶接觸壓力減小,從而使得彈體輸入混凝土靶介質的流線動能減小。

    因此,結合圖14 可知,由于流線強度能的增大,彈體輸入流線的動能減小,流線動能克服強度能后的剩余動能減小,隨著彈體頭部形狀系數(shù)的增大,飛濺區(qū)減小,動態(tài)開坑作用時間以及飛濺區(qū)范圍達到穩(wěn)定的時間逐漸增大,飛濺區(qū)范圍達到穩(wěn)定的時間占動態(tài)開坑作用的時間比例減?。粡楏w頭部形狀系數(shù)從2 增大到6 時,飛濺區(qū)直徑減小了36%,深度減小了28%,而飛濺區(qū)范圍達到穩(wěn)定的時間增加了36%,動態(tài)開坑作用時間增加了60%。彈體頭部長度是影響整個開坑作用時間和飛濺區(qū)范圍達到穩(wěn)定時間的主要因素,隨著彈體頭部形狀系數(shù)的增大,彈頭長度增大,開坑作用時間增大[31],飛濺區(qū)范圍達到穩(wěn)定的時間也增大。與初始撞擊速度對彈體侵徹混凝土靶動態(tài)開坑作用過程的影響相比,彈體頭部形狀系數(shù)的改變對動態(tài)開坑作用過程的影響較小。

    3.3 混凝土強度的影響

    基于Forrestal 等[14]的試驗彈體結構質量參數(shù)(r=38.1 mm,J=3,M=13 kg,v0=600 m/s),開展彈體侵徹單軸抗壓強度分別為20、40、60 和80 MPa 混凝土靶的理論計算。圖15(a)為彈體侵徹單軸抗壓強度為40、80 MPa 混凝土靶的動態(tài)開坑作用過程,圖15(b)~(c)為彈體侵徹不同強度混凝土靶動態(tài)開坑過程中飛濺區(qū)直徑和深度隨侵徹深度的變化過程。

    圖15 彈體撞擊不同強度混凝土靶的動態(tài)開坑作用過程Fig.15 Cratering process of concrete with different compressive strength of target

    彈體撞擊速度相同時,混凝土強度越高,彈體侵徹混凝土靶的粉碎區(qū)越小。由式(12)、(16)可知,混凝土強度增大,開坑作用過程中彈頭表面對應的流線強度能和彈靶接觸面壓力也隨之增大。兩者的增大對開坑飛濺區(qū)大小的影響相反,流線強度能越大,飛濺區(qū)越小,而彈靶接觸面壓力越大,粉碎介質獲得的動能越大,飛濺區(qū)越大。由于開坑粉碎區(qū)的減小,且流線強度能增大效應大于流線輸入動能增大效應,結合圖15(b)~(c)可知,隨著混凝土強度的增大,飛濺區(qū)減小,動態(tài)開坑作用時間以及飛濺區(qū)范圍達到穩(wěn)定的時間縮短,而飛濺區(qū)范圍達到穩(wěn)定的時間占動態(tài)開坑作用時間的比例增大?;炷翉姸葟?0 MPa 增大到80 MPa 時,飛濺區(qū)直徑減小了49%,深度減小了39%,而飛濺區(qū)范圍達到穩(wěn)定的時間縮短了22%,動態(tài)開坑作用時間縮短了6%。因此,混凝土強度改變對于飛濺區(qū)大小以及飛濺區(qū)范圍達到穩(wěn)定的時間影響較大,對動態(tài)開坑作用時間的影響較小。此外,結合3.1~3.2 節(jié)分析可知,撞擊速度對彈體侵徹混凝土動態(tài)開坑作用過程影響最明顯,其次為彈體頭部形狀系數(shù)和混凝土強度。

    4 結 論

    結合彈體頭部形狀函數(shù)、Z 模型流線場分布以及法向膨脹理論,建立了考慮混凝土飛濺過程影響的開坑阻力計算模型,分析了典型彈靶參數(shù)對混凝土靶動態(tài)開坑過程的影響規(guī)律,主要結論如下。

    (1) 侵徹過程中,混凝土靶粉碎飛濺對彈靶接觸面積影響較大,進而影響開坑阻力,大塊層裂主要發(fā)生在彈頭通過開坑區(qū)以后,對彈體開坑過程影響小。

    (2) 與開坑阻力線性增長模型、彈頭簡單嵌入模型相比,考慮混凝土飛濺過程影響的開坑阻力模型在能夠準確計算開坑阻力的同時,還能進一步描述彈體動態(tài)開坑過程中飛濺區(qū)大小的變化過程。

    (3) 彈體開坑飛濺區(qū)范圍隨彈體頭部形狀系數(shù)和混凝土強度的增大而減小;飛濺區(qū)范圍達到穩(wěn)定的時間和動態(tài)開坑作用時間隨初速和混凝土強度的增大而縮短,隨彈體頭部形狀系數(shù)的增大而增大;相較于彈體頭部形狀系數(shù)和混凝土強度,初始撞擊速度對動態(tài)開坑作用過程影響更為顯著。

    猜你喜歡
    彈頭彈體流線
    熱身管下小口徑槍彈彈頭殼材料對其膛內運動的影響規(guī)律
    兵工學報(2022年9期)2022-10-11 01:25:48
    尾錐角對彈體斜侵徹過程中姿態(tài)的影響研究
    橢圓截面彈體斜侵徹金屬靶體彈道研究*
    爆炸與沖擊(2022年2期)2022-03-17 07:28:44
    幾何映射
    導彈突防概率解析與仿真計算結果差異分析
    任意夾角交叉封閉邊界內平面流線計算及應用
    STOPAQ粘彈體技術在管道施工中的應用
    上海煤氣(2018年6期)2018-03-07 01:03:22
    旋轉彈控制系統(tǒng)結構與彈體靜穩(wěn)定特性研究
    大型綜合交通樞紐流線組織設計
    CO2混相驅五點井網(wǎng)流線模擬
    韩国高清视频一区二区三区| 男女国产视频网站| 免费高清在线观看视频在线观看| 国产成人午夜福利电影在线观看| 亚洲人成网站在线播| 99热这里只有是精品在线观看| 一区二区三区免费毛片| 少妇猛男粗大的猛烈进出视频| 女人精品久久久久毛片| 欧美bdsm另类| 国产成人免费观看mmmm| 国产国拍精品亚洲av在线观看| 秋霞在线观看毛片| 美女xxoo啪啪120秒动态图| 亚洲精品乱久久久久久| 国产一区二区在线观看av| 亚洲精品日韩av片在线观看| 欧美精品一区二区免费开放| 2022亚洲国产成人精品| 两个人的视频大全免费| 岛国毛片在线播放| 久久精品国产亚洲av涩爱| 一区二区日韩欧美中文字幕 | 国产极品粉嫩免费观看在线 | 国产视频首页在线观看| 色婷婷av一区二区三区视频| 国产精品无大码| 久久久久久久久久久免费av| 国产成人一区二区在线| 尾随美女入室| 国产精品人妻久久久久久| 黄色欧美视频在线观看| 免费黄网站久久成人精品| 欧美亚洲日本最大视频资源| 51国产日韩欧美| 成年人午夜在线观看视频| 国产老妇伦熟女老妇高清| 91精品伊人久久大香线蕉| 亚洲高清免费不卡视频| 亚洲高清免费不卡视频| 嫩草影院入口| 一二三四中文在线观看免费高清| 久久精品久久精品一区二区三区| 多毛熟女@视频| 多毛熟女@视频| 欧美另类一区| 日本av免费视频播放| 五月伊人婷婷丁香| 亚洲av成人精品一二三区| 26uuu在线亚洲综合色| 欧美精品一区二区免费开放| 91国产中文字幕| 一区在线观看完整版| 超色免费av| 亚洲精品自拍成人| 91久久精品电影网| 亚洲人成77777在线视频| 黑丝袜美女国产一区| 亚洲成色77777| 国产免费现黄频在线看| 国产av精品麻豆| 岛国毛片在线播放| 国产成人精品在线电影| 国产男女超爽视频在线观看| 高清午夜精品一区二区三区| 视频在线观看一区二区三区| 搡女人真爽免费视频火全软件| 777米奇影视久久| 制服人妻中文乱码| 久久99热这里只频精品6学生| 国产片内射在线| 黑人欧美特级aaaaaa片| 91久久精品国产一区二区成人| 久久 成人 亚洲| av国产久精品久网站免费入址| 午夜老司机福利剧场| 一二三四中文在线观看免费高清| 夜夜骑夜夜射夜夜干| 国产免费视频播放在线视频| 日韩电影二区| 亚洲国产精品国产精品| 妹子高潮喷水视频| 人妻系列 视频| 日本黄色片子视频| 亚洲人与动物交配视频| 最近的中文字幕免费完整| 久久精品夜色国产| 欧美亚洲日本最大视频资源| 国产一区有黄有色的免费视频| 中国美白少妇内射xxxbb| 国产亚洲av片在线观看秒播厂| 夜夜看夜夜爽夜夜摸| 美女福利国产在线| 精品一区二区三区视频在线| 亚洲人与动物交配视频| 欧美97在线视频| 国产精品无大码| 久久99蜜桃精品久久| 国产成人精品无人区| 22中文网久久字幕| 91久久精品电影网| 亚洲精品乱码久久久久久按摩| 人妻夜夜爽99麻豆av| 成人毛片a级毛片在线播放| av.在线天堂| 日本免费在线观看一区| 久久久久国产网址| 人人妻人人澡人人爽人人夜夜| 成年美女黄网站色视频大全免费 | 亚洲精品日韩av片在线观看| 国产精品成人在线| 最新的欧美精品一区二区| 欧美亚洲 丝袜 人妻 在线| 亚洲av中文av极速乱| 亚洲精品乱码久久久久久按摩| 国产黄色免费在线视频| 天天躁夜夜躁狠狠久久av| 最新的欧美精品一区二区| 国产一区二区在线观看av| 国产一区二区在线观看av| 亚州av有码| 色哟哟·www| 国产高清有码在线观看视频| 亚洲高清免费不卡视频| 18禁在线无遮挡免费观看视频| 国产精品久久久久久精品电影小说| 日韩亚洲欧美综合| 欧美3d第一页| 纯流量卡能插随身wifi吗| 一级爰片在线观看| 少妇的逼水好多| 在线精品无人区一区二区三| 欧美日本中文国产一区发布| 国产成人精品在线电影| 亚洲第一区二区三区不卡| 欧美精品人与动牲交sv欧美| 午夜av观看不卡| 精品国产一区二区久久| 亚州av有码| 18禁在线播放成人免费| 一区二区三区四区激情视频| 一级毛片我不卡| 最近的中文字幕免费完整| 九九爱精品视频在线观看| 国产日韩欧美亚洲二区| 国产色婷婷99| 亚洲精华国产精华液的使用体验| 18+在线观看网站| 色网站视频免费| 男的添女的下面高潮视频| 久久99蜜桃精品久久| 色吧在线观看| 亚洲av电影在线观看一区二区三区| 91久久精品国产一区二区成人| 欧美精品亚洲一区二区| 亚洲国产成人一精品久久久| 精品人妻一区二区三区麻豆| 99久久精品一区二区三区| 伦精品一区二区三区| 天堂中文最新版在线下载| 亚洲熟女精品中文字幕| 精品久久国产蜜桃| 高清欧美精品videossex| 特大巨黑吊av在线直播| 久久国产精品男人的天堂亚洲 | 在线观看www视频免费| 久久久a久久爽久久v久久| av视频免费观看在线观看| 国产精品.久久久| 这个男人来自地球电影免费观看 | 婷婷色综合大香蕉| 一本—道久久a久久精品蜜桃钙片| 搡老乐熟女国产| 国语对白做爰xxxⅹ性视频网站| 丝袜脚勾引网站| 婷婷色av中文字幕| 午夜福利视频在线观看免费| 午夜福利,免费看| 两个人的视频大全免费| 五月伊人婷婷丁香| 一级a做视频免费观看| 丝袜美足系列| 三级国产精品片| 久久久久久久亚洲中文字幕| 国产精品不卡视频一区二区| 成人亚洲精品一区在线观看| 边亲边吃奶的免费视频| 国产一区二区在线观看日韩| 2018国产大陆天天弄谢| 美女福利国产在线| 秋霞在线观看毛片| 天天影视国产精品| 亚洲在久久综合| 男的添女的下面高潮视频| 亚洲欧洲国产日韩| 国产国拍精品亚洲av在线观看| 少妇猛男粗大的猛烈进出视频| 国产乱来视频区| 日本av手机在线免费观看| 91在线精品国自产拍蜜月| 99久国产av精品国产电影| 婷婷色av中文字幕| 国产日韩欧美在线精品| 国产一区二区三区综合在线观看 | 黄色一级大片看看| 中文天堂在线官网| 亚洲天堂av无毛| 777米奇影视久久| 大片免费播放器 马上看| 蜜桃国产av成人99| a级片在线免费高清观看视频| 欧美日韩一区二区视频在线观看视频在线| 日本与韩国留学比较| 18禁裸乳无遮挡动漫免费视频| 午夜福利视频精品| 在线免费观看不下载黄p国产| 欧美 日韩 精品 国产| 国产成人精品一,二区| 国产一区二区三区综合在线观看 | 国产免费一级a男人的天堂| 久久99蜜桃精品久久| 777米奇影视久久| 能在线免费看毛片的网站| 国产亚洲一区二区精品| 精品久久久噜噜| 大香蕉久久网| 国产极品天堂在线| 我的女老师完整版在线观看| 制服人妻中文乱码| 热99国产精品久久久久久7| 国产深夜福利视频在线观看| 国产成人一区二区在线| 成人毛片60女人毛片免费| 日韩中字成人| 97超视频在线观看视频| 99九九在线精品视频| 欧美3d第一页| 免费看不卡的av| 制服丝袜香蕉在线| 九九久久精品国产亚洲av麻豆| 9色porny在线观看| 国产 精品1| 在线观看一区二区三区激情| 五月开心婷婷网| 国产成人精品久久久久久| 精品久久久久久电影网| 日本爱情动作片www.在线观看| 久久人人爽人人片av| 精品一品国产午夜福利视频| 最近中文字幕2019免费版| 日韩精品有码人妻一区| 黄片播放在线免费| 色婷婷久久久亚洲欧美| 十八禁高潮呻吟视频| 久久精品人人爽人人爽视色| 观看av在线不卡| 综合色丁香网| 99热6这里只有精品| 日韩电影二区| 亚洲欧美一区二区三区黑人 | 啦啦啦啦在线视频资源| 精品久久久噜噜| 久久人妻熟女aⅴ| 精品久久蜜臀av无| 丝袜美足系列| 中文字幕亚洲精品专区| 精品人妻熟女毛片av久久网站| 久久99精品国语久久久| 色哟哟·www| 免费日韩欧美在线观看| 新久久久久国产一级毛片| 女性生殖器流出的白浆| av免费观看日本| 亚洲av.av天堂| 免费观看无遮挡的男女| 晚上一个人看的免费电影| 亚洲精品美女久久av网站| 一区在线观看完整版| 午夜免费鲁丝| 欧美亚洲日本最大视频资源| 久久亚洲国产成人精品v| 男女啪啪激烈高潮av片| 大码成人一级视频| 一级毛片我不卡| 日本91视频免费播放| 美女脱内裤让男人舔精品视频| 中国三级夫妇交换| 99久久精品国产国产毛片| 在线观看一区二区三区激情| 久久精品人人爽人人爽视色| 视频区图区小说| 国产成人免费无遮挡视频| 18禁观看日本| 夫妻性生交免费视频一级片| 国产欧美日韩综合在线一区二区| 久久国产亚洲av麻豆专区| 九九爱精品视频在线观看| 免费av不卡在线播放| 看非洲黑人一级黄片| 国产成人免费无遮挡视频| 九九爱精品视频在线观看| 国产精品.久久久| 亚洲国产av影院在线观看| 99九九在线精品视频| 国产日韩欧美亚洲二区| 久久久久久久久久人人人人人人| 极品人妻少妇av视频| 热re99久久国产66热| 99久久人妻综合| 黑丝袜美女国产一区| 色吧在线观看| 一本久久精品| 亚洲国产毛片av蜜桃av| 永久网站在线| 国产亚洲午夜精品一区二区久久| 人妻少妇偷人精品九色| 免费看不卡的av| 中文字幕精品免费在线观看视频 | 丝袜脚勾引网站| 日韩电影二区| 亚洲精品视频女| 免费av中文字幕在线| 国产日韩欧美视频二区| av专区在线播放| 男女啪啪激烈高潮av片| 搡老乐熟女国产| 久久99热6这里只有精品| 女性生殖器流出的白浆| av免费在线看不卡| 人妻人人澡人人爽人人| 又黄又爽又刺激的免费视频.| 我的老师免费观看完整版| 男女边吃奶边做爰视频| 老司机亚洲免费影院| 亚洲av成人精品一区久久| 亚洲成人av在线免费| 99久久中文字幕三级久久日本| 久久久久网色| 高清视频免费观看一区二区| 日本色播在线视频| 亚洲欧美日韩卡通动漫| 香蕉精品网在线| 国产一区二区三区综合在线观看 | 国产 一区精品| av不卡在线播放| 亚洲精品一二三| 天天躁夜夜躁狠狠久久av| 国产成人精品婷婷| 狂野欧美激情性bbbbbb| 考比视频在线观看| 欧美亚洲 丝袜 人妻 在线| 久久久久人妻精品一区果冻| 国产一级毛片在线| 国产毛片在线视频| 国产一区二区在线观看日韩| av在线播放精品| 成年av动漫网址| 亚洲av二区三区四区| 国产精品国产三级国产av玫瑰| 一本久久精品| 国产免费又黄又爽又色| 国产永久视频网站| 丰满迷人的少妇在线观看| 久久久久久久大尺度免费视频| 伦理电影大哥的女人| 国产极品粉嫩免费观看在线 | 久久精品国产亚洲av天美| 国产熟女欧美一区二区| 免费久久久久久久精品成人欧美视频 | 99热这里只有是精品在线观看| 中国美白少妇内射xxxbb| 婷婷色综合大香蕉| 王馨瑶露胸无遮挡在线观看| 久久久国产欧美日韩av| 91午夜精品亚洲一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 午夜日本视频在线| 中文字幕最新亚洲高清| 婷婷色av中文字幕| .国产精品久久| 国产亚洲欧美精品永久| 99热全是精品| 少妇高潮的动态图| 国产男人的电影天堂91| 欧美xxxx性猛交bbbb| 91精品一卡2卡3卡4卡| 亚洲伊人久久精品综合| 亚洲欧美清纯卡通| av专区在线播放| 观看av在线不卡| 高清欧美精品videossex| 中文字幕久久专区| 街头女战士在线观看网站| 国产成人免费无遮挡视频| 一级爰片在线观看| 这个男人来自地球电影免费观看 | 搡老乐熟女国产| 亚洲精品国产av蜜桃| 国产免费福利视频在线观看| 久久久久人妻精品一区果冻| 99热这里只有精品一区| 中文字幕人妻丝袜制服| 一级a做视频免费观看| a级毛色黄片| 看十八女毛片水多多多| 国产亚洲欧美精品永久| 超碰97精品在线观看| 国产av精品麻豆| 青春草亚洲视频在线观看| 婷婷成人精品国产| .国产精品久久| 精品一区在线观看国产| 18禁观看日本| 国产成人精品久久久久久| 亚洲av电影在线观看一区二区三区| 丝瓜视频免费看黄片| 美女大奶头黄色视频| 日日摸夜夜添夜夜爱| av在线app专区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日本免费在线观看一区| 五月伊人婷婷丁香| 亚洲国产色片| 精品亚洲成国产av| 久久国产精品大桥未久av| 97在线视频观看| 国产欧美日韩一区二区三区在线 | 久久 成人 亚洲| 天天影视国产精品| 国产黄色免费在线视频| 男女高潮啪啪啪动态图| 午夜精品国产一区二区电影| 免费观看在线日韩| 狂野欧美激情性xxxx在线观看| 亚洲人成网站在线播| 纵有疾风起免费观看全集完整版| 亚洲精品久久午夜乱码| 亚洲三级黄色毛片| 免费人妻精品一区二区三区视频| 毛片一级片免费看久久久久| 国产色爽女视频免费观看| 婷婷色麻豆天堂久久| 卡戴珊不雅视频在线播放| 久久久久久久久久久丰满| 欧美精品人与动牲交sv欧美| 日日爽夜夜爽网站| 欧美97在线视频| 最黄视频免费看| 久久久久网色| 少妇高潮的动态图| 一级毛片电影观看| 99热这里只有精品一区| 伊人亚洲综合成人网| 卡戴珊不雅视频在线播放| 日韩一区二区视频免费看| 国产欧美日韩综合在线一区二区| 曰老女人黄片| 国产成人精品福利久久| 日韩三级伦理在线观看| a级毛色黄片| 狠狠婷婷综合久久久久久88av| 亚洲国产日韩一区二区| 亚洲国产毛片av蜜桃av| 久久99精品国语久久久| 成人无遮挡网站| 日日摸夜夜添夜夜添av毛片| 十分钟在线观看高清视频www| 人妻系列 视频| 欧美日韩综合久久久久久| 免费av中文字幕在线| 涩涩av久久男人的天堂| 亚洲av国产av综合av卡| 久久97久久精品| 午夜免费男女啪啪视频观看| 人人妻人人澡人人看| 亚洲国产精品成人久久小说| 成人免费观看视频高清| 中国国产av一级| 欧美97在线视频| 26uuu在线亚洲综合色| 久久久久网色| 日韩熟女老妇一区二区性免费视频| 亚州av有码| 香蕉精品网在线| 欧美精品一区二区大全| 91精品国产国语对白视频| 亚洲丝袜综合中文字幕| 看非洲黑人一级黄片| 日韩一区二区三区影片| 国产熟女欧美一区二区| 久久久久人妻精品一区果冻| 久久久国产欧美日韩av| av免费观看日本| 狂野欧美激情性bbbbbb| 欧美xxⅹ黑人| 精品卡一卡二卡四卡免费| 国产毛片在线视频| 色吧在线观看| 人妻系列 视频| 岛国毛片在线播放| 国产黄频视频在线观看| 毛片一级片免费看久久久久| 欧美97在线视频| 视频在线观看一区二区三区| 搡女人真爽免费视频火全软件| 黄片播放在线免费| 久久99热6这里只有精品| 国产又色又爽无遮挡免| 你懂的网址亚洲精品在线观看| 成人国产麻豆网| 国产亚洲精品久久久com| 五月开心婷婷网| 欧美另类一区| 一级爰片在线观看| 国产一区二区在线观看av| 婷婷色av中文字幕| av在线观看视频网站免费| 亚洲国产精品成人久久小说| 熟女电影av网| 国产精品久久久久久av不卡| 亚洲久久久国产精品| 亚洲精品国产色婷婷电影| 免费少妇av软件| 精品久久久久久电影网| 中文字幕免费在线视频6| 视频在线观看一区二区三区| 亚洲四区av| 亚洲人成77777在线视频| 亚洲,一卡二卡三卡| 亚洲,欧美,日韩| 国产片内射在线| 亚洲三级黄色毛片| 免费不卡的大黄色大毛片视频在线观看| 久热久热在线精品观看| 日韩一本色道免费dvd| 亚洲国产色片| 我的老师免费观看完整版| 国产精品嫩草影院av在线观看| 美女脱内裤让男人舔精品视频| 国产成人免费无遮挡视频| kizo精华| 一边亲一边摸免费视频| 欧美 日韩 精品 国产| 国产精品久久久久久久电影| 最后的刺客免费高清国语| 国产精品.久久久| 久久 成人 亚洲| 久久久欧美国产精品| 午夜久久久在线观看| 少妇精品久久久久久久| 草草在线视频免费看| 国产欧美日韩一区二区三区在线 | 国产精品99久久99久久久不卡 | 美女中出高潮动态图| 国产不卡av网站在线观看| 欧美bdsm另类| 国产老妇伦熟女老妇高清| 日韩欧美一区视频在线观看| 日本午夜av视频| 老司机亚洲免费影院| 欧美激情极品国产一区二区三区 | 免费观看在线日韩| 亚洲情色 制服丝袜| 久久久a久久爽久久v久久| 超色免费av| 99久久精品国产国产毛片| 国产精品99久久久久久久久| 国产深夜福利视频在线观看| 成年美女黄网站色视频大全免费 | av线在线观看网站| 男女啪啪激烈高潮av片| 精品卡一卡二卡四卡免费| 亚洲人成网站在线观看播放| 国产一区二区在线观看av| 麻豆乱淫一区二区| 人人妻人人澡人人爽人人夜夜| 久久久精品94久久精品| 97在线人人人人妻| 欧美人与性动交α欧美精品济南到 | 亚洲色图 男人天堂 中文字幕 | 色婷婷av一区二区三区视频| 国产无遮挡羞羞视频在线观看| 这个男人来自地球电影免费观看 | 卡戴珊不雅视频在线播放| 国产精品人妻久久久久久| 久久影院123| 男人添女人高潮全过程视频| 十分钟在线观看高清视频www| 99热这里只有精品一区| 久久久久久久久大av| 日本91视频免费播放| 亚洲av.av天堂| 亚洲av在线观看美女高潮| 久久99蜜桃精品久久| 少妇人妻精品综合一区二区| 91国产中文字幕| 中文字幕制服av| 十分钟在线观看高清视频www| 精品久久久噜噜| 国产亚洲欧美精品永久| 亚洲精华国产精华液的使用体验| 天堂8中文在线网| 免费高清在线观看日韩| 在线 av 中文字幕| 黑丝袜美女国产一区| 最近手机中文字幕大全| 插逼视频在线观看| av电影中文网址| 99视频精品全部免费 在线| 国模一区二区三区四区视频| 看十八女毛片水多多多| 国产av国产精品国产| 九草在线视频观看| 一级毛片我不卡| 亚洲第一区二区三区不卡|