水下爆炸过程的高精度数值计算

采用三阶精度的Parabolic Piecewise Method(PPM)方法和Volume Of Fluid(VOF)方法相结合，运用Lagrange-Remapping算法，数值计算爆轰产物、水、空气多种介质之间的相互作用过程，该方法可以用来处理界面两边高密度比可压缩流动以及强剪切滑移运动等问题。     

炸药爆轰驱动水的初期过程

采用高速转镜分幅相机,利用阴影照相技术研究了PBX-01炸药爆轰产物在水中的初期演变过程, 获得了爆轰产物与水界面在不同时刻的图像。利用LS-DYNA程序对该实验进行了数值模拟,得到了水中冲 击波迹线、炸药爆轰产物与水界面随时间膨胀的位置、速度。数值计算结果与实验结果具有较好的一致性,所 得参数可为进一步设计改进实验提供参考。理论分析表明,水在爆轰产物作用下的最高瞬间温度约342K, 未到达水在实验压力条件下的相变汽化温度。     

水中钢板爆破水介质对装药量的影响

为分析水介质对水中钢板爆破装药量的影响,理论推导了钢板背衬水介质条件下,钢板爆破最小装药量与空气中钢板爆破最小装药量的倍数关系:最小倍数关系为3.76。数值计算了钢板背衬空气介质和背衬水介质情形下钢板爆破的最小装药量,其倍数关系在3.5左右,与理论结果相近,表明钢板背衬水介质或空气介质是决定装药量大小的关键因素。数值计算了装药在空气介质中及水介质中钢板爆破爆轰产物对钢板的冲量大小,结果接近,表明水介质对炸药爆轰产物的约束作用是影响水中钢板装药量的次要因素。     

空气自由场中强爆炸冲击波传播二维数值模拟

编写了适用于模拟具有高密度比、高压力比的强激波问题的二维柱对称多介质流体计算程序。利用有限体积方法求解流体的Euler方程组,采用level set方法捕捉爆炸产物与空气的运动界面,并通过求解物质界面两侧Riemann问题的精确解来计算爆炸产物与空气之间的数值通量。研制了三角形网格自适应技术来实现网格的自动加密和粗化,在保证捕捉激波峰值的前提下有效地提高了计算效率。利用计算程序对1 kt TNT当量的空气自由场强爆炸问题进行数值模拟,计算得到的峰值超压、冲击波到达时间等物理参数与点爆炸理论结果基本一致。     

PBX-01炸药水中爆轰产物状态方程研究

提出通过水中实验确定炸药的水中爆轰产物JWL状态方程参数的方法;选择PBX-01高能炸药进行水中实验,利用ANSYS/LS-DYNA程序建立炸药的水中实验模型,将实验结果与数值计算结果进行对比,确定PBX-01炸药水中爆轰产物的JWL状态方程参数。研究结果显示,圆筒实验确定的JWL参数在反映炸药水中爆轰产物的膨胀状态时有所不足,水中实验确定的JWL状态方程参数能够更准确地描述PBX-01炸药水中爆轰产物的膨胀过程,因此对水中爆炸的研究需要通过水中爆炸实验建一套状态方程参数。     

带隔板装药爆轰波马赫反射理论研究和数值模拟

基于三波理论和Whitham方法对带隔板装药爆轰波相互作用后发生的正规反射和非正规反射进行了理论分析,给出了爆轰波发生马赫反射时临界入射角和马赫杆增长角等参数的变化规律,提出了马赫杆高度的计算模型。基于凝聚炸药爆轰Jones-Wilkins-Lee(JWL)模型和冲击起爆的Lee-Tarver模型,利用有限元计算软件对带隔板装药爆轰波的传播过程进行了数值模拟。结果表明,发生马赫反射后,随着爆轰波的传播,马赫杆的高度不断增加。数值模拟结果与理论计算结果吻合较好,说明本文中采用的理论模型和数值模拟方法能够较准确地描述带隔板装药爆轰波马赫反射的传播过程。     

大体积混凝土水下接触爆炸破坏分区特征分析

为探究大体积混凝土水下接触爆炸破坏分区特征,基于水中爆炸冲击波与混凝土的相互作用过程分析,建立了综合考虑爆炸冲击波冲击破碎和爆轰产物准静态压拉致裂机制的混凝土爆炸破坏分区计算方法,并与有限元数值模拟和试验实测数据开展对比验证。结果表明：与空中接触爆炸相比,水对爆轰产物膨胀起抑制作用,使得爆炸荷载持时增加、作用于周围介质的冲量增大;采用建议的环向压碎判据计算破碎区,并将开裂区分为动态压裂、准静态压裂和准静态拉裂区的计算方法能够很好地预测混凝土水下接触爆炸破坏分区范围;炸药类型和起爆水深一定时,混凝土的抗拉强度和抗压强度比对开裂区范围起重要影响。     

爆轰波平掠惰性气体界面及其解耦现象的数值研究

采用基元化学反应模型和迎风TVD格式,数值研究了爆轰波平掠惰性气体界面时的物理现象及其作用机制,并用点隐方法克服化学反应源项引起的计算刚性。数值结果显示,当爆轰波平掠过惰性气体界面时,形成了爆轰波、界面、透射激波以及稀疏波相互作用的现象。在高N2比例稀释的可燃混合气体情况下,当爆轰波平掠过特定惰性气体界面时,它与惰性气体界面相互作用产生的稀疏波可以导致爆轰波的解耦。     

水及盛水容器对近距离爆炸载荷影响的实验研究

为研究用水包围炸药的方式对爆炸载荷的作用机理,在两端开口的钢筒内进行了水直接包覆炸药的爆炸实验,利用光纤位移干涉仪获取了钢筒外壁的径向速度和动态变形。结果表明:不同于无水爆炸,爆轰产物通过水的“裂缝”在空气中形成冲击波,造成该冲击波出现时间更晚、强度更低、持续时间更长,并要求相应的数值模拟采用二维以上的计算模型;盛水结构的材料密度越低、厚度越小对爆炸载荷的影响越小。     

射孔弹爆轰波间干扰的有限元分析

通过有限元显式动力分析程序的显式算法,对石油射孔弹顺序起爆后弹间爆轰冲击波干扰的过程进行了2维数值模拟,并考虑了弹体对弹间干扰的影响,对模拟结果进行了讨论;分析了影响高孔密射孔技术的弹间干扰产生的原因,主要是爆轰冲击波的压力场受到不同程度的干扰,致使初始爆轰波形破坏.计算结果与已有的物理现象和规律吻合,说明该计算模型和...     

A NEW EFFICIENT METHOD TO BOUNDARY VALUE PROBLEM FOR BALLISTIC ROCKET GUIDANCE

IntroductionThe problem of motion control for the mass center of ballistic rocket pertained typicallyto the boundary value problem for two points with the non-linear and multi-variables,namely,given the initial states of rocket,in order to satisfy the res…     

全机绕流Euler方程多重网格分区计算方法

全机三维复杂形状绕流数值求解只能采用分区求解的方法,本文采用可压缩Euler方程有限体积方法以及多重网格分区方法对流场进行分区计算。数值方法采用改进的van Leer迎风型矢通量分裂格式和MUSCL方法,基于有限体积方法和迎风型矢通量分裂方法,建立一套处理子区域内分界面的耦合条件。各个子区域之间采用显式耦合条件,区域内部采用隐式格式和局部时间步长等,以加快收敛速度。计算结果飞机表面压力分布等气动力特性与实验值进行了比较,二者基本吻合。计算结果表明采用分析“V”型多重网格方法,能提高计算效率,加快收敛速度达到接近一个量级。根据全机数值计算结果和可视化结果讨论了流场背风区域旋涡的形成过程。     

多体系统动力学求解与动画输出的并行处理

在介绍了多体系统动力学的建模方法及数值仿真的原理之后，着重介绍如何利用局域网的技术，在多台计算机之间实现数值求解和动画计算与输出的并行处理。用一个柔性索动力学仿真算例，论证了该技术的可行性，为多体系统动力学仿真的并行处理打下了基础。     

Calculation of coupled problem between temperature and phase transformation during gas quenching in high pressure

The gas quenching is a modern, effective processing technology. On the basis of nonlinear surface heat-transfer coefficient obtained by Cheng during the gas quenching, the coupled problem between temperature and phase transformation during gas quenching in high pressure was simulated by means of finite element method. In the numerical calculation, the thermal physical properties were treated as the functions of temperature and the volume fraction of phase constituents. In order to avoid effectual "oscillation" of the numerical solutions under smaller time step, the Norsette rational approximate method was used.     

爆轰产物压力势能缓释方法研究

为扩大脉冲爆轰技术的应用领域 ,抑制其振动、噪声等不利因素 ,提高能源转换效率的可能性 ,本文通过实验和数值计算 ,研究了单次爆轰产物通过孔盘将压力势能缓慢释放到相邻腔体的现象。实验在长 2m、内径 40mm的爆轰管中进行 ,数值模拟采用自适应的有限体积方法和有限速率化学反应模型。实验和数值模拟分析的结果均表明 ,在爆轰管与相邻腔体之间放置孔径适当的隔离孔盘 ,可以使爆轰产物的压力势能逐渐地传递到腔体内 ,造成随时间逐渐上升的压力分布。     

弹性接触问题的常刚度迭代方法及其加速

通过引入接触间隙荷载，在荷载项中反映接触状态变化的影响，避免了重新形成接触刚度矩阵，构造了弹性接触问题的常刚度迭代格式。然后借助于不动点迭代的Steffensen加速技术给出了一种求解弹性接触问题的常刚度迭代加速算法。文中最后给出两个算例表明了本文方法的有效性。     

基于单位分解法的无网格数值流形方法

在数值流形方法和单位分解法的基础上，提出了无网格数值流形方法. 无网格数值流形 方法在分析时采用了双重覆盖系统，即数学覆盖和物理覆盖. 数学覆盖提供的节点形成求解 域的有限覆盖和单位分解函数；而物理覆盖描述问题的几何区域及其域内不连续性. 与原有 的数值流形方法相比，无网格数值流形方法的数学覆盖形状更加灵活，可以用一系列节点的 影响域来建立数学覆盖和单位分解函数，具有无网格方法的特性，从而摆脱了传统的数值流 形方法中网格所带来的困难. 与无网格方法相比，由于采用了有限覆盖技术，试函数的构造 不受域内不连续的影响，克服了原有的无网格方法在处理不连续问题时所遇到的困难. 详细推导了无网格数值流形方法的试函数和求解方程，最后给出了算例，验证了该方法的正 确性.     

长杆弹对混凝土的侵爆效应

对适用于长杆弹在混凝土介质中侵彻和爆炸全过程的三维数值模拟方法和技术进行了研究。描述了进行相应数值模拟的有关方法和关键技术。确定了靶体C30混凝土材料所使用的本构模型及其相应的参数。对卵形头长杆弹在C30混凝土中侵彻到一定深度再爆炸的全物理过程进行了三维数值模拟,分别给出了C30混凝土靶体在侵彻和爆炸作用下的破坏效应图像。将侵彻计算图像与实验结果进行了比较,两者定性符合。     

凸轮机构系统磨损及可靠性寿命的数值仿真

基于离散数学理论和计算机技术提出了平底从动件凸轮机构系统磨损及可靠性寿命的数值仿真模型，引入了位置矢量和磨损步长的概念，使连续的磨损过程离散化，对凸轮的磨损寿命及其可靠度进行了预测，并对从动件的运动规律进行了动态跟踪．结果表明，复杂的磨损过程可用数值法进行模拟仿真，从而摆脱了传统的仅依赖实验的相关研究方法，通过一系列离散性准静态模型解决了经典微积分数学方法无法解决的动态和非线性磨损问题．所建立的方法具有良好的工程应用前景．     

PROGRESS IN MEASUREMENT TECHNOLOGY FOR DRILL STRING VIBRATION

随着深井、超深井数量的快速增加, 钻柱振动导致的钻具失效问题更加突出. 针对钻柱振动的研究可分为理论与数值模拟、测量分析两种方法. 由于钻柱振动具有复杂的非线性特征, 使得理论与数值模拟研究受到了很大限制, 因此井下振动测量技术的研究显得尤为重要. 本文比较全面和系统地介绍了国内外钻柱振动测量技术研究现状和进展情况, 并对几个应用较为成熟的国外测量系统的工作原理、分析方法及应用技术进行了详细的综述. 所得结论可为我国井下振动测量技术的发展提供重要参考.