Numerical Simulation of Penetration in Concrete Sheet Based on SPH Method

QIANG Hongfu; SUN Xinya; WANG Guang; CHEN Fuzhen

doi:10.11858/gywlxb.20180634

Volume 33 Issue 2

Apr 2019

Turn off MathJax

Article Contents

Abstract

References

Chinese Journal of High Pressure Physics > 2019 > 33(2): 024101.

CHEN Chuan, GONG Zizheng, YANG Wulin, LI Ming, YU Qian. Influence of Geometry of Space Debris on Laser Ablation Impulse[J]. Chinese Journal of High Pressure Physics, 2018, 32(4): 040101. doi: 10.11858/gywlxb.20180508

Citation:

QIANG Hongfu, SUN Xinya, WANG Guang, CHEN Fuzhen. Numerical Simulation of Penetration in Concrete Sheet Based on SPH Method[J]. Chinese Journal of High Pressure Physics, 2019, 33(2): 024101. doi: 10.11858/gywlxb.20180634

Citation:

PDF( 16376 KB)

Numerical Simulation of Penetration in Concrete Sheet Based on SPH Method

doi: 10.11858/gywlxb.20180634

Missile Engineering College, Rocket Force University of Engineering, Xi’an 710025, China

Received Date: 14 Sep 2018
Rev Recd Date: 24 Sep 2018

Abstract

Abstract

As the continuous improvement of the strength of concrete structure, more and more protective fortifications have chosen concrete as the main building material. On the basis of the smooth particle hydrodynamics method, the TCK-HJC composite constitutive model is proposed to simulate the deformation damage of the concrete target during the penetration process of the rigid conical projectile, and the pseudo-fluid model is used to describe the failed concrete fragments. The deformation process of the conical projectile penetrating into the concrete target under different penetration angles (0°, 60°) were analyzed, and the pressure on the surface of the target was also obtained. The scattering angles of the broken concrete fragments were obtained, and the comparison with experiment shows the simulation result is robust and accurate, which lays a technical foundation for further simulation of the mechanical properties of brittle materials.
- concrete,
- TCK-HJC composite constitutive model,
- smoothed particle hydrodynamics method,
- brittle materials

FullText(HTML)

空间碎片是人类航天活动产生的太空垃圾。随着航天活动的日益频繁，空间环境日趋恶化，空间碎片对航天器安全和航天员生命构成了严重的现实威胁，是空间长期可持续发展面临的严重挑战。按照目前空间碎片的增长趋势推算，在未来50年左右空间碎片密度将达到一个临界值，引发碎片间的链式撞击效应(凯斯勒效应)，之后近地空间将彻底不可用^[1]。当前已经采用的空间碎片减缓、碰撞规避、防护等措施都不能从根本上遏制空间碎片数量的快速增长。为了彻底扭转空间碎片环境不断恶化的趋势，保持空间碎片环境态势的稳定，乃至达到清洁空间的目的，对空间碎片进行主动移除已经成为国际航天界的当务之急和唯一选择^[2]。在众多的空间碎片主动移除技术中，激光烧蚀驱动技术具有操作简便、反应迅速、成本低、可靠性高以及可多次使用的特点，不仅能够高效移除尺寸为1~10cm的危险空间碎片，也能够防止空间碎片与航天器或碎片之间的相互撞击，可有效地减少空间碎片的数量，成为当前空间碎片主动移除的研究热点之一^[2-4]。

激光烧蚀驱动移除空间碎片技术的原理是：高功率激光辐照碎片表面，使光斑区材料熔化、汽化、等离子体化，高温高压等离子体和气体物质飞散形成反喷等离子体羽流，反喷羽流与空间碎片的动量交换使碎片获得一个与其相反的冲量，从而驱动碎片变轨，使低轨碎片逐步落入大气层烧毁^[5]。激光烧蚀驱动碎片原理如图 1所示^[2]。

图 1 激光驱动碎片原理

Figure 1. Schematic of driving debris by laser

下载: 全尺寸图片幻灯片

目前，国内外对激光与空间碎片相互作用的冲量耦合规律和冲量耦合系数变化规律都已经进行了大量深入研究^[6-7]。但这些研究主要是针对静止状态下理想的平面目标开展的，而真实的空间碎片的几何外形绝大多数都是不规则的，不规则外形会对激光辐照驱动效果造成很大影响。为了实现激光烧蚀对空间碎片的精准驱动，必须研究其不规则的几何外形对激光驱动冲量的影响。目前，由于尚缺乏不规则外形目标冲量耦合系数的实验测量数据，因此普遍采用的方法是将其离散为多个已有冲量耦合系数数据的平面目标，分别计算驱动效果后再进行累加。通过该思路金星等^[8]、张品亮等^[9]分别研究了覆盖式光斑和点光斑作用于立方体、圆柱体等典型几何外形目标的驱动效果。美国劳伦斯利弗莫尔国家实验室的Liedahl等^[10-11]研究了激光作用于外形不规则碎片的一般计算方法，但是该方法不仅对每一特定形状目标都必须根据其外形做针对性的分析计算，而且对目标外形信息的准确度、完整度要求甚高，导致该方法使用起来非常复杂，不具备普适性和通用性，基本无法真正在工程中进行应用。

为解决这一问题，本研究提出了一种基于平面目标实验结果，通过对几何形状不规则目标进行表面三角化重构进而简化计算激光驱动效果的方法。根据复杂外形目标表面的顶点信息，基于在实验中测得的同一平面目标材料的冲量耦合系数，就可对其激光驱动效果进行精确计算。以几种典型的外形规则目标和不规则目标为例，应用该方法，研究目标几何形状对激光驱动冲量的影响规律。

1. 计算方法

理想的激光辐照烧蚀驱动平面目标的基本公式为^[12]

$m\Delta v = {C_{\rm{m}}}E$

(1)

式中：m为目标质量(单次脉冲作用过程中，烧蚀质量可忽略不计^[13])，Δv为目标产生的速度增量，C_m为冲量耦合系数，E为单次激光脉冲辐照到目标表面的能量。

对于外形不规则目标，在计算激光烧蚀驱动行为时，首先需要分别计算每个受照射表面元的冲量，然后相加获得目标总冲量，如图 2所示。根据国内外大量的实验研究，无论激光入射方向如何，等离子体羽流喷射方向始终沿着烧蚀区域法向^{[10, 14]}, 因此，有计算公式如下^[7]

$m\Delta \mathit{\boldsymbol{v}} = - {C_{\rm{m}}}I\sum\limits_i {{S_i}} \left| {\mathit{\boldsymbol{k}}\cdot{\mathit{\boldsymbol{n}}_i}} \right|{\mathit{\boldsymbol{n}}_i}$

(2)

图 2 激光照射外形不规则目标

Figure 2. Schematic of irregularly shaped target under laser irradiation

下载: 全尺寸图片幻灯片

式中：S_i为第i个面元面积，k为照射激光单位向量，n_i为第i个表面元法向量，I为激光能量密度。

对于存在曲面的外形不规则目标，则需筛选出该曲面被照射部分区域，对其进行积分，获得该面产生的冲量，如图 3所示，计算公式如下：

${{(m\Delta \mathit{\boldsymbol{v}})}_{i}}=-{{C}_{\text{m}}}I\iint{\left| \mathit{\boldsymbol{k}}\cdot {{\mathit{\boldsymbol{n}}}_{i}} \right|{{\mathit{\boldsymbol{n}}}_{i}}}\text{d}S$

(3)

图 3 激光照射曲面

Figure 3. Schematic of curved surface under laser irradiation

下载: 全尺寸图片幻灯片

因为不规则目标表面形状各异、大小不一，且平面与曲面混杂，各表面面积和法向量难以计算获取，不规则曲面方程尤其难以得到，因而以上公式在实际中难以应用。为解决这一问题，我们提出的新思路是：根据不规则目标表面顶点坐标信息，将其表面三维重构成多个三角面组合，对每个三角面建立通用化的计算方法进行计算，最后对所有三角面进行求和。具体算法如下。

(1) 不规则目标表面三维重构

计算过程中输入的目标信息为不规则目标表面顶点信息。根据拓扑学中关于墨卡托投影Delaunay三角剖分的基本理论，对于任何一个凸多面体，经过拓扑变换后都可以等价于球形，其表面点和凸多面体的顶点一一对应，且表面顶点之间的相互位置关系并不随着拓扑变换而改变^[15]。因此，对其三角化可以等效为对所对应的球面上的点进行三角化。

具体处理步骤为：第1步，将目标表面所有顶点坐标变换到以其几何中心为原点的直角坐标系下，并投影到单位球面上；第2步，筛选出单位球面的南、北极极点，以最小范围划定南、北极圈，对其内的顶点进行三角化(三角扇形)；第3步，将剩余球面通过墨卡托投影展开成二维平面，并对其内定点采用Bowyer-Watson算法进行Delaunay三角化；第4步，将获得的三角化后的点面关系信息返回到原几何体，即获得其表面三角化重构信息。这一过程如图 4所示。

图 4 三角化三维重构过程

Figure 4. Schematic of three-dimensional triangulation reconstruction process

下载: 全尺寸图片幻灯片

(2) 基本参数计算

每个三角面的3个顶点和不规则目标质心构成一个四面体，将不规则目标分割为若干个四面体的组合，以每个四面体为一个基本计算单元，如图 5中四面体OABC即为基本计算单元。分别计算每个计算单元的体积、质量及针对质心的转动惯量，并通过累加获得不规则目标的总体积、质量及针对质心的转动惯量。

图 5 基本计算单元

Figure 5. Basic computing unit

下载: 全尺寸图片幻灯片

(3) 驱动效果计算

针对不规则目标的每个三角面，筛选出被激光照射的面，图 5中三角面ABC即为基本激光照射面，n_i为该面法向量，k为入射激光向量。根据(1)式及该表面材料的冲量耦合系数实验值，计算每个被照射表面所获得的冲量，并将各表面冲量分解为沿x、y、z 3轴的冲量和以过质心分别与x、y、z 3轴平行的直线为转轴的冲量矩，进一步累加计算后获得该次激光脉冲作用后目标的速度和转动角速度改变量。

本研究提出的方法仅需获得目标表面顶点信息即可通过三维重构对激光照射驱动效果进行计算。在实际中通过雷达光学等探测手段返回的目标不规则表面信息也多为表面顶点坐标形式，因此对于各种复杂的不规则目标，其表面顶点信息较容易获得，便于工程应用。

2. 典型几何形状目标算例

将表面三角化重构计算方法编写为C语言程序，并选取立方体、球体和圆柱体3个典型几何形状目标进行计算，通过对比本研究方法计算结果与已有公式计算结果，验证该方法的准确性。

2.1 立方体

立方体因其定点明确、边界规整，可直接进行三角化，也可在其边界和表面获得插值点后进行计算，其不同三角化效果如图 6所示。

图 6 不同顶点数量下立方体三角化重构

Figure 6. Triangulation reconstruction of cube with different vertex numbers

下载: 全尺寸图片幻灯片

对于边长为a的立方体，根据(2)式有

$m\Delta \mathit{\boldsymbol{v}}=-{{C}_{\text{m}}}I{{a}^{2}}\sum\limits_{i}{\left| \mathit{\boldsymbol{k}}\cdot {{\mathit{\boldsymbol{n}}}_{i}} \right|}{{\mathit{\boldsymbol{n}}}_{i}}$

(4)

具体到三维直角坐标系中各边分别与3个坐标轴平行的立方体，在受到单位矢量为(cos θ sin α, cos θ cos α, sin θ)、能量密度为I的激光覆盖照射时，单次脉冲所获得冲量为

$\left\{ \begin{align} & {{P}_{x}}=\text{cos}\theta \text{sin}\alpha {{C}_{\text{m}}}I{{a}^{2}} \\ & {{P}_{y}}=\text{cos}\theta \text{sin}\alpha {{C}_{\text{m}}}I{{a}^{2}} \\ & {{P}_{z}}=\text{sin}\theta {{C}_{\text{m}}}I{{a}^{2}} \\ \end{align} \right.$

(5)

以密度为2.7g/cm³的典型铝合金材料为例，其冲量耦合系数为0.000 06N·s·J^-1。对于边长为1cm的铝合金立方体，当其受到方向矢量为(-1，-1，-1)、到靶表面能量密度为5.7J/cm²的激光覆盖照射时，其公式计算结果与三角化计算结果对比如表 1所示。

表 1 不同顶点数量下立方体三角化计算结果

Table 1. results of cube with different vertex numbers

Vertex number	Mass/ g	Moment of inertia/ (10^-4g·m²)			Impulse/ (mN·s)			Velocity/ (mm·s^-1)			Angular velocity/ (rad·s^-1)
Vertex number	Mass/ g	J_x	J_y	J_z	P_x	P_y	P_z	Δv_x	Δv_y	Δv_z	ω_x	ω_y	ω_z
8	2.7	1.31	1.31	1.31	1.96	1.96	1.96	-2.493	-2.493	-2.493	0	0	0
14	2.7	1.31	1.31	1.31	1.96	1.96	1.96	-2.493	-2.493	-2.493	0	0	0
50	2.7	1.31	1.31	1.31	1.96	1.96	1.96	-2.493	-2.493	-2.493	0	0	0
130	2.7	1.31	1.31	1.31	1.96	1.96	1.96	-2.493	-2.493	-2.493	0	0	0
Formula	2.7	1.31	1.31	1.31	1.96	1.96	1.96	-2.493	-2.493	-2.493	0	0	0

下载: 导出CSV

| 显示表格

由表 1中计算结果可看出，对于立方体，三角化计算结果与公式计算的结果完全一致，且不受插值点数量和三角面划分的影响。

2.2 球体

因球体表面都是曲面并没有顶点，所以需先对其表面插值获得顶点信息才能通过三角化进行计算，不同插值点数量下三角化如图 7所示。

图 7 不同顶点数量下球体三角化

Figure 7. Triangulation reconstruction of sphere with different vertex numbers

下载: 全尺寸图片幻灯片

对于半径为R的球体，根据(3)式有

$m\Delta \mathit{\boldsymbol{v}}=\frac{2}{3}{{C}_{\text{m}}}I\text{ }\!\!\pi\!\!\text{ }{{R}^{2}}\mathit{\boldsymbol{k}}$

(6)

对密度为2.7g/cm³、半径为2.5cm的铝合金球体，冲量耦合系数为0.000 06N·s·J^-1，当其受到方向矢量为(-1，-1，-1)、到靶表面能量密度为5.7J/cm²的激光覆盖照射时，其公式计算结果与三角化计算结果对比如表 2所示。

表 2 不同顶点数量下球体三角化计算结果

Table 2. Calculation results of sphere with different vertex numbers

Vertex number	Mass/ g	Moment of inertia/ (g·m²)			Impulse/ (mN·s)			Velocity/ (mm·s^-1)			Angular velocity/ (10^-7rad·s^-1)
Vertex number	Mass/ g	J_x	J_y	J_z	P_x	P_y	P_z	Δv_x	Δv_y	Δv_z	ω_x	ω_y	ω_z
72	158.926 603 1	0.109	0.109	0.105	2.693	2.691	2.534	5.820	5.820	5.477	2 220	121	10 200
222	170.834 427 5	0.122	0.122	0.121	2.744	2.746	2.715	5.518	5.521	5.459	6.74	6.98	3 790
382	172.723 740 5	0.124	0.124	0.123	2.766	2.766	2.718	5.501	5.500	5.406	2.22	4.95	3 250
762	174.891 641 2	0.126	0.126	0.126	2.769	2.768	2.764	5.438	5.437	5.429	2.02	5.55	2 820
1 452	175.755 229 0	0.127	0.127	0.127	2.775	2.775	2.771	5.425	5.424	5.415	1.71	3.79	217
2 452	176.064 045 8	0.128	0.128	0.128	2.779	2.779	2.773	5.422	5.422	5.407	1.42	3.46	192
Formula	176.625	0.129	0.129	0.129	2.780	2.780	2.780	5.406	5.406	5.406	0	0	0

下载: 导出CSV

| 显示表格

从表 2中可看出，当在球表面插值72点(经线8点、纬线10点)时，计算结果与公式结果有较大偏差。几何体参数方面：如图 7(a)所示，在球表插值72点划分方法下，三角化重构的几何体与球体差距较大，导致计算所得球体质量小于公式结果，相对偏差为10.02%；由于丢失的质量都位于球体表面附近，进一步导致了更大的转动惯量损失，对x(y)轴转动惯量的相对偏差为15.50%，对z轴转动惯量的相对偏差为18.61%；且由于经线和纬线插值点数目不同以及z轴与x、y轴投影方向表面三角化划分方式不同，导致z轴和其他两轴的参数值不同。激光作用效果方面：由于球体表面插值72点三角化重构的表面与原球面的偏差，使激光作用点位置和法向量角度均发生改变，导致计算所得冲量也与公式计算结果不同，x(y)轴方向冲量的相对偏差为3.13%，z轴方向冲量的相对偏差为8.85%；x(y)轴方向速度增量相对偏差为7.66%，z轴方向速度增量相对偏差为1.31%。

对同一球体采用更密集的插值点划分方式：222点(经线13点、纬线20点，见图 7(b))、382点(经线21点、纬线20点，见图 7(c))、762点(经线21点、纬线40点，见图 7(d))、1 452点(经线31点、纬线50点，见图 7(e))、2 452点(经线51点、纬线50点, 见图 7(f))。

从表 2中计算结果可看出，插值点的加密显著降低了与公式计算结果的偏差，随着插值点数的增加，偏差缩小。2 452个插值点时：质量相对偏差为0.32%；转动惯量相对偏差(x、y、z轴)均为0.78%；x(y)轴冲量的相对偏差为0.36%, z轴冲量的相对偏差为0.25%；x(y)速度的相对偏差为0.30%，z轴速度的相对偏差为0.02%。可见，当插值点达到一定数量时，三角化计算值已与公式计算值十分接近。

据此可看出，在对球体采用本研究方法计算时，插值点数量对计算精度有较大影响。插值点过少时，粗糙的三角化将导致较大的计算偏差，插值点的增加能显著提高计算精度。

2.3 圆柱体

因圆柱体侧面是曲面并没有顶点，所以需先对其底圆和母线插值获得顶点信息才能通过三角化重构进行计算，不同插值点数量下三角化如图 8所示。

图 8 不同顶点数量下圆柱体三角化

Figure 8. Triangulation reconstruction of cylinder with different vertex numbers

下载: 全尺寸图片幻灯片

对半径为R、高为H的圆柱体，受到单位矢量为(cos θsin α, cos θcos α, sin θ)、能量密度为I的激光覆盖照射时，根据(2)式、(3)式有

$\left\{ \begin{align} & {{P}_{x}}=\frac{\text{ }\!\!\pi\!\!\text{ }}{2}\text{sin}\alpha \text{cos}\theta {{C}_{\text{m}}}IHR \\ & {{P}_{y}}=\frac{\text{ }\!\!\pi\!\!\text{ }}{2}\text{cos}\alpha \text{cos}\theta {{C}_{\text{m}}}IHR \\ & {{P}_{z}}=\text{sin}\theta {{C}_{\text{m}}}I\text{ }\!\!\pi\!\!\text{ }{{R}^{2}} \\ \end{align} \right.$

(7)

对底面半径为2.5 cm、高为7 cm的铝合金圆柱体，冲量耦合系数为0.000 06 N·s·J^-1，当其受到方向矢量为(-1，-1，-1)、到靶表面能量密度为5.7 J/cm²的激光覆盖照射时，其公式计算结果与三角化计算结果对比如表 3所示。

表 3 不同顶点数量下圆柱体三角化计算结果

Table 3. Calculation results of cylinder with different vertex numbers

Vertex number	Mass/ g	Moment of inertia/ (g·m²)			Impulse/ (mN·s)			Velocity/ (mm·s^-1)			Angular velocity/ (10^-8rad·s^-1)
Vertex number	Mass/ g	J_x	J_y	J_z	P_x	P_y	P_z	Δv_x	Δv_y	Δv_z	ω_x	ω_y	ω_z
32	347.290 1	0.560 5	0.560 5	0.295 7	5.305	5.302	3.603	5.249	5.249	3.565	3.72×10^-4	4.75×10^-4	2.52
62	365.152 7	0.597 2	0.597 2	0.326 7	5.369	5.369	3.788	5.053	5.053	3.565	1.21×10^-8	5.23×10^-8	2.50
122	369.618 3	0.606 7	0.606 7	0.334 8	5.386	5.386	3.835	5.006	5.006	3.565	3.52×10^-16	2.76×10^-16	2.39
162	369.618 3	0.606 7	0.606 7	0.334 8	5.386	5.386	3.835	5.006	5.006	3.565	3.52×10^-16	2.75×10^-16	2.39
202	369.618 3	0.606 7	0.606 7	0.334 8	5.386	5.386	3.835	5.006	5.006	3.565	3.52×10^-16	2.78×10^-16	2.39
242	370.305 3	0.608 5	0.608 5	0.336 4	5.388	5.388	3.844	4.997	4.997	3.565	2.01×10^-16	1.92×10^-16	2.40
Formula	370.992 4	0.609 9	0.609 9	0.337 6	5.389	5.389	3.849	4.988	4.988	3.563	0	0	0

下载: 导出CSV

| 显示表格

从表 3中结果可看出，当在圆柱表面插值32点(底圆10点、母线3点)时，计算结果与已有公式结果有一定差距。几何体参数方面：如图 8(a)所示，该划分方法下，三角化重构的几何体与圆柱体差距较大，导致计算所得圆柱体质量小于公式结果，相对偏差为6.39%；进而导致了更大的转动惯量损失，x、y轴转动惯量的相对偏差为8.10%, z轴转动惯量的相对偏差为12.41%。激光作用效果方面：由于重构后底面与圆柱体底圆差距过大，一方面导致三角化重构表面与原圆柱体侧面曲面产生偏差，使激光作用点位置和法向量角度均发生改变；另一方面导致底面受辐照面积减小，使计算所得x、y、z轴冲量均与公式计算结果不同，x、y轴方向冲量的相对偏差为1.56%、z轴方向冲量的相对偏差为6.39%。在x、y轴方向，由于质量和冲量两者的共同偏差，导致速度增量出现5.23%的相对偏差；z轴方向则由于三角化导致的激光辐照表面损失和质量损失等比例，圆柱体三角化后的面质比保持不变，因而该方向速度相对偏差很小，为0.056%。

对同一圆柱体采用更密集的插值点划分方式：62点(底圆20点、母线3点，见图 8(b))、122点(底圆40点、母线3点，见图 8(c))、162点(底圆40点、母线4点，见图 8(d))、202点(底圆80点、母线3点，见图 8(e))、242点(底圆60点、母线4点，见图 8(f))。其中，62点、122点方案均为细化底圆曲线上的插值，从表 3中可看出，随着插值点的增加，计算精度显著提高。162点和202点分别细化了母线上和底圆平面中的插值，计算结果与122点时完全一致，说明对直线插值点加密以及对平面三角化的细化无法提高计算精度，误差来源于曲线和曲面。242点是在162点的基础上进一步加密了底圆的插值点，偏差进一步缩小：质量相对偏差为0.19%；x、y轴转动惯量相对偏差为0.23%，z轴转动惯量相对偏差为0.36%；x、y轴冲量相对偏差为0.020%，z轴冲量相对偏差为0.13%；x、y轴速度相对偏差为0.18%，z轴速度相对偏差为0.056%。

2.4 讨论

根据以上对立方体、圆柱体、球体计算结果的对比和分析，可见：(1)对于立方体这类直线棱划分的平面包裹的几何体，三角化重构计算方法准确度很高，其结果与公式计算结果完全一致; (2)直线(如立方体的棱、圆柱的母线)上插值点的划分密度对计算结果没有影响，可直接采用几何体天然顶点; (3)对于圆柱体侧面和球体表面等曲面的三角化精度和曲线上插值点的密度会大大影响计算精度，曲线上插值点越密集，曲面三角化后表面越精确，计算结果越准确。

在实际应用中，不规则形状目标的曲面和曲线信息往往是未知的，无法像以上规则几何体一样根据顶点信息和曲线方程进行插值加密，其计算精度就取决于初始获得的目标表面顶点信息。初始获得的目标表面顶点越多、越密集、越准确，三角化重构出来的目标表面就越接近目标真实情况，在此基础上进行的激光驱动计算结果也就越精确。

通过该方法还可以反过来研究激光烧蚀驱动碎片技术中对目标探测系统精度的要求。首先，选取一个已知准确外形信息且包含复杂表面细节的不规则几何体，通过三角化计算获得精确的作用效果。然后，通过逐渐删除顶点逐渐减少信息的完整度和准确度来模拟探测精度下降对目标细节的丢失，并依次将计算结果与精确结果对比，获得不同精度下的误差变化规律。最后，根据碎片移除、主动防御的不同操控精度需求可获得不同的目标信息精度需求，进而提出对探测系统的不同精度要求。

3. 激光辐照外形规则目标产生的冲量规律

在激光辐照驱动目标过程中，一般关注的驱动效果主要有两个：一是驱动产生冲量的方向，二是驱动产生冲量的大小。一般期待的冲量方向为沿激光辐照方向，而目标外形会导致实际冲量方向偏离激光照射方向，并使沿激光辐照方向的冲量减小。因此对目标外形因素导致的冲量方向与激光方向的偏角β (见图 9)和目标外形对沿激光辐照方向冲量大小的影响进行研究。

图 9 冲量偏角β

Figure 9. Diagram of impulse angle β

下载: 全尺寸图片幻灯片

3.1 立方体

根据(5)式可知，单位矢量为(cos θsin α, cos θcos α, sin θ)的激光覆盖照射立方体时冲量方向矢量为(cos θsin α, cos θcos α, sin θ)，与作用激光矢量一致，因此其几何外形不影响激光照射后产生冲量的方向。

立方体沿激光作用方向的合冲量大小始终为C_mIa²，等于其单个表面接收到同强度激光垂直照射时获得的冲量大小。

3.2 球体

由于球体的对称性，任意角度激光覆盖式辐照方向必然过球心，所以其冲量方向始终和激光照射方向相同，其几何外形不影响激光照射后产生的冲量方向。

球体沿激光作用方向的合冲量大小始终为 $\frac{2}{3}{C_{\rm{m}}}I{\rm{ \mathsf{ π} }}{R^2}$ ，等于相同半径平面目标受到同等强度激光垂直照射时获得的冲量大小的2/3。

3.3 圆柱体

激光照射方向与x轴的夹角为γ，如图 10所示。计算得到不同γ下沿激光方向的冲量占激光照射同截面平面目标冲量的比例如图 11所示，不同γ下冲量与激光方向的偏角β如图 12所示。

图 10 激光照射圆柱体

Figure 10. Schematic of cylinder under laser irradiation

下载: 全尺寸图片幻灯片

图 11 圆柱体在不同角度激光照射下产生的冲量占平面目标冲量的比例

Figure 11. Impulse percentage of cylinder to planner target irradiated by laser at different angles

下载: 全尺寸图片幻灯片

图 12 圆柱体在不同角度激光照射下产生的冲量偏角

Figure 12. Impulse angle of cylinder irradiated by laser at different angles

下载: 全尺寸图片幻灯片

从图 11中可看出：激光从正侧面照射圆柱体时冲量损失比例为22%，冲量损失最小；激光从顶部和底部照射时冲量损失比例为为44%，冲量损失最大。从图 12中可看出，激光从圆柱体对称轴方向入射时冲量偏角为0°，离对称轴越远冲量偏角越大，最大为9.6°。综合两者结果，从圆柱体侧边垂直方向照射能在获得最大冲量的同时使冲量方向与激光方向一致。

3.4 圆锥体

对高10 cm、底面半径2.5 cm的圆锥体，激光照射方向与x轴夹角为γ，如图 13所示。计算得到不同γ下沿激光方向的冲量占激光照射同截面平面目标冲量的比例如图 14所示，不同γ下产生的冲量与激光方向的偏角β如图 15所示。

图 13 激光照射圆锥体

Figure 13. Schematic of cone under laser irradiation

下载: 全尺寸图片幻灯片

图 14 圆锥体在不同角度激光照射下产生的冲量占平面目标冲量的比例

Figure 14. Impulse percentage of cone to planner target irradiated by laser at different angles

下载: 全尺寸图片幻灯片

图 15 圆锥体在不同角度激光照射下产生的冲量偏角

Figure 15. Impulse angle of cone irradiated by laser at different angles

下载: 全尺寸图片幻灯片

从图 14中可看出：激光从圆锥底部和侧边方向入射时冲量损失比例均较小，分别为22%和16%；激光从圆锥顶点方向入射时冲量损失比例最大，为80%。从图 15中可看出：激光除了从对称轴方向入射以外，从垂直于侧面方向和底角所对方向入射时也有两个偏角接近0°极小值；在两个侧面方向偏角最大，为48.1°。综合两者结果，从圆锥体底面垂直方向照射能在获得最大冲量的同时使冲量方向与激光方向一致。

4. 激光辐照外形不规则目标产生的冲量规律

这里以小行星为例。外形不规则目标选择编号为101955的小行星贝努，如图 16所示，是一个典型的外形不规则目标。该小行星发现于1999年9月，NASA对其化学、物理和运动学特性进行了大量的测量，获得了丰富的数据，并于2016年9月8日发射了“奥西里斯”探测器，计划在2019年与其交会期间对其进行环绕探测和接触采样。贝努是一颗对地球有威胁的“潜在危险天体”(Potentially Hazardous Asteroid，PHA)，最小轨道交会距离约0.002 AU(1 AU=1.495 978 7×10¹¹ m)，其撞击地球的危险性目前排在PHA第2位(在2182年)。激光烧蚀驱动是抵御小行星撞击地球的可能手段之一，因此，以贝努为算例除了可以检验本研究所提方法外，同时也有重要的实际价值。

图 16 行星贝努^[16]

Figure 16. Asteroid Bennu^[16]

下载: 全尺寸图片幻灯片

贝努的平均直径为492 m，赤道尺寸为565 m×535 m，自旋周期为4.297 d(角速度1.69×10^-5 rad/s)。从NASA网站上的小行星数据库可获得包含其顶点坐标信息的几何外形数据^[16]，其三角化后结果如图 17所示。

图 17 贝努表面三角化重构

Figure 17. Triangulation reconstruction of the Bennu

下载: 全尺寸图片幻灯片

贝努为碳质小行星，表面存在富含碳元素的风化层，因此取其冲量耦合系数为0.000 086 N·s·J^-1，据估算贝努密度为1.181 g/cm³。选取激光方向矢量为(-1，-1，-1)，到靶表面能量密度为5.7 J/cm²。在其坐标xy平面内改变激光照射角度γ, 计算得到不同γ下沿激光辐照方向的冲量占激光照射同截面平面目标及同半径球体冲量的比例如图 18所示，不同γ下冲量与激光方向的偏角β如图 19所示。

图 18 贝努在不同角度激光照射下产生的冲量占平面/球体目标的比例

Figure 18. Impulse percentage of Bennu to plane/sphere target irradiated by laser at different angles

下载: 全尺寸图片幻灯片

图 19 贝努在不同角度激光照射下产生的冲量偏角

Figure 19. Impulse angle of Bennu irradiated by laser at different angles

下载: 全尺寸图片幻灯片

根据结果可知：由于贝努近似一个扁球体，不同角度激光辐照下产生的冲量大小受外形影响较小，在90°照射角附近冲量损失最小，相比同直径规则球体损失了17%，相比同截面平面目标损失了43%；其外形对冲量偏角的绝对值影响较小，最大值为8°，但不同方向上差别较大，在90°和180°附近冲量偏角达到最小值。

5. 结论

基于目标表面三角化三维重构提出了一种可以精确计算激光辐照外形不规则目标产生的冲量大小及方向的方法，以立方体、球体和圆柱体3个典型的外形规则目标为对象，验证了方法的计算精度。利用该方法研究了激光辐照典型外形规则和不规则目标产生的冲量规律，结果显示：对于立方体和球体，由于其自身的对称性，其外形对激光辐照产生的冲量方向无影响，但是会导致冲量大小发生改变；对于圆柱体和圆锥体，其外形对激光辐照产生的冲量大小和方向均有影响，影响大小随激光入射角度的不同而改变，在特定入射角度下，圆柱体和圆锥体均能在不影响激光作用冲量方向的情况下获得最大冲量。对于小行星“贝努”这种类似球体的不规则目标，其外形会使激光作用冲量的大小产生较大损失，而对冲量方向的影响相对较小，但受照射角度影响明显。这一方法及研究结论对深化激光烧蚀驱动移除空间碎片技术及其工程应用有重要价值。

References(33)

References

[1]	HOLMQUIST T J, TEMPLETON D W, BISHNOI K D. Constitutive modeling of aluminum nitride for large strain, high-strain rate, and high-pressure applications [J]. International Journal of Impact Engineering, 2001, 25(3): 211–231. doi: 10.1016/S0734-743X(00)00046-4
[2]	孙炜海, 鞠桂玲, 杨班权. 平头弹丸侵彻B₄C陶瓷/金属复合靶板的数值模拟 [J]. 装甲兵工程学院学报, 2014, 28(2): 45–48 SUN W H, JU G L, YANG B Q. Numerical simulation on penetration of B₄C ceramic/metal composite targets struck by flat-ended projectiles [J]. Journal of Academy of Armored Force Engineering, 2014, 28(2): 45–48
[3]	BACKMAN M E, GOLDSMITH W. The mechanics of penetration of projectiles into targets [J]. International Journal of Engineering Science, 1978, 16(1): 1–99. doi: 10.1016/0020-7225(78)90002-2
[4]	钱伟长. 穿甲力学 [M].北京: 高等教育出版社, 1984. QIAN W C. Penetration mechanics [M]. Beijing: High Education Press, 1984.
[5]	CORBETT G G, REID S R, JOHNSON W. Impact loading of plates and shells by free-flying projectiles: a review [J]. International Journal of Impact Engineering, 1996, 18(2): 141–230. doi: 10.1016/0734-743X(95)00023-4
[6]	FORRESTAL M J, FREW D J, HANCHAK S J, et al. Penetration of grout and concrete targets with ogive-nose steel projectiles [J]. International Journal of Impact Engineering, 1996, 18(5): 465–476. doi: 10.1016/0734-743X(95)00048-F
[7]	LI Q M, REID S R, WEN H M, et al. Local impact effects of hard missiles on concrete targets [J]. International Journal of impact engineering, 2005, 32(1): 224–284. doi: 10.1016/j.ijimpeng.2005.04.005
[8]	陈小伟. 穿甲/侵彻问题的若干工程研究进展 [J]. 力学进展, 2009, 39(3): 316–351 doi: 10.3321/j.issn:1000-0992.2009.03.006 CHEN X W. Advances in the penetration/perforation of rigid projectiles [J]. Advances in Mechanics, 2009, 39(3): 316–351 doi: 10.3321/j.issn:1000-0992.2009.03.006
[9]	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
[10]	GEBBEKEN N, GREULICH S, PIETZSCH A. Hugoniot properties for concrete determined by full-scale detonation experiments and flyer-plate-impact tests [J]. International Journal of Impact Engineering, 2006, 32(12): 2017–2031. doi: 10.1016/j.ijimpeng.2005.08.003
[11]	梁斌. 动能攻坚战斗部对混凝土靶侵爆效应研究 [D].绵阳:中国工程物理研究院, 2009. LIANG B. Research on invading explosion effect of concrete target by kinetic energy [D]. Mianyang: China Academy of Engineering Physics, 2009.
[12]	黄民荣.刚性弹体对混凝土靶的侵彻与贯穿机理研究 [D].南京:南京理工大学, 2011. HUANG M R. Penetration and perforation mechanism of rigid projectile into the concrete target [D]. Nanjing: Nanjing University of Science & Technology, 2011.
[13]	邓佳杰, 张先锋, 葛贤坤, 等. 基于局部相互作用理论的侵彻弹头部形状优化及仿真 [J]. 爆炸与冲击, 2017, 37(4): 611–620 doi: 10.11883/1001-1455(2017)04-0611-10 DENG J F, ZHANG X F, GE X K, et al. Nose-shape optimization and simulation of projectiles penetrating into concrete target based on local interaction theory [J]. Explosion and Shock Waves, 2017, 37(4): 611–620 doi: 10.11883/1001-1455(2017)04-0611-10
[14]	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
[15]	王成, 王万军, 宁建国. 聚能装药对混凝土靶板的侵彻研究 [J]. 力学学报, 2015, 47(4): 672–684 WANG C, WANG W J, NING J G. Investigation on shaped charge penetrating into concrete targets [J]. Chinese Journal of Theoretical and Applied Mechanics, 2015, 47(4): 672–684
[16]	薛建锋, 沈培辉, 王晓鸣. 弹体斜侵彻混凝土靶的实验研究及其数值模拟 [J]. 爆炸与冲击, 2017, 37(3): 536–543 doi: 10.11883/1001-1455(2017)03-0536-08 XUE J F, SHEN P H, WANG X M. Experimental study and numerical simulation of projectile obliquely penetrating into concrete target [J]. Explosion and Shock Waves, 2017, 37(3): 536–543 doi: 10.11883/1001-1455(2017)03-0536-08
[17]	KONG X Z, WU H, FANG Q, et al. Projectile penetration into mortar targets with a broad range of steiking velocities: test and analyses [J]. International of Journal of Impact Engineering, 2017, 106: 18–29. doi: 10.1016/j.ijimpeng.2017.02.022
[18]	GUPTA N K, AITMAN B S, CARGILE J D, et al. Normalimpact of ogive nosed projectiles on thin plates [J]. International of Journal of Impact Engineering, 2001, 25(7): 641–660. doi: 10.1016/S0734-743X(01)00003-3
[19]	WARREN T L, FOSSUM A F, FREW D J. Penetration into low-strength (23 MPa) concrete: target characterization andsimulations [J]. International of Journal of Impact Engineering, 2004, 30(5): 477–503. doi: 10.1016/S0734-743X(03)00092-7
[20]	董军, 邓国强, 杨科之, 等. 弹丸对混凝土薄板的冲击破坏效应 [J]. 岩石力学与工程学报, 2005, 24(4): 713–720 doi: 10.3321/j.issn:1000-6915.2005.04.029 DONG J, DENG G Q, YANG K Z, et al. Damage effect of thin concrete slabs subjected to projectile impact [J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(4): 713–720 doi: 10.3321/j.issn:1000-6915.2005.04.029
[21]	强洪夫, 范树佳, 陈福振, 等. 基于SPH 方法的聚能射流侵彻混凝土靶板数值模拟 [J]. 爆炸与冲击, 2016, 36(4): 516–524 doi: 10.11883/1001-1455(2016)04-0516-09 QIANG H F, FAN S J, CHEN F Z, et al. Numerical simulation on penetration of concrete target by shaped charge jet with SPH method [J]. Explosion and Shock Waves, 2016, 36(4): 516–524 doi: 10.11883/1001-1455(2016)04-0516-09
[22]	张馨予, 吴艳青, 黄风雷. PBX装药弹体侵彻混凝土薄板的数值模拟 [J]. 含能材料, 2018, 26(1): 101–108 ZHANG X Y, WU Y Q, HUANG F L. Numerical simulation on the dynamic damage of PBX charges filled in projectiles during penetrating thin concrete targets [J]. Chinese Journal of Energetic Materials, 2018, 26(1): 101–108
[23]	武海军, 黄风雷, 付跃升, 等. 钢筋混凝土中爆炸破坏效应数值模拟分析 [J]. 北京理工大学学报, 2007, 27(3): 200–204 doi: 10.3969/j.issn.1001-0645.2007.03.004 WU H J, HUANG F L, FU Y S, et al. Numerical simulation of reinforced concrete breakage under internal blast loading [J]. Transactions of Beijing Institute of Technology, 2007, 27(3): 200–204 doi: 10.3969/j.issn.1001-0645.2007.03.004
[24]	SWEGLE J W, ATTAWAY S W. On the feasibility of using smoothed particle hydrodynamics for underwater explosion calculations [J]. Computational Mechanics, 1995, 17(3): 151–168. doi: 10.1007/BF00364078
[25]	LIBERSKY L D, PETSCHEK A G, CARNEY T C, et al. High strain Lagrangian hydrodynamics: a three-dimensional SPH code for dynamic material response [J]. Journal of Computational Physics, 1993, 109(1): 67–75. doi: 10.1006/jcph.1993.1199
[26]	MONAGHAN J J. Smoothed particle hydrodynamics [J]. Annual Review of Astronomy and Astrophysics, 1992, 30(1): 543–574. doi: 10.1146/annurev.aa.30.090192.002551
[27]	JOHNSON G R, COOK W. A constitutive model and data for metals subjected to large strains, high strain rates and high temperatures [C]// Proceedings of the Seventh International Symposium on Ballistics. The Hague, 1983: 541–547.
[28]	QIANG H, WANG K, GAO W. Numerical simulation of shaped charge jet using multi-phase SPH method [J]. Transactions of Tianjin University, 2008, 14(1): 495–499.
[29]	强洪夫, 孙新亚, 陈福振, 等. 钻地弹侵爆多层混凝土建筑物的SPH 数值模拟 [J/OL]. 爆炸与冲击(2018–06–02) [2018–09–14].http://kns.cnki.net/kcms/detail/51.1148.O3.20180531.1019.034.html. QIANG H F, SUN X Y, CHEN F Z, et al. Numerical simulation of earth penetrating weapon penetration and explosion multi-layer concrete structures with SPH method [J]. Explosion and Shock Waves (2018–06–02) [2018–09–14]. http://kns.cnki.net/kcms/detail/51.1148.O3.20180531.1019.034.html.
[30]	陈福振, 强洪夫, 高巍然. 气粒两相流传热问题的光滑离散颗粒流体动力学方法数值模拟 [J]. 物理学报, 2014, 63(23): 230206 doi: 10.7498/aps.63.230206 CHEN F Z, QIANG H F, GAO W R. Numerical simulation of heat transfer in gas-particle two-phase flow with smoothed discrete particle hydro dynamics [J]. Acta Physica Sinica, 2014, 63(23): 230206 doi: 10.7498/aps.63.230206
[31]	强洪夫, 高巍然. 完全变光滑长度SPH法及其实现 [J]. 计算物理, 2008, 25(5): 569–575 doi: 10.3969/j.issn.1001-246X.2008.05.008 QIANG H F, GAO W R. SPH method with fully variable smoothing lengths and implementation [J]. Chinese Journal of Computational Physics, 2008, 25(5): 569–575 doi: 10.3969/j.issn.1001-246X.2008.05.008
[32]	强洪夫.光滑例子流体动力学方法及应用 [M]. 北京: 科学出版社, 2017: 246–287.
[33]	金乾坤. 混凝土动态损伤与失效模型 [J]. 兵工学报, 2006, 27(1): 10–14 doi: 10.3321/j.issn:1000-1093.2006.01.003 JIN Q K. Dynamic damage and failure model for concrete materials [J]. Acta Armamentarii, 2006, 27(1): 10–14 doi: 10.3321/j.issn:1000-1093.2006.01.003

Relative Articles

[1]	LI Jinlin, JIANG Jianwei, MEN Jianbing, WANG Shuyou, LI Mei. Numerical Simulation of the Structure of Composite Liner to Enhance After-Effect[J]. Chinese Journal of High Pressure Physics, 2022, 36(1): 015102. doi: 10.11858/gywlxb.20210785
[2]	LÜ Chao, ZHANG Xuping, WANG Guiji, LUO Binqiang, LUO Ning, WU Heng’an, TAN Fuli, ZHAO Jianheng, LIU Cangli, SUN Chengwei. Micro-Scale Response Characteristics of Ni₅₂Ti₄₈ Alloy under Shock Loading[J]. Chinese Journal of High Pressure Physics, 2021, 35(4): 040107. doi: 10.11858/gywlxb.20210769
[3]	WANG Yansheng, LI Weibing, HUANG Xuanning, WANG Xiaoming. Orthogonal Design of the Liner Structure in Dual-Mode Charge Warhead[J]. Chinese Journal of High Pressure Physics, 2020, 34(4): 045102. doi: 10.11858/gywlxb.20200537
[4]	GUO Tengfei, LI Weibing, LI Wenbin, HONG Xiaowen. Controlling Effect of Tantalum Liner's Structural Parameters on EFP Formation and Penetration Performance[J]. Chinese Journal of High Pressure Physics, 2018, 32(3): 035104. doi: 10.11858/gywlxb.20170667
[5]	GONG Bailin, LI Ming, CHU Zhe, WU Haijun, KE Ming. Penetration Performance of Depleted Uranium Alloys Liner[J]. Chinese Journal of High Pressure Physics, 2018, 32(3): 035102. doi: 10.11858/gywlxb.20170647
[6]	HAO Yan-Jun, LIU Jin-Xu, LI Jian-Chong, FENG Wei, LIU Xu-Zhen, CHU Dong-Shan, LI Shu-Kui. Performance of Titanium Alloy Shaped Charge Liner Penetrating Pure Copper Target and Carbon Steel Target[J]. Chinese Journal of High Pressure Physics, 2017, 31(5): 535-540. doi: 10.11858/gywlxb.2017.05.005
[7]	GAO Yong-Hong, PANG Rui-Qiang, WANG Feng-Ying, GU Xiao-Hui. Effect of Material Particle Size on the Properties of Powder Liner[J]. Chinese Journal of High Pressure Physics, 2014, 28(4): 455-460. doi: 10.11858/gywlxb.2014.04.011
[8]	YU Chuan, WANG Wei, CHEN Hao, YU De-Shui, XIE Gang, ZHANG Zhen-Tao. Design of Explosively Formed Projectile Liner with Small Radius and Experiment of Penetrating Multi-Layer Steel Target[J]. Chinese Journal of High Pressure Physics, 2014, 28(1): 69-72. doi: 10.11858/gywlxb.2014.01.011
[9]	ZHAO Chang-Xiao, LONG Yuan, JI Chong, LI Yu-Chun, XIE Quan-Min, LU Liang. Distribution Law of Pressure on Liner Surface under Multi-Point Initiation[J]. Chinese Journal of High Pressure Physics, 2013, 27(1): 83-89. doi: 10.11858/gywlxb.2013.01.012
[10]	ZHENG Yu, WANG Xiao-Ming, LI Wen-Bin, LI Wei-Bing. Effects of Liner Curvature Radius on Formation of Double-Layered Spherical Segment Charge Liner into Tandem Explosively Formed Projectile (EFP)[J]. Chinese Journal of High Pressure Physics, 2009, 23(3): 229-235 . doi: 10.11858/gywlxb.2009.03.011
[11]	YU Jin-Quan, TANG Xiao-Jun, YUAN Shuai, HU Hai-Bo. Shock Induced Phase Transition of Pure Iron and Fe-Based Alloy Cuneal Samples[J]. Chinese Journal of High Pressure Physics, 2008, 22(1): 89-94 . doi: 10.11858/gywlxb.2008.01.019
[12]	SUN Shu-Hua, LI Jie, ZHAO Hai-Li, XU Rui, LIU Ri-Ping. Microstructural Characteristics in Al-Based Alloys Solidified under High Pressure[J]. Chinese Journal of High Pressure Physics, 2008, 22(4): 434-438 . doi: 10.11858/gywlxb.2008.04.017
[13]	YANG Jing-Ru, LIU Lin, REN Yan-Jun, LIU Jian-Hua, ZHANG Rui-Jun. Effect of High Pressure and Heat Treatment on Microstructure of Aluminum Bronze[J]. Chinese Journal of High Pressure Physics, 2008, 22(3): 318-322 . doi: 10.11858/gywlxb.2008.03.017
[14]	WANG Wu-Xiao, JIE Wan-Qi, YUAN Sen. Chemical Reaction of CuO-Al System in Preparation of Al2O3/Cu Composite[J]. Chinese Journal of High Pressure Physics, 2004, 18(4): 353-358 . doi: 10.11858/gywlxb.2004.04.011
[15]	ZHOU Yu-Song, WU Xi-Jun, LI Bing-Han, XU Guo-Liang. Synthesis of Bulk Samples of Nanorystalline Tungsten[J]. Chinese Journal of High Pressure Physics, 2000, 14(3): 219-223 . doi: 10.11858/gywlxb.2000.03.011
[16]	WANG Tie-Fu, WANG Lei, RUAN Wen-Jun, ZHAO Tong-Hu. The Effect of the Grain Size of a Liner on the Performance of Shaped Charge Jets[J]. Chinese Journal of High Pressure Physics, 1996, 10(4): 291-298 . doi: 10.11858/gywlxb.1996.04.009
[17]	YU Ri-Cheng, XU Da-Peng, SU Wen-Hui, PANG Yan-Xin. The Formation of Decagonal Quasicrystal and a New Crystalline Approximant in Al70Co15Tb5Ni10 Alloy[J]. Chinese Journal of High Pressure Physics, 1994, 8(2): 141-145 . doi: 10.11858/gywlxb.1994.02.008
[18]	YU Ri-Cheng, XU Da-Peng, SU Wen-Hui, PANG Yan-Xin. A Study on New Orthorhombic Approximants of the Al70Co15Tb5Ni10 Decagonal Quasi-crystal[J]. Chinese Journal of High Pressure Physics, 1993, 7(3): 208-213 . doi: 10.11858/gywlxb.1993.03.007
[19]	MAO Zi-Li, CHEN Hong, WANG Wen-Kui. Formation of Bulk Metallic Glass Zr60Ni20Al20 by High Pressure Quenching[J]. Chinese Journal of High Pressure Physics, 1992, 6(3): 212-216 . doi: 10.11858/gywlxb.1992.03.008
[20]	ZHANG Shu-Hua. Ultrasonic Equation of State for Isotropic Polycrystalline Al, Cu and Ta[J]. Chinese Journal of High Pressure Physics, 1991, 5(1): 57-61 . doi: 10.11858/gywlxb.1991.01.009

Supplements(0)

Cited By

Proportional views

Proportional views

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

Figures(11) / Tables(1)

Get Citation

PDF

XML

Article Metrics

Article views(10021) PDF downloads(56)