Dynamic crushing and energy absorption of regular,irregular and functionally graded cellular structures

The in-plane dynamic crushing of two dimensional honeycombs with both regular hexagonal and irregular arrangements was investigated using detailed finite element models. The energy absorption of honeycombs made of a linear elastic-perfectly plastic material with constant and functionally graded density were estimated up to large crushing strains. Our numerical simulations showed three distinct crushing modes for honeycombs with a constant relative density: quasi-static, transition and dynamic. Moreover, irregular cellular structures showed to have energy absorption similar to their counterpart regular honeycombs of same relative density and mass. To study the dynamic crushing of functionally graded cellular structures, a density gradient in the direction of crushing was introduced in the computational models by a gradual change of the cell wall thickness. Decreasing the relative density in the direction of crushing was shown to enhance the energy absorption of honeycombs at early stages of crushing. The study provides new insight into the behavior of engineered and biological cellular materials, and could be used to develop novel energy absorbent structures.     

具有恒定冲击载荷的梯度泡沫金属材料设计

多胞材料可通过大变形大量地吸收冲击能量，引入密度梯度可进一步提高其耐撞性。梯度多胞材料的宏观力学响应对材料密度分布极为敏感，不同类型的细观构型的影响也极为不同。已有的研究工作主要局限在对给定的密度梯度分析其动态响应，较少对耐撞性设计方法进行研究。本文针对梯度闭孔泡沫金属材料，基于非线性塑性冲击波模型发展了耐撞性反向设计方法，以维持冲击物受载恒定为目标，运用级数法获得了简化模型和渐近解。利用变胞元尺寸法构建了连续梯度变化的三维Voronoi细观有限元模型，并利用ABAQUS/Explicit有限元软件对理论设计进行数值验证。结果表明，反向设计理论简化模型的渐近解对于梯度闭孔泡沫金属材料的耐撞性设计是有效的，所提出的耐撞性设计方法在控制冲击吸能过程和冲击物受载方面具有指导意义。     

梯度多胞牺牲层的抗爆炸分析

运用一维非线性塑性冲击波模型和细观有限元模型对密度梯度多胞牺牲层的抗爆炸性能进行了分析。基于率无关的刚性-塑性硬化模型,建立了描述冲击波在多胞牺牲层中传播的控制方程,分别给出了正、负密度梯度多胞材料在指数型爆炸载荷作用下的响应特性。研究了可正好吸收爆炸能量的梯度多胞牺牲层的临界厚度与载荷强度、覆盖层质量、多胞材料的密度梯度等参数之间的关系,给出了以临界厚度和支撑端应力峰值为指标的密度梯度设计图。运用二维细观有限元模型验证了基于非线性塑性冲击波模型的抗爆炸分析的有效性。     

负梯度闭孔泡沫金属的力学性能分析

运用三维Voronoi技术生成闭孔梯度泡沫模型,结合有限元分析方法模拟负梯度闭孔泡沫金属在不同冲击速度下的力学行为。结果表明,随着冲击速度的提高,得到了与均匀泡沫一样的三种变形模式：准静态模式,过渡模式和冲击模式。通过对名义应力应变曲线和变形模式的研究,提出了一种新的定义局部密实化应变的方法,并研究了相对密度和密度梯度对它的影响。分别建立了相对密度和密度梯度与冲击速度的变形模式图。通过引入密实化因子,确定了三种变形模式对应的临界冲击速度。最后讨论了不同冲击速度下,密度梯度大小对泡沫材料能量吸收能力的影响。结果表明,在高速冲击的变形初期,密度梯度的绝对值越大,泡沫材料的能量吸收能力越强。     

Closed-form solution for shock wave propagation in density-graded cellular material under impact

Density-graded cellular materials have tremendous potential in structural applications where impact resistance is required. Cellular materials subjected to high impact loading result in a compaction type deformation, usually modeled using continuum-based shock theory. The resulting governing differential equation of the shock model is nonlinear, and the density gradient further complicates the problem. Earlier studies have employed numerical methods to obtain the solution. In this study, an analytical closed-form solution is proposed to predict the response of density-graded cellular materials subjected to a rigid body impact. Solutions for the velocity of the impinging rigid body mass, energy absorption capacity of the cellular material, and the incident stress are obtained for a single shock propagation. The results obtained are in excellent agreement with the existing numerical solutions found in the literature. The proposed analytical solution can be potentially used for parametric studies and for effectively designing graded structures to mitigate impact.     

多孔材料吸能行为对相对密度和冲击速度的依赖性

多孔材料是一种优异的吸能缓冲材料,但由于其变形模式的非单一性以及动态应力应变曲线的难获取性,其吸能行为对相对密度和冲击速度的依赖性关系还并不完全明朗。本文基于不需要提前作本构假定的波传播法,开展了多孔材料的吸能行为研究。采用多孔材料的细观有限元模型进行Taylor冲击虚拟实验,获取全场质点速度时程曲线,结合Lagrange分析法得到多孔材料的局部应力应变信息,进而探讨了动态吸能性能对材料相对密度和冲击速度的依赖性。研究结果表明多孔材料的吸能行为可依据变形模式分为三个阶段。在冲击模式下,多孔材料单位体积吸能与相对密度成线性增加关系,此时惯性起主导作用;在过渡模式下,惯性的主导作用减弱,单位体积吸能量的增加速率随相对密度的增加而减弱;在准静态模式下,多孔材料只能发生微小的变形,其吸能很少。本文进一步获得了区别于多孔材料准静态应力-应变曲线的动态应力-应变状态曲线,并考察了其与相对密度之间的关系。结果表明：随着相对密度的增加,多孔材料的动态压实应变将变小,而动态塑性平台应力将提高。     

冲击载荷下径向密度排布对泡沫金属力学性能影响的研究

参照层状密度梯度泡沫模型实现方法，利用3D-Voronoi技术设计了新型径向密度梯度泡沫模型，并用有限元软件，对它在不同冲击载荷下的力学行为进行数值模拟。研究冲击速度、密度梯度和平均相对密度对金属泡沫冲击端、支撑端应力和能量吸收能力的影响，发现:径向正梯度泡沫与层状正、负梯度泡沫相比，其两端的应力值均较小，可同时保护冲击端、支撑端物体；径向负梯度泡沫两端应力变化幅度较小，能够保证物体受力稳定；几种泡沫金属的能量吸收能力在不同冲击速度下发生交替变化。对于径向梯度泡沫，能量吸收能力对密度梯度大小不敏感，对梯度方向敏感，径向负梯度泡沫的能量吸收能力始终大于径向正梯度泡沫；平均相对密度越大，径向正、负梯度泡沫两端应力越大、吸能效果越好。     

多孔泡沫牺牲层的动态压溃及缓冲吸能机理研究

本文对强动载荷下多孔泡沫牺牲层的动态压溃行为及缓冲吸能机理进行了研究. 基于刚性-理想塑性-锁定(R-PP-L)及刚性-塑性硬化(R-PH)两类多孔泡沫材料本构, 建立了强动载荷下多孔泡沫牺牲层动态响应的理论分析模型, 分析了一维冲击波在多孔泡沫牺牲层中的传播规律; 利用Voronoi方法建立了多孔泡沫牺牲层的二维细观有限元模型, 获得了冲击载荷下多孔泡沫牺牲层的变形模式和动态响应曲线, 讨论了多孔泡沫材料的层间界面效应对多孔泡沫牺牲层缓冲吸能的影响. 研究结果表明, 考虑多孔泡沫材料塑性硬化影响的理论分析模型能够预测入射波在远端的反射及对多孔泡沫牺牲层的二次压缩过程和端部应力增强现象; 相比较存在界面的多孔泡沫牺牲层, 连续设计的多孔泡沫牺牲层可增强其缓冲吸能能力, 但在界面处增加设计刚性面板则能够降低界面胞元不完整对缓冲吸能的影响; 相同冲量载荷下, 端部应力峰值随冲击能量增大而增大, 而端部冲击波的反射可能是端部应力增强的主要诱因.      

Propagation of compaction waves in metal foams exhibiting strain hardening

One-dimensional models for compaction of cellular materials exhibiting strain hardening are proposed for two different impact scenarios. The models reveal the characteristic features of deformation under the condition of decreasing velocity during the impact event. It was established that an unloading plastic wave of strong discontinuity propagates in the foam and it has a significant dynamic effect on the foam compaction and energy absorption. The proposed models are based on the actual experimentally derived stress strain curves. The compaction mechanism in three aluminium based foam materials, two of them with relatively low density – Alporas and Cymat with 9% and 9.3% relative density, respectively and a higher density Cymat foam with 21% relative density, is analysed. Numerical simulations were carried out to verify the proposed models.The predictions of the proposed models are compared with published analytical models of compaction of cellular materials which assume a predefined densification strain. It is shown that the approximation of a cellular material with significant strain hardening by the Rigid Perfectly-Plastic-Locking (RPPL) model can lead to an overestimation of the energy absorption capacity for the observed stroke due to the non-uniform strains along the compacted zone of the actual material in contrast to the predefined constant densification strain in the RPPL model. The assumption of a constant densification strain leads also to an overestimation of the maximum stress, which occurs under impact.     

多孔金属夹层板在冲击载荷作用下的动态响应

借助两种有限元软件ABAQUS和LS_DYNA, 模拟和分析了两种厚度不同的泡沫铝合金夹层板(三明治板)、方孔蜂窝形夹层板和波纹形夹层板在冲击载荷下的动态响应. 4种夹层板的单位面积密度相同,冲击载荷分别用泡沫铝子弹与不锈钢子弹模拟. 讨论了泡沫金属夹层板和格构式夹层板在不同冲击载荷作用下的变形机制,重点在于对夹层板的吸能特性及板内各部分吸能变化规律的探讨.研究结果表明: 在泡沫子弹冲击下,夹层板主要是通过自身变形来消耗子弹动能,并转化为自身内能. 厚度为22\,mm的泡沫金属夹层板吸收能量最多,底面变形最小,是结构性能最优的夹层板;在刚性子弹高速冲击穿透过程中,格构式夹层板的吸能性能比单位面积密度相同的泡沫金属夹层板的吸能性能更好. 波纹形夹层板的能量吸收能力在4种板中最高.      

应用CBR确定面心立方金属的表面弹性参量

首先从能量变分出发基于同时考虑应变梯度效应和表面效应的跨尺度力学理论, 推导出表面能和表面弹性本构等基本关系, 然后基于简单的准连续Cauchy-Born法则(CBR)建立一种确定表面能密度以及表面弹性参量的方法.进一步以面心立方(face-centre-cubic,FCC)金属为例, 系统地获得了常用FCC金属表面弹性参量的数值, 结果与他人应用分子动力学计算得到的结果相吻合.     

砂墙吸能效应对化爆相似律影响的数值研究

同介质中同能量密度的化学炸药爆炸满足爆炸相似律,在此类化爆实验中,缩比实验作为一种经济有效的实验方法被广泛采用。一些研究人员进行化爆缩比实验时,为了保证实验安全,通常会在实验中加入砂墙,起到吸收爆炸能量,消减冲击波的作用。为了研究砂墙消波吸能效应对化爆相似律的影响,建立了同介质中同能量密度化学炸药爆炸的一维计算模型,并利用两相流程序对三种缩比条件下(1∶1/1∶5/1∶10)的物理模型进行了数值模拟。数值计算结果显示:由于砂墙的吸能效应,化爆不再满足相似律。     

Investigation of the Nitrogen Shock Wave Structure on the Basis of Trajectory Calculations of the Molecular Interaction

The shock wave structure in a diatomic gas is investigated using the direct statistical simulation (Monte-Carlo) method. The energy exchange between translational and rotational degrees of freedom (TR-exchange) is calculated by solving the dynamic problem of the interaction between rigid-rotator molecules within the framework of classical mechanics. The density profiles calculated are compared with the experimental data and on this basis the nitrogen rotational relaxation time is estimated. The possibility of using simplified intermolecular interaction models, namely, the variable-diameter sphere model employed together with a phenomenological consideration of the TR-exchange, is studied. Gasdynamic parameter profiles in the shock wave are analyzed. Simple approximations of the velocity gradient and translational and rotational temperature profiles are obtained on the basis of a parametric calculation of the shock wave structure. This makes it possible approximately to describe the gasdynamic parameter profiles in terms of elementary functions.     

The quasi-stationary structure of radiating shock waves. I. The one-temperature fluid

We calculate the quasi-stationary structure of a radiating shock wave propagating through a spherically symmetric shell of cold gas by solving the time-dependent equations of radiation hydrodynamics on an implicit adaptive grid. We show that this code successfully resolves the shock wave in both the subcritical and supercritical cases and, for the first time, we have reproduced all the expected features – including the optically thin temperature spike at a supercritical shock front – without invoking analytic jump conditions at the discontinuity. We solve the full moment equations for the radiation flux and energy density, but the shock wave structure can also be reproduced if the radiation flux is assumed to be proportional to the gradient of the energy density (the diffusion approximation), as long as the radiation energy density is determined by the appropriate radiative transfer moment equation. We find that Zel'dovich and Raizer's (1967) analytic solution for the shock wave structure accurately describes a subcritical shock but it underestimates the gas temperature, pressure, and the radiation flux in the gas ahead of a supercritical shock. We argue that this discrepancy is a consequence of neglecting terms which are second order in the minimum inverse shock compression ratio [, where is the adiabatic index] and the inaccurate treatment of radiative transfer near the discontinuity. In addition, we verify that the maximum temperature of the gas immediately behind the shock is given by , where is the gas temperature far behind the shock. Received 21 September 1998/ Accepted 2 February 1999     

多孔金属及其夹芯结构力学性能的研究进展

高孔隙率多孔金属及其夹芯复合结构是一种物理功能与结构一体化的新型、轻质高强材料/结构,具有高比强度、高比刚度和优良的吸能和缓冲性能等多种功能,引起了学术界和工程界众多研究者的极大关注. 本文概述了轻质多孔金属及其夹芯结构的制备方法、多功能特性及其应用,介绍了多孔金属夹芯结构元件(梁、板、壳)遭受准静态和动态冲击载荷下的理论、实验和模拟方面的国内外研究现状,分析和讨论了多孔金属及其夹芯结构力学行为研究中的研究手段和基本问题,重点关注了多孔金属夹芯结构的变形/失效、动态响应和能量吸收.     

On the energy conservation and critical velocities for the propagation of a ＂steady-shock＂ wave in a bar made of cellular material

The propagation of shock waves in a cellular bar is systematically studied in the framework of continuum solids by adopting two idealized material models, viz. the dynamic rigid, perfectly plastic, locking （D-R-PP-L） model and the dynamic rigid, linear hardening plastic, locking （D-R-LHP-L） model, both considering the effects of strain-rate on the material properties. The shock wave speed relevant to these two models is derived. Consider the case of a bar made of one of such material with initial length L 0 and initial velocity v i impinging onto a rigid target. The variations of the stress, strain, particle velocity, specific internal energy across the shock wave and the cease distance of shock wave are all determined analytically. In particular the ＂energy conservation condition＂ and the ＂kinematic existence condition＂ as proposed by Tan et al. （2005） is re-examined, showing that the ＂energy conservation condition＂ and the consequent ＂critical velocity＂, i.e. the shock can only be generated and sustained in R-PP-L bars when the impact velocity is above this critical velocity, is incorrect. Instead, with elastic deformation, strain-hardening and strain-rate sensitivity of the cellular materials being considered, it is appropriate to redefine a first and a second critical impact velocity for the existence and propagation of shock waves in cellular solids. Starting from the basic relations for shock wave propagating in D-R-LHP-L cellular materials, a new method for inversely determining the dynamic stress-strain curve for cellular materials is proposed. By using e.g. a combination of Taylor bar and Hopkinson pressure bar impact experimental technique, the dynamic stress-strain curve of aluminum foam could bedetermined. Finally, it is demonstrated that this new formulation of shock theory in this one-dimensional stress state can be generalized to shocks in a one-dimensional strain state, i.e. for the case of plate impact on cellular materials, by simply making proper replacements of the elastic and plastic constants.     

航行体梯度密度式头帽结构设计及降载性能分析

针对航行体在以大于100 m/s的速度高速入水过程中承受巨大的冲击载荷可能导致的结构损坏、弹道失控等现象, 而现有的缓冲措施降载能力有限的难题, 本文设计了一种航行体高速入水梯度密度式缓冲头帽, 确保航行体能够高速安全入水, 并给出了详细的设计过程. 同时基于ALE (arbitrary Lagrangian-Eulerian)算法建立了航行体带缓冲头帽高速入水数值计算模型, 且数值计算的结果与试验测试数据具有较好的一致性. 然后在此基础上, 开展了航行体带梯度密度式缓冲头帽高速入水降载特性的数值研究, 探究了双层缓冲件不同分层厚度、正负密度梯度排列以及层间密度差等重要参数对缓冲头帽能量吸收以及缓冲降载效果的影响规律, 并进行了大尺度模型高速入水冲击测试试验, 根据航行体模型干模态分析时的二阶弯曲模态固有频率对试验数据进行滤波处理. 研究结果表明, 在本文所研究的范围内, 分层的缓冲件相比较于不分层的缓冲件表现出更强的冲击能量吸收效果, 且缓冲件吸收的冲击能量随着分层数的增加而增加; 负密度梯度排列的缓冲件其缓冲能力强于正密度梯度的缓冲件; 当层间密度差越大时, 冲击能量的损耗也将越大, 缓冲头帽的降载效果越好.      

内爆炸载荷下梯度泡沫铝夹芯管的动态响应

基于3D-Voronoi技术构建了泡沫铝芯层的三维细观有限元模型，对梯度泡沫铝夹芯管在内爆炸载荷下的动态响应进行了数值模拟。分析讨论了夹芯管结构内外管的壁厚、泡沫芯层的相对密度、芯层梯度分布等参数对夹芯管结构的抗爆性能与吸能性能的影响，并与无芯层的双层圆管进行了对比。结果表明:泡沫材料的相对密度可通过改变泡沫胞元大小和胞元壁厚进行调控，利用两种方式构建的夹芯管计算结果一致；保持内、外圆管总质量不变，增大内管壁厚可以有效减小外管的塑性变形，但会影响泡沫芯层的能量耗散；泡沫芯层的填充可以有效降低内管的塑性变形，正梯度泡沫铝夹芯管的抗爆性能优于均匀泡沫及负梯度泡沫夹芯管。     

高温非平衡流动中的氧分子振动态精细分析

高超声速流动在头激波压缩后常处于高 温条件下的热化学非平衡状态. 本文采用态-态方法和双温度模型计算分析了一维正激波后和高超声速钝体绕流驻点线上的氧气热化学非平衡流动. 态-态方法将氧气的每个振动能级当成独立的组分,通过耦合 Euler 方程或驻点线上的降维 Navier-Stokes 方程,数值求解得 到了高温流动中的精细热化学非平衡状态. 而双温度模型假设氧气的振动能级服从 Boltzmann 分布,通过求解振动能方程得到振动温度. 一维正激波后热化学松弛过程的计算结果表明,态-态计算预测的温度分布和氧原子浓度分布较好地吻合了文 献中的实验结果,而经典的双温度模型的预测结果误差较大,且不同双温度模型的计算结果比较发散. 态-态方法详细地给出了所有振动能级的变化过程. 无论是正激波还是脱体激波后的流场,都是高振动能级首先得到激发;但是数密度大的低振动能级先达到热平衡,而高能级 分子要经过很长距离后才能达到热平衡. 在驻点附近,复合反应生成的氧气分子处于高振动能级,导致高振动能级分子数密度显著高于平衡分布. 计算还发现,经典双温度模型的离解反应速率明显偏离态-态计算结果,无法准确体现振动离解耦合效应对离解反应 速率的影响,但是 Park 双温度模型将离解失去的振动能取为 0.3$\sim 高超声速流动在头激波压缩后常处于高 温条件下的热化学非平衡状态. 本文采用态-态方法和双温度模型计算分析了一维正激波后和高超声速钝体绕流驻点线上的氧气热化学非平衡流动. 态-态方法将氧气的每个振动能级当成独立的组分,通过耦合 Euler 方程或驻点线上的降维 Navier-Stokes 方程,数值求解得 到了高温流动中的精细热化学非平衡状态. 而双温度模型假设氧气的振动能级服从 Boltzmann 分布,通过求解振动能方程得到振动温度. 一维正激波后热化学松弛过程的计算结果表明,态-态计算预测的温度分布和氧原子浓度分布较好地吻合了文 献中的实验结果,而经典的双温度模型的预测结果误差较大,且不同双温度模型的计算结果比较发散. 态-态方法详细地给出了所有振动能级的变化过程. 无论是正激波还是脱体激波后的流场,都是高振动能级首先得到激发;但是数密度大的低振动能级先达到热平衡,而高能级 分子要经过很长距离后才能达到热平衡. 在驻点附近,复合反应生成的氧气分子处于高振动能级,导致高振动能级分子数密度显著高于平衡分布. 计算还发现,经典双温度模型的离解反应速率明显偏离态-态计算结果,无法准确体现振动离解耦合效应对离解反应 速率的影响,但是 Park 双温度模型将离解失去的振动能取为 0.3$\sim $0.5 倍分子离解能是比较合理的.     

Determinations of both length scale and surface elastic parameters for fcc metals

In the present research, a systematical study of trans-scale mechanics theory is performed. The surface/interface energy density varying with material deformation is considered, and the general surface/interface elastic constitutive equations are derived. New methods to determine the material length scale parameter and the surface elastic parameters based on a simple quasi-continuum method, i.e. the Cauchy–Born rule, are developed and applied to typical fcc metals. In the present research, the material length parameters will be determined through an equivalent condition of the strain energy density calculated by adopting the strain gradient theory and by adopting the Cauchy–Born rule, respectively. Based on the surface constitutive equations obtained in the present research, the surface elastic parameters are calculated by using the Gibbs definition of surface energy density and the Cauchy–Born rule method.