爆炸载荷下脆性颗粒体系破碎特性的数值研究

试验发现,以球形TNT为中心爆源,球形玻璃珠构成的颗粒和球壳中发生破碎的颗粒体积分数随当量比(颗粒球壳的质量与TNT炸药的质量比)的增加呈现指数衰减规律。采用有限元与离散元耦合的连续非连续数值方法,揭示了中心炸药起爆后颗粒环壳内爆炸波的传播衰减和在环壳外界面反射后的稀疏卸载过程。由于爆炸波的短脉冲特性,颗粒内部应力场始终处于应力非均衡状态,采用应力均衡状态下颗粒破碎强度的Weibull分布会得到远高于试验测得的破碎颗粒体积分数。因此采用破坏波传播特征时间内的平均诱导应力而非瞬时诱导应力作为颗粒破碎强度的应力指标,并通过试验结果确定破坏波传播特征时间。考虑了应力传播的非均匀性对于颗粒破碎的影响,得到了平均诱导应力峰值的概率分布随比例距离的变化规律,结合修正后的颗粒破碎强度Weibull分布建立了破碎颗粒体积分数随比例距离的变化模型。     

滑坡冲击铲刮变量的计算方法研究

高速远程滑坡运动过程中,冲击铲刮效应不仅增加滑坡的体积与规模,而且会增大滑坡成灾范围,风险预测与判断出现明显误差,导致灾难性事件的发生。目前对于高速滑坡的铲刮深度、铲刮范围和铲刮体积等变量计算往往是采用基于经验的铲刮率算法,其是通过体积增量来反算铲刮变量的数学方法,而实际情况中高速滑坡的冲击铲刮是滑体冲击力和地表可铲刮材料之间的力学屈服破坏作用的结果。本文基于接触力学、弹塑性力学和岩土力学,提出了地表可铲刮层在附加冲击荷载作用下铲刮变量的理论计算方法,认为竖向冲击和切向剪切是滑体冲击铲刮过程的两种主要作用方式。结合实际岩土体材料参数,发现不考虑切向荷载时,铲刮层塑性区主要以竖向破裂区为主,考虑切向荷载时,铲刮层塑性区沿切向迁移,直至塑性边界贯通至地表,符合实际情况,实现了对铲刮变量的定量化计算。     

基于FDM-DEM耦合的冲击损伤大理岩静态断裂力学特征研究

为探究循环冲击损伤后大理岩的静态断裂力学特征,基于有限差分（finite difference method,FDM）-离散元（discrete element method,DEM）耦合的建模技术构建了三维分离式霍普金森压杆（split Hopkinson pressure bar,SHPB）数值模型,其中杆件系统和岩石试件分别采用FLAC3D和PFC3D程序建模。利用该模型对中心直切槽半圆盘（NSCB）试样进行了恒定子弹速度下的循环冲击,随后对受损试样进行静态三点弯曲断裂实验。通过编写Fish程序,提取试样断裂面数据,对断裂面进行重构并定量计算表面粗糙度。通过与相关室内实验结果的对比分析,验证了本文数值分析的合理性与可靠性。模拟结果表明,随着循环冲击次数的增加,试样内部微裂纹、破碎颗粒均增加。连接力场分布混乱,部分力链发生断裂。力链的变化是试样力学性能劣化的根本原因。在静态三点弯曲断裂实验中,冲击5次后试样的静态断裂韧度较天然试样产生一定程度的降低。试样在静载过程中产生的微裂纹和碎块的数量随循环冲击次数的增加而增加,断裂面粗糙度随循环冲击次数的增加而增加。     

[100]LiF单晶在60 GPa以内冲击加载下的高压屈服强度特性

获取光学窗口自身的高压强度特性是开展材料高压高应变率冲击响应行为精密测量和数据反演的重要基础。利用平板撞击和双屈服面法,通过冲击-卸载、冲击-再加载原位粒子速度剖面精细测量和数据反演,获得了约60 GPa范围内[100]LiF屈服强度特性随冲击压力的变化规律。结果表明:在实验压力范围内,[100]LiF的屈服强度随加载压力的提高而显著提高,压力硬化效应显著;同时,LiF在冲击加载下的屈服强度高于磁驱准等熵加载结果,应变率硬化效应强于热软化效应。采用Huang-Asay模型确定了可描述冲击加载[100]LiF强度特性的本构模型参数,为LiF在强度、相变、层断裂等加窗测量实验中的深入应用和数据准确解读提供了重要支撑。     

一种水下爆炸冲击波压力调控方法

起爆位置和装药形状对水下爆炸冲击波压力有较为显著的影响,这使得利用小当量装药在局部方向形成与大当量装药一定程度等效的冲击波成为可能。为了能够在小当量装药条件下开展舰船结构及设备抗水下爆炸冲击实验,基于细长装药结构和参数优化设计,设计了一种冲击波压力幅值和持续时间可调的装药方法。首先,基于简单波理论给出了水下爆炸冲击波压力调控的原理,以及装药参数优化设计的目标函数和约束条件;然后,采用自主数值模拟软件研究了细长装药的水下爆炸能量输出规律,通过实验验证了数值模拟的置信度,研究发现起爆位置和装药形状对水下爆炸冲击波压力峰值和持续时间的影响是显著的,在炸药爆速一定的情况下,长药柱水下爆炸冲击波压力的持续时间可通过几何近似确定;最后,为了进一步考察该方法的有效性,以1000 kg TNT和100 m爆距的水下爆炸冲击波压力-时间曲线作为原型,设计了2种与该原型冲击波压力等效的装药方案,并通过数值模拟进行了验证。研究结果表明:设计的装药能够在预定的持续时间内,在装药起爆端一侧形成与原型等效的冲击波压力-时间曲线。由于没有考虑对气泡载荷的等效,因此该调控方法仅适用于中远场爆炸冲击问题。     

Al/PTFE活性材料冲击载荷作用下响应特性研究

为了研究铝（Al）/聚四氟乙烯（polytetrafluoroethylene, PTFE）活性材料冲击载荷作用下响应特性,制备了具有反应活性的Al/PTFE块体材料,设计了拉氏实验,采用不同厚度的铝隔板控制入射冲击波幅值,利用锰铜压阻计测量了冲击波在材料中传播过程压力演化过程。同时,基于AUTODYN有限元软件,采用Lee-Tarver三项式点火模型对Al/PTFE活性材料拉氏实验进行数值模拟,并探讨了冲击波在500 mm长的Al/PTFE活性材料中长距离传播行为。研究结果表明,冲击波压力在Al/PTFE活性材料内短距离传播过程中存在明显的衰减,但是,当冲击波传播到远距离时,冲击波压力幅值和冲击波速度趋于稳定,分别为1.3 GPa和2 180 m/s;同时,距离铝隔板越远的材料,其反应度越低并最终趋于0.17。正是由于材料化学反应释能,导致了冲击波压力传播过程最终趋于稳定状态。     

纯铁和铁基合金楔形样品中的冲击相变

2005年，Koflov等对滑移爆轰加载下的楔形纯铁DT-2样品进行了研究，样品仅在斜面部分发生了局部层裂，材料中明显存在不同的几个区域：材料无破坏区域、层裂开始和发展区域、加载表面-高能炸药接触面、满足α-ε相变条件的加载表面邻近区域和初始状态为硼发生高应变率变形的材料区域。但是滑移爆轰加载波系作用非常复杂，无法准确给出各区域宏细观变形破坏所对应的加载应力状态。     

用作VISAR窗口的LiF晶体折射率变化修正因子

近几年，因为样品／窗口界面处理工艺的完善和窗体寄生干涉问题的解决，加窗VISAR测试技术在冲击波物理研究中得到了愈来愈广泛的应用。但窗口材料在冲击作用（如压缩、拉伸、加热等）下，由于折射率变化将引入附加多普勒频移，从而对最终的测速结果引入修正项。此修正项与窗口材料的折射率变化特性直接相关，而且对某些窗口材料，此修正项值还比较大。因此，为获取加窗测试中正确的速度剖面测量数据，必须确定窗口材料在冲击作用下的折射率变化修正因子。在目前所用VISAR窗口中，LiF晶体因为具有中等阻抗，     

Compton散射下强激光等离子体空气冲击波波前的传输特性

 根据球对称激光等离子体空气冲击波在自由空间中传输的约束条件，对多光子非线性Compton散射的强激光等离子体空气冲击波波前的传输特性进行了研究，给出了散射下空气冲击波波前的运动方程，并进行了数值模拟。结果表明：该冲击波的衰减过程不仅与爆炸源和爆炸过程的特性、释放总能量、空气的弹性有关，而且还与散射有关，散射效应使冲击波的初始半径增大，衰减过程加快，能量转移率提高；数值模拟与实验结果符合得很好。     

立方非线性产品包装损坏边界的三维曲面表示

用三次Ｂ样条函数配点法，求得在矩形波冲击下具有立方非线性产品包装的冲击谱。讨论了软、硬非线性及外界激励幅值对产品损坏边界的影响。提出用损坏边界曲面来描述非线性产品包装的安全区域，并作出了相应的损坏边界曲面．