SPH算法基础上程序设计及其运用于成型装药战斗部之方法探讨.pdfVIP

SPH算法基础上程序设计及其运用于成型装药战斗部之方法探讨

第一章绪论

研究背景和意义

成型装药技术广泛地运用于弹药武器设计同工业生产,比如破甲弹设计,石

油开采,线形切割等。现如今，工程问题的计算机数值计算方法的出现给

理论和实验间的相结合起到了纽带的作用。计算机数值仿真能为纯理论给

予实际的验证，同时也能提供一个了解物理本质、解释物理现象和发现新

现象的方法。目前，所流行的数值计算方法依照是否依赖有限兀网格为基

础，主要分为两大类，一种是有网格方法，例如:有限差分方法，有限单

兀法(FEM)等;另一种就是无网格法，例如:粒子法等。有网格方法已在工

程科学的各个领域得到了广泛的应用，例如:计算流体力学，计算固体力

学等，并在这些方面取得了相当大的成功。

尽管目前有网格方法已经在国防技术领域以及民用技术中取得了显著的成

果，但并不代表它在各个方面都是无可挑剔的。有网格法在很多的工程应

用上还存在问题，并目_在很大程度上这些问题限制了它在工程应用中的

进一步发展。综合分析这些困难，有网格方法存在困难的根源是以网格为

背景，因为用这种方法来解决问题的时候必须要尽量得到高质量的网格，

而有些问题得到高质量的网格是很难做到的「0，例如:在自由面问题、不

规则边界问题、移动接触面问题、大变形问题上都无法得到一个高质量的

网格，这样就不能保证精度较高的结果。因此，需要采用其他的数值计算

方法来取代有网格方法，这就是无网格方法，而无网格法中最为典型的方

法就是基于核函数的光滑粒子流体动力学方法(简称SPH方法)。

SPH方法是一种粒子方法，核心思想是插值理论，类似于质点网格方法CPi

C-Particlein

Cell，但根本不同点在于SPH方法中计算空间导数时不需要任何网格，即

通过一个称为“核函数”的积分核进行核函数估值近似，将流体力学方程

或固体控制方程组转化为数值计算用的SPH方程组。在整个物理场中，物

质被离散为若干有限个“粒子”，这些粒子身上负载着所有的物理场量，

并能按照流体规律任意流动，这样就不需要考虑网格的缠绕和扭曲变形等

问题，因此该算法应用限制较少，理论上可用于任何变形问题。

本课题通过对成型装药战一斗部基于SPH算法的仿真软件的开发和应用，

为科研工作者提供了一条方便快捷的通道。避免了因网格引起的建模困难

、计算过程易出错等难以解决的问题，并目_能得到模型每时每刻各个点

的物理量，为更深入详细的研究提供了有力条件，同时该软件为大型通用

软件使用不熟练的工程设计人员提供了快捷方便的途径。

2成型装药战斗部数值模拟的研究现状

成型装药战一斗部在国防建设中具有重要的地位，是当代战争必不可少的

军事装备，主要包括聚能射流((ShapedChargeJet

SC)战一斗部、爆炸成型弹丸(ExplosivelyFormedProjectile

EFP)战一斗部及杆式侵彻体(JettingProjectile

Charge)二大类。它们都具有成型复杂、初速大、高动态等特点。较其他

常规战一斗部，虽然在应用方面较为突出，但同时在设计制造方面也更为

困难。3.2.1函数及其导数的积分26-28

3.2.3粒子近似法28-29

3.3射流形成过程SPH相关问题29-39

3.3.1核函数的构造29-31

3.3.2相互作用粒子对31-33

3.3.3其它重要问题及33-39

第四章基于SPH算法的成型装药39-58

4.1程序设计思想39

4.2自编制SPH程序结构39-42

4.3自编制SPH程序主要模块42-51

4.3.1前后处理接口模块42-44

4.3.2初始化参数模块44

4.3.3核函数计算模块44-45

4.3.4粒子对搜索模块45-46

4.3.5光滑长度更新模块46-47

4.3.6人工粘度模块47-48

4.3.7内力计算模块48

4.3.8