压力对纯的和含硫化亚铁的橄榄石电导率影响的实验研究

刘长财; 胡海英; 代立东; 孙文清

doi:10.11858/gywlxb.20180674

压力对纯的和含硫化亚铁的橄榄石电导率影响的实验研究

doi: 10.11858/gywlxb.20180674

刘长财^{1, 2,},
胡海英^1,*, ,,
代立东¹,
孙文清¹

1.
中国科学院地球化学研究所地球内部物质高温高压重点实验室，贵州贵阳　550081
2.
中国科学院大学，北京　100049

基金项目: 国家自然科学基金（41772042，41774099，41474078）；中国科学院前沿科学重点研究项目（QYZDB-SSW-DQC009）；国家重点研究发展项目（2016YFC0601101）

详细信息

作者简介:
刘长财（1992－），男，硕士研究生，主要从事高温高压下矿物岩石电学性质研究. E-mail：liuchangcai@mail.gyig.ac.cn

通讯作者:
胡海英（1983－），女，博士，副研究员，主要从事高温高压下矿物岩石物理性质研究. E-mail：huhaiying@vip.gyig.ac.cn

中图分类号: P319.2
计量
- 文章访问数: 6490
- HTML全文浏览量: 3311
- PDF下载量: 23
出版历程
- 收稿日期: 2018-11-05
- 修回日期: 2019-01-09

Experimental Study on the Effect of Pressure on the Electrical Conductivity of Pure and Iron Sulfide-Bearing Olivine

LIU Changcai^{1, 2
,},
HU Haiying^{1,*
, ,},
DAI Lidong¹,
SUN Wenqing¹

1.
Key Laboratory of High-Temperature and High-Pressure Study of the Earth’s Interior, Institute of Geochemistry, Chinese Academy of Sciences, Guiyang 550081, China
2.
University of Chinese Academy of Sciences, Beijing 100049, China

摘要

摘要: 在YJ-3000 t紧装式六面顶大腔体压机上，用Solartron-1260阻抗/增益-相位分析仪，在1～3 GPa、723～1273 K的条件下，原位测量了纯的和含15%(质量分数)FeS的橄榄石电导率。实验结果表明：在实验温度范围内，含15% FeS的橄榄石电导率比纯橄榄石的电导率高2～3个数量级，且电导率值在0.1～10 S/m范围内；纯的和含15% FeS的橄榄石电导率都随着温度的增加而增大，但是纯的橄榄石电导率对温度的敏感性更强；纯的和含15% FeS的橄榄石电导率随压力变化表现出相反的特性，随着压力的升高，纯橄榄石电导率微弱地降低，而含15% FeS的橄榄石电导率显著地增加。由含15% FeS的橄榄石电导率对温度、压力的效应以及实验获得的活化焓可知，15% FeS在橄榄石中形成了相互连通的网络，主导着橄榄石的导电过程。
- 橄榄石 /
- FeS /
- 压力 /
- 电导率
Abstract: We performed in situ electrical conductivity measurements on pure and iron FeS-bearing olivine in a multi–anvil apparatus using the impedance spectroscopy technique under the condition of 1–3 GPa and 723–1273 K. The experimental results indicated that the electrical conductivities of 15% (mass fraction) FeS–bearing olivine, in the range of 0.1–10 S/m, are 2 to 3 orders of magnitude higher than that of pure olivine in the experimental temperature range. The electrical conductivities of pure and 15% FeS-bearing olivine increase with increasing temperature. The dependence of the electrical conductivity of pure olivine on temperature is much stronger. The effect of pressure on the electrical conductivity of pure and iron FeS-bearing olivine is different. With the rise of pressure, the electrical conductivity of pure olivine slightly decreases, whereas the electrical conductivity of the 15% FeS-bearing olivine increases significantly. Based on the experimental results including the Arrhenius parameters, it is proposed that the 15% FeS can form an interconnected network in olivine, which dominates the conduction process of olivine.
- olivine /
- FeS /
- pressure /
- electrical conductivity

HTML全文

随着人们对安全要求的不断提高和人为爆炸事故的频繁发生，发生爆炸时工程的安全可靠性已逐渐成为人们关注的重要问题，深入了解爆炸载荷下结构的动态响应具有重要意义。与圆管相比，填充可压缩泡沫芯层的夹芯管能够通过冲击衰减和塑性功耗散来吸收大量爆炸能量，具有更好的抗爆性能。因此，泡沫夹芯结构作为抗爆吸能装置，引起了研究者们的广泛关注^[1-3]。

Fan等^[4]研究了短空心圆管和夹芯管在横向冲击载荷下的动态响应，发现结构的能量耗散主要通过塑性铰的形成来实现。Hall等^[5]研究了横向冲击载荷下泡沫铝、泡沫黄铜和泡沫钛填充管的能量吸收性能，结果表明，泡沫铝填充管的能量吸收特性最佳。Song等^[6]提出了横向爆炸载荷作用下钢圆管塑性变形的解析模型，预测了钢管的挠度和变形角，发现理论预测结果与实验结果吻合较好。Yuen等^[7]对横向爆炸载荷作用下的圆管进行了实验和数值模拟分析，并与Wierzbicki等^[8]提出的理论模型进行了对比，发现实验结果与数值模拟结果吻合良好，但理论模型预测的最终挠度偏低。Jing等^[9-11]研究了爆炸载荷下泡沫铝夹芯壳结构的动态响应和吸能特性，结果表明，载荷强度和几何构型影响结构的变形/破坏模式和吸能特性，增加面板厚度或泡沫铝芯层的密度可提高其抗爆性能，并建立了流固耦合、芯层压缩和整体弯曲拉伸的三阶段框架理论模型预测结构的动态响应。Hoo Fatt等^[12]考虑了泡沫芯层的弹塑性行为，对外爆炸载荷作用下的复合材料夹芯壳进行了应力波分析，结果表明，高峰值、低持续时间的冲击波作用下结构能量吸收最高，低峰值、长持续时间的冲击波反之。Liu等^[13]通过实验和数值模拟研究了两端固支的泡沫铝填充管在横向外爆炸载荷作用下的动力响应和能量吸收性能，结果表明，泡沫铝芯层的引入减小了结构跨中挠度，提高了其能量吸收能力。上述研究表明，填充泡沫的管状结构具有良好的抗冲击性能，是一种理想的吸能装置。

多目标优化可以平衡各个目标，使目标之间作出协调和折中处理，最终总体的目标尽可能达到最优。Baroutaji等^[14-15]通过多目标优化设计方法对准静态横向压缩下的圆管和泡沫铝夹芯管进行了优化，发现圆管在最大壁厚和最小直径下耐撞性最好，而夹芯管在最大芯层厚度和最小内管直径下耐撞性最好。Qi等^[16]对爆炸载荷下的泡沫铝夹芯板进行了多目标优化，通过非支配排序遗传算法（non-dominated sorting genetic algorithm, NSGA-Ⅱ）对结构的几何尺寸和芯层的相对密度进行优化，得到了最佳参数，平衡了结构的吸能能力和结构变形量。

综上所述，已有的研究主要集中于横向冲击载荷下薄壁圆管、填充管和夹芯壳结构的动态响应和能量吸收性能的研究。但是，关于泡沫金属夹芯管结构在横向爆炸载荷作用下的变形失效模式、吸能机制及抗冲击性能等方面的研究比较有限，已有的研究也局限于带曲率的夹芯壳结构，针对爆炸载荷下夹芯管结构的多目标优化的研究尚不多见。本研究基于动力显式有限元法，探究泡沫铝夹芯管在横向爆炸载荷下的动态响应，并对其抗爆性能进行多目标优化。建立横向爆炸作用下泡沫铝夹芯管的有限元模型，研究泡沫铝夹芯管的几何参数、芯层的相对密度和爆炸加载条件对内外管跨中挠度和吸能性能的影响；以内管、外管壁厚和芯层相对密度为设计变量，以芯层比吸能和外管跨中挠度为优化指标，采用响应面法（response surface method, RSM）和NSGA-Ⅱ算法实现泡沫铝夹芯管抗爆性能多目标优化设计，为抗爆结构设计提供技术指导。

1. 有限元模拟

1.1 有限元模型

泡沫铝夹芯管由内、外两个同心的薄壁圆管和泡沫铝芯层3部分组成，内、外薄壁圆管采用钢材，泡沫铝芯层采用相对密度分别为10%、15%和20%的泡沫铝材料。图1为泡沫铝夹芯管横截面的示意图，其中：D为外管外径，H为外管壁厚，d为内管外径，h为内管壁厚。试件长度为L。

图 1 泡沫铝夹芯管

Figure 1. Aluminum foam-filled sandwich tube

下载: 全尺寸图片幻灯片

采用有限元软件ABAQUS/EXPLICIT建立泡沫铝夹芯管的有限元模型，如图2所示。泡沫铝夹芯管的中线与z轴重合，内管和外管采用常规的四节点壳单元（S4R），泡沫铝芯层采用八节点线性实体单元（C3D8R），采用中性轴算法划分网格。内、外管与泡沫芯层之间定义为绑定约束，试件两端约束全部自由度来模拟两端固支的边界条件，分析过程采用动态显式算法。

图 2 有限元模型

Figure 2. Finite element model

下载: 全尺寸图片幻灯片

内外管采用各向同性硬化双线性本构模型，密度ρ=7 830 kg/m³，杨氏模量E=193 GPa，泊松比μ=0.25，屈服应力σ₀=205 MPa，在最大应变ε_UTS=0.4时，最大拉伸强度σ_UTS =520 MPa，应变率敏感性参数C和P分别取100 s⁻¹和10^[17]。泡沫铝芯层采用各向同性硬化塑性模型定义，用可压缩泡沫材料模拟泡沫铝芯层。不同相对密度泡沫铝的材料参数如表1^[18]所示，其中：ρ^*为相对密度，μ_c、μ_p分别为压缩屈服应力比和塑性泊松比。

表 1 泡沫铝材料参数^[18]

Table 1. Material parameters for aluminum foam^[18]

ρ^*/%	ρ/(g·cm⁻³)	E/MPa	μ	μ_c	μ_p
10	0.27	48	0	1.41	0.17
15	0.40	71	0	1.12	0.29
20	0.54	143	0	1.10	0.30

下载: 导出CSV

| 显示表格

1.2 爆炸载荷的假设

球形炸药作用在试件上形成了脉冲载荷，假定作用在泡沫铝夹芯管上的脉冲载荷分布如图3所示，其脉冲近似方程可写为

图 3 爆炸载荷假设^[13]

Figure 3. Blast load assumption^[13]

下载: 全尺寸图片幻灯片

$p(x,\theta ,t) = \left\{ \begin{array}{*{20}{l}} {p_0}\left( {1 - t/\tau } \right){{\rm{cos}}^2}(\text{π} x/L){\rm{sin}}\,\theta & \;\;\;t \leqslant \tau \\ 0&\;\;\; t > \tau \\ \end{array} \right.$

(1)

式中：τ为爆炸载荷响应时间，取τ=50 μs^[17]；L为试件总长，L=280 mm；p₀为作用在泡沫铝夹芯管迎爆面处的超压峰值。

根据TNT球形炸药在无限大气中的超压峰值计算公式，结合实验中的实际炸药质量和起爆距离，超压峰值^[19]（单位：10⁵ Pa）可以表示为

${p_0} = \frac{{0.085}}{Z}+\frac{{0.3}}{{{Z^2}}}+\frac{{0.8}}{{{Z^3}}}$

(2)

$Z = R/\sqrt[{{3}}]{{{W_{\rm{TNT}}}}}$

(3)

式中：W_TNT为TNT炸药质量，kg；R为起爆距离，m。实验中所采用的乳化炸药与TNT炸药的换算关系为

${W_{\rm{TNT}}} = 0.7{W_{\rm{s}}}$

(4)

式中：W_s为乳化炸药质量。结合式(1)～式(4)，可得出作用在试件上的总冲量I_T为

${I_{\rm{T}}} = 4\int_0^{L/2} {\int_0^{{\text{π}} /2} {\int_0^\tau {p(x,\theta ,t)\frac{D}{2}{\rm{d}}x{\rm{d}}\theta {\rm{d}}t = } } } \frac{1}{4}{p_0}DL\tau$

(5)

表2列出了横向爆炸载荷下两端固支的两组泡沫铝填充管试件的理论冲量与文献[13]中实验冲量的对比，其中m为泡沫铝填充管的质量，冲量I_E为实验结果，超压峰值p₀由式(2)计算得到，总冲量I_T由式(5)给出。由表2可知，理论计算冲量与实验冲量吻合较好，验证了爆炸载荷假设的合理性。

表 2 理论冲量与实验冲量的对比

Table 2. Comparison of impulse between theoretical solution and experiment

Test piece	W_s/g	R/mm	D/mm	H/mm	m/kg	p₀/MPa	I_E/(N·s)	I_T/(N·s)	$\dfrac{ {I{_{ \rm{T } } }-I{_{ \rm{E} } } } }{ {I{_{ \rm{E} } } } }\Big/$ %
No.1	30	40	76	0.7	96.8	27.7	7.1	7.4	4.2
No.2	50	40	76	0.7	96.7	45.8	12.9	12.2	−5.4

下载: 导出CSV

| 显示表格

1.3 数值模拟验证及泡沫铝夹芯管网格敏感性验证

为了验证本研究中有限元模型的有效性，根据文献[13]中的实验模型，利用ABAQUS建立有限元模型，爆炸载荷以式(1)的形式施加在泡沫铝填充管上。图4给出了试件1的数值模拟和实验的最终变形模式，其中泡沫铝夹芯管截面中的4.1 mm为最终变形区域的径向厚度尺寸。由图4可知，结构的整体变形模式与有限元结果吻合较好。表2中试件1的外管跨中挠度的实验值与模拟值分别为14.0和12.6 mm，试件2的外管跨中挠度的实验值和模拟值分别为22.1和29.0 mm，数值模拟结果与实验测试结果吻合较好，表明所建立的有限元模型是合理的。

图 4 实验^[13]与数值模拟结果的对比

Figure 4. Comparison of experimental^[13] and numerical simulation results

下载: 全尺寸图片幻灯片

为研究横向爆炸载荷下泡沫铝夹芯管的动态响应，对其有限元模型的网格敏感性进行了验证。采用乳化炸药的质量为30 g，起爆位置距离试件跨中40 mm，载荷以式(1)给出的脉冲载荷函数形式施加在泡沫铝夹芯管上。试件的几何参数：外管外径为101 mm，壁厚为0.7 mm；内管外径为54 mm，壁厚为0.7 mm，泡沫铝芯层的相对密度为15%。泡沫铝芯层和内、外管分别采用1.5、2.0、2.5和3.0 mm共4种不同的单元网格尺寸。图5给出了在横向爆炸载荷下4种不同单元网格尺寸泡沫铝夹芯管试件内外管跨中挠度时程曲线，其中虚线和实线分别代表内、外管挠度曲线，外管的挠度大于内管挠度，内外管跨中挠度随时间上下波动，取第一个波谷点后的挠度曲线的平均值作为泡沫铝夹芯管的最终挠度。当网格尺寸为1.5、2.0和2.5 mm时，三者显示了较好的一致性，而网格尺寸为3.0 mm时，数值模拟与其他网格尺寸结果相差较大。综合考虑计算精度和效率，数值模拟中采用2.0 mm的网格尺寸。

图 5 网格敏感性验证

Figure 5. Grid sensitivity verification

下载: 全尺寸图片幻灯片

2. 模拟结果与分析

为研究爆炸载荷下泡沫铝夹芯管结构的动态响应和能量吸收特性，对结构的几何参数、爆炸加载条件进行分析，研究结构的内外管直径与壁厚、泡沫铝芯层的相对密度、炸药起爆距离及炸药质量对泡沫铝夹芯管动态响应的影响。考虑9组不同几何参数的泡沫铝夹芯管试件在横向爆炸载荷下的抗爆性能，其中泡沫芯层的相对密度为15%，试件详细的几何参数如表3所示。

表 3 试件几何参数

Table 3. Geometric parameter of specimens

Test piece	D/mm	d/mm	H/mm	h/mm	L/mm	m/g
NS01	101	76	0.7	0.7	280	1 242.62
NS02	101	65	0.7	0.7	280	1 326.02
NS03	101	54	0.7	0.7	280	1 388.13
NS04	90	54	0.7	0.7	280	1 150.28
NS05	79	54	0.7	0.7	280	933.72
NS06	101	76	0.8	0.7	280	1 312.19
NS07	101	76	0.9	0.7	280	1 381.75
NS08	101	76	0.7	0.8	280	1 294.97
NS09	101	76	0.7	0.9	280	1 347.31

下载: 导出CSV

| 显示表格

泡沫铝夹芯管对于爆炸冲击波能量的吸收主要靠内外管的变形与泡沫铝芯层的坍塌塑性变形来实现，定义比吸能E_a为单位质量泡沫铝夹芯管所吸收的能量

$E_{\mathrm{a}}=\frac{E_{\mathrm{t}}}{m}$

(6)

式中：E_t为泡沫铝夹芯管吸收的总能量，在数值模拟中取2 ms时的塑性能量耗散值为结构整体的吸能，m为泡沫铝夹芯管的质量。

2.1 外管直径和壁厚的影响

为研究泡沫铝夹芯管结构中外管直径和壁厚对泡沫铝夹芯管抗爆性能的影响，对2组试样共计6个试件进行了横向爆炸载荷下的动态响应和能量吸收性能研究，其中：第一组试样NS03、NS04和NS05的内外管壁厚和内管直径保持不变，仅改变外管的直径；另一组试样NS01、NS06和NS07的内外管直径和内管壁厚保持不变，仅改变外管的壁厚。加载条件为：乳化炸药质量30 g，起爆位置距离试件跨中40 mm。表4给出了6个不同几何参数试件吸收的总能量、比吸能和内、外管跨中挠度（w_i和w_o）。由表4可知，随着外管直径的减小，试件的总吸能、比吸能和外管跨中挠度减小，内管跨中挠度增大。原因是随着外管直径的减小，泡沫铝夹芯管迎爆面积减小，结构承受的爆炸载荷降低，外管塑性变形减小，因而结构的总吸能减小；而内管跨中挠度增大是由于泡沫铝芯层厚度减小导致芯层密实化使内管吸收了较多的能量。随着外管壁厚的增大，试件整体结构刚度增加，内外管的塑性变形和芯层的压缩量减小，因此试件的总吸能、比吸能和内、外管跨中挠度均减小。

表 4 外管直径和壁厚的影响

Table 4. Effect of outer tube diameter and wall thickness

Test piece	E_t/J	m/g	E_a/(J·g⁻¹)	w_o/mm	w_i/mm
NS03	539.878	1 388.13	0.389	10.78	1.95
NS04	442.605	1 150.28	0.385	9.06	2.21
NS05	354.622	933.72	0.380	8.88	2.87
NS01	542.923	1 242.62	0.437	10.06	2.45
NS06	456.645	1 312.19	0.348	7.61	2.16
NS07	388.472	1 381.75	0.281	5.81	1.76

下载: 导出CSV

| 显示表格

图6给出了不同外管直径和壁厚的泡沫铝夹芯管试件内外管的跨中挠度时程曲线，其中虚线和实线分别代表内、外管跨中挠度。显然，外管跨中挠度大于内管跨中挠度。外管壁厚对外管跨中挠度的影响程度大于外管直径的影响；随着外管直径的增大及外管壁厚的减小，泡沫铝夹芯管芯层压缩量增大，结构吸收的总能量增加。

图 6 外管直径和壁厚对泡沫铝夹芯管内外管跨中挠度的影响

Figure 6. Effects of outer tube diameter and wall thickness on mid-span deflection of inner and outer tubes of aluminum foam-filled sandwich tubes

下载: 全尺寸图片幻灯片

2.2 内管直径和壁厚的影响

为研究泡沫铝夹芯管结构中内管直径和壁厚对泡沫铝夹芯管动态响应的影响，对2组试样共计6个试件进行了横向爆炸载荷下的力学性能研究，其中：第一组试样NS01、NS02和NS03的外管直径和内、外管壁厚保持不变，仅改变内管的直径；另一组试样NS01、NS08和NS09的内外管直径和外管壁厚保持不变，仅改变内管的壁厚。试件的加载条件为：乳化炸药质量30 g，起爆位置距离试件跨中40 mm。表5给出了试件的总吸能、比吸能及内、外管跨中挠度。由表5可知，随着内管直径的减小，试件的比吸能及内管跨中挠度减小，外管跨中挠度增大。内管直径的减小对泡沫铝夹芯管迎爆面积无影响，因而对结构吸收总能量的影响很小，芯层厚度的增大使内管吸能减少，内管跨中挠度减小。随着内管壁厚的增大，试件整体结构刚度增加，内外管的塑性变形和芯层的压缩量减小，因此试件的总吸能、比吸能和内、外管跨中挠度均减小。

表 5 内管直径和壁厚的影响

Table 5. Effect of inner tube diameter and wall thickness

Test piece	E_t/J	m/g	E_a/(J·g⁻¹)	w_o/mm	w_i/mm
NS01	542.923	1 242.62	0.437	10.06	2.45
NS02	546.692	1 326.02	0.412	10.39	2.32
NS03	539.878	1 388.13	0.389	10.78	1.95
NS08	542.470	1 294.97	0.419	9.92	2.16
NS09	542.161	1 347.31	0.402	9.81	1.92

下载: 导出CSV

| 显示表格

图7给出了不同内管直径和不同内管壁厚的泡沫铝夹芯管试件内外管跨中挠度时程曲线，其中虚线和实线分别代表内管和外管跨中挠度时程曲线，外管挠度大于内管挠度。由图7(b)可知，内管壁厚对外管跨中挠度的影响程度小于对内管跨中挠度的影响；对比图6，外管的直径和壁厚对内、外管跨中挠度的影响程度大于内管直径和壁厚的影响。

图 7 内管直径和壁厚对泡沫铝夹芯管内外管跨中挠度的影响

Figure 7. Effects of inner tube diameter and wall thickness on mid-span deflection of inner and outer tubes of aluminum foam-filled sandwich tubes

下载: 全尺寸图片幻灯片

2.3 芯层相对密度的影响

为研究泡沫铝夹芯管结构中泡沫铝夹芯层的相对密度对泡沫铝夹芯管动态性能的影响，对以NS01试件为基础分别填充相对密度为10%、15%和20%的泡沫铝芯层的3个试件进行横向爆炸载荷下的力学性能研究，试件的加载条件为：乳化炸药质量30 g，起爆位置距离试件跨中40 mm。表6给出了不同相对密度芯层的泡沫铝夹芯管结构的总能量、比吸能和内外管跨中挠度。由表6可知，随着泡沫铝芯层相对密度的增大，结构的整体刚度增大，内、外管的塑性变形和芯层压缩量减小。因此，结构的总吸能、比吸能和内、外管跨中挠度均减小。图8给出了不同相对密度芯层的泡沫铝夹芯管试件的内外管跨中挠度时程曲线，其中虚线和实线分别代表内、外管跨中挠度随时间的变化规律，外管挠度大于内管挠度。由图8可知，相对密度越小，泡沫铝芯层的压缩量越大，结构的总吸能和比吸能越大。

表 6 泡沫铝相对密度的影响

Table 6. Effect of relative density of aluminum foam

ρ^*/%	p₀/MPa	E_t/J	m/g	E_a/(J·g⁻¹)	w_o/mm	w_i/mm
10	36.8	614.301	1 116.12	0.550	14.77	4.99
15	36.8	542.923	1 242.62	0.437	10.06	2.45
20	36.8	379.921	1 378.86	0.276	4.82	2.35

下载: 导出CSV

| 显示表格

图 8 泡沫铝相对密度对泡沫铝夹芯管内外管跨中挠度的影响

Figure 8. Effect of relative density of aluminum foam on mid-span deflection of inner and outer tubes of aluminum foam-filled sandwich tubes

下载: 全尺寸图片幻灯片

图9给出了填充不同相对密度泡沫铝夹芯管内、外管的轴向变形轮廓，其中虚线和实线分别代表内、外管变形轮廓。由图9可知，内、外管的变形可分为3个区域，包括中心部位的塑性大变形区域、围绕塑性铰移动的刚性截面和未变形边界区域^[17]。

图 9 不同相对密度泡沫铝夹芯管的轴向变形轮廓

Figure 9. Axial deformation profiles of aluminum foam-filled sandwich tubes with different relative densities

下载: 全尺寸图片幻灯片

2.4 爆炸加载条件的影响

由2.3节可知，横向爆炸载荷下泡沫铝夹芯管NS01试件在所研究的试件中抗爆性能最佳，因此以NS01试件为研究对象，在距离试件跨中37、40和43 mm处施加炸药质量为30 g的爆炸载荷，研究炸药起爆距离对试件动态响应和吸能特性的影响。在距离跨中40 mm处分别施加乳化炸药质量为30、35和40 g的爆炸载荷，研究炸药质量对试件动态响应和吸能特性的影响。表7给出了不同爆炸加载条件下试件NS01的总能量、比吸能和内、外管跨中挠度。由表7可知，随着炸药起爆距离的减小和炸药质量的增加，泡沫铝夹芯管承受的爆炸载荷增大，内外管的塑性变形区域增大。因此，结构的总吸能、比吸能和内、外管跨中挠度均增大。

表 7 爆炸加载条件的影响

Table 7. Effect of blast loading conditions

R/mm	W_s/g	p₀/MPa	E_t/J	E_a/(J·g⁻¹)	w_o/mm	w_i/mm
37	30	34.9	935.594	0.753	14.96	4.34
40	30	27.7	542.923	0.437	10.06	2.45
40	35	32.3	776.974	0.625	12.60	3.62
40	40	36.8	1058.910	0.852	16.20	4.98
43	30	22.4	312.212	0.251	5.92	1.54

下载: 导出CSV

| 显示表格

图10给出了炸药起爆距离与炸药质量对泡沫铝夹芯管试件内外管跨中挠度的影响，其中虚线和实线分别代表内、外管跨中挠度时程曲线，外管的挠度大于内管的挠度。由图10可知，炸药起爆距离和炸药药量对泡沫铝夹芯管的内、外管跨中挠度影响显著。

图 10 炸药起爆距离和炸药质量对泡沫铝夹芯管内外管跨中挠度的影响

Figure 10. Effect of explosive initiation distance on mid-span deflection of inner and outer tubes of aluminum foam-filled sandwich tubes

下载: 全尺寸图片幻灯片

泡沫铝夹芯管各组成部分的吸能占比如表8所示。由表8可知，外管吸收的能量最多，吸能百分比均超过了50%，其次是泡沫铝芯层，吸能最少的是内管。原因是外管作为迎爆面，首先承受爆炸载荷，之后外管压缩芯层，芯层部分密实化后与外管共同压缩内管变形。随着外管壁厚的减小，外管的吸能百分比减小，泡沫铝芯层的吸能百分比增大；内管壁厚对各部分吸能百分比的影响较小；当夹芯管外管直径或内管直径一定时，随着芯层厚度的增加，内管的吸能百分比减小；芯层的相对密度对泡沫铝夹芯管的吸能百分比影响显著，随着芯层相对密度的减小，芯层和内管的吸能百分比减小，外管的吸能百分比增加。

表 8 泡沫铝夹芯管组成部分的吸能百分比

Table 8. Energy absorption percentage of aluminum foam-filled sandwich tube components

Test piece	$\rho ^*$ /%	R/mm	W_s/g	p₀/MPa	Energy absorption percentage/%
Test piece	$\rho ^*$ /%	R/mm	W_s/g	p₀/MPa	Outer tube	Foam core	Inner tube
NS01	15	40	30	27.7	59.9	32.0	8.1
NS02	15	40	30	27.7	62.6	30.5	6.9
NS03	15	40	30	27.7	61.8	33.7	4.5
NS04	15	40	30	27.7	62.4	30.7	6.9
NS05	15	40	30	27.7	59.7	32.3	8.0
NS06	15	40	30	27.7	64.0	28.1	7.9
NS07	15	40	30	27.7	66.3	25.6	8.2
NS08	15	40	30	27.7	59.5	32.4	8.0
NS09	15	40	30	27.7	59.2	32.9	7.9
NS01	10	40	30	27.7	79.8	17.5	2.7
NS01	20	40	30	27.7	52.0	32.2	15.8
NS01	15	37	30	34.9	61.7	31.0	7.3
NS01	15	43	30	22.4	61.5	29.6	8.9
NS01	15	40	35	32.3	62.0	30.5	7.5
NS01	15	40	40	36.8	61.5	31.2	7.3

下载: 导出CSV

| 显示表格

3. 抗爆性能的多目标优化

3.1 抗爆性能的定义及优化方法

由参数分析可知，内、外管壁厚和泡沫铝芯层的相对密度对泡沫铝夹芯管的抗爆性能具有一定影响，期望芯层比吸能最大化的同时，尽可能将外管跨中挠度控制在合理的范围内。因此，将最大化芯层比吸能E_CSA和最小化外管跨中挠度w_o作为优化目标，通过多目标优化方法研究泡沫铝夹芯管的抗爆性能。以泡沫铝夹芯管NS01内、外管壁厚以及芯层的相对密度作为设计参数（壁厚单位：mm），则多目标优化问题表示为

$\left\{ \begin{gathered} {{f}}\left( x \right){{ = {\rm{min}}}}\left( {{1 \mathord{\left/ {\vphantom {1 {{E_{{\rm{CSA}}}}(H,h,{\rho ^*})}}} \right. } {{E_{{\rm{CSA}}}}(H,h,{\rho ^*})}},{w_{\text{o}}}(H,h,{\rho ^*})} \right) \\ {\rm{s.t.}}\;{\text{ }}10{\text{%}} \leqslant {\rho ^*} \leqslant 20{\text{%}} {\text{ }} \\ \;\;\;\;0.7 \leqslant H \leqslant 1.78 \\ \;\;\;\;0.7 \leqslant h \leqslant 1.78 \\ \end{gathered} \right.$

(7)

多目标优化方案如图11所示。首先，采用正交试验设计方法，在设计变量的范围内选取样本点；然后，采用响应面法（RSM），根据有限元模拟得到样本点数据建立E_CSA和w_o的代理模型，并对代理模型进行评估；最后，基于非支配排序遗传算法（NSGA-Ⅱ）对参数进行优化，得到多目标优化问题的Pareto前沿解集。

图 11 优化方法流程

Figure 11. Flowchart of the optimization method

下载: 全尺寸图片幻灯片

采用正交试验设计的方法选取试验点^[20]，将外管与内管壁厚、泡沫铝芯层的相对密度作为3个影响因子，内、外管壁厚采用4个水平，内、外管壁厚水平分别采用0.70、1.06、1.42、1.78 mm，相对密度采用3个水平，相对密度水平采用10%、15%和20%，按照正交试验设计方法选取48个样本点进行数值计算，乳化炸药的质量为30 g，起爆位置距离试件跨中40 mm，分别得到不同样本点的外管跨中挠度和泡沫铝芯层的比吸能，并以此建立响应面模型。

响应面法具有简单易行、效果良好的特点，在建立高度非线性的代理模型中得到了广泛应用，使用三次多项式构建代理模型，得到48个样本点对应的E_CSA和w_o，通过数据拟合和最小二乘法确定响应面代理模型的系数。由设计变量H、h、 $\,\rho^*$ 构造的E_CSA和w_o三次多项式响应面函数为

$\begin{split} 1/{E_{\rm{CSA}}} =& - 6.477\,3+49.519\,7H - 23.364\,1h+122.008\,0{\rho^ {*} } - 8.864\,1{H^2} - 4.141\,7{h^2} - 481.066\,6{\rho ^ {*2} }+18.808\,4Hh - \\ &739.643\,0H{\rho ^ {*} }+301.309\,9h{\rho ^ {*} }+0.389\,1{H^3} - 0.228\,9{h^3} - 224.429\,9{\rho ^ {*3} }- 4.478\,3{H^2}h+116.427\,8{H^2}{\rho ^ {*} } + \\ & 2.412\,4H{h^2}+21.482\,4{h^2}{\rho^ {*} }+2\,533.473\,6H{\rho ^ {*2} }- 847.899\,4h{\rho ^ {*2} } - 127.966\,0Hh{\rho ^ * } \\[-15pt]\end{split}$

(8)

$\begin{split} {w_{\text{o}}} = & 74.447\,9 - 95.263\,9H - 19.720\,6h - 162.469\,6{\rho ^ * }+44.208\,1{H^2}+3.002\,8{h^2} - 324.652\,3{\rho ^ {*2} }+11.999\,2Hh+ \\ & 151.621\,9H{\rho ^ * }+97.555\,6h{\rho ^ * } - 6.028\,2{H^3} - 0.044\,7{h^3} - 135.309\,7{\rho ^ {*3} }- 2.002\,2{H^2}h - 75.762\,0{H^2}{\rho ^ * } - \\ &1.305\,7H{h^2} - 4.677\,9{h^2}{\rho ^ * }+430.972\,2H{\rho ^ {*2} } - 201.805\,6h{\rho ^ {*2} } - 18.248\,5Hh{\rho ^ * } \end{split}$

(9)

为了评估代理模型的准确性，采用拟合优度R²和相对误差的绝对值 $\left| \delta \right|$ 两种评价指标对多项式响应面进行误差分析，这两种指标的计算方法为

${R^2} = 1 - \frac{{\displaystyle\sum\limits_{i=1}^n {{{({y_i} - {{\hat y}_i})}^2}} }}{{\displaystyle\sum\limits_{i=1}^n {{{({y_i} - \bar y)}^2}} }}$

(10)

$\left| \delta \right| = \frac{{\displaystyle\sum\limits_{i=1}^n {\left| {{y_i} - {{\hat y}_i}} \right|} }}{{\displaystyle\sum\limits_{i=1}^n {\left| {{y_i} - \bar y} \right|} }}$

(11)

式中： ${y_i}$ 和 ${\hat y_i}$ 分别为第i个样本点的有限元模拟和多项式响应面的响应值； $\bar y$ 为 ${y_i}$ 的平均值；n为样本点的个数；R²为评价响应面函数与数值计算数据拟合准确度的指标， $\left| \delta \right|$ 用于评估模型的整体精度， ${R^2} > 0.9$ 、 $\left| \delta \right|$ <0.2表示代理模型有足够高的精确度。结果显示：E_CSA和w_o的R²分别为0.992 6和0.993 8，均大于0.99，接近1，E_CSA和w_o的 $\left| \delta \right|$ 分别为0.096 9和0.084 7，均小于0.1，表明建立的多项式代理模型精度满足要求。

3.2 优化结果分析

最终得到的Pareto前沿如图12所示，空心点为NSGA-Ⅱ得到的数据。当E_CSA增大时，w_o也增大，两个优化目标相互冲突。由图12可知，所有的解集中w_o都小于初始设计点。表9给出了设计点的结构参数，并将有限元模拟结果与响应面函数得到的结果进行比较，误差均小于9.5%，验证了代理模型的准确性。当仅要求芯层比吸能最大化时，最佳参数设计点为Op1，与初始设计点相比，外管跨中挠度减小了28.9%，芯层比吸能提高了11.7%。当仅要求外管跨中挠度最小化时，最佳参数设计点为Op2，与初始设计点相比，外管跨中挠度减小了95.7%，芯层比吸能降低了87.5%。为使E_CSA和w_o达到平衡，采用最小距离选择法设计Op3，最小距离选择法的公式为

图 12 Pareto前沿

Figure 12. Pareto front

下载: 全尺寸图片幻灯片

表 9 最佳设计点

Table 9. Optimum design point

Design point	H/mm	h/mm	$\,\rho ^ *$ /%	w_o			E_CSA
Design point	H/mm	h/mm	$\,\rho ^ *$ /%	FEM/mm	RSM/mm	Error/%	FEM/(J·g⁻¹)	RSM/(J·g⁻¹)	Error/%
Op1	0.70	0.76	13.7	0.457	0.505	−9.5	10.50	11.14	−5.7
Op2	1.78	1.75	18.8	0.051	0.056	−8.9	0.64	0.59	8.5
Op3	0.84	1.76	20.0	0.194	0.205	−5.4	2.85	2.83	0.7

下载: 导出CSV

| 显示表格

$\min S = \sqrt {{{\left[ {\frac{{\left( {1/{E_{\rm{CSA}}}\left( {H,h,{\rho ^*}} \right)} \right) - {{\left( {1/{E_{\rm{CSA}}}\left( {H,h,{\rho ^*}} \right)} \right)}_{\min }}}}{{{{\left( {1/{E_{\rm{CSA}}}\left( {H,h,{\rho ^*}} \right)} \right)}_{\max }} - {{\left( {1/{E_{\rm{CSA}}}\left( {H,h,{\rho ^*}} \right)} \right)}_{\min }}}}} \right]}^2}+{{\left[ {\frac{{{w_{\rm{o}}}\left( {H,h,{\rho ^*}} \right) - {w_{\rm{o}}}{{\left( {H,h,{\rho ^*}} \right)}_{\min }}}}{{{w_{\rm{o}}}{{\left( {H,h,{\rho ^*}} \right)}_{\max }} - {w_{\rm{o}}}{{\left( {H,h,{\rho ^*}} \right)}_{\min }}}}} \right]}^2}}$

(12)

式中：S为Pareto解集中每个非劣解到理想最优解（ideal point）的距离， $\left( {{1 \mathord{\left/ {\vphantom {1 {{E_{\rm{CSA}}}\left( {H,h,{\rho ^*}} \right)}}} \right. } {{E_{{\rm{CSA}}}}\left( {H,h,{\rho ^*}} \right)}}} \right)$ 和 ${w_{\rm{o}}}\left( {H,h,{\rho ^*}} \right)$ 为Pareto解集中每个非劣解对应芯层比吸能的倒数和外管跨中挠度的值， ${\left( {{1 \mathord{\left/ {\vphantom {1 {{E_{{\rm{CSA}}}}\left( {H,h,{\rho ^*}} \right)}}} \right. } {{E_{{\rm{CSA}}}}\left( {H,h,{\rho ^*}} \right)}}} \right)_{\max }}$ 和 ${\left( {{1 \mathord{\left/ {\vphantom {1 {{E_{{\rm{CSA}}}}\left( {H,h,{\rho ^*}} \right)}}} \right. } {{E_{{\rm{CSA}}}}\left( {H,h,{\rho ^*}} \right)}}} \right)_{\min }}$ 分别为所有Pareto解集中芯层比吸能的倒数的最大值和最小值， ${w_{\rm{o}}}{\left( {H,h,{\rho ^*}} \right)_{\max }}$ 和 ${w_{\rm{o}}}{\left( {H,h,{\rho ^*}} \right)_{\min }}$ 分别为所有Pareto解集中外管跨中挠度的最大值和最小值。最佳参数设计点为Op3，与初始设计点相比，外管跨中挠度减小了80.7%，芯层比吸能降低了52.6%。

4. 结　论

采用数值模拟方法研究了泡沫铝夹芯管在横向爆炸载荷下的动态响应，对比分析了泡沫铝夹芯管的几何参数、芯层的相对密度和爆炸加载条件对泡沫铝夹芯管的总吸能、比吸能和内、外管跨中挠度的影响。以泡沫铝夹芯管内管壁厚、外管壁厚和芯层相对密度为设计变量进行多目标优化，以提高泡沫铝夹芯管结构的抗爆性能，得到以下主要结论。

(1) 泡沫铝夹芯管结构的几何参数和芯层的相对密度对其抗爆性能存在显著影响，相对于内管的直径和壁厚，外管直径和壁厚对内外管跨中挠度的影响更显著。随着外管直径的增大和内管直径的减小，泡沫铝夹芯管内管跨中挠度减小，外管跨中挠度增大，比吸能随内外管直径的减小而减小；随着内管和外管壁厚的增大，泡沫铝夹芯管的比吸能和内外管跨中挠度减小。泡沫铝芯层的相对密度越大，泡沫铝夹芯管的比吸能和内、外管跨中挠度越小。

(2) 随着炸药质量的增加和起爆距离的减小，泡沫铝夹芯管的变形区域增大，比吸能和内、外管跨中挠度增大。分析了不同工况下泡沫铝夹芯管各组成部分的吸能情况，其中外管吸收的能量最多，均超过了50%，其次是泡沫铝芯层，吸能最少的是内管；内管壁厚对各部分吸能百分比的影响较小；芯层的厚度和相对密度对泡沫铝夹芯管的吸能百分比影响显著。

(3) 以内管壁厚、外管壁厚和芯层相对密度为设计变量，以芯层的比吸能和外管跨中挠度为优化指标，采用响应面法和非支配排序遗传算法可以实现泡沫铝夹芯管抗爆性多目标优化设计，得到了Pareto前沿图，发现外管跨中挠度取最小值时，对应的芯层比吸能最小。因此，当确定了外管跨中挠度和芯层比吸能的偏好时，可确定泡沫铝夹芯管结构的内外管壁厚和芯层的相对密度。

图 1 实验样品组装图

Figure 1. Experimental setup for electrical conductivity measurements

下载: 全尺寸图片幻灯片

图 2 (a)纯橄榄石和(b)含15% FeS的橄榄石的背散射图

Figure 2. Backscattered electron images of (a) pure olivine and (b)15% FeS–bearing olivine

下载: 全尺寸图片幻灯片

图 3 不同温度压力条件下纯橄榄石和含15% FeS的橄榄石复阻抗谱的对比

Figure 3. Comparison of complex impedance spectra of pure olivine and 15% FeS–bearing olivine under different temperature and pressure conditions

下载: 全尺寸图片幻灯片

图 4 纯的和含15% FeS橄榄石的电导率对数与温度倒数的关系

Figure 4. Logarithm of electrical conductivity reciprocal temperature for pure olivine and 15% FeS–bearing olivine at different pressures

下载: 全尺寸图片幻灯片

图 5 活化焓与压力的关系

Figure 5. Activation enthalpy versus pressure

下载: 全尺寸图片幻灯片

表 1 样品电导率的Arrhenius关系拟合参数

Table 1. Fitted parameters of Arrhenius relation for the electrical conductivity of samples

Sample	p/GPa	T/K	lg $\sigma_0$	$\sigma_0$ /(S∙m⁻¹)	${\Delta H}/{\rm{eV}}$	${\gamma^2}$
	1	723–1273	2.53	338.84	1.18	0.9913
Olivine	2	723–1273	1.96	91.20	1.11	0.9781
	3	723–1273	1.05	11.22	1.03	0.9867
	1	723–1273	0.47	2.95	0.20	0.9953
FeS-bearing olivine	2	723–1173	1.73	53.70	0.29	0.9802
	3	723–1073	2.78	602.56	0.39	0.9845
Note: ${\gamma}^2$ —adjust R-square.

下载: 导出CSV

参考文献(19)

[1]	RINGWOOD A E. Composition and petrology of the Earth’s mantle [M]. New York: McGraw-Hill, 1975.
[2]	IRIFUNE T, RINGWOOD A E. Phase-transformations in a harzburgite composition to 26 GPa: implication for dynamical behavior of the subducting slab [J]. Earth and Planetary Science Letters, 1987, 86(2/3/4): 365–376.
[3]	YOSHINO T, WALTER M J, KATSURA T. Core formation in planetesimals triggered by permeable flow [J]. Nature, 2003, 422(6928): 154–157. doi: 10.1038/nature01459
[4]	YOSHINO T, WALTER M J, KATSURA T. Connectivity of molten Fe alloy in peridotite based on in situ electrical conductivity measurements:implications for core formation in terrestrial planets [J]. Earth and Planetary Science Letters, 2004, 222(2): 625–643. doi: 10.1016/j.jpgl.2004.03.010
[5]	WANG D J, KARATO S I, JIANG Z T. An experimental study of the influence of graphite on the electrical conductivity of olivine aggregates [J]. Geophysical Research Letters, 2013, 40(10): 2028–2032. doi: 10.1002/grl.50471
[6]	BAGDASSAROV N, GOLABEK G J, SOLFERINO G, et al. Constraints on the Fe-S melt connectivity in mantle silicates from electrical impedance measurements [J]. Earth and Planetary Science Letters, 2009, 177(3/4): 139–146.
[7]	WATSON H C, ROBERTS J J, TYBURCZY J A. Effect of conductive impurities on electrical conductivity in polycrystalline olivine [J]. Geophysical Research Letters, 2010, 37: L02302.
[8]	WATSON H C, ROBERTS J J. Connectivity of core forming melts: experimental constraints from electrical conductivity and X-ray tomography [J]. Physics of the Earth and Planetary Interiors, 2011, 186(3/4): 172–182.
[9]	ZHANG Z, POMMIER A. Electrical investigation of metal-olivine systems and application to the deep interior of mercury [J]. Journal of Geophysical Research–Planets, 2017, 122(12): 2702–2718. doi: 10.1002/2017JE005390
[10]	OMURA K, KURITA K, KUMAZAWA M. Experimental study of pressure dependence of electrical conductivity of olivine at high temperatures [J]. Physics of the Earth and Planetary Interiors, 1989, 57(3/4): 291–303.
[11]	XU Y S, SHANKLAND T J, DUBA A G. Pressure effect on electrical conductivity of mantle olivine [J]. Physics of the Earth and Planetary Interiors, 2000, 118(1/2): 149–161.
[12]	DAI L D, KARATO S I. The effect of pressure on the electrical conductivity of olivine under the hydrogen-rich conditions [J]. Physics of the Earth and Planetary Interiors, 2014, 232: 51–56. doi: 10.1016/j.pepi.2014.03.010
[13]	DAI L D, HU H Y, LI H P, et al. Influence of temperature, pressure, and oxygen fugacity on the electrical conductivity of dry eclogite and geophysical implications [J]. Geochemistry Geophysics Geosystems, 2016, 17(6): 2394–2407. doi: 10.1002/2016GC006282
[14]	HU H Y, DAI L D, LI H P, et al. Influence of dehydration on the electrical conductivity of epidote and implications for high–conductivity anomalies in subduction zones [J]. Journal of Geophysical Research–Solid Earth, 2017, 122(4): 2751–2762. doi: 10.1002/2016JB013767
[15]	SHI C Y, ZHANG L, YANG W G, et al. Formation of an interconnected network of iron melt at Earth’s lower mantle conditions [J]. Nature Geosciences, 2013, 6(11): 971–975. doi: 10.1038/ngeo1956
[16]	ROBERTS J J, TYBURCZY J A. Impedance spectroscopy of single and polycrystalline olivine: evidence for grain boundary transport [J]. Physics and Chemistry of Miners, 1993, 20(1): 19–26.
[17]	HIRAGA T, ANDERSON I M, KOHLSTEDT D L. Grain boundaries as reservoirs for incompatible elements in the Earth’s mantle [J]. Nature, 2004, 427(6976): 699–703. doi: 10.1038/nature02259
[18]	DAI L D, LI H P, HU H Y, et al. Experimental study of grain boundary electrical conductivities of dry synthetic peridotite under high temperature, high-pressure, and different oxygen fugacity conditions [J]. Journal of Geophysical Research–Solid Earth, 2008, 113(B12): B12211. doi: 10.1029/2008JB005820
[19]	TERASAKI H, FROST D J, RUBIE D C, et al. The effect of oxygen and sulphur on the dihedral angle between Fe–O–S melt and silicate minerals at high pressure: implications for Martian core formation [J]. Earth and Planetary Science Letters, 2005, 232(3/4): 379–392.