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

    六棱柱形戰(zhàn)斗部預(yù)制破片驅(qū)動(dòng)的數(shù)值模擬與試驗(yàn)

    2016-07-15 01:59:41武敬博茍瑞君鄭俊杰王學(xué)雷趙燕兵
    火炸藥學(xué)報(bào) 2016年3期

    武敬博,茍瑞君,鄭俊杰,王學(xué)雷,趙燕兵

    (中北大學(xué)化工與環(huán)境學(xué)院, 山西 太原030051)

    ?

    六棱柱形戰(zhàn)斗部預(yù)制破片驅(qū)動(dòng)的數(shù)值模擬與試驗(yàn)

    武敬博,茍瑞君,鄭俊杰,王學(xué)雷,趙燕兵

    (中北大學(xué)化工與環(huán)境學(xué)院, 山西 太原030051)

    摘要:為研究六棱柱形戰(zhàn)斗部預(yù)制破片的殺傷效果,建立了正六棱柱戰(zhàn)斗部的三維模型,并以等高、等外徑的普通圓柱形戰(zhàn)斗部為對(duì)照組,利用LS-DYNA軟件分別模擬了二者在端面中心起爆與偏心兩線(xiàn)起爆時(shí)對(duì)破片的驅(qū)動(dòng)過(guò)程,分析了破片的速度與密度增益,并設(shè)計(jì)了實(shí)彈試驗(yàn)對(duì)模擬結(jié)果進(jìn)行了驗(yàn)證。模擬結(jié)果表明,偏心兩線(xiàn)起爆時(shí),六棱柱形戰(zhàn)斗部和圓柱形戰(zhàn)斗部破片速度的增益分別為9.2%和12.2%,與試驗(yàn)值的誤差均在10%以?xún)?nèi)。試驗(yàn)結(jié)果表明,與端面中心起爆的圓柱形戰(zhàn)斗部相比,六棱柱形戰(zhàn)斗部在端面中心起爆和偏心兩線(xiàn)起爆時(shí)分別可使破片密度提高53.6%和74.1%,且使破片在較遠(yuǎn)距離處仍有較優(yōu)的聚集效果。

    關(guān)鍵詞:爆炸力學(xué);六棱柱形戰(zhàn)斗部;密度增益;LS-DYNA軟件;預(yù)制破片;偏心起爆

    引 言

    通常評(píng)價(jià)戰(zhàn)斗部預(yù)制破片殺傷能力的指標(biāo)主要是破片速度與密度。不同結(jié)構(gòu)、裝藥比的戰(zhàn)斗部,速度的增益能力不同[1-2]。單點(diǎn)偏心起爆能顯著增加目標(biāo)方位的破片初速,但對(duì)目標(biāo)方向的密度提升作用較小。破片的飛散方向取決于戰(zhàn)斗部的形狀、起爆方式等,早期爆炸成形戰(zhàn)斗部通過(guò)削弱目標(biāo)部分的外殼,使破片更易向目標(biāo)方向飛散[3],但輔裝藥對(duì)主裝藥的影響卻難以控制;后期采取整體變形的思路,比如“膠囊”戰(zhàn)斗部[4]和D型戰(zhàn)斗部[5]等,這兩種戰(zhàn)斗部能顯著提升目標(biāo)方向的破片密度,若進(jìn)一步提高其爆炸控制系統(tǒng)的可靠性和準(zhǔn)確性,此類(lèi)戰(zhàn)斗部將在未來(lái)戰(zhàn)場(chǎng)上顯示出更大威力。聚焦戰(zhàn)斗部通過(guò)構(gòu)造內(nèi)腰鼓型的結(jié)構(gòu),縮小了破片的飛散角,形成密集的殺傷帶[6],以增強(qiáng)毀傷性能。為滿(mǎn)足戰(zhàn)斗部互換性與專(zhuān)業(yè)性的要求,需要一種結(jié)構(gòu)簡(jiǎn)單、易于實(shí)現(xiàn)的戰(zhàn)斗部。LI Yuan等[7]在研究戰(zhàn)斗部結(jié)構(gòu)與起爆點(diǎn)對(duì)破片殺傷效應(yīng)的影響時(shí),指出了六棱柱結(jié)構(gòu)的戰(zhàn)斗部破片在速度、動(dòng)能以及密度方面的優(yōu)勢(shì)。

    本研究對(duì)六棱柱形戰(zhàn)斗部和傳統(tǒng)圓柱形戰(zhàn)斗部的破片驅(qū)動(dòng)過(guò)程進(jìn)行了數(shù)值模擬,對(duì)比分析了二者破片速度和飛散情況;并通過(guò)實(shí)彈試驗(yàn),采用測(cè)速和統(tǒng)計(jì)靶板彈孔分布的方法驗(yàn)證了模擬結(jié)果,以期為戰(zhàn)斗部殺傷元增益技術(shù)的發(fā)展提供參考。

    1數(shù)值模擬

    1.1戰(zhàn)斗部結(jié)構(gòu)

    本研究旨在考察棱柱形戰(zhàn)斗部定向方向的破片增益效果,為提高效率,仿真建模和試驗(yàn)時(shí)破片只附著1/2戰(zhàn)斗部。設(shè)置圓柱組與六棱柱組以作對(duì)比,戰(zhàn)斗部的具體尺寸參數(shù)如圖1(a)和(b)所示,除角度外圖中單位均為mm。戰(zhàn)斗部主要由預(yù)制破片、襯筒、主裝藥3部分構(gòu)成,破片為單枚質(zhì)量2g、直徑6mm的鎢制小球,襯筒為2mm厚的2A12鋁制殼體,炸藥為B炸藥,裝藥密度1.69g/cm3。根據(jù)上述參數(shù)計(jì)算得出,棱柱形戰(zhàn)斗部和圓柱形戰(zhàn)斗部的主裝藥質(zhì)量分別約為3062g和3520g,前者約輕16%。

    為比較不同工況下破片的飛散情況,分別設(shè)置上端面中心起爆與偏心線(xiàn)起爆兩種方式。偏心線(xiàn)起爆為沿戰(zhàn)斗部豎直方向的兩列同步起爆,兩起爆線(xiàn)的夾角為60°,每列設(shè)4個(gè)等距的起爆點(diǎn),共8個(gè)起爆點(diǎn)(如圖1(c)所示),定義此時(shí)兩列起爆點(diǎn)所夾面的正對(duì)方向?yàn)槎ㄏ蚍较?,設(shè)定向區(qū)的中線(xiàn)為0°位,則定向區(qū)為-30°~30°的范圍。

    圖1 戰(zhàn)斗部尺寸參數(shù)及偏心起爆示意圖Fig.1 Warhead size parameters and schematic diagramof eccentric initiation

    1.2算法分析及有限元模型的建立

    采用流固耦合算法模擬炸藥與襯筒、破片的作用,炸藥周?chē)畛淇諝饨橘|(zhì),二者劃分為Euler網(wǎng)格,破片、襯筒劃分為L(zhǎng)agrange網(wǎng)格,放置在Euler網(wǎng)格中。襯筒與破片之間設(shè)置自動(dòng)面面接觸,對(duì)空氣域的外邊界施加壓力輸出條件以消除邊界效應(yīng);鎢球之間定義為自動(dòng)單面接觸,單位采用cm-g-μs,建立全尺寸模型進(jìn)行計(jì)算,如圖2所示。

    1.3材料參數(shù)

    炸藥采用高能炸藥計(jì)算模型(*MAT_HIGH_EXPLOSIVE_BURN),狀態(tài)方程使用*EOS_JWL狀態(tài)方程,具體參數(shù)如表1所示,其中ρ、D和pCJ分別為炸藥的密度、爆速和C-J壓力,A、B、R1、R2、ω為與炸藥性質(zhì)有關(guān)的常數(shù)。襯筒、破片采用*MAT_JOHNSON_COOK材料模型和*EOS_GRUNEISEN狀態(tài)方程。通過(guò)材料模型中自帶的失效參數(shù)的控制,在單元畸變過(guò)大時(shí)刪除,參數(shù)見(jiàn)表2,其中ρ、E、μ、G分別為材料的密度、楊氏模量、泊松比和剪切模量,A為歸一化內(nèi)聚強(qiáng)度,D1為失效參數(shù)。空氣采用*MAT_NULL模型和線(xiàn)性多項(xiàng)式狀態(tài)方程,參數(shù)見(jiàn)表3,其中ρ為空氣密度,C0~C5為空氣相關(guān)的多項(xiàng)式方程系數(shù)。部分參數(shù)取自文獻(xiàn)[8-9]。

    表2 破片與鋁殼的模型參數(shù)

    表3 空氣模型與狀態(tài)方程參數(shù)

    1.4模擬結(jié)果及分析

    炸藥起爆后,爆轟波以球面波形式傳播,襯筒殼體開(kāi)始膨脹,在極短的時(shí)間內(nèi)殼體破裂后,破片向外飛散,傳播過(guò)程如圖3所示。由圖3可知,150μs時(shí)鎢球的速度已趨于穩(wěn)定,按一定角度飛散。由于受到端面稀疏波的影響,破片速度呈現(xiàn)兩端小、中間大的變化梯度,反映到圖3中,不同速度的破片形成了向外的類(lèi)似弧形面的結(jié)構(gòu)。

    圖3 各起爆方式下兩種戰(zhàn)斗部在不同時(shí)刻的破片飛散情況Fig.3 Fragment dispersions of the two kinds of warheadsat different moments under various initiation modes

    1.4.1速度對(duì)比分析

    通過(guò)端面中心起爆和偏心兩線(xiàn)8點(diǎn)起爆的方式,考察六棱柱形戰(zhàn)斗部的破片速度增益。在LSPP后處理中,提取兩種起爆方式下六棱柱與圓柱形戰(zhàn)斗部的破片平均速度,結(jié)果如圖4所示。從圖4可以看出,在端面中心起爆時(shí),六棱柱形戰(zhàn)斗部破片平均速度約為970m/s,而圓柱形戰(zhàn)斗部破片速度約為998m/s;偏心兩線(xiàn)起爆時(shí),六棱柱形戰(zhàn)斗部破片平均速度約為1089m/s,較其端面中心起爆約提高9.2%,圓柱形戰(zhàn)斗部破片速度約為1120m/s,較其端面中心起爆提高約12.2%。在兩種起爆方式下,與圓柱形戰(zhàn)斗部相比,六棱柱形戰(zhàn)斗部的破片速度都僅小2.8%左右,考慮到其主裝藥量少約16%,可認(rèn)為六棱柱形戰(zhàn)斗部對(duì)破片速度有較優(yōu)的增益能力。

    圖4 各起爆方式下兩種戰(zhàn)斗部破片的平均速度Fig.4 Average velocity of fragments of the two kinds ofwarheads under various initiation modes

    (1)

    端面起爆時(shí),距起爆端不同距離處的破片初速可由式(2)[10]求得

    vx=[1-0.3615exp(-1.111x/dx)]×

    {1-0.1925exp[-3.03(l-x)/dx]}×vM

    (2)

    式中:dx為距離起爆端x距離處的裝藥直徑。

    將戰(zhàn)斗部相關(guān)參數(shù)帶入式(1)、式(2),求得六棱柱形戰(zhàn)斗部端面中心起爆時(shí)的破片平均速度約為1185m/s,比模擬結(jié)果大20%。由于格尼公式是采用如產(chǎn)物均勻膨脹、能量全部作用于破片飛散等假設(shè)推導(dǎo)出來(lái)的,并沒(méi)有考慮厚度不同的殼體膨脹速度的差距,且爆炸能量并非全部作用于破片飛散,因此求得的破片速度偏大,同時(shí)本研究所述的六棱柱形戰(zhàn)斗部的爆轟產(chǎn)物泄出十分明顯,這也限制了格尼公式的使用。因此仍需更多的實(shí)驗(yàn)數(shù)據(jù)來(lái)進(jìn)一步修正可適用于六棱柱形戰(zhàn)斗部的工程計(jì)算公式。

    1.4.2破片飛散對(duì)比分析

    從圖3可以看出,圓柱形戰(zhàn)斗部的破片沿圓周均勻分散,而六棱柱形戰(zhàn)斗部的破片在速度穩(wěn)定之后,呈現(xiàn)出明顯的緊密聚集的現(xiàn)象,類(lèi)似于“殺傷帶”的破片群在命中目標(biāo)后,將造成更有效的殺傷作用。

    未爆時(shí),定向區(qū)內(nèi)的破片排布在與戰(zhàn)斗部豎直鉛垂面(圖2中YOZ平面)夾角為-30°~30°的范圍內(nèi),共計(jì)40行,每行9個(gè)鎢珠。取戰(zhàn)斗部赤道面的破片進(jìn)行分析,依次編號(hào)為1~9。提取150μs時(shí)1~9號(hào)破片的徑向速度分量與合速度,計(jì)算破片飛散方向偏離YOZ平面正向的夾角θ,如圖5所示。

    從圖5可以看出,六棱柱形戰(zhàn)斗部各個(gè)破片偏離的程度明顯較小,偏離角的極差約為30°,且相鄰破片之間的差異亦較小,約為3°,即定向方向上的破片集中在30°張角的范圍內(nèi),且飛散至更遠(yuǎn)處時(shí)

    仍能保持較好的聚集效果。

    圖5 各起爆方式下兩種戰(zhàn)斗部破片的偏離角度Fig.5 Fragment drift-off angles of the two kinds ofwarheads under various initiation modes

    為更直觀地說(shuō)明,將150μs時(shí)定向區(qū)全部破片的飛散方向偏離YOZ面的角度分布情況進(jìn)行統(tǒng)計(jì),如圖6所示,橫坐標(biāo)為破片合速度的方向偏離豎直面的夾角θ,縱坐標(biāo)為該對(duì)應(yīng)角度下的破片數(shù)量。

    圖6 各起爆方式下不同徑向偏離角對(duì)應(yīng)的破片數(shù)量統(tǒng)計(jì)

    從圖6中可以看出,與圓柱形戰(zhàn)斗部破片的松散分布相比,六棱柱形戰(zhàn)斗部幾乎所有破片偏離豎直面的角度均在約-15°~15°的區(qū)間內(nèi),破片密集程度很高。同時(shí)可發(fā)現(xiàn),偏心兩線(xiàn)起爆時(shí),小角度區(qū)間內(nèi)的破片數(shù)量增多,破片的密集度提高。偏心兩線(xiàn)起爆的情況下,爆轟波發(fā)生碰撞,波陣面的方向改變,傳播至破片時(shí),各個(gè)破片所受驅(qū)動(dòng)力方向(即切線(xiàn)的法向方向)的張角變小,使得破片能以較密集的趨勢(shì)向外飛散。

    2試驗(yàn)驗(yàn)證

    2.1實(shí)彈參數(shù)

    為驗(yàn)證模擬結(jié)果,針對(duì)前述兩種形狀的戰(zhàn)斗部分別進(jìn)行端面中心和偏心兩線(xiàn)兩種方式下的靜爆試驗(yàn),實(shí)彈參數(shù)如表4所示。樣彈的規(guī)格、各部分的材料與前文一致,僅環(huán)繞1/2破片。在頂面作孔,引爆孔中的起爆藥柱實(shí)現(xiàn)端面中心起爆;側(cè)表面上做8個(gè)等徑的圓形槽(同圖1(c)所示),放置同規(guī)格的起爆藥柱,利用一入八出的爆炸邏輯網(wǎng)絡(luò)實(shí)現(xiàn)同步偏心兩線(xiàn)起爆。

    表4 實(shí)彈參數(shù)

    2.2現(xiàn)場(chǎng)布置

    根據(jù)扇形靶試驗(yàn)方法,在戰(zhàn)斗部的定向方向排布靶板,現(xiàn)場(chǎng)布置實(shí)況及示意圖如圖7所示。試驗(yàn)時(shí),樣彈立式放置,起爆端朝上,彈體幾何中心與靶板中心高度平齊;距彈體10m處設(shè)置靶板,靶板材質(zhì)為Q235鋼,規(guī)格為2m×1m×8mm(高×寬×厚),編號(hào)為1~7;考慮到戰(zhàn)斗部的對(duì)稱(chēng)性,為減少誤差,重點(diǎn)統(tǒng)計(jì)爆心發(fā)散30°范圍內(nèi)的破片分布,即只統(tǒng)計(jì)1~5號(hào)靶板的彈孔數(shù)目,而將6、7號(hào)作為參考靶;在13.5m處設(shè)置一塊獨(dú)立靶,以考察更遠(yuǎn)處六棱柱形戰(zhàn)斗部所驅(qū)動(dòng)破片的數(shù)量。分別在距爆心11.0m和7.5m處設(shè)置兩個(gè)激光測(cè)速儀,通過(guò)激光發(fā)射系統(tǒng)、反射膜、激光接收系統(tǒng)等形成矩形激光光幕,通過(guò)檢測(cè)破片穿越光幕時(shí)引發(fā)的光通量變化時(shí)間,可計(jì)算得出該破片的即時(shí)速度。

    圖7 試驗(yàn)現(xiàn)場(chǎng)布置情況及示意圖Fig.7 Setup and schematic diagram of the test site

    2.3試驗(yàn)結(jié)果與分析

    2.3.1 速度增益分析

    依照球形破片速度的衰減公式[3],由測(cè)速點(diǎn)的

    速度反推計(jì)算得出破片初速,去除明顯偏小值,獲得平均速度,結(jié)果見(jiàn)表5。從表5可知,模擬值均較實(shí)測(cè)值略小,其原因在于采用ALE算法模擬爆轟產(chǎn)物膨脹擴(kuò)大時(shí),產(chǎn)物的前沿面較為平滑,破片的速度曲線(xiàn)沒(méi)有明顯的突躍變化;模擬時(shí),鋁殼襯筒在變形過(guò)大的情況下失效刪除,部分能量提前逸散,但實(shí)際中鋁殼經(jīng)歷了膨脹炸裂的過(guò)程,同時(shí)破碎的殼體可能穿過(guò)測(cè)速光幕,這也會(huì)影響速度平均值的精確度。模擬值與實(shí)測(cè)值的誤差在10%以?xún)?nèi),可認(rèn)為模擬值是可信的。

    表5 破片速度實(shí)測(cè)值與模擬值對(duì)比

    2.3.2 密度增益分析

    破片穿靶效果如圖8所示。從圖8可以看出,大多數(shù)破片均能穿透靶板且形成圓形通孔,少數(shù)破片未穿透而嵌在靶板上,這些均為有效數(shù)據(jù);此外還有極少數(shù)的微小孔和大型不規(guī)則孔,微小孔可能是爆炸沖擊波吹起的碎石子或其他雜物擊中靶板留下的痕跡,而大型不規(guī)則孔是由撕裂的金屬雜物造成的,這些均會(huì)對(duì)統(tǒng)計(jì)結(jié)果造成影響,應(yīng)予以剔除。統(tǒng)計(jì)各靶板的彈孔數(shù)目列于表6。

    圖8 破片著靶效果實(shí)物圖Fig.8 Picture of target plate penetrated by fragments

    戰(zhàn)斗部形狀起爆方式彈孔數(shù)目/個(gè)1號(hào)靶2號(hào)靶3號(hào)靶4號(hào)靶5號(hào)靶總計(jì)破片密度/(個(gè)·m-2)密度增益/%六棱柱端面中心334246292217217.253.6六棱柱偏心兩線(xiàn)124071472519519.574.1圓柱端面中心262913261811211.20圓柱偏心兩線(xiàn)373825343416816.850.0

    由表6可知,不同起爆方式下,六棱柱形戰(zhàn)斗部的破片密度均高于圓柱形戰(zhàn)斗部。以圓柱端面中心起爆為基準(zhǔn),單點(diǎn)起爆和兩線(xiàn)起爆分別可使破片密度提高53.6%和74.1%。此外試驗(yàn)結(jié)果也證明了偏心兩線(xiàn)起爆能增大目標(biāo)方向破片密度的推論。

    分析3號(hào)靶板可發(fā)現(xiàn),在13.5m處,六棱柱形戰(zhàn)斗部形成的破片數(shù)目要遠(yuǎn)遠(yuǎn)多于圓柱形戰(zhàn)斗部,該現(xiàn)象證明了六棱柱形戰(zhàn)斗部具有較好的破片聚集能力,在較遠(yuǎn)距離處仍能保證充足的殺傷元數(shù)量。

    六棱柱形戰(zhàn)斗部的殼體存在轉(zhuǎn)折角,該處有應(yīng)力集中的現(xiàn)象,爆轟產(chǎn)物傳播至此處時(shí),會(huì)造成轉(zhuǎn)折角處優(yōu)先產(chǎn)生裂痕,而此時(shí)產(chǎn)物對(duì)破片的驅(qū)動(dòng)加速過(guò)程尚未結(jié)束。爆轟波靠近裂痕處的部分向外溢出,而主體仍繼續(xù)推動(dòng)破片運(yùn)動(dòng),破片受到的合力近似地垂直于棱柱側(cè)平面,從而使其徑向分散程度變小,破片以更密集的姿態(tài)飛向目標(biāo)。

    與破片軸向飛散戰(zhàn)斗部、聚焦戰(zhàn)斗部等情況類(lèi)似,六棱柱形戰(zhàn)斗部某個(gè)方向的破片聚集,必然會(huì)造成兩個(gè)破片密集區(qū)域之間存在一定的殺傷元稀疏區(qū),六棱柱棱角處的破片很少,在沒(méi)有外殼的模擬條件下,相鄰破片密集區(qū)之間有近30°的死角,這一問(wèn)題的解決有賴(lài)于精確制導(dǎo)和引戰(zhàn)配合技術(shù)的發(fā)展。六棱柱形戰(zhàn)斗部能夠顯著提高目標(biāo)方向的破片密度,同時(shí)結(jié)構(gòu)簡(jiǎn)單,容易實(shí)現(xiàn),便于和偏心起爆技術(shù)相結(jié)合,同時(shí)在保證殺傷能力的基礎(chǔ)上減輕了質(zhì)量,對(duì)戰(zhàn)斗部的發(fā)展有積極作用。

    3結(jié)論

    (1) 采用數(shù)值模擬得到六棱柱形戰(zhàn)斗部對(duì)預(yù)制破片的驅(qū)動(dòng)過(guò)程,所得的破片初速、觀測(cè)區(qū)的破片分布情況與試驗(yàn)結(jié)果基本吻合,模擬結(jié)果可信。

    (2) 相同裝藥高度和外徑時(shí),與圓柱形戰(zhàn)斗部相比,六棱柱形戰(zhàn)斗部裝藥質(zhì)量少16%,但在同一起爆方式下,二者的破片速度僅相差2.8%,而偏心兩線(xiàn)起爆可使二者破片的速度分別提升9.2%和12.2%,差距較小。綜合考慮,可認(rèn)為六棱柱形戰(zhàn)斗部對(duì)破片速度有較優(yōu)的增益能力。

    (3) 以圓柱形戰(zhàn)斗部端面中心起爆為基準(zhǔn),六棱柱形戰(zhàn)斗部的端面中心和偏心兩線(xiàn)起爆、圓柱形戰(zhàn)斗部?jī)删€(xiàn)起爆,分別可使破片密度提高53.6%、74.1%和50.0%。偏心兩線(xiàn)起爆方式有助于提高破片密度。六棱柱形戰(zhàn)斗部生成的破片群有良好的聚集性,能有效增加破片著靶的數(shù)量,且在較遠(yuǎn)距離處仍能保證良好的殺傷密度。

    參考文獻(xiàn):

    [1]葉小軍, 韓玉, 陳慶寶. 偏心起爆戰(zhàn)斗部速度增益的數(shù)值模擬及實(shí)驗(yàn)[J].火炸藥學(xué)報(bào),2009,32(3):29-34.YE Xiao-jun, HAN Yu, CHEN Qing-bao. Numerical simulation and experiment of velocity gains on the non-central detonation warhead [J]. Chinese Journal of Explosives & Propellants(Huozhayao Xuebao),2009, 32(3):29-34.

    [2]Wang Ma-fa, Lu Fang-yun, Li Xiang-yu.et al. A formula for calculating the velocities of fragments from velocity enhanced warhead [J] .Propellants, Explosives, Pyrotechnics, 2013(38):232-237.

    [3]盧芳云,李翔宇,林玉亮.戰(zhàn)斗部結(jié)構(gòu)與原理[M].北京:科學(xué)出版社,2009.

    [4]滕璽,米雙山.定向戰(zhàn)斗部的現(xiàn)狀分析與發(fā)展[J]. 飛航導(dǎo)彈,2014(4):89-94.

    TENG Xi, MI Shuang-shan. Study on status and development of directional warhead [J]. Aerodynamic Missile Journal, 2014(4):89-94.

    [5]龔柏林,盧芳云,李翔宇.D型預(yù)制破片戰(zhàn)斗部破片飛散過(guò)程的數(shù)值模擬[J]. 彈箭與制導(dǎo)學(xué)報(bào), 2010, 30(1):88-90.

    GONG Bai-lin, LU Fang-yun, LI Xiang-yu. Simulation and study on the fragment ejection process of premade D-shape warhead [J].Journal of Projectiles, Rockets, Missiles and Guidance,2010, 30(1):88-90.

    [6]王娟娟,郭雙鋒,袁寶慧,等. 聚焦戰(zhàn)斗部破片聚焦性能的數(shù)值計(jì)算和試驗(yàn)研究[J]. 火炸藥學(xué)報(bào),2013,36(2):87-90.WANG Juan-juan, GUO Shuang-feng, YUAN Bao-hui, et al. Numerical simulation and experiment study on directional-focused fragmentation warhead [J]. Chinese Journal of Explosive & Propellants(Huozhayao Xuebao), 2013, 36(2):87-90.

    [7]LI Yuan, WEN Yu-quan. Simulation of damage effectiveness of hexagonal prism aimable warhead with multi-point initiations [J].Journal of Beijing Institute of Technology, 2014(23):1-7.

    [8]趙海鷗. LS-DYNA動(dòng)力分析指南[M].北京:兵器工業(yè)出版社,2003.

    [9]楊亞?wèn)|,李向東,葉小軍,等. 可變形定向戰(zhàn)斗部變形控制參數(shù)的匹配[J].火炸藥學(xué)報(bào),2012,35(3):61-65,79.YANG Ya-dong, LI Xiang-dong YE Xiao-jun, et al. Matching of deformation control parameters of deformable warheads[J]. Chinese Journal of Explosive & Propellants(Huozhayao Xuebao), 2012, 35(3): 61-65, 79.

    [10]張寶平, 張慶明, 黃風(fēng)雷.爆轟物理學(xué)[M].北京:兵器工業(yè)出版社,2009: 317-319.

    Numerical Simulation and Experiment of Premade Fragments Droved by Hexagonal Prism Shaped Warhead

    WU Jing-bo, GOU Rui-jun, ZHENG Jun-jie, WANG Xue-lei, ZHAO Yan-bing

    (School of Chemical Engineering and Environment, North University of China, Taiyuan 030051, China)

    Abstract:To study the damage effects of premade fragment formed by hexagonal prism shaped warhead, a 3D model of hexagonal prism shaped warhead was established and an ordinary cylinder warhead with same height and same outer diameter was set as control group. The driven processes of fragments under end-surface initiation and eccentric two lines initiation for both were simulated by LS-DYNA software, respectively, and the simulated gains of fragment velocity and density were analyzed. An actual test was designed to verify the simulated results. The simulated results show that the gain of fragments velocity of hexagonal prism shaped warhead and cylindrical warhead is 9.2% and 12.2% respectively under eccentric two lines initiation, the errors with experiment values are all less than 10%. The experimental results show that compared with cylinder warhead under end-surface initiation, the hexagonal prism shaped warhead makes the fragment density improve by 53.6% and 74.1% under end-surface initiation and eccentric two lines initiation respectively and the fragment still have a better focus effect in a long distance.

    Keywords:explosion mechanics; hexagonal prism shaped warhead; density gain; LS-DYNA software; premade fragment;eccentric initiation

    DOI:10.14077/j.issn.1007-7812.2016.03.018

    收稿日期:2016-01-16; 2016-01-31

    基金項(xiàng)目:武敬博(1992-),男, 碩士研究生,從事戰(zhàn)斗部增益設(shè)計(jì)研究。E-mail: zbwujingbo@163.com

    作者簡(jiǎn)介:茍瑞君(1968-),女,教授,從事武器系統(tǒng)對(duì)抗技術(shù)和現(xiàn)代爆炸技術(shù)研究。E-mail: grjzsh@163.com

    中圖分類(lèi)號(hào):TJ414

    文獻(xiàn)標(biāo)志碼:A

    文章編號(hào):1007-7812(2016)03-0089-06

    国产一级毛片七仙女欲春2 | 美国免费a级毛片| 国产高清videossex| 亚洲中文日韩欧美视频| 欧美+亚洲+日韩+国产| 美女免费视频网站| 琪琪午夜伦伦电影理论片6080| 欧美日韩乱码在线| 亚洲一码二码三码区别大吗| 嫩草影院精品99| 精品国产美女av久久久久小说| 亚洲国产欧美日韩在线播放| 狂野欧美激情性xxxx| 国产成人影院久久av| 女生性感内裤真人,穿戴方法视频| 日韩欧美一区二区三区在线观看| 欧美日韩精品网址| 欧美不卡视频在线免费观看 | 久久久国产成人精品二区| 这个男人来自地球电影免费观看| 香蕉国产在线看| 亚洲一码二码三码区别大吗| 香蕉国产在线看| 老司机午夜十八禁免费视频| 日本在线视频免费播放| 午夜久久久在线观看| 午夜日韩欧美国产| 男男h啪啪无遮挡| 日韩一卡2卡3卡4卡2021年| 一边摸一边抽搐一进一小说| 欧美黄色淫秽网站| 亚洲第一欧美日韩一区二区三区| 桃色一区二区三区在线观看| 国产高清videossex| 天堂影院成人在线观看| 日日摸夜夜添夜夜添小说| 成人欧美大片| 中文资源天堂在线| 手机成人av网站| 身体一侧抽搐| 成人av一区二区三区在线看| 亚洲精品在线美女| 一级a爱视频在线免费观看| 又紧又爽又黄一区二区| 999久久久精品免费观看国产| 天天躁夜夜躁狠狠躁躁| 此物有八面人人有两片| 色在线成人网| 久久天堂一区二区三区四区| 国产黄a三级三级三级人| 一个人免费在线观看的高清视频| 女生性感内裤真人,穿戴方法视频| 国产一区二区三区在线臀色熟女| 性色av乱码一区二区三区2| 母亲3免费完整高清在线观看| 久久精品成人免费网站| 欧美日韩一级在线毛片| 午夜免费激情av| 国内精品久久久久久久电影| 级片在线观看| 在线观看舔阴道视频| 在线观看舔阴道视频| 亚洲精品一卡2卡三卡4卡5卡| 999精品在线视频| 国产精品影院久久| 精品电影一区二区在线| 最好的美女福利视频网| 久久国产精品男人的天堂亚洲| 在线观看免费视频日本深夜| 侵犯人妻中文字幕一二三四区| 欧美乱色亚洲激情| 老熟妇乱子伦视频在线观看| 欧美成人午夜精品| 午夜福利欧美成人| 一进一出好大好爽视频| 亚洲午夜精品一区,二区,三区| 两性午夜刺激爽爽歪歪视频在线观看 | 成年人黄色毛片网站| 国产成人精品久久二区二区91| 亚洲一区二区三区不卡视频| 欧美久久黑人一区二区| 露出奶头的视频| 欧美日韩中文字幕国产精品一区二区三区| 丝袜美腿诱惑在线| 女警被强在线播放| 精品国产乱子伦一区二区三区| 超碰成人久久| 亚洲国产欧美日韩在线播放| 黄色女人牲交| 狠狠狠狠99中文字幕| 日本一区二区免费在线视频| 日韩欧美国产在线观看| 亚洲专区国产一区二区| 国产真人三级小视频在线观看| 精品国产乱码久久久久久男人| 老司机午夜十八禁免费视频| 男男h啪啪无遮挡| 999精品在线视频| 搡老岳熟女国产| 啦啦啦 在线观看视频| 亚洲国产欧美一区二区综合| 校园春色视频在线观看| 黑人欧美特级aaaaaa片| 欧美成人免费av一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 日日爽夜夜爽网站| 精品电影一区二区在线| 亚洲aⅴ乱码一区二区在线播放 | 一级毛片高清免费大全| www.自偷自拍.com| 国产一区二区激情短视频| 午夜老司机福利片| 在线十欧美十亚洲十日本专区| 午夜a级毛片| 欧美成人一区二区免费高清观看 | 这个男人来自地球电影免费观看| 色综合婷婷激情| 成人18禁在线播放| 国产亚洲av嫩草精品影院| 久久人妻福利社区极品人妻图片| 神马国产精品三级电影在线观看 | 一级毛片女人18水好多| 日本在线视频免费播放| 97超级碰碰碰精品色视频在线观看| 久久中文字幕一级| 男人舔女人下体高潮全视频| 国产午夜福利久久久久久| 亚洲最大成人中文| 老司机福利观看| 欧美性猛交╳xxx乱大交人| 村上凉子中文字幕在线| 久久久久久久久免费视频了| 精品欧美一区二区三区在线| 成熟少妇高潮喷水视频| 亚洲精品av麻豆狂野| 黄网站色视频无遮挡免费观看| 在线观看一区二区三区| 亚洲熟妇熟女久久| 精品久久久久久久毛片微露脸| 母亲3免费完整高清在线观看| 热99re8久久精品国产| 免费在线观看亚洲国产| 国产一区二区三区在线臀色熟女| 此物有八面人人有两片| 岛国在线观看网站| 日韩大尺度精品在线看网址| 成年女人毛片免费观看观看9| 一级毛片女人18水好多| 一区二区三区激情视频| 中文字幕人成人乱码亚洲影| 99国产综合亚洲精品| 99国产综合亚洲精品| 欧美日韩乱码在线| 日本精品一区二区三区蜜桃| 波多野结衣av一区二区av| 久久国产精品影院| 熟女少妇亚洲综合色aaa.| 精品久久久久久成人av| 久久久久久久午夜电影| 日韩精品青青久久久久久| 成人三级做爰电影| 亚洲,欧美精品.| 欧美性猛交黑人性爽| 国产高清有码在线观看视频 | 婷婷精品国产亚洲av在线| 999久久久精品免费观看国产| 精品一区二区三区视频在线观看免费| 麻豆成人av在线观看| 久久国产精品人妻蜜桃| 精品无人区乱码1区二区| 亚洲国产欧洲综合997久久, | 国产麻豆成人av免费视频| 午夜老司机福利片| 91麻豆av在线| 日本在线视频免费播放| 免费在线观看黄色视频的| 美女 人体艺术 gogo| 欧美色视频一区免费| 十八禁人妻一区二区| 美女 人体艺术 gogo| 麻豆成人午夜福利视频| 制服诱惑二区| 欧美日韩一级在线毛片| tocl精华| 久久久国产成人免费| 视频在线观看一区二区三区| 母亲3免费完整高清在线观看| 最近最新免费中文字幕在线| 国产高清有码在线观看视频 | 中文字幕另类日韩欧美亚洲嫩草| 国产单亲对白刺激| www日本在线高清视频| 美女免费视频网站| 成人午夜高清在线视频 | 两性午夜刺激爽爽歪歪视频在线观看 | 一进一出抽搐动态| 两人在一起打扑克的视频| 人妻久久中文字幕网| 制服人妻中文乱码| 日日爽夜夜爽网站| 又黄又粗又硬又大视频| 亚洲午夜精品一区,二区,三区| 两个人看的免费小视频| 一级毛片女人18水好多| 亚洲国产欧美网| 美国免费a级毛片| 激情在线观看视频在线高清| 亚洲专区国产一区二区| 国产高清视频在线播放一区| 国产精品,欧美在线| 欧美日韩亚洲综合一区二区三区_| 亚洲成av片中文字幕在线观看| 亚洲精品中文字幕在线视频| 久久精品国产亚洲av香蕉五月| 一个人观看的视频www高清免费观看 | 99久久综合精品五月天人人| 看黄色毛片网站| av在线播放免费不卡| 禁无遮挡网站| 欧美最黄视频在线播放免费| 亚洲 欧美 日韩 在线 免费| 欧美黑人欧美精品刺激| 在线观看免费午夜福利视频| 精品卡一卡二卡四卡免费| 亚洲国产毛片av蜜桃av| 久久精品成人免费网站| 欧美久久黑人一区二区| ponron亚洲| 亚洲人成77777在线视频| 热99re8久久精品国产| 欧美乱码精品一区二区三区| 在线视频色国产色| av在线天堂中文字幕| 黄片大片在线免费观看| 国产aⅴ精品一区二区三区波| 黄网站色视频无遮挡免费观看| 亚洲中文av在线| 欧美日韩一级在线毛片| 人人澡人人妻人| av天堂在线播放| 国产亚洲欧美在线一区二区| 日韩高清综合在线| 热99re8久久精品国产| 久久青草综合色| 精品国产亚洲在线| 亚洲精品粉嫩美女一区| 好男人在线观看高清免费视频 | 久久久久久九九精品二区国产 | 国产精品香港三级国产av潘金莲| 美女高潮喷水抽搐中文字幕| 国产高清有码在线观看视频 | 久久人人精品亚洲av| 国产野战对白在线观看| 叶爱在线成人免费视频播放| xxxwww97欧美| 中出人妻视频一区二区| 日韩欧美一区二区三区在线观看| 亚洲熟妇中文字幕五十中出| 亚洲av熟女| 中文字幕人妻熟女乱码| 成人亚洲精品av一区二区| 久久 成人 亚洲| 变态另类成人亚洲欧美熟女| 大型黄色视频在线免费观看| 午夜福利18| 国产av一区二区精品久久| av欧美777| 国产私拍福利视频在线观看| 免费人成视频x8x8入口观看| 制服诱惑二区| 亚洲第一欧美日韩一区二区三区| 亚洲一区中文字幕在线| 亚洲国产精品sss在线观看| 国产精品av久久久久免费| 精品久久久久久久久久久久久 | 亚洲国产中文字幕在线视频| 国产精品电影一区二区三区| 一进一出好大好爽视频| 欧美日韩中文字幕国产精品一区二区三区| 后天国语完整版免费观看| 午夜免费鲁丝| 成人特级黄色片久久久久久久| 国产国语露脸激情在线看| 老司机午夜十八禁免费视频| 久久精品亚洲精品国产色婷小说| 草草在线视频免费看| 在线看三级毛片| 在线观看66精品国产| 19禁男女啪啪无遮挡网站| 最近最新中文字幕大全电影3 | 精品免费久久久久久久清纯| 91九色精品人成在线观看| 长腿黑丝高跟| 久久青草综合色| 国产精品,欧美在线| www.精华液| 久久久国产欧美日韩av| 国产极品粉嫩免费观看在线| 99久久久亚洲精品蜜臀av| 久久久久九九精品影院| 色婷婷久久久亚洲欧美| 国产黄色小视频在线观看| 看片在线看免费视频| 天堂动漫精品| 啦啦啦观看免费观看视频高清| 午夜免费观看网址| 国产高清videossex| 搡老岳熟女国产| 97超级碰碰碰精品色视频在线观看| 国产精品综合久久久久久久免费| 好男人在线观看高清免费视频 | 亚洲人成伊人成综合网2020| www.自偷自拍.com| 国内精品久久久久久久电影| 午夜免费成人在线视频| 很黄的视频免费| 国产爱豆传媒在线观看 | 我的亚洲天堂| 天堂动漫精品| 精品国内亚洲2022精品成人| 亚洲久久久国产精品| 高清在线国产一区| 一进一出好大好爽视频| 最新美女视频免费是黄的| 国产蜜桃级精品一区二区三区| 亚洲av电影不卡..在线观看| 日韩 欧美 亚洲 中文字幕| 性欧美人与动物交配| 亚洲 国产 在线| 男人操女人黄网站| 国产精品综合久久久久久久免费| 国产精品美女特级片免费视频播放器 | 精品日产1卡2卡| 精品久久久久久,| 国内精品久久久久久久电影| 2021天堂中文幕一二区在线观 | 巨乳人妻的诱惑在线观看| aaaaa片日本免费| 大型黄色视频在线免费观看| 非洲黑人性xxxx精品又粗又长| 久久久久久人人人人人| 国产麻豆成人av免费视频| 亚洲成人久久性| 久久性视频一级片| 麻豆久久精品国产亚洲av| 欧美成狂野欧美在线观看| 女同久久另类99精品国产91| 禁无遮挡网站| 12—13女人毛片做爰片一| 好男人在线观看高清免费视频 | 在线av久久热| 久久国产精品男人的天堂亚洲| 欧美黄色淫秽网站| 亚洲中文日韩欧美视频| 一级黄色大片毛片| 免费在线观看完整版高清| 成年人黄色毛片网站| 精品久久久久久久久久免费视频| 亚洲第一av免费看| 精品一区二区三区四区五区乱码| 999久久久精品免费观看国产| a级毛片a级免费在线| 亚洲av电影不卡..在线观看| 亚洲专区中文字幕在线| 中亚洲国语对白在线视频| 人人妻人人看人人澡| 亚洲精品美女久久久久99蜜臀| 特大巨黑吊av在线直播 | 桃红色精品国产亚洲av| 中文字幕精品免费在线观看视频| 国产一卡二卡三卡精品| 日韩一卡2卡3卡4卡2021年| 亚洲最大成人中文| 90打野战视频偷拍视频| 国产在线精品亚洲第一网站| 亚洲欧美日韩无卡精品| 久久性视频一级片| 超碰成人久久| 午夜福利免费观看在线| 制服人妻中文乱码| 日本黄色视频三级网站网址| 欧美日韩亚洲综合一区二区三区_| 欧美色欧美亚洲另类二区| 亚洲三区欧美一区| 琪琪午夜伦伦电影理论片6080| e午夜精品久久久久久久| 国产精品亚洲一级av第二区| 国产私拍福利视频在线观看| 97超级碰碰碰精品色视频在线观看| 欧美日韩精品网址| 亚洲精品国产精品久久久不卡| 人人妻人人澡人人看| 国产精品亚洲一级av第二区| 91成年电影在线观看| 成人18禁高潮啪啪吃奶动态图| 97碰自拍视频| 老鸭窝网址在线观看| 丁香六月欧美| 天天躁夜夜躁狠狠躁躁| 国产成人精品无人区| 精品国内亚洲2022精品成人| 人人妻人人澡欧美一区二区| 制服诱惑二区| 成人一区二区视频在线观看| 一区二区三区高清视频在线| 亚洲一码二码三码区别大吗| √禁漫天堂资源中文www| 99久久久亚洲精品蜜臀av| 中亚洲国语对白在线视频| av电影中文网址| 久久精品国产清高在天天线| x7x7x7水蜜桃| 亚洲色图av天堂| 国产一区二区在线av高清观看| 婷婷丁香在线五月| 人人妻人人澡欧美一区二区| 一本一本综合久久| 一级黄色大片毛片| 日韩欧美免费精品| 日日摸夜夜添夜夜添小说| 99re在线观看精品视频| 99热6这里只有精品| 在线观看66精品国产| 俄罗斯特黄特色一大片| 男女视频在线观看网站免费 | 成人午夜高清在线视频 | 黑丝袜美女国产一区| xxx96com| 男人舔女人下体高潮全视频| 国内精品久久久久久久电影| 欧美三级亚洲精品| 色播亚洲综合网| 香蕉久久夜色| 国产精品久久久人人做人人爽| 欧美色欧美亚洲另类二区| 亚洲av成人不卡在线观看播放网| 亚洲国产精品成人综合色| 成人精品一区二区免费| 99精品欧美一区二区三区四区| 午夜激情福利司机影院| 丁香六月欧美| 国产又色又爽无遮挡免费看| 黄色片一级片一级黄色片| 午夜两性在线视频| 亚洲aⅴ乱码一区二区在线播放 | 91麻豆av在线| 婷婷精品国产亚洲av在线| 在线观看午夜福利视频| 亚洲一区二区三区不卡视频| 久久香蕉国产精品| 国产又黄又爽又无遮挡在线| 亚洲精品一区av在线观看| www.www免费av| 女性被躁到高潮视频| 日韩视频一区二区在线观看| www.熟女人妻精品国产| 久久久国产欧美日韩av| 热99re8久久精品国产| 夜夜看夜夜爽夜夜摸| 露出奶头的视频| 最新美女视频免费是黄的| 色在线成人网| 操出白浆在线播放| 亚洲 欧美 日韩 在线 免费| 亚洲精品久久国产高清桃花| 欧美性长视频在线观看| 啦啦啦 在线观看视频| 亚洲欧洲精品一区二区精品久久久| av中文乱码字幕在线| 亚洲精品久久国产高清桃花| 免费无遮挡裸体视频| 嫩草影院精品99| 色播在线永久视频| 亚洲熟妇熟女久久| 亚洲国产中文字幕在线视频| 亚洲第一av免费看| 亚洲免费av在线视频| 天天躁夜夜躁狠狠躁躁| 在线观看免费日韩欧美大片| 亚洲一区中文字幕在线| 99国产综合亚洲精品| 精品国产国语对白av| 91成年电影在线观看| 一级毛片精品| 深夜精品福利| 精品第一国产精品| 1024香蕉在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲第一av免费看| 首页视频小说图片口味搜索| www.自偷自拍.com| 操出白浆在线播放| 久久婷婷成人综合色麻豆| 一级片免费观看大全| 国产免费男女视频| 国产三级黄色录像| 精品久久久久久久久久免费视频| av欧美777| 久久久久国产一级毛片高清牌| 99国产极品粉嫩在线观看| 久99久视频精品免费| 熟女电影av网| 99久久无色码亚洲精品果冻| 成人永久免费在线观看视频| 露出奶头的视频| 国产激情欧美一区二区| 亚洲久久久国产精品| 最新美女视频免费是黄的| 他把我摸到了高潮在线观看| 美女扒开内裤让男人捅视频| 精品国产一区二区三区四区第35| 成人国语在线视频| 国产男靠女视频免费网站| 日本撒尿小便嘘嘘汇集6| 成人国产一区最新在线观看| 99在线视频只有这里精品首页| 妹子高潮喷水视频| 亚洲aⅴ乱码一区二区在线播放 | netflix在线观看网站| 亚洲人成77777在线视频| 中文字幕高清在线视频| 90打野战视频偷拍视频| 亚洲av片天天在线观看| 久久久精品欧美日韩精品| 久久久久免费精品人妻一区二区 | 久久精品国产99精品国产亚洲性色| 国产野战对白在线观看| av片东京热男人的天堂| av中文乱码字幕在线| 黄色片一级片一级黄色片| 久久精品人妻少妇| 很黄的视频免费| 亚洲自拍偷在线| www国产在线视频色| 亚洲成人精品中文字幕电影| 成年免费大片在线观看| 久久人人精品亚洲av| 免费在线观看亚洲国产| 午夜福利在线在线| 成人国语在线视频| 99re在线观看精品视频| 熟女少妇亚洲综合色aaa.| 男女下面进入的视频免费午夜 | 在线观看66精品国产| 欧美成人一区二区免费高清观看 | 在线观看免费日韩欧美大片| 国产一区二区三区视频了| 国产真实乱freesex| 精品人妻1区二区| 十分钟在线观看高清视频www| 亚洲成人精品中文字幕电影| 色综合欧美亚洲国产小说| 男人舔奶头视频| 精品一区二区三区四区五区乱码| 婷婷精品国产亚洲av在线| 首页视频小说图片口味搜索| 国产又爽黄色视频| 欧洲精品卡2卡3卡4卡5卡区| 色综合亚洲欧美另类图片| 久久国产精品人妻蜜桃| 国产成年人精品一区二区| 久久久久久久久免费视频了| 久久香蕉国产精品| 香蕉久久夜色| a级毛片a级免费在线| 又黄又爽又免费观看的视频| 这个男人来自地球电影免费观看| 亚洲国产日韩欧美精品在线观看 | 2021天堂中文幕一二区在线观 | 久久久久久人人人人人| 欧美日韩福利视频一区二区| 自线自在国产av| 最新在线观看一区二区三区| 日韩免费av在线播放| 国产97色在线日韩免费| 免费看a级黄色片| e午夜精品久久久久久久| 国产又爽黄色视频| 在线播放国产精品三级| 日韩欧美三级三区| 国产精品电影一区二区三区| 成人免费观看视频高清| 亚洲自偷自拍图片 自拍| cao死你这个sao货| 妹子高潮喷水视频| 精品一区二区三区四区五区乱码| 日韩欧美免费精品| 激情在线观看视频在线高清| 成年人黄色毛片网站| 午夜免费成人在线视频| 精品久久蜜臀av无| 久久天堂一区二区三区四区| 黄色a级毛片大全视频| or卡值多少钱| 亚洲国产精品合色在线| 国产精品久久视频播放| 动漫黄色视频在线观看| 国产欧美日韩一区二区精品| 99国产极品粉嫩在线观看| 可以免费在线观看a视频的电影网站| 男人舔女人下体高潮全视频| 国产在线精品亚洲第一网站| 亚洲熟妇熟女久久| 亚洲专区字幕在线| 国产成人影院久久av| 国产精品免费视频内射| 麻豆成人午夜福利视频| 国产男靠女视频免费网站| 狂野欧美激情性xxxx| 2021天堂中文幕一二区在线观 | 无遮挡黄片免费观看| 美女扒开内裤让男人捅视频| 久久天堂一区二区三区四区| 看片在线看免费视频| 露出奶头的视频|