斜波压缩下锡的相变和层裂行为

获取不同热力学路径下锡的动态响应实验数据,是深入研究其相变和损伤物理过程的基础.利用小型磁驱装置CQ-4完成了金属锡的斜波加载实验,获取了锡含有相变和层裂损伤物理信息的实验数据.实验结果显示,在加载段锡依次经历了弹塑性转变和β-γ相变两种物理过程,屈服强度约0.194 GPa,相变压力随着锡厚度的增加从7.54 GPa减小到7.14 GPa.在卸载段出现了明显的层裂损伤,层裂强度约1.1 GPa,与相同加载压力下冲击实验结果有巨大差异,层裂片厚度约0.38 mm.结合由锡的多相Helmholtz自由能计算的多相状态方程、Hayes相变动力学方程和损伤度理论,对斜波压缩实验过程进行一维流体动力学数值模拟,计算结果可以很好描述锡的弹塑性转变、相变和层裂三个物理过程.     

基于“钉床型”飞片的斜波加载技术及应用

为了实现斜波加载,设计了一种“钉床型”广义波阻抗梯度飞片,即在基座上密排叠加许多小圆锥,简称“钉床型”飞片。该飞片采用激光选区熔化金属增材制造技术进行制备。利用一级轻气炮加载装置和全光纤激光位移干涉测试系统,开展不同工况下“钉床型”飞片高速击靶压缩实验和层裂实验,重点讨论小圆锥高度和撞击速度对斜波压缩加载波形的影响规律,以及斜波加载对不锈钢靶板层裂特性的影响。实验结果显示:（1）“钉床型”飞片对靶板产生的压缩是逐步的,从自由面速度剖面上观察到压缩波上升前沿时间被显著延长,形成了斜波波阵面,明显不同于冲击压缩的陡峭波阵面;（2）在飞片击靶速度近似恒定条件下,斜波波阵面的上升沿时间、平台速度峰值都明显依赖于“钉床型”飞片上的小圆锥高度,随着小圆锥高度增大,上升沿时间呈线性增大,而平台速度峰值呈线性减小;（3）在“钉床型”飞片的几何尺寸保持不变的条件下,斜波波阵面的上升沿时间随着飞片击靶速度的增大而线性减小,平台速度峰值则线性增大;（4）与冲击加载相比,“钉床型”飞片产生的斜波加载不会对材料的层裂强度产生明显影响,但对材料内部损伤演化速率有一定的影响。     

ICE自由面台阶靶数据处理的直接计算方法

斜波压缩台阶靶实验中，不同厚度界面粒子速度历史与材料压缩特性参数存在联系。然而利用普遍采用的数据处理方法无法直接获得该联系。本文中借助特征线理论在建立上述关联的基础上，实现未知EOS下斜波压缩流场的直接计算过程。经数值计算表明，该方法不仅在无强度效应数据处理中能够准确计算理论值，而且在含有强度效应数据处理中也能够较好地逼近理论值，并可在真实实验数据处理中获得与文献符合较好的结果。该研究可为探索具有较完整理论的强度效应数据处理方法提供新途径。     

HMX粒度与爆轰波传播性能的关系研究

采用楔形装药,用临界截面厚度来表征爆轰波传播特性,分别对两种密度条件下HMX粒度对其自身和以HMX为主体的混合炸药HMX/F2641（wHMX/wF2641=95∶5）爆轰波传播性能的影响作了实验研究。结果表明:HMX粒度对爆轰波传播特性有显著的影响,随着HMX粒度的减小,临界截面厚度变小,即爆轰波传播性能增强;同一HMX粒度下, HMX/F2641爆轰波传播性能优于HMX;装药密度的增加有利于爆轰波传递。     

泡沫金属在冲击载荷下的动态压缩行为

基于微CT扫描影像信息,建立泡沫金属材料二维细观有限元模型,考虑不规则胞孔的不均匀分布,根据实验结果拟合孔壁材料的弹塑性本构参数。研究了泡沫金属在不同加载速度下的压缩变形机理,重点讨论泡沫金属中弹塑性波的传播、惯性效应和从冲击端传递到静止端的应力变化特征。对于相对密度为0.3的泡沫铝,弹性波速约为5 km/s,与孔壁材料的弹性波速相当,塑性波速表现为随着加载速度的增大而增大。在加载速度为50~100 m/s间变形模式从准静态模式转变为动态模式,未发现明显的临界速度,动态锁死应变随着加载速度的增大而增大。由于塑性波发生反射,试件会发生二次压缩过程,相应地,静止端产生二次应力平台。受惯性作用的影响,二次应力平台也随着加载速度的增大而提高。     

斜波压缩实验数据的正向Lagrange处理方法研究

提出了一种联合使用Lagrange方法和转换函数方法来处理斜波压缩实验数据的新途径，和传统的数据处理方法相比，该方法更适合处理材料的复杂力学响应并具有较好的精度，同时该方法使用过程中对样品材料参数的初猜值精度要求更低。分析了转换函数方法的可靠性和健壮性，讨论了转换函数和Lagrange方法联合在斜波压缩强度实验数据处理中的应用，获得了可靠的结果。     

动态压缩下Zr基非晶合金失效释能机理

为研究Zr基非晶合金动态压缩条件下的失效释能机理，采用力学试验机、霍普金森杆、高速摄影、差示扫描量热分析（differential scanning calorimetry, DSC）、扫描电镜（scanning electron microscope，SEM）等，得到了材料应力应变曲线、高速摄影图像、失效式样微观形貌及DSC曲线，根据实验数据计算了材料的晶化激活能，并拟合了材料的JH-2（Johnson-Holmquist II）模型，对材料动态失效过程进行有限元数值模拟。实验结果表明，压缩条件下材料为脆性断裂，断口处观察到典型的脉状纹样及液滴状结构，材料失效过程伴随着释能现象；数值模拟结果表明，材料裂纹局部的瞬时内能大于材料晶化激活能。动态压缩下材料的失效释能机理即为材料破碎释放储存的弹性势能，并导致材料局部晶化释能，释能强度与应变率成正相关。     

磁驱动平面准等熵加载装置、实验技术及应用研究新进展

利用脉冲大电流装置产生随时间变化平滑上升的磁压力,实现对平面、柱面等不同结构样品的磁驱动准等熵（斜波）压缩,为极端条件下材料动力学研究提供了一种偏离Hugoniot状态热力学路径的加载手段。本文从磁驱动准等熵加载装置、实验技术、数据处理方法等方面综述了磁驱动准等熵加载技术研究近十年的新进展,评述了利用磁驱动准等熵加载技术和方法开展极端条件下材料高压状态方程、高压强度与本构关系、相变与相变动力学等方面研究的进展情况,展望了磁驱动准等熵加载技术发展及其在材料动力学、武器物理和高能量密度物理等方面的应用前景。     

两参数轴向冲击载荷作用下圆柱壳弹塑性动力屈曲

研究圆柱壳在两参数轴向冲击载荷下的弹塑性动力屈曲问题，基本控制方程由弹塑性连续介质中关于加速度的最小原理获得，本构关系采用增量理论。研究表明：屈曲过程可划分为两相，两相之间由临界时间ｔ表征，并分别讨论了应力波对屈曲的影响，压缩波与弯曲波的相互作用及几何尺寸，材料参数，初始缺陷，载荷峰值及持续时间等诸多因素与动力屈曲的关系。     

平面应变各向异性本构关系及在应力波传播模拟中的应用

为了研究平面应变条件下各向异性材料中应力波传播的特点,利用各向异性弹性Hooke定律、 Tsai-Hill屈服准则、经典塑性流动理论,引入修正的物态方程计及高压下的体积压缩非线性,建立了平面应 变条件下正交各向异性复合材料的弹塑性本构关系,并且分析了二维问题中材料变形引起的主轴旋转及客 观应力率修正问题。最后采用动态显式有限元方法自行编写程序模拟某种纤维增强复合材料碰撞过程中平 面应力波的传播,模拟结果显示,在平面应变条件下应力波在该材料的传播过程中表现出明显的二维效应、各 向异性特点及弹塑性特点。     

Phase transition and dynamics of iron under ramp wave compression

The ramp wave compression experiments of iron with different thicknesses were performed on the magnetically driven ramp loading device CQ-4. Numerical simulations of this process were done with Hayes multi-phase equation of state (H-MEOS) and dynamic equations of phase transition. The calculated results of H-MEOS are in good agreement with those of shock phase transition, but are different from those under ramp wave compression. The reason for this is that the bulk modulus of the material in the Hayes model and the wave velocity are considered constant. Shock compression is a jump from the initial state to the final state, and the sound speed is related to the slope of the Rayleigh line. However, ramp compression is a continuous process, and the bulk modulus is no longer a constant but a function of pressure and temperature. Based on Murnaghan equation of state, the first-order correction of the bulk modulus on pressure in the Hayes model was carried out. The numerical results of the corrected H-MEOS agree well with those of pure iron in both ramp and shock compression phase transition experiments. The calculated results show that the relaxation time of iron is about 30 ns and the phase transition pressure is about 13 GPa. There are obvious differences between the isentropic and adiabatic process in terms of pressure–specific volume and temperature–pressure. The fluctuation of the sound speed after 13 GPa is caused by the phase transition.     

基于分子动力学模拟的单晶硅冲击压缩相变研究

晶体硅具有复杂的相变机制,在相图研究中受到广泛关注,其在动载荷下的变形机制是当前研究热点。为揭示晶体硅在强动加载下的变形和相变行为特征,基于分子动力学方法,采用平板冲击加载方式,模拟研究了单晶硅在初始环境温度为300 K时分别沿[001]、[110]和[111]晶向的不同强度下的冲击压缩行为,冲击粒子速度为0.3～3.2 km/s。研究发现,随着冲击粒子速度的增加,单晶硅剪切应力在逐渐增加后由于结构相变发生急剧下降,相变阈值和相变机制均呈现各向异性。其中,沿[001]晶向冲击压缩下观察到多种固-固相变以及固-液相变,并观察到与最新文献的实验高度一致的固-液共存现象。研究结果可为动加载下晶体硅的相变研究提供纳米尺度的结果支撑。     

Theoretical latent hardening in crystals—I. General equations for tension and compression with application to F.C.C. crystals in tension

Equations for latent strengths in single slip, based upon the simple theory of finite distortional crystal hardening introduced by K.S. Havner and A.H. Shalaby (1977), are derived for both tensile and compression tests without restriction as to crystal class. Detailed comparisons between theoretical results and the experiments of P.J. Jackson and Z.S. Basinski (1967) on copper crystals in tension are presented. There is good qualitative agreement between theory and experiment regarding the diversity of anisotropic hardening among slip systems. Moreover, there is satisfactory quantitative agreement between the theory and the extrapolated experimental data in the stage III, large-strain range. It is suggested that further experimental investigation of latent hardening at large prestrains would be desirable.The simple theory predicts anisotropic hardening and the perpetuation of single slip in axial loading of cubic crystals initially oriented for single slip, but predicts symmetric, isotropic hardening of specimens initially oriented in positions of 4, 6 or 8-fold multiple-slip. These predictions are in general accord with experimental observations from tests of f.c.c. and b.c.c. crystals.     

非均质结构PBX炸药的动态压缩过程模拟

采用离散元方法,构建了考虑PBX炸药晶体颗粒与黏结剂非均匀性的计算模型,通过在炸药试样两端采用相向飞片加载、将原SHPB实验的应力边界替代为速度边界的方法,开展了PBX炸药的动态压缩数值模拟研究,再现了考虑PBX炸药非均质特性的SHPB实验过程,获得了不同加-卸载路径下炸药应力应变曲线以及相应的损伤破坏图像。模拟结果表明:在动态压缩过程中,虽然PBX炸药处于整体应力平衡,由于PBX高度非均质性,内部应力分布并不均匀,晶体间应力以应力桥形式传递;应力曲线为试样整体平均,试样内局部所受真实应力可能高于曲线峰值应力;出现卸载回滞曲线的试样发生局部损伤破坏,而完全软化曲线的试样发生整体失稳破坏。     

PMN-PT夹层弹性半空间结构中SH波的传播Propagation of SH waves in sandwich structures consisting of PMN-PT single crystal layer and elastic half-spaces

分析了弹性上下半空间和PMN‐PT单晶层组成的夹层结构中SH波的传播性质,PMN‐PT单晶沿[011]c方向极化,宏观上呈mm2对称,且晶体沿角度θ方向切割。基于正交各向异性压电材料和各向同性弹性材料的基本方程,得到了夹层结构中SH波传播时行列式形式的频散方程。通过对数值算例进行分析可以看出,PMN‐PT单晶的切割角度和弹性材料属性对结构中的相速度有很大影响,因此波的某些传播性能可以通过材料的设计以及晶体切割的方向来实现,这些结论为声表面波器件的开发和应用提供了理论依据。     

双屈服面法测量金属材料动高压屈服强度的若干改进

在J.R.Asay等提出的双屈服面法（简称AC方法）测量动态屈服强度的基础上,对实验设计和数据处理做了部分改进。在实验设计方面,在满足实验物理设计的前提下,将 AC方法的卸载实验装置(双层飞片和双台阶样品)简化为单飞片和单台阶样品,以简化波系作用,并减少了声速测量不确定度的引入因素。在数据处理方面,发现由Lagrange声速和工程应变关系曲线外推求预冲击态的体积声速容易引入较大的不确定度,而通过Lagrange声速和粒子速度关系曲线可以更清晰地得到塑性卸载起始点,使图解法和积分法得到的屈服强度值趋于一致。利用改进后的AC方法,得到了低温退火后LY12铝在20.0 GPa预冲击压力下的屈服强度为0.6 GPa。     

Thermal and mechanical analysis of material response to non-steady ramp and steady shock wave loading

Ramp wave experiments on the Sandia Z accelerator provide a new approach to study the rapid compression response of materials at pressures, temperatures and stress or strain rates not attainable in conventional shock experiments. Due to its shockless nature, the ramp wave experiment is often termed as an isentropic (or quasi-isentropic) compression experiment (ICE). However, in reality there is always some entropy produced when materials are subjected to large amplitude compression even under shockless loading. The entropy production mechanisms that cause deformation to deviate from the isentropic process can be attributed to mechanical and thermal dissipations. The former is due to inelasticity associated with various deformation mechanisms and the rate effect that is inherent in all the deformation processes and the latter is due to irreversible heat conduction. The main purpose of the current study is to gain insights into the effects of ramp and shock loading on the entropy production and thermomechanical responses of materials. Another purpose is to investigate the role of heat conduction in the material response to both the non-steady ramp wave and steady shock.Numerical simulations are used to address the aforementioned research objectives. The thermomechanical response associated with a steady shock wave is investigated first by solving a set of nonlinear ordinary differential equations. Using the steady wave solutions as the reference, the material responses under non-steady ramp waves are then studied with numerical wave propagation simulation. It is demonstrated that the material response to ramp and shock loading is essentially a manifestation of the interaction between the time scale associated with the loading and the intrinsic time scales associated with mechanical deformation and heat transfer. At lower loading rates as encountered in ramp loading, the loading path is closer to an isentrope and results in lower entropy production. The reasonable ramp rate to obtain a quasi-isentropic state depends on the intrinsic time scales of the dissipation mechanisms which are strongly material dependent. Thus shockless loading does not necessarily produce an isentropic response. Between two equilibrium states, heat conduction was shown to have significant effect on the temperature history but it contributes little to the overall temperature change if the specific heat remains constant. It also affects the history of entropy, but only the irreversible part of heat conduction contributes to the net entropy change. The various types of thermomechanical responses of materials would manifest themselves more significantly in terms of the thermal history than the mechanical history. Thus temperature measurement appears to be an important experimental tool in distinguishing the various mechanisms for the thermomechancial responses of the materials.