• 
    

    
    

      99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

      步槍彈侵徹明膠的表面受力模型

      2014-03-01 06:54:02莫根林吳志林馮杰劉坤
      兵工學(xué)報(bào) 2014年2期
      關(guān)鍵詞:彈頭明膠質(zhì)心

      莫根林,吳志林,馮杰,劉坤

      (南京理工大學(xué) 機(jī)械工程學(xué)院,江蘇 南京210094)

      0 引言

      彈道明膠被廣泛地用作人體組織的替代物,用于人體創(chuàng)傷的評(píng)估。它具有透明性質(zhì),利用高速攝影可以方便地獲得彈頭在其中的位置和方向信息。實(shí)驗(yàn)表明步槍彈在明膠中的翻滾角度能夠大于90°,此時(shí)小攻角條件下的彈道方程不再適用,需要建立能夠描述彈頭在大攻角條件下運(yùn)動(dòng)的彈道模型[1]。

      文獻(xiàn)[2]中:假設(shè)彈頭阻力和速度呈正比,建立了彈頭的水平運(yùn)動(dòng)方程;假設(shè)彈頭豎直方向的受力和彈頭的轉(zhuǎn)矩僅和彈頭的攻角有關(guān),建立了彈頭的豎直運(yùn)動(dòng)方程和繞質(zhì)心的轉(zhuǎn)動(dòng)方程;根據(jù)實(shí)驗(yàn)數(shù)據(jù)擬合得到了阻力系數(shù)和攻角的關(guān)系。Weinacht 等[3]在文獻(xiàn)[2]的基礎(chǔ)上獲得了彈頭在任意起始攻角下轉(zhuǎn)動(dòng)方程的解析解。Peters 等[4]采用修正的伯努利方程建立了彈頭質(zhì)心的水平運(yùn)動(dòng)方程,在彈頭的平面轉(zhuǎn)動(dòng)方程中引入了動(dòng)力矩系數(shù)和靜力矩系數(shù),建立了力矩系數(shù)和彈頭外形的聯(lián)系。文獻(xiàn)[5 -6]假設(shè)彈頭質(zhì)心的速度方向水平不變、阻力系數(shù)CD和法向力系數(shù)CN為偏角φ 的函數(shù),利用實(shí)驗(yàn)數(shù)據(jù)對(duì)模型假設(shè)進(jìn)行了驗(yàn)證;劉坤等[7]在彈頭的平面彈道模型中區(qū)分了攻角和偏角的概念,從而平面模型中考慮了彈頭豎直方向的運(yùn)動(dòng)。上述建立的模型均為平面彈道模型,沒有考慮彈頭的運(yùn)動(dòng)平面和拍攝平面的關(guān)系,也沒有考慮彈頭轉(zhuǎn)動(dòng)對(duì)侵徹過程的影響。

      為了獲得彈頭運(yùn)動(dòng)規(guī)律,首先需要確定彈頭的受力。在靶板的侵徹動(dòng)力學(xué)中,彈頭表面和侵徹介質(zhì)的相互作用力通常用彈頭的幾何尺寸和動(dòng)力學(xué)參數(shù)表示[8]。將接觸應(yīng)力表達(dá)為微元速度和靶標(biāo)強(qiáng)度的方法分別在鋁板、混凝土靶、巖石靶和鋼板的侵徹實(shí)驗(yàn)中得到了有效驗(yàn)證[9-11]。本文將侵徹明膠彈頭表面的微元應(yīng)力表達(dá)為微元速度和明膠強(qiáng)度的函數(shù),通過步槍彈侵徹明膠的實(shí)驗(yàn),對(duì)該模型進(jìn)行間接驗(yàn)證。

      1 模型

      1.1 現(xiàn)象描述

      現(xiàn)有步槍彈彈頭侵徹明膠過程中會(huì)發(fā)生失穩(wěn)翻滾,步槍彈在30 cm 明膠塊中的位置和姿態(tài)變化示意圖如圖1所示。

      圖1 彈頭在明膠中的運(yùn)動(dòng)過程Fig.1 Penetration process of rifle bullet in gelatin

      從圖中可以看出,彈頭侵徹過程中質(zhì)心在水平方向和豎直方向均有位移變化,彈頭姿態(tài)的變化是空間的而不是平面的,彈頭攻角能夠達(dá)到幾十度。為了描述彈頭在明膠中的空間運(yùn)動(dòng),必須確定彈頭任意時(shí)刻受到的合力和合力矩。

      1.2 模型假設(shè)

      為簡(jiǎn)化模型,提出如下假設(shè):

      1)彈頭為剛性體,不考慮彈頭侵徹過程中的彈塑性變形;

      2)彈頭為旋轉(zhuǎn)體,不考慮其表面劃痕等因素的影響;

      3)彈頭在明膠中的運(yùn)動(dòng)時(shí)間一般小于1 ms,所以不考慮重力影響;

      4)根據(jù)文獻(xiàn)[9 -11],假設(shè)微元面侵徹明膠時(shí)的受力為F,可以表示為

      式中:a、c 為力學(xué)系數(shù);v 為面元的速度;dS 表示面元的面積;n 為面元的外法線方向。

      1.3 幾何和動(dòng)力學(xué)參數(shù)

      旋轉(zhuǎn)體彈頭可以看作是一條曲線繞中心線旋轉(zhuǎn)形成的實(shí)體。由于曲線可以近似看為微小直線段的組合,因此彈頭表面可以看作是很多圓錐面的組合,對(duì)彈頭表面的受力分析可以簡(jiǎn)化為對(duì)錐面的受力分析。

      為了分析錐面上任意位置微元的受力。建立慣性坐標(biāo)系O'x'y'z'和連體坐標(biāo)系Oxyz,如圖2所示。令O'x'y'z'的基矢量為i'、j'和k',Oxyz 的基矢量為i、j 和k.采用歐拉角描述連體坐標(biāo)系相對(duì)慣性坐標(biāo)系的轉(zhuǎn)動(dòng),其中φ、ψ 和θ 分別為自轉(zhuǎn)角、進(jìn)動(dòng)角和章動(dòng)角。θ 為Oz 軸與O'z'的夾角,原點(diǎn)O 位于彈頭的質(zhì)心位置,Oz 與彈頭的對(duì)稱軸重合。zt和zb分別為微元頂面和底面的z 坐標(biāo),R2和R1分別為微元頂面和底面圓面的半徑。錐面上任意z 位置處A 點(diǎn)速度為

      式中:vc為質(zhì)心運(yùn)動(dòng)速度;ω 為彈頭轉(zhuǎn)動(dòng)角速度;rA為A 點(diǎn)相對(duì)矢徑。

      圖2 錐面和連體坐標(biāo)系Oxyz 與固定坐標(biāo)系O'x'y'z'關(guān)系Fig.2 Schematic of a frustum of a cone with body coordinate system Oxyz in fixed coordinate system O'x'y'z'

      轉(zhuǎn)動(dòng)角速度ω 可以表示為

      式中:ωx、ωy和ωz為角速度在x、y 和z 方向的分量,表達(dá)式為

      矢徑r'A可表示為

      式中:r'c為彈頭質(zhì)心的矢徑;rA可以表示為

      ρ 和δ 分別為rA在Oxy 平面的投影向量的模和投影向量與x 軸的夾角。

      采用坐標(biāo)z 和δ 表示微元的面積dS,可得

      A 點(diǎn)的外法線方向n 可以表示為

      式中:

      將(2)式~(9)式帶入(1)式即可求得任意位置點(diǎn)出的受力F.將F 對(duì)質(zhì)心O 點(diǎn)取矩得

      1.4 運(yùn)動(dòng)方程

      根據(jù)高等動(dòng)力學(xué)理論[12],彈頭的質(zhì)心運(yùn)動(dòng)方程可以表示為

      式中:m 表示彈頭質(zhì)量;x'c、y'c和z'c表示彈頭質(zhì)心在固定坐標(biāo)系中的坐標(biāo);∑Fx'、∑Fy'和∑Fz'為彈頭在i'、j'和k'方向所受的合力。

      彈頭的歐拉運(yùn)動(dòng)方程可以表示為

      式中:Jx、Jy和Jz分別為彈頭繞x 軸、y 軸和z 軸的轉(zhuǎn)動(dòng)慣量;∑Mx、∑My和∑Mz為彈頭繞x 軸、y軸和z 軸的力矩分量。

      對(duì)(4)式、(11)式和(12)式整理可得x'c、y'c、z'c、φ、ψ 和θ 表示的剛體的空間方程為

      式中:k1、k2和k3分別為

      利用數(shù)值積分,將(1)式和(10)式得到的∑Fx'、∑Fy'、∑Fz'、∑Mx、∑My和∑Mz帶入(13)式和(14)式。在初始條件已知的情況下,通過4 階龍格-庫(kù)塔法即可得到彈頭位置和姿態(tài)的變化規(guī)律。

      2 實(shí)驗(yàn)驗(yàn)證

      實(shí)驗(yàn)中選用10%的彈道明膠,明膠塊尺寸為30 cm×30 cm×30 cm,明膠塊在使用前在4 ℃保溫箱中保溫24 h.彈道槍距離明膠塊50 m 發(fā)射7.62 mm步槍彈侵徹明膠塊,通過光電靶測(cè)量明膠前1 m 處彈頭速度,采用1.5 MHz 幀頻的高速攝影從側(cè)面拍攝彈頭在明膠中的運(yùn)動(dòng)姿態(tài)。圖3為高速攝影和明膠塊在固定坐標(biāo)系中相對(duì)位置示意圖。圖4為彈頭在明膠出口位置的實(shí)驗(yàn)圖像,可以看出彈頭在出口時(shí)幾乎沒有塑性變形,因此可以將侵徹過程中的彈頭看作剛性體。

      圖3 高速攝像機(jī)和明膠塊在固定坐標(biāo)系O'x'y'z'中的位置Fig.3 Positions of the high-speed camera and gelatin block in O'x'y'z'

      圖4 明膠塊出口處的7.62 mm 彈頭Fig.4 7.62 mm bullet at exit of the gelatin block

      模型中的初始參數(shù)如表1和表2所示,下標(biāo)“0”表示相應(yīng)物理量的初始值,其中彈頭的角速度分量通過理論計(jì)算得到[13]。由于高速攝影獲得的平面圖像無法確定三維空間的歐拉角,所以初始?xì)W拉角由試算得到。實(shí)驗(yàn)表明空氣中彈頭歐拉角的取值必須保證彈頭的攻角比較小[13-14],因此彈頭侵徹明膠的初始?xì)W拉角必須是一個(gè)較小值。根據(jù)表1和表2可以得到彈頭攻角為3.5°,可以認(rèn)為給定的初始值是合理的。明膠的力學(xué)參數(shù)根據(jù)鋼球侵徹明膠的實(shí)驗(yàn)數(shù)據(jù)擬合得到[15]。

      表1 7.62 mm 彈頭的初始位移和速度分量Tab.1 Initial displacements and initial velocity components of 7.62 mm projectile

      表2 7.62 mm 彈頭的初始?xì)W拉角和初始角速度Tab.2 Initial Eulerian angles and initial angular velocities of 7.62 mm projectile

      理論計(jì)算結(jié)果和實(shí)驗(yàn)值的比較如圖5~圖7所示。圖5為彈頭質(zhì)心z'c坐標(biāo)隨時(shí)間的變化,最大實(shí)驗(yàn)誤差7.8 mm,平均實(shí)驗(yàn)誤差3.1 mm;圖6為y'c隨z'c的變化,最大實(shí)驗(yàn)誤差2.3 mm,平均誤差2.1 mm;圖7為彈頭在O'y'z'平面投影轉(zhuǎn)角隨z'c坐標(biāo)的變化,最大實(shí)驗(yàn)誤差8.9°,平均誤差2.3°.可見,實(shí)驗(yàn)值和理論值之間的誤差較小。

      圖5 z'c與t 的關(guān)系Fig.5 z'c vs.t

      圖6 y'c與z'c的關(guān)系Fig.6 y'c vs.z'c

      圖7 O'y'z'平面投影轉(zhuǎn)角φx'與z'c的關(guān)系Fig.7 Projection angle φx' vs.z'c on O'y'z'

      3 討論分析

      3.1 章動(dòng)角和平面投影轉(zhuǎn)角

      理論計(jì)算得到的彈頭章動(dòng)角θ 和彈頭在O'y'z'平面投影轉(zhuǎn)角φx'的關(guān)系如圖8所示。它們的關(guān)系可表示為

      彈頭的平面運(yùn)動(dòng)模型中,若假設(shè)彈頭在高速攝影拍攝平面內(nèi)運(yùn)動(dòng)[4-7],此時(shí)其進(jìn)動(dòng)角為0°,從(15)式可知,θ 和φx'的大小相同;在空間運(yùn)動(dòng)中,當(dāng)進(jìn)動(dòng)角變化或不為0°時(shí),θ 和φx'將有很大的不同。因此平面模型的實(shí)驗(yàn)驗(yàn)證需要確定彈頭的運(yùn)動(dòng)平面。

      圖8 投影角φx'和θ 的比較Fig.8 Comparison between φx' and θ

      彈頭外形在彈尖和彈尾處變化比較劇烈,與該處接觸的明膠必然處于應(yīng)力集中狀態(tài)。明膠在軸對(duì)稱膨脹過程中,該處也就較容易產(chǎn)生徑向裂紋。顯然裂紋方向的變化是由彈頭的進(jìn)動(dòng)角變化引起的,所以彈頭在明膠中進(jìn)動(dòng)角的變化可以從明膠中的裂紋方向看出。圖9為第2 節(jié)中7.62 mm 彈頭進(jìn)動(dòng)角(理論值)隨侵徹距離的變化。圖10為明膠中裂紋方向變化形成的裂紋曲面,和圖9中的曲線變化相似,明膠中的裂紋方向在0.25 m 前基本不變,在0.25 m 后發(fā)生劇烈變化。圖10進(jìn)一步表明了模型假設(shè)的合理性。

      圖9 進(jìn)動(dòng)角隨z'c的變化Fig.9 ψ vs.z'c

      3.2 力學(xué)系數(shù)a 和c

      圖10 明膠中的裂紋曲面Fig.10 Curved surface of cracks in gelatin block

      明膠是一種應(yīng)變率相關(guān)材料。Kwon 等[16]利用單軸壓縮實(shí)驗(yàn)獲得了10%明膠應(yīng)變率在2 000 ~3 200 s-1之間時(shí),應(yīng)力與應(yīng)變的實(shí)驗(yàn)數(shù)據(jù);Cronin 等[17]通過單軸壓縮實(shí)驗(yàn)得到了應(yīng)變率為0.01 s-1、10 s-1、117 s-1和1 552 s-1時(shí),應(yīng)力與應(yīng)變的實(shí)驗(yàn)數(shù)據(jù);Moy 等[18]采用拉伸實(shí)驗(yàn)得到了應(yīng)變率為0.001 s-1、0.01 s-1和1 s-1時(shí),應(yīng)力與應(yīng)變的實(shí)驗(yàn)數(shù)據(jù)。上述實(shí)驗(yàn)均表明明膠的力學(xué)性質(zhì)與應(yīng)變率具有較大相關(guān)性。對(duì)于明膠的本構(gòu)模型,文獻(xiàn)[19 -20]建立了準(zhǔn)靜態(tài)條件下的明膠超彈性本構(gòu)模型和粘彈性固體模型。

      Shepherd 等[21]通過實(shí)驗(yàn)獲得了沖擊速度75 ~860 m/s 時(shí),20% 明膠的雨貢紐曲線。文獻(xiàn)[16 -21]中的實(shí)驗(yàn)數(shù)據(jù)是在明膠不存在剪切破壞的條件下獲得的。在破片侵徹實(shí)驗(yàn)中可以觀察到彈道周圍的明膠存在明顯的剪切破壞,這種破壞對(duì)侵徹過程的影響程度尚不清楚。因此不能直接從上述文獻(xiàn)中得到10%明膠的力學(xué)參數(shù)a 和c.

      Mo 等[15]采用本文的面元積分方法,對(duì)直徑3 mm、4 mm 和4.8 mm 球形破片侵徹明膠的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了擬合,間接確定了模型中的力學(xué)系數(shù)。本文力學(xué)系數(shù)a 和c 的取值與文獻(xiàn)[15]相同。實(shí)驗(yàn)中破片撞擊速度在600 ~800 m/s 范圍。對(duì)于低速100 ~300 m/s 情況下,明膠的力學(xué)系數(shù)是否保持不變;明膠力學(xué)系數(shù)是否和破片外形的關(guān)系都有待進(jìn)一步研究。將在實(shí)驗(yàn)中增加對(duì)柱形和菱形破片的侵徹研究,并采用不同發(fā)射速度研究速度對(duì)力學(xué)系數(shù)的影響。

      4 結(jié)論

      本文在彈頭表面受力模型的基礎(chǔ)上,分析了旋轉(zhuǎn)體彈頭的受力,運(yùn)用4 階龍格-庫(kù)塔法通過求解剛體的空間運(yùn)動(dòng)方程,得到了7.62 mm 步槍彈侵徹明膠過程的質(zhì)心位置、空間姿態(tài)變化規(guī)律。理論計(jì)算值和實(shí)驗(yàn)值誤差較小,從而間接證明了表面受力模型的合理性。模型中力學(xué)系數(shù)是從球形破片侵徹明膠的實(shí)驗(yàn)中間接得到的,它們隨沖擊速度、破片外形的影響有待進(jìn)一步研究。

      References)

      [1] Kneubuehl B P,Coupland R M,Rothschild M A,et al.Wound ballistics:basics and applications[M].Berlin:Springer,2011:95 -100.

      [2] Roecker E T.On yawing motion,velocity decay and path deflection of projectile penetrating cavitating media,ADA168220[R].Maryland:US Army Ballistic Research Laboratory Aberdeen Proving Ground,1986.

      [3] Weinacht P,Cooper G R.Analytical solutions for the linear and nonlinear yawing motion in dense media[C]∥23rd International Symposiumon Ballistics.Tarragona, Spain: PAP Congresos,2007:1339 -1346.

      [4] Peters C E,Sebourn C L,Crowder H L.Wound ballistics of unstable projectiles.part I:projectile yaw growth and retardation[J].Journal of Trauma,1996,40(3S):10S-15S.

      [5] Kneubuehl B,Sellier K.Woundballistics:a new understanding of the behavior of a bullet in a dense medium[C]∥Proceedings 13th International Symposium of Ballistics.Stockholm:IBC,1992:445 - 462.

      [6] Sellier K,Kneubuehl B P.Woundballlistics and the scientific background[M].New York:Elsevier,1994:154 -168.

      [7] 劉坤,吳志林,徐萬和,等.彈頭侵徹明膠的運(yùn)動(dòng)模型[J].爆炸與沖擊,2012,32(6):616 -622.LIU Kun,WU Zhi-lin,XU Wan-he,et al.A motion model for bullet penetrating gelatin[J].Explosion and Shock Waves,2012,32(6):616 -622.(in Chinese)

      [8] Ben-Dor G,Dubinsky A,Elperin T.Ballistic impact:recent advances in analytical modeling of plate penetration dynamics-a review[J].Applied Mechanics Reviews,2005,58(6):355-371.

      [9] Warren T L,Tabbara M R.Spherical cavity-expansion forcing function in PRONTO 3D for application to penetration problems,SAND-97-1174[R].Albuquerque,NM:Sandia National Lab,1997.

      [10] Warren T L.Simulations of the penetration of limestone targets by ogive-nose 4340 steel projectiles[J].International Journal of Impact Engineering,2002,27(5):475 -496.

      [11] Satapathy S,Bless S.Response of long rods to moving lateral pressure pulse:numerical evaluation in DYNA3D[J].International Journal of Impact Engineering,2001,26(1):663 -674.

      [12] 張延教.高等動(dòng)力學(xué)[M].南京:華東工程學(xué)院,1984:56 -280.ZHANG Yan-jiao.Advanced dynamics[M].Nanjing:East China Institute of Engineering,1984:56 -280.(in Chinese)

      [13] 高乃同.自動(dòng)武器彈藥學(xué)[M].北京:國(guó)防工業(yè)出版社,1990:334 -336.GAO Nai-tong.Automatic weapons and ammunition[M].Beijing:National Defense Industry Press,1990:334 -336.(in Chinese)

      [14] Knudsen P J T,S?rensen O H.The initial yaw of some commonly encountered military rifle bullets[J].Int J Legal Med,1994,107(3):141 -146.

      [15] Mo G L,Wu Z L,Liu K.Drag force analysis of spheres penetrating gelatin based on surface integral[J].Applied Mechanics and Materials,2012,105:561 -565.

      [16] Kwon J,Subhash G.Compressive strain rate sensitivity of ballistic gelatin[J].Journal of Biomechanics,2010,43(3):420 -425.

      [17] Cronin D S,F(xiàn)alzon C.Dynamic characterization and simulation of ballistic gelatin[C]∥2009 SEM Conference & Exposition on Experimental & Applied Mechanics.Albuquerque,NM:SEM,2009:1 -4.

      [18] Moy P,Tusit W T,Gunnarsson C L.Tensile deformation of ballistic gelatin as a function of loading rate[C]∥Proceedings of the 6th International Congress and Exposition.Orlando,F(xiàn)L:SEM,2008:2 -5.

      [19] Ottone M L,Deiber J A.Modeling the rheology of gelatin gels for finite deformations.part 1.elastic rheological model[J].Polymer,2005,46(13):4928 -4937.

      [20] Ottone M L,Peirotti M B,Deiber J A.Modeling the rheology of gelatin gels for finite deformations.part 2.viscoelastic solid model[J].Polymer,2005,46(13):4938 -4949.

      [21] Shephert C J,Appleby-Thomas G J,Hazell P J,et al.The dynamic behavior of ballistic gelatin[J].AIP Conference Proceedings,2009,1195:1399 -1402.

      猜你喜歡
      彈頭明膠質(zhì)心
      熱身管下小口徑槍彈彈頭殼材料對(duì)其膛內(nèi)運(yùn)動(dòng)的影響規(guī)律
      重型半掛汽車質(zhì)量與質(zhì)心位置估計(jì)
      基于GNSS測(cè)量的天宮二號(hào)質(zhì)心確定
      導(dǎo)彈突防概率解析與仿真計(jì)算結(jié)果差異分析
      超高壓明膠理化性質(zhì)分析
      一種海洋測(cè)高衛(wèi)星質(zhì)心在軌估計(jì)算法
      航天器工程(2014年5期)2014-03-11 16:35:53
      復(fù)凝法制備明膠微球
      河南科技(2014年22期)2014-02-27 14:18:07
      宋海新:“問題”明膠生產(chǎn)者
      法人(2014年2期)2014-02-27 10:42:17
      淺談明膠
      均質(zhì)半圓盤質(zhì)心計(jì)算的微元選取及討論
      物理與工程(2010年1期)2010-03-25 10:01:48
      新巴尔虎左旗| 遂平县| 体育| 湖州市| 亳州市| 龙陵县| 理塘县| 天门市| 思南县| 斗六市| 清镇市| 长武县| 青铜峡市| 黄大仙区| 岗巴县| 潮安县| 肇庆市| 古丈县| 诸城市| 姜堰市| 五原县| 武邑县| 民勤县| 乐业县| 米泉市| 南川市| 濉溪县| 曲周县| 汶川县| 大安市| 纳雍县| 响水县| 昌江| 陆川县| 边坝县| 武夷山市| 北川| 上饶市| 白水县| 白河县| 盐池县|