《Materials》:Dynamic Impact Characteristics of Airdrop Cushioning Materials and a C ? σm Curve-Based Cushioning Pad Design Method
为减轻空投物资着陆过程中的冲击损伤,本研究提出了一种基于C ? σm(缓冲系数–最大应力)曲线的缓冲垫设计方法,旨在在吸能效率与轻量化需求之间实现平衡。研究建立了一个中型空投冲击仿真模型,并通过跌落冲击试验进行了验证;随后对三种具有代表性的缓冲材料——蜂窝纸板、聚氨酯泡沫(PU)和铝泡沫——开展了系统的动态冲击分析。三者的缓冲特性曲线比较表明,这三种材料均呈现凹形C ? σm分布特征(先下降后上升),且分别对应不同的空投缓冲最优应力范围:铝泡沫适用于高应力(≥500 kPa),聚氨酯泡沫适用于中等应力(350–450 kPa),蜂窝纸板适用于低应力(≤200 kPa)。吸能潜力随最优应力阈值升高而降低,而缓冲层厚度与空投载荷适用范围呈正相关。在低应力阶段,最大应力与能量密度之间存在显著函数依赖关系,因此厚度对吸能的影响可忽略。在材料脆值约束条件下,基于C ? σm曲线的图解法能够准确确定缓冲垫的最优厚度与承载面积。在设计案例3中,通过优化承载面积,所需缓冲厚度由100.5 cm降至25.0 cm,显著减小了缓冲体积。上述结果为空投货物冲击防护设计提供了可靠的材料层面认识与理论支持,并为空投缓冲材料的选择提供了基于材料本征力学响应的明确指导。
该文发表于《Materials》,围绕中型空投货物着陆冲击防护中的材料选择与缓冲结构定量设计问题展开研究。空投在应急救援、物流运输和特种物资投送中具有重要作用,但货物着陆瞬间的冲击载荷往往直接决定任务成败及载荷完好性。现有重载空投系统多依赖安装于平台底部的气囊缓冲系统,虽然能够承受较高冲击能量,但存在成本高、易倾覆、对复杂气象和高原环境适应性不足等问题,且对轻量化或中小型货物的冲击加速度控制不够精确。对于质量相对较低的中小型空投对象,采用材料缓冲吸收着陆能量是一条可行路径,但不同缓冲材料在空投场景中的系统比较仍不充分,尤其缺乏一种能够依据材料本构响应、货物质量和脆值要求来选择材料并同步确定缓冲垫厚度与承载面积的系统化理论方法。因此,研究人员开展本研究,目的在于建立面向中型空投工况的冲击仿真与材料筛选框架,并提出一种基于C ? σ
m曲线的定量设计方法,以弥补传统经验设计和大量跌落试验选材的不足。
研究人员首先从空投包装构型与冲击能量转换机理入手,分析托盘式空投货物自高度H下落至着陆时的重力势能、缓冲垫单位体积弹性变形能、最大应力σ
m与缓冲系数C之间的关系,建立了由动态冲击响应映射到工程设计参数的理论基础。文中指出,C ? σ
m曲线是将材料动态本构行为转化为工程缓冲设计曲线的核心工具,可用于评价单位体积吸能能力与最大应力之间的平衡关系。基于此,研究将空投缓冲问题从传统的经验试验主导模式,推进为可计算、可比较、可图解设计的定量化模式。论文的核心意义在于,不仅比较了三类典型轻质缓冲材料——蜂窝纸板、聚氨酯泡沫、铝泡沫——在相同空投相关边界条件下的动态冲击特征,还提出了一种可根据目标脆值和载荷水平反推缓冲垫参数的材料无关型设计原则,为中型空投防护提供了新的理论与方法学支撑。
研究采用的主要技术方法包括:建立中型空投托盘—货箱—缓冲垫—地面的三维数值模型,并在ANSYS/LS-DYNA R15.0中进行跌落冲击仿真;基于蜂窝纸板动态压缩试验结果及文献中的聚氨酯泡沫、铝泡沫应力–应变数据构建材料模型;通过实际跌落试验对蜂窝纸板空投组件仿真模型进行验证;在统一接触面积、厚度和跌落速度条件下,系统提取峰值加速度G
m、静应力σ
s、最大应力σ
m、能量密度u与缓冲系数C,并绘制σ
s ? G
m曲线、u ? σ
m曲线及C ? σ
m曲线,最终结合案例完成图解式缓冲设计。样本来源方面,模型验证所用实验对象为固定于托盘上的蜂窝纸板缓冲空投箱组件。
在“空投缓冲包装”部分,论文首先说明了传统中型空投系统的结构配置,包括空投平台、货物箱体、约束带、导向滑轮和系留环等。现实作业中,投送物通常直接捆扎在平台表面,而平台本身缺乏一体化缓冲设计,因此落地后冲击载荷几乎无衰减地传递至货物。研究据此明确指出,在平台与货物之间设置合适缓冲材料,是提高空投安全性的必要措施。随后,研究人员从重力势能向缓冲垫应变能转化的角度建立了分析框架,推导出缓冲系数C及其与最大加速度、最大应力之间的关系,并进一步说明利用C ? σ
m曲线可根据货物质量、冲击加速度和跌落高度求取所需缓冲厚度。这一部分为全文后续材料分析和设计方法建立了统一的力学基础。
在“空投仿真模型”部分,研究建立了中型空投缓冲包装的几何模型,并对货箱、托盘、绑带、地面与缓冲材料分别赋予相应力学参数。蜂窝纸板采用A型蜂窝芯,结合6 m/s动态压缩试验得到其应力–应变响应,并使用分段线性塑性模型描述其应变率相关行为。结果显示,蜂窝纸板压缩过程可分为弹性、屈服、平台和密实化四个阶段,仿真曲线与试验曲线吻合较好,说明模型能够反映其主要吸能机制。对于聚氨酯泡沫和铝泡沫,研究采用LS-DYNA中的可压溃泡沫模型(MAT_63),并依据已有文献数据输入应力–应变曲线。论文特别指出,在本研究涉及的低速冲击范围内,聚氨酯泡沫的高应变率敏感性较弱,铝泡沫的应变率硬化效应也可忽略,因此采用相应静态或文献压缩数据能够满足仿真精度要求。随后,研究通过不同网格密度的网格无关性分析证明参考网格的峰值加速度相对误差仅1.5%,并设置了接触算法、Rayleigh阻尼、时间步长缩放和hourglass控制等求解参数,以保证数值稳定性。最后,研究以300 kg空投组件进行跌落验证,仿真得到的峰值冲击过载为54.6 g,试验值为52.8 g,相对误差为3.4%;不同质量条件下,多数组别的相对误差控制在10%以内,证明该模型在冲击时程响应和蜂窝纸板压缩形态上均具有较好可信度。
在“4.1 蜂窝纸板缓冲特性及曲线”中,研究人员对不同质量空投箱在蜂窝纸板缓冲下的落地过程进行了仿真,并根据峰值加速度计算σ
m与C。结果表明,蜂窝纸板的过载G
m随货物质量增加呈“先下降后上升”趋势。轻载时,材料仅经历短暂弹性压缩,压缩行程极小,导致峰值力显著高于平台应力区间,从而出现较大过载;中等质量时,蜂窝结构稳定压溃并延长减速过程,G
m降至最优;重载继续增加后,蜂窝单元逐步密实化,吸能效率下降,G
m再次升高。通过结合相邻载荷段过载增长率的突变和最大应力陡增判据,研究确定150 mm蜂窝纸板的吸能饱和阈值为580 kg,对应临界最大应力155.65 kPa、临界过载33 g。由此说明蜂窝纸板在一定低应力范围内具有较好的缓冲能力,但超过饱和阈值后其塑性吸能难以维持,缓冲性能迅速劣化。对应的C ? σ
m曲线呈典型凹形,最低点位于540 kg工况,表明存在明确的最优缓冲应力区间。
在“4.2 不同厚度缓冲特性与应力–能量法”中,研究比较了50 mm、100 mm和150 mm蜂窝纸板在相同跌落速度下的σ
s ? G
m曲线和u ? σ
m曲线。结果显示,在相同静应力下,厚度越大,峰值加速度越低,说明更大的缓冲厚度提供了更长的压缩行程和更充分的能量耗散。三种厚度的σ
s ? G
m曲线均为U形,并各自对应一组最优静应力和最小G
m。另一方面,在低应力阶段,三种厚度的u ? σ
m曲线基本重合并近似线性增加,说明此阶段单位体积吸能能力主要由应力水平决定,厚度影响可忽略;超过临界拐点后,曲线分化明显,薄板在高应力阶段表现出更高单位体积能量密度。研究对低应力区数据进行了二次多项式拟合,拟合优度R
2为0.962,证明蜂窝纸板在σ
m ≤ 150 kPa范围内可由统一本构关系较准确预测。这一结果意味着,蜂窝纸板适合低静应力、轻载空投场景,但其平台应力和承载能力有限,不适于重载设备防护。
在“4.3 不同材料缓冲特性分析”中,研究将铝泡沫和聚氨酯泡沫引入相同边界条件下进行比较。两种材料在不同厚度下的σ
s ? G
m曲线与蜂窝纸板相同,均呈现先降后升的凹形特征,说明三类材料都存在最优静应力工作点。对比可见,随着厚度增加,铝泡沫和聚氨酯泡沫的峰值加速度均明显降低,最优静应力持续升高,且这种承载能力的提升幅度显著大于蜂窝纸板,表明二者更适于通过增加厚度来适配重载空投。以150 mm为例,铝泡沫最小G
m为26.95 g,聚氨酯泡沫为30.82 g,蜂窝纸板为28.75 g,说明三者在各自最优条件下都可提供良好防护,但适用负载范围不同。进一步比较指出,在低静应力轻载区,铝泡沫和聚氨酯泡沫因初始刚度更高而表现出更高G
m,因此不如蜂窝纸板适用;而在高静应力重载区,二者因平台应力更高、单位体积吸能更强,显著优于蜂窝纸板。基于150 mm数据,研究归纳出材料匹配范围:铝泡沫适用于20–25 kPa高静应力重载场景,聚氨酯泡沫适用于10–15 kPa中等静应力中载场景,蜂窝纸板适用于2–5 kPa低静应力轻载场景。与此同时,聚氨酯泡沫和铝泡沫的u ? σ
m曲线也显示出与蜂窝纸板一致的阶段性:低应力阶段近似统一本构,厚度效应不显著;超过屈服/密实化转折点后,塑性塌陷主导吸能。二者在低应力段的二次拟合同样具有很高拟合度,从而为工程设计提供了预测依据。
在“4.4 缓冲系数特性分析”中,论文从C ? σ
m曲线角度对三种材料进行综合比较。结果表明,三者的曲线均为典型凹形,反映了材料从压缩不足、吸能效率低,到进入最佳吸能区,再到过度密实化、性能退化的全过程。铝泡沫曲线最低点为σ
m = 611.20 kPa、C = 2.21,在最优点处具有三者中最高的最大应力和最低的缓冲系数,且最小点之后曲线斜率平缓,说明其在高应力范围内仍能维持稳定吸能效率。蜂窝纸板最低点位于σ
m = 125.65 kPa、C = 2.34,对应最低的最优工作应力,但最小点后曲线迅速上升,显示其在中高应力区能量吸收效率迅速恶化。聚氨酯泡沫最低点为σ
m = 399.40 kPa、C = 2.51,整体性能介于前两者之间。综合σ
s ? G
m、u ? σ
m和C ? σ
m曲线,研究确认三种材料可构成具有互补性的载荷覆盖体系:铝泡沫对应高载荷,聚氨酯泡沫对应中载荷,蜂窝纸板对应轻载荷。
在“4.5 空投应用案例分析”中,研究提出了基于C ? σ
m曲线的缓冲设计流程。该流程先收集货物质量m、允许最大加速度G、冲击速度v
0和缓冲面积A等输入参数,再计算静应力σ
s和最大冲击应力σ
m,随后将σ
m投影至三种材料的C ? σ
m曲线以确定适宜材料及对应C值,最后利用厚度公式求得缓冲层厚度h。案例1中,500 kg机械零件、脆值120 g,在全表面缓冲条件下计算得到σ
m为486.45 kPa,适合选用铝泡沫或聚氨酯泡沫,所需厚度分别为4.0 cm和4.5 cm;当落地速度由6 m/s升至8 m/s时,σ
s和σ
m不变,但厚度需相应增至7.0 cm和8.0 cm,说明速度影响厚度需求而不改变应力匹配。案例2中,300 kg仪器仪表、脆值50 g,计算得到σ
m接近蜂窝纸板最优区域,因此选择蜂窝纸板,厚度约9.7 cm。案例3中,300 kg精密部件、脆值20 g,若采用全表面缓冲,则σ
m仅48.6 kPa,远离最优工作区,导致C高达10.92、所需厚度达100.5 cm;当缓冲面积缩减至0.5 m
2后,σ
m提升至117.7 kPa,进入蜂窝纸板高效缓冲区,C降至2.72,厚度降至25.0 cm。该案例清楚表明,通过调节承载面积可提升材料工作应力,使其进入最佳吸能区,从而在满足脆值约束的同时显著降低缓冲体积。进一步仿真验证显示,三个案例的设计结果与模拟峰值加速度相对误差均在10%以内,且模拟值均低于设计脆值,证实该设计方法具有良好工程可行性。
讨论部分的核心在于,研究并未停留于“不同材料刚度不同”这一经验性认识,而是通过统一仿真条件下的系统比较和蜂窝纸板试验验证,建立了材料工作应力区间、吸能效率与几何设计变量之间的定量联系。研究表明,中型空投缓冲设计不宜仅依据单一材料性能指标或传统经验进行,而应同时考虑材料在目标应力区的C值、低应力阶段的统一本构关系以及厚度与承载面积对工作应力的调控作用。尤其是对高脆值精密件而言,单纯增加厚度并非最佳策略,合理缩小承载面积、提高材料工作应力进入最佳C ? σ
m区间,能够更有效地降低总体缓冲体积。论文同时指出,蜂窝纸板模型经过物理试验验证,而聚氨酯泡沫和铝泡沫部分目前主要基于文献材料数据和数值预测,因此其结果体现为在本研究框架下的可靠模拟结论。总体上,本文在中型空投缓冲领域提供了材料筛选、参数计算和结构协同优化的一体化思路。
研究结论部分可译为:本研究利用ANSYS/LS-DYNA R15.0对蜂窝纸板、聚氨酯泡沫和铝泡沫的缓冲性能进行了仿真,并通过跌落冲击试验对数值模型进行了验证,获得了可用于中型空投缓冲设计的可靠C ? σ
m曲线。结果表明:(1)三种材料均呈现具有明显最优效率点的凹形C ? σ
m曲线。铝泡沫适用于高应力范围,具有稳定效率和高承载能力;聚氨酯泡沫适用于中等载荷,性能较为均衡;蜂窝纸板则在低应力应用中具有轻质、低成本优势。三者共同覆盖了中型空投的大部分载荷范围,并形成系统化材料选择框架。(2)最优静应力随缓冲垫厚度增加而升高,厚度减小则使阈值降低。铝泡沫最优应力随厚度增加提升133%,表现出最大的适应性增益。最小峰值加速度随厚度增加而降低。在饱和前的低应力阶段,最大应力与能量密度遵循统一本构关系(R
2 > 0.96),可据此准确预测该范围内任意跌落高度–缓冲厚度比条件下的缓冲性能。(3)研究建立了基于C ? σ
m曲线的定量设计方法。通过调节承载面积,可将最大应力调控至最优缓冲效率区间,从而在满足脆弱货物防护要求的同时显著减小缓冲体积。该原理不依赖于特定材料,可扩展应用于其他缓冲材料与载荷条件。整体而言,本研究贡献了一个面向中型空投的C ? σ
m曲线缓冲设计框架、一套在一致空投工况下对三类代表性材料的系统比较表征方法,以及一种用于协同优化缓冲厚度与承载面积的材料无关型应力调控原则。