聚苯乙烯泡沫材料的动态压缩特性

雷忠琦; 姚小虎; 龙舒畅; 常建虎; 王海波

doi:10.11858/gywlxb.20180655

聚苯乙烯泡沫材料的动态压缩特性

doi: 10.11858/gywlxb.20180655

1.
华南理工大学土木与交通学院力学系，广东广州　510640
2.
合肥美的电冰箱有限公司，安徽合肥　230601

基金项目: 国家自然科学基金（11672110）

详细信息

作者简介:
雷忠琦（1993－），男，硕士，主要从事冲击动力学研究. E-mail: lzq000@foxmail.com

通讯作者:
姚小虎（1974－），男，博士，教授，博士生导师，主要从事冲击动力学及纳米力学研究. E-mail: yaoxh@scut.edu.cn

中图分类号: O347.1
计量
- 文章访问数: 7857
- HTML全文浏览量: 5149
- PDF下载量: 122
出版历程
- 收稿日期: 2018-10-15
- 修回日期: 2018-11-12

Dynamic Compression Characteristics of Polystyrene Foam Materials

1.
Department of Engineering Mechanics, School of Civil Engineering and Transportation, South China University of Technology, Guangzhou 510640, China
2.
Hefei Midea Refrigerator Division, Hefei 230601, China

摘要

摘要: 通过万能材料试验机和落锤式冲击试验装置，对发泡聚苯乙烯泡沫材料进行准静态和动态压缩试验，探讨密度和加载速率对材料动态压缩特性的影响。基于落锤试验数据，考虑密度相关性，通过修正得到恒定应变率下材料动态本构关系的经验公式。基于LS-DYNA中的MAT57和MAT163以及ABAQUS中的Low Density Foam和Crushable Foam 4种材料模型，建立了适用于有限元仿真的发泡聚苯乙烯泡沫本构模型。通过模拟落锤冲击过程，对比试验结果发现：MAT163和Crushable Foam模型能较好地预测材料动态响应和能量吸收性能，验证了动态本构模型的可靠性，并且这两种特定的材料模型在模拟发泡聚苯乙烯泡沫冲击碰撞时具有良好的适用性。
- 聚苯乙烯泡沫 /
- 应变率 /
- 密度 /
- 动态压缩特性 /
- 本构模型
Abstract: The quasi-static and dynamic compression tests of expanded polystyrene (EPS) foam material were carried out through a universal material tester and a drop-weight impact machine. The density and loading rate effects of the dynamic compression characteristics for the polystyrene foam material were discussed. By considering the density correlation, the empirical formula of dynamic constitutive relation under constant strain rate was modified based on the drop-weight test data. The constitutive model of EPS foam was established based on the material finite element models of MAT57, MAT163 in LS-DYNA, and Low Density Foam and Crushable Foam in ABAQUS. By simulating the impact process of drop-weight and comparing with the test results, it shows that MAT163 and Crushable Foam model can predict the dynamic response and energy absorption performance better. The results verify the reliability of the dynamic constitutive model. Meanwhile, it shows that these two specific material models have good applicability in simulating the impact problem of EPS foam.
- polystyrene foam /
- strain rate /
- density /
- dynamic compression characteristics /
- constitutive model

HTML全文

火炮发射时，后效期高温高压的火药气体从膛口瞬时流出，对火炮主体形成很强的后坐力，对炮架产生强冲击作用，影响火炮射击精度和机动性能，制约高性能火炮的发展。炮口制退器作为一种安装在炮口部位的排气装置，通过控制后效期火药气体的流量分配和气流速度，减小射击时火药气体作用于后坐部分的冲量，并为炮身提供制退力，减小火炮后坐动能和炮架的射击载荷，因而成为一种广泛应用的反后坐技术^[1]。然而，炮口制退器的使用具有一定的负面效应。一方面，炮口制退器造成炮口区域强激波，加剧后效期火药气体流场的复杂性及膛口焰现象，加大火炮发射初速扰动及膛口冲击波超压，对火炮身管周围设备及操作人员造成不良影响；另一方面，加装制退器后火炮身管本身质量及挠度的增大也对射击精度造成不利影响。炮口制退器自身的质量和结构决定了其制退性能及负面危害表现，通过选材及结构的优化设计实现高性能、低危害的综合性能成为当前炮口制退器研究与发展的方向^[2–6]。

按照结构形式的不同，传统的炮口制退器分为冲击式、反作用式及冲击反作用式。由于工艺、重量的限制，高效率与低负面效应往往难以兼顾。近年来，增材制造技术的发展为复杂异型结构的加工制造以及钛合金等难加工轻金属材料的成型制备提供了有力支撑^[7]，也为传统装置结构的优化与创新设计提供了更大的自由空间。随着计算机技术和计算流体力学的发展，数值模拟已成为低成本研究膛口流场及膛口装置的重要手段，例如：张焕好等^[8–9]、代淑兰等^[10]利用数值模拟方法，对不同类型炮口制退器的膛口流场波系、膛口激波及二次焰特征进行了详细研究，得到了制退器装置效率；Lei等^[11]、Chaturvedi等^[12]通过将流场仿真与流固耦合相结合，实现了炮口制退器结构应力与变形响应研究及性能评价。本研究利用增材制造技术优势，针对机载火炮小口径炮口制退器综合性能的设计需求，提出一种叠加冲击式与反作用式传统制退器结构特征与优点的钛合金新型炮口制退器结构方案，对膛口流场发展过程及特征进行数值模拟分析，并以传统的高效率冲击式炮口制退器为参照，对新型结构制退器的综合性能予以评价。

1. 制退器结构

炮口制退器的基本构型是后端与炮管封闭相接，前端收缩，中间拥有膨胀腔室、内部挡板和侧孔的圆筒形结构。冲击式炮口制退器有较大的腔室直径和侧孔面积，火药气体膨胀较为充分，大部分经侧向流出并冲击反射挡板，制退效率高，但炮口冲击波危害也较大。反作用式炮口制退器的直径一般较小，拥有多排小面积侧孔，没有或有很小的反射挡板，制退效率较低，但膛口火药气体流场相对简单，对弹丸初速的扰动小，不会影响大威力脱壳穿甲弹发射。基于上述两种制退器的结构特征和性能的经验规律，本研究以低冲击波危害和高制退效率为导向，在三维建模平台上设计一种面向机载航炮的30 mm小口径新型结构炮口制退器，如图1所示。

图 1 新型结构炮口制退器

Figure 1. New-structure muzzle brake

下载: 全尺寸图片幻灯片

新型制退器内层为冲击式内壁结构，见图1，其上设置编号为1～3的3组阵列式大面积后倾侧孔，旨在引导火药气体出膛后的侧向分流，以实现类似于冲击式炮口制退器的高制退效率。此外，3组内层侧孔的后倾角度及面积采用非均等设计，参照一般冲击式炮口制退器的结构参数，对三维设计模型直接设定或导出得到内层侧孔的非均等结构参数，见表1。其中：为了限制火药气体在临近膛口位置的向后偏转角度及炮口后方冲击波的强度，第1组侧孔的后倾角度及面积均小于其他两组；根据火药气体出膛以后沿轴向流量逐渐衰减的特点，第2组侧孔的面积设置为大于其他两组。新型结构制退器外层为胡椒瓶式套筒结构，其上设置成组分布的阵列式小面积圆孔，各组圆孔与内层侧孔周向重叠，旨在借鉴反作用式炮口制退器的特点，对火药气体进行分散和疏导，从而抑制侧向强激波的产生，实现较低的冲击波危害效益。新型制退器弹孔处为收缩式前端面，并设有分流孔，以在出口处进一步减小中央弹孔流量及火药气体在制退器出口部位的膨胀扰动。内层结构与外层结构通过布置在孔间的阵列式筋板相连，实现整体结构的一体化，便于增材制造成型。同时，结构中不存在过于复杂或与轴向夹角过大的成型曲面，满足增材制造高质量成型工艺要求。设计方案选材为3D打印钛合金粉末，材料参数^[13]见表2。为评价新型炮口制退器的综合性能，以图2所示的传统冲击式炮口制退器为对比分析对象。

表 1 新型炮口制退器结构参数

Table 1. Structural parameters of new-structure muzzle brake

Structure	Backward angle/(°)	Area/mm²
Inner side hole 1	5	392.96
Inner side hole 2	10	470.34
Inner side hole 3	10	392.96
Outer side hole	0	78.50
Front hole		146.54

下载: 导出CSV

| 显示表格

表 2 3D打印钛合金材料参数

Table 2. Material parameters of titanium alloy

Material	Density/(g·cm⁻³)	Young’s modulus/GPa	Yield strength/MPa	Tensile strength/MPa	Poisson’s ratio
TC4	4.43	118	944	1058	0.3

下载: 导出CSV

| 显示表格

图 2 冲击式炮口制退器

Figure 2. Impulsive type muzzle brake

下载: 全尺寸图片幻灯片

2. 数值方法与模型参数

采用无黏三维Euler方程描述带制退器的膛口流场气体流动。考虑到计算目的及效率，对炮口气流进行适当简化，忽略气固两相性、多组分、化学反应以及运动弹丸的影响，将火药气体与外界大气视为理想气体，且满足气体状态方程。有限元仿真的封闭方程组为

$\frac{\partial {U}}{\partial t}+\frac{\partial {F}}{\partial x}+\frac{\partial {G}}{\partial y}+\frac{\partial {H}}{\partial z}=0$

(1)

$e=\frac{p}{\gamma -1}+\frac{1}{2}\rho \left({u}^{2}+{v}^{2}+{w}^{2}\right)$

(2)

$pV=nRT$

(3)

式中：e为单位质量气体的总能量， $p$ 为火药气体压强，V为气体体积，T为体系温度， ${U}=[\rho,\rho u,\rho v,\rho w,e{]}^{\rm{T}}$ ， ${F}=[\rho u,\rho {u}^{2}+p,\rho uv,\rho uw,(e+p)u{]}^{\rm{T}}$ ， ${G}=[\rho v,\rho uv,\rho {v}^{2}+p,\rho vw,(e+p)v{]}^{\rm{T}}$ ， ${H}= [\rho w,\rho uw,\rho vw,\rho {w}^{2}+p,(e+p)w{]}^{\rm{T}}$ ， $\,\rho$ 为火药气体密度，u、v、w分别为x、y、z方向的速度分量， $\gamma$ 为理想气体绝热指数，R为通用气体常数，n为气体的物质的量。方程的离散采用有限体积法，时间推进采用二阶精度Runge-Kutta法，对流项采用二阶精度Roe格式，湍流模型采用单方程Spalart-Allmaras模型，数值计算在ANSYS FLUENT 19.0中进行。

后效期火药气体排空过程数值模拟的计算域设置为以炮口端面中心为原点、直径为3 m的球形区域，以使火药气体充分发展。计算初始时刻为弹丸飞离炮口时刻，以膛底位置为x=0，此时膛内气体沿炮轴方向任意一点的压力p_x、速度v_x和温度T_x的初始分布状态表示为

${p}_{x}={p}_{\mathrm{d}}\left[1+\frac{{m}_{\mathrm{w}}}{2{\psi }_{1}{m}_{\mathrm{q}}}\left(1-\frac{{x}^{2}}{{L}^{2}}\right)\right]$

(4)

${v}_{x}=\frac{{v}_{\mathrm{d}}}{L}x$

(5)

${T}_{x}=\frac{{p}_{x}{M}_{\mathrm{g}}}{{\rho }_{\mathrm{g}}R}$

(6)

${\rho }_{\mathrm{g}}=\frac{{m}_{\mathrm{w}}}{{V}_{\mathrm{c}}+LS}$

(7)

式中： $\,\rho$ _g为膛内火药气体的平均密度，m_w、m_q、M_g分别为装药量、弹丸质量、火药气体摩尔质量，V_c、L、S分别为药室容积、身管行程长度、身管横截面积， ${\psi }_{1}$ 为内弹道次要功系数，v_d为弹底火药气体速度，p_d为膛口处弹底压力。本算例中，未加装炮口制退器时火炮身管质量M₀=7.09 kg，初始时刻弹底火药气体速度v_d=780 m/s，p_d=27.9 MPa，膛口外为常温常压条件。计算时间步长取0.01 ms，当膛内气压降低至约0.2 MPa时视为后效期结束，终止计算。

3. 结果分析

3.1 膛口冲击波

图3、图4所示为新型制退器轴向和径向截面处膛口冲击波在不同时刻的压力云图，为突出显示流场发展过程，图中隐去流场超高压及部分负压区域。图5为火药气体流场充分发展时轴向和径向截面处速度矢量图。可以看出，后效期内火药气体在流入和穿出新型制退器的过程中，与内外双层结构依次产生交互作用，形成独特的流场结构与激波形貌，并构成制退效益的基础。后效期膛口冲击波的形成过程可以分为3个阶段。（1）火药气体膨胀加速进入制退器腔室后，依次流经内层侧孔和外层小孔喷出，形成侧孔冲击波，见图3（a）和图4（a）。（2）随着火药气体的向前推进和持续补充，沿炮轴方向各外层圆孔依次形成各自的侧孔冲击波，并且伴有追赶、相交和叠加过程，最终在炮口一侧合并成为一个大的冲击波。由于膨胀速率不同，在此过程中出现了后排侧孔冲击波追赶并超过前排侧孔冲击波的现象，见图3（b）和图3（c）；与此同时，周向上不同角度的三维侧孔冲击波也逐渐相交合并，见图4（b）和图4（c）。（3）火药气体流出弹孔以后，在出口处形成迅速膨胀的弹孔冲击波，并与侧孔冲击波相交，形成具有各向叠加特征的典型膛口波系结构，见图3（d）、图4（d）和图5。

图 3 轴向截面膛口压力云图

Figure 3. Pressure nephogram at axial plane

下载: 全尺寸图片幻灯片

图 4 径向截面膛口压力云图

Figure 4. Pressure nephogramat at radial cross section

下载: 全尺寸图片幻灯片

图 5 膛口流场速度矢量图

Figure 5. Velocity vector diagram of muzzle flow field

下载: 全尺寸图片幻灯片

3.2 超压分布

图6和图7为炮轴中心线、侧孔位置径向轴线处不同时刻压力分布曲线。可以看出，火药气体到达弹孔之前，不断提升新型制退器腔室内的压力水平。t=0.65 ms时，火药气体从中央弹孔流出并迅速提升当地压力，直至t=1.00 ms时弹孔处压力峰值超过制退器腔室内压力最大值，见图6（a）和图6（b）。与此同时，径向轴线方向的压力分布呈现随时间剧烈波动的高瞬态特征，见图7（a）。t=1.00 ms之后，火药气体流场逐渐衰减，腔室内外压力水平均逐渐下降（见图6（c）和图7（b）），同时径向轴线方向压力水平也呈现出规则的逐级下降趋势。

图 6 不同时刻炮轴中心线压力分布曲线

Figure 6. Pressure distributions along central axis of gun at different moments

下载: 全尺寸图片幻灯片

图 7 不同时刻侧孔径向轴线压力分布曲线

Figure 7. Pressure distributions along radial axis of side holes at different moments

下载: 全尺寸图片幻灯片

根据机载小口径航炮的实际工况，以炮口端面为中心、炮膛轴线为正方向，在0.5和1.0 m 的周向距离上，炮口区域侧向0°～150°范围内设置若干膛口流场超压监测点，如图8所示，以精确分析加装制退器后炮口区域的超压特征，评估炮口冲击波对机载平台的危害效应。得到的后效期内各点超压峰值监测结果与差值情况如表3所示。其中，超压值代表压力监测值高于环境压力的部分，加装新型制退器与传统冲击式炮口制退器情况下的超压峰值与光膛口情况下的超压峰值差值分别记为Δ_new、Δ_impact，负值代表产生减小超压效果，正值代表超压增大。同时，根据超压峰值，绘制了不同情况下炮口超压分布，如图9所示。

表 3 不同制退器的超压峰值监测结果

Table 3. Overpressure peak values of different muzzle brakes

Monitor point	Overpressure peak/kPa			Δ_new/kPa	Δ_impact/kPa
Monitor point	Smooth muzzle	New-structure muzzle brake	Impact muzzle brake	Δ_new/kPa	Δ_impact/kPa
P₀₁	117.84	362.54	200.83	244.70	82.99
P₀₂	46.31	65.46	53.20	19.15	6.89
P₁₁	83.80	79.03	44.58	−4.76	−39.21
P₁₂	28.07	23.61	19.56	−4.46	−8.51
P₂₁	48.84	29.28	30.90	−19.56	−17.93
P₂₂	19.76	11.45	15.40	−8.31	−4.36
P₃₁	22.29	18.85	50.87	−3.45	28.57
P₃₂	10.54	6.99	15.30	−3.55	4.76
P₄₁	28.27	26.45	41.75	−1.82	13.48
P₄₂	11.55	9.54	9.73	−2.01	−1.82
P₅₁	34.75	44.68	40.09	9.93	5.34
P₅₂	15.71	14.79	13.58	−0.91	−2.13

下载: 导出CSV

| 显示表格

图 8 膛口流场超压监测点阵

Figure 8. Overpressure monitoring points of muzzle flow filed

下载: 全尺寸图片幻灯片

图 9 膛口超压峰值分布对比

Figure 9. Comparison of overpressure peak distribution

下载: 全尺寸图片幻灯片

由表3和图9可知，相比于光膛口情况，两种类型制退器的膛口冲击波超压特点可总结如下：（1）炮口正前0°方向，0.5和1.0 m两个距离上两种制退器均造成超压峰值的明显增大，并且新型制退器的增压幅度更大；（2）自炮口侧向30°方向向后，新型制退器膛口冲击波超压峰值呈U形分布，冲击型制退器呈W形分布；（3）30°～120°侧向范围内，新型制退器均产生减小超压效果，各点超压峰值降低幅度为1.82～19.56 kPa，正侧方90°处超压峰值最小，该方位0.5和1.0 m距离处分别为18.85和6.99 kPa；（4）30°～120°侧向范围内，冲击型制退器对应的超压峰值高低不等，在正侧方90°方向0.5和1.0 m距离上均形成增大超压效果，超压峰值分别为50.87和15.30 kPa；（5）150°方向0.5 m处，新型制退器对应的超压峰值略高于冲击型制退器及光膛口情况，1.0 m处同样略高于冲击型制退器但低于光膛口情况。

综上可知，新型制退器的膛口冲击波表现出高超压集中于炮口正前方、低超压集中于正侧方（炮口与机载平台直线距离最近位置处）的分布特征，同时其后侧方的整体超压水平也较好地控制在低于或近似光膛口情况，相比传统冲击型炮口制退器，对于缓解因制退器安装导致的机身及炮口装置的冲击侵害更具优势。同时应该注意到，为了追求高制退效率，新型制退器设置了较大的侧孔后倾角度及向后偏转面积，因此在炮口后侧150°方向0.5 m距离处超压出现一定的增大，说明新型制退器的结构参数仍存在改进空间，以实现高制退效率和更低冲击波危害的进一步平衡。

3.3 制退效率及结构强度

为分析后效期内新型制退器受载变化及制退力产生机制，沿炮膛轴线在新型炮口制退器腔室内层侧孔1～3位置及前端面侧孔处分别设置压力监测点，如图10所示。图11为后效期初期阶段各监测点压力随时间的变化曲线。图12和图13分别为后效期初期制退器受力、后效期全程火炮身管各部分受力随时间的变化曲线。

图 10 制退器腔内压力监测点

Figure 10. Monitoring points inside muzzle brake

下载: 全尺寸图片幻灯片

图 11 腔内各点压力随时间变化

Figure 11. Variation of overpressure inside muzzle brake with time

下载: 全尺寸图片幻灯片

图 12 后效期初期制退器受力曲线

Figure 12. Curve of muzzle brake force in the early gas ejection period

下载: 全尺寸图片幻灯片

图 13 膛底、制退器及总受力曲线

Figure 13. Curves of chamber bottom force, muzzle brake force and total thrust

下载: 全尺寸图片幻灯片

结合图11和图12可知，火药气体的流出和偏转作用不断产生对火炮身管的后坐力和向前的制退力，后效期初期制退器受力曲线呈现阶梯式上升趋势。图12中点Ⅰ～点Ⅳ处的抬升显示了当火药气体依次流经制退器内层侧孔时逐级产生的冲撞和制退效果。火药气体流出中央弹孔之前，对制退器前端面的冲击作用使制退器受力曲线继续抬升，直至达到峰值点Ⅴ，此后由于火药气体对制退器外壁面产生反向冲击作用，又使制退器的正向受力水平略降至峰值以下。

依据仿真监测结果，首先由总受力F随时间的变化曲线计算得到后效期火药气体对身管的总冲量I，再结合动量守恒定律及冲量定律，计算得到后效期起始时刻身管自由后坐速度v、后效期结束时刻身管最大后坐速度v_max以及火炮后坐动能E，进而依据定义得到制退器的效率 $\eta$ ，具体计算公式为

$v=\frac{{m}_{\mathrm{q}}+0.5{m}_{\mathrm{w}}}{M+{m}_{\mathrm{w}}+{m}_{\mathrm{q}}}{v}_{\mathrm{d}}$

(8)

$I={\int }_{0}^{\tau }F\mathrm{d}t=M({v}_{\mathrm{m}\mathrm{a}\mathrm{x}}-v)$

(9)

$E=\frac{M{v}_{\mathrm{m}\mathrm{a}\mathrm{x}}^{2}}{2}$

(10)

$\eta =\frac{{E}_{0}-E}{{E}_{0}}\times 100\text{%}$

(11)

式中： $\tau$ 为后效期持续时间，E₀为光膛口时火炮后坐动能，M为加装不同制退器后身管总质量。本算例中，设光膛口时火炮身管质量为7.09 kg，各种情况下的计算结果如表4所示，其中负号代表火炮发射相反方向。可见，新型制退器的制退效率达到45.62%，满足小口径炮口制退器的使用要求，并略超过以高效率为特点的传统冲击式炮口制退器性能水平。不过，新型炮口制退器情况下火药气体对膛口装置的反作用冲量相比传统冲击式炮口制退器的提升不明显，说明新型炮口制退器在内部结构参数上还有优化空间，可以进一步提升新结构形式的制退性能。

表 4 制退器效率对比

Table 4. Comparison of efficiency for different muzzle brakes

Structure	M/kg	I/(N·s)	v/(m·s⁻¹)	v_max/(m·s⁻¹)	E/J	$\eta$ /%
Smooth muzzle	7.09	−123.68	−47.00	−64.45	14719.0
New-structure muzzle brake	8.37	−29.40	−40.23	−43.74	8003.2	45.62
Impact muzzle brake	8.10	−26.15	−41.47	−44.70	8091.3	45.03

下载: 导出CSV

| 显示表格

制退器与后效期火药气体相互作用过程中，其本身所承受的动态载荷作用通过流固耦合技术由流体计算结果传递至制退器各壁面，并通过瞬态动力学有限元计算方法求得变形与等效应力随时间的变化曲线，如图14和图15所示。后效期初期，高温高压的火药燃气出膛后迅速膨胀，在剧烈变化的流体压力动态载荷作用下，制退器所承受的等效应力及其变形也在短时间内迅速抬升达到峰值；约2 ms以后，随着流场的充分发展并衰减，制退器结构的应力及变形转为逐渐下降趋势。

图 14 制退器结构的等效应力随时间变化曲线

Figure 14. Variation of equivalent stress with time of muzzle brake

下载: 全尺寸图片幻灯片

图 15 制退器结构的变形随时间变化曲线

Figure 15. Variation of deformation with time of muzzle brake

下载: 全尺寸图片幻灯片

图16和图17分别为新型制退器的等效应力和变形在其峰值时刻的分布云图。t=0.25 ms时，火药气体首先对炮口端面处的侧孔产生冲击，并在该时刻造成距离炮口端面第2排外层圆孔处出现全局最大范式等效应力，为224.37 MPa。t=1.10 ms时，火药气体已从中央弹孔流出，火药气体对制退器前端面的冲击作用使得最大变形出现在前端面与中央弹孔的结合处，约为14 μm，其他部位的变形呈现从炮口端面沿炮口正向递减趋势。由以上结果可知，后效期内制退器结构的最大应力小于钛合金材料许用应力，满足结构强度的使用要求。

图 16 0.25 ms时制退器结构的等效应力云图

Figure 16. Equivalent stress distribution of the muzzle brake at 0.25 ms

下载: 全尺寸图片幻灯片

图 17 1.10 ms时制退器结构的变形云图

Figure 17. Deformation distribution of the muzzle brake at 1.10 ms

下载: 全尺寸图片幻灯片

4. 结　论

提出了一种叠加冲击式内壁与反作用式外孔特征的新型小口径钛合金炮口制退器结构方案，基于膛口流场数值仿真、流固耦合技术等对其性能进行分析，并与传统冲击型炮口制退器进行对比。结果表明：新型结构方案满足小口径炮口制退器性能要求，制退效率相比传统冲击型制退器略有提高；在机载航炮应用背景下，其膛口冲击波超压分布特征及峰值水平相比传统冲击型制退器具有更好的机身与炮口装置防护优势和低冲击波危害效益，同时结构强度满足使用要求。新型结构方案对于新材料和新制造技术背景下炮口制退器的创新设计和发展具有参考意义。本研究所展示的新型制退器结构参数在高效率和更低负面危害的平衡上具有进一步改进与优化空间。研究过程中膛口流场数值模拟存在一定简化，仿真结果与实际情况存在一定偏差。另外，由于后效期膛口火药气体具有高温属性，一般借助防烧蚀涂层等技术强化钛合金装置的高温使用性能，因此未在结构强度分析中考虑热耦合效应。

图 1 试样变形前后对比

Figure 1. Comparison of deformation before and after test

下载: 全尺寸图片幻灯片

图 2 在不同应变率准静态压缩下EPS28泡沫的应力-应变曲线

Figure 2. Stress-strain curve of EPS28 foam under quasi-static compression test with different strain rates

下载: 全尺寸图片幻灯片

图 3 不同密度EPS泡沫在准静态（0.001 s^–1）压缩下的应力-应变曲线

Figure 3. Stress-strain curves of EPS foam with different densities under quasi-static compression test (0.001 s^–1)

下载: 全尺寸图片幻灯片

图 4 EPS25在不同速度冲击下的应力-应变曲线

Figure 4. Stress-strain curve of EPS25 under different impact velocities

下载: 全尺寸图片幻灯片

图 5 不同密度EPS泡沫在3.960 m/s冲击下的应力-应变曲线

Figure 5. Stress-strain curve of EPS foam with different densities under the impact velocity of 3.960 m/s

下载: 全尺寸图片幻灯片

图 6 不同密度EPS的应变率与CDIF的关系

Figure 6. Strain rate vs. CDIF for EPS with different densities

下载: 全尺寸图片幻灯片

图 7 恒定应变率下EPS28的动态压缩本构曲线

Figure 7. Dynamic compression constitutive curve of EPS28 under constant strain rate

下载: 全尺寸图片幻灯片

图 8 率相关加/卸载应力-应变曲线

Figure 8. Rate-dependent loading/unloading stress-strain curves

下载: 全尺寸图片幻灯片

图 9 各项同性硬化模型的屈服面演化

Figure 9. Evolution of yield surface for isotropic hardened models

下载: 全尺寸图片幻灯片

图 10 率相关加载应力-应变曲线及弹性卸载行为

Figure 10. Rate-dependent loading stress-strain curves and elastic unloading behavior

下载: 全尺寸图片幻灯片

图 11 EPS泡沫有限元模型变形前后对比

Figure 11. Comparison of EPS foam FEA model before and after deformation

下载: 全尺寸图片幻灯片

图 12 H=0.2 m时EPS28仿真结果与试验结果对比

Figure 12. Simulated and experimental results comparison of H=0.2 m for EPS28

下载: 全尺寸图片幻灯片

图 13 H=0.5 m时EPS28仿真结果与试验结果对比

Figure 13. Simulated and experimental results comparison of H=0.5 m for EPS28

下载: 全尺寸图片幻灯片

图 14 H=0.8 m时EPS28仿真结果与试验结果对比

Figure 14. Simulated and experimental results comparison of H=0.8 m for EPS28

下载: 全尺寸图片幻灯片

表 1 EPS落锤试验工况

Table 1. Experimental conditions of EPS

Size of specimen/（mm×mm×mm）	Density of EPS/(kg·m^–3)	Mass of drop-hammer/kg	Drop height/m
100×100×20	20，25，28	9.8648	0.2，0.5，0.8

下载: 导出CSV

表 2 EPS在20%应变处的应力与压缩动态增长因子

Table 2. Stress and CDIF of EPS at 20% strain

EPS28			EPS25			EPS20
Strain rate/s^–1	Stress/MPa	CDIF	Strain rate/s^–1	Stress/MPa	CDIF	Strain rate/s^–1	Stress/MPa	CDIF
0.001	0.212	1.000	0.001	0.178	1.000	0.001	0.144	1.000
0.01	0.244	1.147	0.01	0.197	1.102	0.01	0.159	1.100
76.50	0.282	1.325	80.25	0.235	1.319	87.00	0.170	1.180
140.50	0.314	1.478	144.00	0.251	1.408	148.50	0.182	1.262
184.75	0.352	1.657	188.25	0.266	1.492	191.00	0.199	1.379

下载: 导出CSV

表 3 EPS在30%应变处的应力与压缩动态增长因子

Table 3. Stress and CDIF of EPS at 30% strain

EPS28			EPS25			EPS20
Strain rate/s^–1	Stress/MPa	CDIF	Strain rate/s^–1	Stress/MPa	CDIF	Strain rate/s^–1	Stress/MPa	CDIF
0.001	0.235	1.000	0.001	0.201	1.000	0.001	0.162	1.000
0.01	0.269	1.144	0.01	0.222	1.104	0.01	0.177	1.095
54.75	0.305	1.298	63.25	0.261	1.302	76.00	0.195	1.201
129.00	0.335	1.427	134.75	0.273	1.358	142.00	0.207	1.280
174.75	0.371	1.580	181.50	0.285	1.422	185.75	0.221	1.362

下载: 导出CSV

表 4 基于Crushable Foam模型的EPS28材料参数

Table 4. EPS28 material parameters based on crushable foam model

Elastic parameters			Plastic parameters
Poisson’s ratio	Modulus/ MPa	Density/ (kg·m^–3)	Plastic Poisson’s ratio	Stress/ MPa	Plastic strain	Yield stress ratio	Strain rate/s^–1
0	6.213	28	0	0.168	0	1	0.001
				0.219	0.227	1.204	1
				0.265	0.498	1.273	10
				0.341	0.871	1.343	112
				0.446	1.186	1.514	150
				0.599	1.465	1.685	200

下载: 导出CSV

表 5 基于LS-DYNA MAT57模型的EPS28材料参数

Table 5. EPS28 material parameters based on LS-DYNA MAT57 model

Density/ (kg·m^–3)	Modulus/ MPa	Poisson’s ratio	Tension cut-off stress/MPa	Viscous coefficient	Shape factor	Hysteretic unloading factor	E_d/MPa	β₁
28	6.213	0	1	0.1	10	0.1	0.36	169.23

下载: 导出CSV

表 6 MAT163和Crushable Foam预测的吸收能量与试验结果对比

Table 6. MAT163 and Crushable Foam predicted results compared with the test results for absorbed energy

H/m	Absorbed energy/J			Relative error/%
H/m	Test	MAT163	Crushable Foam	MAT163	Crushable Foam
0.2	19.192	17.949	18.416	6.48	4.04
0.5	46.172	45.087	46.479	2.35	0.66
0.8	71.514	71.308	72.410	0.29	1.25

下载: 导出CSV

参考文献(19)

[1]	OUELLET S, CRONIN D, WORSWICK M. Compressive response of polymeric foams under quasi-static, medium and high strain rate conditions [J]. Polymer Testing, 2006, 25(6): 731–743. doi: 10.1016/j.polymertesting.2006.05.005
[2]	CRONIN D S, OUELLET S. Low density polyethylene, expanded polystyrene and expanded polypropylene: strain rate and size effects on mechanical properties [J]. Polymer Testing, 2016, 53: 40–50. doi: 10.1016/j.polymertesting.2016.04.018
[3]	HORVATH J S. Expanded polystyrene (EPS) geofoam: an introduction to material behavior [J]. Geotextiles and Geomembranes, 1994, 13(4): 263–280. doi: 10.1016/0266-1144(94)90048-5
[4]	GIBSON L, ASHBY M. Cellular solids: structure and properties [M]. Cambridge: Cambridge University Press, 1997.
[5]	SAINT-MICHEL F, CHAZEAU L, CAVAILLÉ J Y, et al. Mechanical properties of high density polyurethane foams: I. effect of the density [J]. Composites Science and Technology, 2006, 66(15): 2700–2708. doi: 10.1016/j.compscitech.2006.03.009
[6]	CHEN W, HAO H, HUGHES D, et al. Static and dynamic mechanical properties of expanded polystyrene [J]. Materials & Design, 2015, 69: 170–180.
[7]	王志亮, 诸斌. EPS泡沫冲击压缩和吸能特性试验研究 [J]. 建筑材料学报, 2013, 16(4): 630–636 doi: 10.3969/j.issn.1007-9629.2013.04.014 WANG Z L, ZHU B. Experimental study on impact compression and energy-absorbing property of expanded polystyrene foam [J]. Journal of Building Materials, 2013, 16(4): 630–636 doi: 10.3969/j.issn.1007-9629.2013.04.014
[8]	LING C, IVENS J, CARDIFF P, et al. Deformation response of EPS foam under combined compression-shear loading. Part I: experimental design and quasi-static tests [J]. International Journal of Mechanical Sciences, 2018, 144: 480–489. doi: 10.1016/j.ijmecsci.2018.06.014
[9]	LING C, IVENS J, CARDIFF P, et al. Deformation response of EPS foam under combined compression-shear loading. Part II: high strain rate dynamic tests [J]. International Journal of Mechanical Sciences, 2018, 145: 9–23. doi: 10.1016/j.ijmecsci.2018.06.015
[10]	LING C, CARDIFF P, GILCHRIST M D. Mechanical behaviour of EPS foam under combined compression-shear loading [J]. Materials Today Communications, 2018, 16: 339–352. doi: 10.1016/j.mtcomm.2018.07.001
[11]	SHAH Q H, TOPA A. Modeling large deformation and failure of expanded polystyrene crushable foam using LS-DYNA [J]. Modelling and Simulation in Engineering, 2014: 1.
[12]	OZTURK U E, ANLAS G. Finite element analysis of expanded polystyrene foam under multiple compressive loading and unloading [J]. Materials & Design, 2011, 32(2): 773–780.
[13]	姚小虎, 任会兰, 林荣, 等. 聚合物泡沫材料动态力学性能及其能量吸收研究 [J]. 高压物理学报, 2012, 26(5): 531–536 YAO X H, REN H L, LIN R, et al. Study on dynamic mechanical properties and energy absorption of polymeric foams [J]. Chinese Journal of High Pressure Physics, 2012, 26(5): 531–536
[14]	郭伟国, 李玉龙, 索涛. 应力波基础简明教程 [M]. 西安: 西北工业大学出版社, 2007: 120–154.
[15]	NAGY A, KO W L, LINDHOLM U S. Mechanical behavior of foamed materials under dynamic compression [J]. Journal of Cellular Plastics, 1974, 10(3): 127–134. doi: 10.1177/0021955X7401000306
[16]	ZHANG J, KIKUCHI N, LI V, et al. Constitutive modeling of polymeric foam material subjected to dynamic crash loading [J]. International Journal of Impact Engineering, 1998, 21(5): 369–386. doi: 10.1016/S0734-743X(97)00087-0
[17]	ABAQUS analysis user’s manual version 2017 [M]. ABAQUS Inc., 2017.
[18]	OZTURK U E. Mechanical behavior of low density polymeric foams under multiple loading and unloading [D]. Turkey: Bogazici University, 2008.
[19]	LS-DYNA theory manual version Vol. R7.1 [M]. Livermore: Livermore Software Technology, 2006.