湯成莉，劉蘇蘇，楊玉萍，房佳恒

豬皮對(duì)球形彈丸阻滯作用的仿真研究

湯成莉，劉蘇蘇，楊玉萍，房佳恒

(南通大學(xué)機(jī)械工程學(xué)院，江蘇 南通 226019)

為了明確皮膚對(duì)球形彈丸的阻滯作用機(jī)理，建立了模擬皮膚的3層結(jié)構(gòu)模型。基于拉伸及壓縮力學(xué)特性曲線，利用MATLAB對(duì)皮膚模型進(jìn)行參數(shù)識(shí)別?；跊_擊實(shí)驗(yàn)得到的皮膚穿透閾值速度驗(yàn)證了皮膚模型的準(zhǔn)確性。利用LS-DYNA仿真模擬了4.5 mm鉛球侵徹皮膚過程，從球形彈丸速度衰減、皮膚鼓包方面分析了皮膚對(duì)球形彈丸的阻滯作用。仿真模擬了相同動(dòng)能不同密度條件下，直徑為4.5 mm的彈丸侵徹皮膚的過程和相同動(dòng)能不同直徑條件下，鉛球侵徹皮膚的過程，結(jié)果表明，建立的皮膚3層結(jié)構(gòu)模型可以用于模擬皮膚的真實(shí)傷害，作為皮膚的數(shù)值仿真模型；在初始動(dòng)能相同的情況下，球形彈丸密度越小，體積越大，能量耗散比越大，皮膚阻滯作用越明顯，為后續(xù)研究投射物對(duì)真實(shí)人體的傷害奠定了基礎(chǔ)。

球形彈丸；阻滯作用；速度衰減；能量耗散；皮膚

目前，創(chuàng)傷彈道學(xué)研究中主要采用肌肉仿生靶標(biāo)或肌肉-骨骼復(fù)合仿生靶標(biāo)。溫垚珂等[1]利用明膠代替肌肉，分析了明膠應(yīng)力應(yīng)變場(chǎng)、壓力波和能量轉(zhuǎn)換的規(guī)律。黃玲等[2]建立了9 mm手槍彈侵徹肌肉-骨骼復(fù)合仿真模型，分析了彈丸速度、不同章動(dòng)角等因素對(duì)復(fù)合模型致傷效應(yīng)的影響。ZHANG等[3]建立了7.62 mm步槍彈侵徹明膠-骨骼復(fù)合仿真模型，通過數(shù)值模擬分析了空腔膨脹情況對(duì)骨骼間接損傷的影響。皮膚對(duì)投射物阻滯作用的研究較少。DIMAIO等[4]通過向人體下肢皮膚射擊不同口徑的子彈進(jìn)行皮膚阻滯作用實(shí)驗(yàn)，得到了3種子彈侵徹皮膚的穿透閾值速度。JUSSILA等[5]計(jì)算并修正了4.5 mm鉛球侵徹皮膚的穿透閾值速度，并通過實(shí)驗(yàn)進(jìn)行驗(yàn)證。BREEZE和CLASPER[6]對(duì)投射物穿透皮膚的穿透閾值速度進(jìn)行綜述性研究，提出“未來的皮膚仿真材料模型必須與穿透閾值速度相匹配”。

數(shù)值模擬能夠直觀地觀察更多細(xì)節(jié)，比實(shí)驗(yàn)更有助于闡明皮膚阻滯作用。然而皮膚有限元模型的應(yīng)用主要停留在臨床醫(yī)學(xué)方向，普遍采用超彈性模型作為皮膚的材料模型。FLYNN等[7]采用Mooney-Rivlin超彈性模型模擬人體皮膚，建立了體內(nèi)張力作用下的人臉有限元模型，用于模擬人臉面部表情。SWAIN和GUPTA[8]使用Varga超彈性模型建立了皮膚傷口愈合階段的生物力學(xué)模型，強(qiáng)調(diào)在愈合過程中皮膚會(huì)出現(xiàn)應(yīng)力和皺紋。LI和LUO[9]基于Odgen超彈性模型，提出了一種新的皮膚損傷模型。在投射物侵徹皮膚過程中，其產(chǎn)生大應(yīng)變行為，表現(xiàn)出應(yīng)變硬化現(xiàn)象[10]；投射物完全穿透皮膚后，其產(chǎn)生塑性變形，留下彈道傷口，而單一的超彈性模型不能描述沖擊過程中皮膚的應(yīng)變硬化現(xiàn)象和彈道傷口。本文建立了可以應(yīng)用于沖擊研究的皮膚3層結(jié)構(gòu)模型和球形彈丸侵徹皮膚有限元模型的過程，為后續(xù)建立更真實(shí)的肢體(皮膚、肌肉、骨骼、血管等)復(fù)合仿真模型奠定了基礎(chǔ)。

1 皮膚模型建立及參數(shù)識(shí)別

1.1 皮膚模型建立

皮膚分為表皮層、真皮層和皮下組織層[11]。真皮層的厚度決定了皮膚的厚度。膠原纖維是真皮層可以經(jīng)受大應(yīng)力、大變形的主要原因，但膠原纖維只占真皮的60%~80%[12]。本文設(shè)想皮膚為3層：上下兩層采用雙線性各向同性硬化模型(bilinear isotropic hardening constants)，可以描述皮膚的大應(yīng)變行為和在大應(yīng)變情況下的硬化現(xiàn)象，占皮膚材料模型的40%；中間一層采用Mooney-Rivlin超彈性模型，可以描述材料在30%~200%的應(yīng)變行為，占皮膚材料模型的60%，皮膚的3層材料模型如圖1所示。

圖1 皮膚材料模型

Mooney-Rivlin超彈性模型的應(yīng)變勢(shì)能為

初始d和Mooney-Rivlin模型的楊氏模量1可表示為

其中，10，01，d為由實(shí)驗(yàn)確定的材料常數(shù)；1，2為Green張量不變量；為泊松比。

1.2 參數(shù)識(shí)別

ANKERSEN等[13]研究了豬腹部皮膚比背部皮膚更貼近于人腹部皮膚；SHERGOLD等[14]通過對(duì)比發(fā)現(xiàn)豬腹部皮膚(平行于脊柱方向)與人體腹部皮膚的力學(xué)特性曲線相似。因此，本文選取健康長(zhǎng)白豬的腹部皮膚(平行于脊柱)作為實(shí)驗(yàn)材料，用手術(shù)刀小心地剔除皮下脂肪組織，并保持豬皮的平滑性，通過多次測(cè)量取平均值得到豬皮厚度為2.5 mm。分別對(duì)長(zhǎng)方體豬皮(長(zhǎng)×寬：40 mm× 15 mm)和圓柱狀豬皮(20 mm)進(jìn)行拉伸及壓縮實(shí)驗(yàn)，得到皮膚在0.10/s和0.01/s應(yīng)變率下的拉伸及壓縮力學(xué)特性曲線，如圖2(a)所示，求得BISO模型所需材料參數(shù)2(楊氏模量，彈性區(qū)間的斜率)、σ(屈服應(yīng)力，應(yīng)力-應(yīng)變曲線的彈性區(qū)域及塑性區(qū)域曲線外推的交點(diǎn))、tan(切線模量，屈服極限和強(qiáng)度極限之間的斜率)。結(jié)合壓縮力學(xué)特性曲線，如圖2(b)所示，運(yùn)用MATLAB對(duì)曲線進(jìn)行基于最小二乘法的非線性擬合得出Mooney-Rivlin模型的參數(shù)10和01，皮膚的材料參數(shù)見表1和表2。

圖2 不同應(yīng)變率下皮膚力學(xué)特性曲線及參數(shù)識(shí)別

表1 BISO模型參數(shù)

表2 Mooney-Rivlin模型參數(shù)

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

2.1 沖擊實(shí)驗(yàn)

利用9 mm一級(jí)輕氣炮進(jìn)行皮膚沖擊實(shí)驗(yàn)，炮管中用彈托固定4.5 mm鉛球。選取健康長(zhǎng)白豬的腹部皮膚，其彈道沖擊的有效尺寸為100.0 mm× 100.0 mm×2.5 mm，利用夾具將豬皮固定在靶室中(圖3)。通過測(cè)速光幕靶來測(cè)量鉛球速度、數(shù)據(jù)采集處理系統(tǒng)來記錄鉛球速度，如圖4所示。經(jīng)過多次沖擊實(shí)驗(yàn)統(tǒng)計(jì)后得到皮膚的穿透閾值速度約為94 m/s。

圖3 沖擊實(shí)驗(yàn)

圖4 實(shí)驗(yàn)裝置示意圖

2.2 結(jié)果分析

在建立4.5 mm鉛球侵徹皮膚(長(zhǎng)×寬×高為100.0 mm×100.0 mm×2.5 mm)有限元模型時(shí)，鉛球采用JOHNSON_COOK模型，材料參數(shù)見表3。彈丸與皮膚之間采用面面侵蝕接觸(*CONTACT_ ERODING_SURFACE_TO_SURFACE)。對(duì)彈丸賦予垂直于皮膚表面的初速度TH為94 m/s，對(duì)皮膚模型的邊緣面給與固定約束。根據(jù)有限元網(wǎng)格劃分的一般原則，對(duì)彈丸侵徹皮膚的附近地區(qū)進(jìn)行網(wǎng)格加密，保證仿真結(jié)果的精度和可信度。而對(duì)距侵徹中心較遠(yuǎn)的地區(qū)的網(wǎng)格進(jìn)行放大，以減小網(wǎng)格數(shù)來提高計(jì)算速度。對(duì)于皮膚模型，將半徑為20 mm的侵徹中心圓分為40段，侵徹中心圓以外分為逐漸稀疏的40段，厚度方向分為25段，如圖5(a)所示。彈丸沿半徑方向分別劃分為8段，通過設(shè)置粘性的沙漏系數(shù)為0.01，避免沙漏對(duì)數(shù)值仿真結(jié)果可靠性的影響。建立的彈丸侵徹皮膚有限元模型如圖5(b)所示。

表3 鉛球材料模型

仿真模擬了初速度為94 m/s的鉛球侵徹皮膚過程，得到彈丸速度衰減曲線，如圖6所示。與實(shí)驗(yàn)數(shù)據(jù)對(duì)比發(fā)現(xiàn)，仿真結(jié)果與實(shí)驗(yàn)結(jié)果基本吻合，驗(yàn)證了皮膚模型的正確性。

圖5 彈丸侵徹皮膚有限元模型

圖6 鉛球速度衰減曲線

3 皮膚阻滯作用機(jī)理研究

3.1 侵徹過程分析

由圖7可以看出，隨著鉛球的逐漸深入，由于鉛球大部分能量傳遞給皮膚，使得最大應(yīng)力集中在鉛球侵徹的前方。在110 μs時(shí)，皮膚所受應(yīng)力超過皮膚第1層的抗拉強(qiáng)度，達(dá)到第1個(gè)應(yīng)力峰值，皮膚第1層斷裂；在240 μs時(shí)，皮膚所受Von Mises應(yīng)力達(dá)到最大值，皮膚3層全部斷裂，最大應(yīng)力開始集中在傷口處；在500 μs時(shí)，彈丸侵徹完成，留下星型傷口。

結(jié)合圖6速度曲線圖可發(fā)現(xiàn)，在0~90 μs之間，鉛球處于侵徹初期，速度快速衰減。在90 μs(皮膚第1層開始斷裂)時(shí)，速度進(jìn)入衰減緩慢期，皮膚鼓包的直徑只有14 mm。在110 μs時(shí)，鉛球侵徹皮膚的第2層，由于皮膚第2層沒有發(fā)生嚴(yán)重變形及第3層幾乎沒有變形，速度衰減慢。直到120 μs時(shí)，皮膚第2，3層發(fā)生明顯變形，速度進(jìn)入快速衰減階段。在110~120 μs內(nèi)，速度衰減較慢，皮膚鼓包直徑幾乎沒有變化，近乎為16.7 mm。隨著鉛球深入侵徹皮膚模型，在240 μs時(shí)皮膚3層全部斷裂，皮膚形成的鼓包范圍增大至39.1 mm。240 μs后鉛球速度衰減降低直至500 μs停止衰減，此時(shí)皮膚鼓包直徑為58.2 mm。

圖7 典型時(shí)刻鉛球侵徹皮膚過程圖

3.2 彈丸不同因素對(duì)皮膚阻滯作用的影響

為了研究彈丸不同因素對(duì)皮膚阻滯作用的影響規(guī)律，仿真模擬了相同動(dòng)能不同密度下4.5 mm鉛球、鋼球和鋁球侵徹皮膚的過程和相同動(dòng)能不同直徑鉛球侵徹皮膚的過程。通過仿真模擬得到彈丸能量耗散比(耗散能量與初始能量之間的比值)與初始能量之間的定量關(guān)系曲線，如圖8所示。隨著初始動(dòng)能的增加，彈丸能量耗散比快速降低，但是皮膚的阻滯作用是有限的，隨著沖擊動(dòng)能的增大，其耗散比接近于一個(gè)定值。在初始動(dòng)能相同的情況下，球形彈丸密度越小，尺寸越大，能量耗散比越大，受到的皮膚阻滯作用越明顯，如當(dāng)初始動(dòng)能為32.517 J時(shí)，相同大小的鉛球能量耗散比為8.76%，鋼球?yàn)?7.80%，鋁球?yàn)?4.16%；當(dāng)初始動(dòng)能為10.619 J，相同密度條件下，直徑為4.5 mm鉛球能量耗散比為20.19%，直徑為5.0 mm鉛球?yàn)?1.95%，直徑為6.0 mm鉛球?yàn)?1.33%。

圖8 彈丸能量耗散比曲線

4 結(jié)束語

(1) 本文建立了應(yīng)用于沖擊的皮膚3層結(jié)構(gòu)模型和球形彈丸侵徹皮膚有限元模型，并通過沖擊實(shí)驗(yàn)驗(yàn)證了皮膚模型的正確性。

(2) 基于4.5 mm鉛球侵徹皮膚有限元模型，發(fā)現(xiàn)鉛球速度在皮膚第1層結(jié)構(gòu)破裂時(shí)由快速衰減逐漸轉(zhuǎn)變?yōu)榫徛p，在皮膚第2，3層結(jié)構(gòu)發(fā)生大變形后再次進(jìn)入快速衰減階段，在皮膚3層結(jié)構(gòu)全部破裂后衰減越來越緩慢直至侵徹完成時(shí)停止；皮膚鼓包范圍從皮膚鼓包出現(xiàn)時(shí)增大，直到侵徹完成時(shí)停止，皮膚留下星型傷口。

(3) 模擬了相同動(dòng)能不同密度下4.5 mm鉛球、鋼球和鋁球侵徹皮膚的過程和相同動(dòng)能不同直徑鉛球侵徹皮膚的過程，得出在初始動(dòng)能相同的情況下，球形彈丸密度越小，體積越大，能量耗散比越大，皮膚阻滯作用越明顯。

[1] 溫垚珂, 徐誠, 陳愛軍, 等. 球形破片高速侵徹明膠靶標(biāo)的數(shù)值模擬[J]. 彈道學(xué)報(bào), 2012, 24(3): 25-30.

[2] 黃玲, 劉艷, 徐誠, 等. 手槍彈對(duì)骨骼間接作用的數(shù)值模擬及致傷力學(xué)響應(yīng)研究[J]. 彈道學(xué)報(bào), 2013, 25(2): 95-99.

[3] ZHANG X Y, XU C, WEN Y K, et al. The experimental and numerical study of indirect effect of a rifle bullet on the bone[J]. Forensic Science International, 2015, 257: 473-480.

[4] DIMAIO V J M, COPELAND A R, BESANT- MATTHEWS P E, et al. Minimal velocities necessary for perforation of skin by air gun pellets and bullets[J]. Journal of Forensic Sciences, 1982, 27(4): 894-898.

[5] JUSSILA J, LEPP?NIEMI A, PARONEN M, et al. Ballistic skin simulant[J]. Forensic Science International, 2005, 150(1): 63-71.

[6] BREEZE J, CLASPER J C. Determining the velocity required for skin perforation by fragment simulating projectiles: a systematic review[J]. Journal of the Royal Army Medical Corps, 2013, 159(4): 265-270.

[7] FLYNN C, STAVNESS I, LLOYD J, et al. A finite element model of the face including an orthotropic skin model underin vivotension[J]. Computer Methods in Biomechanics and Biomedical Engineering, 2015, 18(6): 571-582.

[8] SWAIN D, GUPTA A. Interfacial growth during closure of a cutaneous wound: stress generation and wrinkle formation[J]. Soft Matter, 2015, 11(32): 6499-6508.

[9] LI W G, LUO X Y. An invariant-based damage model for human and animal skins[J]. Annals of Biomedical Engineering, 2016, 44(10): 3109-3122.

[10] NICOLLE S, DECORPS J, FROMY B, et al. New regime in the mechanical behavior of skin: strain-softening occurring before strain-hardening[J]. Journal of the Mechanical Behavior of Biomedical Materials, 2017, 69: 98-106.

[11] WILDNAUER R H, BOTHWELL J W, DOUGLASS A B. Stratum corneum biomechanical properties I. influence of relative humidity on normal and extracted human stratum corneum[J]. Journal of Investigative Dermatology, 1971, 56(1): 72-78.

[12] 盧天健, 徐峰. 皮膚的力學(xué)性能概述[J]. 力學(xué)進(jìn)展, 2008, 38(4): 393-426.

[13] ANKERSEN J, BIRKBECK A E, THOMSON R D, et al. Puncture resistance and tensile strength of skin simulants[J]. Proceedings of the Institution of Mechanical Engineers, Part H: Journal of Engineering in Medicine, 1999, 213(6): 493-501.

[14] SHERGOLD O A, FLECK N A, RADFORD D. The uniaxial stress versus strain response of pig skin and silicone rubber at low and high strain rates[J]. International Journal of Impact Engineering, 2006, 32(9): 1384-1402.

Simulation study on the blocking effect of pig skin on spherical projectile

TANG Cheng-li, LIU Su-su, YANG Yu-ping, FANG Jia-heng

(School of Mechanical Engineering, Nantong University, Nantong Jiangsu 226019, China)

In order to clarify the blocking mechanism of skin on spherical bullet, a three-layer structural model of simulated skin was established. Based on the mechanical curves of tensile and compression experiments, the parameters of the skin model were identified by MATLAB. Based on the skin penetration threshold velocity obtained from the impact experiment, the accuracy of the skin model was verified. Then, the penetration process of4.5 mm shot into the skin was simulated by LS-DYNA simulation, and the blocking effect of skin on spherical projectile was analyzed from the aspects of velocity attenuation and skin bulging. The penetration process of4.5 mm projectile into the skin under the condition of the same kinetic energy and different density, together with the process of shot put into the skin with the same kinetic energy and different diameter, was numerically simulated. The results show that the three-layer structural model of skin can be taken as a numerical simulation model to simulate the real damage of skin. Furthermore, with the same initial kinetic energy, the smaller the density of spherical projectile is, the larger the volume is, the larger the energy dissipation ratio is, and the more obvious the skin blocking effect is. Such finding slay a foundation for the future study of the harm of the projectile to the real human body.

spherical bullet; blocking effect; velocity attenuation; energy dissipation; skin

TP 391

10.11996/JG.j.2095-302X.2020020283

A

2095-302X(2020)02-0283-05

2019-10-09；

2019-11-22

南通市2018-基礎(chǔ)科學(xué)研究項(xiàng)目(JC2018098)

湯成莉(1995-)，女，安徽六安人，碩士研究生。主要研究方向?yàn)閯?chuàng)傷彈道學(xué)。E-mail：1109584472@qq.com

劉蘇蘇(1988-)，男，江蘇南通人，講師，博士，碩士生導(dǎo)師。主要研究方向?yàn)閯?chuàng)傷彈道學(xué)。E-mail：liususu1006@139.com