Processing math: 100%

高压下HBT晶体的弹性性质

李佐 张凤玲 廖大麟

尹华伟, 蒋轲, 张料, 黄亮, 王陈凌. 钢纤维混凝土板在冲击与爆炸荷载下的K&C模型[J]. 高压物理学报, 2020, 34(3): 034205. doi: 10.11858/gywlxb.20190853
引用本文: 李佐, 张凤玲, 廖大麟. 高压下HBT晶体的弹性性质[J]. 高压物理学报, 2020, 34(4): 041301. doi: 10.11858/gywlxb.20190823
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: LI Zuo, ZHANG Fengling, LIAO Dalin. Elastic Properties of HBT Crystal under High Pressure[J]. Chinese Journal of High Pressure Physics, 2020, 34(4): 041301. doi: 10.11858/gywlxb.20190823

高压下HBT晶体的弹性性质

doi: 10.11858/gywlxb.20190823
基金项目: 贵州省科技厅、毕节市、贵州工程应用技术学院联合基金(LH[2014]7525);贵州省教育厅青年科技人才成长项目(KY[2018]398)
详细信息
    作者简介:

    李 佐(1983-),男,硕士,讲师,主要从事高压下含能材料物性研究. E-mail:lizuo212@163.com

  • 中图分类号: O521.21

Elastic Properties of HBT Crystal under High Pressure

  • 摘要: 采用基于密度泛函理论的第一性原理方法,研究了HBT晶体在常压和高压下的晶格常数和弹性性质及各向异性,利用3种不同理论模型研究了高压下HBT晶体的各向异性性质。研究结果表明:高压下HBT晶体的弹性常数和弹性模量显著增加,晶体表现出高压韧性;同时,高压下HBT晶体具有较大的弹性模量和力学各向异性,随着压强的增大,HBT晶体的各向异性程度减弱。此外,热力学性质计算结果表明,HBT晶体具有较高的德拜温度,并且德拜温度随着压强的增大而升高。

     

  • 普通混凝土具有显著的脆性特征,钢纤维的引入极大地提高了混凝土的延性,并可减少结构部件的剥落飞射,避免了大量的人员伤亡,因此被广泛应用于各类防护结构。近几十年来,随着计算能力的高速发展,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模型,能够准确表征钢纤维混凝土板受冲击和爆炸荷载作用时的动力响应。

  • 图  优化后的HBT晶体几何结构

    Figure  1.  Geometric structure of optimized HBT crystal

    图  理论和实验得到的压强-体积比关系

    Figure  2.  Pressure as a function of volume ratio V/V0 from theoretical and experimental data

    图  弹性常数随压强的变化

    Figure  3.  Variations of elastic constant with pressure

    图  体弹、剪切和杨氏模量随压强的变化关系

    Figure  4.  Variations of bulk modulus B, shear modulus G and Young’s modulus E with pressure

    图  体弹模量和剪切模量之比与压强的关系

    Figure  5.  Pressure dependence of the ratio of bulk modulus and shear modulus

    图  不同压强下HBT晶体的声速和德拜温度

    Figure  6.  Compressional wave velocity, shear wave velocity and averaged wave velocity of HBT crystal under different pressures

    图  模量分数比与压强的关系

    Figure  7.  Pressure dependence of the fractional ratio of modulus

    图  AUAC的各向异性关系

    Figure  8.  Anisotropy diagram of the AU and AC

    图  AL随压强的变化关系

    Figure  9.  Pressure as a function of anisotropy index AL

    图  10  不同压强下剪切模量各向异性的二维投影

    Figure  10.  Anisotropic two-dimensional projection of shear modulus at different pressures

    表  1  HBT晶体晶格参数的计算值和实验值

    Table  1.   Calculated lattice parameters of HBT crystal along with experimental data

    Methoda/nmb/nmc/nmβ/(°)V/nm3
    Exp. [7]1.240 10.551 30.983 5115.570.606 69
    Exp. [21]1.241 50.551 70.984 2115.630.607 79
    This work1.211 20.697 01.038 8113.650.634 54
    下载: 导出CSV
  • [1] 范开敏, 杨莉, 孙庆强, 等. 六角相ErA x (A = H, He)体系弹性性质的第一性原理研究 [J]. 物理学报, 2013, 62(11): 116201. doi: 10.7498/aps.62.116201

    FAN K M, YANG L, SUN Q Q, et al. First-principles study on elastic properties of hexagonal phase ErAx(A = H, He) [J]. Acta Physica Sinica, 2013, 62(11): 116201. doi: 10.7498/aps.62.116201
    [2] ANIRBAN P, CATALIN R P. Contribution of molecular flexibility to the elastic–plastic properties of molecular crystal α-RDX [J]. Modelling and Simulation in Materials Science and Engineering, 2017, 25(1): 015006. doi: 10.1088/1361-651X/25/1/015006
    [3] RYKOUNOV A A. Investigation of the pressure dependent thermodynamic and elastic properties of 1,3,5-triamino-2,4,6-trinitrobenzene using dispersion corrected density functional theory [J]. Journal of Applied Physics, 2015, 117(21): 215901. doi: 10.1063/1.4921815
    [4] LIAN D, LU L Y, WEI D Q. High-pressure behavior of β-HMX crystal studied by DFT-LDA [J]. Chinese Physics Letters, 2008, 25(3): 899–902. doi: 10.1088/0256-307X/25/3/026
    [5] WU Q, ZHU W H, XIAO H M. DFT study on crystalline 1,1-diamino-2,2-dintroethylene under high pressures [J]. Journal of Molecule Modeling, 2013, 19(9): 4039–4047. doi: 10.1007/s00894-013-1931-8
    [6] XU X J, ZHU W H, XIAO H M. DFT studies on the four polymorphs of crystalline CL-20 and the influences of hydrostatic pressure on epsilon-CL-20 crystal [J]. The Journal of Physical Chemistry B, 2007, 111(8): 2090–2097. doi: 10.1021/jp066833e
    [7] KLAPÖTKE T M, SABATE C M. Bistetrazoles: nitrogen-rich, high-performing, insensitive energetic compounds [J]. Chemistry of Materials, 2008, 20(11): 3629–3637. doi: 10.1021/cm703657k
    [8] KLAPÖTKE T M, SABATE C M. 5, 5′-Hydrazinebistetrazole: an oxidation-stable nitrogen-rich compound and starting material for the synthesis of 5, 5′-Azobistetrazolates [J]. Zeitschrift für Anorganische und Allgemeine Chemie, 2007, 633(15): 2671–2677. doi: 10.1002/zaac.200700271
    [9] KARAGHIOSOFF K, KLAPOTKE T M, SABATE C M. Nitrogen-rich compounds in pyrotechnics: alkaline earth metal salts of 5, 5′-Hydrazine-1, 2-diylbis(1H-tetrazole) [J]. European Journal of Inorganic Chemistry, 2009, 2: 238–250.
    [10] EBESPACHER M, KLAPOTKE T M, SABATE C M. Nitrogen-rich alkali metal 5, 5′-Hydrazinebistetrazolate salts: environmentally friendly compounds in pyrotechnic mixtures [J]. New Journal of Chemistry, 2009, 33: 517–527. doi: 10.1039/B818927G
    [11] DE LUCIA F C, GOTTFRIED J L. Characterization of a series of nitrogen-rich molecules using laser induced breakdown spectroscopy [J]. Propellants Explosives Pyrotechnics, 2010, 35(3): 268–277. doi: 10.1002/prep.201000009
    [12] GIANNOZZI P, ANDREUSSI O, BRUMME T, et al. Advanced capabilities for materials modelling with Quantum ESPRESSO [J]. Journal of Physics: Condensed Matter, 2017, 29(46): 465901. doi: 10.1088/1361-648X/aa8f79
    [13] PERDEW J P, BURKE K, ERNZERHOF M. Generalized gradient approximation made simple [J]. Physical Review Letter, 1996, 77(18): 3865–3868. doi: 10.1103/PhysRevLett.77.3865
    [14] PERDEW J P, BURKE K, ERNZERHOF M. Errata: generalized gradient approximation made simple [J]. Physical Review Letter, 1997, 78(7): 1396.
    [15] HAMANN D R. Optimized norm-conserving vanderbilt pseudopotentials [J]. Physical Review B, 2013, 88(8): 085117. doi: 10.1103/PhysRevB.88.085117
    [16] FLETCHER R. Practical methods of optimization [M]. New York: Wiley, 1980: 126−127.
    [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] CORSO D A. Pseudopotentials periodic table: from H to Pu [J]. Computational Materials Science, 2014, 95: 337–350. doi: 10.1016/j.commatsci.2014.07.043
    [19] CORSO D A. Elastic constants of beryllium: a first-principles investigation [J]. Journal of Physics: Condensed Matter, 2016, 28(7): 075401. doi: 10.1088/0953-8984/28/7/075401
    [20] MOUHAT F, COUDERT F X. Necessary and sufficient elastic stability conditions in various crystal systems [J]. Physical Review B, 2014, 90(22): 224104. doi: 10.1103/PhysRevB.90.224104
    [21] CIEZAK-JENKINS J A, JENKINS T A. Shear induced weakening of the hydrogen bonding lattice of the energetic material 5,5′-Hydrazinebistetrazole at high-pressure [J]. Journal of Molecular Structure, 2017, 1129: 313–318. doi: 10.1016/j.molstruc.2016.09.084
    [22] BORN M. On the stability of crystal lattices. I [J]. Mathematical Proceedings of the Cambridge Philosophical Society, 1940, 36(2): 160–172. doi: 10.1017/S0305004100017138
    [23] HILL R. The elastic behaviour of a crystalline aggregate [J]. Proceedings of the Physical Society (Section A), 1952, 65(5): 349–354. doi: 10.1088/0370-1298/65/5/307
    [24] PUGH S F. XCII Relations between the elastic moduli and the plastic properties of polycrystalline pure metals [J]. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 1954, 45(367): 823–843. doi: 10.1080/14786440808520496
    [25] 范航, 聂福德, 龙瑶, 等. 钝感高能炸药三氨基三硝基苯高温高压下热力学性质的分子动力学模拟研究 [J]. 物理学报, 2016, 65(6): 066201. doi: 10.7498/aps.65.066201

    FAN H, NIE F D, LONG Y, et al. A molecular dynamics simulation of thermodynamic properties of 1, 3, 5-triamino-2, 4, 6-trinitrobenzene under high pressure and high temperature [J]. Acta Physica Sinica, 2016, 65(6): 066201. doi: 10.7498/aps.65.066201
    [26] ANDERSON O L. A simplified method for calculating the Debye temperature from elastic constants [J]. Journal of Physics and Chemistry of Solids, 1963, 24(7): 909–917. doi: 10.1016/0022-3697(63)90067-2
    [27] 陈志谦, 林兴香, 李春梅. 硅化锆(ZrSi)弹性及其各向异性研究 [J]. 西南大学学报(自然科学版), 2016, 38(3): 1–8.

    CHEN Z Q, LIN X X, LI C M. Investigations of elasticity anisotropy of ZrSi [J]. Journal of Southwest University (Natural Science Edition), 2016, 38(3): 1–8.
    [28] CHUNG D H, BUESSEM W R. The elastic anisotropy of crystals [J]. Journal of Applied Physics, 1967, 38(5): 2010–2012. doi: 10.1063/1.1709819
    [29] RANGANATHAN S I, OSTOJA-STARZEWSKI M. Universal elastic anisotropy index [J]. Physical Review Letter, 2008, 101(5): 055504. doi: 10.1103/PhysRevLett.101.055504
    [30] TOHER C, PLATA J J, LEVY O, et al. High-throughput computational screening of thermal conductivity, Debye temperature, and Grüneisen parameter using a quasiharmonic Debye model [J]. Physical Review B, 2014, 90(17): 174107. doi: 10.1103/PhysRevB.90.174107
    [31] GAILLAC R, PULLUMBI P, COUDERT F X. ELATE: an open-source online application for analysis and visualization of elastic tensors [J]. Journal of Physics: Condensed Matter, 2016, 28(27): 275201. doi: 10.1088/0953-8984/28/27/275201
  • 加载中
图(10) / 表(1)
计量
  • 文章访问数:  6053
  • HTML全文浏览量:  2535
  • PDF下载量:  33
出版历程
  • 收稿日期:  2019-08-19
  • 修回日期:  2019-09-20

目录

/

返回文章
返回