Loading [MathJax]/jax/output/SVG/jax.js

固体氢在极端压强下的超导性能

杜昱 孙莹 王彦超 钟鑫

尹华伟, 蒋轲, 张料, 黄亮, 王陈凌. 钢纤维混凝土板在冲击与爆炸荷载下的K&C模型[J]. 高压物理学报, 2020, 34(3): 034205. doi: 10.11858/gywlxb.20190853
引用本文: 杜昱, 孙莹, 王彦超, 钟鑫. 固体氢在极端压强下的超导性能[J]. 高压物理学报, 2024, 38(2): 020105. doi: 10.11858/gywlxb.20230722
YIN Huawei, JIANG Ke, ZHANG Liao, HUANG Liang, WANG Chenling. K&C Model of Steel Fiber Reinforced Concrete Plate under Impact and Blast Load[J]. Chinese Journal of High Pressure Physics, 2020, 34(3): 034205. doi: 10.11858/gywlxb.20190853
Citation: DU Yu, SUN Ying, WANG Yanchao, ZHONG Xin. Superconductivity of Solid Hydrogen under Extreme Pressure[J]. Chinese Journal of High Pressure Physics, 2024, 38(2): 020105. doi: 10.11858/gywlxb.20230722

固体氢在极端压强下的超导性能

doi: 10.11858/gywlxb.20230722
基金项目: 国家自然科学基金(52288102,52090024);科技部重点研发计划(2022YFA1402304);中国科学院战略性先导科技专项(XDB33000000)
详细信息
    作者简介:

    杜 昱(1999-),女,博士研究生,主要从事极端高压下的计算凝聚态物理研究.E-mail:duyu21@mails.jlu.edu.cn

    通讯作者:

    钟 鑫(1988-), 女,博士,副教授,主要从事高压下新型高温超导材料设计与超导电性研究.E-mail:zx777@jlu.edu.cn

  • 中图分类号: O521.2

Superconductivity of Solid Hydrogen under Extreme Pressure

  • 摘要: 氢元素在常压下具有最简单的晶体结构和物理性质。随着压强增加,氢单质发生相变,由绝缘体转变为金属,被称为金属氢。数值模拟表明,金属氢具有高温超导电性,因此,金属氢研究也被称为高压物理领域的“圣杯”课题。利用基于密度泛函理论的第一性原理计算方法,对固体氢在极端高压(0.5~5.0 TPa)下的结构和超导电性开展了系统研究。研究结果表明:固体氢的高压相变序列为I41/amd→oC12→cI16;对于同一种结构,随着压强增加,电声耦合系数减小,费米面处电子态密度减小,特征振动频率增加,超导转变温度发生小幅变化;在2.0 TPa压强下,固体氢的超导转变温度高达418 K(库伦赝势经验值μ*=0.10)。研究工作将为金属氢及其超导电性的后续理论和实验研究提供参考。

     

  • 普通混凝土具有显著的脆性特征,钢纤维的引入极大地提高了混凝土的延性,并可减少结构部件的剥落飞射,避免了大量的人员伤亡,因此被广泛应用于各类防护结构。近几十年来,随着计算能力的高速发展,LS-DYNA显式有限元动力分析软件被广泛应用于工程结构的动态行为研究。

    K&C模型常用于分析冲击荷载与爆炸荷载作用下混凝土结构的动力响应,研究人员只需输入混凝土的单轴抗压强度,其他参数可自动生成,然而该模型存在以下不足:(1)自动生成的强度面参数适用于普通混凝土,但不适用于三轴强度显著变化的钢纤维混凝土(Steel fiber reinforced concrete,SFRC);(2)K&C模型的损伤机制能够表征普通混凝土的力学行为,但对于韧性、延性、能量吸收能力显著增强的SFRC,不能准确表征其应变软化行为;(3)SFRC的应变率效应在很大程度上依赖于钢纤维的体积分数、抗拉强度、几何形状和混凝土的基体强度,因此普通混凝土的动力增强因子模型不适用于SFRC。

    本研究将针对以上不足,对K&C模型进行改进,并运用LS-DYNA显式有限元动力分析软件模拟高速弹丸侵彻SFRC靶板[1]和爆炸荷载下SFRC板的动力响应[2],利用模拟结果论证改进后的K&C模型,以期能够准确描述SFRC的动力响应。

    K&C模型采用压子午线定义了3个独立的失效强度面

    Δσm=a0+pa1+a2p,Δσy=a0y+pa1y+a2yp,Δσr=pa1f+a2fp (1)

    式中:Δσm为最大强度面;Δσy为初始屈服强度面;Δσr为残余强度面;aiaiyaifi = 0, 1, 2)为混凝土三轴压缩试验得到的一组参数;p为静水压力,p=(σ1+σ2+σ3)/3,其中σ1σ2σ3分别为第1、第2、第3主应力;失效强度面偏应力Δσ=[(σ1σ2)2+(σ2σ3)2+(σ3σ1)2]/2。原K&C模型失效强度面参数为

    a0=0.296fc,a1=0.446,a2=0.081/fc (2)
    a0y=0.223fc,a1y=0.625,a2y=0.258/fc (3)
    a1f=0.442,a2f=0.118/fc (4)

    式中:fc为单轴抗压强度。

    国内外学者研究发现,SFRC的三轴压缩强度与素混凝土有显著差异,且依赖于钢纤维体积分数[3-11]。为了明确钢纤维体积分数对SFRC三轴压缩强度的影响,本研究在原最大强度面参数模型的基础上,根据大量实验数据,采用最小二乘法建立了新的最大强度面参数模型

    a0=(0.232+0.133Vf)fc,a1=0.378+0.103Vf,a2=(0.1490.047Vf)/fc (5)

    式中:Vf为钢纤维体积分数,0 < Vf < 3.0%。

    图1给出了不同钢纤维体积分数SFRC的三轴压缩实验数据及根据式(2)、式(5)模型计算得到的预测值。由图1可知,式(5)与不同钢纤维体积分数的SFRC三轴压缩实验数据[3-11]的吻合程度优于式(2)。

    图  1  钢纤维体积分数不同时的三轴压缩数据以及式(2)、式(5)的预测值
    Figure  1.  Triaxial compression data of different fiber contents as well as predictions from Eq.(2) and Eq.(5)

    为了确定初始屈服面,必须首先确定a0ya1ya2y。通常认为初始屈服面近似为三轴受压应力路径上Δσ=0.45Δσm处点的轨迹。最大破坏面上的点(p,Δσm)与屈服面上对应的点(p,Δσy)的关系可表示为

    Δσy=0.45Δσm,p=p0.55Δσm/3 (6)

    根据式(6)的假定及已提出的最大强度面模型,初始屈服面参数模型可表示为

    a0y=(0.194+0.068Vf)fc,a1y=0.524+0.202Vf,a2y=(0.4520.142Vf)/fc (7)

    残余强度面参数的确定需要用到三轴压缩实验的残余强度数据,但相关的实验数据难以得到。根据Marlvar等[12]的建议,残余强度面参数可表示为

    a1f=0.374+0.102Vf,a2f=(0.2180.068Vf)/fc (8)

    改进的失效强度面模型应按式(5)、式(7)、式(8)进行修正,当fc = 175.3 MPa,Vf = 1.0%时,SFRC的失效强度面模型如图2所示,该模型将应用于2.1节的讨论。

    图  2  改进的失效强度面
    Figure  2.  Modified failure strength surface

    在K&C模型中,当偏应力超过初始屈服面,但未达到最大强度面时,当前多轴应力状态为两个强度面的线性插值;当偏应力超过最大强度面时,当前多轴应力状态取为最大强度面与残余强度面的线性插值,其关系式表示为

    Δσ={η(ΔσmΔσy)+Δσyλ<λmη(ΔσmΔσr)+Δσrλλm (9)

    式中:λ为有效塑性应变,η为损伤比例因子,损伤演化函数η(λ)由输入的13组(η, λ)数据定义。演化全过程分为两个阶段:(1) λ < λm时为硬化段,λ从零开始逐渐积累到λmη逐渐增大到1,材料的应变硬化发展到峰值,此时偏应力点位于最大强度面;(2)λλm时,材料进入应变软化阶段,λm持续积累,η逐渐减小为零,此时偏应力点位于残余强度面。

    ηλ的关系描述了应力-应变全过程曲线,必须谨慎定义。原K&C模型的损伤演化函数η(λ)描述的损伤演化速度过快,不能正确表征SFRC的延性、韧性及断裂应变能。为了精准描述由塑性剪切应变和塑性体积应变引起的损伤演变过程,需要对损伤演化函数进行修正。

    本研究使用分段函数[13-14]建立新的损伤演化模型,其关系式表示为

    η(λ)={α(λλm)+(32α)(λλm)2+(α2)(λλm)3λ<λmλ/λmαc(λ/λm1)αd+λ/λmλλm (10)

    式中:α为控制应变硬化段的常量,αcαd为控制应变软化段的常量,λm为有效塑性应变峰值。在LS-DYNA软件中建立单个单元模型,单元尺寸为2 mm,约束底面节点自由度,对顶面节点施加速度矢量。模拟SFRC的单轴压缩试验,通过反复试验法确定ααcαdλm,直到模拟的应力-应变曲线与典型的SFRC应力-应变关系相符。基于大量试验研究,建议α = 3,αc = 0.381,αd = 1.90,λm = 9.5 × 10–5。改进的η(λ)将用于接下来的讨论,与原损伤演化函数的比较如图3所示。

    图  3  损伤演化函数η(λ)的比较
    Figure  3.  Comparisons of damage evolution function η(λ)

    虽然本节对η(λ)进行了修正,但是K&C模型中材料的软化行为由η(λ)与损伤参数b1b2协同控制,因此要准确描述SFRC的软化段特性,还需对损伤参数进行校准,1.3节将介绍校准损伤参数b1b2的简便方法,并进一步论证本研究所建立损伤演化函数η(λ)的有效性。

    K&C模型中,有效塑性应变积累函数为

    λ={ˉεp0dˉεpγf(1+p/γfft)b1p0ˉεp0dˉεpγf(1+p/γfft)b2p<0 (11)

    式中:ˉεp为有效塑性应变增量;γf为应变率增强因子;ft为单轴抗拉强度;b1b2为损伤参数,分别控制单轴压缩与拉伸软化段的力学行为。

    Malvar等[12]建议损伤参数b1的取值应与单轴压缩试验(UUC)下的应力-应变曲线相匹配,建立单个单元模型,模拟fc = 175.3 MPa、Vf = 1.0%时SFRC的单轴压缩试验[1]。原K&C模型自动生成的损伤演化函数、损伤参数(b1 = 0)与修正后的损伤演化函数、损伤参数(b1 = 0.75)对应力-应变曲线的描述如图4所示。由图4可知,原K&C模型不能准确描述SFRC应变软化段的力学行为,而修正后的损伤模型较好地预测了SFRC的峰后特性,且与试验数据的吻合程度较高,进一步论证了修正后损伤演化函数η(λ)的有效性。

    图  4  单轴压缩试验的应力-应变曲线与K&C模型的预测结果比较
    Figure  4.  Comparisons of experiment at stress-strain curve of UUC test and prediction of K&C models

    当应力点从最大强度面向残余强度破坏面移动时,材料进入软化阶段,参数b2控制了受拉应力-应变曲线的软化行为。Bazant等[15]认为,只有当SFRC的软化行为受到局部或特征长度控制时,网格细化后的应力-应变曲线才会客观,否则它们将依赖于网格尺寸。他们提出消除网格依赖性的方法是强迫单轴拉伸应力-应变曲线下方的面积为Gf/hw,其中,Gf为断裂应变能,hw为裂缝正面宽度或裂缝局部宽度,通常hw为最大骨料直径的1~6倍。Malvar等[12]归纳总结出b2取值的表达式为

    b2=(0.09h2w0.98hw+3.06)(10.004f2c+0.097fc0.484) (12)

    众所周知,当钢纤维掺入混凝土中时,SFRC的断裂能呈阶跃式增长,且随着钢纤维体积分数的增长而增长。因此式(12)对于普通混凝土是适宜的,却不能正确反映SFRC对能量的吸收能力。根据董振英等[16]的实验得到的SFRC的断裂应变能481 N/m,建立单个单元模型,以模拟单轴拉伸试验(UUT),采用Bazant建议的方法,强迫单轴拉伸应力-应变曲线下方的面积为Gf/hw,校准SFRC的拉伸损伤参数(b2 = 0.2),由式(12)自动生成的拉伸损伤参数b2 = 2.1。原模型与校准损伤参数及修正损伤演化函数后的K&C模型对单轴拉伸应力-应变曲线的描述如图5所示。由图5可知,修正后的K&C模型的应力-应变曲线软化段表现出更高的韧性、延性和断裂应变能,且与实验数据的吻合程度较高。

    图  5  单轴拉伸试验的应力-应变曲线与K&C模型的预测结果比较
    Figure  5.  Comparisons of experiment at stress-strain curve of UUT test and prediction of K&C models

    SFRC的骨料强度大于砂浆与交界面强度,在高应变率下,断裂机理发生改变,裂纹扩展极为迅速,裂纹的扩展可能直接通过骨料区域,而不是骨料与砂浆的交界面,混凝土的单轴受压强度增长可达100%以上,而单轴受拉强度增长可达600%以上。目前,已有几种应变率模型用于预测混凝土在动态压缩荷载下的强度增强效应。被广泛接受的模型由CEB-FIP[17]模式规范提出,定义动态强度和静态强度的比值与应变率的函数关系,即受压动力增强因子(CDIF)表示为

    CDIF=fc,dfc,s={(˙ε˙εs)1.026δ˙ε<30s1γs(˙ε˙εs)1/3˙ε30s1 (13)

    式中:fc,d为动态压缩强度;fc,s为静态压缩强度;˙ε为当前应变率;˙εs为基准静态应变率,˙εs = 3×105s1;lg γs = 6.156δ – 2,δ = fc,0/(5fc,0 + 9fc,s),fc,0 = 10 MPa。

    Wang等[18]开展了SFRC的霍普金森压杆试验研究,在CEB-FIP公式的基础上考虑了纤维效应,修正了过渡应变,提出了CDIF模型

    CDIF={(˙ε˙εs)1.026δ˙ε<(30+23i)s1η(˙ε˙εs)(1+0.05i)/3˙ε(30+23i)s1 (14)

    式中:i为考虑纤维效应的参数,i等于0或1,分别代表素混凝土及纤维增强混凝土;η = γs(1 – 0.339 2i),其他变量同CEB-FIP模型。

    CEB-FIP模型能够较好地预测普通混凝土的DIF值,但对SFRC的DIF值估计过高,如图6所示。众多学者的研究数据[19-23]表明,过渡应变率之后,SFRC的应变率敏感性低于普通混凝土,原因在于:(1)根据ACI446[24]的研究结果,高应变率下混凝土强度的增长速率随着裂缝扩展速度的增大而增大,而SFRC的高密实性和钢纤维的阻裂作用使得裂缝扩展速率与素混凝土相比较慢;(2)钢纤维的应变率敏感性低于混凝土基体。因此需提出适用于SFRC的CDIF模型。

    图  6  不同模型得到的CDIF值与实验数据的比较
    Figure  6.  Comparisons of CDIF values obtained from different models and experiments

    由于侧向惯性效应,在过渡应变率之后,DIF曲线的斜率急剧增大,表明混凝土的物理性能发生了变化。钢纤维的约束效应有效地减少了混凝土的横向膨胀,使得惯性效应更加显著,需要更高的过渡应变率触发SFRC的物理性能变化,因此有必要提出适用于SFRC的过渡应变率模型。

    Wang等[18]提出的模型虽与众多学者的试验数据边界大致相符,但没有给出具体的钢纤维体积分数所对应的曲线,对过渡应变率模型的描述也不够精确。本研究在式(14)的基础上,引入了钢纤维体积分数Vf,并修正过渡应变率,通过最小二乘法回归分析大量实验数据,给出CDIF值,即

    CDIF={(˙ε˙εs)1.026δ˙ε<εts1γVfγs(˙ε˙εs)(1+0.018Vf)/3˙εεts1 (15)

    式中:εt为过渡应变率,εt= 41.087 + 7.677 Vf – 0.746 V2fγVf= 0.607exp( –Vf/4.972) + 0.294,γVf为钢纤维体积分数因子;其他变量同CEB-FIP模型。如图7所示,式(15)得到的CDIF值与实验数据[19-23]吻合较好。

    图  7  新模型得到的CDIF值与实验数据的比较
    Figure  7.  Comparisons of CDIF values obtained from new models and experiments

    混凝土的动态抗拉强度表现出比抗压强度更高的应变率敏感性。目前研究人员在各自实验研究[17, 25-28]的基础上,提出了数种SFRC受拉动力增强因子模型(TDIF),如图8所示。由图8可以看出,TDIF模型存在较大差异,可能是因为混凝土的基体强度、纤维含量、形状、长径比、抗拉强度对TDIF值的影响较大。由于现有实验数据的分歧较大,尚不具备提出新的TDIF模型的条件,因此本研究中采用被广泛接受的CEB-FIP公式作为受拉增强因子模型

    图  8  不同模型得到的TDIF比较
    Figure  8.  Comparisons of TDIF values obtained from different models
    TDIF=ft,dft,s={(˙ε˙εs)1.016δ˙ε<30s1β(˙ε˙εs)1/3˙ε30s1 (16)

    式中:ft,d为动态抗拉强度;ft,s为静态抗拉强度;˙εs = 3 × 10–6 s–1δ = fc,0/(10fc,0 + 6fc),fc,0 = 10 MPa;lg β = 7.11δ – 2.33。

    为了论证第1节中改进的失效强度面参数模型、损伤模型、受压动力增强因子模型的有效性与可靠性,对高速弹丸侵彻SFRC板[1]进行仿真分析。数值模拟所用SFRC板的钢纤维体积分数为1%,尺寸设定为310 mm × 310 mm × 50 mm,支撑于钢框架之上,自由跨度为240 mm。模拟所用的平头弹丸直径为25 mm,长度为75 mm,质量为297 g;弹丸材料为SKH-51工具钢,洛氏硬度HRC设定为60,侵彻过程中塑性变形较小。弹丸分别以58.2、76.0和104.0 m/s的速度对SFRC板进行模拟冲击实验。

    SFRC板的1/4计算模型如图9所示,对称面上施加对称约束。采用SOLID164单元对弹体和SFRC板进行离散,单元尺寸为2 mm,采用单点高斯积分拉格朗日算法。SFRC板采用改进的K&C材料模型,耦合*EOS_TABULATED_COMPACTION状态方程。最大强度面、初始屈服强度面、残余强度面参数值分别按式(5)、式(7)、式(8)确定,损伤演化函数按式(10)进行修正,损伤比例因子按1.3节的方法校准,DIF模型为式(15)和式(16)。SFRC板的主要材料参数列于表1。其中,Lw为裂缝局部宽度,ρ为密度,ν为泊松比,ω为膨胀系数。表1中未给出的参数均采用自动生成值。

    表  1  改进的钢纤维混凝土K&C模型参数
    Table  1.  Parameters of the modified K&C model of SFRC
    Strength surface
    a0/MPa a1 a2/MPa−1 a0y/MPa a1y a2y/MPa−1 a1f a2f/MPa−1
    64.0 0.481 5.82×10−4 45.93 0.726 1.77×10−3 0.476 8.56×10−4
    Damage Others
    b1 b2 b3 α αc αd λm Lw/mm fc/MPa ft/MPa ρ/(kg·m−3) ν ω
    0.75 −1.50 1.15 3.00 0.381 1.90 9.5×10−5 24 175.3 13.8 2 600 0.19 0.5
    下载: 导出CSV 
    | 显示表格
    图  9  SFRC靶板在平头弹丸冲击下的有限元模型
    Figure  9.  Finite element model of flat-ended projectile impacting SFRC plate

    弹丸采用*MAT_JOHNSON_COOK材料模型。该模型适用于应变率变化范围较大及材料软化引起的绝热温度升高问题,是弹丸冲击模拟研究中常用的模型,需耦合*EOS_GRÜNEISEN状态方程。弹丸的主要材料参数如表2所示。其中:G为剪切模量,ABNCM为输入的材料常数,Tm为熔化温度,TR为室内温度,ε为基准静态应变率,cp为定压比热容,SPALL为剥落处理方式,IT为塑性应变迭代选项,D1D2D3D4D5均为失效破坏参数。表2中未给出的参数均采用自动生成值。

    表  2  平头弹丸的材料参数
    Table  2.  Materials parameters of the flat ended projectile
    ρ/(g·cm−3)G/GPaA/MPaB/MPaNCMTm/KTR/K
    7.832107925100.260.0141.031 793294
    ε/(μs−1)cp/(J·kg−1·K−1SPALLITD1D2D3D4D5
    1.0 × 10−44.77 × 10−53.00.04.000.000.000.000.00
    下载: 导出CSV 
    | 显示表格

    弹丸与靶板之间采用侵蚀接触算法*CONTACT_ERODING_SURFACE_TO_SURFACE,动摩擦因子与静摩擦因子均为0.15。采用关键字*INITIAL_VELOCITY_GENERATION产生弹丸初速度。沙漏控制方式采用LS-DYNA软件的标准黏性方式,计算时间步长比例因子为0.6。

    K&C材料模型未包含侵蚀失效准则。为了模拟SFRC板的破坏失效模式,还需引入侵蚀失效准则*MAT_ADD_EROSION,选取失效主应变εmax作为破坏判据,一旦单元主应变大于εmax,该单元将被立即删除。侵蚀失效值依赖于网格尺寸,且因材料而异,并适用于特定研究环境。本研究在选择侵蚀失效值时,主要考虑避免单元的大量缺失和保持质量守恒。如图4中SFRC的峰值应变为0.003,考虑到在应变软化阶段,断裂应变为压缩应变的2倍,还需考虑应变率效应和纤维约束效应等其他影响(建议为10倍),因此失效主应变取为0.06。

    实验结果[1]及改进前后K&C模型预测的SFRC板动力响应的对比如图10所示。由图10可知,原K&C模型的预测结果与实验结果的破坏模式吻合较差,不能预测到实验中的破坏主裂缝,而改进的K&C模型的预测结果与实验结果的破坏模式吻合良好,准确描述了剪切裂缝和拉伸裂缝的发展。SFRC板在弹丸侵彻下的结痂直径与实验结果的比较如表3所示,原K&C模型预测的结痂直径与实验结果偏差较大,而改进的K&C模型精准地预测了弹丸速度为58.2、76.0和104.0 m/s时的结痂直径,与实验结果的相对偏差分别为12.0%、4.0%和6.4%。综上所述,原K&C模型的模拟结果表明,SFRC的脆性破坏特征显著,而改进的K&C模型则准确表征了SFRC的韧性、延性及能量吸收能力,说明改进的K&C模型能更准确地预测SFRC板受冲击荷载时的动力响应。

    表  3  实验数据与数值模拟结痂弹坑直径比较
    Table  3.  Comparisons of experimental data and numerical simulation of the scabbed crater diameter
    Projectile velocity/(m·s−1)Scabbed crater diameter/mmModel error/%
    ExperimentOriginal K&CModified K&COriginal K&CModified K&C
    58.2120.5132106 9.512.0
    76.0119.214812424.2 4.0
    104.0120.3 2812876.7 6.4
    下载: 导出CSV 
    | 显示表格
    图  10  不同撞击速度下SFRC板破坏模式的实验结果[1]与数值模拟结果对比
    Figure  10.  Comparison of the damage patterns of SFRC plate from experimental results[1] and numerical simulation under different impact velocities

    K&C模型常被用来分析爆炸荷载作用下混凝土结构的动力响应,赵春风等[29]采用K&C模型分析安全壳穹顶钢筋混凝土板的抗爆性能。为进一步验证改进的K&C模型在爆炸荷载作用下预测SFRC板动力响应的有效性和可靠性,对SFRC板在爆炸荷载作用下的动力响应[2]进行数值模拟。SFRC板的尺寸设为550 mm × 550 mm × 50 mm,支撑于钢框架之上,自由跨度为460 mm。引爆圆柱形炸药对钢纤维体积分数分别为0.5%、1.0%的SFRC板进行冲击加载,TNT炸药等效当量为0.488 kg,炸药质心距SFRC板表面0.272 5 m,比例距离为0.35 m·kg–1/3

    SFRC的1/4模型如图11所示,在对称面施加对称约束,空气域的非对称面采用关键字*BOUNDARY_NON_REFLECTING施加无反射边界条件,模拟无约束爆炸条件。数值计算模型由钢纤维混凝土板、TNT炸药和空气域组成。SFRC板采用拉格朗日算法,TNT炸药和空气域单元均采用单点多物质ALE算法。另需采用关键字*ALE_MULTI_MATERIAL_GROUP将TNT炸药与空气定义为ALE材料组,它代表一个独特的流体,可与模型中的拉格朗日SFRC结构产生相互作用。采用流固耦合机制*CONSTRAINED_LAGRANGE_IN_SOLID计算SFRC板与空气和炸药之间的相互作用。TNT材料模型采用*MAT_HIGH_EXPLOSIVE_BURN,并联立*EOS_JWL状态方程。空气域采用*MAT_NULL材料模型,并联立*EOS_LINEAR_POLYNOMIAL状态方程。钢纤维混凝土材料模型为改进的K&C模型,材料参数的校准方法及其他环境变量同2.1节。

    图  11  SFRC板在爆炸荷载作用下的有限元模型
    Figure  11.  Finite element model of SFRC plate under blast load

    当模拟爆炸冲击波到达SFRC板时,前表面材料处于三轴压缩应力状态,可导致弹坑状的材料局部破坏。另一部分冲击波透射到SFRC板后表面时,后表面材料处于多轴拉伸应力状态,可导致混凝土的剥落和碎片弹射。实验结果[2]与改进的K&C模型模拟结果的对比如图12所示。由图12可知:改进的K&C模型较好地预测了SFRC板背面混凝土的剥落现象;数值模拟与实验的破坏模式均为典型的对角线破坏,且吻合良好;随着钢纤维体积分数的增加,SFRC板损伤破坏的程度减轻。

    图  12  不同纤维体积分数时SFRC板挠曲破坏模式的实验结果[2]与数值预测对比
    Figure  12.  Comparisons of the damage patterns of SFRC plate from simulation predictions and experimental results[2] under different fiber volume fractions

    通过对原K&C混凝土模型进行改进,研究了钢纤维混凝土板在冲击与爆炸荷载作用下的动力响应,得出以下主要结论。

    (1)根据大量SFRC三轴压缩实验数据,在原失效强度面模型的基础上,考虑了钢纤维体积分数对三轴强度的影响,采用最小二乘法建立了适用于SFRC的失效强度面模型。该模型能准确描述SFRC的失效强度面。

    (2)通过LS-DYNA软件模拟单元轴向拉、压试验,采用反复试验法,建立了SFRC的损伤演化模型,介绍了校准拉、压损伤参数的简便方法,模拟结果表明,新的损伤机制能准确描述SFRC的应变软化段力学特征。

    (3)根据大量高应变率下SFRC的单轴压缩实验数据,在CEB-FIP受压动力增强因子模型的基础上,考虑了钢纤维体积分数对DIF值及过渡应变率的影响,建立了SFRC的CDIF模型,该模型能准确表征SFRC受压时的应变率效应。

    (4)通过LS-DYNA软件模拟了高速弹丸侵彻钢纤维混凝土靶板与爆炸荷载下钢纤维混凝土板的动力响应,模拟结果验证了改进的K&C模型,能够准确表征钢纤维混凝土板受冲击和爆炸荷载作用时的动力响应。

  • 图  (a) 极端高压下各结构相对于I41/amd相固体氢的焓差随压强的变化,(b) 考虑零点振动能后极端高压下各结构相对于oC12相固体氢的焓差随压强的变化

    Figure  1.  (a) Calculated enthalpies of various structures as a function of pressure with respect to I41/amd structure; (b) calculated enthalpies of various structures with the inclusion of ZPE as a function of pressure with respect to the oC12 structure

    图  3.0 TPa下I41/amd相(a)、oC12相(b)(不同颜色的原子代表不同层)、cI16相(c)固体氢的结构和oC12相电子局域函数(d)(均使用VESTA绘制[28]

    Figure  2.  I41/amd structure (a), oC12 structure (b) (spheres of different colors indicate different layers), and cI16 structure (c), and ELF profile of the oC12 structure (d) at 3.0 TPa (All structures and the ELF profile were drawn by using VESTA[28])

    图  不同压强下各结构谱函数α2F(ω)及电声耦合系数λ(ω)

    Figure  3.  Éliashberg spectral function α2F(ω) together with the electron-phonon integral λ(ω) at different pressures

    图  不同压强下各结构的声子谱

    Figure  4.  Phonon dispersion curves of each structure at different pressures

    图  各结构的超导转变温度随压强的变化(μ*=0.10)

    Figure  5.  Superconducting critical temperatures Tc of each structure with μ*=0.10

    表  1  优化后的结构参数

    Table  1.   Optimized structural parameters

    Structure p/TPa Lattice parameters Atomic coordinates (fractional)
    I41/amd 1.0 a=b=1.067 9 Å, c=2.993 9 Å,
    α=β=γ=90°
    H1(4a) 0 1 0
    oC12(Cmcm) 3.0 a=0.947 9 Å, b=2.812 3 Å, c=2.319 0 Å,
    α=β=γ=90°
    H1(8f) 0.500 0 0.140 5 0.580 7
    H2(4c) 0.500 0 0.595 6 1.250 0
    cI16(I¯43d) 4.0 a=b=c=1.929 3 Å, α=β=γ=90° H1(16c) 0.955 5 0.455 5 0.044 5
    下载: 导出CSV

    表  2  不同压强下各结构的超导转变温度及超导计算中的重要参数

    Table  2.   Calculated superconducting parameters and the superconducting transition temperature for I41/amd, oC12, cI16 under different pressures

    Structure p/TPa λ ωlog/K NEF/
    (states·Ry−1/atom)
    Tc/K
    μ*=0.10 μ*=0.13 μ*=0.17
    I41/amd 0.5 1.94 1 784.5 0.23 339 307 280
    1.0 1.87 2 039.0 0.20 382 354 320
    1.5 1.53 2 539.5 0.19 385 348 310
    oC12 2.0 2.65 1 580.2 0.18 418 391 356
    3.0 1.77 2 324.7 0.16 398 356 318
    4.0 1.48 2 822.3 0.14 388 353 312
    cI16 3.0 1.73 2 303.0 0.16 356 345 307
    4.0 1.42 2 782.6 0.14 392 362 322
    5.0 1.27 3 159.0 0.13 295 258 217
    下载: 导出CSV
  • [1] WIGNER E, HUNTINGTON H B. On the possibility of a metallic modification of hydrogen [J]. The Journal of Chemical Physics, 1935, 3(12): 764–770. doi: 10.1063/1.1749590
    [2] PICKARD C J, NEEDS R J. Structure of phase Ⅲ of solid hydrogen [J]. Nature Physics, 2007, 3(7): 473–476. doi: 10.1038/nphys625
    [3] MCMAHON J M, CEPERLEY D M. Ground-state structures of atomic metallic hydrogen [J]. Physical Review Letters, 2011, 106(16): 165302. doi: 10.1103/PhysRevLett.106.165302
    [4] EREMETS M I, DROZDOV A P, KONG P P, et al. Semimetallic molecular hydrogen at pressure above 350 GPa [J]. Nature Physics, 2019, 15(12): 1246–1249. doi: 10.1038/s41567-019-0646-x
    [5] DIAS R P, SILVERA I F. Observation of the Wigner-Huntington transition to metallic hydrogen [J]. Science, 2017, 355(6326): 715–718. doi: 10.1126/science.aal1579
    [6] LOUBEYRE P, OCCELLI F, DUMAS P. Synchrotron infrared spectroscopic evidence of the probable transition to metal hydrogen [J]. Nature, 2020, 577(7792): 631–635. doi: 10.1038/s41586-019-1927-3
    [7] ASHCROFT N W. Metallic hydrogen: a high-temperature superconductor? [J]. Physical Review Letters, 1968, 21(26): 1748–1749. doi: 10.1103/PhysRevLett.21.1748
    [8] CUDAZZO P, PROFETA G, SANNA A, et al. Ab initio description of high-temperature superconductivity in dense molecular hydrogen [J]. Physical Review Letters, 2008, 100(25): 257001. doi: 10.1103/PhysRevLett.100.257001
    [9] YAN Y, GONG J, LIU Y H. Ab initio studies of superconductivity in monatomic metallic hydrogen under high pressure [J]. Physics Letters A, 2011, 375(9): 1264–1268. doi: 10.1016/j.physleta.2011.01.045
    [10] MAKSIMOV E G, SAVRASOV D Y. Lattice stability and superconductivity of the metallic hydrogen at high pressure [J]. Solid State Communications, 2001, 119(10/11): 569–572.
    [11] SZCZȨS̀NIAK R, JAROSIK M W. The superconducting state in metallic hydrogen under pressure at 2 000 GPa [J]. Solid State Communications, 2009, 149(45/46): 2053–2057.
    [12] MCMAHON J M, CEPERLEY D M. High-temperature superconductivity in atomic metallic hydrogen [J]. Physical Review B, 2011, 84(14): 144515. doi: 10.1103/PhysRevB.84.144515
    [13] LIU H Y, WANG H, MA Y M. Quasi-molecular and atomic phases of dense solid hydrogen [J]. The Journal of Physical Chemistry C, 2012, 116(16): 9221–9226. doi: 10.1021/jp301596v
    [14] KRESSE G, FURTHMÜLLER J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set [J]. Physical Review B, 1996, 54(16): 11169–11186. doi: 10.1103/PhysRevB.54.11169
    [15] BLÖCHL P E. Projector augmented-wave method [J]. Physical Review B, 1994, 50(24): 17953–17979. doi: 10.1103/PhysRevB.50.17953
    [16] PERDEW J P, BURKE K, ERNZERHOF M. Generalized gradient approximation made simple [J]. Physical Review Letters, 1996, 77(18): 3865–3868. doi: 10.1103/PhysRevLett.77.3865
    [17] MONKHORST H J, PACK J D. Special points for Brillouin-zone integrations [J]. Physical Review B, 1976, 13(12): 5188–5192. doi: 10.1103/PhysRevB.13.5188
    [18] TOGO A, CHAPUT L, TADANO T, et al. Implementation strategies in phonopy and phono3py [J]. Journal of Physics: Condensed Matter, 2023, 35(35): 353001. doi: 10.1088/1361-648X/acd831
    [19] TOGO A. First-principles phonon calculations with phonopy and phono3py [J]. Journal of the Physical Society of Japan, 2023, 92(1): 012001. doi: 10.7566/JPSJ.92.012001
    [20] GIANNOZZI P, BARONI S, BONINI N, et al. QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials [J]. Journal of Physics: Condensed Matter, 2009, 21(39): 395502. doi: 10.1088/0953-8984/21/39/395502
    [21] MIGDAL A B. Interaction between electrons and lattice vibrations in a normal metal [J]. Soviet Physics JETP, 1958, 34(7): 996–1001.
    [22] ÉLIASHBERG G M. Interactions between electrons and lattice vibrations in a superconductor [J]. Soviet Physics JETP, 1960, 11(3): 696–702.
    [23] SCALAPINO D J, SCHRIEFFER J R, WILKINS J W. Strong-coupling superconductivity.Ⅰ [J]. Physical Review, 1966, 148(1): 263–279. doi: 10.1103/PhysRev.148.263
    [24] VIDBERG H J, SERENE J W. Solving the Eliashberg equations by means of N-point Padé approximants [J]. Journal of Low Temperature Physics, 1977, 29(3/4): 179–192.
    [25] HANFLAND M, SYASSEN K, CHRISTENSEN N E, et al. New high-pressure phases of lithium [J]. Nature, 2000, 408(6809): 174–178. doi: 10.1038/35041515
    [26] GREGORYANZ E, LUNDEGAARD L F, MCMAHON M I, et al. Structural diversity of sodium [J]. Science, 2008, 320(5879): 1054–1057. doi: 10.1126/science.1155715
    [27] MCMAHON M I, GREGORYANZ E, LUNDEGAARD L F, et al. Structure of sodium above 100 GPa by single-crystal X-ray diffraction [J]. Proceedings of the National Academy of Sciences of the United States of America, 2007, 104(44): 17297–17299.
    [28] MOMMA K, IZUMI F. VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data [J]. Journal of Applied Crystallography, 2011, 44(6): 1272–1276. doi: 10.1107/S0021889811038970
    [29] ALLEN P B, DYNES R C. Transition temperature of strong-coupled superconductors reanalyzed [J]. Physical Review B, 1975, 12(3): 905–922. doi: 10.1103/PhysRevB.12.905
    [30] DROZDOV A P, EREMETS M I, TROYAN I A, et al. Conventional superconductivity at 203 kelvin at high pressures in the sulfur hydride system [J]. Nature, 2015, 525(7567): 73–76. doi: 10.1038/nature14964
    [31] MA L, WANG K, XIE Y, et al. High-temperature superconducting phase in clathrate calcium hydride CaH6 up to 215 K at a pressure of 172 GPa [J]. Physical Review Letters, 2022, 128(16): 167001. doi: 10.1103/PhysRevLett.128.167001
    [32] KONG P P, MINKOV V S, KUZOVNIKOV M A, et al. Superconductivity up to 243 K in the yttrium-hydrogen system under high pressure [J]. Nature Communications, 2021, 12(1): 5075. doi: 10.1038/s41467-021-25372-2
    [33] DROZDOV A P, KONG P P, MINKOV V S, et al. Superconductivity at 250 K in lanthanum hydride under high pressures [J]. Nature, 2019, 569(7757): 528–531. doi: 10.1038/s41586-019-1201-8
    [34] SOMAYAZULU M, AHART M, MISHRA A K, et al. Evidence for superconductivity above 260 K in lanthanum superhydride at megabar pressures [J]. Physical Review Letters, 2019, 122(2): 027001. doi: 10.1103/PhysRevLett.122.027001
  • 加载中
图(5) / 表(2)
计量
  • 文章访问数:  265
  • HTML全文浏览量:  78
  • PDF下载量:  56
出版历程
  • 收稿日期:  2023-08-21
  • 修回日期:  2023-10-18
  • 录用日期:  2024-02-26
  • 网络出版日期:  2024-03-25
  • 刊出日期:  2024-04-09

目录

/

返回文章
返回