【作者】 文尚刚；

【导师】 孙承纬；

【作者基本信息】 中国工程物理研究院 ， 工程力学， 2001， 博士

【摘要】 本文首先回顾了爆轰传播理论发展历史，介绍了几种基本的理论模型及与爆轰波传播有关的其它内容。然后对等位线方法（Level Set，LS）作了一个简单的介绍，并根据该方法的原理，建立三维爆轰波传播的LS方程，研究该LS方程的数值求解方法，编制了三维爆轰波传播计算程序LS3D，提出了处理边界条件和爆轰波相互作用的方法，进行爆轰波相互作用多点起爆实验设计，完成各种起爆条件下爆轰波发展过程的数值模拟，并与实验结果作了比较。 本文提出了在三维直角坐标系中爆轰波传播LS方程离散化方法，建立了相应的差分方程，研究了对不同形状炸药在多种起爆形式（点、线、面、散心、聚心）下初值条件的设定，以及各种边界条件的处理方法。根据本文建立的计算编码LS3D，可以确定不同时刻爆轰波阵面的位置，再通过计算编码LS2DDRAW和LS3DDRAW画出某一截面的波阵面曲线图，或在整个计算空间内的三维波阵面曲面图。 本文的实验工作设计了两种多点起爆装置，采用了转镜式高速分幅照相和扫描照相技术及SVR相机高速摄影技术，观测了几种多点起爆装置中爆轰波在炸药柱（DG01A、B炸药和JB9014）底端面的发展过程及爆轰波相互作用后波形的演变情况。 本文最后根据所建立的三维爆轰波传播的LS方程的基本理论模型和该方程数值解法，编写了计算爆轰波传播的三维计算编码LS3D，并使用该编码对不同形状炸药在不同起爆方式（点起爆、线起爆、面起爆、发散爆轰、聚心爆轰等）下理想爆轰波的传播过程进行了计算，并且对钝感炸药JB9014几种多点起爆实验及圆柱形和圆弧片JB9014的平面起爆条件下爆轰波在炸药中的传播进行了数值模拟。并与实验结果进行比较，计算结果与实验结果符合得还比较好。另外，本文还给出了其它几个算例，主要是聚心爆轰、爆轰波在多介质炸药的平面波透镜中的传播及爆轰波的相互作用等。 本文计算结果表明，计算编码LS3D可以对不同起爆方式下不同形状炸药中三维爆轰波的传播进行较准确的计算，适于处理不同形状、多介质的炸药装置中 三维爆更波传播的匕方法研究爆轰传播问题，是工程计算的一种有力工具。但应指出的是，这类计算结果的精度在很大程度上取决于所采用的爆速曲率关系是否准确，以及所采用的爆轰波在炸药与空气或其它介质边界处的夹角是否合适等。这些关键参数的测量是近代爆轰传播实验研究的主要课题。 三维爆轰波传播的数值模拟，采用流体力学编码计算是十分困难且很不准确的。目前LS方法是一种较为现实和可靠的方法。本文的研究工作在这个重要方向上取得了较大进展，提出了该方法关键问题的一些初步的解决方法，明确了为发展和完善该方法应进一步探讨的难点。更多还原

【Abstract】 The development history of detonation wave propagation is reviewed first in this thesis. Some basic theory model and other contents about detonation wave propagation are introduced. Then the level set method is introduced. By means of this method, the properties of three dimensional detonation waves are studied. The main contents of this thesis includes: the LS equations of three dimensional detonation wave propagation; the numerical solution method of LS equations; the experimental design of multi-point initiation; the numerical simulation results of detonation wave propagating on different initiation conditions. In this paper, the numerical methods for the LS equation of three-dimensional detonation wave are introduced. Detail are introduced as following: transfering the LS equation to difference equation, setting the initial conditions for different shape explosive at different initiation way (point, line, plane, divergent or convergent); the dealing method of different edge conditions; how to draw the curve graph in two dimensional space or surface graph in three dimensional space with the calculated data. Two kinds of initiation way are designed in this paper. By means of high-speed frame camera, streak camera and SVR electrical camera, the propagation process and the interaction of the detonation wave are recorded for several kinds of explosive (DGO IA, B and JB90 14) at different imitation way. Finally, based on the basic theory model of LS method applied to three dimensional detonation wave propagation, we compile the calculation code LS3D. By employing the code LS3D, the propagation processes of ideal detonation wave at different shape explosives initiated by different way (point, line, plane) are calculated. The detonation propagation problems of insensitive high explosive JB9014 with different shape (cylinder, arc, pie) initiated from one plane are numerical simulated. The calculated results agree with the experimental results at some degree. In addition, the problems of convergent detonation and detonation wave propagating in planewave lens with multi-medium and the interaction of detonation waves are calculated. All the calculated results show that: the code LS3D can be used for calculating the detonation wave propagation in different explosive initiated by different way. It can deal with different shape, multi-medium explosive. The veracity of calculated results is depended on the relation of detonation speed with the curvature of detonation front and the edge angle of explosive with air or other medium. Measurement for these key parameters is an important study aspect of detonation wave propagation. It is very difficult and not veracious to simulate the three dimensional detonation wave propagation by means of fluid mechanics code. But the Level Set method is a practicable technology. In this paper, some primary methods of resolution for some key problems of LS method are given, and some propose are given for improving this method.更多还原