激波与堆积粉尘相互作用的数值模拟

基于双流体模型和测定的堆积粉尘的本构方程 ,利用AUSM+ 格式 ,对激波与堆积粉尘的相互作用进行了数值模拟。计算所反映的流场结构与实验图像一致。此外还对激波强度 ,颗粒材料密度等对流场的影响进行了讨论。     

On the response of rubbers at high strain rates—III. Effect of hysteresis

In this series of papers, we examine the propagation of waves of finite deformation in rubbers through experiments and analysis; in the present paper, Part III, the effect of hysteretic material behavior on the free retraction of prestretched rubber is considered. A rubber strip stretched to many times its initial length is released at one end and the resulting unloading is examined. A high-speed video camera was used to monitor the motion and to determine the evolution of strain and particle velocity in rubber strips. Simple waves as well as shock waves are observed in these unloading experiments. The measurements are modeled using a power-law model of the material behavior. The hysteretic material response and the formation of shocks are characterized.     

Full-Field Deformation Measurements of Aluminum Plates Under Free Air Blast Loading

Assessment of the structural dynamic response caused by an explosion is complex due to the high velocity impact, the transient nature of the deformation and the interactions between the structure and the pressure wave. This paper deals with full-field measurements of aluminum plates under free air blast loading conditions. Forty grams of explosive material C4 is detonated at a stand-off-distance of 250?mm and two synchronized high-speed cameras in a stereoscopic setup are used to capture the plate response with an inter frame rate of 6,000 fps. The transient deformation fields are calculated using a three-dimensional digital image correlation technique. The observations appeared to be interesting and somewhat counter intuitively. Results show that a free air blast load induces a highly localized, rapid material response which can be essentially divided in two different stadia. First, when the shock impulse occurs, all particles are forced to move out-of-plane and provided with initial velocities. Secondly, when the pressure wave has vanished, the deformation is further driven by this imparted momentum. This paper shows that a 3D high-speed DIC system is a powerful tool for the assessment of the dynamic response of a structure subjected to extreme loading conditions such as explosions and that this system is capable of accurately measuring surface displacement and deformation data at high rates. Moreover, a free air blast load makes it possible to load a plate specimen at different strain rates in different zones. This makes the test suitable for future material identification using inverse methods, which profit from heterogeneous displacement and strain fields.     

Particle response to shock waves in solids: dynamic witness plate/PIV method for detonations

Studies using transparent, polymeric witness plates consisting of polydimethlysiloxane (PDMS) have been conducted to measure the output of exploding bridge wire (EBW) detonators and exploding foil initiators (EFI). Polymeric witness plates are utilized to alleviate particle response issues that arise in gaseous flow fields containing shock waves and to allow measurements of shock-induced material velocities to be made using particle image velocimetry (PIV). Quantitative comparisons of velocity profiles across the shock waves in air and in PDMS demonstrate the improved response achieved by the dynamic witness plate method. Schlieren photographs complement the analysis through direct visualization of detonator-induced shock waves in the witness plates.     

机动方舱核爆炸冲击波动力响应和安定性分析

本文以核爆炸当量、高度、距离为参数，采用等超压圆给出了普通方舱和防核加固方舱结构在核冲击波下的安全域。研究了限制舱体倾斜的压力条件。本文还以冲击波压力的时间函数为激励，求解了方舱的动力响应。从试验结果推论和计算模拟两方面，分析了方舱承受冲击波的能力。     

Shock tube study of the dynamical behavior of granular materials

The purpose of the present study is to clarify both the compression phenomenon and the gas filtration effect that take place inside a granular medium when it is dynamically loaded by a shock wave. In order to measure the pore pressure and the total stress at different locations along the granular medium, pressure transducers were placed along the side-wall and at the end-wall of the shock tube test section, which was filled with the granular material. In order to elucidate the gas filtration effect, the results of two experiments with identical granular media but with and without filtration were compared. The gas filtration was eliminated by means of a thin plastic film, which was placed at the front edge of the granular medium. Based on this comparison quantitative information on the gas filtration and its role in the stress formation inside granular media of different material and length was obtained. Furthermore, curves of the dynamic compression and the Young moduli of the granular medium for the range of the operating conditions were reconstructed.     

轴、侧向同卸荷下砂岩力学特性影响的试验研究

本文以实际岩体工程为背景,利用WDT-1500 仪器开展了轴向、侧向同时卸荷条件下砂岩的三轴试验. 结果表明:轴、侧向同卸荷这种卸荷路径下,砂岩试样破坏时并没有出现应力峰值,为了定义试样的破坏强度,将最大与最小主应力差随最小主应力的变化关系曲线上应力跌落的拐点处的应力值定义为破坏强度. 砂岩变形初始段发生应力跌落和轴向应变回弹,破坏前无明显的弹性和屈服阶段;试验的过程中,砂岩的侧向变形明显大于轴向变形,其体积应变一直处于膨胀状态;相对于砂岩的常规三轴试验结果,试样破坏时的强度在轴向、侧向同时卸荷条件下有所降低. 初始轴压和初始围压对试样的力学特征有十分显著的影响,但围压的卸荷速率却并不显著. 砂岩的破坏特征主要是以张-拉为主的混合张剪的破坏.      

小不耦合系数装药爆破孔壁压力峰值计算方法

基于轮廓爆破孔壁压力峰值计算方法的相关研究,充分考虑空气冲击波的传播与爆轰产物膨胀的过程,理论分析了小不耦合系数装药爆破过程中空气冲击波与炮孔壁的相互作用,建立了三维空气介质径向不耦合装药单孔爆破有限元模型,研究了工程爆破中常用的多种小不耦合系数装药组合工况下,炸药单点起爆后的炮孔壁压力峰值,并获得了相应工况下的孔壁压力峰值较爆生气体准静态等熵膨胀压力的压力增大倍数。结果表明:小不耦合系数装药爆破过程中,爆轰产物参数会对空气冲击波波后物质参数产生显著影响,揭示了小不耦合系数装药爆破与轮廓爆破在孔壁压力峰值计算方法上的本质差异;柱状装药结构爆轰波沿轴向传播使得空气冲击波撞击炮孔壁时存在叠加效应,孔壁压力峰值也相应增大,通过统计分析不同炸药类型、不同岩石类型工况下压力增大倍数与不耦合系数的关系,发现压力增大倍数随不耦合系数的增大近似呈线性增长;基于理论推导结果及常用爆破孔壁压力峰值计算形式,综合考虑炸药性能、孔壁岩石介质条件、不耦合装药系数对空气冲击波撞击炮孔壁后压力增大倍数的影响,提出了不耦合系数较小时爆破孔壁压力峰值计算方法。     

Compression and rarefaction waves in shock-compressed metals

The behavior of duralumin and copper is studied under conditions of specimen loading by two successive shock waves and during unloading after the shock compression. The amplitude of the first shock wave was 150–250 kbar. Direct measurements were performed of the difference in main stresses behind the shock front in duralumin. The results obtained do not agree with existing concepts of the behavior of solids under dynamic loading. Possible causes of this divergence are considered.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 2, pp. 146–153, March–April, 1976.In conclusion, the authors thank G. A. Savel'ev for his aid in the preparation and performance of the experiments.     

颗粒材料中致密波结构研究

采用一维两相流模型与相应颗粒构形应力函数,研究了致密波的形成及其结构．用简化两相流模型系统地讨论致密波对有关因素的依赖关系．分析指出：小于基体材料音速的致密波仅能在非理想颗粒材料中存在,从波前到波后,所有状态物理量光滑过渡．大于基体材料音速的致密波,波头可能存在间断．应力函数与致密粘性确定后,致密波速度决定致密波结构、宽度、终态压实度．采用一维两相流模型模拟了活塞驱动颗粒床形成致密波这一动态过程．用线方法（MOL）对该方程组求数值解．计算表明,经过短暂的非稳态过程,颗粒床中形成一稳态致密波．分析了活塞速度与初始孔隙率对致密波结构的影响,并对简化两相流模型与两相流模型的计算结果进行了对比．     

Analytical Solution for Isothermal Flow in a Shock Tube Containing Rigid Granular Material

Analytical solution of shock wave propagation in pure gas in a shock tube is usually addressed in gas dynamics. However, such a solution for granular media is complex due to the inclusion of parameters relating to particles configuration within the medium, which affect the balance equations. In this article, an analytical solution for isothermal shock wave propagation in an isotropic homogenous rigid granular material is presented, and a closed-form solution is obtained for the case of weak shock waves. Fluid mass and momentum equations are first written in wave and (mathematical) non-conservation forms. Afterwards by redefining the sound speed of the gas flowing inside the pores, an analytical solution is obtained using the classical method of characteristics, followed by Taylor’s series expansion based on the assumption of weak flow which finally led to explicit functions for velocity, density and pressure. The solution enables plotting gas velocity, density and pressure variations in the porous medium, which is of high interest in the design of granular shock isolators.     

氢氧燃烧及爆轰驱动激波管

分析并观察了沿驱动段轴向分布多火塞燃烧驱动段的性能．提出主膜处同一管截面均匀分布三火花塞引燃的点火方法．用这种点火方法驱动产生的入射激波强度重复性较高,激波后气流速度、温度和压力的定常性亦大大改善,可满足气动试验实际要求．提出在驱动段尾端串接卸爆段来消除爆轰波反射高压,从而可使反向爆轰驱动段用来产生高焓高密度试验气流．这种反向爆轰驱动产生的入射激波重复性高,激波衰减弱．在主膜处的收缩段产生的反射波可缓解爆轰波后跟随的稀疏波的不利影响,从而使前向爆轰驱动具有实用性．在产生的入射激波强度相同条件下,前向爆轰驱动所需的爆轰驱动段可爆混合气初始压力可较反向爆轰低近一个量级．     

散体介质SHPB被动围压试验体应力计算的理论修正方法

SHPB被动围压试验为探究散体介质在爆炸和冲击荷载作用下的力学行为提供了一个行之有效的方法。针对相关试验设计和计算中存在的弊端和不足,借助经典板壳理论将SHPB被动围压试验中用于约束散体介质的刚性套筒简化为受均匀带状内压作用的圆柱形壳体。理论计算了套筒径向位移、环向应变与均匀带状内压及套筒几何、力学参数的关系,得到了套筒径向位移、环向应变沿其轴向的分布规律;分析了套筒长度、厚度、内外径以及均匀带状内压宽度之间等无量纲几何参数对计算结果的影响;将理论计算结果与试验和数值模拟结果进行对比,验证了理论计算结果的正确性。本文中提出的理论修正方法可为指导散体介质SHPB被动围压试验提供参考。     

Effects of chamber temperature and pressure on the characteristics of high speed diesel jets

This study is an investigation into the effects of temperature and pressure within a test chamber on the dynamic characteristics of injected supersonic diesel fuel jets. These jets were generated by the impact of a projectile driven by a horizontal single stage powder gun. A high speed video camera and a shadowgraph optical system were used to capture their dynamic characteristics. The test chamber had controlled air conditions of temperature and pressure up to 150 °C and 8.2 bar, respectively. It was found experimentally that, at the highest temperature, a maximum jet velocity of around 1,500 m/s was obtained. At this temperature, a narrow pointed jet appeared while at the highest pressure, a thick, blunt headed jet was obtained. Strong shock waves were generated in both cases at the jet head. For analytical prediction, equations of jet tip velocity and penetration from the work of Dent and of Hiroyasu were employed to describe the dynamic characteristics of the experiments at a standard condition of 1 bar, 30 °C. These analytical predictions show reasonable agreement to the experimental results, the experimental trend differing in slope because of the effect of the pressure, density fluctuation of the injection and the shock wave phenomena occurring during the jet generation process.     

点火能对丙烷-空气预混气体爆炸过程及管壁动态响应的影响

在长12 m的无缝不锈钢直管中,通过改变初始点火能量,探究了点火能对封闭管道内丙烷-空气混合气体爆炸传播特性和激波对管壁动态加载的影响。结果表明,初始点火能对预混气体爆炸火焰传播规律以及管壁的动态响应有显著影响:点火能越大,爆炸越剧烈,爆炸压力峰值压力和管壁最大应变就越大,且压力波和管壁应变的发展一致。火焰在传播过程中受到管道末端反射波的作用会发生短暂熄灭和复燃;管壁承受冲击波加载,应变信号主要分布在0~781.25 Hz,管壁最大应变率大于10-3 s-1,实验工况下管壁应变属动态响应。     

Mesoscale calculations of the dynamic behavior of a granular ceramic

Mesoscale calculations have been conducted in order to gain further insight into the dynamic compaction characteristics of granular ceramics. The primary goals of this work are to numerically determine the shock response of granular tungsten carbide and to assess the feasibility of using these results to construct the bulk material Hugoniot. Secondary goals include describing the averaged compaction wave behavior as well as characterizing wave front behavior such as the strain rate versus stress relationship and statistically describing the laterally induced velocity distribution. The mesoscale calculations were able to accurately reproduce the experimentally determined Hugoniot slope but under predicted the zero pressure shock speed by 12%. The averaged compaction wave demonstrated an initial transient stress followed by asymptotic behavior as a function of grain bed distance. The wave front dynamics demonstrate non-Gaussian compaction dynamics in the lateral velocity distribution and a power-law strain rate–stress relationship.     

DEVELOPMENT OF EXPERIMENTAL METHODS FOR IMPACT TESTING BY COMBINING HOPKINSON PRESSURE BAR WITH OTHER TECHNIQUES

The split Hopkinson pressure bar（SHPB） technique and the wave propagation inverse analysis（WPIA） technique are both extensively used to experimentally investigate the impact behavior of materials, although neither of them alone provides a fully satisfactory analysis. In the present paper, attention is given to new experimental techniques by incorporating a damagemodified constitutive model into the SHPB technique and combining the Hopkinson pressure bar（HPB） technique with WPIA. First, to distinguish the response due to dynamic constitutive behavior and the response due to dynamic damage evolution, the SHPB method incorporating a damage-modified constitutive model is developed, including an explicit damage-modified Zhu–Wang–Tang model and an implicit damage-modified constitutive model. Second, when the SHPB results become invalid, a method of combining new Lagrange inverse analyses with the HPB technique is developed, including cases of the HPB arranged in front of a long specimen and behind the specimen. As examples of these new methods, typical results are given for nonlinear viscoelastic polymers and concretes considering damage evolution, a super-elastic Ti–Ni alloy with phase transformation and an aluminum foam with shock waves propagating within it.     

爆轰加载下金属锡层裂破碎数值模拟

对爆轰加载下低熔点金属锡的层裂破碎问题开展了数值模拟。在利用实验数据对所采用数值方法和材料模型开展对比验证的基础上,通过对样品内部物理量时间及空间分布演化对比分析,剖析了冲击加-卸载中样品内部应力波与材料相互作用过程。此外,通过对比分析不同厚度锡样品在爆轰加载下的动态行为特征,进一步认识了自由面反射稀疏波、边侧稀疏波和入射稀疏波共同作用下层裂破碎演化机制。结果表明,当样品较薄时,层裂破碎行为由反射稀疏波主导;随着样品厚度的增大,反射稀疏波主导区缩小,入射稀疏波和边侧稀疏波主导区逐渐增大。     

Evolution of blast wave profiles in simulated air blasts: experiment and computational modeling

Shock tubes have been extensively used in the study of blast traumatic brain injury due to increased incidence of blast-induced neurotrauma in Iraq and Afghanistan conflicts. One of the important aspects in these studies is how to best replicate the field conditions in the laboratory which relies on reproducing blast wave profiles. Evolution of the blast wave profiles along the length of the compression-driven air shock tube is studied using experiments and numerical simulations with emphasis on the shape and magnitude of pressure time profiles. In order to measure dynamic pressures of the blast, a series of sensors are mounted on a cylindrical specimen normal to the flow direction. Our results indicate that the blast wave loading is significantly different for locations inside and outside of the shock tube. Pressure profiles inside the shock tube follow the Friedlander waveform fairly well. Upon approaching exit of the shock tube, an expansion wave released from the shock tube edges significantly degrades the pressure profiles. For tests outside the shock tube, peak pressure and total impulse reduce drastically as we move away from the exit and majority of loading is in the form of subsonic jet wind. In addition, the planarity of the blast wave degrades as blast wave evolves three dimensionally. Numerical results visually and quantitatively confirm the presence of vortices, jet wind and three-dimensional expansion of the planar blast wave near the exit. Pressure profiles at 90° orientation show flow separation. When cylinder is placed inside, this flow separation is not sustained, but when placed outside the shock tube this flow separation is sustained which causes tensile loading on the sides of the cylinder. Friedlander waves formed due to field explosives in the intermediate-to far-field ranges are replicated in a narrow test region located deep inside the shock tube.     

Effective simulation of flexible lateral boundaries in two- and three-dimensional DEM simulations

Discrete element method （DEM） models to simulate laboratory element tests play an important role in advancing our understanding of the mechanics of granular material response, including bonded or cemented, particulate materials. Comparisons of the macro-scale response observed in a real physical test and a ＂virtual＂ DEM-simulated test can calibrate or validate DEM models. The detailed, particle scale information provided in the DEM simulation can then be used to develop our understanding of the material behaviour. It is important to accurately model the physical test boundary conditions in these DEM simulations. This paper specifically considers triaxial tests as these tests are commonly used in soil mechanics. In a triaxial test, the test specimen of granular material is enclosed within a flexible latex membrane that allows the material to deform freely during testing, while maintaining a specified stress condition. Triaxial tests can only be realistically simulated in 3D DEM codes, however analogue, 2D, biaxial DEM simulations are also often considered as it is easier to visualize particle interactions in two dimensions. This paper describes algorithms to simulate the lateral boundary conditions imposed by the latex membrane used in physical triaxial tests in both 2D and 3D DEM simulations. The importance of carefully considering the lateral boundary conditions in DEM simulations is illustrated by considering a 2D biaxial test on a specimen of frictional unbonded disks and a 3D triaxial test on a bonded （cemented） specimen of spheres. The comparisons indicate that the lateral boundary conditions have a more significant influence on the local, particle-scale response in comparison with the overall macro-scale observations.