• ISSN 1001-1455  CN 51-1148/O3
  • EI Compendex、CA收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊
高级检索 E-mail Alert

ZrCuNiAlAg块体非晶合金JH-2模型研究

石永相 施冬梅 李文钊 余志统 尚春明

引用本文:
Citation:

ZrCuNiAlAg块体非晶合金JH-2模型研究

    作者简介: 石永相(1992- ),男,硕士,助理工程师,2536973141@qq.com;
    通讯作者: 施冬梅, 13383013059@163.com
  • 中图分类号: O345

Study on JH-2 model of the ZrCuNiAlAg bulk amorphous alloy

    Corresponding author: SHI Dongmei, 13383013059@163.com ;
  • CLC number: O345

  • 摘要: 为了更好的进行ZrCuNiAlAg块体非晶合金药型罩的爆炸成形及侵彻仿真研究,首要就是建立其材料模型。本文结合ZrCuNiAlAg块体非晶合金力学性能试验结果计算得到了材料的JH-2模型参数,研究确定了ZrCuNiAlAg块体非晶合金JH-2模型。为了验证ZrCuNiAlAg块体非晶合金JH-2模型的准确性,采用Autodyn建立了平板撞击试验有限元模型,模拟了ZrCuNiAlAg块体非晶合金材料在高压、高应变率等环境条件下的变形过程,仿真计算得到的靶板背面自由面粒子速度与试验结果相比,速度平均偏差均在3%以内,表明ZrCuNiAlAg块体非晶合金JH-2模型能很好的描述该材料在大变形、高应变率、高压等环境条件下的力学行为,验证了ZrCuNiAlAg块体非晶合金JH-2模型准确性。
  • 图 1  压缩试样断裂前后

    Figure 1.  Before and after fracture of compressed specimens

    图 2  真实应力-应变曲线

    Figure 2.  Real stress-strain curves

    图 3  动态应力应变曲线

    Figure 3.  Dynamic stress-strain curve

    图 4  飞片及靶板

    Figure 4.  Flyer and target

    图 5  自由面粒子速度历史

    Figure 5.  Velocity history of free surface particles

    图 6  D-u曲线

    Figure 6.  D-u curve

    图 7  p-μ曲线

    Figure 7.  p-μ curve

    图 8  σ1/σ2-p的变化曲线

    Figure 8.  σ1/σ2-p curve

    图 9  有限元模型

    Figure 9.  Finite element model

    图 10  A点速度时间曲线

    Figure 10.  The velocity-time curve of point A

    表 1  SHPB试验结果

    Table 1.  SHPB test results

    气压/MPa应变率/s−1原始尺寸/mm整形片整形片尺寸/mm压缩强度/GPa
    0.176003.84×3.6110×10×0.251.06
    0.2222143.81×3.6410×10×0.251.25
    0.3027833.85×3.6110×10×0.251.46
    0.3531293.82×3.6310×10×0.251.59
    下载: 导出CSV

    表 2  压缩力学性能试验结果

    Table 2.  Test results of compressive mechanical properties

    试验编号$\dot \varepsilon $/s−1σ/GPap/GPaσ*p*
    CMT#14×10−21.230.4070.2110.071
    CMT#24×10−21.240.4090.2130.072
    CMT#34×10−21.240.4160.2130.072
    SHPB#12 2141.250.4180.2150.073
    SHPB#22 7831.460.4890.2510.086
    SHPB#33 1291.590.5290.2730.093
    SHPB#43 1321.640.5710.2820.100
    下载: 导出CSV

    表 3  块体非晶合金碎片准静态压缩试验结果

    Table 3.  Quasi-static compression test results for bulk amorphous alloy fragments

    试验编号$\dot \varepsilon $/s−1σ/GPap/GPaσ*p*
    11×10−21.260.420.2160.074
    22×10−21.410.430.2420.075
    下载: 导出CSV

    表 4  ZrCuNiAlAg块体非晶合金材料模型参数

    Table 4.  Material model parameters

    ρ/(g·cm−3)A1/kPaA2/kPaA3/kPaD1NC
    6.5812.303×1094.716×10108.873×10110.0051.1530.094
    T1/kPaABσHEL/kPaD2MG/kPa
    2.303×1090.2961.035.82×10610.3833.704×106
    下载: 导出CSV

    表 5  计算与试验结果对比

    Table 5.  Comparisons between calculated results and experimental results

    试验编号试验速度/(m·s−1)计算速度/(m·s−1)误差/%
    13503402.9
    23903792.8
    34394272.7
    45024892.6
    55505352.7
    下载: 导出CSV
  • [1] INOUE A, TAKEUCHI A. Recent progress in bulk glassy, nanoquasicrystalline and nanocrystalline alloys [J]. Materials Science & Engineering A, 2004, 43(8): 16–30. DOI: 10.1016/j.msea.2003.10.159.
    [2] SHENG H W, LUO W K, ALAMGIR F M, et al. Atomic packing and short-to-medium-range order in metallic glasses [J]. Nature, 2006, 439(7075): 419–425. DOI: 10.1038/nature04421.
    [3] MIRACLE D B. A structural model for metallic glasses [J]. Nature Materials, 2004, 3(10): 692–702. DOI: 10.1038/nmat1205.
    [4] LI Yaoqi, SONG Min, HE Yuehui. Effect of quenching mode on the mechanical properties of a Zr64Al10Ni15Cu11 metallic glass [J]. Materials and Design, 2010, 31(7): 3555–3558. DOI: 10.1016/j.matdes.2010.01.040.
    [5] 仇在宏, 杨元政, 赵德强, 等. 锆基块体非晶合金的力学性能研究进展 [J]. 金属功能材料, 2004, 11(5): 25–27. DOI: 10.3969/j.issn.1005-8192.2004.05.008.
    QIU Zaihong, YANG Yuanzheng, ZHAO Deqiang, et al. Progress in mechanical properties of Zr-based bulk amorphous alloys [J]. Metallic Functional Materials, 2004, 11(5): 25–27. DOI: 10.3969/j.issn.1005-8192.2004.05.008.
    [6] INOUE A, YOKOYAMA Y, SHINOHARA Y, et al. Preparation of bulk Zr-base amorphous alloy by a zone melting method [J]. Mater Trans, JIM, 1994, 35(2): 923–926. DOI: 10.2320/matertrans1989.35.923.
    [7] 蒋维科. Zr基块体非晶合金的制备及力学性能研究[D]. 兰州: 兰州理工大学, 2013.
    [8] 寇生中, 黄文军, 郑宝超, 等. Y元素添加对Zr基非晶形成和力学性能的影响 [J]. 铸造工程, 2010, 31(11): 1430–1431.
    KOU Shengzhong, HUANG Wenjun, ZHEN Baochao, et al. Glass forming ability and mechanical property of Zr-based alloy with the addition of yttrium [J]. Foundry Technology, 2010, 31(11): 1430–1431.
    [9] 武晓峰, 张海峰, 李宏, 等. Zr基大块非晶合金的微区变形及力学性能 [J]. 中国有色金属学报, 2003, 13(6): 1368–1371. DOI: 10.3321/j.issn:1004-0609.2003.06.008.
    WU Xiaofeng, ZHANG Haifeng, LI Hong, et al. Deformation and mechanical properties of Zr-based bulk metallic glasses under nanoindenter [J]. The Chinese Journal of Nonferrous Metals, 2003, 13(6): 1368–1371. DOI: 10.3321/j.issn:1004-0609.2003.06.008.
    [10] MAO J, ZHANG H F, HUA M F, et al. Effects of casting temperature on mechanical properties of Zr-based metallic glasses [J]. Materials Science and Engineering, 2010, A527: 981–985. DOI: 10.1016/j.msea.2009.09.040.
    [11] 邢秋玮. 锆基非晶合金力学性能影响因素的研究[D]. 兰州: 兰州理工大学, 2013. DOI: 10.7666/d.Y2445918.
    [12] 李杳奇. Zr基块体非晶合金压缩力学性能的研究[D]. 长沙: 中南大学, 2011. DOI: 10.7666/d.y1915234.
    [13] MA G, ZHU Z, WANG Z, et al. Deformation behavior of the Zr53.5Cu26.5Ni5Al12Ag3 bulk metallic glass over a wide range of strain rate and temperatures [J]. Journal of Materials Science & Technology, 2015, 31(9): 941–945. DOI: 10.1016/j.jmst.2015.06.001.
    [14] WANG Weihua, WEI Ping, WANG Limin, et al. Equation of state of bulk metallic glasses studied by an ultrasonic method [J]. Applied Physics Letters, 2001, 79(24): 3947–3949. DOI: 10.1063/1.1426272.
    [15] PAN M X, WANG W H, ZHAO D Q, et al. The equation of state and potential function of Zr41Ti14Cu12.5Ni10Be22.5 bulk metallic glass [J]. Journal of Physics: Condensed Matter, 2002(14): 5665–5671. DOI: 10.1088/0953-8984/14/23/302.
    [16] 潘念侨. Zr基非晶合金材料动态本构关系及其释能效应研究[D]. 南京: 南京理工大学: 2016.
    [17] 石永相, 施冬梅. ZrCuNiAlAg块体非晶合金动静态力学性能研究 [J]. 热加工工艺, 2019, 48(6): 83–86. DOI: 10.14158/j.cnki.1001-3814.2019.06.020.
    SHI Yongxiang, SHI Dongmei. Study on dynamic and static mechanical properties of ZrCuNiAlAg bulk amorphous alloy [J]. Hot Working Technology, 2019, 48(6): 83–86. DOI: 10.14158/j.cnki.1001-3814.2019.06.020.
    [18] 罗兴柏, 张玉令, 丁玉葵.爆炸力学理论教程[M]. 北京:国防工业出版社, 2016: 210.
    [19] 杨震琦, 庞宝君, 王立闻, 等. JH-2模型及其在Al2O3陶瓷低速撞击数值模拟中的应用 [J]. 爆炸与冲击, 2010, 30(5): 464–468. DOI: 10.11883/1001-1455(2010)05-0463-09.
    YANG Zhenqi, PANG Baojun, WANG Liwen, et al. JH-2 model and its application to numericals imulation on Al2O3 ceramic under low-velocity impact [J]. Explosion and Shock Waves, 2010, 30(5): 464–468. DOI: 10.11883/1001-1455(2010)05-0463-09.
  • [1] 杨震琦庞宝君王立闻迟润强 . JH-2模型及其在Al2O3陶瓷低速撞击数值模拟中的应用. 爆炸与冲击, 2010, 30(5): 463-471. doi: 10.11883/1001-1455(2010)05-0463-09
    [2] 蒋东李永池于少娟邓世春 . 钨合金长杆弹侵彻约束AD95陶瓷复合靶. 爆炸与冲击, 2010, 30(1): 91-95. doi: 10.11883/1001-1455(2010)01-0091-05
    [3] 俞鑫炉董新龙潘顺吉 . 不同爆炸载荷下TA2钛合金圆管膨胀破坏过程. 爆炸与冲击, 2018, 38(1): 148-154. doi: 10.11883/bzycj-2017-0014
    [4] 范志庚陈常青万强 . 泡沫铝率相关性能的有限元模拟. 爆炸与冲击, 2014, 34(6): 742-747. doi: 10.11883/1001-1455(2014)06-0742-06
    [5] 谢凡张涛陈继恩刘土光 . 应力三轴度的有限元计算修正. 爆炸与冲击, 2012, 32(1): 8-14. doi: 10.11883/1001-1455(2012)01-0008-07
    [6] 张晓天贾光辉黄海 . 基于节点分离Lagrange有限元方法的超高速碰撞碎片云数值模拟. 爆炸与冲击, 2010, 30(5): 499-504. doi: 10.11883/1001-1455(2010)05-0499-06
    [7] 郑波王安稳 . 直杆碰撞刚性壁弹塑性动力后屈曲有限元分析. 爆炸与冲击, 2007, 27(2): 126-130. doi: 10.11883/1001-1455(2007)02-0126-05
    [8] 王吉王肖钧卞梁 . 光滑粒子法与有限元的耦合算法及其在冲击动力学中的应用. 爆炸与冲击, 2007, 27(6): 522-528. doi: 10.11883/1001-1455(2007)06-0522-07
    [9] 董杰李永池陈学东 . 微孔洞唯象损伤力学模型及其应用. 爆炸与冲击, 2008, 28(5): 443-447. doi: 10.11883/1001-1455(2008)05-0443-05
    [10] 李杰李猛深李宏施存程 . 弹体侵彻干砂的数值模型. 爆炸与冲击, 2015, 35(5): 633-640. doi: 10.11883/1001-1455(2015)05-0633-08
    [11] 李月强衣娜席丰 . 钢柱抗爆响应分析单自由度模型适用性评估. 爆炸与冲击, 2017, 37(5): 957-963. doi: 10.11883/1001-1455(2017)05-0957-07
    [12] 施泽彬朱哲明汪小梦王雄 . Ⅰ型裂纹中低速冲击荷载下起裂韧度测试新方法. 爆炸与冲击, 2018, 38(6): 1247-1254. doi: 10.11883/bzycj-2017-0132
    [13] 刘晶波姚玲 . 有限元畸变单元离散模型中的SH波动. 爆炸与冲击, 1993, 13(1): 9-17.
    [14] 敬霖韩亮亮周彭滔 . 基于SPH方法铁路车轴遭受道砟撞击的数值模拟. 爆炸与冲击, 2018, 38(3): 603-615. doi: 10.11883/bzycj-2016-0265
    [15] 金龙文雷彬李治源张倩 . 轨道炮刨削形成机理分析及数值模拟. 爆炸与冲击, 2013, 33(5): 537-543. doi: 10.11883/1001-1455(2013)05-0537-07
    [16] 王喜军宗智聂春藤井睦雄 . 激光刻痕聚丙烯塑料失效破坏的参数识别. 爆炸与冲击, 2012, 32(5): 457-462. doi: 10.11883/1001-1455(2012)05-0457-06
    [17] 宋顺成才鸿年王富耻 . 关于弹体冲击和贯穿混凝土的三维数值模拟. 爆炸与冲击, 2006, 26(1): 1-6. doi: 10.11883/1001-1455(2006)01-0001-06
    [18] 王军姚熊亮杨棣 . 浮动冲击平台冲击环境对设备响应的影响. 爆炸与冲击, 2015, 35(2): 236-242. doi: 10.11883/1001-1455(2015)02-0236-07
    [19] 李猛深李杰李宏施存程张宁 . 爆炸荷载下钢筋混凝土梁的变形和破坏. 爆炸与冲击, 2015, 35(2): 177-183. doi: 10.11883/1001-1455(2015)02-0177-07
    [20] 李强李鹏辉赵君官茹站勇 . 冲击载荷下自增强身管残余应力变化规律. 爆炸与冲击, 2011, 31(6): 635-640. doi: 10.11883/1001-1455(2011)06-0635-06
  • 加载中
图(10)表(5)
计量
  • 文章访问数:  66
  • HTML全文浏览量:  169
  • PDF下载量:  15
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-06-20
  • 录用日期:  2019-07-04
  • 刊出日期:  2019-09-01

ZrCuNiAlAg块体非晶合金JH-2模型研究

    作者简介:石永相(1992- ),男,硕士,助理工程师,2536973141@qq.com
    通讯作者: 施冬梅, 13383013059@163.com
  • 1. 西北核技术研究所,陕西 西安 710024
  • 2. 陆军工程大学石家庄校区,河北 石家庄 050003

摘要: 为了更好的进行ZrCuNiAlAg块体非晶合金药型罩的爆炸成形及侵彻仿真研究,首要就是建立其材料模型。本文结合ZrCuNiAlAg块体非晶合金力学性能试验结果计算得到了材料的JH-2模型参数,研究确定了ZrCuNiAlAg块体非晶合金JH-2模型。为了验证ZrCuNiAlAg块体非晶合金JH-2模型的准确性,采用Autodyn建立了平板撞击试验有限元模型,模拟了ZrCuNiAlAg块体非晶合金材料在高压、高应变率等环境条件下的变形过程,仿真计算得到的靶板背面自由面粒子速度与试验结果相比,速度平均偏差均在3%以内,表明ZrCuNiAlAg块体非晶合金JH-2模型能很好的描述该材料在大变形、高应变率、高压等环境条件下的力学行为,验证了ZrCuNiAlAg块体非晶合金JH-2模型准确性。

English Abstract

  • Zr基块体非晶合金因为其良好的性能及应用前景已经引起了极大的研究兴趣[1-5],它具备较高的非晶形成能力[6],因而制备出的Zr基非晶合金具备高强度、高硬度、低杨氏模量、高弹性应变和高弹性极限,且还具备良好的超塑性和良好的断裂韧性。

    Zr基块体非晶合金力学性能的影响因素众多,包括组成元素、合成温度、试样规格、加载速率等。例如在ZrCuAl体系中,Fe和Ti微量金属元素的添加降低了体系的力学性能,而Co元素的添加反而提高了体系的力学性能并呈现出一定的塑性[7]。在Zr65-xCu17.5Ni10Al7.5Yx合金中,随着Y含量的增加,合金的热稳定性和力学性能都有所降低[8]。Zr基非晶合金制备过程中,冷却速度越大,维氏硬度HV和弹性模量越小,当退火温度低于玻璃转变温度时,硬度随退火温度的提高而提高[9]。Zr基非晶合金的制备还存在一个合适的温度,一定条件下温度的过高或者过低都会影响其抗压强度[10-11]。试样规格对Zr基非晶合金力学性能的影响主要是指试样尺寸和高径比的影响。一般而言随着试样规格的减小,块体非晶合金的塑性增强,如:李杳奇[12]研究了试样尺寸和高径比对Zr56Al10.9Ni4.6Cu27.8Nb0.7块体非晶合金力学性能的影响,结果表明合金的塑性随着试样尺寸和高径比的减小而增大。试验温度、加载速率等试验条件也会影响Zr基非晶合金的力学性能;Ma等[13]对Zr53.5Cu26.5Ni5Al12Ag3块体非晶合金的力学性能开展了研究,结果表明在应变率为1.0×10−4 s−1时,随着温度从室温升高到593 K时,该非晶合金的断裂强度从1911 MPa降低到了1 690 MPa,当温度增加到673 K时,非晶合金则呈现出均匀变形行为。加载速率的不同会导致Zr基非晶合金处于不同的应变率范围,应变率会对Zr基块体非晶合金的力学性能产生影响。

    为了更好的将该材料应用于武器弹药中,掌握其冲击载荷下的动态力学性能,研究其状态方程、本构模型等就显得特别重要。对于Zr基块体非晶合金状态方程和本构模型的研究已经开展了很多:Wang等[14]研究表明ZrTiCuNiBe块体非晶合金的状态方程适用于Murnaghan形式;Pan等[15]研究表明Zr41Ti14Cu12.5Ni10Be22.5块体非晶合金的P-V曲线适合Bridgman方程;潘念侨[16]采用分离式霍普金森压杆(split Hopkinson pressure bar,SHPB)研究了ZrTiNiCuBe块体非晶合金的本构关系,拟合出了其JH-2模型和Johnson-cook模型,对比发现JH-2模型能更好的描述材料的力学行为。

    本文主要采用万能试验机、SHPB及轻气炮开展ZrCuNiAlAg块体非晶合金材料的动静态力学性能和JH-2模型研究,并应用Autodyn有限元软件模拟ZrCuNiAlAg块体非晶合金平板撞击试验过程,进行JH-2模型的验证。

    • 对于ZrCuNiAlAg块体非晶合金材料的准静态压缩试验前期已经进行过总结,试验采用的试样压缩前后如图1所示,试验通过改变其压缩速度得到不同应变率下准静态压缩试验结果,试验过程及结果详见文献[17]。

      图  1  压缩试样断裂前后

      Figure 1.  Before and after fracture of compressed specimens

      采用MATLAB将万能试验机采集到的力和位移的曲线转换为真实应力和真实应变曲线,如图2所示。

      图  2  真实应力-应变曲线

      Figure 2.  Real stress-strain curves

      ZrCuNiAlAg块体非晶合金的准静态压缩试验中当压缩速度在0.1~10 mm/min变化时,合金的压缩强度变化范围为1.15~1.37GPa,最大与最小强度相差0.22 GPa,平均压缩强度为1.29 GPa。图2中,曲线顶点对应的横坐标值为断裂应变,对应的纵坐标值为断裂应力,曲线斜率为弹性模量,可以发现各个应变率下,应变均随应力的增大而增大,直至试样发生断裂。上述试验结果表明在室温和低应变率条件下ZrCuNiAlAg块体非晶合金具备脆性断裂特征,且无应变率效应。

    • 分离式霍普金森压杆试验采用的试样与准静态压缩试样相同,试验过程中考虑了韧铜整形及马氏体时效钢对压杆端面的防护。分离式霍普金森压杆试验通过改变气室压强来控制撞击杆的撞击速度进而实现不同的加载应变率,不同气压下SHPB试验结果如表1所示。采用Matlab和Origin对原始数据进行处理后得到不同应变率条件下材料的动态应力应变曲线,如图3所示。

      气压/MPa应变率/s−1原始尺寸/mm整形片整形片尺寸/mm压缩强度/GPa
      0.176003.84×3.6110×10×0.251.06
      0.2222143.81×3.6410×10×0.251.25
      0.3027833.85×3.6110×10×0.251.46
      0.3531293.82×3.6310×10×0.251.59

      表 1  SHPB试验结果

      Table 1.  SHPB test results

      图  3  动态应力应变曲线

      Figure 3.  Dynamic stress-strain curve

      整个SHPB试验结果表明:

      (1)不同应变率条件下,随着应变的增加,应力迅速提高直至试件发生断裂,且压缩变形过程中无明显屈服阶段,试件呈脆性断裂特征;

      (2)随着应变率的增大,断裂应力即压缩强度逐渐增大,表明ZrCuNiAlAg块体非晶合金在动态应变率条件下具有明显的应变率强化效应。

    • 采用单级轻气炮开展ZrCuNiAlAg块体非晶合金的平板撞击试验研究,试验原理如图4所示。其中飞片材料为无氧高导电性铜(OFHC),而靶板由靶环、环氧树脂及ZrCuNiAlAg块体非晶合金,尺寸为24 mm×24 mm×4 mm。

      图  4  飞片及靶板

      Figure 4.  Flyer and target

      试验过程中通过采用长短探针和PDV测试技术分别测量飞片的撞击速度及靶板后自由表面速度,而后采用阻抗匹配法,依据冲击波速度与质点速度之间的线性关系[18],将已知飞片材料Hugoniot参数作为标准,将待测材料的冲击压缩特性与标准材料进行比较,获得待测材料的Hugoniot参数。

      试验测得的自由面粒子速度时间曲线如图5所示,由实测的自由面速度历史,可以得到不同撞击速度下材料的Hugoniot弹性极限σHEL,即

      图  5  自由面粒子速度历史

      Figure 5.  Velocity history of free surface particles

      ${\sigma _{\rm{HEL}}} = \rho Cu_{{\rm{fs}}}^{{\rm{HEL}}}/2$

      式中:ρ为密度,C为纵波波速,$u_{{\rm{fs}}}^{\rm {HEL}}$为弹性段自由面速度幅值。

      ZrCuNiAlAg块体非晶合金平板撞击试验结果详见文献[17],其中u为质点速度,D为冲击波速度,拟合得到其D-u型Hugoniot曲线方程为$D = 3 971 + 4.077u$,如图6所示,其中C=3 971 m/s,S=4.077,CS均为材料Hugoniot参数。

      图  6  D-u曲线

      Figure 6.  D-u curve

    • 试验中所采用的ZrCuNiAlAg块体非晶合金试样的平均密度ρ0为6.581 g/cm3,通过准静态试验获得其弹性模量E为100.36 GPa,泊松比为0.35,采用$C = \sqrt {\displaystyle\frac{{\left( {1 - \nu } \right)E}}{{\rho \left( {1 + \nu } \right)\left( {1 - 2\nu } \right)}}} $[18]计算求得纵波波速C为4 947.02 m/s,剪切模量G为37.04 Gpa,体积模量k1为111.67 GPa。

    • 状态方程主要用来描述材料的压力、体积应变和内能之间的关系,JH-2模型中的多项式(Polynomial)状态方程形式为:

      $ p = {A_1}\mu + {A_2}{\mu ^2} + {A_3}{\mu ^3} + ({B_0} + {B_1}\mu ){\rho _0}e \;\;\;\;\;\;\left( {\mu \text{>}0,{\text{压缩}}} \right) $

      $ p = {T_1}\mu + {T_2}{\mu ^2} + {B_0}{\rho _0}e\;\;\;\;\;\;\;\;\;\;\;\quad\quad\;\;\quad\quad\;\left( {\mu \text{<}0,{\text{拉伸}}} \right) $

      式中:μ=ρ/ρ0-1,为体积应变;p为压力;ρ0为平均密度;ρ为当前密度;e为能量;A0A1A2B0B1T1T2均为材料常数。

      考虑到ZrCuNiAlAg块体非晶合金主要应用于高压、高应变率条件下,且不考虑绝热温升导致的热软化等,即忽略内能e的影响,因此适合该合金的多项式(polynomial)状态方程可以简化为:

      ${p_{\rm{r}}} = {A_1}\mu + {A_2}{\mu ^2} + {A_3}{\mu ^3}\quad\quad (\mu \text{>}0) $

      高压下材料的p-μ曲线可通过平板撞击试验获得。平板撞击试验中一维应变下产生的冲击波满足

      ${\rho _0}(D - {u_0}) = \rho (D - u)$

      设初始条件下质点速度${u_0} = 0$,得:

      $\mu = \frac{u}{{D - u}}$

      利用式(5)计算得到μ,描绘(pμ)关系曲线,利用最小二乘法拟合(pμ)的多项式方程,如图7所示。

      图  7  p-μ曲线

      Figure 7.  p-μ curve

      即状态方程为:

      $\begin{split} {p_{\rm{r}}} = &2.303 \times {10^3}\mu + 4.716 \times {10^4}{\mu ^2} +\\ &8.873 \times {10^5}{\mu ^3}\quad \left( {\mu \text{>}0} \right)\end{split} $

      则与之对应的状态方程参数A0A1A2的值分别为2.303×1012、4.716×1013、8.873×1014 Pa。

    • JH-2本构模型描述的等效应力与其等效接触、断裂应力和损伤因子相关,其定义式为[19]

      ${\sigma ^*} = \sigma _{\rm{I}}^* - D(\sigma _{\rm{I}}^* - \sigma _{\rm{F}}^*)\quad\quad(0 \text{≤} D\text{≤} 1) $

      式中:σI*表示等效接触应力;D表示损伤因子(0≤D≤1);σF*表示等效断裂应力。

      当材料未发生损伤(D=0)时,等效接触应力表示为:

      $\sigma _{\rm{I}}^ * = A{({p^ * } + T_{}^ * )^N}[1 + C\ln (\dot \varepsilon / {{\dot\varepsilon _0}} )]$

      当材料完全损伤(D=1)时,等效断裂应力表示为:

      $\sigma _{\rm{F}}^ * = B{({p^ * })^M}[1 + C\ln (\dot \varepsilon / {{\dot\varepsilon _0}} )]$

      式中:ACNBM为材料常数;${\sigma ^ * } = \sigma /{\sigma _{\rm{HEL}}}$${p^ * } = p/{p_{\rm{HEL}}}$${T^ * } = T/{p_{\rm{HEL}}}$σ为真实应力;p为静水压力;T为最大静水拉应力;p*表示为等效静水压力;T*表示为等效最大静水拉应力;σHLE为材料的Hugoniot弹性极限;pHEL为材料处于Hugoniot弹性极限时的压力分量;$\dot\varepsilon$表示为材料真实应变率;$ {{\dot\varepsilon _0}} $表示为参考应变率。

      材料未发生损伤时(D=0),JH-2本构模型对应的材料状态方程为:

      ${p_{\rm{r}}} = 2.303 \times {10^3}\mu + 4.716 \times {10^4}{\mu ^2} + 8.873 \times {10^5}{\mu ^3}$

      根据平板撞击试验结果计算得到材料的σHLE为5.82 GPa。材料的Hugoniot弹性极限σHLE、偏应力张量SHEL满足

      ${\sigma _{\rm{HEL}}} = {p_{\rm{HEL}}} + \frac{2}{3}{S_{\rm{HEL}}}$

      ${S_{\rm{HEL}}} = 2G\frac{{{\mu _{\rm{HEL}}}}}{{1 + {\mu _{\rm{HEL}}}}}$

      联立式(10)~(12)可得式

      ${\sigma _{\rm{HEL}}} = 2.303 \times {10^3}{\mu _{\rm{HEL}}} + 4.716 \times {10^4}\mu _{\rm{HEL}}^2 + 8.873 \times {10^5}\mu _{\rm{HEL}}^3 + \frac{4}{3}G\left(\frac{{{\mu _{\rm{HEL}}}}}{{1 + {\mu _{\rm{HEL}}}}}\right)$

      由式(13)可以求解出${\mu _{\rm{HEL}}} = 0.002\;36$,带入式(10)解得pHEL=5.71 GPa。

      表2为ZrCuNiAlAg块体非晶合金压缩力学性能试验结果。

      试验编号$\dot \varepsilon $/s−1σ/GPap/GPaσ*p*
      CMT#14×10−21.230.4070.2110.071
      CMT#24×10−21.240.4090.2130.072
      CMT#34×10−21.240.4160.2130.072
      SHPB#12 2141.250.4180.2150.073
      SHPB#22 7831.460.4890.2510.086
      SHPB#33 1291.590.5290.2730.093
      SHPB#43 1321.640.5710.2820.100

      表 2  压缩力学性能试验结果

      Table 2.  Test results of compressive mechanical properties

      材料未发生损伤时模型中的应变率效应可以解耦。取3.1×103 s−1和0.04 s−1应变率条件下的(σp)数据,拟合出σ-p关系σ1σ2,将函数σ1σ2相除,得到σ1/σ2p的变化曲线,如图8所示。可以发现当p超过10 GPa时,σ1/σ2函数值接近于常数1.08,将上述参数代入式(8)得到:

      图  8  σ1/σ2-p的变化曲线

      Figure 8.  σ1/σ2-p curve

      $\frac{{{\sigma _1}}}{{{\sigma _2}}} = \frac{{1 + C\ln 3129}}{{1 + C\ln 0.04}} = 1.08$

      计算得到应变率相关系数C=0.094。由计算过程可知,ZrCuNiAlAg块体非晶合金的JH-2模型更适用于高压状态下,压力越高,其材料参数C越准确。

      对于式(9)中的未知数AN,还需要两个独立方程。选取应变率为2 214 s−1,压缩强度为1.25 GPa;另一个选取σHEL,应变率为5 000 s−1,此外假设材料所受的最大静水拉力为0.4 GPa。将数据代入式(9)得到参数A=0.296,N=1.153。

      因此,等效接触应力表示为:

      $\sigma _{\rm{I}}^ * {\rm{ = 0}}{\rm{.296}}{\left( {{p^ * } + {T^ * }} \right)^{1.153}}[1 + 0.094\ln \left( {\dot\varepsilon / {{\dot\varepsilon _0}} } \right)]$

      等效断裂应力表达式中BM材料参数可通过平板撞击后ZrCuNiAlAg块体非晶合金碎片的强度与压力关系来获得。取平板撞击后块体非晶合金碎片开展准静态压缩试验,获得其强度与压力的关系如表3所示。

      试验编号$\dot \varepsilon $/s−1σ/GPap/GPaσ*p*
      11×10−21.260.420.2160.074
      22×10−21.410.430.2420.075

      表 3  块体非晶合金碎片准静态压缩试验结果

      Table 3.  Quasi-static compression test results for bulk amorphous alloy fragments

      表3两组数据带入式(10)可得参数M=0.383、B=1.03。因此,当材料完全损伤(D=1)时,等效接触应力表示为:

      $\sigma _{\rm{F}}^ * {\rm{ = 1}}{\rm{.03}}{\left( {{{{p}}^ * }} \right)^{0.383}}[1 + 0.094\ln \left( {\dot \varepsilon / {{\dot\varepsilon _0}} } \right)]$

      JH-2本构模型中损伤变量D满足:

      $D = \Sigma (\Delta {\varepsilon _{\rm{p}}}/\varepsilon _{\rm{p}}^f)$

      式中$\varepsilon _{\rm{p}}^f $满足:

      $\varepsilon _{\rm{p}}^f = {D_1}({p^ * } + \sigma _{{\rm{t}},{\rm{m}}}^ * )/{D_2}$

      上式D1D2为材料的损伤度函数,由于ZrCuNiAlAg块体非晶合金呈脆性特性,所以采取试验的方法取得D1D2的值难度系数较大,故参考玻璃陶瓷等脆性材料的取值将D1D2分别取为0.005和1[19]

    • 基于ZrCuNiAlAg块体非晶合金的平板撞击试验,采用AUTODYN开展其数值模拟研究,验证ZrCuNiAlAg块体非晶合金的JH-2模型的准确性。为了简化整个数值模拟,忽略固定靶板的靶环及环氧树脂,仅模拟飞片撞击ZrCuNiAlAg块体非晶合金的过程,试验过程中在ZrCuNiAlAg块体非晶合金背面设置可移动高斯点A,如图9所示。

      图  9  有限元模型

      Figure 9.  Finite element model

      有限元模型中,飞片采用AUTODYN材料库中定义的CU-OFHC材料模型,其状态方程和本构模型分别采用SHOCK和Steinberg Guinan模型。靶板材料为ZrCuNiAlAg块体非晶合金,材料模型为自定义的Polynomial状态方程、JH-2模型,其模型参数如表4所示。

      ρ/(g·cm−3)A1/kPaA2/kPaA3/kPaD1NC
      6.5812.303×1094.716×10108.873×10110.0051.1530.094
      T1/kPaABσHEL/kPaD2MG/kPa
      2.303×1090.2961.035.82×10610.3833.704×106

      表 4  ZrCuNiAlAg块体非晶合金材料模型参数

      Table 4.  Material model parameters

      改变飞片的撞击速度,模拟计算得到不同飞片撞击速度下A点速度时间曲线如图10所示,取冲击波过后第一个平台仿真计算结果与其试验结果进行对比,如表5所示。

      图  10  A点速度时间曲线

      Figure 10.  The velocity-time curve of point A

      试验编号试验速度/(m·s−1)计算速度/(m·s−1)误差/%
      13503402.9
      23903792.8
      34394272.7
      45024892.6
      55505352.7

      表 5  计算与试验结果对比

      Table 5.  Comparisons between calculated results and experimental results

      图10可以看出不同速度下仿真计算得到的速度时间曲线变化趋势较为接近,前期上下波动较大,后期逐渐趋为稳定。从表5可以看出计算与试验速度偏差均在3%以内,说明仿真计算得到的结果是准确的,进一步说明了ZrCuNiAlAg块体非晶合金JH-2模型参数的准确性,表明JH-2模型能很好的描述该材料在外界作用下的响应。

    • (1)开展了ZrCuNiAlAg块体非晶合金的动静态力学性能试验研究,研究表明ZrCuNiAlAg块体非晶合金在准静态和动态条件下呈脆性断裂特征,且动态条件下具有明显的应变率强化效应。

      (2)结合力学性能试验结果,计算得到了ZrCuNiAlAg块体非晶合金的JH-2模型参数,研究确定了ZrCuNiAlAg块体非晶合金JH-2模型。

      (3)采用AUTODYN仿真计算得到的靶板背面粒子速度与试验结果相比,平均偏差均在3%以内,表明ZrCuNiAlAg块体非晶合金JH-2模型能很好的描述该材料在大变形、高应变率、高压等环境条件下的力学行为,验证了ZrCuNiAlAg块体非晶合金JH-2模型的准确性。

参考文献 (19)

目录

    /

    返回文章
    返回