多组分PBX炸药细观结构冲击点火数值模拟

建立了炸药颗粒自由堆积三维计算模型,考虑了炸药颗粒尺寸和位置的随机分布、粘结剂对炸药颗粒的包覆和不同组分炸药颗粒间的级配;采用非线性有限元方法,对炸药颗粒由自由堆积到密实装药的压药过程进行模拟,构建了PBX炸药细观结构模型;对多元PBX炸药（HMX+TATB+Estane）细观结构冲击点火过程进行了计算,考虑冲击作用下炸药内部热力耦合作用和自热反应,计算不同组分比例混合炸药的冲击点火性能,分析了炸药组分对冲击点火的影响。研究表明TATB含量增加,混合炸药冲击感度降低。

     

HMX颗粒炸药低速撞击点火实验研究

改造落锤撞击仪器,在落锤与撞击底座中增加光学系统和压力测试部件,利用高速摄像仪,拍摄了 HMX颗粒炸药受撞击加载过程中微秒量级时间的变形-熔化-点火-燃烧过程图像。研究了下落高度、颗粒松 散程度对变形和点火时间长短及燃烧剧烈程度的影响。结果表明,落锤下落高度越高,塑性扩展持续时间越 短,易发生点火而不易发生喷射现象;较低的落锤下落高度易导致剧烈的喷射,而不易发生点火。颗粒排列越 松散的样品,塑性扩展时间越短,更容易形成多个点火位置,从而使热点聚合导致燃烧的机会越大。     

射弹倾斜撞击带盖板炸药引发爆轰的条件

用二级轻气炮发射圆柱形、球形钢射弹以不同的角度撞击带不同厚度钢盖板的 ＴＮＴ／ＲＤＸ（４０／６０）炸药,得到了不同条件下引发炸药爆轰的阈值射弹速度。可以用ｖｄ~　（１／２）＝（１ ＋ｋ）［Ａ＋Ｂｈ／（ｄｃｏｓ）］描述临界引爆条件。对于带钢盖板的ＴＮＴ／ＲＤＸ（４０／６０）炸药,Ａ＝ ３．３３．Ｂ＝５．３４;圆柱形平头射弹撞击,ｋ／７５,球形射弹撞击,ｋ０．５＋０．２（１／ｃｏｓ－１）。由 此,本研究将Ｊａｃｏｂｓ 引爆判据推广到了斜碰撞条件。     

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

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

PBX9502炸药爆轰产物的状态方程

采用vanderWaals等效单组分流体模型和Ross硬球微扰理论软球修正模型,计算了爆轰气相产 物的状态方程;用石墨相、金刚石相、类石墨液相和类金刚石液相4种相态描述凝聚成分,由Gibbs自由能最 小确定了不同状态下的凝聚产物相态。对爆轰产物混合系统采用自由能最小原理,通过化学平衡方程组求解 了炸药爆轰产物系统的平衡组分。使用该理论计算了高含碳炸药PBX9502Chapman-Jouguet(CJ)点的爆轰 参数,计算值与实验值符合很好;同时计算了3条等温线,并与Sesame库比较,发现温度超过1000 ℃时,计 算值与Sesame库的计算结果比较接近。在计算的5802K 等温线上发现了一个拐点,分析发现是由于在此 处游离态的碳发生了相变。     

PBX炸药的拉伸断裂损伤分析

采用巴西实验作为间接拉伸加载手段,研究了某PBX(polymer-bonded explosive)炸药试样拉伸作用下的断裂损伤特性.发展了PBX炸药的光学制备技术,获得了试样在光学显微镜下的细观表面形貌和断裂损伤形貌,结合高速摄影和数字相关分析技术获得了试样的形变和破坏过程.实验结果表明该PBX炸药的静态拉伸强度低于...     

PBX 炸药本构关系的实验研究

在SHPB装置上开展了某典型PBX炸药的单轴压缩、间接拉伸实验。利用入射波整形技术,结合半导体应变片和石英晶体联合测试的方法,实现了PBX炸药的应力平衡和恒应变率加载,得到了不同应变率下某PBX炸药的压缩应力应变曲线,初步建立了该炸药的修正Sargin唯象模型;开展了3种PBX炸药的动态巴西实验,通过高速摄影获得了炸药的破坏过程,结合数字相关技术,获得了试样裂纹附近的应变场分布,初步建立了描述3种炸药动态拉伸行为的修正Johnson-Cook模型。模型曲线与实验结果符合较好。     

孔隙度对PBX炸药冲击起爆影响的实验研究

为了研究孔隙度(装药密度)对PBX炸药冲击起爆爆轰成长的影响,采用炸药冲击起爆锰铜压阻一维拉格朗日实验测试系统,测量了不同孔隙度的PBXC03炸药(HMX的质量分数为87%,TATB的质量分数为7%,黏结剂的质量分数为6%)冲击起爆过程不同拉格朗日位置的压力-时间历史。结果显示:在本文装药范围和加载条件下,孔隙度对PBX炸药冲击起爆爆轰过程的影响不单调,中等密度的炸药冲击起爆和爆轰成长最快,这是热点点火过程与燃烧反应过程共同作用的结果。     

含能破片引爆带壳炸药过程的数值模拟

将引爆模拟战斗部等效为带铝壳炸药,设计了一种新型含能破片作为毁伤元,利用非线性有限元LS -DYNA软件对该含能破片侵彻、引爆带壳炸药的作用过程进行了数值模拟。用“升降法”得到了该“含能破片”对不同盖板厚度带壳炸药的引爆速度,同时与普通破片引爆同规格带壳炸药进行了对比,并进行了实验验证。结果表明,通过控制含能破片的撞击速度和含能物质的延迟起爆时间,可有效引爆盖板厚度为8~16mm的带铝壳炸药。

     

垫层对破片冲击起爆带壳炸药影响的数值模拟

将战斗部简化为带钢壳和垫层的炸药组件,用3种不同厚度的材料作垫层,利用非线性有限元软件AUTODYN对高速破片侵彻、引爆带壳炸药的作用过程进行了数值模拟,得到了高速破片冲击起爆带壳带垫层炸药的波后压力剖面。对不同厚度多种材料作为垫层时带壳炸药的冲击到爆轰距离进行了对比分析,结果表明,通过控制垫层的材料和厚度,可以对破片冲击起爆带壳炸药进行有效抑制。研究结果对战斗部破片起爆的防护设计具有指导意义。

     

冲击起爆弹丸内装药延迟起爆数值模拟

对穿爆燃弹的穿靶及释能过程进行数值模拟,解释了SPH-FEM耦合方法在模拟冲击起爆过程的优越性。通过对不同弹芯头部形状、直径及材料的对比分析,得到了不同工况下装药热点成长的压力-时间曲线,计算结果表明:在弹芯直径不变的情况下,弹芯头部直角尖刺越短,装药的热点成长时间越短;在头部尖刺长度不变的情况下,减小弹芯直径,热点生成时间也缩短;选用钢弹芯比钨合金弹芯有一定靶后释能优势。模拟穿靶效果与真实穿靶效果符合较好,这种方法可以为穿爆类弹丸设计提供依据。     

复合材料T型接头冰雹高速撞击损伤的数值模拟

利用高速空气炮进行了冰雹撞击复合材料T型接头的实验,研究了不同撞击速度下结构的损伤情况。同时,采用光滑质点流体动力学方法与黏聚区模型相结合的方法,对冰雹撞击复合材料T型接头进行了数值模拟。通过与实验结果的对比,验证了数值模拟模型的有效性。在此基础上,利用数值模型研究了影响复合材料T型接头冰雹撞击损伤的各种因素。计算结果表明：冰雹撞击对复合材料T型接头造成的损伤主要为分层损伤;冰雹的撞击速度、尺寸、入射角等,都对T型接头的损伤程度有很大影响;T型接头的分层面积与冰雹的撞击能量之间呈近似线性关系,分层面积随着撞击能量的增大而增大;冰雹的入射角越大,分层面积与撞击力峰值也越大。

     

高温高压气体在低含水率介质中迁移的数值模拟

针对气液两相非等温渗流模型高度非线性的特点,发展了适宜的数值离散方法。根据相态转换准则和控制方程的性质,采用最低饱和度法简化算法。空间离散方面,使用有限体积法;时间离散方面,设计了一套包含合理求解顺序的Picard迭代法,解决了方程组强耦合的问题。利用上述数值方法对高温高压气体的迁移行为进行数值模拟,证明了气体在低含水率介质和等效孔隙度的干燥介质内的运动基本一致,并分析了空腔内的气液相态转变过程。在此基础上,研究了多孔介质孔隙度和渗透率对气体压强演化和示踪气体迁移的影响。研究表明,孔隙度越小（相同渗透率）、渗透率越高（相同孔隙度）,示踪气体的迁移距离越远,并给出了估算不同孔隙度和渗透率下迁移距离的半经验公式。     

冲击作用下易碎式密封盖数值仿真及实验分析

确定了易碎式密封盖的材料和结构,运用有限元法进行了火箭弹和燃气流冲击的数值仿真计算,通过冲破实验进行了间接验证,从而掌握了密封盖的破碎效果和性能。结果表明,密封盖能沿预置沟槽方向破裂,能满足引信体的安全性和可靠性,不影响火箭弹的正常飞行;燃气流冲击力对密封盖影响不大,不会对密封盖造成破坏。     

Numerical simulation of Taylor impact tests

The present paper discusses some aspects in the numerical simulation of Taylor impact tests. A phenomenological internal variable theory is presented and restricted by the second law of thermodynamics. The constitutive model consists of a thermo-hyperelastic equation for the stresses and an evolution equation for the plastic internal variable. Specific thermo-plastic constitutive functions are proposed and corresponding material parameters are identified. Taylor impact tests are numerically simulated using the explicit finite element program LS-DYNA augmented by an user-defined material subroutine and the effect of variation of selected model parameters is discussed. Numerical results allow new interpretations of experimental observations and test data and gives advice on identification of material parameters in rate-dependent inelastic constitutive models.     

高应变率下断裂韧性实验的数值模拟

采用有限元软件ANSYS/LS-DYNA程序对静态和冲击荷载作用下的含裂纹半圆弯曲(SCB)实验进行了数值模拟。根据静态实验的模拟结果,提出了适合复合型加载的Ⅰ型应力强度因子拟合公式,采用该公式计算应力强度因子的最大误差不超过10%。动态实验的模拟结果表明:对于纯Ⅰ型加载的SCB实验,动态应力强度因子随着试样半径、支座间距以及相对裂纹长度的变化呈现规律性变化;当试样半径小于60mm、相对支座间距为1.2、相对裂纹长度在0.1~0.4范围内时,惯性效应的影响较小,采用静态拟合公式计算裂尖的动态应力强度因子的误差约10%;对于复合型加载的SCB实验,当相对裂纹长度为0.2~0.4、裂纹倾角在10°~40°范围内时,采用静态拟合公式计算裂尖的动态应力强度因子的误差小于10%。     

充满液体弹塑性圆柱壳轴向冲击屈曲的数值模拟

基于对充满液体的弹塑性圆柱壳轴向冲击屈曲的实验观察 ,建立了相应的计算模型 ,借用 LS-DYNA软件包 ,在 HP-C36 0工作站上完成了数值模拟 ,得到了动画显示 ,给出的屈曲模态及一些特征量的时程曲线与实验结果相当一致。文中对壳的几何、物理参数对屈曲性能的影响进行了较详细的讨论     

DAMAGE PROGRESSIVE MODEL OF COMPRESSION OF COMPOSITE LAMINATES AFTER LOW VELOCITY IMPACT

IntroductionCompressivepropertiesofcompositelaminatesafterlowvelocityimpactareoneofthemostseriouscircumstanceswhichmustbetakenintoaccountindamagetolerancedesignofcompositestructures[1].Impactdamagegenerallymanifestsintheformsofdelaminations,matrixcracksandfibrefracture[2 ,3].Sofaralotofresearcheshavebeentakeninlowvelocityimpactbehaviorofcompositelaminatesandtheirpost_impactcompressiveproperties.Compressivefailuremechanismsofcompositelaminatesafterlowvelocityimpacthavenotbeenknownclearly .Ther…     

Numerical modeling of ice behavior under high velocity impacts

In this work a constitutive relation for ice at high strain rates and an algorithm for its numerical integration are developed. This model is based on the Drucker–Prager plasticity criteria, which allows a different behavior in tension and in compression. In addition a failure criteria, based on pressure cut-offs, is implemented to describe the ice damage. In order to validate the constitutive model, numerical simulations were compared with experimental results, in which ice cylinders were impacted against a steel plate, allowing the measurement of the contact load. Three different numerical solvers are used in order to analyze its performance to appropriately modeling the ice behavior.     

镁单晶低温压缩形变机制的原子模拟

应用分子动力学方法模拟镁单晶在单向压缩作用下的力学性能和微观结构演化过程.在镁单晶低温(5K)压缩的模拟结果中可以同时观察π2面上两种弗兰克柏氏矢量的滑移,即1/6[0223]和1/6[2203].两种滑移的同时启动使得晶体在完成滑移后又恢复到完整的密排六方(hcp)晶体结构,只在滑移区的端部形成内部缺陷和表面形成滑移台阶.