Dynamic limit analysis formulation for impact simulation of structural members

This paper introduces an extended concept of limit analysis to deal with the dynamic equilibrium condition considering the inertia and strain-rate effect for dynamic behavior of structures. The conventional limit analysis method has been applied to only static collapse analysis of structures without consideration of dynamic effects in the structural behavior. A dynamic formulation for the limit analysis has been derived for incremental analysis dealing with time integration, strain and stress evaluation, strain hardening, strain-rate hardening and thermal softening. The time dependent term in the governing equation is integrated with the WBZ-α method. The dynamic material behavior is described by the Johnson–Cook model in order to consider strain-rate hardening and thermal softening as well as strain hardening. Simulations have been carried out for impact analysis of a Taylor bar and an S-rail and their numerical results are compared with elasto-plastic explicit analysis results by LS-DYNA3D. Comparison demonstrates that the dynamic finite element limit analysis can predict the crashworthiness of structural members effectively with less effort and computing time than the commercial code compared. The crashworthiness of a structure with the rate-dependent constitutive model is also compared to that with the quasi-static constitutive relation in order to investigate the dynamic effect on deformation of structures.     

Strain-rate effect on initial crush stress of irregular honeycomb under dynamic loading and its deformation mechanism

The seemingly contradictory understandings of the initial crush stress of cellular materials under dynamic loadings exist in the literature, and a comprehensive analysis of this issue is carried out with using direct information of local stress and strain. Local stress/strain calculation methods are applied to determine the initial crush stresses and the strain rates at initial crush from a cell-based finite element model of irregular honeycomb under dynamic loadings. The initial crush stress under constant-velocity compression is identical to the quasi-static one, but less than the one under direct impact, i.e. the initial crush stresses under different dynamic loadings could be very different even though there is no strain-rate effect of matrix material. A power-law relation between the initial crush stress and the strain rate is explored to describe the strain-rate effect on the initial crush stress of irregular honeycomb when the local strain rate exceeds a critical value, below which there is no strain-rate effect of irregular honeycomb. Deformation mechanisms of the initial crush behavior under dynamic loadings are also explored. The deformation modes of the initial crush region in the front of plastic compaction wave are different under different dynamic loadings.     

随机缺陷对蜂窝结构动态行为影响的有限元分析

通过对胞壁随机移除的蜂窝结构动态变形过程的有限元模拟,分析了随机缺陷对蜂窝结构变形模式的影响,得到蜂窝结构在两个加载方向上的变形模式图及不同模式间转换的临界速度. 对含缺陷蜂窝结构平台应力的研究发现,当变形模式为过渡模式或动态模式时结构平台应力与冲击速度的平方成线性关系. 相同密度下,低缺陷蜂窝结构的平台应力在由过渡模式向动态模式转变的临界速度附近高于规则蜂窝结构,较高的随机缺陷则使蜂窝结构的平台应力在由准静态模式向过渡模式转变的临界速度附近显著下降.关键词:多孔材料,蜂窝,缺陷,平台应力,有限元分析      

沉痛悼念许兵先生

通过对胞壁随机移除的蜂窝结构动态变形过程的有限元模拟,分析了随机缺陷对蜂窝 结构变形模式的影响,得到蜂窝结构在两个加载方向上的变形模式图及不同模式间转换的临 界速度. 对含缺陷蜂窝结构平台应力的研究发现,当变形模式为过渡模式或动态模式时结构 平台应力与冲击速度的平方成线性关系. 相同密度下,低缺陷蜂窝结构的平台应力在由过渡 模式向动态模式转变的临界速度附近高于规则蜂窝结构,较高的随机缺陷则使蜂窝结构的平 台应力在由准静态模式向过渡模式转变的临界速度附近显著下降. 关键词：多孔材料,蜂窝,缺陷,平台应力,有限元分析     

In-plane crushing behavior and energy absorption design of composite honeycombs

Theoretical analysis and numerical simulation methods were used to study the in-plane crushing behavior of single-cell structures and regular and composite honeycombs. Square, hexagonal, and circular honeycombs were selected as honeycomb layers to establish composite honeycomb models in the form of composite structures and realize the complementary advantages of honeycombs with type I and type II structures. The effects of honeycomb layer arrangement, plastic collapse strength, relative density, and crushing velocity on the deformation mode, plateau stress, load uniformity, and energy absorption performance of the composite honeycombs were mainly considered. A semi-empirical formula for plateau stress and energy absorption rate per unit mass for the composite honeycombs was developed. The results showed that the arrangement mode of honeycomb layers is an important factor that affects their mechanical properties. Appropriately selecting the arrangement of honeycomb layers and the proportion of honeycomb layers with different structures in a composite honeycomb can effectively improve its load uniformity and control the magnitude of plateau stress and energy absorption capacity.     

Strain-rate history and temperature effects on the torsional-shear behavior of a mild steel

Previous investigations on the effects of strain-rate and temperature histories on the mechanical behavior of steel are briefly reviewed. A study is presented on the influence of strain rate and strain-rate history on the shear behavior of a mild steel, over a wide range of temperature Experiments were performed on thin-walled tubular specimens of short gage length, using a torsional split-Hopkinson-bar apparatus adapted to permit quasi-static as well as dynamic straining at different temperatures. The constant-rate behavior was first measured at nominal strain rates of 10^?3 and 10³ s^?1 for ?150, ?100, ?50, 20, 200 and 400°C. Tests were then carried out, at the same temperatures, in which the strain rate was suddenly increased during deformation from the lower to the higher rate at various large values of plastic strain. The increase in rate occurred in a time of the order of 20 μs so that relatively little change of strain took place during the jump. The low strain-rate results show a well-defined elastic limit but no yield drop, a small yield plateau is found at room temperature. The subsequent strain hardening shows a maximum at 200°C, when serrated flow occurs and the ductility is reduced. The high strain-rate results show a considerable drop of stress at yield. The post-yield flow stress decreases steadily with increasing temperature, throughout the temperature range investigated. At room temperature and below, the strain-hardening rate becomes negative at large strains. The adiabatic temperature rise in the dynamic tests was computed on the assumption that the plastic work is entirely converted to heat. This enabled the isothermal dynamic stress-strain curves to be calculated, and showed that considerable thermal softening took place. The initial response to a strain-rate jump is approximately elastic, and has a magnitude which increases with decrease of testing temperature; it is little affected by the amount of prestrain. At 200 and 400° C, a yield drop occurs after the initial stress increment. The post-jump flow stress is always greater than that for the same strain in a constant-rate dynamic test, the strain-hardening rate becoming negative at large strains or low testing temperature. This observed effect of strain-rate history cannot be explained by the thermal softening accompanying dynamic deformation. These and other results concerning total ductility under various strain-rate and temperature conditions show that strain-rate history strongly affects the mechanical behavior of the mild steel tested and, hence, should be taken into account in the formulation of constitutive equations for that material.     

Dynamic crushing behavior of honeycomb structures with irregular cell shapes and non-uniform cell wall thickness

Dynamic crushing responses of honeycomb structures having irregular cell shapes and non-uniform cell wall thickness are studied using the Voronoi tessellation technique and the finite element (FE) method. FE models are constructed for such honeycomb structures based on Voronoi diagrams with different degrees of cell shape irregularity and cell wall thickness non-uniformity. The plateau stress, the densification strain energy and the initiation strain are determined using the FE models. Simulation results reveal that the “X” and “V” shaped deformation modes evident in a perfectly ordered honeycomb at low or moderate impact velocities are disrupted as cell shapes become irregular and/or cell wall thickness gets non-uniform. The “I” shaped deformation mode is clearly seen in all honeycomb structures at high impact velocities. Both the plateau stress and the densification strain energy are found to decrease as the degree of cell shape irregularity or the degree of cell wall thickness non-uniformity increases, with the weakening effect induced by the presence of non-uniform cell wall thickness being more significant. When the two types of imperfections co-exist in a honeycomb structure, the interaction between them is seen to exhibit a complicated pattern and to have a nonlinear effect on both the plateau stress and the densification strain energy. It is also found that stress waves propagate faster in a honeycomb structure having irregular cell shapes and slower in a honeycomb structure having non-uniform cell wall thickness than in a perfectly ordered honeycomb. Finally, the strain hardening of the cell wall material is seen to have a strengthening effect on the plateau stress, which is more significant for perfectly ordered honeycombs than for imperfect honeycomb structures.     

含随机填充孔圆形蜂窝结构的面内冲击性能

研究多孔材料细观结构与宏观力学性能之间的关系, 建立具有固定相对密度的含随机固体填充孔的圆形蜂窝结构模型。在此模型的基础上具体讨论了不同孔洞填充比和冲击速度对圆形蜂窝结构变形模式、动态冲击平台应力以及能量吸收性能的影响。研究结果表明:填充孔在蜂窝变形过程中有局部牵制作用, 蜂窝材料变形模式仍为准静态模式、过渡模式、动态模式; 当变形模式为过渡模式或动态模式时, 结构的平台应力与速度的平方成线性关系, 存在明显的速度效应; 高速冲击下, 含固体填充孔的蜂窝结构单位质量吸收的能量高于规则蜂窝结构。研究结果可为蜂窝材料的研究和设计提供参考。     

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

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

铝蜂窝的动态力学性能及影响因素

采用显式动力有限元方法研究了具有不同胞元结构的六角形铝蜂窝在冲击荷载下的力学性能,讨 论了铝蜂窝的变形模式和动态承载力以及影响因素。通过改变胞壁夹角得到5种不同的胞元结构,计算采用 了3种冲击速度。结果表明,在准静态变形模式下,胞元的几何因素对铝蜂窝的承载力起主导作用;一旦蜂窝 的变形呈现动态模式后,惯性效应显著,对蜂窝承载力起决定作用,胞元几何因素的影响不再明显;在过渡模 式下,惯性效应与几何因素共同主导蜂窝的动态承载力,并且冲击速度越高,惯性效应的影响越大。     

Dynamic Crushing Strength Analysis of Auxetic Honeycombs

The in-plane dynamic crushing behavior of re-entrant honeycomb is analyzed and compared with the conventional hexagon topology.Detailed deformation modes along two orthogonal directions are examined,where a parametric study of the effect of impact velocity and cell wall aspect ratio is performed.An analytical formula of the dynamic crushing strength is then deduced based on the periodic collapse mechanism of cell structures.Comparisons with the finite element results validate the effectiveness of the proposed analytical method.Numerical results also reveal higher plateau stress of re-entrant honeycomb over conventional hexagon topology,implying better energy absorption properties.The underlying physical understanding of the results is emphasized,where the auxetic effect(negative Poisson's ratio) induced in the re-entrant topology is believed to be responsible for this superior impact resistance.     

受轴向冲击薄壁圆管的几何畸变相似律研究

对于受轴向冲击载荷作用的薄壁圆管动态响应的相似律问题,由于圆管的薄壁特性导致厚度无法与高度和半径按相同的比例进行结构缩放,从而产生模型的几何畸变,此时传统的相似律已无法描述原型与畸变模型之间的动态响应规律。基于薄壁圆管轴向冲击问题的控制方程,通过能量守恒和量纲分析,推导了考虑几何畸变条件下轴向冲击载荷作用的理想弹塑性薄壁圆管动态响应的相似律。通过在给定应变与应变率区间上建立比例模型预测的流动屈服应力与原型流动屈服应力的最佳逼近关系,将几何畸变相似律进一步推广至包含应变率和应变硬化的材料。通过数值方法验证了提出的几何畸变模型相似律的适用性。分析结果表明,提出的考虑厚度畸变的受轴向冲击薄壁圆管的相似律可用于预测原型结构的冲击动态响应,并显著降低比例模型与原型结构平均载荷和能量的偏差。     

基于十四面体模型的闭孔泡沫材料动态力学性能的有限元分析

应用有限元方法分析了基于十四面体模型的三维闭孔泡沫材料的动态力学性能。计算中所有十四面体具有相同的尺寸,主要研究了不同初始冲击速度、不同相对密度以及组成泡沫的机体材料的应变强化对泡沫材料的变形模态、平台力及密实化应变能的影响,尽可能全面地描述了泡沫材料的能量吸收特性。数值结果表明:冲击速度对泡沫模型的模态影响较大,特别受到高速冲击时,冲击端泡沫形成I形然后向支撑端传播;相对密度对能量吸收能力的贡献较大,密实化应变能随相对密度呈二次曲线变化;冲击速度、相对密度及机体材料的应变强化分别与坪应力呈线性关系。     

High-strain-rate compressive behavior of a rigid polyurethane foam with various densities

The dynamic compressive stress-strain behavior of a rigid polyurethane foam with four values of density (78, 154, 299, and 445 kg/m³) has been determined in the strain-rate range of 1000–5000 s⁻¹. A pulse shaping technique was used with a split Hopkinson pressure bar to ensure homogeneous deformation in the foam specimens under dynamic compression. Dynamic stress equilibrium in the specimen was monitored during each experiment using piezoelectric force transducers mounted close to the specimen end-faces. Quasi-static experiments were also performed to demonstrate rate effects. Experimental results show that both the quasistatic and the dynamic stress-strain curves of the foam exhibit linear elasticity at small strains until a peak is reached. After the peak, the stress-strain curves have a plateau region followed by a densification region. The peak stress is strain-rate sensitive and depends on the square of the foam density.     

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

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

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.     

Why does necking ignore notches in dynamic tension?

Recent experimental work has revealed that notched tensile specimens, subjected to dynamic loading, may fail by growing a neck outside of the notched region. This apparent lack of sensitivity to a classical stress concentration case was reported but not explained or modeled.The present paper combines experimental and numerical work to address this issue. Specifically, it is shown that the dynamic tensile failure locus is dictated by both the applied velocity boundary condition and the material mechanical properties, specifically strain-rate sensitivity and strain-rate hardening.It is shown that at sufficiently high impact velocities, the flows stress in the notch vicinity becomes quite higher than in the rest of the specimen, so that while the former resists deformation, it transfers the load to the latter. The result will be the formation of a local neck and failure away from the notch.This effect is shown to be active when the material properties are perturbed only at the local level, as in the case of machining of the notch, which in itself may again be sufficient to stabilize the structure under local failure until a neck forms elsewhere.While the physical observations are quite counterintuitive with respect to the engineering views of stress concentrator's effect, the present work rationalizes those observations and also provides information for the designers of dynamically tensioned structures that may contain notches or similar flaws.     

Coupled viscoplasticity damage constitutive model for concrete materials

A coupled viscoplasticity damage constitutive model for concrete materials is developed within the framework of irreversible thermodynamics.Simultaneously the Helmholtz free energy function and a non-associated flow potential function are given, which include the internal variables of kinematic hardening,isotropic hardening and dam- age.Results from the numerical simulation show that the model presented can describe the deformation properties of the concrete without the formal hypotheses of yield criterion and failure criteria,such as the volume dilatancy under the compression,strain-rate sen- sitivity,stiffness degradation and stress-softening behavior beyond the peak stress which are brought by damages and fractures.Moreover,we could benefit from the application of the finite element method based on this model under complex loading because of not having to choose different constitutive models based on the deformation level.     

Dynamic pre-strain and inertia effects on the fracture of metals

A combined experimental and theoretical study is described which examines the influence of strain-rate and dynamic pre-strain on the ductile fracture of thin cylinders. The thin-cylinder configuration is particularly important in this case because it allows inertia terms to be directly incorporated into the theory of plastic instability. A series of quasi-static and dynamic tests is conducted on three materials with differing degrees of strain-rate sensitivity and strain-hardening. The experimental observation that fracture is inhibited at high strain-rates is in accord with the theory when inertia can no longer be considered insignificant. It is also shown that dynamic pre-strain has little or no effect on the flow stress or the strain at fracture in materials which-are essentially strain-rate insensitive, but does reduce the fracture strain in the strain-rate sensitive materials.     

Modeling of Temperature and Strain-Rate Effects in Metals Using an Internal Variable Model

The thermo-viscoplastic behavior of three metals is characterized in a large range of loading conditions by using a new phenomenological constitutive model. The flow stress is decomposed into the sum of an effective stress with an internal stress depending upon an internal parameter which describes the strain hardening effect. The evolution of the internal stress is sensitive to the history of strain-rate and temperature. A systematic method is used for determining the model’s parameters. The model predictions show a good correlation with experimental data. Temperature history effects are especially analyzed.