地球内核超离子态铁合金及其效应

何宇; 孙士川; 李和平

doi:10.11858/gywlxb.20240707

地球内核超离子态铁合金及其效应

doi: 10.11858/gywlxb.20240707

中国科学院地球化学研究所, 中国科学院地球内部物质高温高压重点实验室, 贵州贵阳　550081

基金项目: 国家自然科学基金（42074104）；中国科学院青年交叉团队项目（JCTD-2022-16）

详细信息

作者简介:
何　宇（1985－），男，博士，研究员，主要从事高温高压下地球内部物质性质的计算模拟研究.E-mail：heyu@mail.gyig.ac.cn

中图分类号: O521.2; P313
计量
- 文章访问数: 298
- HTML全文浏览量: 174
- PDF下载量: 51
出版历程
- 收稿日期: 2024-01-09
- 修回日期: 2024-02-29
- 网络出版日期: 2024-05-10
- 刊出日期: 2024-06-03

Superionic Iron Alloys in Earth’s Inner Core and Their Effects

Key Laboratory of High-Temperature and High-Pressure Study of the Earth’s Interior, Institute of Geochemistry, Chinese Academy of Sciences, Guiyang 550081, Guizhou, China

摘要

摘要: 超离子态介于固态与液态之间，被认为广泛存在于地球和行星内部。计算研究发现，在地球内核温度压力条件下，铁-氢、铁-碳、铁-氧合金处于超离子态，表现为氢、碳、氧等元素在固态铁合金中像液体一样快速流动。流动的轻元素导致铁合金软化及地震波速降低，与地球物理观测到的内核密度亏损和低剪切波速的特征一致。内核超离子态铁-氢合金可以与地磁场发生相互作用，在偶极地磁场的驱动下形成定向排列组构，从而解释了内核的各向异性结构成因。内核超离子态铁-轻元素合金的发现更新了人们对内核物态的认知，对掌握地球内核的结构、组成和演化以及内核结构与地球磁场的关系等具有重要意义。
- 超离子态 /
- 地球内核 /
- 铁合金 /
- 各向异性 /
- 第一性原理计算
Abstract: Under the conditions of high temperature and high pressure, a series of materials transform into superionic states, which fall between the solid and liquid states and are widely believed to exist in the interior of Earth and exoplanets. Computational research has found that under the temperature and pressure of the Earth’s inner core, iron-hydrogen, iron-carbon, and iron-oxygen alloys transform to superionic states, manifested as elements such as hydrogen, carbon, and oxygen flowing rapidly like liquids in solid iron alloys. The flowing light elements cause softening of Fe alloys and a decrease in seismic wave velocities, explaining the characteristics of core density deficient and low shear wave velocity observed in geophysics. The superionic iron-hydrogen alloy in the core can interact with the geomagnetic field, forming a lattice preferred orientation fiber driven by a dipole geomagnetic field, explaining the origin of the anisotropic structure in the inner core. The discovery of superionic iron-light-element alloys in the inner core has updated our understanding of the state of the inner core, and is of great significance for understanding the structure, composition, and evolution of Earth’s inner core, as well as the relationship between the inner core structure and the Earth’s magnetic field.
- superionic state /
- Earth’s inner core /
- Fe alloy /
- anisotropic structure /
- first-principles calculations

HTML全文

软质高分子聚合物材料因生产成本低、质量轻、具有良好的抗冲击性能等优点在军工、汽车、航空、电子^[1-4]等行业得到广泛应用，如用于汽车碰撞实验的假人皮肤材料^[5]、航空航天仿生材料^[6]等。在上述领域，构件/材料常受动态冲击载荷作用，明确其动态力学性能，构建描述其黏弹性力学特征的本构模型，可为软质高分子聚合物材料的设计开发、性能优化及可靠性分析提供理论模型和方法支撑。

学者们曾利用分离式霍普金森压杆（Split Hopkinson pressure bar, SHPB）装置对PVC弹性体等软质高分子聚合物材料开展了动态压缩力学性能研究^[7-8]。与SHPB相比，分离式霍普金森拉杆（Split Hopkinson tensile bar, SHTB）装置采用拉伸加载形式，便于开展材料的动态冲击损伤、变形、断裂等行为研究。近年来，研究人员已对丁氰橡胶^[9]、聚脲^[10]、半晶聚合物高密度聚乙烯（HDPE）^[11]、猪皮^[12]、聚氨酯弹性体^[13]等低阻抗软材料进行了静、动态拉伸力学性能测试。同时，自动网格法应变测量^[14]、脉冲整形^[15]等技术以及高/低温环境箱^[16]、扫描电镜^[17-18]、高速摄像机^[19]等装置也被用于软材料的拉伸变形、裂纹萌生及扩展机理分析。通过静、动态拉伸实验发现：软质高分子聚合物材料的拉伸强度会随着应变率对数的增大而线性增加，且不同分子链结构的软质高分子聚合物在同一高拉伸应变率下有不同的拉伸强度和拉伸应变^[20]。

但是，软质高分子聚合物材料动态拉伸实验过程中仍存在试样连接方式难以确定、胶黏剂黏接强度低、加载应变率不恒定等问题。为揭示PVC弹性体材料的静、动态拉伸力学性能，利用Instron-5943万能材料试验机和改进型SHTB实验装置对PVC弹性体进行静、动态直接拉伸实验，分析应变率对PVC弹性体材料静、动态拉伸力学性能的影响规律。

1. 实验材料及测试方法

1.1 准静态拉伸实验

实验材料为PVC弹性体。准静态拉伸试样根据GB/T 528–2009^[21]中的3型试样标准设计由模具冲压而成，具体尺寸见图1，试样呈哑铃状，厚度为2 mm。由于弹性体材料动态拉伸实验试样的设计尚无统一标准，因此本研究中动态拉伸实验试样将依据动态拉伸实验的具体效果进行设计并制备，详见2.3.1节。

图 1 准静态拉伸试样尺寸（单位：mm）

Figure 1. Dimension of quasi-static tensile specimen (Unit：mm)

下载: 全尺寸图片幻灯片

PVC弹性体的准静态拉伸力学性能实验是在Instron-5943万能材料试验机上进行的。实验时，将试验机的加载速率设置为96 mm/min，使其应变率为0.1 s⁻¹。拉伸过程中的载荷（F）和位移（ΔL）数据通过传感器输入至计算机中，利用式（1）、式（2）将其转化为应力、应变数据

$\sigma {\rm{ = }}\frac{F}{A_{\rm s}}$

(1)

$\varepsilon {\rm{ = }}\frac{{\Delta L}}{{{L\rm_s}}}$

(2)

式中： $\sigma$ 为应力，A_s为试样截面积， $\varepsilon$ 为应变，L_s为试样长度。

1.2 动态拉伸实验

采用合肥姜水动态力学实验技术有限公司第一实验室的分离式霍普金森拉杆（SHTB）实验系统开展聚氯乙烯（PVC）弹性体材料在高应变率下的动态拉伸实验，其装置原理及实物见图2。

图 2 分离式霍普金森拉杆装置示意图

Figure 2. Schematic of split Hopkinson tensile bar setup

下载: 全尺寸图片幻灯片

调节储气瓶阀向储气室充气，待储气室内气压到达预定值时，关闭储气瓶阀，触发发射按钮，使储气室内高压气体流入发射通道，推动管状套筒子弹在入射杆上加速至一定速度，并以该速度撞击入射杆顶端的法兰盘，在入射杆中产生一个压缩应力波，其持续时间由子弹长度决定。套筒子弹与法兰盘撞击瞬间形成的压缩应力波分为两部分：一部分压缩应力波被吸收杆捕获，最终被能量吸收器（阻尼器）吸收；另一部分压缩应力波在入射杆自由端面处反射形成入射杆中的拉伸应力波。当该波传递至试样端面时，由于杆系与试样间的波阻抗不匹配，一部分拉伸应力波会被反射回入射杆形成反射波，另一部分则会穿过试样并传递至透射杆形成透射波。

根据一维弹性应力波的传播理论，由入射杆和透射杆上的半导体应变片分别获取入射波 ${\varepsilon \rm_i}\left( t \right)$ 、反射波 ${\varepsilon \rm_r}\left( t \right)$ 和透射波 ${\varepsilon \rm_t}\left( t \right)$ 的应变脉冲信号，可计算出试样两端面的力 ${F_1}$ 、 ${F_2}$ 及位移 ${u_1}$ 、 ${u_2}$

${F_1} = {E_0}{A_0}\left[ {{\varepsilon \rm_i}\left( t \right) + {\varepsilon \rm_r}\left( t \right)} \right]$

(3)

${F_2} = {E_0}{A_0}{\varepsilon \rm_t}\left( t \right)$

(4)

${u_1} = {C_0}\int {\left[ {{\varepsilon \rm_i}\left( t \right) - {\varepsilon \rm_r}\left( t \right)} \right]} {\rm{d}}t$

(5)

${u_2} = {C_0}\int {{\varepsilon \rm_t}} \left( t \right){\rm{d}}t$

(6)

当试样内部满足应力均匀性假设^[22]时，材料的应力、应变、应变率具有以下关系

${\sigma \rm_s} = \frac{{{F_1} + {F_2}}}{{2{A\rm_s}}}$

(7)

${\varepsilon \rm_s} = \frac{{{u_1} - {u_2}}}{{2{L\rm_s}}}$

(8)

${\dot \varepsilon \rm_s} = \frac{{{\rm{d}}\varepsilon }}{{{\rm{d}}t}}$

(9)

由于二波法和三波法比较时，三波法的可信度较高，更能反映真实的测试结果，故使用三波法处理实验所得数据。整理式（3）～式（9）可得数据处理公式

${\dot \varepsilon \rm_s} = \frac{{{C_0}}}{{{L\rm_s}}}\left[ {{\varepsilon \rm_i}(t) - {\varepsilon \rm_r}(t) - {\varepsilon \rm_t}(t)} \right]$

(10)

${\varepsilon \rm_s} = \frac{{{C_0}}}{{{L_{\rm{s}}}}}\int_0^t {\left[ {{\varepsilon \rm_i}(t) - {\varepsilon \rm_r}(t) - {\varepsilon \rm_t}(t)} \right]} {\rm{d}}t$

(11)

${\sigma \rm_s} = \frac{{{A_0}{E_0}}}{{2{A\rm_s}}}\left[ {{\varepsilon \rm_i}\left( t \right) + {\varepsilon \rm_r}\left( t \right) + {\varepsilon \rm_t}\left( t \right)} \right]$

(12)

式中： ${\dot \varepsilon \rm_s}$ 、 ${\varepsilon \rm_s}$ 和 ${\sigma \rm_s}$ 分别为动态拉伸实验中试样的应变率、应变和应力， ${\varepsilon \rm_i}\left( t \right)$ 、 ${\varepsilon \rm_r}\left( t \right)$ 和 ${\varepsilon \rm_t}\left( t \right)$ 分别为应变片所获取的电压信号经转换之后的入射应变、反射应变和透射应变， ${C_0}$ 、 ${A_0}$ 和 ${E_0}$ 分别为杆的弹性波波速、横截面积和弹性模量， $t$ 为应力波在试样内的传播时间。

1.3 动态拉伸实验关键问题及解决方法

SHTB实验装置的拉杆均为直径20 mm的实心铝杆，杆系总长度为9500 mm，其中，入射杆长度为5000 mm，透射杆长度为3000 mm，吸收杆长度为1500 mm。法兰盘安装于入射杆上，管状铝质套筒子弹长度为600 mm。软质高分子聚合物材料的SHTB实验与SHPB实验不同，实验过程存在动态拉伸试样连接方式及胶黏剂优选、恒应变率加载与入射波完全卸载等问题，可联合波形分析和高速摄像等方法加以解决。

1.3.1 试样连接方式的确定及胶黏剂优选

目前，SHTB实验中试样的连接方式主要有3种：螺纹连接、挂接以及黏接。由于PVC弹性体材料自身特性不宜采用螺纹连接，因而本研究首先尝试了较常用的挂接方式（使用特殊夹具夹持固定试样）。

挂接夹具和哑铃状试样尺寸如图3所示，由于聚氯乙烯（PVC）弹性体材料的波速较低，难以达到应力均匀，故试样不可过长；同时若试样长度过短，试样两端夹持部分会影响实验结果的准确性。根据橡胶材料测试经验^[23]，本研究动态拉伸试样有效拉伸长度设置为2 mm。

图 3 挂接夹具及哑铃状试样尺寸（单位：mm）

Figure 3. Dimensions of clamps and dumbbell-shaped specimens (Unit：mm)

下载: 全尺寸图片幻灯片

挂接所得波形见图4。由图4可知：入射脉冲为矩形方波，且在100 μs内便可达到最大应变；透射脉冲较宽且透射信号在0～0.2 μs之间有明显下降点，但始终未下降至水平基线位置。

图 4 挂接波形

Figure 4. Mounting waveform

下载: 全尺寸图片幻灯片

通过分析发现，出现这一现象主要有两方面原因：一是挂接中使用的夹具对波形传播造成了干扰；二是透射杆长度过短，使得通过应变片的透射波沿透射杆返回，与尚未完全通过应变片的透射波在应变片处发生了叠加。

为解决上述问题，对原有SHTB实验装置进行了改进。改进后的SHTB杆系总长度为11500 mm，其中，入射杆长度为5000 mm，透射杆长度为5000 mm，吸收杆长度为1500 mm。试样的连接方式改为黏接（将试样用合适强度的胶黏剂黏接于入射杆与透射杆之间），试样尺寸为 ${\varnothing}$ 20 mm× 2 mm。

由于胶黏剂完全固化后才可达到最大黏接强度，为加快实验速度，首先选用了强度高、完全固化时间短（约10 min）的502瞬间强力胶进行了实验，试样黏接示意图及加载波形见图5。

图 5 改进后的SHTB黏接试样及加载波形

Figure 5. Modified SHTB bonded sample and loading waveform

下载: 全尺寸图片幻灯片

由图5可知，入射波仍为矩形方波，且可在100 μs内达到最大应变；透射脉冲在极短时间内上升至最大应变后迅速下降且回到基线位置。为明确透射波上升至最大应变后立即卸载的原因，采用高速摄像机（V1212，深圳约克科技有限公司）记录了动态拉伸实验过程，发现波传递至试样与入射杆黏接面的瞬间，502瞬间强力胶发生强烈抖动，导致试样与入射杆黏接面脱胶，致使透射波上升至最大应变后立即卸载，高速摄像照片见图6。

图 6 502瞬间强力胶黏接试样拉伸对比

Figure 6. Tensile comparisons of 502 superglue bonded samples

下载: 全尺寸图片幻灯片

为解决脱胶问题，对比了多种类型的胶黏剂（如YH-818专用橡胶胶水、固特灵401胶水、环氧树脂AB胶、JL-330橡胶专用胶等）并进行了胶黏剂优选。研究发现，环氧树脂AB胶和JL-330橡胶专用胶可以将试样牢固地黏接于入射杆与透射杆之间。利用SHTB装置对这两种胶黏剂黏接的试样进行了动态拉伸实验，以确定这两种胶黏剂的适用性，实验所得波形见图7。

图 7 不同胶黏剂黏接典型波形

Figure 7. Typical waveforms of different glues bonded samples

下载: 全尺寸图片幻灯片

由图7可知，采用两种胶黏剂所得入射脉冲均为矩形方波，且可在100 μs内达到最大应变，其中，环氧树脂AB胶所得透射脉冲上升至最大应变后迅速下降并可回到基线，JL-330橡胶专用胶所得透射波加载时间变长（约在200 μs左右），且达到最大应变后缓慢下降，说明JL-330橡胶专用胶可在一定程度上缓解透射波上升至最大应变后立即卸载的问题，故后续实验将尝试用JL-330橡胶专用胶进行聚氯乙烯弹性体动态拉伸实验。

1.3.2 恒应变率加载

当SHTB实验采用直接加载时，入射波为方波且在极短时间内即可达到最大应变，难以保证整个加载过程为恒应变率加载。为确保SHTB实验结果的准确性，在法兰盘端面（套筒子弹撞击端面）上放置脉冲整形器，以获得所需的入射波特征（通常保证入射波上升沿所对应的加载时间大于加载脉冲在试样中3个往复所需时间）并实现恒应变率加载。在本SHTB实验中，脉冲整形器材料采用铜版纸，并加工成外径为28.4 mm、内径为20 mm的环形片，涂抹少量凡士林将其黏附在法兰盘的端面上，脉冲整形器的安装见图8。

图 8 脉冲整形器安装示意图

Figure 8. Schematic of pulse shaper installation

下载: 全尺寸图片幻灯片

试样中应力波波速约为70 m/s，达到基本均匀需要约240 μs才能使波在试样中传播3个来回。经过多次实验发现，实验中加入整形器后所得入射波的上升沿为270～300 μs，试样达到应力平衡状态的时间充足，其波形如图9所示。

图 9 整形后波形

Figure 9. Waveform after shaping

下载: 全尺寸图片幻灯片

1.3.3 入射波完全卸载

SHTB实验过程中，常将吸收杆与入射杆间留出一定空隙（2～5 mm），以避免入射杆与吸收杆之间发生二次撞击，但在PVC弹性体的SHTB实验过程中发现，入射杆与吸收杆间的空隙会导致入射波偏离基线，即未完全卸载，见图10（a）。对实验装置进行分析发现，法兰盘厚度远小于子弹长度，使得子弹与法兰盘一起运动，这是引发入射波偏离基线的主要原因。为了解决这一问题，在实验时应消除入射杆与吸收杆之间的空隙，以卸载入射波，进而获得理想实验波形，见图10（b）。

图 10 空隙对入射波卸载情况的影响

Figure 10. Effect of gap on unloading of incident wave

下载: 全尺寸图片幻灯片

2. 实验结果与分析

2.1 准静态拉伸应力-应变曲线分析

PVC弹性体在准静态加载下的拉伸载荷-位移曲线如图11所示。由图11可以看出，PVC弹性体在准静态下的拉伸力学性能基本呈线性增加趋势。从载荷-位移曲线中可以计算得出，PVC弹性体在0.1 s⁻¹应变率下的拉伸模量为27 MPa。

图 11 准静态拉伸应力-应变曲线

Figure 11. Quasi-static tensile stress-strain curve

下载: 全尺寸图片幻灯片

2.2 动态拉伸应力-应变曲线分析

对PVC弹性体动态拉伸实验数据进行处理，得到应变率为400～1850 s⁻¹时的拉伸应力-应变曲线，如图12所示。由图12可以看出，PVC弹性体在高应变率拉伸加载条件下的应力-应变曲线可分为弹性、塑性以及卸载3个阶段，呈现出明显的非线性特征，且在不同的高应变率载荷作用下具有不同的拉伸弹性模量、峰值应力、峰值应变，表现出一定的应变率敏感性。图12中，应变为0.1时的对应点与原点连线，用其斜率表征杨氏模量E_g，称之为割线模量。为进一步说明PVC弹性体的应变率效应，取图12中应变为0.1处对应的割线模量、峰值应力、峰值应变进行比较，其具体变化情况见表1。

图 12 动态拉伸应力-应变曲线

Figure 12. Dynamic tensile stress-strain curves

下载: 全尺寸图片幻灯片

表 1 动态参数变化情况

Table 1. Variations of dynamic properties

Strain rate/s⁻¹	Secant modulus/MPa	Peak stress/MPa	Peak strain/%
400	89.78	9.19	14.12
700	115.62	13.25	16.07
950	134.82	17.99	26.85
1 850	166.81	27.63	53.19

下载: 导出CSV

| 显示表格

由表1可知，当应变率范围为400～1850 s⁻¹时，割线模量、峰值应力及峰值应变均随应变率的增加而增大。

3. ZWT非线性黏弹性模型

高分子聚合物的动态力学性能可由一个非线性弹簧和两个Maxwell单元组成的朱-王-唐（ZWT）本构模型（见图13）描述，其表达式为

图 13 ZWT模型示意图

Figure 13. Schematic of ZWT model

下载: 全尺寸图片幻灯片

$\sigma = {E_0}\varepsilon + \alpha {\varepsilon ^2} + \beta {\varepsilon ^3} + {E_1}\int_0^t {\dot \varepsilon \exp \left( - \frac{{t - \tau }}{{{\theta _1}}}\right){\rm{d}}\tau } + {E_2}\int_0^t {\dot \varepsilon \exp \left( - \frac{{t - \tau }}{{{\theta _2}}}\right){\rm{d}}\tau }$

(13)

式中：前3项用于表征材料非线性弹性响应且与应变率无关，后两项分别代表低、高应变率下的黏弹性响应，E₀、 $\alpha$ 和 $\;\beta$ 为试验确定的弹性常数，E₁、 $\theta$ ₁和E₂、 $\theta$ ₂分别为试验确定的低、高应变率下的弹性常数和松弛时间。当材料受到低应变率加载时，高应变率所对应的Maxwell单元始终处于松弛状态；而当材料受到高应变率加载时，低应变率所对应的Maxwell单元则来不及松弛。相应的模型表达式可简化为

$\begin{array}{l} \sigma = {E_0}\varepsilon + \alpha {\varepsilon ^2} + \beta {\varepsilon ^3} + {E_1}\displaystyle\int_0^t {\dot \varepsilon \exp \left( - \displaystyle\frac{{t - \tau }}{{{\theta _1}}}\right){\rm{d}}\tau } \end{array}$

(14)

$\begin{array}{l} \sigma = \left( {{E_0} + {E_1}} \right)\varepsilon + \alpha {\varepsilon ^2} + \beta {\varepsilon ^3} + {E_2}\displaystyle\int_0^t {\dot \varepsilon \exp \left( - \displaystyle\frac{{t - \tau }}{{{\theta _2}}}\right){\rm{d}}\tau } \\ \end{array}$

(15)

本实验的加载率可近似看作恒应变率加载，故式（14）和式（15）可写为

$\begin{array}{l} \sigma = {E_0}\varepsilon + \alpha {\varepsilon ^2} + \beta {\varepsilon ^3} + {E_1}{\theta _1}\dot \varepsilon \left[ {{\rm{1 - }}\exp \left( { - \displaystyle\frac{\varepsilon }{{{\theta _1}\dot \varepsilon }}} \right)} \right] \end{array}$

(16)

$\begin{array}{l} \sigma = \left( {{E_0} + {E_1}} \right)\varepsilon + \alpha {\varepsilon ^2} + \beta {\varepsilon ^3} + {E_{\rm{2}}}{\theta _{\rm{2}}}\dot \varepsilon \left[ {{\rm{1 - }}\exp \left( { - \displaystyle\frac{\varepsilon }{{{\theta _{\rm{2}}}\dot \varepsilon }}} \right)} \right] \\ \end{array}$

(17)

采用式（16）和式（17）拟合PVC弹性体在高应变率下的实验数据，相应的拟合参数见表2。

表 2 ZWT模型拟合参数值

Table 2. ZWT model fitting parameter values

Strain rates/s⁻¹	E₀	E₁	$\alpha$	$\;\beta$	E₂	${\theta}$ ₂
400	−7.364	6.930	0.166	0.009	8.414	0.002 0
700	−0.350	0.630	0.076	0.004	5.074	0.002 0
950	51.249	−51.151	0.051	0.002	5.482	0.001 7
1 850	2.392	−2.550	0.027	0	8.647	0.001 0

下载: 导出CSV

| 显示表格

图14为实验结果与本构模型拟合结果的对比。可以看出，两者吻合较好，说明该本构能够较好地描述PVC弹性体在高应变率范围内的动态力学性能。

图 14 ZWT拟合曲线与实验曲线

Figure 14. Fitting curves of ZWT model and experimental curves

下载: 全尺寸图片幻灯片

4. 结　论

开展了PVC弹性体材料的动态拉伸实验，联合波形分析和高速摄像等方法对动态拉伸试样连接方式及胶黏剂进行了优选，实现了恒应变率加载与入射波完全卸载，研究了PVC弹性体在高应变率（400～1850 s⁻¹）拉伸载荷作用下的力学性能，并构建了黏弹性本构模型。通过实验结果分析发现：

（1）提出的基于SHTB装置的动态拉伸技术可用于软质高分子聚合物材料，能满足霍普金森实验的两个基本假设，测试所得结果可反映材料的动态力学性能，且装置结构简单，便于操作和推广；

（2）聚氯乙烯（PVC）弹性体在低应变率（0.1 s⁻¹）拉伸载荷作用下具有明显的非线性弹性特征，PVC弹性体在高应变率（400～1850 s⁻¹）载荷作用下具有不同的拉伸弹性模量、峰值应力、峰值应变，具有明显的应变率敏感性；

（3）PVC弹性体的动态拉伸力学性能的黏弹性特征明显，用ZWT本构模型描述特征时误差较小。

研究结果可为软质高分子聚合物材料的静、动态力学性能研究提供参考。为研究温度对PVC弹性体静、动态拉伸力学性能的影响，后期将用高/低温环境箱开展不同环境温度下的PVC弹性体材料的静、动态拉伸力学性能测试实验。

图 1 常压下H和C在金属中的扩散系数随温度的变化关系^[39]

Figure 1. Diffusion coefficients of H and C in different metals versus temperature at normal pressure^[39]

下载: 全尺寸图片幻灯片

图 2 超离子态和液态下Fe-C、Fe-H、Fe-O合金中Fe和轻元素（C、H、O）的MSD随模拟时间的变化关系

Figure 2. MSD-time curves of light elements (C, H, and O) in superionic and liquid Fe-C, Fe-H and Fe-O alloys

下载: 全尺寸图片幻灯片

图 3 260和360 GPa下C、H、O在hcp铁合金中的扩散系数随温度的变化关系

Figure 3. Diffusion coefficients of C, H, and O in hcp-Fe versus temperature at 260 and 360 GPa

下载: 全尺寸图片幻灯片

图 4 360 GPa下超离子态hcp相FeH_0.25、FeO_0.0625和FeC_0.0625的纵波波速（v_P）和横波波速（v_S）随温度的变化（文献中的计算结果：hcp-Fe的波速和预熔化效应（蓝色）^[52]、6600 K时bcc-Fe的波速（橙色）^[58]、5000 K时bcc-Fe₁₃Si₃的Voigt模型波速（青色）以及Voigt-Ruess-Hill（VRH）波速（紫色）^[59] ）

Figure 4. Compressional (v_P) and shear (v_S) wave velocities in hcp-FeH_0.25, hcp-FeO_0.0625 and hcp-FeC_0.0625 as a function of temperature at 360 GPa upon superionic transition (In comparison with previous calculation results of pure hcp-Fe with the effect of pre-melting (blue symbols)^[52], seismic velocities in bcc-Fe at 6600 K (orange symbols)^[58], and seismic velocities in bcc-Fe₁₃Si₃ at 5000 K with Voigt and Voigt-Ruess-Hill (VRH) model, respectively (cyan and purple symbols)^[59].)

下载: 全尺寸图片幻灯片

图 5 (a) 内核中地磁场示意图以及随深度变化的各向异性构造变化^[36]，(b) 蓝色区域的超离子态波速模型与观测值的对比^[66]，(c) 最内核（绿色区域）模型波速与观测值的对比^{[22, 27, 67]}

Figure 5. (a) Schematic diagram of poloidal and toroidal geomagnetic field in the inner core (IC) and depth-dependent anisotropic texture change^[36]; (b) the superionic model at the depth of blue area in comparison with geophysical observation data^[66]; (c) the superionic model of innermost inner core (green area) in comparison with geophysical observation data^{[22, 27, 67]}

下载: 全尺寸图片幻灯片

表 1 H、O、C在超离子态铁-轻元素合金中的扩散活化焓和超离子态转变温度

Table 1. Diffusion enthalpy of H, O, and C in superionic Fe-light element alloys and their superionic transition temperatures

Alloy	Pressure/GPa	D₀/(cm²·s⁻¹)	ΔH/eV	T_s/K
FeH_0.25	260	0.148	1.58	2000
FeH_0.25	360	0.317	2.14	2230
FeO_0.0625	260	0.079	1.62	2350
FeO_0.0625	360	0.067	1.72	2630
FeC_0.0625	260	0.087	1.83	2600
FeC_0.0625	360	0.074	1.60	2360

下载: 导出CSV

参考文献(67)

[1]	FARADAY M. Ⅶ. experimental researches in electricity: twelfth series [J]. Philosophical Transactions of the Royal Society of London, 1838, 128: 83–123. doi: 10.1098/rstl.1838.0008
[2]	TUBANDT C, LORENZ E. Molekularzustand und elektrisches leitvermögen kristallisierter salze [J]. Zeitschrift für Physikalische Chemie, 1914, 87U(1): 513–542. doi: 10.1515/zpch-1914-8737
[3]	DEMONTIS P, LESAR R, KLEIN M L. New high-pressure phases of ice [J]. Physical Review Letters, 1988, 60(22): 2284–2287. doi: 10.1103/PhysRevLett.60.2284
[4]	MILLOT M, COPPARI F, RYGG J R, et al. Nanosecond X-ray diffraction of shock-compressed superionic water ice [J]. Nature, 2019, 569(7755): 251–255. doi: 10.1038/s41586-019-1114-6
[5]	CAVAZZONI C, CHIAROTTI G L, SCANDOLO S, et al. Superionic and metallic states of water and ammonia at giant planet conditions [J]. Science, 1999, 283(5398): 44–46. doi: 10.1126/science.283.5398.44
[6]	REDMER R, MATTSSON T R, NETTELMANN N, et al. The phase diagram of water and the magnetic fields of Uranus and Neptune [J]. Icarus, 2011, 211(1): 798–803. doi: 10.1016/j.icarus.2010.08.008
[7]	LIU C, GAO H, WANG Y, et al. Multiple superionic states in helium-water compounds [J]. Nature Physics, 2019, 15(10): 1065–1070. doi: 10.1038/s41567-019-0568-7
[8]	KIMURA T, MURAKAMI M. Fluid-like elastic response of superionic NH₃ in Uranus and Neptune [J]. Proceedings of the National Academy of Sciences of the United States of America, 2021, 118(14): e2021810118.
[9]	BINNS J, HERMANN A, PEÑA-ALVAREZ M, et al. Superionicity, disorder, and bandgap closure in dense hydrogen chloride [J]. Science Advances, 2021, 7(36): eabi9507. doi: 10.1126/sciadv.abi9507
[10]	LI H F, OGANOV A R, CUI H X, et al. Ultrahigh-pressure magnesium hydrosilicates as reservoirs of water in early Earth [J]. Physical Review Letters, 2022, 128(3): 035703. doi: 10.1103/PhysRevLett.128.035703
[11]	LI J W, LIN Y H, MEIER T. Silica-water superstructure and one-dimensional superionic conduit in Earth’s mantle [J]. Science Advances, 2023, 9(35): eadh3784. doi: 10.1126/sciadv.adh3784
[12]	HOU M Q, HE Y, JANG B G, et al. Superionic iron oxide-hydroxide in Earth’s deep mantle [J]. Nature Geoscience, 2021, 14(3): 174–178. doi: 10.1038/s41561-021-00696-2
[13]	HIROSE K, WOOD B, VOČADLO L, et al. Light elements in the Earth’s core [J]. Nature Reviews Earth & Environment, 2021, 2(9): 645–658. doi: 10.1038/s43017-021-00203-6
[14]	刘锦, 吕超甲, 赵超帅. 矿物高压稳定性与深部挥发分循环过程 [J]. 矿物岩石地球化学通报, 2022, 41(2): 245–259. doi: 10.19658/j.issn.1007-2802.2022.41.011 LIU J, LYU C J, ZHAO C S. High-pressure stability of minerals and volatiles cycling in the deep Earth [J]. Bulletin of Mineralogy, Petrology and Geochemistry, 2022, 41(2): 245–259. doi: 10.19658/j.issn.1007-2802.2022.41.011
[15]	BIRCH F. Elasticity and constitution of the Earth’s interior [J]. Journal of Geophysical Research: Solid Earth, 1952, 57(2): 227–286. doi: 10.1029/JZ057i002p00227
[16]	TATENO S, HIROSE K, OHISHI Y, et al. The structure of iron in Earth’s inner core [J]. Science, 2010, 330(6002): 359–361. doi: 10.1126/science.1194662
[17]	TURNEAURE S J, SHARMA S M, GUPTA Y M. Crystal structure and melting of Fe shock compressed to 273 GPa: in situ X-ray diffraction [J]. Physical Review Letters, 2020, 125(21): 215702. doi: 10.1103/PhysRevLett.125.215702
[18]	KRAUS R G, HEMLEY R J, ALI S J et al. Measuring the melting curve of iron at super-Earth core conditions [J]. Science, 2022, 375(6577): 202–205. doi: 10.1126/science.abm1472
[19]	VOČADLO L, ALFÈ D, GILLAN M J, et al. Possible thermal and chemical stabilization of body-centred-cubic iron in the Earth’s core [J]. Nature, 2003, 424(6948): 536–539. doi: 10.1038/nature01829
[20]	BELONOSHKO A B, AHUJA R, JOHANSSON B. Stability of the body-centred-cubic phase of iron in the Earth’s inner core [J]. Nature, 2003, 424(6952): 1032–1034. doi: 10.1038/nature01954
[21]	BELONOSHKO A B, LUKINOV T, FU J, et al. Stabilization of body-centred cubic iron under inner-core conditions [J]. Nature Geoscience, 2017, 10(4): 312–316. doi: 10.1038/ngeo2892
[22]	DZIEWONSKI A M, ANDERSON D L. Preliminary reference Earth model [J]. Physics of the Earth and Planetary Interiors, 1981, 25(4): 297–356. doi: 10.1016/0031-9201(81)90046-7
[23]	SONG X D. Anisotropy of the Earth’s inner core [J]. Reviews of Geophysics, 1997, 35(3): 297–313. doi: 10.1029/97RG01285
[24]	DEUSS A. Heterogeneity and anisotropy of Earth’s inner core [J]. Annual Review of Earth and Planetary Sciences, 2014, 42: 103–126. doi: 10.1146/annurev-earth-060313-054658
[25]	MORELLI A, DZIEWONSKI A M, WOODHOUSE J H. Anisotropy of the inner core inferred from PKIKP travel times [J]. Geophysical Research Letters, 1986, 13(13): 1545–1548. doi: 10.1029/GL013i013p01545
[26]	TANAKA S, HAMAGUCHI H. Degree one heterogeneity and hemispherical variation of anisotropy in the inner core from PKP(BC)-PKP(DF) times [J]. Journal of Geophysical Research: Solid Earth, 1997, 102(B2): 2925–2938. doi: 10.1029/96JB03187
[27]	ISHII M, DZIEWOŃSKI A M. The innermost inner core of the Earth: evidence for a change in anisotropic behavior at the radius of about 300 km [J]. Proceedings of the National Academy of Sciences of the United States of America, 2002, 99(22): 14026–14030. doi: 10.1073/pnas.172508499
[28]	SONG X D, RICHARDS P G. Seismological evidence for differential rotation of the Earth’s inner core [J]. Nature, 1996, 382(6588): 221–224. doi: 10.1038/382221a0
[29]	BADDING J V, HEMLEY R J, MAO H K, et al. High-pressure chemistry of hydrogen in metals: in situ study of iron hydride [J]. Science, 1991, 253(5018): 421–424. doi: 10.1126/science.253.5018.421
[30]	PÉPIN C M, DEWAELE A, GENESTE G, et al. New iron hydrides under high pressure [J]. Physical Review Letters, 2014, 113(26): 265504. doi: 10.1103/PhysRevLett.113.265504
[31]	LI Y G, VOČADLO L, SUN T, et al. The Earth’s core as a reservoir of water [J]. Nature Geoscience, 2020, 13(6): 453–458. doi: 10.1038/s41561-020-0578-1
[32]	YUAN L, STEINLE-NEUMANN G. Strong sequestration of hydrogen into the Earth’s core during planetary differentiation [J]. Geophysical Research Letters, 2020, 47(15): e2020GL088303. doi: 10.1029/2020GL088303
[33]	TAGAWA S, SAKAMOTO N, HIROSE K, et al. Experimental evidence for hydrogen incorporation into Earth’s core [J]. Nature Communications, 2021, 12(1): 2588. doi: 10.1038/s41467-021-22035-0
[34]	HE Y, KIM D Y, STRUZHKIN V V, et al. The stability of FeH_x and hydrogen transport at Earth’s core mantle boundary [J]. Science Bulletin, 2023, 68(14): 1567–1573. doi: 10.1016/j.scib.2023.06.012
[35]	HE Y, SUN S C, KIM D Y, et al. Superionic iron alloys and their seismic velocities in Earth’s inner core [J]. Nature, 2022, 602(7896): 258–262. doi: 10.1038/s41586-021-04361-x
[36]	SUN S C, HE Y, YANG J Y, et al. Superionic effect and anisotropic texture in Earth’s inner core driven by geomagnetic field [J]. Nature Communications, 2023, 14(1): 1656. doi: 10.1038/s41467-023-37376-1
[37]	WANG W Z, LI Y G, BRODHOLT J P, et al. Strong shear softening induced by superionic hydrogen in Earth’s inner core [J]. Earth and Planetary Science Letters, 2021, 568: 117014. doi: 10.1016/j.jpgl.2021.117014
[38]	YANG H, DOU P X, XIAO T T, et al. The geophysical properties of FeH_x phases under inner core conditions [J]. Geophysical Research Letters, 2023, 50(22): e2023GL104493. doi: 10.1029/2023GL104493
[39]	FUKAI Y, SUGIMOTO H. Diffusion of hydrogen in metals [J]. Advances in Physics, 1985, 34(2): 263–326. doi: 10.1080/00018738500101751
[40]	SHIBAZAKI Y, OHTANI E, FUKUI H, et al. Sound velocity measurements in dhcp-FeH up to 70 GPa with inelastic X-ray scattering: implications for the composition of the Earth’s core [J]. Earth and Planetary Science Letters, 2012, 313/314: 79–85. doi: 10.1016/j.jpgl.2011.11.002
[41]	MASHINO I, MIOZZI F, HIROSE K, et al. Melting experiments on the Fe-C binary system up to 255 GPa: constraints on the carbon content in the Earth’s core [J]. Earth and Planetary Science Letters, 2019, 515: 135–144. doi: 10.1016/j.jpgl.2019.03.020
[42]	OZAWA H, HIROSE K, TATENO S, et al. Phase transition boundary between B1 and B8 structures of FeO up to 210 GPa [J]. Physics of the Earth and Planetary Interiors, 2010, 179(3/4): 157–163. doi: 10.1016/j.pepi.2009.11.005
[43]	OISHI Y, KAMEI Y, AKIYAMA M, et al. Self-diffusion coefficient of lithium in lithium oxide [J]. Journal of Nuclear Materials, 1979, 87(2/3): 341–344. doi: 10.1016/0022-3115(79)90570-1
[44]	KARKI B B, STIXRUDE L, CLARK S J, et al. Structure and elasticity of MgO at high pressure [J]. American Mineralogist, 1997, 82(1/2): 51–60. doi: 10.2138/am-1997-1-207
[45]	KARKI B B, STIXRUDE L, WENTZCOVITCH R M. High-pressure elastic properties of major materials of Earth’s mantle from first principles [J]. Reviews of Geophysics, 2001, 39(4): 507–534. doi: 10.1029/2000RG000088
[46]	VOIGT W. Lehrbuch der kristallphysik [M]. Leipzig: Teubner Verlag, 1928.
[47]	REUSS A. Berechnung der fließgrenze von mischkristallen auf grund der plastizitätsbedingung fűr einkristalle [J]. Journal of Applied Mathematics and Mechanics, 1929, 9(1): 49–58. doi: 10.1002/zamm.19290090104
[48]	HILL R. The elastic behaviour of a crystalline aggregate [J]. Proceedings of the Physical Society: Section A, 1952, 65(5): 349–354.
[49]	ANDERSON D L. Theory of the Earth [M]. Boston: Blackwell Scientific Publications, 1989.
[50]	HULL S, FARLEY T W D, HAYES W, et al. The elastic properties of lithium oxide and their variation with temperature [J]. Journal of Nuclear Materials, 1988, 160(2/3): 125–134. doi: 10.1016/0022-3115(88)90039-6
[51]	HE Y, SUN S C, LI H P. Ab initio molecular dynamics investigation of the elastic properties of superionic Li₂O under high temperature and pressure [J]. Physical Review B, 2021, 103(17): 174105. doi: 10.1103/PhysRevB.103.174105
[52]	MARTORELL B, VOČADLO L, BRODHOLT J, et al. Strong premelting effect in the elastic properties of hcp-Fe under inner-core conditions [J]. Science, 2013, 342(6157): 466–468. doi: 10.1126/science.1243651
[53]	甘波, 李俊, 蒋刚, 等. Fe高压熔化线的实验研究进展 [J]. 高压物理学报, 2021, 35(6): 060101. doi: 10.11858/gywlxb.20210859 GAN B, LI J, JIANG G, et al. A review of the experimental determination of the melting curve of iron at ultrahigh pressures [J]. Chinese Journal of High Pressure Physics, 2021, 35(6): 060101. doi: 10.11858/gywlxb.20210859
[54]	ANZELLINI S, DEWAELE A, MEZOUAR M, et al. Melting of iron at Earth’s inner core boundary based on fast X-ray diffraction [J]. Science, 2013, 340(6131): 464–466. doi: 10.1126/science.1233514
[55]	LI J, WU Q, LI J B, et al. Shock melting curve of iron: a consensus on the temperature at the Earth’s inner core boundary [J]. Geophysical Research Letters, 2020, 47(15): e2020GL087758. doi: 10.1029/2020GL087758
[56]	HOU H Q, LIU J, ZHANG Y J, et al. Melting of iron explored by electrical resistance jump up to 135 GPa [J]. Geophysical Research Letters, 2021, 48(20): e2021GL095739. doi: 10.1029/2021GL095739
[57]	ZHANG Y J, WANG Y, HUANG Y Q, et al. Collective motion in hcp-Fe at Earth’s inner core conditions [J]. Proceedings of the National Academy of Sciences of the United States of America, 2023, 120(41): e2309952120.
[58]	BELONOSHKO A B, SIMAK S I, OLOVSSON W, et al. Elastic properties of body-centered cubic iron in Earth’s inner core [J]. Physical Review B, 2022, 105(18): L180102. doi: 10.1103/PhysRevB.105.L180102
[59]	WU Z Q, WANG W Z. Shear softening of Earth’s inner core as indicated by its high Poisson ratio and elastic anisotropy [J]. Fundamental Research. DOI: 10.1016/j.fmre.2022.08.010.
[60]	BELONOSHKO A B, SKORODUMOVA N V, DAVIS S, et al. Origin of the low rigidity of the Earth’s inner core [J]. Science, 2007, 316(5831): 1603–1605. doi: 10.1126/science.1141374
[61]	KARATO S I. Inner core anisotropy due to the magnetic field: induced preferred orientation of iron [J]. Science, 1993, 262(5140): 1708–1711. doi: 10.1126/science.262.5140.1708
[62]	YOSHIDA S, SUMITA I, KUMAZAWA M. Growth model of the inner core coupled with the outer core dynamics and the resulting elastic anisotropy [J]. Journal of Geophysical Research: Solid Earth, 1996, 101(B12): 28085–28103. doi: 10.1029/96JB02700
[63]	FROST D A, LASBLEIS M, CHANDLER B, et al. Dynamic history of the inner core constrained by seismic anisotropy [J]. Nature Geoscience, 2021, 14(7): 531–535. doi: 10.1038/s41561-021-00761-w
[64]	SINGH S C, TAYLOR M A J, MONTAGNER J P. On the presence of liquid in Earth’s inner core [J]. Science, 2000, 287(5462): 2471–2474. doi: 10.1126/science.287.5462.2471
[65]	何宇, 孙士川, 徐云帆, 等. 地球内核各向异性结构成因: 矿物学模型和动力学机制 [J/OL]. 矿物岩石地球化学通报 (2023-09-28)[2024-01-09]. https://kns.cnki.net/kcms2/article/abstract?v=z5VdU6XQV3X2PXbTT1OgcuOSUGALw3UfEUzDIb8cGLjV1OFDYBL1x9lftDrN9bJg3zFaPXciDx4-jtQR4v9o1dgio2ra8UjlPAeskyFgNML5Wt6L8hM4_CvCzP0zQsJS-1kFXp9zKdE=&uniplatform=NZKPT&language=CHS. DOI: 10.19658/j.issn.1007-2802.2023.42.100. HE Y, SUN S C, XU Y F, et al. The origin of anisotropic structure of the Earth’s inner core from mineralogical model to dynamic mechanism [J/OL]. Bulletin of Mineralogy, Petrology and Geochemistry (2023-09-28)[2024-01-09]. https://kns.cnki.net/kcms2/article/abstract?v=z5VdU6XQV3X2PXbTT1OgcuOSUGALw3UfEUzDIb8cGLjV1OFDYBL1x9lftDrN9bJg3zFaPXciDx4-jtQR4v9o1dgio2ra8UjlPAeskyFgNML5Wt6L8hM4_CvCzP0zQsJS-1kFXp9zKdE=&uniplatform=NZKPT&language=CHS. DOI: 10.19658/j.issn.1007-2802.2023.42.100.
[66]	BRETT H, DEUSS A. Inner core anisotropy measured using new ultra-polar PKIKP paths [J]. Geophysical Journal International, 2020, 223(2): 1230–1246. doi: 10.1093/gji/ggaa348
[67]	SUN X L, SONG X D. The inner inner core of the Earth: texturing of iron crystals from three-dimensional seismic anisotropy [J]. Earth and Planetary Science Letters, 2008, 269(1/2): 56–65. doi: 10.1016/j.jpgl.2008.01.049

施引文献

资源附件(0)

访问统计

点击查看大图

图(5) / 表(1)

计量

文章访问数: 298
HTML全文浏览量: 174
PDF下载量: 51

1. 实验材料及测试方法
1.1 准静态拉伸实验
1.2 动态拉伸实验
1.3 动态拉伸实验关键问题及解决方法
2. 实验结果与分析
2.1 准静态拉伸应力-应变曲线分析
2.2 动态拉伸应力-应变曲线分析
3. ZWT非线性黏弹性模型
4. 结　论

地球内核超离子态铁合金及其效应

doi: 10.11858/gywlxb.20240707

作者简介: 何 宇（1985－），男，博士，研究员，主要从事高温高压下地球内部物质性质的计算模拟研究.E-mail：heyu@mail.gyig.ac.cn

计量

出版历程

Superionic Iron Alloys in Earth’s Inner Core and Their Effects

1. 实验材料及测试方法

1.1 准静态拉伸实验

1.2 动态拉伸实验

1.3 动态拉伸实验关键问题及解决方法

1.3.1 试样连接方式的确定及胶黏剂优选

1.3.2 恒应变率加载

1.3.3 入射波完全卸载

2. 实验结果与分析

2.1 准静态拉伸应力-应变曲线分析

2.2 动态拉伸应力-应变曲线分析

3. ZWT非线性黏弹性模型

4. 结 论

计量

出版历程

目录

1. 实验材料及测试方法

1.1 准静态拉伸实验

1.2 动态拉伸实验

1.3 动态拉伸实验关键问题及解决方法

2. 实验结果与分析

2.1 准静态拉伸应力-应变曲线分析

2.2 动态拉伸应力-应变曲线分析

3. ZWT非线性黏弹性模型

4. 结 论

作者简介:
何　宇（1985－），男，博士，研究员，主要从事高温高压下地球内部物质性质的计算模拟研究.E-mail：heyu@mail.gyig.ac.cn

4. 结　论

4. 结　论