非等比例缩比侵彻/贯穿相似规律研究

刘源; 皮爱国; 杨荷; 冯吉奎; 黄风雷

doi:10.11883/bzycj-2019-0086

非等比例缩比侵彻/贯穿相似规律研究

doi: 10.11883/bzycj-2019-0086

刘源^1,,
皮爱国^1, ,,
杨荷²,
冯吉奎¹,
黄风雷^1, ,

1.
北京理工大学爆炸科学与技术国家重点实验室，北京 100081
2.
中国工程物理研究院电子工程研究所，四川绵阳 621999

详细信息

作者简介:
刘　源（1994－　），男，硕士研究生，lymoxing@163.com

通讯作者:
皮爱国（1977－　），男，博士，副教授，aiguo_pi@bit.edu.cn

黄风雷（1965－　），男，博士，教授，huangfl@bit.edu.cn

中图分类号: O385
计量
- 文章访问数: 5168
- HTML全文浏览量: 1628
- PDF下载量: 112
- 被引次数: 10
出版历程
- 收稿日期: 2019-03-25
- 修回日期: 2019-04-27
- 刊出日期: 2020-03-01

Study on similarity law of non-proportionally scaled penetration/perforation test

LIU Yuan^1
,,
PI Aiguo^{1
, ,},
YANG He²,
FENG Jikui¹,
HUANG Fenglei^{1
, ,}

1.
State Key Laboratory of Explosive Science and Technology, Beijing Institute of Technology, Beijing 100081, China
2.
Institute of Electronic Engineering, China Academy of Engineering Physics, Mianyang 621999, Sichuan, China

摘要

摘要: 高速侵彻弹体的弹载部件/关键元器件的生存性与可靠性考核是引战系统研制领域的热点与难点问题，受原型试验的成本限制，利用缩比弹体搭载原型引信部件开展非等比例缩比试验研究是可行途径。针对传统等比例缩比方案无法满足弹体刚体过载相似性要求的情况，研究了非等比例缩比侵彻/贯穿相似规律，提出了非等比例缩比侵彻试验设计方法。数值计算结果表明：侵彻半无限厚混凝土靶条件下，非等比例缩比弹刚体过载的脉宽、幅值均可实现与原型弹刚体过载一致的加载条件；贯穿多层薄靶的条件下，通过调节靶板布置及弹体初速等试验工况，合理设计缩比弹体结构，可使非等比例缩比试验的弹体刚体过载峰值和脉宽覆盖原型试验。通过缩比模型试验得到的刚体过载特性可以为弹体及引信部件抗过载防护设计提供可靠的参考依据。
- 非等比例缩比 /
- 侵彻/贯穿 /
- 刚体过载 /
- 弹载引信 /
- 数值计算
Abstract: Survivability and reliability assessment of components/key components on high-speed penetrating projectiles is a hot and difficult issue in the field of EPW development. Due to the cost limitation of prototype test, it is feasible to carry out non-proportionally scale experimental research by carrying prototype fuze components on scaled projectiles. Through the analysis of the process mechanism of a projectile penetrating concrete target, the analytic solution of rigid-body deceleration when the projectile penetrating the semi-infinite thick concrete target and the multi-layer thin concrete target are discussed respectively. From the point of view of similarity of rigid-body deceleration, the non-proportionally reduced-scale criterion of projectile is proposed when the traditional scaling scheme can not meet the requirements of similarity. The numerical results show that under the condition of penetrating semi-infinite thick concrete target, the rigid-body deceleration of the non-proportionally reduced-scale projectile can achieve the same conditions as which of the prototype projectile from the point of view of pulse width and amplitude; under the condition of penetrating multi-layered thin target, through reasonably setting the scale factor and adjusting the layout of the target plate and the initial velocity of the projectile. The pulse width and amplitude of the rigid-body deceleration in the reduced scale test can cover them in the prototype test. The rigid body deceleration characteristics obtained from scaled model test can provide reliable overload environment reference for missile projectile design.
- non-proportionally scaled /
- penetration/perforation /
- rigid-body deceleration /
- rocket-borne detonator /
- numerical simulation

HTML全文

在包含运动边界的非定常流动模拟中，流场物理空间随时间变化，边界运动导致网格也随时间变化，典型应用包括多体分离、机翼抖振颤振、舱门启闭等^[1]。CFD计算一般采用运动嵌套网格、笛卡尔自适应网格和变形动网格等技术处理网格变化。在采用上述技术处理网格变化时，均涉及流场信息的插值问题^[1-2]，即利用旧网格上的流动信息插值获得重构区域流动参数。此外，在爆轰流场模拟中，爆轰计算网格与冲击波传播计算网格差两个量级^[3-4]，为提高计算效率，往往需要将物理量由爆轰计算时使用的密网格插值传递到冲击波传播模拟中的疏网格。

插值过程一般应具备守恒性和单调性等基本特性。对爆轰计算等守恒性要求高的问题，一般的插值不能保证计算区域内物理量的守恒，会降低计算精度，且无法捕捉激波的正确位置^[5]，需要发展守恒插值方法。守恒插值方法主要可分为两类，基于面通量(SFB/DC，simplified face-based donor-cell)的方法和基于单元相交(CIB/DC, cell-intersection-based donor-cell)的方法^[6]。

SFB/DC方法^[6]效率较高，但要求新、旧网格具有相同的拓扑结构，不适用于网格重构等拓扑发生变化的情况。CIB/DC方法思想简单直接，对计算网格的类型、结构不进行任何假设，适用范围很广，并且具有保号性，插值过程中不会产生新的极值^[6]。但该方法推广至三维情况时，由于三维网格单元重叠区域切割计算非常复杂，给这种方法的应用带来了很大困难。为了实现严格意义下的三维CIB/DC算法，P.E.Farrell等^[7]进行了一系列卓有成效的工作，在2011年进一步发展为“局部超网格”(Local Supermesh)方法，成功实现了非结构网格下的三维CIB/DC算法^[8]。

本文中在P.E.Farrell等和S.Menon等^[9]工作的基础上，构造了应用于二维/三维混合网格的CIB/DC算法，实现了基于精确网格切割的积分守恒插值。验证结果表明，该方法能够保证插值过程中计算域内物理量的严格守恒，具有比常规二阶插值更高的精度。

1. 守恒插值基本思想

守恒插值要求在插值过程中保证守恒量的积分值保持不变，假设计算域为Ω，网格T^a是计算域的一种划分，则在网格T^a中，Ω被划分为N_a个网格单元T^a=(T₁^a, T₁^a, …, T_{N_a}^a)，并满足如下条件：

$\mathit{\Omega } = \bigcup\limits_i^{{{{N}}_{\rm{a}}}} {T_i^a} , \quad T_i^a \cap T_j^b = \emptyset \quad i \ne j$

在网格T_i^a中，物理量ϕ的分布函数为：

${\phi ^a}(\mathit{\boldsymbol{x}}) = \phi _i^a(\mathit{\boldsymbol{x}})\quad {\rm{ if }}\mathit{\boldsymbol{x}} \in T_i^a$

式中：ϕ_i^a(x)可以指定多种形式，通常为多项式分布。

网格T^b是计算域Ω的另一种形式的划分，划分形式与要求同T^a，则同样有T^b上的物理量分布ϕ^b(x)。如果对于Ω的任意子区域D，满足如下条件：

$\int_D {{\phi ^a}} (\mathit{\boldsymbol{x}}){\rm{d}}V = \int_D {{\phi ^b}} (\mathit{\boldsymbol{x}}){\rm{d}}V\quad \forall D \subset \mathit{\Omega }$

(1)

则称ϕ^b(x)是ϕ^a(x)的守恒插值结果。

对于一般的情形，要使式(1)满足，必须有T^b≡T^a，ϕ^b(x)≡ϕ^a(x)，但这就失去了在两套网格间进行插值的意义。为此，对D进行限制，不能是Ω的任意子区域，而必须是网格T^b中的某个网格：

$\int_D {{\phi ^a}} (\mathit{\boldsymbol{x}}){\rm{d}}V = \int_D {{\phi ^b}} (\mathit{\boldsymbol{x}}){\rm{d}}V\quad \forall D \subset \left\{ {T_i^b|i = 1, 2, \cdots , {N_b}} \right\}$

(2)

上式就是一般意义上的守恒插值关系式。

对式(2)，按照基于单元相交的思想，对D进行分解：

$D = T_i^b = \sum\limits_{j = 1}^{{N_a}} {\left( {T_i^b \cap T_j^a} \right)} = \sum\limits_{j = 1}^{{N_a}} {{V_{ij}}}$

式中：V_ij=T_i^b∩T_j^a。将上式代入式(2)，则有：

$\int_D {{\phi ^b}} (\mathit{\boldsymbol{x}}){\rm{d}}V = \int_D {{\phi ^a}} (\mathit{\boldsymbol{x}}){\rm{d}}V = \sum\limits_{j = 1}^{{N_a}} {\int_{{V_{ij}}} {{\phi ^a}} } (\mathit{\boldsymbol{x}}){\rm{d}}V$

(3)

式(3)就是基于单元相交守恒插值方法(CIB/DC)的基本依据。从式(2)到式(3)未引入任何近似。因此，在给定分布ϕ(x)的情况下，CIB/DC方法是严格守恒的。

从式(3)还可以看出，如果已知网格T^a上的物理量分布ϕ^a(x)，如何计算网格单元间的相交区域V_ij便成为CIB/DC方法的关键。

2. “超网格”基本思想

图 1为文献[8]给出的超网格示意图，超网格是由T^a和T^b所有结点以及它们边界交点构成的网格，具有如下特性：

图 1 超网格示意图

Figure 1. Schematic of the supermesh

下载: 全尺寸图片幻灯片

(1) 每一个超网格不再和T^a或T^b中任何网格单元相交，即它总可以在T^a中找到一个网格包含或等于它，也总可以在T^b中找到一个网格包含或等于它；

(2) T^a的每一个网格总可以明确表示为一个或多个超网格的集合，T^b的每一个网格总可以明确表示为一个或多个超网格的集合。

这样一来，对于T^b任意网格D可以写出它所包含的超网格，每一个超网格就是V_ij=T_i^b∩T_j^a，可以在T^a上找到具体位置和对应的物理量，采用ϕ^b(x)≡ϕ^a(x)，按照D所包含的超网格进行求和，就得到D的物理量，从而实现了从T^a到T^b的信息传递。由于T^a和T^b在超网格的交集完全等价，使传递的物理量严格守恒，实现精确守恒插值。

3. 单元相交算法

由上节可知，获得超网格是实现精确守恒插值的关键。下面给出计算平面凸多边形和三维凸多面体交集的具体算法。

3.1 平面凸多边形相交算法

如图 2所示，将目标三角形A₁A₂A₃视为切割多边形，使用切割三角形B₁B₂B₃对目标三角形进行切割，每次切割时只保留目标多边形中与切割多边形在切割线同侧的部分。图中依次使用△B₁B₂B₃的3条边B₁B₂、B₂B₃、B₃B₁切割△A₁A₂A₃，切割后保留的多边形依次为CA₃A₂G、CDEA₂G、CDEFB₁，多边形CDEFB₁即为△A₁A₂A₃和△B₁B₂B₃的交集。

图 2 二维网格切割示意图

Figure 2. Schematic of the two-dimensional grid cutting

下载: 全尺寸图片幻灯片

在使用直线去切割平面凸多边形时，以图 3为例，切割直线经过点O，其法向单位矢量为n，切割直线将平面划分为两部分，定义法向n所指的半平面为正半平面，反方向的半平面为负半平面。如果切割之后目标多边形仅保留正半平面的部分，那么切割操作等效于求目标多边形和正半平面的交集。定义平面上任意一点P到切割直线的距离为：

$d = {\mathit{\boldsymbol{x}}_{OP}} \cdot \mathit{\boldsymbol{n}} = \left( {{\mathit{\boldsymbol{x}}_P} - {\mathit{\boldsymbol{x}}_O}} \right) \cdot \mathit{\boldsymbol{n}}$

图 3 直线切割平面凸多边形

Figure 3. Schematic of a straight line cuttingplane convex polygon

下载: 全尺寸图片幻灯片

式中：x_P和x_P分别为点P和点O坐标。

根据上述定义，正半平面的点到切割直线的距离均为正数，负半平面的点到切割直线的距离均为负数，切割直线上的点到切割直线的距离为零。如果凸多边形各顶点到切割直线的距离均为非负数(图 3(a))，则切割结果即为原多边形；如果凸多边形各顶点到切割直线的距离均为非正数，则切割结果为空集；除上述两种情况外，凸多边形与切割直线必然有两个交点(图 3(b))。具体如下：

(1) 设目标多边形为N，其顶点依次为N₁, N₂, …, N_N；切割直线为l，切割结果为R；

(2) 计算N的各个顶点到l的距离，记顶点N_i到l的距离为d_i；

(3) 初始化一个空的有序点集P；

(4) 从顶点N₁起，依次对多边形N的边N₁N₂, N₂N₃, …, N_NN₁进行判定。

对线段N_iN_j的判定过程为：如果d_id_j＜0，则线段与切割直线有新的交点，记交点为N_new；依次将N_i、N_new、N_j中距离为非负数的顶点加入点集P；如果d_id_j≥0，则依次将N_i、N_j中距离为非负数的顶点加入点集P；

(5) 完成所有边的判定后，如果P为空集，则R为空集；若P为非空集，则依次连接P中的点形成的多边形就是所求的切割结果R。

3.2 三维凸多面体相交算法

在三维流体计算中，常用的网格包括四面体、四棱锥、三棱柱、六面体等多种类型。其中，四棱锥、三棱柱、六面体这3种网格都具有四边形表面，网格生成过程中不能严格保证四边形表面的四个顶点位于同一平面内，会给多面体表面的定义带来歧义。

为了简化算法的复杂性，同时消除歧义，在计算两个任意凸多面体的交集时，首先将凸多面体分解为多个小四面体，如图 4所示，其中点N₆为点N₂、N₃、N₄、N₅的重心；然后通过计算这些小四面体的交集，来获得多面体的交集。计算多面体P_a和P_b的交集的具体步骤为：

图 4 三维网格分解示意图

Figure 4. Schematic of the three dimensional mesh decomposition

下载: 全尺寸图片幻灯片

(1) 将P_a分解为一组四面体{T_aⁱ}，将P_b分解为一组四面体{T_b^j}，并且满足：

$\begin{array}{l} \bigcup\limits_{i = 1}^{{N_a}} {T_a^i} = {P_a}, {\rm{ }}T_a^i \cap T_a^j = \emptyset {\rm{ }}\;\forall 1 \le i, {\rm{ }}j \le {N_a}\;\;且i \ne j\\ \bigcup\limits_{j = 1}^{{N_b}} {T_b^j} = {P_b}, {\rm{ }}T_b^i \cap T_b^j = \emptyset {\rm{ }}\;\forall 1 \le i, {\rm{ }}j \le {N_b}\;\;且i \ne j \end{array}$

(2) 计算{T_aⁱ}和{T_bⁱ}中任意两个四面体的交集：

$T_c^{ij} = T_a^i \cap T_b^j\;\;1 \le i \le {N_a}{\rm{, }}1 \le j \le {N_b}$

并记录交集T_c^ij的体积和重心分别为V_c^ij，x_c^ij。

(3) 多面体P_a和P_b交集可表示为：

${P_a} \cap {P_b} = \left( {\bigcup\limits_{i = 1}^{{N_a}} {T_a^i} } \right) \cap \left( {\bigcup\limits_{j = 1}^{{N_b}} {T_b^j} } \right) = \bigcup\limits_{i = 1}^{{N_a}} {\bigcup\limits_{j = 1}^{{N_b}} {\left( {T_a^i \cap T_b^j} \right)} } = \bigcup\limits_{i = 1}^{{N_a}} {\bigcup\limits_{j = 1}^{{N_b}} {T_b^j \cap T_c^{ij}} }$

交集的体积和重心分别是：

${V_c} = \sum\limits_{i, j} {V_c^{ij}} , \quad {x_c} = \frac{1}{{{V_c}}}\sum\limits_{i, j} {V_c^{ij}} \mathit{\boldsymbol{x}}_c^{ij}$

该方法的核心在于如何计算两个四面体的交集T_c^ij。相对于平面多边形的情形，即使四面体是最简单的多面体，其交集的计算也是十分复杂的：两个四面体相交，最多可能产生一个12顶点的8面体。本文中使用切割法计算两个四面体的交集：将一个四面体看做目标四面体，另一个四面体为切割四面体，依次使用切割四面体的4个表面去切割目标四面体，最终留下的多面体就是这两个四面体的交集。

4. 相交单元的搜索算法

相对于通常的流场插值方法，守恒插值的计算效率较低，主要有两个原因：(1)计算两个单元交集的算法较为复杂；(2)对同一个新网格单元，普通插值方法只需用到一个旧网格单元(称为贡献单元)的信息，而守恒插值需要用到多个贡献单元的信息。前文已对单元相交算法进行了详细的描述，本节将介绍守恒插值中多个贡献单元的搜索算法。

在贡献单元的搜索算法中，首先需要快速定位出新网格单元在旧网格中的位置。对此，可采用四叉树^[11]或八叉树^[12](三维)方法加速搜寻，如图 5(a)所示，完全填充网格为旧网格，其上的非填充网格为新网格中的某个单元，通过四叉树方法，可快速定位出新网格单元的格心所在的旧网格单元(图中填充单元)，这是新网格单元的第1个贡献单元；然后，逐个搜索各个贡献单元的直接相邻单元，如果相邻单元与新网格单元相交，则将该单元加入贡献单元的列表，如图 5(b)中浅色填充单元所示；当所有贡献单元的相邻单元(图 5(c)与(b)中不同的填充单元)都不与新网格单元相交时，搜索结束。

图 5 相交单元的搜索方法

Figure 5. Search method of intersection elements

下载: 全尺寸图片幻灯片

5. 守恒插值验证算例

设计了3个算例来验证本文发展的守恒插值方法。第1个算例主要考核相同拓扑结构下网格点位置变化的影响；第2个算例主要考核网格形状的影响；第3个算例主要考核三维情况下的应用情况。

算例 1

考虑如下的测试函数：

$\phi (x, y) = 1 + \sin (2{\rm{ \mathsf{ π} }}x)\sin (2{\rm{ \mathsf{ π} }}y)$

计算区域为[0, 1]×[0, 1]，计算网格为一组拓扑结构相同的非结构三角形网格，分别表示为M₀、M₁、…、M_N，其中M₀如图 6所示。网格M_N中的第(i, j)个网格点的坐标由下式给出：

$\left\{ {\begin{array}{*{20}{l}} {x(\xi , \eta , t) = (1 - \alpha (t))\xi + \alpha (t){\xi ^2}}\\ {y(\xi , \eta , t) = (1 - \alpha (t))\eta + \alpha (t){\eta ^2}} \end{array}} \right.$

图 6 数值计算建模图

Figure 6. Geometry of the specimenused for simulation

下载: 全尺寸图片幻灯片

式中：α(t)=0.5sin(4πt), ξ=(i-1)/(i_max-1), η=(j-1)/(j_max-1)t=n/T, 其中i=1, 2, …, i_max; j=1, 2, …, j_max; n=0, 1, …, N。

测试过程中，在网格M₀上给定测试函数的准确分布，然后依次插值到网格M₁、M₂、…、M_N上，共进行N次插值，最后比较网格M_N上的函数分布与准确值之间的误差。在本文计算中，取i_max=j_max=65，N=200，T=20。

图 7给出了计算误差随插值次数的变化情况，其中纵坐标代表计算区域内各单元函数分布插值后计算值与理论值之差的1范数。两种插值方法中，插值误差都随插值次数的增加而增大，但守恒型插值方法的误差始终较小，不到常规二阶插值误差的一半。

图 7 计算结果的误差比较

Figure 7. Different errors of computational results

下载: 全尺寸图片幻灯片

图 8为两种方法的物理量积分值比较，其中纵坐标代表每个单元计算值总和与理论值的比，从图中明显可以看出，守恒插值保证了计算域内物理量的守恒。

图 8 计算结果的积分比较

Figure 8. Different integrals of computational results

下载: 全尺寸图片幻灯片

算例2

考虑如下的测试函数：

$\phi (r) = 2 + \cos ({\rm{ \mathsf{ π} }}r/L)$

计算区域是边长为1的正方形，L为正方形的对角线长度，r是与正方形中心的距离。测试中使用两种不同的网格，网格1是均匀分布的结构网格，网格2是非结构网格，并与网格1具有相同的边界网格点分布，如图 9所示。测试中，首先在网格1上给出测试函数的准确分布，然后将其插值到网格2上，再插值回网格1，如此反复，进行200次插值后，将物理量分布与初始的准确分布进行比较，统计插值过程中的误差。分别采用常规二阶插值和守恒型二阶插值进行测试，并使用不同疏密的网格以分析插值方法的精度。测试中使用的4组网格的网格量如表 1所示。

图 9 数值计算建模图

Figure 9. Geometry of the specimen used for simulation

下载: 全尺寸图片幻灯片

表 1 计算网格参数

Table 1. Parameters of the computational grid

序号	结构网格	非结构网格
序号	结构网格	节点数	单元数
1	33×33	695	1 260
2	65×65	2 593	4 928
3	129×129	10 044	19 574
4	257×257	39 687	78 348

下载: 导出CSV

| 显示表格

图 10是计算误差与网格步长之间的关系，其中横坐标是表 1中4个不同结构网格计算单元的边长，纵坐标代表计算区域内各单元函数分布插值后计算值与理论值之差的1范数。可以看出在相同的网格步长下，守恒型插值的误差更小，而且误差与网格步长之间的变化斜率更大，说明其具有更高阶的精度。计算结果表明，常规二阶插值精度约为1.4阶，而守恒型二阶插值精度约为1.87阶。表 2是不同算例的守恒误差比较，守恒插值方法的守恒误差接近机器零，明显优于普通二阶插值。

图 10 插值方法的精度比较

Figure 10. Comparison of precision with different interpolation methods

下载: 全尺寸图片幻灯片

表 2 不同算例的守恒误差

Table 2. Conservation errors indifferent numerical cases

序号	二阶插值	守恒型二阶插值
1	3.724×10^-3	0
2	-1.447×10^-3	0
3	-3.985×10^-4	1×10^-14
4	5.636×10^-5	-2×10^-14

下载: 导出CSV

| 显示表格

算例3

考虑如下的测试函数：

$\varphi (x, y, z) = 1 + \sin (2{\rm{ \mathsf{ π} }}x)\sin (2{\rm{ \mathsf{ π} }}y)\sin (2{\rm{ \mathsf{ π} }}z)$

计算区域是边长为1的正方体。每次测试使用两套网格，首先在网格1上给出测试函数的准确分布，然后将其插值到网格2上，再插值回网格1，如此反复，进行10次插值后，将物理量分布与初始的准确分布进行比较，统计插值过程中的误差。分别采用普通二阶插值和守恒型二阶插值进行测试，并使用不同疏密的网格以分析插值方法的精度。测试中使用的4组网格的网格量如表 3所示，其中网格步长由下式得到：

${\rm{d}}x = {\left( {\sqrt {{N_{E1}}{N_{E2}}} } \right)^{1/{N_D}}}$

表 3 计算网格参数

Table 3. Parameters of the computational grid

序号	网格1		网格2		dx
序号	节点数	单元数	节点数	单元数	dx
1	914	4 033	676	2 711	0.067
2	2 222	10 632	1 424	6 075	0.050
3	6 527	32 913	3 733	16 681	0.035
4	26 099	138 542	12 600	59 347	0.022

下载: 导出CSV

| 显示表格

式中：N_E1和N_E2分别为两套网格的单元个数，N_D=3为网格的维数。

图 11是测试函数在三维空间中的分布云图，图 12是计算误差与网格步长之间的关系，其中横坐标是表 3中dx的值，纵坐标代表计算区域内各单元函数分布插值后计算值与理论值之差的1范数。在三维条件下，守恒型插值依然具有更小的计算误差和更高的计算精度。

图 11 测试函数分布云图

Figure 11. Nephogram of different test functions

下载: 全尺寸图片幻灯片

图 12 插值方法的精度比较

Figure 12. Comparison of precision with different interpolation methods

下载: 全尺寸图片幻灯片

6. 结论

采用在新旧网格相交单元剖分构建出来的超级网格作为网格之间信息传递的基础，采用四叉树(二维)或八叉树(三维)方法搜寻算法提高计算效率，基于以上方法提出了网格之间流场物理量的守恒插值算法。该算法理论上可以实现网格之间物理量的精确守恒，二维和三维算例也表明产生的误差接近机器零，明显优于常用的二阶插值。本文中提出的方法可以解决爆轰流场模拟中网格尺度和物理量相差较大情况下的插值难题。

图 1 典型尖卵形侵彻弹体尺寸

Figure 1. Sizes of a typical ogiven projectile

下载: 全尺寸图片幻灯片

图 2 弹体SNL-00-02侵彻23 MPa混凝土时的刚体过载曲线比较

Figure 2. Comparison of rigid-body deceleration curves for projectiles (SNL-00-02) penetrating 23 MPa concrete

下载: 全尺寸图片幻灯片

图 3 弹径缩比系数 ${\lambda _D}$ 与弹长缩比系数 ${\lambda _b}$ 、空腔长度缩比系数 ${\lambda _e}$ 的对应关系

Figure 3. Corresponding relationships among the diameter scaling coefficient ${\lambda _D}$ , the length scaling coefficient ${\lambda _b}$ and the cavity length scaling coefficient ${\lambda _e}$

下载: 全尺寸图片幻灯片

图 4 弹体侵彻薄靶的偏转过程示意图

Figure 4. Deflection trajectory of a projectile penetrating into a thin target

下载: 全尺寸图片幻灯片

图 5 弹身长度缩比系数λ_b与弹径缩比系数λ_D、空腔长度缩比系数λ_e的对应取值

Figure 5. The corresponding range of the length reduction coefficient of projectile body λ_l, the diameter reduction coefficient λ_D and the cavity length reduction coefficient λ_e

下载: 全尺寸图片幻灯片

图 6 原型试验和缩比弹试验弹体侵深

Figure 6. Penetration depth of prototype projectile and reduced scale projectile for two types projectile

下载: 全尺寸图片幻灯片

图 7 原型试验和缩比弹试验刚体过载对比

Figure 7. Rigid-body decelerations of prototype projectiles and reduced scale projectiles for two types projectiles

下载: 全尺寸图片幻灯片

图 8 1 000 kg级原型弹与缩比弹侵彻10层混凝土靶板的刚体过载

Figure 8. Rigid-body decelerations of 1 000 kg prototype projectiles and reduced scale projectiles penetrating ten-layer concrete target

下载: 全尺寸图片幻灯片

图 9 1 000 kg级原型弹与缩比弹的有限元模型对比

Figure 9. Models of 1 000 kg class prototype projectile and reduced scale projectile

下载: 全尺寸图片幻灯片

图 10 缩比弹自由状态下初始形状与第12阶拉伸模态振型的位移云图对比

Figure 10. Tensile mode shapes of the 12th order modes of scale projectiles in free state

下载: 全尺寸图片幻灯片

图 11 弹体侵彻十层靶板1/4有限元模型

Figure 11. 1/4 finite element model of a projectile penetrating a ten-layer target plates

下载: 全尺寸图片幻灯片

图 12 原型弹和缩比弹贯穿10层薄靶过程刚体过载时程曲线

Figure 12. Rigid-body deceleration curve of a prototype projectile and a scale projectile while penetrating 10 layers of thin target

下载: 全尺寸图片幻灯片

表 1 模型和原型物理量缩比因子

Table 1. Scaling factors of physical quantities between model and prototype

变量	量纲	原型与模型参量比	缩比因子
特征尺寸	L	l_p/l_m	n
时间	T	t_p/t_m	n
质量	M	M_p/M_m	n³
密度	ML⁻³	ρ_p/ρ_m	1
应力	ML⁻¹T⁻²	σ_p/σ_m	1
过载	LT⁻²	a_p/a_m	1/n
速度	LT⁻¹	V_p/V_m	1

下载: 导出CSV

表 2 原型弹的弹靶参数

Table 2. Expemrimental parameters of prototype projectiles and targets

弹型	D/mm	M/kg	l/mm	η_CRH	v₀/(m∙s⁻¹)	ρ_t/(kg∙m^-3)	f_c/MPa
Ⅰ	250	500	2 000	3.5	450	2 420	45
Ⅱ	380	1 200	2 500	3.5	850	2 420	45

下载: 导出CSV

表 3 缩比系数表

Table 3. Scaling parameters of prototype projectiles

缩比弹序号	λ_D	λ_M	λ_b	λ_e
1	1.50	2.25	1.20	1.20
2	2.00	4.00	1.30	1.20
3	2.50	6.25	1.50	1.20
4	3.00	9.00	1.60	1.20

下载: 导出CSV

表 4 Ⅱ型（1 000 kg级）弹原型试验及对应的缩比试验参数

Table 4. Prototype test and corresponding scaling test parameters of type II (1 000 kg) projectile

量	原型试验	缩比试验	相关缩比系数
M/kg	1 200	80	λ_M=15
l/mm	2 500	890	λ_l=2.8
D/mm	380	150	${\lambda _D}{\rm{ = }}2.5$
η_CRH	3	3
v₀/(m∙s⁻¹)	850	850	${\lambda _{{v_0}}} = 1$
θ₀/(°)	0/10/15	0/10/15
H/mm	300+180×9	300+140×9	λ_H=1和1.35
z/mm	3 000	2 900	${\lambda _z}{\rm{ = }}0.97$
σ_bc/MPa	40	30	λ_fc=0.75
注：z为靶间距；σ_bc为靶材料抗压强度。

下载: 导出CSV

表 5 1 000 kg级原型弹及其缩比弹自由状态下前20阶模态频率

Table 5. Top 20 modal frequencies of 1 000 kg class prototype projectile and reduced scale projectile in free state

模态序号	模态频率/Hz		模态序号	模态频率/Hz
模态序号	原型弹	缩比弹	模态序号	原型弹	缩比弹
7	335.67	1 159.2	14	1 416.2	4 670.5
8	335.67	1 159.2	15	1 453.2	4 670.5
9	716.72	2 132.5	16	2 012.8	6 291.2
10	833.99	2 708.3	17	2 012.8	6 291.2
11	833.99	2 708.3	18	2 200.1	6 463.9
12	1 085.9	3 196.4	19	2 236.0	7 044.5
13	1 416.2	4 563.3	20	2 633.9	8 244.7

下载: 导出CSV

参考文献(30)

[1]	王树有, 顾晓辉, 赵有守. 混凝土侵彻试验相似准则验证分析 [J]. 南京理工大学学报(自然科学版), 2005, 29(5): 43–46. DOI: 10.3969/j.issn.1005-9830.2005.05.011. WANG S Y, GU X H, ZHAO Y H. Experimental analysis of simularity criteria for concrete penetration [J]. Journal of Nanjing University of Science and Technology (Natural Science), 2005, 29(5): 43–46. DOI: 10.3969/j.issn.1005-9830.2005.05.011.
[2]	王世虎. 硬目标侵彻中的加速度信号研究[D]. 北京: 北京理工大学, 2010: 5−30.
[3]	FORRESTAL M J, FREW D J, HICKERSON J P, et al. Penetration of concrete targets with deceleration-time measurements [J]. International Journal of Impact Engineering, 2003, 28(5): 479–497. DOI: 10.1016/S0734-743X(02)00108-2.
[4]	FREW D J, FORRESTAL M J, Cargile J D. The effect of concrete target diameter on proticle deceleration and penetration depth [J]. International Journal of Impact Engineering, 2006, 32(10): 1584–1594. DOI: 10.1016/j.ijimpeng.2005.01.012.
[5]	XU Y, KEER L M, LUK V K. Elastic-cracked model for penetration into un-reinforced concrete targets with ogival nose projectiles [J]. International Journal Solids Structures, 1997, 34(12): 1479–1491. DOI: 10.1016/S0020-7683(96)00099-6.
[6]	LUNDGREN. A strain gage based projectile health monitor and salvage indicating circuit for kinetic energy penetrating projectiles [C] // 53th NDIA Fuze Conferenc. Florida: NDIA, 2009.
[7]	刘小虎, 刘吉, 王乘, 等. 弹丸低速垂直侵彻无钢筋混凝土的实验研究 [J]. 爆炸与冲击, 1999, 19(4): 323–328. LIU X H, LIU J, WANG C, et al. Experimental studies on the projectile penetrating normally into a plain concrete [J]. Explosion and Shock Waves, 1999, 19(4): 323–328.
[8]	杨明, 杨志刚, 林祖森. 加速度计输出信号的信息组成 [J]. 中北大学学报(自然科学版), 2000, 21(1): 44–46. DOI: 10.3969/j.issn.1673-3193.2000.01.012. YANG M, YANG Z G, LIN Z S. Information compositions of the output signal of accelerometers [J]. Journal of North China Institute of Technology (Natural Science Edtition), 2000, 21(1): 44–46. DOI: 10.3969/j.issn.1673-3193.2000.01.012.
[9]	ZHANG W, CHEN L, XIONG J, et al. Ultra-high g deceleration-time measurement for the penetration into steel target [J]. International Journal of Impact Engineering, 2007, 34(3): 436–447. DOI: 10.1016/j.ijimpeng.2006.01.008.
[10]	徐鹏, 祖静, 范锦彪. 高速动能弹侵彻硬目标加速度测试技术研究 [J]. 振动与冲击, 2007, 26(11): 118–122. DOI: 10.3969/j.issn.1000-3835.2007.11.028. XU P, ZU J, FAN J B. Study on acceleration test technique of high velocity kinetic energy projectile penetrating into hard target [J]. Journal of Vibration and Shock, 2007, 26(11): 118–122. DOI: 10.3969/j.issn.1000-3835.2007.11.028.
[11]	何丽灵, 高进忠, 陈小伟, 等. 弹体高过载硬回收测量技术的实验探讨 [J]. 爆炸与冲击, 2013, 33(6): 608–612. DOI: 10.11883/1001-1455(2013)06-0608-0612. HE L L, GAO J Z, CHEN X W, et al. Experimental study on measurement technology for projectile deceleration [J]. Explosion and Shock Waves, 2013, 33(6): 608–612. DOI: 10.11883/1001-1455(2013)06-0608-0612.
[12]	周宁, 任辉启, 沈兆武, 等. 侵彻钢筋混凝土过程中弹丸过载特性的实验研究 [J]. 实验力学, 2006, 21(5): 572–578. DOI: 10.3969/j.issn.1001-4888.2006.05.005. ZHOU N, REN Q H, SHEN Z W, et al. Experimental study on overload characteristics of projectile penetrating reinforced concrete [J]. Journal of Experimental Mechanics, 2006, 21(5): 572–578. DOI: 10.3969/j.issn.1001-4888.2006.05.005.
[13]	黄家蓉, 刘瑞朝, 何翔, 等. 侵彻过载测试信号的数据处理方法 [J]. 爆炸与冲击, 2009, 29(5): 555–560. DOI: 10.11883/1001-1455(2009)05-0555-0560. HUANG J R, LIU R C, HE X, et al. A new data processing technique for measured penetration overloads [J]. Explosion and Shock Waves, 2009, 29(5): 555–560. DOI: 10.11883/1001-1455(2009)05-0555-0560.
[14]	赵生伟, 初哲, 李明. 抗侵彻过载战斗部装药安定性实验研究 [J]. 兵工学报, 2010, 31(S1): 284–287. ZHAO S W, CHU Z, LI M. Experiment investigation on stability of explosive in anti-overload warhead [J]. Acta Armamentarii, 2010, 31(S1): 284–287.
[15]	张会锁, 罗旭, 张远高. 弹体过载记录仪安装方式对侵彻过载峰值的影响分析 [J]. 中北大学学报(自然科学版), 2014, 35(3): 252–257. DOI: 10.3969/j.issn.1673-3193.2014.03.006. ZHANG H S, LUO X, ZHANG Y G. Research on the influence of the missile overload recorder installation style on penetration acceleration peak value [J]. Journal of North University of China (Natural Science Edition), 2014, 35(3): 252–257. DOI: 10.3969/j.issn.1673-3193.2014.03.006.
[16]	李计林, 徐文峥, 王晶禹, 等. 侵彻过程中弹载火工品过载特性数值模拟 [J]. 火工品, 2009(6): 30–34. DOI: 10.3969/j.issn.1003-1480.2009.06.010. LI J L, XU W J, WANG J Y, et al. The numerical simulation for overload characteristic of initiating explosive device on missile during penetration [J]. Initiators and Pyrotechnics, 2009(6): 30–34. DOI: 10.3969/j.issn.1003-1480.2009.06.010.
[17]	HALDAR A, HAMIEH H A. Local effect of solid missiles on concrete structures [J]. Journal of Structural Engineering, 1984, 110(5): 948–960. DOI: 10.1061/(ASCE)0733-9445(1984)110:5(948).
[18]	徐建波. 长杆射弹对混凝土的侵彻特性研究[D]. 湖南: 国防科学技术大学, 2001: 23−56.
[19]	杨超, 赵宝荣, 付克勤, 等. 缩比件弹体侵彻混凝土过程相似律研究 [J]. 兵器材料科学与工程, 2003(5): 3–7. DOI: 10.3969/j.issn.1004-244X.2003.05.001. YANG C, ZHAO B R, FU K Q, et al. Research on similarity law of penetration concrete of equiscale projectile [J]. Ordnance Material Science and Engineering, 2003(5): 3–7. DOI: 10.3969/j.issn.1004-244X.2003.05.001.
[20]	武海军, 黄风雷, 陈利, 等. 动能弹侵彻钢筋混凝土相似性分析 [J]. 兵工学报, 2007, 28(3): 276–280. DOI: 10.3321/j.issn:1000-1093.2007.03.005. WU H J, HUANG F L, CHEN L, et al. Similarity law analyses of penetration behavior in reinforced concrete [J]. Acta Armamentarii, 2007, 28(3): 276–280. DOI: 10.3321/j.issn:1000-1093.2007.03.005.
[21]	陈小伟, 张方举, 杨世全, 等. 动能深侵彻弹的力学设计(III): 缩比实验分析 [J]. 爆炸与冲击, 2006, 26(2): 105–114. DOI: 10.11883/1001-1455(2006)02-0105-0114. CHEN X W, ZHANG F J, YANG S Q, et al. Mechanics of structural design of epw (Ⅲ): investigations on the reduced-scale tests [J]. Explosion And Shock Waves, 2006, 26(2): 105–114. DOI: 10.11883/1001-1455(2006)02-0105-0114.
[22]	MEBAR Y. A method for scaling ballistic penetration phenomena [J]. International Journal of Impact Engineering, 1997, 19(9): 821–829. DOI: 10.1016/S0734-743X(97)00020-1.
[23]	YOUNG C W. Equations for predicting earth penetration by projectiles: an update: SAND88-0013 [R]. USA: Sandia National Laboratories, 1988.
[24]	LI Q M, CHEN X W. Dimensionless formulae for penetration depth of concrete target impacted by a non-deformable projectile [J]. International Journal of Impact Engineering, 2003, 28(1): 93–116. DOI: 10.1016/S0734-743X(02)00037-4.
[25]	朱彤. 结构动力模型相似问题及结构动力试验技术研究[D]. 大连: 大连理工大学, 2004.
[26]	MAI Y W, ATKINS A G. Crack propagation in non-proportionally scaled elastic structures [J]. International Journal of Mechanical Sciences, 1978, 20(7): 437–449. DOI: 10.1016/0020-7403(78)90033-4.
[27]	FORRESTAL M J, TZOU D Y. A spherical cavity-expansion penetration model for concrete targets [J]. International. Journal. Solids Structure, 1997, 34(31-32): 4127–4146. DOI: 10.1016/S0020-7683(97)00017-6.
[28]	MA Z F, DUAN Z P, OU Z C, et al. The experimental and theoretical research on attitude of projectile obliquely penetrating into thin concrete target [J]. Acta Armamentarii, 2015, 36(S1): 248–252. DOI: 10.15918/j.tbit1001-0645.2016.10.003.
[29]	王冰. 大长径比弹体侵彻混凝土的数值模拟研究[D]. 北京: 北京理工大学, 2012: 5−30.
[30]	PENG Y, WU H, FANG Q, et al. Residual velocities of projectiles after normally perforating the thin ultra-high performance steel fiber reinforced concrete slabs [J]. International Journal of Impact Engineering, 2016, 97(11): 1–9. DOI: 10.1016/j.ijimpeng.2016.06.006.

施引文献

期刊类型引用(7)

1.	余昕，吴乘胜，金奕星. 面向交界面网格的超网格算法研发及CFD应用. 中国造船. 2023(03): 233-246 . 百度学术
2.	王子维，刘东健，陈逖. 单级跨声速压气机旋转失速模拟方法. 气体物理. 2021(01): 30-37 . 百度学术
3.	常思源，白晓征，崔小强，刘君. 一种改进的非定常激波装配算法. 航空学报. 2020(02): 163-176 . 百度学术
4.	刘耘臻，万志强，杨超. 飞行载荷外部气动力的二次规划等效映射方法. 北京航空航天大学学报. 2020(03): 541-547 . 百度学术
5.	孙伟，陈泽栋. 风载荷对飞行器尾罩分离影响的数值模拟分析. 弹箭与制导学报. 2019(06): 102-106 . 百度学术
6.	崔鹏程，唐静，李彬，马明生，邓有奇. 基于超网格的重叠网格守恒插值方法. 航空学报. 2018(03): 23-35 . 百度学术
7.	靳晨晖，王刚，王泽汉. 子母弹多体分离过程的非定常CFD/RBD数值仿真. 气体物理. 2018(04): 47-63 . 百度学术

其他类型引用(3)

资源附件(0)

访问统计

图(12) / 表(5)

计量

文章访问数: 5168
HTML全文浏览量: 1628
PDF下载量: 112
被引次数: 10

1. 守恒插值基本思想
2. “超网格”基本思想
3. 单元相交算法
3.1 平面凸多边形相交算法
3.2 三维凸多面体相交算法
4. 相交单元的搜索算法
5. 守恒插值验证算例
6. 结论

非等比例缩比侵彻/贯穿相似规律研究

doi: 10.11883/bzycj-2019-0086

作者简介: 刘 源（1994－ ），男，硕士研究生，lymoxing@163.com

通讯作者: 皮爱国（1977－ ），男，博士，副教授，aiguo_pi@bit.edu.cn 黄风雷（1965－ ），男，博士，教授，huangfl@bit.edu.cn

计量

出版历程

Study on similarity law of non-proportionally scaled penetration/perforation test

1. 守恒插值基本思想

2. “超网格”基本思想

3. 单元相交算法

3.1 平面凸多边形相交算法

3.2 三维凸多面体相交算法

4. 相交单元的搜索算法

5. 守恒插值验证算例

6. 结论

期刊类型引用(7)

其他类型引用(3)

计量

出版历程

目录

1. 守恒插值基本思想

2. “超网格”基本思想

3. 单元相交算法

3.1 平面凸多边形相交算法

3.2 三维凸多面体相交算法

4. 相交单元的搜索算法

5. 守恒插值验证算例

6. 结论

作者简介:
刘　源（1994－　），男，硕士研究生，lymoxing@163.com

通讯作者:
皮爱国（1977－　），男，博士，副教授，aiguo_pi@bit.edu.cn

黄风雷（1965－　），男，博士，教授，huangfl@bit.edu.cn