Shock enhancement of cellular structures under impact loading: Part II analysis

Numerical simulations of two distinct testing configurations using a Hopkinson bar (pressure bar behind/ahead of the shock front) are performed with an explicit finite element code. It allows us to confirm the observed test data such as velocity and force time histories at the measurement surface. A comparison of the simulated local strain fields during shock front propagation with those measured by image correlation provides an additional proof of the validity of such simulations.Very simple rate insensitive phenomenological constitutive model are used in such simulations. It shows that the shock effect is captured numerically with a basic densification feature. It means that strength enhancement due to shock should not be integrated in the constitutive model of foam-like materials used in industrial FE codes.In order to separate shock enhancement from entire strength enhancement, an improvement of an existing model with easily identifiable parameters for shock enhancement prediction is proposed. For a quick estimate of the shock enhancement level, a simple power law densification model is proposed instead of the classical RPPL model proposed by Reid and co-workers [Tan et al., 2005. Dynamic compressive strength properties of aluminium foams. Part I—experimental data and observations. J. Mech. Phys. Solids 53, 2174-2205]. It is aimed at eliminating the parameter identification uncertainty of the RPPL model. Such an improved model is easily identifiable and gives a good prediction of the shock enhancement level.     

On the behaviour characterization of metallic cellular materials under impact loading

This paper reviews the common mechanical features of the metallic cellular material under impact loading as well as the characterization methods of such behaviours. The main focus is on the innovations of various testing methods at impact loading rates.Following aspects were discussed in details.（1） The use of soft nylon Hopkinson/Kolsky bar for an enhanced measuring accuracy in order to assess if there is a strength enhancement or not for this class of cellular materials under moderate impact loading;（2） The use of digital image correlations to determine the strain fields during the tests to confirm the existence of a pseudo-shock wave propagation inside the cellular material under high speed impact： （3） The use of new combined shear compression device to determine the loading envelop of cellular materials under impact multiaxial loadings.     

Shock wave interaction with cellular materials

The equations governing the head-on collision of a planar shock wave with a cellular material and a numerical scheme for solving the set of the governing equations were outlined. In addition, the condition for the transmitted compression waves to transform into a shock wave, inside the cellular material was introduced. It was proved analytically that a cellular material cannot be used as a means of reducing the pressure load acting on the end-wall of the shock tube. In subsequent papers, the interaction of planar shock waves with specific cellular materials, e.g., foams and honeycombs will be described in detail.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Dynamic fracture toughness of high strength metals under impact loading: increase or decrease

An elusive phenomenon is observed in previous investigations on dynamic fracture that the dynamic fracture toughness(DFT) of high strength metals always increases with the loading rate on the order of TPa.m 1 /2.s 1.For the purpose of verification,variation of DFT with the loading rate for two high strength steels commonly used in the aviation industry,30CrMnSiA and 40Cr,is studied in this work.Results of the experiments are compared,which were conducted on the modified split Hopkinson pressure bar(SHPB) apparatus,with striker velocities ranging from 9.2 to 24.1 m/s and a constant value of 16.3 m/s for 30CrMnSiA and 40Cr,respectively.It is observed that for 30CrMnSiA,the crack tip loading rate increases with the increase of the striker velocity,while the fracture initiation time and the DFT simultaneously decrease.However,in the tests of 40Cr,there is also an increasing tendency of DFT,similar to other reports.Through an in-depth investigation on the relationship between the dynamic stress intensity factor(DSIF) and the loading rate,it is concluded that the generally increasing tendency in previous studies could be false,which is induced from a limited striker velocity domain and the errors existing in the experimental and numerical processes.To disclose the real dependency of DFT on the loading rate,experiments need to be performed in a comparatively large striker velocity range.     

Effect of Compressible Foam Properties on Pressure Amplification During Shock Wave Impact

A comprehensive study is made of the influence of the physical properties of compressible open-cell foam blocks exposed to shock-wave loading, and particularly on the pressure distribution on the shock tube walls. Seven different foams are used, with three different shock Mach numbers, and three different slab lengths. Foam properties examined include permeability, density, stiffness, tortuosity and cell characteristics. The investigations concentrate on both side-wall and back-wall pressures, and the peak pressures achieved, as well as collapse velocities of the front face and the strength and nature of the reflected shock wave. The consequences of deviations from one-dimensionality are identified; primarily those due to wall friction and side-wall leakage. The results presented are the most comprehensive and wide ranging series conducted in a single facility and are thus a significant resource for comparison with theoretical and numerical studies. The different foams show significant differences in behavior, both in terms of peak pressure and duration, depending primarily on their density and permeability.This paper was based on work presented at the 2nd International Symposium on Interdisciplinary Shock Wave Research, Sendai, Japan on March 1–3, 2005.     

高温后钢管活性粉末混凝土的动态力学性能

采用霍普金森压杆装置对高温后钢管活性粉末混凝土(reactive powder concrete-filled steel tube,RPC-FST)进行冲击压缩实验,分析了应变率效应及温度效应对试件动态力学性能的影响。结果表明:高温(200、300 ℃)后RPC-FST仍具有较好的抗冲击能力、延性和完整性;冲击荷载作用下,RPC-FST的应变率效应明显弱于RPC的应变率效应;随着过火温度的提高,RPC-FST的峰值应力逐渐增大,变形能力增强,抗冲击能力提高。动力提高系数随过火温度的提高而增大,说明高温后RPC-FST的应变率效应更显著。     

A technique for rapid two-stage dynamic tensile loading of polymers

A method for rapid two-stage dynamic-dynamic tensile loading of polymers, based on a tensile Hopkinson bar apparatus, is established. In this technique, the initial incident wave and its reflection are used to load a specimen in quick succession. Consequently, the specimen is stressed, momentarily unloaded, then reloaded until fracture. By adopting appropriate assumptions, a procedure to obtain the associated stress-strain curves for such double-stage loading is formulated. These assumptions are examined experimentally and analytically to substantiate their validity. To verify the proposed approach, a relatively rate-insensitive material, LEXAN 141 polycarbonate, was subjected to two-stage dynamic tension. The stress-strain curves obtained via the procedure established were compared with results from static loading. Favorable correlation between the two indicates that the proposed technique can be applied to the study of load history effects on the dynamic behavior of polymeric materials.     

平板撞击下C30混凝土中冲击波的传播特性

采用1级气炮加载技术和锰铜应力计多点测试技术,开展了C30混凝土在平板撞击条件下的冲击压缩实验研究。基于锰铜应力计实测的应力波形,研究了混凝土中冲击波的传播特性,结果显示冲击波的应力峰值随传播距离呈现明显的衰减特性,衰减过程可分为2个阶段。在早期阶段,卸载波没有赶上前面传播的冲击波,冲击波应力峰值衰减较慢,主要是混凝土材料的本构粘性效应所引起的;而后期阶段应力峰值的快速衰减则归因于混凝土材料的本构粘性效应、后续的来自飞片自由面的反射波追赶卸载、边侧稀疏波卸载及波传播的几何弥散效应的共同作用;另外,冲击波在混凝土中传播的升时也随着传播距离逐渐增大,即由强间断波逐渐转化为弱间断波。     

Numerical study of rate-dependent strength behavior under ramp and shock wave loading

The main objective of the current study was to gain a detailed understanding on the rate-dependent strength behavior under ramp and shock wave loading. A forward, numerical-simulation-based cause and effect analysis was used to address the research objective. The apparent strength associated with shock and ramp wave loadings with different risetimes and shapes was investigated. It was shown that intrinsic material strength could vary with pressure, temperature, and deformation history, but the apparent strength, which was larger than the intrinsic strength, was a result of the interaction between the rate sensitivity of the strength and the rate of the external loading. The degree of interaction led to different levels of mechanical and thermal dissipations and their partition, which was manifested by different temperature, stress, and deformation histories.     

Search for conditions of compressive fracture of hard brittle ceramics at impact loading

In this paper we discuss three different experimental configurations to diagnosing the modes of inelastic deformation and to evaluating the failure thresholds at shock compression of hard brittle solids. One of the manifestations of brittle material response is the failure wave phenomenon, which has been previously observed in shock-compressed glasses. However, based on the measurements from our “theory critical” experiments, both alumina and boron carbide did not exhibit this phenomenon. In experiments with free and pre-stressed ceramics, while the Hugoniot elastic limit (HEL) in high-density B₄C ceramic was found to be very sensitive to the transverse stress, it was found relatively less sensitive in Al₂O₃, implying brittle response of the boron carbide and ductile behavior of alumina. To further investigate the effects of stress states on the shock response of brittle materials, a “divergent flow or spherical shock wave” based plate impact experimental technique was employed to vary the ratio of longitudinal and transversal stresses and to probe conditions for compressive fracture thresholds. Two different experimental approaches were considered to generate both longitudinal and shear waves in the target through the impact of convex flyer plates. In the ceramic target plates, the shear wave separates a region of highly divergent flow behind the decaying spherical longitudinal shock wave and a region of low-divergent flow. Experiments with divergent shock loading of alumina and boron carbide ceramic plates coupled with computer simulations demonstrated the validity of these experimental approaches to develop a better understanding of fracture phenomena.     

On the local nature of the strain field calculation method for measuring heterogeneous deformation of cellular materials

A strain field calculation method based on the optimal local deformation gradient technique has been developed to calculate the ‘local’ strain tensor of cellular materials using cell-based finite element models. The local nature and accuracy of this method may be strongly dependent on the cut-off radius, which is introduced to collect the effective nodes for determining the optimal local deformation gradient of a node. Two different schemes are first analyzed to determine the suitable cut-off radius by characterizing the heterogeneous deformation of Voronoi honeycombs under uniaxial compression and we suggest that in Scheme 1, the cut-off radius defined based on the reference configuration is about 1.5 times the average cell radius; in Scheme 2, the cut-off radius defined based on the current configuration is about 0.5 times the average cell radius. Then, Scheme 3, a combined scheme of the two former schemes, is further suggested. It is demonstrated that the optimal cut-off radius in Scheme 3 characterizes the local strain reasonable well whether the compression rate is low or high. Finally, the strain field calculation method with the optimal cut-off radius is applied to reveal the evolution of the heterogeneous deformation of two different configurations of double-layer cellular cladding under a linear decaying blast load. The 2D fields and the 1D distributions of local engineering strain are calculated. These results interpret the shock wave propagation mechanisms in both claddings and provide useful understanding in the design of a double-layer cellular cladding.     

Shock pattern in the plume of rocket nozzles: needs for design consideration

For ideal nozzles, basically two different types of shock structures in the plume may appear for overexpanded flow conditions, a regular shock reflection or a Mach reflection at the nozzle centreline. Especially for rocket propulsion, other nozzle types besides the ideal nozzles are often used, including simple conical, thrust-optimized or parabolic contoured nozzles. Depending on the contour type, another shock structure may appear: the so-called cap-shock pattern. The exact knowledge of the plume pattern is of importance for mastering the engine operation featuring uncontrolled flow separation inside the nozzle, appearing during engine start-up and shut-down operation. As consequence of uncontrolled flow separation, lateral loads may be induced. The side-load character strongly depends on the nozzle design, and is a key feature for the nozzle’s mechanical structure layout. It is shown especially for the VULCAIN and VULCAIN 2 nozzle, how specific shock patterns evolve during transients, and how - by the nozzle design - undesired flow phenomena can be avoided.     

Behavior of the surface of a bubbly liquid after detonation wave impact

In safety engineering, one position of interest inside heterogeneous systems of the type liquid–gas is the contact surface between these two phases. Under certain conditions, e.g. shock wave impact, phenomena can take place at this position that can have a significant influence on the explosion behavior of the system. In this work an investigation is presented about the existence of such phenomena on the surface of liquid cyclohexane with or without the existence of oxygen containing bubbles. The observations have been performed during the time before, as well as after, a detonation wave reflection on that surface. High-speed pressure and optical measurements have been applied. Apart from the experimental observations, also a theoretical analysis and discussion is presented in this contribution, which contains the comparison between calculated and experimental values.     

A phenomenological model of damage in articular cartilage under impact loading

Damage progression in high-strain rate and impact tests on articular cartilage is considered. A new type of kinetic damage evolution law is proposed and used to draw implications about the accumulated damage and the coefficient of restitution. Based on the developed damage model, a new fracture criterion is introduced.     

Impact behavior of honeycombs under combined shear-compression. Part I: Experiments

This paper presents a combined shear-compression impact test for soft cellular materials designed in order to investigate their behavior under impact multiaxial loadings. A large-diameter Nylon Split Hopkinson Pressure Bar system (SHPB) with beveled ends of different angles is used to apply the desired shear-compression combinations. The data processing methods are studied and validated by virtual testing data generated with FEM simulations. A series of experiments on an aluminum honeycomb were performed at the impact velocity of about 15 m/s with the loading angles ranging from 0° (corresponding to the pure compression) to 60°. It shows a strong effect of the additional shear loading because both the initial peak and the crush strength decrease with increasing loading angles. The quasi-static shear-compression experiments were also performed using the same beveled ends on a universal INSTRON machine and a notable strength enhancement under impact loading is observed. Images captured during quasi-static and impact tests permit for the determination of the two co-existing deforming patterns under combined shear-compression and reveal the influence of the loading rate on the occurrence of these two patterns.     

冲击载荷作用下黑砂岩动态断裂参数的分形修正

基于分形理论研究了偏折裂纹扩展路径对动载荷作用下黑砂岩的动态断裂力学参数的测试误差影响作用,采用传统的分离式霍普金森压杆（split Hopkinson pressure bar, SHPB）实验装置对修正侧开单裂纹半孔板（improved single cleavage semi-circle specimen, ISCSC）试样进行动态冲击实验,随后采用裂纹扩展计进行裂纹起裂时间与裂纹扩展速度等动态断裂力学参数测试,采用分形理论对测试的裂纹扩展速度与动态应力强度因子进行修正,利用实验-数值法对黑砂岩的动态断裂韧度进行计算。研究结果表明,ISCSC构型构件能够有效应用于岩石材料动态裂纹扩展行为的研究,并发生了止裂现象,经分形修正的裂纹扩展速度与动态断裂韧度更接近实际裂纹动态扩展情况,修正前后得到黑砂岩材料的裂纹扩展速度误差为33.51%,动态断裂韧度最大误差为7.68%,说明利用分形理论对动态断裂韧度等动态断裂参数计算更合理。

     

纵向冲击压缩下LY12铝合金圆环的塑性失稳

通过对铝合金圆环的纵向冲击压缩研究发现,一定条件下在试件的宏观塑性硬化阶段会出现明显的应力降过程。为揭示此应力降的发生机制,对润滑、细磨、粗磨3种端面粗糙条件下,外径、内径和高度比值为6:3:2的LY12铝合金圆环进行系统的Hopkinson压杆纵向冲击实验。结果表明:应力降主要发生在较大的应变和较高的应变率条件。进一步对实验样品的金相观察发现:应力降产生的内在机制为绝热剪切带的形成和发展,此现象是一种动态塑性失稳的过程。以上结果为金属材料在冲击条件下绝热剪切带产生的研究提供了参考。     

混凝土类材料动态压缩强度在多维应力状态下的应变率效应

对混凝土类材料动态压缩应变率效应研究的发展及问题进行了概述,对比不同应力状态下混凝土类材料动态压缩应变率效应的表现特征,揭示了不同加载路径下实测动态强度提高系数的显著差异。研究表明,在高应变率下,基于初始一维应力加载路径的试件将因横向惯性效应导致的侧向围压而演化至多维应力状态,传统霍普金森杆技术无法获得高应变率下基于真实一维应力路径的动态强度提高系数,在强度模型中直接应用实测数据将过高估计材料的动态强度。鉴于应变率效应的加载路径依赖性,将仅包含应变率的强度提高系数模型扩展至同时计及应变率和应力状态的多维应力状态模型,并结合Drucker-Prager准则在强度模型中给予了实现。针对具有自由和约束边界试件开展的数值霍普金森杆实验表明,多维应力状态下的应变率效应模型可以考虑应变率效应随应力状态改变的特点,从而准确预测该类材料的动态压缩强度。研究结果可为正确应用霍普金森杆技术确定脆性材料的动态压缩强度提供参考。

     

On axisymmetric motion of an incompressible elastic medium under impact loading

The method of matched asymptotic expansions was employed to obtain approximate solutions to the one-dimensional boundary-value problems of nonlinear dynamic elasticity theory of impact loading on the surface of a cylindrical cavity of an incompressible medium that causes antiplane motion or torsion of the medium. The expansion of the solution in the near-front region is based on solutions of evolution equations different from the equations for quasi-simple waves. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 6, pp. 144–151, November–December, 2006.