冲击载荷下分层梯度泡沫材料中的应力波传播特性

赵转 李世强 刘志芳

引用本文:
Citation:

冲击载荷下分层梯度泡沫材料中的应力波传播特性

    作者简介: 赵 转(1984-),女,硕士,讲师,主要从事冲击动力学研究. E-mail:zhaozhuan001@163.com;
    通讯作者: 李世强, lishiqiang@tyut.edu.cn
  • 中图分类号: O347

Stress Waves Propagation in Layered Graded CellularMaterials under Dynamic Crush Loading

    Corresponding author: LI Shiqiang, lishiqiang@tyut.edu.cn ;
  • CLC number: O347

  • 摘要: 以一维应力波传播理论为基础,建立了3层金属泡沫材料受到刚性块撞击时的理论模型,研究了刚性块撞击3层泡沫圆杆时的动力响应过程,从理论上给出了刚性块在撞击过程中的速度衰减规律数值解。利用ANSYS/LS-DYNA分析了受刚性块撞击时塑性应力波在3层泡沫材料中的传播过程,比较了刚性块以及层间界面处节点速度的变化规律。通过对比有限元结果与理论结果发现:理论模型能够较好地预测冲击载荷下分层泡沫材料各界面的速度衰减规律;3层梯度泡沫材料比相同质量的单层均质泡沫材料具有更加高效的吸能和缓冲能力。由于理论假设忽略了反射波以及泡沫材料应变硬化效应的影响,理论解与有限元模拟结果之间存在一定的误差。
  • 图 1  刚性质量块撞击一端固定的3层泡沫材料

    Figure 1.  Impact of a rigid mass against a triple-layered foam material

    图 2  3层泡沫铝圆杆受刚性块撞击的有限元模型及3种Cymat泡沫铝压缩应力-应变曲线

    Figure 2.  Triple layer foam rod impacted by a rigid mass and stress-strain curves of three Cymat foam

    图 3  各节点速度变化规律对比

    Figure 3.  Comparison of velocity profiles

    表 1  泡沫杆参数

    Table 1.  Parameters of the foam rod

    L1/mmL2/mmL3/mm${\rho _{01}}$/(kg·m–3)${\rho _{02}}$/(kg·m–3)${\rho _{03}}$/(kg·m–3)
    304060570410253
    下载: 导出CSV
  • [1] MA G W, YE Z Q. Energy absorption of double-layer foam cladding for blast alleviation [J]. International Journal of Impact Engineering, 2007, 34(2): 329–347. doi: 10.1016/j.ijimpeng.2005.07.012
    [2] GARDNER N, SHUKLA A. The blast response of sandwich composites with a functionally graded core and polyurea interlayer [M]//Dynamic Behavior of Materials. New York: Springer, 2011: 215–223.
    [3] WANG E, GARDNER N, SHUKLA A. The blast resistance of sandwich composites with stepwise graded cores [J]. International Journal of Solids and Structures, 2009, 46(18/19): 3492–3502.
    [4] LI S, LI X, WANG Z, et al. Sandwich panels with layered graded aluminum honeycomb cores under blast loading [J]. Composite Structures, 2017, 173: 242–254. doi: 10.1016/j.compstruct.2017.04.037
    [5] LI S, LI X, WANG Z, et al. Finite element analysis of sandwich panels with stepwise graded aluminum honeycomb cores under blast loading [J]. Composites Part A: Applied Science and Manufacturing, 2016, 80: 1–12. doi: 10.1016/j.compositesa.2015.09.025
    [6] LIU X, TIAN X, LU T J, et al. Blast resistance of sandwich-walled hollow cylinders with graded metallic foam cores [J]. Composite Structures, 2012, 94(8): 2485–2493. doi: 10.1016/j.compstruct.2012.02.029
    [7] ZHANG L, HEBERT R, WRIGHT J T, et al. Dynamic response of corrugated sandwich steel plates with graded cores [J]. International Journal of Impact Engineering, 2014, 65: 185–194. doi: 10.1016/j.ijimpeng.2013.11.011
    [8] REID S R, PENG C. Dynamic uniaxial crushing of wood [J]. International Journal of Impact Engineering, 1997, 19(5/6): 531–570.
    [9] TAN P J, REID S R, HARRIGAN J J, et al. Dynamic compressive strength properties of aluminium foams. Part II—‘shock’ theory and comparison with experimental data and numerical models [J]. Journal of the Mechanics and Physics of Solids, 2005, 53(10): 2206–2230. doi: 10.1016/j.jmps.2005.05.003
    [10] KARAGIOZOVA D. Velocity attenuation and force transfer by a single-and double-layer claddings made of foam materials [J]. International Journal of Protective Structures, 2011, 2(4): 417–437. doi: 10.1260/2041-4196.2.4.417
    [11] GURUPRASAD S, MUKHERJEE A. Layered sacrificial claddings under blast loading Part I—analytical studies [J]. International Journal of Impact Engineering, 2000, 24(9): 957–973. doi: 10.1016/S0734-743X(00)00004-X
  • [1] 徐明利张若棋张光莹 . 用元胞法研究层合材料中的应力波. 高压物理学报, doi: 10.11858/gywlxb.2000.03.005
    [2] 周南 . 脉冲辐照热激波与压缩应力波. 高压物理学报, doi: 10.11858/gywlxb.1994.03.006
    [3] 马棋棋熊迅郑宇轩周风华 . 无机玻璃动态压缩破坏的离散元模拟. 高压物理学报, doi: 10.11858/gywlxb.20190719
    [4] 任善良文鹤鸣周琳 . 平头弹穿透接触式双层金属板的理论研究. 高压物理学报, doi: 10.11858/gywlxb.20170695
    [5] 孙宇新张进李永池胡时胜董杰 . 内爆加载下热塑性管壳的应力波演化与层裂效应研究. 高压物理学报, doi: 10.11858/gywlxb.2005.04.006
    [6] 张越李世强王志华吴桂英 . 二阶层级自相似四边形蜂窝动力压缩行为数值模拟. 高压物理学报, doi: 10.11858/gywlxb.2017.04.002
    [7] 夏明黄正祥顾晓辉王叶中 . 磁爆加载下强磁体的磁场研究. 高压物理学报, doi: 10.11858/gywlxb.2012.04.012
    [8] 罗鑫许金余李为民白二雷 . 变截面杆共轴撞击数值分析及其应用前景. 高压物理学报, doi: 10.11858/gywlxb.2012.06.018
    [9] 薛建锋沈培辉王晓鸣 . 弹体侵彻混凝土开坑阶段阻力的计算. 高压物理学报, doi: 10.11858/gywlxb.2016.06.010
    [10] 王鹏李世强于国际吴桂英 . 轴向冲击载荷作用下梯度圆管的动力学行为. 高压物理学报, doi: 10.11858/gywlxb.2017.06.013
    [11] 王春奎刘小苹 . LY-12铝高温凝聚态动力学性质研究高温弹性模量的测定. 高压物理学报, doi: 10.11858/gywlxb.1991.01.005
    [12] 雷忠琦姚小虎龙舒畅常建虎王海波 . 聚苯乙烯泡沫材料的动态压缩特性. 高压物理学报, doi: 10.11858/gywlxb.20180655
    [13] 庞宝君杨震琦王立闻迟润强 . 橡胶材料的动态压缩性能及其应变率相关的本构模型. 高压物理学报, doi: 10.11858/gywlxb.2011.05.005
    [14] 高光发 . 混凝土材料动态压缩强度的应变率强化规律. 高压物理学报, doi: 10.11858/gywlxb.2017.03.007
    [15] 喻寅李媛媛贺红亮王文强 . 脆性材料动态断裂的介观格子模型. 高压物理学报, doi: 10.11858/gywlxb.20190707
    [16] 李雪艳李志斌张舵 . 考虑温度效应的泡沫铝准静态压缩本构模型. 高压物理学报, doi: 10.11858/gywlxb.20170642
    [17] 姚小虎任会兰林荣张晓晴 . 聚合物泡沫材料动态力学性能及其能量吸收研究. 高压物理学报, doi: 10.11858/gywlxb.2012.05.008
    [18] 刘占芳常敬臻唐录成姚国文李建鹏 . 平面冲击波压缩下氧化铝陶瓷的动态强度. 高压物理学报, doi: 10.11858/gywlxb.2007.02.003
    [19] 王少恒杨震琦管公顺庞宝君陈海波 . 空心球/铝合金复合材料动态压缩性能与损伤破碎特性. 高压物理学报, doi: 10.11858/gywlxb.2013.03.019
    [20] 谷卓伟金孝刚张清福孙悦 . 材料预加热冲击压缩实验技术及高温下不锈钢的动态响应. 高压物理学报, doi: 10.11858/gywlxb.1998.03.005
  • 加载中
图(3)表(1)
计量
  • 文章访问数:  328
  • 阅读全文浏览量:  317
  • PDF下载量:  12
出版历程
  • 收稿日期:  2019-01-29
  • 录用日期:  2019-02-21
  • 网络出版日期:  2019-09-18

冲击载荷下分层梯度泡沫材料中的应力波传播特性

    作者简介:赵 转(1984-),女,硕士,讲师,主要从事冲击动力学研究. E-mail:zhaozhuan001@163.com
    通讯作者: 李世强, lishiqiang@tyut.edu.cn
  • 1. 运城职业技术学院建筑工程系,山西 运城 044000
  • 2. 太原理工大学机械与运载工程学院应用力学研究所,山西 太原 030024

摘要: 以一维应力波传播理论为基础,建立了3层金属泡沫材料受到刚性块撞击时的理论模型,研究了刚性块撞击3层泡沫圆杆时的动力响应过程,从理论上给出了刚性块在撞击过程中的速度衰减规律数值解。利用ANSYS/LS-DYNA分析了受刚性块撞击时塑性应力波在3层泡沫材料中的传播过程,比较了刚性块以及层间界面处节点速度的变化规律。通过对比有限元结果与理论结果发现:理论模型能够较好地预测冲击载荷下分层泡沫材料各界面的速度衰减规律;3层梯度泡沫材料比相同质量的单层均质泡沫材料具有更加高效的吸能和缓冲能力。由于理论假设忽略了反射波以及泡沫材料应变硬化效应的影响,理论解与有限元模拟结果之间存在一定的误差。

English Abstract

  • 多孔材料因其压缩应力-应变曲线包含一段较长的平台区域,在受到冲击载荷时能够吸收较多的能量,起到缓冲减震的功能,在航空航天、交通运输等领域具有广泛的应用。然而受工业生产中金属发泡工艺本身的限制,泡沫材料的厚度、平台应力以及密实化应变有限,压缩效率较低,单层泡沫材料作为吸能填充层在爆炸冲击中的缓冲和能量吸收能力受到制约[1]。分层泡沫材料在用作防护层时没有厚度限制,可以根据载荷工况进行梯度及分层优化设计[2],合理的分层泡沫组合能够对冲击波产生有效的衰减作用,提高结构的抗冲击能力[3-4]。与均质多孔材料相比,梯度和分层设计可以便捷地对不同性质的材料进行有序的组合,并结合不同材料的特点,充分发挥材料的吸能特性,使之具有更加优越的性能,因此分层梯度金属多孔材料及结构冲击动力学行为引起了广泛的关注[5-7]

    本研究以一维应力波传播理论为基础,建立了3层梯度金属泡沫材料受到刚性块撞击时的理论模型,分析各层界面速度的衰减规律,并利用ANSYS/LS-DYNA分析受刚性块撞击时应力波在3层泡沫材料中的传播过程,再与理论结果进行对比。

    • Reid等[8]最早运用波阵面上的动量守恒条件,研究了应力波在木头圆杆中的传播特性,提出了有关应力波在多孔材料中传播的刚性理想塑性锁定模型(RPPL模型)。Tan等[9]通过对冲击载荷下金属泡沫铝的动态压缩性能实验研究,提出了有关泡沫铝材料在冲击作用下的冲击模型。Karagiozova[10]针对刚性质量块撞击单层和双层梯度泡沫圆杆的动力响应及应力波传播问题,利用能量守恒原理,结合波阵面上的质量守恒及动态连续条件,分析了质量块速度变化及衰减特征。

      本研究以一维应力波传播理论为基础,建立质量为M的刚性质量块撞击一端固定的3层泡沫圆杆的理论模型,设泡沫圆杆的截面积为A0,时间为t,撞击过程如图1所示。假设整个过程分为3个相互耦合的阶段:第1阶段($0 < t \leqslant {t_1}$)为塑性冲击波通过第1层泡沫之前;第2阶段 (${t_1} < t \leqslant {t_2}$)为塑性冲击波通过第1层泡沫至通过第2层泡沫之前;第3阶段(${t_2} < t \leqslant {t_3}$)为塑性冲击波通过第2层泡沫至到达固定端之前。由于各层之间的波阻抗发生间断,压缩过程中前驱波可能导致各层之间的耦合压缩作用,这种耦合作用通过假设层间压缩区(Assumption of Compression Zone,如图1所示)来描述。理论模型中忽略泡沫材料界面间的反射作用,并且假设每层泡沫都符合RPPL模型[8]

      图  1  刚性质量块撞击一端固定的3层泡沫材料

      Figure 1.  Impact of a rigid mass against a triple-layered foam material

      假设刚性质量块的初始速度为V0,3层泡沫圆杆的初始密度分别为${\rho _{01}}$${\rho _{02}}$${\rho _{03}}$,初始长度分别为L1L2L3,初始平台屈服应力分别为$\sigma _1^{\rm{Y}}$$\sigma _2^{\rm{Y}}$$\sigma _3^{\rm{Y}}$,泊松比均假设为零。根据波阵面上的质量守恒及动量守恒定律[8, 11],得到如下结果。

      (1)第1阶段。当$0 \leqslant t \leqslant {t_1}$时,初始条件为:各界面速度${V_1}(0) = {V_0}$${V_2}(0) = 0$${V_3}(0) = 0$;各压缩区长度${h_1}(0) = 0$${h_2}(0) = 0$${h_3}(0) = 0$;各界面位移${u_1}(0) = 0$${u_2}(0) = 0$${u_3}(0) = 0$;压缩过程中假设自由端速度V1=du1/dt,第1界面处速度为V2=du2/dt,第2界面处速度为V3=du3/dt。对刚性质量块和第1层泡沫杆压缩区域(压缩长度为h1,如图1(b)所示)列动量守恒方程

      $\begin{split} & \{ M/{A_0} + {\rho _{01}}[{u_1} + {\text{δ}}{u_1} + {h_1} + {\text{δ}}{h_1} - ({u_2} + {\text{δ}}{u_2})]\left. {} \right\}({V_1} + {\text{δ}}{V_1})- \\ \, &\;\;\;\;\;\; [M/{A_0} + {\rho _{01}}({u_1} + {h_1} - {u_2})]{V_1} - {\rho _{01}}({\text{δ}}{u_1} + {\text{δ}}{h_1} - {\text{δ}}{u_2}){V_2} = - \sigma _1^{\rm{Y}}{\text{δ}}t \end{split}$

      $[M/{A_0} + {\rho _{01}}({u_1} + {h_1} - {u_2})]{\text{δ}}{V_1} + {\rho _{01}}({V_1} - {V_2})({\text{δ}}{u_1} + {\text{δ}}{h_1} - {\text{δ}}{u_2}) = - \sigma _1^{\rm{Y}}{\text{δ}}t $

      由于$\dfrac{{{\rm{d}}{h_1}}}{{{\rm{d}}t}} = {G_1} - ({V_1} - {V_2})$,冲击波波速${G_1}(\varepsilon ) = \dfrac{{{V_1} - {V_2}}}{{{\varepsilon _1}}}$${\varepsilon _1}$为第1层泡沫材料的密实化应变,所以当${\text{δ}}t \to 0$时,(2)式可写为

      $\begin{split} & [M/{A_0} + {\rho _{01}}({u_1} + {h_1} - {u_2})]\frac{{{\rm{d}}{V_1}}}{{{\rm{d}}t}} + {\rho _{01}}({V_1} - {V_2})\frac{{{\rm{d}}{u_1} + {\rm{d}}{h_1} - {\rm{d}}{u_2}}}{{{\rm{d}}t}} = [M/{A_0} + {\rho _{01}}({u_1} + {h_1} - {u_2})]\frac{{{\rm{d}}{V_1}}}{{{\rm{d}}t}}+\\ \, &\;\;\;\;{\rho _{01}}({V_1} - {V_2})\left\{ {} \right.[{G_1} - ({V_1} - {V_2})] + {V_1} - {V_2}\left. {} \right\} = [M/{A_0} + {\rho _{01}}({u_1} + {h_1} - {u_2})]\frac{{{\rm{d}}{V_1}}}{{{\rm{d}}t}} + {\rho _{01}}\frac{{{{({V_1} - {V_2})}^2}}}{{{\varepsilon _1}}} = - \sigma _1^{\rm{Y}} \end{split} $

      于是

      $\frac{{{\rm{d}}{V_1}}}{{{\rm{d}}t}} = - \left[\sigma _1^{\rm{Y}} + {\rho _{01}}\frac{{{{({V_1} - {V_2})}^2}}}{{{\varepsilon _1}}}\right]\frac{1}{{M/{A_0} + {\rho _{01}}{A_0}({u_1} + {h_1} - {u_2})}}$

      在此阶段,可以记②号界面(如图1所示)上的压力为${p_1}$,通过假设右侧第2层泡沫杆压缩段的长度为h2来考虑耦合压缩作用,对第2层泡沫杆压缩段列动量守恒方程

      ${\rho _{02}}[{h_2} + {u_2} + {\text{δ}}{h_2} + {\text{δ}}{u_2} + \left( {{u_3} + {\text{δ}}{u_3}} \right)]({V_2} + {\text{δ}}{V_2}) - {\rho _{02}}({h_2} + {u_2} - {u_3}){V_2} = ({p_1} - \sigma _2^{\rm{Y}}){\text{δ}}t$

      ${\text{δ}}t \to 0$

      ${\rho _{02}}({h_2} + {u_2} - {u_3})\frac{{{\rm{d}}{V_2}}}{{{\rm{d}}t}} + {\rho _{02}}{V_2}\frac{{{\rm{d}}({h_2} + {u_2} - {u_3})}}{{{\rm{d}}t}} = {p_1} - \sigma _2^{\rm{Y}}$

      同样有如下关系:$\dfrac{{{\rm{d}}{h_2}}}{{{\rm{d}}t}} = {G_2} - \left( {{V_2} - {V_3}} \right)$,冲击波波速${G_2}(\varepsilon ) = \dfrac{{{V_2} - {V_3}}}{{{\varepsilon _2}}}$${\varepsilon _2}$为第2层泡沫材料的密实化应变,${V_3} = \dfrac{{{\rm{d}}{u_3}}}{{{\rm{d}}t}}$。代入(6)式,可得

      ${\rho _{02}}({h_2} + {u_2} - {u_3})\frac{{{\rm{d}}{V_2}}}{{{\rm{d}}t}} = {p_1} - (\sigma _2^{\rm{Y}} + {\rho _{02}}{V_2}{G_2})$

      由于惯性效应,压缩区动态应力符合如下关系

      $\sigma _2^{\rm d} = \sigma _2^{\rm{Y}} + {\rho _{02}}\left( {{V_2} - {V_3}} \right){G_2}$

      联合(7)式和(8)式,得到

      ${\rho _{02}}({h_2} + {u_2} - {u_3})\frac{{{\rm{d}}{V_2}}}{{{\rm{d}}t}} = {p_1} - \sigma _2^{\rm d}$

      对第1层泡沫杆未压缩段列动量守恒方程

      ${\rho _{01}}[{L_1} - ({h_1} + {u_1} - {u_2})]{\text{δ}}{V_2} - {\rho _{01}}({V_1} + {\text{δ}}{V_1})({\text{δ}}{h_1} + {\text{δ}}{u_1} - {\text{δ}}{u_2}) + {\rho _{01}}{V_2}({\text{δ}}{h_1} + {\text{δ}}{u_1} - {\text{δ}}{u_2})= (\sigma _1^{\rm d} - {p_1}){\text{δ}}t $

      略去无穷小量$({\text{δ}}{h_1}{\text{δ}}{V_1} + {\text{δ}}{u_1}{\text{δ}}{V_1} - {\text{δ}}{u_2}{\text{δ}}{V_1})$,则

      $ {\rho _{01}}[{L_1} - ({h_1} + {u_1} - {u_2})]\frac{{{\rm{d}}{V_2}}}{{{\rm{d}}t}} - {\rho _{01}}{V_1}\frac{{{\rm{d}}({h_1} + {u_1} - {u_2})}}{{{\rm{d}}t}}+ {\rho _{01}}{V_2}\frac{{{\rm{d}}({h_1} + {u_1} - {u_2})}}{{{\rm{d}}t}} = \sigma _1^{\rm d} - {p_1} $

      式中:$\sigma _1^{\rm d} = \sigma _1^{\rm{Y}} + {\rho _{01}}({V_1} - {V_2}){G_1}$。由此可得

      $\frac{{{\rm{d}}{V_2}}}{{{\rm{d}}t}} = - \left[\sigma _2^{\rm{Y}} + {\rho _{02}}\frac{{{{\left( {{V_2} - {V_3}} \right)}^2}}}{{{\varepsilon _2}}} - \sigma _1^{\rm{Y}}\right] \cdot \frac{1}{{{\rho _{02}}({h_2} + {u_2} - {u_3}) + {\rho _{01}}[{L_1} - ({u_1} + {h_1} - {u_2})]}}$

      同样,在此阶段记③号界面上的压力为${p_2}$,通过假设右侧第3层泡沫杆压缩段的长度为h3来考虑耦合压缩作用,对第3层泡沫杆压缩段列动量守恒方程

      ${\rho _{03}}({h_3} + {u_3} + {\text{δ}}{h_3} + {\text{δ}}{u_3})({V_3}_{} + {\text{δ}}{V_3}) - {\rho _{03}}({h_3} + {u_3}){V_3} = ({p_2} - \sigma _3^{\rm{Y}}){\text{δ}}t$

      ${\text{δ}}t \to 0$

      ${\rho _{03}}({h_3} + {u_3})\frac{{{\rm{d}}{V_3}}}{{{\rm{d}}t}} + {\rho _{03}}{V_3}\frac{{{\rm{d}}({h_3} + {u_3})}}{{{\rm{d}}t}} = {p_2} - \sigma _3^{\rm{Y}}$

      因为$\dfrac{{{\rm{d}}{h_3}}}{{{\rm{d}}t}} = {G_3} - {V_3}$,冲击波波速${G_3}(\varepsilon ) = \dfrac{{{V_3}}}{{{\varepsilon _3}}}$${\varepsilon _3}$为第3层泡沫材料的密实化应变,${V_3} = \dfrac{{{\rm{d}}{u_3}}}{{{\rm{d}}t}}$,则

      $\sigma _3^{\rm d} = \sigma _3^{\rm{Y}} + {\rho _{03}}{V_3}{G_3}$

      联合(14)式和(15)式,有

      ${\rho _{03}}({h_3} + {u_3})\frac{{{\rm{d}}{V_3}}}{{{\rm{d}}t}} = {p_2} - \sigma _3^{\rm d}$

      对第2层泡沫杆未压缩段列动量守恒方程

      ${\rho _{02}}[{L_2} - ({h_2} + {u_2} - {u_3})]{\text{δ}}{V_3} - {\rho _{02}}({V_2} + {\text{δ}}{V_2})({\text{δ}}{h_2} + {\text{δ}}{u_2} - {\text{δ}}{u_3})+ {\rho _{02}}{V_3}({\text{δ}}{h_2} + {\text{δ}}{u_2} - {\text{δ}}{u_3}) = (\sigma _2^{\rm d} - {p_2}){\text{δ}}t $

      略去无穷小量$({\text{δ}}{h_1}{\text{δ}}{V_1} + {\text{δ}}{u_1}{\text{δ}}{V_1} - {\text{δ}}{u_2}{\text{δ}}{V_1})$,则

      ${\rho _{02}}[{L_2} - ({h_2} + {u_2} - {u_3})]\frac{{{\rm{d}}{V_3}}}{{{\rm{d}}t}} - {\rho _{02}}{V_2}\frac{{{\rm{d}}({h_2} + {u_2} - {u_3})}}{{{\rm{d}}t}} + {\rho _{02}}{V_3}\frac{{{\rm{d}}({h_2} + {u_2} - {u_3})}}{{{\rm{d}}t}}= \sigma _2^{\rm d} - {p_2} $

      所以有

      $\frac{{{\rm{d}}{V_3}}}{{{\rm{d}}t}} = - \left(\sigma _3^{\rm{Y}} + {\rho _{03}}\frac{{{V_3^2}}}{{{\varepsilon _3}}} - \sigma _2^{\rm{Y}}\right) \cdot \frac{1}{{{\rho _{03}}({h_3} + {u_3}) + {\rho _{02}}[{L_2} - ({u_2} + {h_2} - {u_3})]}}$

      (2)第2阶段。${t_1}< t \leqslant {t_2}$, ${V_1} = {V_2}$$\dfrac{{{\rm{d}}{h_2}}}{{{\rm{d}}t}} = {G_2} - ({V_2} - {V_3})$${G_2}(\varepsilon ) = \dfrac{{{V_2} - {V_3}}}{{{\varepsilon _2}}}$$\dfrac{{{\rm{d}}{u_2}}}{{{\rm{d}}t}} = {V_2}$$\dfrac{{{\rm{d}}{u_3}}}{{{\rm{d}}t}} = {V_3}$,根据波阵面上的动量守恒条件

      $\begin{split} & \left\{\right.M/{A_0} + {\rho _{01}}{L_1} + {\rho _{02}}[{u_2} + {\text{δ}}{u_2} + {h_2} + {\text{δ}}{h_2} - ({u_3} + {\text{δ}}{u_3})]\left. {} \right\}({V_2} + {\text{δ}}{V_2}) - \\ &\;\;\; [M/{A_0} + {\rho _{01}}{L_1} + {\rho _{02}}({u_2} + {h_2} - {u_3})]{V_2} - {\rho _{02}}({\text{δ}}{u_2} + {\text{δ}}{h_2} - {\text{δ}}{u_3}){V_3} = - {\sigma _{_{{\rm Y}2}}}{\text{δ}}t \end{split} $

      $ [M/{A_0} + {\rho _{01}}({u_1} + {h_1} - {u_2}){\text{δ}}{V_1}+ {\rho _{01}}({V_1} - {V_2})({\text{δ}}{u_1} + {\text{δ}}{h_1} - {\text{δ}}{u_2}) = - {\sigma _{_{{\rm Y}1}}}{\text{δ}}t $

      $\begin{split} & [M/{A_0} + {\rho _{01}}{L_1} + {\rho _{02}}({u_2} + {h_2} - {u_3})]\frac{{{\rm{d}}{V_2}}}{{{\rm{d}}t}} + {\rho _{02}}({V_2} - {V_3})\frac{{{\rm{d}}{u_2} + {\rm{d}}{h_2} - {\rm{d}}{u_3}}}{{{\rm{d}}t}} = \\ \; & [M/{A_0} + {\rho _{01}}{L_1} + {\rho _{02}}({u_2} + {h_2} - {u_3})]\frac{{{\rm{d}}{V_2}}}{{{\rm{d}}t}} + {\rho _{02}}({V_2} - {V_3})\left\{ {} \right.[{G_2} - ({V_2} - {V_3})] + {V_2} - {V_3}\left. {} \right\} = \\ \; & [M/{A_0} + {\rho _{01}}{L_1} + {\rho _{02}}({u_2} + {h_2} - {u_3})]\frac{{{\rm{d}}{V_2}}}{{{\rm{d}}t}} + {\rho _{02}}\frac{{{{({V_2} - {V_3})}^2}}}{{{\varepsilon _2}}} = - {\sigma _{_{{\rm Y}2}}}\end{split} $

      可得

      $\frac{{{\rm{d}}{V_2}}}{{{\rm{d}}t}} = - \left[{\sigma _{_{{\rm Y}2}}}+ {\rho _{02}}\frac{{{{({V_2} - {V_3})}^2}}}{{{\varepsilon _2}}}\right]\frac{1}{{M/{A_0} + {\rho _{01}}{L_1} + {\rho _{02}}{A_0}({u_2} + {h_2} - {u_3})}}$

      $\frac{{{\rm{d}}{V_3}}}{{{\rm{d}}t}} = - \left(\sigma _3^{\rm{Y}} + {\rho _{03}}\frac{{{V_3}^2}}{{{\varepsilon _3}}} - \sigma _2^{\rm{Y}}\right) \cdot \frac{1}{{{\rho _{03}}({h_3} + {u_3}) + {\rho _{02}}[{L_2} - ({u_2} + {h_2} - {u_3})]}}$

      其中初始条件:${V_2}({t_1}) = {V_1}({t_1})$${h_2}({t_1}) = 0$${h_3}({t_1}) = 0$${u_2}({t_1}) = {\bar u_2}$${u_3}({t_1}) = {\bar u_3}$

      (3)第3阶段。${t_2} < t \leqslant {t_3}$${V_1} = {V_2} = {V_3}$

      $\frac{{{\rm{d}}{V_3}}}{{{\rm{d}}t}} = - \left(\sigma _3^{\rm{Y}} + {\rho _{03}}\frac{{{V_3^2}}}{{{\varepsilon _3}}}\right)\frac{1}{{M/{A_0} + {\rho _{01}}{L_1} + {\rho _{02}}{L_2} + {\rho _{03}}({u_3} + {h_3})}}$

      $\frac{{{\rm{d}}{h_3}}}{{{\rm{d}}t}} = {G_3} - {V_3},\;\;\;\;{G_3}(\varepsilon ) = \frac{{{V_3}}}{{{\varepsilon _3}}},\;\;\;\;\frac{{{\rm{d}}{u_3}}}{{{\rm{d}}t}} = {V_3}$

      初始条件:${V_3}({t_2}) = {V_2}({t_2})$${h_3}({t_2}) = 0$${u_3}({t_2}) = {\overline{\overline u} _3}$

      利用Runge-Kutta方法,将各阶段初始条件代入(4)式、(12)式、(19)式、(23)式、(24)式和(25)式,可求得各界面速度随时间变化规律的数值解。

    • 为了验证理论模型的可靠性,建立如图2所示的有限元模型(1/4模型)。刚性块的质量M为370 g,泡沫杆直径为100 mm,3层泡沫杆的参数如表1所示。

      图  2  3层泡沫铝圆杆受刚性块撞击的有限元模型及3种Cymat泡沫铝压缩应力-应变曲线

      Figure 2.  Triple layer foam rod impacted by a rigid mass and stress-strain curves of three Cymat foam

      L1/mmL2/mmL3/mm${\rho _{01}}$/(kg·m–3)${\rho _{02}}$/(kg·m–3)${\rho _{03}}$/(kg·m–3)
      304060570410253

      表 1  泡沫杆参数

      Table 1.  Parameters of the foam rod

      在有限元计算中,泡沫材料通过Crushable Foam模型模拟其压缩变形,图2给出了3种不同密度的Cymat泡沫铝的准静态压缩应力-应变曲线。为了记录压缩过程中各界面的速度变化规律,分别在刚性块与第1层泡沫杆及中间泡沫杆界面处选择节点1、2、3,记录压缩过程中层间界面的速度时程曲线。

      图3对比了有限元模拟结果和理论模型计算结果,并且与相同质量的单层泡沫圆杆(单层泡沫杆长度L1=200 mm,密度$\rho $=253 kg/m3)的理论和有限元结果进行了比较。从图3可以看出,理论结果与有限元模拟结果吻合较好,说明理论模型能够较好地预测冲击载荷下分层泡沫材料各界面的速度衰减规律。

      图  3  各节点速度变化规律对比

      Figure 3.  Comparison of velocity profiles

      图3可知,在整体质量一定的情况下,增加泡沫材料的层数时,理论结果与有限元模拟结果之间的相对偏差增加。这是由于:理论计算中将连续的压缩过程假设为相互耦合的3个阶段,各阶段之间的计算误差不断累积;并且在有限元计算中泡沫铝材料的应力-应变关系为单轴压缩实验曲线,随着应变的增加有明显的硬化现象,而理论计算并未考虑应变硬化的影响,压缩过程中平台应力保持恒定,所以在整个压缩过程中刚性块速度的有限元结果略高于理论结果。从图3中还可以看出,在理论模型中第1阶段和第2阶段对应的时间分别为t1=0.16 ms和t2=0.33 ms,数值模拟中则分别为t1=0.15 ms和t2=0.39 ms。在3层泡沫材料的有限元结果中,节点2的速度与节点1速度相等之后会略微超越节点1的速度,节点3也有类似情况,这是由于:冲击波在截面处发生反射,反射波抑制了波阵面后方节点的运动,而在理论计算中并未考虑界面处反射波的作用。

      比较单层和3层泡沫材料节点1的结果发现:在相同的冲击速度下,3层梯度泡沫材料对刚性块的衰减作用更明显,且相同质量的3层泡沫材料的总长度更短,所以分层梯度泡沫材料具有更高的压缩效率[11]以及更有效的吸能和缓冲效果。

    • 在理论分析刚性质量块撞击分层梯度泡沫材料的动力响应过程的基础上,利用ANSYS/LS-DYNA建立了相应的有限元模型,分析了受刚性块撞击时应力波在3层泡沫材料中的传播过程,比较了层间界面的速度变化规律。根据波阵面上动量守恒关系与一维应力波传播理论,建立了刚性块撞击分层泡沫材料时各界面速度衰减规律的理论模型,利用Runge-Kutta方法给出了理论模型的数值解,理论计算结果与有限元模拟结果吻合较好,说明理论模型能够较好地预测冲击载荷下分层梯度泡沫材料各界面的速度衰减规律。然而,由于在理论计算中忽略了反射波以及泡沫材料应变硬化效应的影响,理论解与有限元模拟结果之间存在一定偏差。从有限元结果还发现,在相同质量和冲击速度下,与单层泡沫材料相比,3层梯度泡沫材料具有更高效的吸能和缓冲能力。

参考文献 (11)

目录

    /

    返回文章
    返回