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

    在橢圓橫截面彈體正侵徹下有限厚鋁靶的破壞模式及響應(yīng)特性*

    2022-12-21 08:30:56劉均偉張先鋒趙瑤瑤魏海洋李鵬程
    爆炸與沖擊 2022年12期
    關(guān)鍵詞:靶體長軸彈體

    劉均偉,張先鋒,趙瑤瑤,魏海洋,3,劉 闖,李鵬程

    (1. 南京理工大學(xué)機(jī)械工程學(xué)院,江蘇 南京 210094;2. 江蘇永豐機(jī)械有限責(zé)任公司,江蘇 淮安 211722;3. 北京航天長征飛行器研究所,北京 100074)

    目前,相較于圓形橫截面彈體,各種非圓橫截面彈體由于具有在高超聲速武器平臺(tái)上適應(yīng)性好、裝填比大、抗彎強(qiáng)度高及平臺(tái)內(nèi)部有效載荷空間利用率高等優(yōu)點(diǎn)已成為研究熱點(diǎn)[1-3]。諸多學(xué)者在異形彈體對目標(biāo)的侵徹與貫穿作用機(jī)理方面已開展了大量研究工作。Ben-Dor 等[4-5]、Yakunina[6-7]通過研究發(fā)現(xiàn),對于彈體在稠密介質(zhì)中的侵徹問題,非圓橫截面彈體(多圓錐體切面)為最優(yōu)的彈體結(jié)構(gòu)形狀,能夠有效降低靶體侵徹阻力。Woo[8]基于傅里葉級數(shù)展開和最小二乘法,研究了以壓力/速度邊界條件控制的任意截面形狀的空腔膨脹時(shí)空腔邊界的應(yīng)力分布特性,通過比較發(fā)現(xiàn)非圓橫截面彈體具有更優(yōu)異的侵徹性能。Bless[9]針對矩形橫截面彈體和圓形橫截面彈體進(jìn)行了高速侵徹實(shí)驗(yàn),分析了這2 種彈體的侵徹效率,發(fā)現(xiàn)非對稱橫截面彈體在臨界速度以下的侵徹性能更優(yōu)異。杜忠華等[10-11]開展了圓形、矩形和三角形3 類橫截面彈體侵徹半無限金屬板的實(shí)驗(yàn)和理論研究,發(fā)現(xiàn)在彈體最大橫截面面積和初始速度相同的情況下,三角形橫截面彈體的侵徹效率最高,矩形橫截面彈體的侵徹效率次之,而圓形橫截面彈體的侵徹效率最低。高光發(fā)等[12]基于杜忠華等[10-11]的實(shí)驗(yàn)結(jié)果,開展了相應(yīng)的數(shù)值模擬,進(jìn)一步驗(yàn)證了杜忠華等[10-11]的發(fā)現(xiàn)。Dong 等[13-14]、Gao 等[15]、高旭東等[16]開展了一系列橢圓橫截面彈體侵徹半無限厚混凝土靶的實(shí)驗(yàn),彈體橫截面長短軸長度比分別為1.5 和2.0,實(shí)驗(yàn)數(shù)據(jù)表明,橢圓橫截面彈體在700~1 000 m/s 的撞擊速度范圍內(nèi)具有良好的侵徹性能和彈道穩(wěn)定性;同時(shí),通過數(shù)值模擬研究了橢圓橫截面彈體侵徹混凝土靶的受力特性,結(jié)果表明,在相同侵徹速度下,彈體頭部的法向應(yīng)力從彈體橫截面長軸方向向短軸方向逐漸升高。在有限厚靶板方面,王浩等[17-18]開展了橢圓橫截面彈體正貫穿加筋板的研究,獲得了橢圓橫截面彈體剩余速度和靶板彈道極限速度的預(yù)測公式,發(fā)現(xiàn)隨著橢圓彈體橫截面長短軸長度比的增大,靶板的彈道極限速度近似線性升高。Landkof 等[19]基于剛塑性模型假設(shè)和能量守恒原理,考慮了裂紋傳播、花瓣彎曲和靶板整體凹陷變形能,建立了錐形彈體正侵徹薄鋁板彈道極限速度計(jì)算模型。張中國等[20]和Chen 等[21]開展了加筋板貫穿實(shí)驗(yàn)研究,討論了加強(qiáng)筋對靶板破壞模式和彈道參數(shù)的影響,分析了靶板吸收能量的情況,根據(jù)實(shí)驗(yàn)數(shù)據(jù)修正了De Marre 貫穿經(jīng)驗(yàn)公式。Song 等[22]和徐雙喜等[23]同樣基于剛塑性模型假設(shè)和能量守恒原理,分別對截卵形彈體和截錐形彈體貫穿加筋鋼板進(jìn)行了理論分析。

    彈靶作用過程是一個(gè)復(fù)雜的動(dòng)力學(xué)問題,在異形橫截面彈體對目標(biāo)的侵徹作用過程研究中,對彈體阻力以及侵徹剩余速度和侵徹深度的研究較廣泛,而關(guān)于彈體侵徹作用下靶體響應(yīng)特性的研究工作報(bào)道較少。研究異形橫截面彈體侵徹作用下靶體破壞模式及響應(yīng)特性對理解靶體非對稱破壞、變形和耗能特性以及建立靶體阻力函數(shù)等具有重要意義。因此,本文中設(shè)計(jì)并開展3 種不同橫截面彈體正侵徹2A12 鋁板的實(shí)驗(yàn)研究,獲得橢圓橫截面彈體對2A12 鋁板的破壞特性及彈體剩余速度;在此基礎(chǔ)上,建立橢圓橫截面彈體正侵徹2A12 鋁板的數(shù)值模擬模型,結(jié)合實(shí)驗(yàn)結(jié)果驗(yàn)證模型的有效性,并從工程應(yīng)用的角度,研究不同橫截面彈體正侵徹鋁板的剩余速度;另外,系統(tǒng)分析彈體橫截面長短軸長度比對靶體破壞情況及響應(yīng)特性的影響,包括靶體兩側(cè)花瓣的大小、形狀和分布規(guī)律,靶體塑性應(yīng)變區(qū)域的范圍以及在彈體頭部貫入階段靶體沿周向各點(diǎn)的徑向位移、徑向應(yīng)力和切向應(yīng)力的變化規(guī)律等。

    1 實(shí) 驗(yàn)

    1.1 彈體和靶體

    為了研究不同橫截面形狀彈體正侵徹鋁板的破壞特性,開展了橢圓橫截面彈體正侵徹鋁板的實(shí)驗(yàn)。彈體材料為30CrMnSiNi2A 高強(qiáng)度鋼,經(jīng)淬火處理后洛氏硬度HRC 為55 ,屈服強(qiáng)度約為1 600 MPa。實(shí)驗(yàn)彈體如圖1 所示,包括3 種彈型:C1 型彈、T1 型彈和T2 型彈。彈體主要參數(shù)見表1,a和b分別為彈體橫截面長、短半軸長,β 為彈體橫截面的長短軸長度比,L為彈體頭部的長度,ψ 為彈體頭部曲徑比,m為彈體的質(zhì)量。橢圓橫截面彈體的頭部曲徑比不是一個(gè)常數(shù),它與方位角密切相關(guān),即沿圓周方向由長軸向短軸逐漸增大[2],本文中所用曲徑比均為短軸方向的曲徑比。

    圖1 實(shí)驗(yàn)彈體Fig. 1 Projectiles used in the experiments

    表1 三種彈體主要參數(shù)Table 1 Main parameters of three projectiles

    為保證侵徹彈體的完整性以研究彈體橫截面形狀對靶體破壞效應(yīng)的影響,選用強(qiáng)度較低(屈服強(qiáng)度約為200 MPa)的2A12 鋁合金作為靶體材料。靶體尺寸設(shè)為400 mm×400 mm×30 mm 以忽略靶體兩側(cè)自由面對靶體響應(yīng)的影響,如圖2 所示。實(shí)驗(yàn)現(xiàn)場布置如圖3 所示,采用高速攝像機(jī)觀察彈體飛行姿態(tài),同時(shí)將2 個(gè)平面鏡分別放置于靶前和靶后彈體運(yùn)動(dòng)軌跡的正下方,與水平面呈45°夾角,通過測量和分析彈體與平面鏡中彈體的位置變化,得到彈體的速度、運(yùn)動(dòng)軌跡和姿態(tài)。

    圖2 實(shí)驗(yàn)靶體Fig. 2 Targets used in the experiments

    圖3 實(shí)驗(yàn)布局Fig. 3 Experimental layout

    1.2 實(shí)驗(yàn)結(jié)果

    共開展了10 發(fā)實(shí)驗(yàn),其中C1 型彈體4 發(fā),T1 和T2 型彈體各3 發(fā),通過調(diào)整發(fā)射藥量控制彈體速度在200~600 m/s。C1-4 彈體飛行姿態(tài)如圖4 所示,根據(jù)高速攝影以及平面鏡成像可以核定彈體的著靶姿態(tài)、初始撞擊速度和剩余速度,實(shí)驗(yàn)結(jié)果見表2,v0為彈體初始撞擊速度,vr為彈體剩余速度, α 為彈體的著靶角, γ 為彈體的偏航角。從表2 中可以看出:3 種彈體的初始撞擊速度基本相同,具有一定的可比性;彈體最大著靶角 α 為2.69°,最大偏航角 γ 為-2.39°,基本滿足正侵徹條件。

    表2 彈體正侵徹鋁靶的實(shí)驗(yàn)結(jié)果Table 2 Experimental results for normal penetration of projectiles into aluminum targets

    圖4 C1-4 彈體飛行姿態(tài)分析Fig. 4 Analysis of the flight attitude of the C1-4 projectile

    C1-4 彈體侵徹鋁靶前后的對比如圖5 所示,彈體長度和形狀無明顯變化,僅表面有一定的摩擦痕跡,因此,彈體正侵徹鋁靶的過程可視為剛體侵徹過程。

    圖5 實(shí)驗(yàn)前后C1-4 彈體對比Fig. 5 Comparison of the C1-4 projectile before and after the experiment

    圖6 為在3 種不同橫截面形狀彈體侵徹下鋁靶的破壞形貌。可知,C1 型彈體貫穿鋁靶后,在鋁靶正面形成近似圓形的彈孔,孔徑略大于彈體的直徑,鋁靶背面形成多個(gè)花瓣,且花瓣大小、分布較均勻;T1 和T2 型彈體與C1 型彈體類似,在鋁靶正面形成近似彈體最大橫截面形狀和尺寸的孔洞,而鋁靶背面的花瓣大小和分布與C1 型彈體侵徹下的有較大區(qū)別,靶體破壞情況呈現(xiàn)不對稱性,表現(xiàn)為靶體孔洞橫截面短軸方向的花瓣大于長軸方向的,部分靶體孔洞橫截面長軸方向的花瓣甚至發(fā)生了崩落。這與文獻(xiàn)[24]的結(jié)論一致,由于橢圓橫截面空腔邊界位移的不均勻性,空腔邊界應(yīng)力分布呈現(xiàn)從長軸至短軸逐漸降低的趨勢,因此長軸方向受力更高。

    圖6 不同橫截面形狀彈體侵徹下鋁靶的破壞形貌Fig. 6 Damage of aluminum targets under normal penetration of the projectiles with different cross-section shapes

    2 數(shù)值模擬

    2.1 有限元模型及材料參數(shù)

    采用LS-DYNA 軟件對鋁靶在不同橫截面形狀的彈體正侵徹下的破壞特性進(jìn)行數(shù)值模擬分析。鋁板和彈體均采用八節(jié)點(diǎn)六面體單元,單元算法為拉格朗日算法。為了節(jié)約CPU 資源、提高數(shù)值模擬計(jì)算效率,整個(gè)模型采用1/4 建模,鋁靶迎彈面的中心區(qū)域網(wǎng)格劃分密集,網(wǎng)格大小為0.25 mm×0.25 mm×0.3 mm,網(wǎng)格從中心沿四周采用3 次蝴蝶網(wǎng)格過渡,中心區(qū)網(wǎng)格密集區(qū)域?yàn)?0 倍彈徑,彈靶有限元模型如圖7 所示。通過關(guān)鍵字BOUNDARY_SPC_SET 設(shè)置對稱面約束,彈體和靶體之間通過關(guān)鍵字CONTACT_ERODING_NODES_TO_SURFACE 設(shè)置侵蝕接觸,在靶板兩側(cè)邊緣面設(shè)置無反射邊界條件,并施加固定約束,保證數(shù)值模擬的初始條件與實(shí)驗(yàn)中初始條件一致。

    圖7 有限元模型Fig. 7 Finite element models

    從實(shí)驗(yàn)回收彈體可以看出彈體正侵徹薄鋁靶的過程可視為剛體侵徹,因此對彈體采用MAT_RIGID剛體模型,而對鋁靶則采用高應(yīng)變率下適用的JOHNSON_COOK (JC)材料模型和GRüNEISEN 狀態(tài)方程共同表征,具體材料參數(shù)如表3 所示。表3 中A、B、C、m和n分別為靜態(tài)屈服強(qiáng)度、硬化系數(shù)、硬化指數(shù)、溫度軟化指數(shù)和應(yīng)變率系數(shù), ρ 、G和μ分別為材料的密度、剪切強(qiáng)度和泊松比。

    表3 材料參數(shù)Table 3 Material parameters

    2.2 數(shù)值模型準(zhǔn)確性驗(yàn)證

    根據(jù)實(shí)驗(yàn)初始條件,對不同橫截面彈體正侵徹鋁靶過程進(jìn)行了模擬。從表4 可以看到,數(shù)值模擬的彈體剩余速度與實(shí)驗(yàn)結(jié)果吻合較好,相對誤差?r最大為9.2%。從圖8 可以發(fā)現(xiàn),數(shù)值模擬得到的3 種橫截面彈體正侵徹鋁靶的破壞形貌與實(shí)驗(yàn)結(jié)果吻合良好,僅部分鋁靶孔洞邊緣處一些尺寸較小的花瓣與實(shí)驗(yàn)結(jié)果有一定差距,這可能是由于單元侵蝕算法將部分畸變嚴(yán)重網(wǎng)格刪除所致。數(shù)值模擬結(jié)果中靶板破壞情況與實(shí)驗(yàn)結(jié)果基本吻合,表明本文中使用的網(wǎng)格尺寸、材料模型及參數(shù)等是合理的,能較好地模擬不同橫截面彈體正侵徹下鋁靶的結(jié)構(gòu)響應(yīng)和破壞模式,后續(xù)可依據(jù)此模型對不同橫截面彈體正侵徹下鋁靶的力學(xué)行為進(jìn)行分析和預(yù)測。

    圖8 實(shí)驗(yàn)靶板破壞形態(tài)與數(shù)值模擬結(jié)果對比Fig. 8 Comparison of failure morphologies of targets between simulation and experiment

    表4 彈體剩余速度模擬結(jié)果與實(shí)驗(yàn)結(jié)果的對比Table 4 Comparison of residual velocities of projectiles between simulation and experiment

    3 橢圓橫截面彈體正侵徹有限厚鋁靶作用特性數(shù)值模擬分析

    基于已驗(yàn)證的數(shù)值模型,系統(tǒng)研究了彈體橢圓橫截面長短軸長度比對靶體破壞情況以及響應(yīng)特性的影響,包括靶體兩側(cè)的花瓣大小、形狀和分布規(guī)律,以及靶體塑性應(yīng)變區(qū)域的范圍等,并從工程應(yīng)用的角度,計(jì)算了不同橫截面形狀彈體正侵徹鋁靶后的剩余速度;同時(shí)在極坐標(biāo)系中研究了在彈體頭部貫入階段靶體沿周向各點(diǎn)的徑向位移、徑向應(yīng)力和切向應(yīng)力的變化規(guī)律。

    3.1 不同橫截面形狀彈體正侵徹鋁靶后的剩余速度

    為研究彈體橢圓橫截面長短軸長度比對彈體剩余速度的影響,根據(jù)C1 和T1 彈體,在保證彈體橫截面面積相同的情況下設(shè)計(jì)了2 種橢圓橫截面彈體(T3 和T4 彈體),這4 種彈體橫截面長短軸長度比β 分別為1.00、1.61、1.25 和2.00,彈體的初始速度v0為200~600 m/s,計(jì)算相應(yīng)的剩余速度vr,如圖9 所示,S為該彈體的最大橫截面面積。通過調(diào)整彈身部分的開孔尺寸,保證各彈體的質(zhì)量和頭部長度一致。

    從圖9 可以看出,C1、T1、T3 和T4 型彈體的剩余速度差異較小,這表明橢圓彈體橫截面長短軸長度比對其剩余速度的影響較弱。因此,在工程應(yīng)用中,可直接使用圓形橫截面彈體正貫穿靶板的工程模型來計(jì)算具有相同最大橫截面面積的橢圓橫截面彈體的剩余速度。T2 和T3 型彈體雖然具有相同的橫截面長短軸長度比,但T2 型彈體的最大橫截面面積大于T3 型彈體的,從兩者的剩余速度曲線可以看出,T2 型彈體的剩余速度低于T3 型彈體的,這意味著彈體最大橫截面面積對其正侵徹靶板后的剩余速度影響較顯著。在設(shè)計(jì)彈體結(jié)構(gòu)時(shí),在保證彈體強(qiáng)度和剛度的情況下,盡量減小彈體的最大橫截面面積以提高彈體穿靶后的剩余速度。

    圖9 不同橫截面形狀彈體剩余速度對比Fig. 9 Comparison of residual velocities of projectiles with different cross-section shapes

    3.2 彈體橫截面長短軸長度比對靶體破壞情況及響應(yīng)特性的影響

    為分析橢圓橫截面彈體長短軸長度比對靶體破壞情況的影響,以長短軸長度比β 分別為1.00、1.25、1.61 和2.00 的4 種彈體(質(zhì)量、長度、最大橫截面面積等均一致)為計(jì)算模型,分析在400 m/s 彈體撞擊速度下靶體破壞情況及響應(yīng)特性的差異,如圖10~11 所示。

    從圖10 可以看出,圓形與橢圓橫截面彈體侵徹作用下靶體的破壞程度有一定差異(圖中長軸方向?yàn)殚L軸平視圖,即彈道邊緣對應(yīng)彈體長軸)。從靶體背部的花瓣數(shù)量來看,β=1.00 的C1 型彈體侵徹造成的花瓣數(shù)量最少(6 個(gè)),而β=2.00 的T4 型彈體侵徹造成的花瓣數(shù)量最多(多于20 個(gè))。從花瓣的形態(tài)分析,前者靶體背部的花瓣大小一致、形狀規(guī)則,且分布均勻,而后者靶體背部的花瓣則大小不一、形狀不規(guī)則,且分布角度差距較大??偟膩碚f,彈體橫截面長短軸長度比越大,侵徹后靶體背部的花瓣數(shù)量越多,尺寸越小,且短軸方向的花瓣數(shù)量和靶體表面隆起高度大于長軸方向的。此外,從塑性應(yīng)變分布可以看出,在花瓣處的應(yīng)變值明顯高于其余位置的,且長軸方向的高塑性應(yīng)變區(qū)域更集中,意味著靶體在長軸方向受力更高,容易產(chǎn)生破壞,這也解釋了實(shí)驗(yàn)中所觀察到長軸方向花瓣發(fā)生的崩落現(xiàn)象。

    圖10 彈體橫截面長短軸長度比對靶體破壞形貌的影響Fig. 10 Influence of major-to-minor axis length ratios of the projectile cross-sections on damage morphologies of the targets

    從圖11 可以看出:圓形橫截面彈體侵徹作用下靶體表面、背部塑性應(yīng)變區(qū)域?yàn)閳A形;其余3 種橢圓橫截面彈體侵徹作用下靶體表面、背部塑性應(yīng)變區(qū)域呈現(xiàn)類橢圓形,同時(shí)長軸方向的塑性應(yīng)變影響范圍明顯小于短軸方向的。

    圖11 靶體應(yīng)變分布Fig. 11 Strain distribution of the targets

    表5 中統(tǒng)計(jì)了在靶體背部長/短軸方向上塑性應(yīng)變不超過0.2 的區(qū)域的最大坐標(biāo)值,定義長/短軸相對增量=(長/短軸最大坐標(biāo)值-長/短軸半軸長)÷(長/短軸半軸長),表征靶體在長/短軸方向塑性應(yīng)變區(qū)域的增大比例。表5 中數(shù)據(jù)表明,隨著長短軸長度比的增大,長軸相對增量逐漸減小,短軸相對增量逐漸增大,塑性應(yīng)變區(qū)域的不對稱性逐漸加劇。

    表5 靶體背部塑性應(yīng)變區(qū)域范圍對比Table 5 Comparison of plastic strain ranges on the back of targets

    3.3 極坐標(biāo)系下彈體橫截面長短軸長度比對靶體響應(yīng)特性的影響

    在彈體的侵徹理論模型中,靶體的徑向/法向應(yīng)力是學(xué)者們更關(guān)注的問題,其與靶體阻力函數(shù)密切相關(guān),直接影響理論模型的準(zhǔn)確性。因此,通過在數(shù)值模擬后處理軟件中建立極坐標(biāo)系的方法,分析不同橢圓橫截面長短軸長度比的彈體侵徹作用下靶體的響應(yīng)特性。忽略靶體自由面效應(yīng)及邊界效應(yīng)對計(jì)算結(jié)果的影響,僅提取彈體頭部貫入階段靶體沿周向各點(diǎn)的各項(xiàng)參數(shù),靶體測量點(diǎn)位置如圖12 所示,在400 m/s 的初始撞擊速度下,彈體頭部貫入階段持續(xù)約80 μs。

    圖12 數(shù)值模擬中靶體測點(diǎn)分布Fig. 12 Layout of measured points of targets in numerical simulation

    圖13~15 給出了不同長短軸長度比的橢圓橫截面彈體侵徹作用下,各測量點(diǎn)的徑向位移、徑向應(yīng)力及切向應(yīng)力的變化情況。從圖13~14 可以看出,圓形與橢圓橫截面彈體侵徹作用下靶體的響應(yīng)特性有顯著區(qū)別。從靶體的徑向位移和徑向應(yīng)力來看,前者在周向方向的變化趨勢基本一致,各測量點(diǎn)提取的數(shù)據(jù)曲線重合,表明靶體的響應(yīng)沿周向是均勻的;而后者在周向的各點(diǎn)表現(xiàn)出較大的差異,主要如下:(1)各測量點(diǎn)的徑向位移和徑向應(yīng)力峰值不同,峰值呈現(xiàn)出從短軸至長軸逐漸降低的趨勢;(2)各測量點(diǎn)的徑向位移和徑向應(yīng)力達(dá)到峰值的時(shí)間不同,長軸方向最先達(dá)到峰值,而在短軸方向最晚;(3)徑向位移和徑向應(yīng)力峰值大小和比值與彈體橫截面長短軸長度比有關(guān)。

    圖13 靶體沿周向方向的徑向位移Fig. 13 Radial displacement of targets along circumferential direction

    從圖15(a)可以發(fā)現(xiàn),靶體在各點(diǎn)的切向應(yīng)力基本為零,這是由于圓形橫截面彈體為旋轉(zhuǎn)對稱結(jié)構(gòu),在理想侵徹條件(彈體剛性、無傾角/攻角/偏航角)下,靶體在各方向的變形一致且僅沿著徑向方向,因此切向方向無變形,切向應(yīng)力為零;而橢圓橫截面彈體則是面對稱結(jié)構(gòu),導(dǎo)致靶體除了在彈體對稱面(長/短軸平面)上僅沿徑向變形,而在其余方向同時(shí)具有徑向變形和切向變形;因此切向應(yīng)力在0°和90°方向?yàn)榱?,?°與90°之間呈現(xiàn)先升高后降低的趨勢,如圖15(b)~(d)所示。

    圖14 靶體沿周向方向的徑向應(yīng)力Fig. 14 Radial stress of targets along circumferential direction

    圖15 靶體沿周向方向的切向應(yīng)力Fig. 15 Tangential stress of targets along circumferential direction

    為了直觀分析彈體橫截面長短軸長度比對靶體響應(yīng)特性的影響規(guī)律,將徑向位移和徑向應(yīng)力的峰值做無量綱化處理,同時(shí)提取切向應(yīng)力的峰值進(jìn)行研究,如圖16~18 所示。從圖中可以看到彈體橫截面長短軸長度比對靶體的徑向位移/應(yīng)力和切向應(yīng)力的影響顯著,長短軸長度比越大,徑向位移/應(yīng)力和切向應(yīng)力的不對稱性越明顯;徑向位移/應(yīng)力從長軸至短軸逐漸增加,同時(shí)變化梯度呈現(xiàn)先增加后減小的趨勢,類似于三角函數(shù)。而靶體的切向應(yīng)力(除圓形橫截面彈體外)則表現(xiàn)先增加后減小的變化規(guī)律,切向應(yīng)力關(guān)于45°方向基本呈現(xiàn)對稱分布;此外,隨著長短軸長度比的增大,切向應(yīng)力的峰值提升顯著。

    圖16 沿周向方向的無量綱徑向位移Fig. 16 Dimensionless radial displacement of targets along circumferential direction

    圖13~18 的數(shù)值模擬結(jié)果表明,在圓形橫截面彈體侵徹作用下,靶體沿周向方向各測量點(diǎn)的響應(yīng)曲線基本一致,靶體各點(diǎn)僅受徑向應(yīng)力作用而切向應(yīng)力為零,因此絕大部分圓形橫截面彈體的侵徹理論模型將彈靶作用過程簡化為一維是合理的;在橢圓橫截面彈體侵徹作用下,靶體沿周向方向各測量點(diǎn)的響應(yīng)曲線則存在明顯的非對稱特性,靶體各點(diǎn)徑向位移和徑向應(yīng)力峰值、峰值到達(dá)時(shí)間、峰值比值等與周向角密切相關(guān)。不考慮靶體反射形成的拉伸波的情況下,靶體也不再是處于簡單的壓縮狀態(tài),而是處于壓縮/剪切的復(fù)雜應(yīng)力狀態(tài),因此針對橢圓橫截面彈體的侵徹理論模型應(yīng)著重考慮靶體復(fù)雜應(yīng)力狀態(tài)對靶體阻力函數(shù)的影響。

    圖17 沿周向方向的無量綱徑向應(yīng)力Fig. 17 Dimensionless radial stress of targets along circumferential direction

    圖18 沿周向方向的最大切向應(yīng)力Fig. 18 The maximum tangential stress of targets along circumferential direction

    4 結(jié) 論

    開展了3 種不同橢圓橫截面彈體(彈體橫截面長短軸長度比β=1.00, 1.25, 1.61)在200~600 m/s 撞擊速度范圍內(nèi)正侵徹2A12 鋁板的實(shí)驗(yàn),獲得了2A12 鋁板的破壞特征及彈體的剩余速度。在此基礎(chǔ)上,建立了相應(yīng)的數(shù)值模擬模型,結(jié)合實(shí)驗(yàn)結(jié)果驗(yàn)證了模型的有效性,并系統(tǒng)分析了彈體橫截面長短軸長度比對靶體破壞情況及響應(yīng)特性的影響,主要結(jié)論如下。

    (1)彈體最大橫截面面積是影響彈體剩余速度的主要因素,彈體橫截面長短軸長度比對其剩余速度的影響較小。在工程應(yīng)用中,可直接使用圓形橫截面彈體正貫穿靶板的工程模型來計(jì)算具有相同最大橫截面面積的橢圓橫截面彈體的剩余速度。

    (2)彈體橫截面長短軸長度比對靶體的破壞程度有較大影響。當(dāng)β=1.00 時(shí),靶體背部形成的花瓣大小和形狀一致、空間分布均勻,而隨著長短軸長度比的增大,花瓣數(shù)量增加,尺寸變小,且短軸方向的花瓣數(shù)量、靶體表面隆起高度大于長軸方向的。

    (3)彈體橫截面長短軸長度比對靶體響應(yīng)特性有顯著影響。在圓形橫截面彈體侵徹作用下,靶體背部的塑性應(yīng)變區(qū)域?yàn)閳A形;而在橢圓橫截面彈體侵徹作用下,靶體背部的塑性應(yīng)變區(qū)域則近似于橢圓形,且沿短軸方向的塑性應(yīng)變區(qū)域明顯大于沿長軸方向的塑性應(yīng)變區(qū)域。

    (4)在圓形橫截面彈體侵徹作用下,靶體的徑向位移、徑向應(yīng)力和切向應(yīng)力與在橢圓橫截面彈體侵徹作用下靶體的徑向位移、徑向應(yīng)力和切向應(yīng)力顯著不同:前者靶體沿周向方向各點(diǎn)的應(yīng)力狀態(tài)變化規(guī)律基本一致,靶體處于簡單的壓縮狀態(tài),切向應(yīng)力為零;后者靶體各點(diǎn)的應(yīng)力狀態(tài)與彈體橫截面長短軸長度比和周向角密切相關(guān),靶體各點(diǎn)同時(shí)受到壓縮和剪切應(yīng)力的耦合作用。

    猜你喜歡
    靶體長軸彈體
    靶體結(jié)構(gòu)對前混合水射流噴丸強(qiáng)化應(yīng)力特性的影響
    尾錐角對彈體斜侵徹過程中姿態(tài)的影響研究
    超高韌性水泥基復(fù)合材料—纖維混凝土組合靶體抗兩次打擊試驗(yàn)研究*
    爆炸與沖擊(2022年3期)2022-04-11 03:14:32
    橢圓截面彈體斜侵徹金屬靶體彈道研究*
    爆炸與沖擊(2022年2期)2022-03-17 07:28:44
    單管立式長軸多級熔鹽泵的研發(fā)及應(yīng)用
    橢圓與兩焦點(diǎn)弦有關(guān)的幾個(gè)重要性質(zhì)及其推論
    2013年山東卷(理)壓軸題的推廣
    STOPAQ粘彈體技術(shù)在管道施工中的應(yīng)用
    上海煤氣(2018年6期)2018-03-07 01:03:22
    彈丸對預(yù)開孔混凝土靶體侵徹的實(shí)驗(yàn)研究
    LP(LT)型立式長軸排水泵的研制及應(yīng)用
    亚洲不卡免费看| 91成人精品电影| 国精品久久久久久国模美| 97超视频在线观看视频| 2022亚洲国产成人精品| 久久国产精品大桥未久av| 国产色爽女视频免费观看| 国产一级毛片在线| 精品国产国语对白av| 国国产精品蜜臀av免费| 桃花免费在线播放| 亚洲精品aⅴ在线观看| 人妻制服诱惑在线中文字幕| videos熟女内射| 天堂俺去俺来也www色官网| 人人妻人人添人人爽欧美一区卜| 亚洲综合色网址| 欧美精品一区二区免费开放| 大陆偷拍与自拍| 一区二区日韩欧美中文字幕 | 蜜桃在线观看..| 老女人水多毛片| 国产乱来视频区| 国产av国产精品国产| 大片电影免费在线观看免费| 如日韩欧美国产精品一区二区三区 | 久久久久精品性色| 丝袜美足系列| 国国产精品蜜臀av免费| 日韩制服骚丝袜av| 久久免费观看电影| 色网站视频免费| 亚洲精品美女久久av网站| 18+在线观看网站| videos熟女内射| 亚洲精品中文字幕在线视频| 99久久精品国产国产毛片| 人人澡人人妻人| 极品人妻少妇av视频| 夜夜看夜夜爽夜夜摸| 成年女人在线观看亚洲视频| 久久精品国产亚洲av天美| 国产亚洲一区二区精品| 少妇人妻精品综合一区二区| 插阴视频在线观看视频| 国产av国产精品国产| 91久久精品国产一区二区三区| 又大又黄又爽视频免费| 在线观看免费日韩欧美大片 | 国产av国产精品国产| 伦理电影免费视频| 日韩av免费高清视频| 成人国产av品久久久| 午夜老司机福利剧场| 免费大片18禁| 久久影院123| 亚洲国产欧美在线一区| 欧美bdsm另类| 最新的欧美精品一区二区| 交换朋友夫妻互换小说| 色网站视频免费| 免费看光身美女| 99国产精品免费福利视频| 国产午夜精品一二区理论片| 成人综合一区亚洲| 亚洲av在线观看美女高潮| 如何舔出高潮| 国产精品蜜桃在线观看| 九色亚洲精品在线播放| 国产精品嫩草影院av在线观看| 丝袜美足系列| 热re99久久国产66热| 色视频在线一区二区三区| av天堂久久9| 久久久精品94久久精品| 成人漫画全彩无遮挡| 国产亚洲午夜精品一区二区久久| xxx大片免费视频| 18禁裸乳无遮挡动漫免费视频| 桃花免费在线播放| 亚洲欧洲国产日韩| 国产精品秋霞免费鲁丝片| 制服诱惑二区| 国产精品 国内视频| 午夜福利网站1000一区二区三区| 久久久久久久久久人人人人人人| 国产白丝娇喘喷水9色精品| 欧美+日韩+精品| 国产精品国产av在线观看| 亚洲三级黄色毛片| 青春草视频在线免费观看| 97在线人人人人妻| 狂野欧美激情性bbbbbb| a级毛色黄片| 成年人午夜在线观看视频| 日韩欧美精品免费久久| 91久久精品国产一区二区成人| av天堂久久9| av天堂久久9| 久久久精品94久久精品| 国产精品欧美亚洲77777| 九色成人免费人妻av| 桃花免费在线播放| 七月丁香在线播放| 国产熟女午夜一区二区三区 | 91久久精品国产一区二区成人| 黄色怎么调成土黄色| 99久国产av精品国产电影| 综合色丁香网| 免费大片18禁| 最近中文字幕高清免费大全6| 观看美女的网站| 18禁在线播放成人免费| 免费观看无遮挡的男女| 性色avwww在线观看| 精品久久久精品久久久| 男人爽女人下面视频在线观看| 高清欧美精品videossex| 性色avwww在线观看| 午夜福利影视在线免费观看| 一区二区三区乱码不卡18| 欧美日韩国产mv在线观看视频| 少妇猛男粗大的猛烈进出视频| 黄色毛片三级朝国网站| 毛片一级片免费看久久久久| 毛片一级片免费看久久久久| 中文字幕精品免费在线观看视频 | 丰满少妇做爰视频| 免费观看性生交大片5| 久久97久久精品| 久久国产亚洲av麻豆专区| 亚洲成人手机| 国产片特级美女逼逼视频| 简卡轻食公司| 水蜜桃什么品种好| 国产精品久久久久久精品电影小说| 久久精品夜色国产| 最近的中文字幕免费完整| 免费黄色在线免费观看| 男人爽女人下面视频在线观看| 精品一区二区三区视频在线| 九色亚洲精品在线播放| 99热这里只有是精品在线观看| 欧美丝袜亚洲另类| 亚洲精品乱码久久久v下载方式| 精品亚洲成国产av| 丝袜喷水一区| 97超碰精品成人国产| 久久久久国产网址| 成人免费观看视频高清| 日韩熟女老妇一区二区性免费视频| 国产成人精品福利久久| 成人亚洲精品一区在线观看| 99精国产麻豆久久婷婷| 欧美最新免费一区二区三区| 国产av码专区亚洲av| 亚洲成人av在线免费| 国产有黄有色有爽视频| 成年美女黄网站色视频大全免费 | 久久韩国三级中文字幕| 一本久久精品| 人妻少妇偷人精品九色| 青春草国产在线视频| 蜜桃国产av成人99| 亚洲不卡免费看| 午夜福利视频在线观看免费| 秋霞伦理黄片| 不卡视频在线观看欧美| 亚洲经典国产精华液单| 国产免费一级a男人的天堂| 亚洲精品日韩在线中文字幕| 久久99精品国语久久久| 久久久久精品久久久久真实原创| 精品一区二区免费观看| 国产精品国产三级国产av玫瑰| 在线观看三级黄色| 国产视频首页在线观看| 91精品伊人久久大香线蕉| 精品人妻熟女av久视频| 少妇被粗大的猛进出69影院 | 久久这里有精品视频免费| 精品熟女少妇av免费看| 一边亲一边摸免费视频| 国产精品嫩草影院av在线观看| av黄色大香蕉| 久久人妻熟女aⅴ| 中国三级夫妇交换| 日韩免费高清中文字幕av| 亚洲丝袜综合中文字幕| 天天躁夜夜躁狠狠久久av| 多毛熟女@视频| 日韩一区二区三区影片| 午夜免费男女啪啪视频观看| 这个男人来自地球电影免费观看 | 久久精品国产自在天天线| 中文精品一卡2卡3卡4更新| 国产精品一区二区在线不卡| 亚洲欧美日韩卡通动漫| 亚洲欧美清纯卡通| 一级黄片播放器| 亚洲美女搞黄在线观看| 久久狼人影院| 免费看光身美女| 色网站视频免费| 三级国产精品欧美在线观看| 2022亚洲国产成人精品| 国产国语露脸激情在线看| xxxhd国产人妻xxx| 亚洲综合精品二区| 国产亚洲一区二区精品| 97在线人人人人妻| 日韩欧美一区视频在线观看| 国产成人精品福利久久| 高清黄色对白视频在线免费看| 亚洲精品视频女| 国产极品粉嫩免费观看在线 | 大片免费播放器 马上看| 曰老女人黄片| 欧美 日韩 精品 国产| 国产免费又黄又爽又色| 麻豆成人av视频| 亚洲国产精品一区二区三区在线| 亚洲激情五月婷婷啪啪| 免费播放大片免费观看视频在线观看| 97在线视频观看| videossex国产| 国产视频内射| 晚上一个人看的免费电影| 乱码一卡2卡4卡精品| 日韩一区二区视频免费看| 免费高清在线观看日韩| 国产精品国产三级专区第一集| 亚洲国产av影院在线观看| 日韩强制内射视频| 精品少妇内射三级| 久久精品久久久久久久性| 国产成人精品一,二区| kizo精华| 欧美精品一区二区大全| 欧美日韩视频精品一区| 亚洲综合色网址| 国产高清有码在线观看视频| 欧美 亚洲 国产 日韩一| 在线亚洲精品国产二区图片欧美 | 日日撸夜夜添| 欧美性感艳星| 色婷婷久久久亚洲欧美| 18禁动态无遮挡网站| 男人添女人高潮全过程视频| 成人免费观看视频高清| a级毛片在线看网站| 如日韩欧美国产精品一区二区三区 | 亚洲精品,欧美精品| 精品久久久久久电影网| 亚洲情色 制服丝袜| 我要看黄色一级片免费的| 少妇人妻精品综合一区二区| 国产欧美日韩一区二区三区在线 | 色吧在线观看| 免费大片18禁| 国产 精品1| 熟妇人妻不卡中文字幕| 青春草视频在线免费观看| 国产毛片在线视频| 边亲边吃奶的免费视频| 少妇人妻久久综合中文| 成人影院久久| 午夜老司机福利剧场| 51国产日韩欧美| 欧美 日韩 精品 国产| 青春草视频在线免费观看| 国产在线视频一区二区| 欧美精品一区二区大全| 国产成人免费观看mmmm| 日本爱情动作片www.在线观看| 国产老妇伦熟女老妇高清| 99久久人妻综合| 熟女av电影| 美女大奶头黄色视频| 十分钟在线观看高清视频www| 国产片内射在线| 夜夜骑夜夜射夜夜干| 亚洲人成网站在线播| 国产又色又爽无遮挡免| 午夜福利网站1000一区二区三区| 五月伊人婷婷丁香| 亚洲欧洲精品一区二区精品久久久 | 麻豆成人av视频| 亚洲一区二区三区欧美精品| 91国产中文字幕| 亚洲美女黄色视频免费看| 精品久久久久久电影网| 成人国语在线视频| 国产av国产精品国产| 久久人人爽人人片av| 成年女人在线观看亚洲视频| 亚洲成色77777| 亚洲av中文av极速乱| 各种免费的搞黄视频| 午夜免费鲁丝| 妹子高潮喷水视频| 国产淫语在线视频| 成人亚洲欧美一区二区av| 欧美人与性动交α欧美精品济南到 | 一区二区av电影网| 精品少妇内射三级| 狠狠精品人妻久久久久久综合| 国产精品不卡视频一区二区| 91久久精品电影网| 久久久久视频综合| 免费观看无遮挡的男女| 99热6这里只有精品| 国产亚洲精品第一综合不卡 | 亚洲欧美日韩另类电影网站| 久久久a久久爽久久v久久| 亚洲人成77777在线视频| 国产精品一国产av| 搡老乐熟女国产| 夜夜骑夜夜射夜夜干| 晚上一个人看的免费电影| 国产午夜精品一二区理论片| 国产成人精品无人区| 国产一区二区在线观看日韩| 国产亚洲最大av| 欧美xxxx性猛交bbbb| 午夜免费男女啪啪视频观看| 天天操日日干夜夜撸| 日韩制服骚丝袜av| 亚洲精品av麻豆狂野| 欧美亚洲日本最大视频资源| 亚洲无线观看免费| 亚洲精品视频女| 免费不卡的大黄色大毛片视频在线观看| 18禁在线播放成人免费| a 毛片基地| 成人综合一区亚洲| 免费观看在线日韩| 色婷婷av一区二区三区视频| 亚洲成人手机| 亚洲图色成人| 日韩,欧美,国产一区二区三区| 欧美日韩成人在线一区二区| 一级毛片电影观看| 简卡轻食公司| 亚洲三级黄色毛片| 熟女电影av网| 日韩大片免费观看网站| 丰满少妇做爰视频| 亚洲少妇的诱惑av| 亚洲av二区三区四区| 久久99精品国语久久久| 丰满乱子伦码专区| 少妇 在线观看| 日韩电影二区| 欧美最新免费一区二区三区| 日本wwww免费看| 国产成人免费无遮挡视频| 内地一区二区视频在线| 久久女婷五月综合色啪小说| 亚洲国产精品专区欧美| 99久久综合免费| 亚洲人成网站在线播| 亚洲久久久国产精品| av免费观看日本| 在线观看美女被高潮喷水网站| 亚洲欧美日韩另类电影网站| 免费看不卡的av| 欧美亚洲 丝袜 人妻 在线| 精品一区在线观看国产| 亚洲欧美中文字幕日韩二区| 日韩电影二区| 久久免费观看电影| 黑人欧美特级aaaaaa片| 欧美xxxx性猛交bbbb| 日韩欧美一区视频在线观看| 人成视频在线观看免费观看| 精品久久久久久久久av| a 毛片基地| 两个人免费观看高清视频| 亚洲av国产av综合av卡| 成人18禁高潮啪啪吃奶动态图 | 午夜视频国产福利| 2022亚洲国产成人精品| 国产精品一区二区在线观看99| 日日撸夜夜添| 亚洲av欧美aⅴ国产| 国产成人freesex在线| 在线亚洲精品国产二区图片欧美 | 中文字幕免费在线视频6| 亚洲欧美色中文字幕在线| 国产精品.久久久| 国产女主播在线喷水免费视频网站| 日本av免费视频播放| 日本猛色少妇xxxxx猛交久久| 国产视频首页在线观看| 美女福利国产在线| 18禁在线无遮挡免费观看视频| 久久精品久久久久久久性| h视频一区二区三区| 精品国产露脸久久av麻豆| 日韩av在线免费看完整版不卡| 欧美亚洲 丝袜 人妻 在线| 国产精品国产三级专区第一集| 26uuu在线亚洲综合色| 亚洲内射少妇av| 国产精品国产三级国产专区5o| 中国美白少妇内射xxxbb| 天天躁夜夜躁狠狠久久av| 人人妻人人澡人人看| av电影中文网址| 99热全是精品| 亚洲国产精品国产精品| av在线播放精品| 超色免费av| 大码成人一级视频| 亚洲国产成人一精品久久久| 国产又色又爽无遮挡免| 男女啪啪激烈高潮av片| 日韩制服骚丝袜av| 一区二区av电影网| 在线观看美女被高潮喷水网站| 中文字幕免费在线视频6| 欧美精品一区二区大全| 免费久久久久久久精品成人欧美视频 | 一区二区三区免费毛片| 午夜激情福利司机影院| 国产爽快片一区二区三区| av不卡在线播放| 亚洲av不卡在线观看| 国产精品久久久久久久久免| 91精品一卡2卡3卡4卡| 久久免费观看电影| av天堂久久9| 婷婷成人精品国产| 日日摸夜夜添夜夜添av毛片| 免费高清在线观看视频在线观看| av一本久久久久| 国产69精品久久久久777片| 欧美日本中文国产一区发布| 亚洲欧美成人综合另类久久久| 美女福利国产在线| 午夜影院在线不卡| 大陆偷拍与自拍| 少妇猛男粗大的猛烈进出视频| 国产成人一区二区在线| 观看美女的网站| 亚洲精品美女久久av网站| 啦啦啦中文免费视频观看日本| 国产白丝娇喘喷水9色精品| 王馨瑶露胸无遮挡在线观看| 欧美三级亚洲精品| 熟女人妻精品中文字幕| 国产精品偷伦视频观看了| 丰满饥渴人妻一区二区三| 哪个播放器可以免费观看大片| 亚洲精品国产av蜜桃| 精品一区二区免费观看| 欧美一级a爱片免费观看看| 色哟哟·www| 91精品一卡2卡3卡4卡| 女性被躁到高潮视频| 久久99一区二区三区| 国产精品久久久久久精品电影小说| 精品人妻偷拍中文字幕| 日本黄色片子视频| 日韩精品有码人妻一区| 日本91视频免费播放| 久久人人爽人人爽人人片va| 久久av网站| 久久97久久精品| 久久国产精品大桥未久av| 美女cb高潮喷水在线观看| 永久免费av网站大全| 久久精品夜色国产| av免费观看日本| 国产亚洲精品第一综合不卡 | 免费观看的影片在线观看| 永久免费av网站大全| 国产视频内射| 国产极品粉嫩免费观看在线 | 国产精品国产三级专区第一集| 国产精品久久久久久精品古装| 九九久久精品国产亚洲av麻豆| 又黄又爽又刺激的免费视频.| 有码 亚洲区| 日日摸夜夜添夜夜添av毛片| 午夜福利视频精品| 边亲边吃奶的免费视频| 亚洲国产精品成人久久小说| 91国产中文字幕| 精品亚洲乱码少妇综合久久| 日韩中文字幕视频在线看片| 国产一区有黄有色的免费视频| 午夜久久久在线观看| 久久亚洲国产成人精品v| 亚洲精品av麻豆狂野| 亚洲第一区二区三区不卡| 51国产日韩欧美| 男人操女人黄网站| 涩涩av久久男人的天堂| 国国产精品蜜臀av免费| 不卡视频在线观看欧美| 一本久久精品| 亚洲精品亚洲一区二区| 特大巨黑吊av在线直播| 免费不卡的大黄色大毛片视频在线观看| 日韩一区二区视频免费看| 9色porny在线观看| 国产精品一区二区在线观看99| 国产精品久久久久久精品电影小说| 亚洲色图 男人天堂 中文字幕 | 亚洲精品久久成人aⅴ小说 | 欧美精品高潮呻吟av久久| 久久久久久久久久久丰满| 亚洲精品色激情综合| 亚洲精华国产精华液的使用体验| 亚洲精品av麻豆狂野| 国产免费福利视频在线观看| 久久99热6这里只有精品| 一区二区av电影网| 欧美变态另类bdsm刘玥| 中文字幕av电影在线播放| 免费观看av网站的网址| 久久久久国产精品人妻一区二区| 精品国产乱码久久久久久小说| 又黄又爽又刺激的免费视频.| 亚洲av电影在线观看一区二区三区| 性高湖久久久久久久久免费观看| 国产精品嫩草影院av在线观看| av在线观看视频网站免费| 纯流量卡能插随身wifi吗| 国产一区亚洲一区在线观看| 熟女av电影| 人人澡人人妻人| 人人妻人人澡人人爽人人夜夜| av黄色大香蕉| 国产亚洲最大av| 成人免费观看视频高清| a级毛片黄视频| 亚洲色图 男人天堂 中文字幕 | 亚洲无线观看免费| 丰满少妇做爰视频| 欧美xxⅹ黑人| 亚洲av综合色区一区| 国产成人免费观看mmmm| 18+在线观看网站| 毛片一级片免费看久久久久| av网站免费在线观看视频| 纯流量卡能插随身wifi吗| 精品久久国产蜜桃| 国产免费现黄频在线看| 一级黄片播放器| 一区二区av电影网| 免费看不卡的av| 好男人视频免费观看在线| 久久免费观看电影| 久久精品国产鲁丝片午夜精品| 97在线视频观看| 日韩熟女老妇一区二区性免费视频| 熟女人妻精品中文字幕| 日韩制服骚丝袜av| 日韩成人伦理影院| av电影中文网址| 国产一区二区三区av在线| 久久精品国产鲁丝片午夜精品| 丰满乱子伦码专区| 免费av不卡在线播放| 免费黄色在线免费观看| 一级毛片 在线播放| a级片在线免费高清观看视频| av国产久精品久网站免费入址| 男女边摸边吃奶| 精品亚洲成a人片在线观看| 国产日韩一区二区三区精品不卡 | 国产黄片视频在线免费观看| 国产成人精品久久久久久| 婷婷色av中文字幕| 亚洲国产欧美在线一区| 一个人看视频在线观看www免费| 毛片一级片免费看久久久久| 2022亚洲国产成人精品| 国产一区有黄有色的免费视频| 最近中文字幕高清免费大全6| 男女免费视频国产| 亚洲国产欧美日韩在线播放| 色哟哟·www| 精品国产一区二区三区久久久樱花| 亚洲五月色婷婷综合| 亚洲国产精品专区欧美| 熟女av电影| 色婷婷久久久亚洲欧美| 欧美日韩精品成人综合77777| 80岁老熟妇乱子伦牲交| 精品少妇内射三级| 久久99精品国语久久久| av在线老鸭窝| 国产黄片视频在线免费观看| 久久久久久久久久久久大奶| 99热这里只有是精品在线观看| 欧美人与性动交α欧美精品济南到 | 天天躁夜夜躁狠狠久久av| 欧美成人午夜免费资源| 9色porny在线观看| 寂寞人妻少妇视频99o| 久久久久久久大尺度免费视频| 有码 亚洲区| 国产av码专区亚洲av| 婷婷色av中文字幕| 色哟哟·www| 日韩,欧美,国产一区二区三区| 丝袜脚勾引网站| 亚洲国产精品专区欧美| 人妻一区二区av| 久久久久网色| 久久久久久伊人网av|