泡沫铝夹层结构抗冲击性能的近场动力学模拟分析

针对某光学舱所采用的泡沫铝夹层防护结构在破片冲击下的抗冲击性能问题,采用Monte-Carlo方法创建了泡沫铝结构的二维细观模型,在常规态型近场动力学理论中引入了Mises屈服准则和线性各向同性强化模型,建立了近场动力学塑性本构的数值计算框架。基于近场动力学计算程序模拟了低速冲击作用下泡沫铝夹层结构的塑性变形以及有机玻璃背板的裂纹扩展形态,分析了泡沫铝芯材孔隙率对该夹层结构抗冲击性能和损伤模式的影响规律。结果表明：泡沫铝夹层结构良好的塑性变形能力是其发挥缓冲与防护作用的主要因素,并且在一定范围内,泡沫铝芯材孔隙率越高,则夹层结构具有更好的抗冲击性能;当泡沫铝孔隙率从0.4提升到0.7时,泡沫铝对冲击物的动能吸收率从90%提高到99%;模拟结果与实验结果具有较好的一致性,验证了模拟结果的准确性和分析结论的有效性。通过数值模拟,预测了有机玻璃背板的裂纹扩展形态,发现提高泡沫铝的孔隙率能获得更好的防护效果。     

Nonlinear deformation of elastomeric foams

The nonlinear deformation of a porous foam-type elastomeric material is studied, both theoretically and experimentally. The elastomer is modeled by the neo-Hookean material. The one-dimensional compressive behavior of the foam is analysed by using certain kinematic assumptions. The stress required to compress the foam is predicted by the model in terms of the porosity of the foam and the single constant in the neo-Hookean stress-strain form. A particular silicone foam is used as a test of the theory. The neo-Hookean constant is evaluated from a test of the homogeneous elastomer. Hence the behavior of the corresponding foam is predicted theoretically and compared with experimental results. The general results are applicable to closed-cells foams of intermediate density.     

泡沫金属夹芯梁在重复冲击下的动态响应

为了研究重复冲击载荷作用下泡沫金属夹芯梁的动态响应,采用Abaqus数值仿真软件,基于可压碎泡沫模型（crushable foam）,建立了泡沫金属夹芯梁遭受楔形质量块冲击的有限元模型。通过将仿真获得的夹芯梁上下面板最终挠度与重复冲击实验结果进行对比,验证仿真方法的准确性。在此基础之上,分析了泡沫金属夹芯梁在楔形质量块重复冲击作用下的变形模式、加卸载过程以及能量耗散特性。结果表明,在重复冲击载荷作用下,夹芯梁的变形不断累积,上面板主要出现局部凹陷和整体弯曲,而芯层则是局部压缩,下面板表现为整体弯曲。在重复加卸载过程中,加卸载刚度随着冲击次数的增加而增大。随着冲击次数的增加,上面板和芯层的能量吸收增量不断减小,而下面板的能量吸收增量不断增加,且最终均趋于稳定。泡沫金属夹芯梁的塑性变形能增量不断减小,而回弹系数随着冲击次数逐渐增加,最后趋于稳定值。     

超高速碰撞下相变效应的交错网格物质点法研究

超高速碰撞下的结构毁伤过程常伴随材料相变、断裂失效及碎片云的产生与演化,其所具有的压力强间断、材料非线性等几何大变形等问题给数值模拟带来了困难。交错网格物质点法SGMP通过背景网格格心积分消除了物质点法MPM跨网格误差,是模拟固体冲击爆炸等极端大变形与材料破坏问题的一种有效数值分析方法。本文将含金属相变的GRAY状态方程及含非线性内聚力断裂的Johnson-Cook修正金属模型引入SGMP中,模拟超高速碰撞单层、多层靶板问题。结果表明,SGMP和修正金属模型可以很好地模拟超高速碰撞问题中的碎片云形貌特征和相变效应。     

Deformation rate effects on failure modes of open-cell Al foams and textile cellular materials

The compressive behavior of open-cell aluminum alloy foam and stainless steel woven textile core materials have been investigated at three different deformation rate regimes. Quasi-static compressive tests were performed using a miniature loading frame, intermediate rates were achieved using a stored energy Kolsky bar, and high strain rate tests were performed using a light gas gun.In agreement with previous studies on foam materials, the strain rate was not found to have a significant effect on the plateau stress of metallic foams. For all the tests, real time imaging of the specimen combined with digital image correlation analysis allowed the determination of local deformation fields and failure modes. For the Kolsky bar tests, the deformations in the foam specimen were found to be more distributed than for the quasi-static test, which is attributed to moderate inertia effects. The differences in failure mode are more dramatic for the gas gun experiments, where a full compaction shock wave is generated at the impact surface. The stresses in front and behind the shock wave front were determined by means of direct and reverse gas gun impact tests, i.e., stationary and launched specimen, respectively. A one-dimensional shock wave model based on an idealized foam behavior is employed to gain insight into the stress history measurements. We show that the predictions of the model are consistent with the experimental observations. Woven textile materials exhibited moderate dependence of strength on the deformation rate in comparison with open-cell foam materials.     

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.     

近似不可压软材料动力分析的完全拉格朗日物质点法

物质点法(MPM)在模拟非线性动力问题时具有很好的效果,其已被广泛应用于许多大变形动力问题的分析中.然而传统的MPM在模拟不可压或近似不可压材料的动力学行为时会产生体积自锁,极大地影响模拟精度和收敛性.本文针对近似不可压软材料的大变形动力学行为,提出一种混合格式的显式完全拉格朗日物质点法(TLMPM).首先基于近似不可压软材料的体积部分应变能密度,引入关于静水压力的方程;之后将该方程与动量方程基于显式物质点法框架进行离散,并采用完全拉格朗日格式消除物质点跨网格产生的误差,提升大变形问题的模拟精度;对位移和压强场采用不同阶次的B样条插值函数并通过引入针对体积变形的重映射技术改进了算法,提升算法的准确性.此外,算法通过实施一种交错求解格式在每个时间步对位移场和压强场依次进行求解.最后,给出几个典型数值算例来验证本文所提出的混合格式TLMPM的有效性和准确性,计算结果表明该方法可以有效处理体积自锁,准确地模拟近似不可压软材料的大变形动力学行为.     

Full-field characterization of mechanical behavior of polyurethane foams

The foam material of interest in this investigation is a rigid closed-cell polyurethane foam PMDI with a nominal density of 20 pcf (320 kg/m³). Three separate types of compression experiments were conducted on foam specimens. The heterogeneous deformation of foam specimens and strain concentration at the foam–steel interface were obtained using the 3-dimensional digital image correlation (3D-DIC) technique. These experiments demonstrated that the 3D-DIC technique is able to obtain accurate and full-field large deformation of foam specimens, including strain concentrations. The experiments also showed the effects of loading configurations on deformation and strain concentration in foam specimens. These DIC results provided experimental data to validate the previously developed viscoplastic foam model (VFM). In the first experiment, cubic foam specimens were compressed uniaxially up to 60%. The full-field surface displacement and strain distributions obtained using the 3D-DIC technique provided detailed information about the inhomogeneous deformation over the area of interest during compression. In the second experiment, compression tests were conducted for cubic foam specimens with a steel cylinder inclusion, which imitate the deformation of foam components in a package under crush conditions. The strain concentration at the interface between the steel cylinder and the foam specimen was studied in detail. In the third experiment, the foam specimens were loaded by a steel cylinder passing through the center of the specimens rather than from its end surface, which created a loading condition of the foam components similar to a package that has been dropped. To study the effects of confinement, the strain concentration and displacement distribution over the defined sections were compared for cases with and without a confinement fixture.     

Large amplitude oscillatory shear of pseudoplastic and elastoviscoplastic materials

We explore the utility of strain-controlled large amplitude oscillatory shear (LAOS) deformation for identifying and characterizing apparent yield stress responses in elastoviscoplastic materials. Our approach emphasizes the visual representation of the LAOS stress response within the framework of Lissajous curves with strain, strain rate, and stress as the coordinate axes, in conjunction with quantitative analysis of the corresponding limit cycle behavior. This approach enables us to explore how the material properties characterizing the yielding response depend on both strain amplitude and frequency of deformation. Canonical constitutive models (including the purely viscous Carreau model and the elastic Bingham model) are used to illustrate the characteristic features of pseudoplastic and elastoplastic material responses under large amplitude oscillatory shear. A new parameter, the perfect plastic dissipation ratio, is introduced for uniquely identifying plastic behavior. Experimental results are presented for two complex fluids, a pseudoplastic shear-thinning xanthan gum solution and an elastoviscoplastic invert-emulsion drilling fluid. The LAOS test protocols and the associated material measures provide a rheological fingerprint of the yielding behavior of a complex fluid that can be compactly represented within the domain of a Pipkin diagram defined by the amplitude and timescale of deformation.     

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

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

高温下轻质泡沫铝动态力学性能实验

对传统的分离式Hopkinson压杆装置加以改进,设计了一种长杆直接撞击Hopkinson杆的实验方案,检测出低波阻抗材料在高温动态加载下的应力均匀性。对轻质泡沫铝材料的实验表明,在同一撞击速度下,温度越高,试件两端的应力均匀性越差,增加温度与提高撞击速度均会导致泡沫铝材料冲击端与支撑端的应力不均匀性。根据高温下应力均匀性的实验结果,确定高温下试件均匀变形对应的冲击速度,再通过传统的分离式Hopkinson压杆实验得出泡沫铝在高温动态下的力学性能。     

煤岩巷道冲击破坏过程相似模拟试验研究

采用相似材料模拟试验研究方法,利用单轴加载系统和数字散斑光测技术,根据相似材料配比制备脆性材料,对标准试块和三维巷道模型进行了相似试验研究。对相似模拟试验结果分析可知,实验室制备脆性材料模拟冲击地压巷道破坏过程是可行的,冲击破坏特征较明显;散斑计算结果说明,试块位移和应变变化与裂缝出现、扩展和贯通具有一致性,且竖向裂缝开展数量多、速度快;相似模拟试验结果也表明巷道顶板受冲击作用和下沉位移较大,两边帮向巷道内部弯曲变形也较大,甚至出现煤岩抛出和巷道垮塌。本试验也为进一步揭示煤岩巷道冲击破坏规律提供了新的途径。     

软质聚氨酯泡沫的动态压缩力学性能和本构模型

采用Instron 9350落锤试验机研究了中低应变率下软质聚氨酯泡沫的动态压缩力学性能,分析了其应力-应变响应特征和应变率敏感性,讨论了应变率对材料应变率敏感性指数和能量吸收特性的影响,并基于实验结果建立了可准确描述其压缩力学响应的率相关本构模型。结果表明,软质聚氨酯泡沫的静动态压缩应力-应变响应具有典型的三阶段特征,且呈现出明显的应变率强化效应。准静态加载下,材料具有较高的吸能效率但能量吸收值较小,应变率对最大吸能效率和比吸能的影响较小;动态加载下,随着应变率的增加,最大吸能效率显著减小而比吸能明显增大。考虑应变率影响的修正Sherwood-Frost模型和修正Avalle模型都能够很好地表征软质聚氨酯泡沫的静动态压缩应力-应变响应,但修正Avalle模型的参数较少,更便于工程应用。研究结果可为软质聚氨酯泡沫抗冲击结构的设计和优化提供指导。     

土体颗粒物流动物质点法模拟的弹塑性和非牛顿流体本构模型比较研究

土体颗粒物流动是一种典型的大变形破坏,具有非牛顿流体的流动特征。准确模拟土体颗粒物的流动及冲击过程,对滑坡和泥石流等地质灾害的防治具有重要意义。物质点法是一种无网格粒子类方法,已在各类大变形问题中得到了广泛应用。以往土体颗粒物流动的模拟,通常采用弹塑性本构模型,但缺乏对非牛顿本构模型的模拟分析。本文引入非牛顿本构模型的模拟分析,旨在为土体颗粒物流动模拟提供一种新的方法与思路。非牛顿本构模型的模拟分析是将非牛顿广义Cross模型引入三维物质点法,通过人工阻尼力模拟颗粒间的摩擦力,对土体颗粒物的坍塌、沿斜面滑动以及冲击障碍物等问题进行了动态模拟,研究了其运动全过程,并与弹塑性本构模型的模拟结果进行了对比验证。结果表明,基于非牛顿流体本构模型的物质点法可以较好地模拟土体颗粒物加速、减速到再次稳定的流动全过程及其对障碍物的冲击效应。     

Phenomenological modeling of viscous electrostrictive polymers

A common usage for electroactive polymers (EAPs) is in different types of actuators, where advantage is taken of the deformation of the polymer due to an electric field. It turns out that time-dependent effects are present in these applications. One of these effects is the viscoelastic behavior of the polymer material. In view of the modeling and simulation of applications for EAP within a continuum mechanics setting, a phenomenological framework for an electro-viscoelastic material model is elaborated in this work. The different specific models are fitted to experimental data available in the literature. While the experimental data used for inherent electrostriction is restricted to small strains, a large strain setting is used for the model in order to account for possible applications where the polymers undergo large deformations, such as in pre-strained actuators.     

基于数字图像相关方法对冲击载荷下泡沫铝全场变形过程的测试

应用分离式霍普金森压杆(以下简称SHPB)和高速摄影装置研究了冲击载荷下泡沫铝试件全场变形的测量方法。使用SHPB对泡沫铝试件进行冲击压缩实验,同时用高速摄影装置对实验过程进行全程跟踪拍摄。将得到的高速摄影图像采用数字图像相关方法进行分析,由此可得到冲击压缩过程中泡沫铝试件全场应变的分布和变化规律。此研究揭示了冲击载荷下泡沫铝试件局部化变形的发展过程,为研究泡沫铝在不同冲击载荷下不同变形模式的内在机制提供了新的可靠的方法。     

基于力热耦合材料模型的Steven试验数值模拟方法

在力热耦合材料模型中,增加炸药自热放能模型,建立了Steven试验的力-热-化耦合的数值模拟方法。数值计算模型中,应力应变关系采用双线性硬化弹塑性模型,炸药受力后的热作用采用各向同性热材料模型,炸药的化学反应采用Arrhenius反应率函数,同时还考虑了升温和熔化对材料力学、热学性能的影响。针对标准Steven试验,通过数值分析得到了靶板的变形情况和炸药点火的速度阈值,将计算结果与实验数据进行了比较,两者符合较好。表明该方法可以较好地模拟Steven试验,而且与以往的分析模型和方法相比,本文的方法不需要增加经验性的点火准则和判据,具有更广泛的适用性,可以为研究低速撞击条件下炸药的力热响应和局域化点火问题提供参考。     

泡沫铝泡孔动态变形特性研究

在SHPB实验中采用应变冻结法将试件冻结于设定的压缩应变,然后观察内部泡孔的变形情况并讨论其变形机理。分别对同一相对密度(40%)、不同基体材料的两种泡沫铝以及两种相对密度(19%和8%)、相同基体材料的泡沫铝进行实验,讨论了基体材料的本构关系、泡沫的相对密度对变形模式的影响。利用ANSYS/LS DYNA软件对泡沫铝在动态压缩时泡孔的变形进行数值模拟计算,进一步说明泡沫铝泡孔变形的模式与基体材料本构关系及相对密度的关系。     

ZrCuNiAlAg块体非晶合金JH-2模型研究

为了更好的进行ZrCuNiAlAg块体非晶合金药型罩的爆炸成形及侵彻仿真研究,首要就是建立其材料模型。本文结合ZrCuNiAlAg块体非晶合金力学性能试验结果计算得到了材料的JH-2模型参数,研究确定了ZrCuNiAlAg块体非晶合金JH-2模型。为了验证ZrCuNiAlAg块体非晶合金JH-2模型的准确性,采用Autodyn建立了平板撞击试验有限元模型,模拟了ZrCuNiAlAg块体非晶合金材料在高压、高应变率等环境条件下的变形过程,仿真计算得到的靶板背面自由面粒子速度与试验结果相比,速度平均偏差均在3%以内,表明ZrCuNiAlAg块体非晶合金JH-2模型能很好的描述该材料在大变形、高应变率、高压等环境条件下的力学行为,验证了ZrCuNiAlAg块体非晶合金JH-2模型准确性。