• ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

载人空投着陆冲击下不同乘员姿态的损伤影响

戴俊超 周云波 张进成 张明 王显会 孙晓旺

肖大武, 邱志聪, 巫祥超, 何立峰. 金属铍的压缩变形行为[J]. 爆炸与冲击, 2016, 36(2): 285-288. doi: 10.11883/1001-1455(2016)02-0285-04
引用本文: 戴俊超, 周云波, 张进成, 张明, 王显会, 孙晓旺. 载人空投着陆冲击下不同乘员姿态的损伤影响[J]. 爆炸与冲击, 2021, 41(1): 015901. doi: 10.11883/bzycj-2020-0073
Xiao Dawu, Qiu Zhicong, Wu Xiangchao, He Lifeng. Compressive deformation behaviors of beryllium[J]. Explosion And Shock Waves, 2016, 36(2): 285-288. doi: 10.11883/1001-1455(2016)02-0285-04
Citation: DAI Junchao, ZHOU Yunbo, ZHANG Jincheng, ZHANG Ming, WANG Xianhui, SUN Xiaowang. Effects of different postures on crew damage under the impact of manned airdrop landing[J]. Explosion And Shock Waves, 2021, 41(1): 015901. doi: 10.11883/bzycj-2020-0073

载人空投着陆冲击下不同乘员姿态的损伤影响

doi: 10.11883/bzycj-2020-0073
基金项目: 国家自然科学基金(11802140)
详细信息
    作者简介:

    戴俊超(1994- ),男,硕士研究生,19910625452@163.com

    通讯作者:

    周云波(1980- ),男,博士,副教授,yunbo31983@163.com

  • 中图分类号: O382; E923.3

Effects of different postures on crew damage under the impact of manned airdrop landing

  • 摘要: 针对某军用车辆在1 m高度进行无缓冲平台空投实验,并建立座椅与乘员的模拟模型。利用实验获取的座椅安装点冲击信号作为模拟模型的输入数据,并通过实验结果与模拟结果的对比验证了该模型的可靠性。借鉴航空工程相关研究,提出了一种将各关键损伤指标加以归一化的权重评价指标—加权损伤准则(weighted injury criteria,WIC)。研究了乘员仰卧角度和大小腿夹角两个姿态参数对乘员损伤的影响,并以WIC为优化目标,利用遗传算法完成参数优化工作。研究发现:对乘员小腿运动进行约束能降低乘员整体损伤响应,乘员对抗着陆冲击的最佳姿态为仰卧角47°~56°、大小腿夹角62°~68°。
  • 金属铍具有中子散射截面大、吸收截面小、硬度高、模量高、比强高、热学性能良好等特性,因此被广泛应用于航空航天、军事工业、医疗设备、焊接技术等多个技术领域,如中子反射层,反应堆第一壁材料、中子慢化剂,航空航天结构部件、精密仪表、光学器件及X射线管窗口等。国外已开展大量金属铍的变形行为研究,而国内开展的相关研究较少,且主要集中在常温静态拉伸性能方面,对其压缩力学行为尤其是动态压缩特性方面报道较少[1-9]。王零森等[1]研究了晶粒尺寸对铍静态拉伸力学性能的影响,发现随着晶粒度逐渐细化,铍材料的强度显著提高,而晶粒过粗或过细,延伸率均下降。许德美等[2-3]研究了组织缺陷对金属铍室温拉伸断裂行为的影响,其分析结果表明铍的“脆性”特征主要来源于杂质、片状晶体疏松和孔洞等初始缺陷,最关键因素是杂质的尺寸、间距和其在材料内部的分布形态。W.R.Blumenthal等[4-5]对不同制备工艺下的铍进行了较为系统的研究。实验结果表明,铍的压缩应力应变响应具有较强的应变率敏感性和一定的热软化效应,并指出孪生是高应变率下铍变形的主要机制。D.W.Brown等[6-8]系统开展了应变率对热压和轧制铍的力学性能和变形机理的研究工作,分析结果表明屈服强度对应变率不敏感,而加工硬化则受织构的影响具有较强的率相关性。T.Nicholas[9]和D.Breithaupt[10]研究了铍在常温102~103 s-1应变率下的动态压缩性能,结果表明铍具有良好的塑性,应变增大至0.25时样品才发生断裂。由此可见,国外开展的相关研究工作重点关注制备工艺、温度、应变率等条件对金属铍滑移及孪晶变形机制的影响研究,获得描述金属铍变形织构行为的本构模型参数。国内开展的研究则主要围绕金属铍静态拉伸应力状态下的“脆性”行为的微观变形机制,对其压缩行为研究工作较少,尤其是动态加载下温度、应变速率对其变形行为的影响未见相关研究报道。

    本文中利用材料实验机及Hopkinson杆装置系统开展了热等静压金属铍在不同温度、应变率下的压缩力学行为研究,获得金属铍压缩载荷下强度、塑性与实验温度、应变率之间的对应关系。并采用Johnson-Cook本构模型对获得的应力应变曲线进行拟合,模型计算结果与实验结果吻合较好。

    铍在机加后表面会有较大的残余应力,为了消除残余应力对测量结果的影响[11],室温力学实验前对样品进行了蚀刻处理,蚀刻剂配方为:H3PO4,750 mL;H2SO4,30 mL;Cr2O3,71 mg;H2O,200 mL。蚀刻方法为将铍试样放入酸洗液约50 s取出,用蒸馏水等清洗干净。

    静态力学实验在CMT5105型材料试验机及其配置的高温真空炉中进行,高温炉温度控制精度为±3 ℃,真空度优于1×10-2 Pa,试样在1 h内加热到规定温度,保温15 min后开始实验,应变率为1.0×10-3 s-1,测试温度范围为室温至800 ℃。动态压缩实验采用∅10 mm的Hopkinson杆装置。试样为∅5 mm×5 mm的圆柱体,应变率范围为0.5×103~2.5×103 s-1, 在常温下进行。

    图 1所示为铍在不同温度下的准静态压缩实验结果。由图 1应力应变曲线可以看出,金属铍在室温至800 ℃的温度范围内压缩变形具有良好的塑性。屈服强度和流动应力随实验温度升高而降低,加工硬化行为也随之降低。图 2所示为不同固定应变下的流动应力随实验温度的变化。由图中可以看到,在室温至200 ℃时,不同固定应变下流动应力均下降较快,高于200 ℃时流动应力下降趋势变缓,呈线性下降特征。当实验温度高于400 ℃时,不同应变下的流动应力值基本一致,这表明此时材料的塑性变形行为趋于理性塑性流动。

    图  1  金属铍在准静态条件下应力应变关系
    Figure  1.  Relation between stress and strain under quasi-static condition
    图  2  金属铍在准静态条件下流动应力随温度变化曲线
    Figure  2.  Relation between flow stress and temprature under quasi-static condition

    图 3所示为铍的动态压缩实验结果。可以看出,铍的屈服强度和加工硬化行为随应变率增大而显著增大,在初始变形阶段,加工硬化行为呈现非线性特征,随变形量增大,转变为线性硬化。由文献[4]可知,准静态和动态加载下,金属铍的塑性变形控制机制有显著区别。与大多数对称性低、滑移系统少的密排六方晶系金属一样,由于晶体的取向不利于发生滑移,孪生成为铍塑性变形的重要方式。在初始变形阶段,变形机制由位错滑移控制,随着变形增大,位错滑移困难,通过孪生协调变形,尤其在动态加载过程中,晶粒内部将产生大量的孪晶,由于滑移与孪生机制的竞争导致了不同应变率、不同应变下金属铍屈服强度和加工硬化行为的显著区别。

    图  3  金属铍的动态压缩力学行为
    Figure  3.  The dynamic compressive behavior of beryllium

    Johnson-Cook模型是目前应用最广泛的本构模型之一,模型中将流动应力表述为应变硬化效应、应变率效应和温度软化效应的乘积,方程的基本形式如下:

    σ=(A+Bεnp)(1+Cln˙ε˙ε0)(1Tm)
    (1)

    式中:σ为Von-Mises流动应力,εp为等效塑性应变,A为准静态下的屈服应力,B为应变硬化系数,n为应变硬化指数,˙ε0为参考应变率(可取准静态应变率),C为应变率敏感系数,T*为温度相关项,具体表达式为(T-Tr)/(Tm-Tr), TrTm分别为参考温度和熔化温度,一般取Tr为300 K,m为热软化系数。由式(1)可见,Johnson-Cook本构模型忽略了材料变形历史的影响,即如果材料服从Johnson-Cook本构模型,则不同应变率下的应力应变曲线是相似的。

    而由图 1~3中的应力应变曲线可以看到,不同温度或应变率下铍的应力应变曲线呈发散趋势,传统的Johnson-Cook本构模型已不适用。因此,本文中采用一个修正的Johnson-Cook本构模型对实验数据进行拟合,在应变硬化项中增加屈服强度温度相关线性函数,同时参考Zerrilli-Armstrong本构模型中描述hcp晶体结构材料变形硬化的函数关系式, 在幂指数应变硬化项中添加应变率指数硬化项和温度指数软化项,分别描述温度、变形历史对材料屈服强度和流动应力的影响,以及流动应力随应变率明显的增加趋势, 其表达式为:

    σ=[A(1A1T)+Bεnpea˙ε(B1+B2eβT)](1+Cln˙ε˙ε0)
    (2)

    和传统Johnson-Cook模型相比,修正模型中增加了4个参数。取准静态应变率10-3 s-1为参考应变率,Tm=1 557 K。本构拟合参数为:A=424 MPa, B=1 010 MPa, A1=1.487, B1=0.107 3, B2=0.885 4, n=0.485, α=0.000 39, β=-13.83, C=0.015。

    采用修正模型计算结果与实验结果对比如图 4所示,实线为采用修正Johnson-Cook本构模型的计算结果。可以看到,模型的计算结果与实验结果符合较好,修正后的Johnson-Cook本构模型能够较好地描述金属铍在不同温度、应变和应变率下的压缩变形行为。

    图  4  修正Johnson-Cook模型计算结果与实验结果对比
    Figure  4.  Comparison of experimental results with calculated results by modified Johnson-Cook model

    本文中研究了较宽温度范围和应变率下热等静压金属铍的压缩力学行为。结果表明铍的屈服强度和加工硬化行为随应变率的提高而显著增大,随温度的升高而降低。常温下其加工硬化行为在初始变形阶段呈现非线性特征,随变形增大转变为线性硬化。温度高于400 ℃时,其变形行为趋于理性塑性流动。考虑温度、变形历史对材料屈服强度和加工硬化的影响,对Johnson-Cook模型进行了修正,修正后的本构模型预测结果和实验结果吻合较好。

  • 图  1  实验装置整体布置

    Figure  1.  Overall arrangement of the experimental setup

    图  2  车内实验假人的姿势及加速度传感器的布置

    Figure  2.  The test dummy posture in vehicle and the acceleration sensor arrangement on the joist

    图  3  座椅安装点Z向加速度时间历程

    Figure  3.  Change of acceleration along Z direction with time at the mounting point on the seat

    图  4  假人的运动状态

    Figure  4.  Motion statuses of the dummy

    图  5  乘员约束系统有限元模型

    Figure  5.  A finite element model of the crew restraint system

    图  6  坐垫靠垫和安全带材料应力应变曲线

    Figure  6.  Stress-strain curves of cushion and seat belt materials

    图  7  实验与模拟假人损伤对比

    Figure  7.  Comparison of experimental and simulated dummy damage

    图  8  乘员约束系统设计

    Figure  8.  Passenger restraint system design

    图  9  小腿约束对乘员响应的影响

    Figure  9.  Effect of calf restraint on passenger response

    图  10  不同仰卧角度对成员身体几个关键部位动态响应的影响

    Figure  10.  Effects of different supine angles on dynamic responses of several key parts of a passenger’s body

    图  11  大小腿的不同夹角对成员身体几个关键部位动态响应的影响

    Figure  11.  Effects of different angles between the calf and the thigh on dynamic responses of several key parts of a passenger’s body

    图  12  第50代优化解集

    Figure  12.  The 50th generation optimization solution set

    图  13  加权损伤准则优化解集

    Figure  13.  Optimized solution set based on weighted injury criteria

    表  1  不同部位伤害标准和限值

    Table  1.   Injury standards and threshold values for different parts

    部位损伤标准注解损伤指标损伤指标限值
    头部GX/g头部X向加速度GX,lim/g26
    GZ/g头部Z向加速度GZ,lim/g20
    颈部Fn,Z/kN颈部压缩轴向力Fn,Z,lim/kN−1.3
    Mn,Y/(N·m)颈部力矩Mn,Y,liml/(N·m)135
    Nij轴向力和弯曲力矩线性合成准则1
    腰椎Flum,Z/kN腰椎轴向力Flum,Z,lim/kN8
    Mlum,Y/(N·m)腰椎Y向弯矩Mlum,Y,lim/(N·m)400
    盆骨Gp/g盆骨Z向加速度
    λdri动态响应指数λdri,lim13.4
    下载: 导出CSV

    表  2  主要结构材料参数

    Table  2.   Main structural material parameters

    材料密度/(g·cm−3)弹性模量/GPa泊松比屈服极限/MPa拉伸失效应变
    Q2357.8210.00.31235.00.24
    泡沫0.228.30.4019.70.50
    A3003H14铝2.770.00.30110.00.10
    下载: 导出CSV

    表  3  实验与模拟数据对比

    Table  3.   Comparison of experimental and simulated data

    方法GZ/gFn,Z/kNMn,Y/(N·m)Flum,Z/kNMlum,Y/(N·m)Gp/g
    实验16.70.7956.74.65126.717.9
    模拟18.80.7663.94.52151.215.6
    相对误差/%12.63.512.72.619.312.8
    下载: 导出CSV

    表  4  设计变量选取

    Table  4.   Selection of design variables

    设计参数采样点
    α/(°)010203040506070
    β/(°)6090120150
    下载: 导出CSV

    表  5  小腿有、无约束模拟数据对比

    Table  5.   Comparison of simulation data between the calf with and without restraint

    模拟条件GX/gGZ/gFn,Z/kNMn,Y/(N·m)Flum,Z/kNMlum,Y/(N·m)IdIw
    小腿无约束10.218.800.76363.94.521573.760.536 0
    小腿有约束8.715.380.71559.33.921223.660.411 6
    相对差值/%14.718.25.77.213.421.82.723.2
    下载: 导出CSV

    表  6  不同仰卧角度的仿真数据对比

    Table  6.   Comparison of simulation data of different supine angles

    α/(°)GX/gGZ/gFn,Z/kNMn,Y/(N·m)Flum,Z/kNMlum,Y/(N·m)IdIw
    08.715.380.5930.3183.92112.83.460.411 6
    109.3114.050.5580.2914.45149.03.180.425 3
    2010.8813.960.5750.2604.28176.23.10.432 6
    3013.2212.330.5650.2083.46223.91.840.410 5
    4015.9511.200.5870.1282.85225.01.820.385 6
    5013.0511.700.5450.1042.27202.01.850.343 6
    6013.8011.890.5640.0801.79239.01.470.346 6
    7013.9010.300.52490.0901.43279.01.490.349 3
    下载: 导出CSV

    表  7  不同大小腿夹角的空投冲击工况模拟数据对比

    Table  7.   Comparison of simulation data among airdrop impact conditions with different angles between the calf and the thigh

    β/(°)GX/gGZ/gFn,Z/kNMn,Y/(N·m)Flum,Z/kNMlum,Y/(N·m)IdIw
    609.811.660.4910.2473.87104.33.550.365 2
    908.715.380.5930.3183.92122.83.660.411 6
    12011.813.20.5660.3164.74127.93.820.435 5
    1508.3913.470.5280.2704.53105.54.420.405 6
    下载: 导出CSV

    表  8  优化结果和仿真结果的对比

    Table  8.   Comparison of optimization results and simulation results

    项目GX/gGZ/gFn,Z/kNMn,Y/(N·m)Flum,Z/kNMlum,Y/(N·m)IdIw
    优化结果0.49870.103 11.702011.350.326
    优化后模拟结果13.310.40.51580.106 01.722291.040.328
    误差/%32.81.213.922.90.6
    下载: 导出CSV

    表  9  优化后方案与原始方案的对比

    Table  9.   Comparison between the optimized solution and the original solution

    项目GX/gGZ/gFn,Z/kNMn,Y/(N·m)Flum,Z/kNMlum,Y/(N·m)IdIw
    初始值10.218.80.720 20.5704.521573.760.536
    优化后模拟结果13.310.40.515 80.1061.722291.040.328
    相对差值/%−3044.720.481.461.9−457238.8
    下载: 导出CSV
  • [1] 刘鑫. 基于气囊缓冲的载人空降乘员防护理论与方法[D]. 长沙: 湖南大学, 2009. DOI: 10.7666/d.d146725.
    [2] 刘鑫, 张志勇. 基于气囊缓冲的载人空降乘员防护装置优化设计 [J]. 机械工程学报, 2010, 48(21): 168–174. DOI: 10.3901/JME.2012.21.168.

    LIU X, ZHANG Z Y. Optimal design of passenger’s protection devices in manned airdrop based on airbag cushion [J]. Journal of Mechanical Engineering, 2010, 48(21): 168–174. DOI: 10.3901/JME.2012.21.168.
    [3] STAPP J P, TAYLOR E R. Space cabin landing impact vector effects on human physiology [J]. Aerospace Medicine, 1964, 35: 1117–1133.
    [4] 成自龙, 韩延方, 曾文艺, 等. 人体坐姿着陆冲击(+G z)耐限区间的研究 [J]. 航天医学与医学工程, 1997, 10(5): 340–343. DOI: 10.16289/j.cnki.1002-0837.1997.05.008.

    CHENG Z L, HAN Y F, ZENG W Y, et al. Human tolerance to landing impact (+G z) in sitting position [J]. Space Medicine and Medical Engineering, 1997, 10(5): 340–343. DOI: 10.16289/j.cnki.1002-0837.1997.05.008.
    [5] 刘炳坤, 王宪民, 王玉兰, 等. 不同体位着陆冲击时人体的动态响应 [J]. 航天医学与医学工程, 2001, 14(2): 121–123. DOI: 10.3969/j.issn.1002-0837.2001.02.010.

    LIU B K, WANG X M, WANG Y L, et al. Human dynamic response to landing impact in selected body orientations [J]. Space Medicine and Medical Engineering, 2001, 14(2): 121–123. DOI: 10.3969/j.issn.1002-0837.2001.02.010.
    [6] 冯宇, 于广春, 吴增斌, 等. 某轮式空降战车模拟着地冲击试验方法研究 [J]. 车辆与动力技术, 2019(2): 62–64. DOI: 10.16599/j.cnki.1009-4687.2019.02.012.

    FENG Y, YU G C, WU Z B, et al. Research on the simulated landing impact test method of a wheeled airborne fighting vehicle [J]. Vehicle and Power Technology, 2019(2): 62–64. DOI: 10.16599/j.cnki.1009-4687.2019.02.012.
    [7] 俞彤, 王显会, 周云波, 等. 爆炸环境下车辆地板加速度对乘员小腿损伤的影响 [J]. 科学技术与工程, 2018, 18(21): 339–344. DOI: 10.3969/j.issn.1671-1815.2018.21.053.

    YU T, WANG X H, ZHOU Y B, et al. The effects of vehicle floor acceleration on the damage of occupant's tibia in explosive environment [J]. Science Technology and Engineering, 2018, 18(21): 339–344. DOI: 10.3969/j.issn.1671-1815.2018.21.053.
    [8] 梁锐, 黄世霖, 杜汇良, 等. 高G值复合着陆冲击下头、颈部损伤成因 [J]. 中国临床康复, 2005, 9(25): 28–30. DOI: 10.3321/j.issn:1673-8225.2005.25.014.

    LIANG R, HUANG S L, DU H L, et al. Causes of head and neck injury under high G landing impact [J]. Chinese Journal of Clinical Rehabilitation, 2005, 9(25): 28–30. DOI: 10.3321/j.issn:1673-8225.2005.25.014.
    [9] North Atlantic Treaty Organization. Procedures for evaluating the protection level of armored vehicles-mine threat: AEP-55: Volume 2[S]. Brussels: Allied Engineering Publication, 2011.
    [10] 通过载人航天器发射、中止、着陆引起的伤害来确定NASA风险准则[J]. 载人航天信息, 2015(2): 1−13.

    Determine NASA risk criteria through injuries caused by manned spacecraft launch, suspension, and landing[J]. Manned Spaceflight Information, 2015(2): 1−13.
    [11] ECE UN. Uniform provision concerning the approval of vehicles with regard to the emission of pollutants according to engine fuel requirements [J]. UN ECE Regulation, 2005(83).
    [12] VIANO D C, AREPALLY S. Assessing the safety performance of occupant restraint systems [M]. DOI: 10.4271/902328.
    [13] BROWN W K, ROTHSTEIN J D, FOSTER P. Human response to predicted Apollo landing impacts in selected body orientations [J]. Aerospace Medicine, 1966, 37(4): 394–398.
    [14] TABIEI A, LAWRENCE C, FASANELLA E L. Validation of finite element crash test dummy models for predicting orion crew member injuries during a simulated vehicle landing: NASA TM-2009-215476 [R]. 2009.
    [15] 王心怡, 刘晓颖, 陈聪, 等. 跌落工况下的人体腰椎损伤风险的几种影响因素 [J]. 汽车安全与节能学报, 2018, 9(2): 178–185. DOI: 10.3969/j.issn.1674-8484.2018.02.008.

    WANG X Y, LIU X Y, CHEN C, et al. Factors for affecting risk of human lumbar spine injuries under dropping conditions [J]. Journal of Automotive Safety and Energy, 2018, 9(2): 178–185. DOI: 10.3969/j.issn.1674-8484.2018.02.008.
    [16] SIMPSON T M, MAUERY T M, KORTE J J, et al. Kriging models for global approximation in simulation-based multidisciplinary design optimization [J]. AIAA Journal, 2001, 39(12): 2233–2241. DOI: 10.2514/2.1234.
  • 期刊类型引用(2)

    1. 宋殿义,蒙朝美,谭清华. 基于数值模拟的蜂窝钢管约束混凝土靶抗侵彻机理. 弹道学报. 2024(03): 82-90 . 百度学术
    2. 李萌,武海军,董恒,任光,张鹏,黄风雷. 基于机器学习的混凝土侵彻深度预测模型. 兵工学报. 2023(12): 3771-3782 . 百度学术

    其他类型引用(3)

  • 加载中
图(13) / 表(9)
计量
  • 文章访问数:  705
  • HTML全文浏览量:  300
  • PDF下载量:  78
  • 被引次数: 5
出版历程
  • 收稿日期:  2020-03-19
  • 修回日期:  2020-06-22
  • 刊出日期:  2021-01-05

目录

/

返回文章
返回