运动装药空中爆炸冲击波特性研究

利用AUTODYN对不同速度装药的空中爆炸冲击波场进行数值计算,定量研究运动装药空中爆炸冲击波场特性。仿真结果表明:装药速度对冲击波场的时空分布有较大影响,爆炸初期冲击波场的移动速度和位移与装药速度正相关;动爆冲击波场为空间非均匀分布,超压值随着与装药运动方向夹角的增大近似呈余弦变化,由大于静爆超压值逐渐降低为小于静爆超压值,且变化幅度均随装药速度的增加而增大。最后,分析了动爆与静爆冲击波场的关联特性,建立了动爆冲击波超压的工程计算模型,模型计算结果与动爆试验结果及仿真计算结果吻合较好。     

反恐中微差控爆破门实验研究

 在解救人质或追捕恐怖分子时,需要利用炸药炸开门体或者墙体开辟一条快速通道,通常采用集团装药或多点药包同时起爆的方法进行爆破,但产生的冲击波超压值过高,很容易对反恐人员造成不必要的伤害。采用多点延时起爆方法进行微差控爆破门,与其它方法相比,在相同药量和距离下可以降低空气冲击波的超压峰值,从而有效避免冲击波对反恐人员的伤害。实验研究结果表明,防盗门中的膨胀螺栓可以作为破门弱点进行毁伤。通过实验,得到了破门的合理药量,利用不同药量下的冲击波超压值计算出了反恐人员的安全距离。     

地下管廊在燃气爆炸作用下的动力响应分析

在燃气通过地下管廊输送过程中,若燃气泄漏进入管廊内部并引起爆炸将会产生严重的后果。以平潭综合试验区环岛路管线工程为背景,借助LS-DYNA非线性动力分析有限元软件,基于流固耦合和ALE多物质算法,建立地下空间燃气爆炸数值计算模型,探讨管廊内燃气爆炸过程中爆炸冲击波对管廊结构的影响。数值计算结果表明:在爆炸过程中,超压峰值从爆源中心转移到波阵面上,导致爆源附近压力低于波阵面压力,爆源附近形成负压区;爆炸荷载作用下,燃气仓内距离爆源最近的内墙上的测点超压值最大,在t=7.8ms时达到最大值为18.65MPa,且在t=10ms左右时,位移和速度达到最大值,分别为10.47mm、3.303m/s;气体爆炸后,管廊燃气仓内墙正负压振荡时间持续较长,振动现象最为明显,最易发生破坏。     

刚性柱附近浅水爆炸荷载特性研究

刚性柱附近浅水爆炸时冲击波传播、气泡射流受多种因素影响。考虑水面、水底、刚性柱与水下爆炸冲击波及气泡的耦合作用,基于LS-DYNA有限元软件,建立浅水爆炸全耦合模型,通过经验公式验证有限元模型的正确性。研究表明:采用炸药直径1/3～1/2中心渐变网格能够较好地保证数值模拟精度。在冲击波传播阶段,刚性柱迎爆区冲击波峰值上升并产生切断现象,冲击波下降段被"截断",而背爆区冲击波峰值衰减约50%,同时正压作用时间增加;在气泡脉动阶段,气泡在收缩阶段产生指向刚性柱的气泡射流,当刚性柱与炸药之间的距离约为一个气泡半径时,刚性柱附近的脉冲荷载增幅最大,脉冲荷载最大测点水深较爆心上移。     

组合装药的撞击安全性与内爆威力试验研究

为了验证某组合装药的撞击安全性与内爆威力,开展了组合装药的大型落锤试验及爆炸罐内爆威力试验研究。结果表明:单一温压炸药在400 kg落锤作用下的临界落高为2.2 m,组合装药分别在2.2、2.3、2.5和2.7 m落高时均未发生爆炸反应,证明组合装药具有更好的撞击安全性;组合装药的初始冲击波超压峰值是单一温压炸药的61.9%,比冲量是单一温压炸药的99.4%,准静态压力峰值是单一温压炸药的94.5%。考虑到炸药在有限空间内爆炸的能量释放特性,将准静态压力峰值作为威力评价标准更为合适。试验结果证明组合装药的内爆威力与单一温压炸药相当。     

水下多点爆炸条件下的冲击波载荷特性

针对水下多爆源起爆的实战背景,开展了两点同时起爆条件下冲击波载荷特性的数值模拟研究。基于自研的多相可压缩流体计算程序,采用高精度的数值格式对流体控制方程进行离散求解。将数值模型计算的自由场水下爆炸的结果与理论结果比较,初步验证了数值模型计算的准确性与可靠性。利用该模型计算了典型工况下水下两点起爆工况,计算结果表明:两爆源对称面上压力相比单爆源线性叠加后的峰值压力增加12%~16%;两爆源垂直截面之间的压力存在双峰现象;而对于两垂直截面之外的测点压力也存在双峰现象,第1个峰值压力与单爆源线性叠加的峰值相等,第2个峰值压力要远低于单爆源线性叠加的峰值,峰值压力下降幅度可高达30%左右。研究结果能够为水下武器防护设计与威胁评估提供参考。     

水下爆炸特性的二维数值模拟研究

 对带有3 mm左右厚铝壳或钢壳的柱形TNT、RS-211、T/Γ（35/65）以及ROT901装药在水下采用一端点起爆后,水中冲击波超压分布进行了二维数值模拟研究。对该四种炸药的每一种,分别考虑了药量为1 kg、3 kg、5 kg和8 kg四种TNT当量的柱形装药情况。柱形装药的长径比为L/D=1.5左右。计算使用的程序为二维欧拉多流体网格法流体力学计算程序MFIC。给出了水中冲击波超压的空间和时间分布;对铝壳和钢壳两种情况的结果进行了分析、比较;同时给出了缩比模型与1∶1模型之间相应物理量之间的缩比关系。     

裸爆和特制半球形结构内爆超压对比实验研究

 爆炸焊接过程会产生爆炸地震、冲击波、有毒气体和噪音等有害效应。设计并制作一个直径为36 m、入口和排烟口直径分别为8 m的半球形结构,用以降低爆炸焊接过程中产生的超压影响。为了研究实际超压降低效果,进行了裸爆和按1/6比例建立的半球形结构内爆炸的超压对比实验。实验结果表明：当超压的传播距离大于20 m时,这种半球形结构能有效降低爆炸焊接过程中产生的超压;对于相同的超压传播距离,切线方向（垂直于半球形结构入口方向）比径向方向（半球形结构入口方向）的超压降低效果好。结合冲击波的传播和反射特点,对超压降低的可能原因进行了分析。     

9.5 m高水平极化有界波电磁脉冲模拟器内场分布特性的初步实验研究

测量了9.5 m高的水平极化有界波电磁脉冲模拟器的内部场,并根据实验测量结果分析了该模拟器内场分布特性,包括一定区域内场均匀性的定量分析及模拟器内部有效测试空间的确定方法,进而对最低位置为距离地面2 m的有效测试空间进行了预估。实验结果表明：位于该模拟器双锥中心正下方且距离该中心5.5～7.5 m的测点场的峰值基本按照测点与双锥中心间距的倒数衰减,且随着测点与双锥中心间距的增大,因锥与极板不连续结构导致的波形变化在时间轴上滞后,而因地面影响导致的波形变化在时间轴上提前;在距离地面比较高的水平面上,两极板之间场的外泄方向场的衰减比双锥中轴线方向场的衰减更慢;该模拟器内部距离地面2 m的水平面上12 m×12 m的区域内所取测点的归一化场平均峰值约为0.678,归一化场平均峰值的标准偏差约为0.068 9,场的均匀性约为2.039 dB。     

同心旋转圆柱间隙流场中纤维取向的动态模拟

本文对两同心旋转圆柱间隙形成的流场以及处于流场中的纤维运动和取向进行了数值研究. 在贴体坐标网格下求解了流场控制方程, 得到了流场中的速度、压力等物理量. 研究了两同心圆柱同速反向旋转以及仅内层圆柱旋转这两种情况下的纤维运动和取向状态. 得到了处于这两种情况下的纤维在流场中从静止到开始运动和取向直至最终达到稳定状态的动态详细过程. 结果表明, 当两个圆柱同速反向旋转时, 纤维运动与取向也相应的呈现两层结构; 而仅内圆柱旋转时, 纤维运动与取向呈单层结构. 在两种情况下, 纤维均沿流线方向运动和取向. 讨论了纤维长径比对纤维取向的影响, 结果表明随着纤维长径比的增加, 纤维沿流线取向的取向度逐渐增强.     

煤矿巷道瓦斯爆炸冲击波与高温气流的关系

 为了给瓦斯爆炸后煤尘二次爆炸的深入研究提供理论依据,应用计算流体力学方法,对煤矿巷道内瓦斯爆炸的瞬态流场进行了数值模拟,得到了冲击波与高温气流流动的时空关系,并借助实验对数值方法进行了验证。研究表明：在瓦斯爆炸后的一定时间内,近场区域和远场的部分区域极有可能引发煤尘二次爆炸。给出了可能发生煤尘二次爆炸的区域随瓦斯区长度的函数关系式,以及远场中峰值温度和峰值超压到达时间的间隔随轴向距离和瓦斯区长度的分布特性。     

泄压口大小对约束空间爆炸准静态超压载荷的影响规律

约束空间爆炸载荷主要包含瞬态冲击波和持续时间较长的准静态超压,为了研究泄压口大小对约束空间爆炸准静态超压载荷的影响规律,采用三阶WENO有限差分格式,编写了约束空间内炸药爆炸冲击波三维数值计算程序。利用Sod激波管、双爆轰波碰撞、激波与熵波相互作用等经典算例验证了数值程序的可靠性。基于验证的程序开展了带有泄压口的约束空间内爆炸载荷数值计算。在分析数值计算结果基础上,提出了描述准静态超压载荷简化理论公式,该理论公式与数值计算结果吻合较好。研究结果可为工程抗爆结构设计提供一定的参考和指导。     

新型结构炮口制退器的膛口冲击波数值研究与性能分析

面向炮口制退器综合性能设计需求,基于增材制造技术优势,提出了一种叠加冲击式内壁与反作用式外孔特征的新型小口径轻金属炮口制退器方案。基于三维无黏Euler方程,建立了后效期火药气体排空过程的有限元仿真模型,通过膛口流场仿真与流固耦合分析,研究了膛口流场发展、膛口冲击波与超压分布特征、制退器效率与其强度性能。结果表明:相比于传统结构,采用同种轻质钛合金材料的新型制退器具有较高的制退效率,对于机载平台具有较优的膛口超压分布特征和低冲击波危害效应,并且结构强度满足使用要求。     

深水爆炸冲击波作用下圆柱壳动态响应影响因素的数值模拟研究

 运用通用有限元程序ABAQUS,对圆柱壳在深水爆炸冲击波作用下的动态响应进行了数值模拟,研究了圆柱壳所处深度、爆心方位和预应力对于圆柱壳动态响应的影响。研究结果表明：在同样的水下爆炸冲击波作用下,随着深度的增大,圆柱壳的毁伤逐渐加重,而且响应过程中的变形逐渐增大,低频响应阶段的速度曲线趋于平坦;爆距不变时,不同的爆心方位在不同深度上的毁伤效果存在一定差异;预应力的存在使圆柱壳的毁伤加重,并且预应力对圆柱壳动态响应的影响与爆心方位有关。研究结果对于潜体结构的生命力评估具有一定参考价值。     

冲击波作用下引信传爆序列殉爆的数值模拟

为研究不敏感弹药或不敏感引信在战备和后勤贮存过程中的殉爆现象,利用装填JH-14C传爆药的某引信传爆序列,开展了冲击波作用下的殉爆数值模拟研究,获得了引信传爆序列爆轰波成长历程、传播规律以及临界殉爆距离,建立了冲击波能量判据,并给出了殉爆条件。结果表明:冲击波由传爆药左上向右下传播,在最右下端面处起爆,引信传爆序列的临界殉爆距离为9.7 mm;当作用冲击波能量大于临界起爆能量时,引信传爆序列发生殉爆。     

旋转圆柱目标的动态散斑统计特性(英)

高斯光束照射以固定角速度旋转的粗糙圆柱体表面,反射空间形成随着被测表面运动变化的动态散斑。研究了照射面积大于圆柱曲面条件下在远场衍射区形成的动态散斑统计特性,得到了不同入射条件及不同圆柱时动态散斑强度起伏空间-时间归一化相关函数。结果表明:随着圆柱半径、旋转速度以及照射光斑有效面积的增大,动态散斑强度起伏相关性快速减小;当圆柱旋转速度恒定,动态散斑光强起伏相关函数的峰值随着探测点之间距离的增大而减小,但峰值的位置随之而增大;在近似点照射情况下随着圆柱半径的增大空间相关长度基本不变,而相关时间性明显增大。     

旋转圆柱目标的动态散斑统计特性(英)

高斯光束照射以固定角速度旋转的粗糙圆柱体表面,反射空间形成随着被测表面运动变化的动态散斑。研究了照射面积大于圆柱曲面条件下在远场衍射区形成的动态散斑统计特性,得到了不同入射条件及不同圆柱时动态散斑强度起伏空间-时间归一化相关函数。结果表明：随着圆柱半径、旋转速度以及照射光斑有效面积的增大,动态散斑强度起伏相关性快速减小;当圆柱旋转速度恒定,动态散斑光强起伏相关函数的峰值随着探测点之间距离的增大而减小,但峰值的位置随之而增大;在近似点照射情况下随着圆柱半径的增大空间相关长度基本不变,而相关时间性明显增大。     

Statistical property of dynamic speckles produced by a rotating cylinder object

高斯光束照射以固定角速度旋转的粗糙圆柱体表面,反射空间形成随着被测表面运动变化的动态散斑.研究了照射面积大于圆柱曲面条件下在远场衍射区形成的动态散斑统计特性,得到了不同入射条件及不同圆柱时动态散斑强度起伏空间-时间归一化相关函数.结果表明:随着圆柱半径、旋转速度以及照射光斑有效面积的增大,动态散斑强度起伏相关性快速减小;当圆柱旋转速度恒定,动态散斑光强起伏相关函数的峰值随着探测点之间距离的增大而减小,但峰值的位置随之而增大;在近似点照射情况下随着圆柱半径的增大空间相关长度基本不变,而相关时间性明显增大.     

串联战斗部不同介质组合的隔爆能力

为了有效提升串联切割战斗部隔爆结构衰减爆炸冲击波的性能,解决前级聚能装药结构与后级随进弹的匹配及隔爆问题,在前级切割器和后级随进弹之间加装隔爆结构,使用有限元分析软件ANSYS/LS-DYNA建立模型,进行不同组合结构隔爆性能的数值模拟,比较隔爆能力。模拟结果表明:前级装药起爆后,爆炸冲击波首先向后级随进弹头靠里区域汇聚,而不是向弹头尖端区域汇聚,因此可以适当减薄外层金属隔爆介质头部尖端区域;将外层金属由硬质钢改成铝时,后端壳体应力峰值的变化很小,故确定外层金属介质为铝;铝-聚脲的隔爆能力优于铝-泡沫铝结构,最终确定"软"隔爆介质为聚脲。通过调整铝和聚脲层的厚度,确定了最佳隔爆参数,能够满足实际应用。     

小行星轨道演化研究中的数值方法

小行星是太阳系最重要的一类小天本,主要分布在两个区域;火星和木星轨道之间的一条主带和近地空间,近地小行星轨道的最大特点是其轨道半长径与地理轨道半长径相近,或近日距离接近甚至小于日地平均距离,其运动可深入到地球轨道的内部,迷将导致该类小行星与地球（还有金星、火星等）十分靠近甚至发生碰撞,这一特征给其轨道演化数值研究带来一些困难,包括天体力学方法中一般消除碰撞奇点的正规化处理以及对定性研究十分成功的辛算法都将在不同程度上失败,通过对几种常用数值方法（包括辛算法）计算效果的比较,根据小天体运动自身的特性,给出了相应处理措施,从而可提高计算结果的可靠性。