用于描述W-Cu粉末混合物冲击压缩行为的 p-αp-λ模型适用性

高明悦 周强

引用本文:
Citation:

用于描述W-Cu粉末混合物冲击压缩行为的 p-αp-λ模型适用性

    作者简介: 高明悦(1993-),女,硕士,主要从事混合粉末冲击压缩行为研究. E-mail: 343702781@qq.com;
    通讯作者: 周强, zqpcgm@gmail.com
  • 中图分类号: O347.1; TG301

p-$\alpha $ and p-$\lambda $ Model for Describing Shock Compressive Behavior of W-Cu Powder Mixture

    Corresponding author: ZHOU Qiang, zqpcgm@gmail.com
  • CLC number: O347.1; TG301

  • 摘要: 研究了3种p-$\alpha $模型和p-$\lambda $模型在预测非均质W-Cu混合粉末冲击压实响应的适用性。利用Mie-Grüneisen方法和Barry等压混合法,基于单质W、Cu粉末的Hugoniot关系预测了同孔隙度的W-Cu混合粉末的Hugoniot线,在高压段与实验结果符合较好,但在低压段与实验偏差较大。分别应用3种p-$\alpha $模型和p-$\lambda $模型对实验结果进行拟合,发现除p-$\alpha $ PL模型外,其他模型均较好地描述W-Cu混合粉末的冲击压缩响应,受经验参数选择的影响,所有模型的压溃强度和压缩路径各不相同,预测功能较差。
  • 图 1  粉末或多孔材料低压区域冲击压缩行为示意图

    Figure 1.  Compressive behavior of powder or porous material

    图 2  零件拆分及装配图

    Figure 2.  Parts split diagram and assembly diagram

    图 3  探针位置

    Figure 3.  Probe location

    图 4  典型的入射波(a)和传播波(b)的DISAR轮廓曲线(显示了结构化波形,以及50%的均衡到达时间和10%、90%的不确定到达时间,用于计算冲击波速度)

    Figure 4.  Typical extracted profile for input (a) and propagated waves (b) that illustrating structured waveforms (Locations of 50% equilibrium time of arrival and 10% and 90% uncertainty arrival times used for calculation of shock velocity are also marked.)

    图 5  计算原理和误差分布

    Figure 5.  Calculation principle and error distribution

    图 6  4种模型拟合结果

    Figure 6.  Fitting results of four models

    表 1  实验相关参数计算结果及误差范围

    Table 1.  Calculated results and errors of the related experimental parameters

    No.${\rho _{00}}$/(g·cm–3)vI/ (km·s–1)uS/(km·s–1)uP/(km·s–1)p/GPa$\rho $/(g·cm–3)
    110.696±0.9100.622±0.0031.254±0.0410.457±0.0196.126±0.43416.826±1.122
    210.484±0.7310.270±0.0010.913±0.0550.213±0.0132.038±0.22313.668±0.536
    310.707±1.2290.784±0.0041.549±0.0890.546±0.0289.056±0.91316.639±1.324
    410.237±1.3700.508±0.0031.107±0.0300.388±0.0164.392±0.21115.752±1.370
    下载: 导出CSV

    表 2  组分和混合物的相关冲击和材料特性

    Table 2.  Relevant shock and material properties of constituents and mixtures

    Material${\rho _0}$/(g·cm–3)C0/(km·s–1)S${\gamma _0}$V0/(cm3·g–1)Y/GPa
    W19.35 4.0641.2041.780.052 01.700
    Cu8.933.9101.5101.970.112 10.065
    W-Cu15.1174.0271.2771.830.066 41.030
    下载: 导出CSV

    表 3  模型拟合参数

    Table 3.  Model fitting parameters

    ModelParameters
    pS/GPapE/GPa${ \alpha _{\rm{E} } }$N
    p-${\alpha} $(MQ)4.600.057 71.3591.952
    p-${\alpha} $(PL)pS/GPan
    4.607.80
    p-${\alpha} $(SS)Y/GPa
    1.273 4
    p-$\lambda $n
    1.51
    下载: 导出CSV
  • [1] ALTSHULER L V, KRUPNIKOV K K, LEDENEV B N, et al. Dynamic compressibility and equation of state for iron under high pressure [J]. Soviet Physics-JETP, 1958: 34.
    [2] MCQUEEN R G, MARSH S P, TAYLOR J W, et al. The equation of state of solids from shock wave studies [M]. Los Alamos, New Mexico: University of California, 1970: 293–417.
    [3] FREDENBURG D A. Shock compaction and impact response of thermit powder mixtures [D]. Georgia: Georgia Institute of Technology, 2010.
    [4] HERRMANN W. Constitutive equation for the dynamic compaction of ductile porous materials [J]. Journal of Applied Physics, 1969, 40(6): 2490–2499. doi: 10.1063/1.1658021
    [5] DAI C D, EAKINS D E, THADHANI N N. Dynamic densification behavior of nanoiron powders under shock compression [J]. Journal of Applied Physics, 2008, 103(9): 093503. doi: 10.1063/1.2908209
    [6] BUTCHER B M, KARNES C H. Dynamic compaction of porous iron [J]. Journal of Applied Physics, 1969, 40(7): 2967–2976. doi: 10.1063/1.1658109
    [7] CARROLL M M, HOLT A C. Static and dynamic pore‐collapse relations for ductile porous materials [J]. Journal of Applied Physics, 1972, 43(4): 1626–1636. doi: 10.1063/1.1661372
    [8] GRADY D, KERLEY E G I, KUHNS L D, et al. Computational modeling and wave propagation in media with inelastic deforming microstructure [J]. Journal de Physique, IV: Proceedings of International Conference, 2000, 10(9): 15–20.
    [9] MEYERS M A. Shock waves: equations of state [M]. John Wiley & Sons, Inc., 2007.
    [10] MCQUEEN R G, MARSH S P. Equation of state for nineteen metallic elements from shock‐wave measurements to two megabars [J]. Journal of Applied Physics, 1960, 31(7): 1253–1269. doi: 10.1063/1.1735815
    [11] ALEKSEEV Y F, AL’TSHULER L V, KRUPNIKOVA V P. Shock compression of two-component paraffin-tungsten mixtures [J]. Journal of Applied Mechanics Technical Physics, 1971, 12(4): 624–627.
    [12] BATSANOV S S. Effects of explosions on materials [M]. New York: Springer, 1994.
    [13] MEYERS M A. Dynamic Behavior of Materials [M]. San Diego: University of California, 1994.
    [14] KRUEGER B R, MUTZ A H, VREELAND T. Correlation of shock initiated and thermally initiated chemical reactions in a 1∶1 atomic ratio nickel‐silicon mixture [J]. Journal of Applied Physics, 1991, 70(10): 5362–5368. doi: 10.1063/1.350217
    [15] WENG J, TAN H, WANG X, et al. Optical-fiber interferometer for velocity measurements with picosecond resolution [J]. Applied Physics Letters, 2006, 89(11): 111101. doi: 10.1063/1.2335948
    [16] FREDENBURG D A, KOLLER D D, RIGG P A, et al. High-fidelity Hugoniot analysis of porous materials [J]. Review of Scientific Instruments, 2013, 84(1): 013903. doi: 10.1063/1.4774394
    [17] MITCHELL A C, NELLIS W J. Shock compression of aluminum, copper, and tantalum [J]. Journal of Applied Physics, 1981, 52(5): 3363–3374. doi: 10.1063/1.329160
    [18] FREDENBURG D A, THADHANI N N. On the applicability of the P- $\alpha $ and P- $\lambda $ models to describe the dynamic compaction response of highly heterogeneous powder mixtures [J]. Journal of Applied Physics, 2013, 113(4): 043507. doi: 10.1063/1.4788700
    [19] BROWN J L, VOGLER T J, GRADY D E, et al. Dynamic compaction of sand [C]//Shock Compression of Condensed Matter-2007, 2007: 1363–1366.
    [20] NEEL C H. Shock compression of a heterogeneous, porous polymer composite [J]. Dissertations & Theses-Gradworks, 2010.
  • [1] 毕延谭华经福谦赵敏光 . 冲击压缩下金属的电导率测量. 高压物理学报, 2003, 17(1): 1-7 . doi: 10.11858/gywlxb.2003.01.001
    [2] 刘高旻杜金梅刘雨生谭华贺红亮 . PZT 95/5铁电陶瓷的冲击压缩Hugoniot特性研究. 高压物理学报, 2008, 22(1): 30-34 . doi: 10.11858/gywlxb.2008.01.007
    [3] 汤文辉张若棋经福谦胡金彪 . 冲击压缩下金属/窗口界面热弛豫过程的热阻模型研究. 高压物理学报, 1993, 7(4): 247-253 . doi: 10.11858/gywlxb.1993.04.002
    [4] 陈源福刘福生张宁超赵北京王军国张明建薛学东 . 瞬态激光拉曼光谱测量系统及其在苯冲击压缩实验中的应用. 高压物理学报, 2013, 27(4): 505-510. doi: 10.11858/gywlxb.2013.04.006
    [5] 孙悦孙东帅俊卿吴国栋施尚春 . 液态水、苯和壬基乳化剂混合物冲击压缩特性的实验研究. 高压物理学报, 2008, 22(1): 62-66 . doi: 10.11858/gywlxb.2008.01.014
    [6] 孙悦孟川民吴国栋杨向东 . 液氮与液态一氧化碳混合物的冲击压缩特性研究. 高压物理学报, 2002, 16(3): 217-223 . doi: 10.11858/gywlxb.2002.03.010
    [7] 胡栋杨向东经福谦王永国李孝昌孙珠妹董石施尚春 . 碳水混合物冲击压缩特性研究. 高压物理学报, 1997, 11(3): 197-202 . doi: 10.11858/gywlxb.1997.03.006
    [8] 古成钢谢盘海王金贵经福谦 . 测量冲击压缩下绝缘材料电阻率技术的一个改进. 高压物理学报, 1988, 2(4): 340-345 . doi: 10.11858/gywlxb.1988.04.008
    [9] 彭克锋潘昊赵凯郑志军虞吉林 . 铜粉末动态压缩行为的多颗粒有限元分析. 高压物理学报, 2019, 33(4): 044102-1-044102-8. doi: 10.11858/gywlxb.20180665
    [10] 汤文辉张若棋胡金彪经福谦 . 金属冲击温度测量中金属/窗口界面表观光谱辐亮度的再研究. 高压物理学报, 1995, 9(4): 264-269 . doi: 10.11858/gywlxb.1995.04.004
    [11] 林华令于万瑞 . 冲击压缩后热传导的理论研究. 高压物理学报, 1994, 8(1): 49-56 . doi: 10.11858/gywlxb.1994.01.008
    [12] 郝龙王翔王青松康强黄金 . 冲击压缩下PMMA的响应和光学特性. 高压物理学报, 2017, 31(5): 579-584. doi: 10.11858/gywlxb.2017.05.011
    [13] 乔良张先锋何勇史安顺张将 . 颗粒金属材料冲击压缩细观数值模拟. 高压物理学报, 2013, 27(6): 863-871. doi: 10.11858/gywlxb.2013.06.011
    [14] 孟川民焦荣珍施尚春董石孙悦杨向东经福谦 . 液氦冲击压缩特性的理论计算. 高压物理学报, 2002, 16(1): 70-74 . doi: 10.11858/gywlxb.2002.01.012
    [15] 胡金彪谭华唐志平杨嘉陵 . 冲击压缩下NaCl单晶光辐射特性的研究. 高压物理学报, 2001, 15(4): 277-284 . doi: 10.11858/gywlxb.2001.04.007
    [16] 施尚春董石黄跃刘福生孙悦 . 液态N2、CO冲击压缩特性研究. 高压物理学报, 1999, 13(4): 295-300 . doi: 10.11858/gywlxb.1999.04.010
    [17] 王金贵李欣竹王为傅秋卫王翔 . 三种钨合金的冲击压缩性. 高压物理学报, 1998, 12(4): 258-263 . doi: 10.11858/gywlxb.1998.04.004
    [18] 李巧燕施尚春杨金文孙悦 . 钕铁硼的冲击压缩特性. 高压物理学报, 2007, 21(2): 210-214 . doi: 10.11858/gywlxb.2007.02.016
    [19] 张修路袁长迎孙悦韦建军董石吴国栋 . 壬基酚聚氧乙烯醚的冲击压缩实验研究. 高压物理学报, 2005, 19(2): 164-168 . doi: 10.11858/gywlxb.2005.02.011
    [20] 李泽仁刘俊程晓锋马如超姚建铨 . 冲击压缩下钨靶对514.5 nm激光的漫反射率. 高压物理学报, 1999, 13(1): 71-75 . doi: 10.11858/gywlxb.1999.01.013
  • 加载中
图(6)表(3)
计量
  • 文章访问数:  256
  • 阅读全文浏览量:  246
  • PDF下载量:  0
出版历程
  • 收稿日期:  2019-05-27
  • 录用日期:  2019-06-11
  • 网络出版日期:  2020-01-17
  • 刊出日期:  2020-02-01

用于描述W-Cu粉末混合物冲击压缩行为的 p-αp-λ模型适用性

    作者简介:高明悦(1993-),女,硕士,主要从事混合粉末冲击压缩行为研究. E-mail: 343702781@qq.com
    通讯作者: 周强, zqpcgm@gmail.com
  • 北京理工大学爆炸科学与技术国家重点实验室,北京 100081

摘要: 研究了3种p-\begin{document}$\alpha $\end{document}模型和p-\begin{document}$\lambda $\end{document}模型在预测非均质W-Cu混合粉末冲击压实响应的适用性。利用Mie-Grüneisen方法和Barry等压混合法,基于单质W、Cu粉末的Hugoniot关系预测了同孔隙度的W-Cu混合粉末的Hugoniot线,在高压段与实验结果符合较好,但在低压段与实验偏差较大。分别应用3种p-\begin{document}$\alpha $\end{document}模型和p-\begin{document}$\lambda $\end{document}模型对实验结果进行拟合,发现除p-\begin{document}$\alpha $\end{document} PL模型外,其他模型均较好地描述W-Cu混合粉末的冲击压缩响应,受经验参数选择的影响,所有模型的压溃强度和压缩路径各不相同,预测功能较差。

English Abstract

  • 多孔材料和粉末材料的冲击压缩一直是重点研究方向。早期研究主要通过多孔或粉末材料的冲击压缩,获得在固体材料中无法获得的热力学状态,以研究完全物态方程,此类研究关注高压区域,忽略粉末材料的自身强度[1-2]。粉末的冲击压实利用冲击波作用于非晶或纳米晶粉末,在瞬态高温高压下得到近乎致密的块体,且仍保持原始晶粒尺寸,是制备亚稳态材料的一种非常有潜力和前景的方法。为了避免高压导致的高温和裂纹问题,往往在1~10 GPa的低压范围内实现压实,因此,粉末在低压区的致密化行为对于粉末冲击压实研究非常重要。

    粉末在低压区的冲击压实是一个复杂的过程,受到多种因素的影响,包括粉末材料的弹塑性、初始颗粒的大小和形状等[3]。为了在较宽的压力范围内更好地描述粉末的压实过程,一般将这一过程分解为弹性段和塑性段,如图1所示,其中:${\alpha _{\rm{E}}} $pE为弹性段结束时的孔隙度和压力,由于在整个压缩过程中,弹性段仅占很小部分,一般都忽略不计,多数模型取${\alpha _{\rm{E}}} ={\alpha _0} $;当压力超过pE后,进入塑性段,随着颗粒形变和孔隙坍缩,孔隙率逐渐降低,压力达到pS时,粉末完全压实,pS也被称为压溃强度。目前已经有较多模型描述粉末压实响应,但存在两个假设前提:(1)忽略了剪切强度;(2)在同等压力、温度条件下,粉末和固体材料的比内能是相等的。假设(1)使得模型只能描述由于孔隙闭合导致的不可逆压实;假设(2)忽略了孔隙和颗粒的表面能,更适合描述颗粒较大的粉末体系。

    图  1  粉末或多孔材料低压区域冲击压缩行为示意图

    Figure 1.  Compressive behavior of powder or porous material

    在众多描述颗粒或粉末材料动态压实行为的模型中,以p-${\alpha} $模型[4-7]p-$\lambda $模型[8]应用最为广泛,是典型的描述型模型,这类模型结合了粉末材料本征特性,通过实验数据拟合获得模型的关键参数。尽管模型的具体形式不同,但都包含一个或者多个经验拟合参数。通常情况下这些参数与可测量的物理参数没有直接关系,只是为了使模型更好地与实验数据匹配,这种灵活性使模型可以很好地描述某种特定粉末系统的冲击压实过程。

    对于由不同性质粉末组成的混合物,粉末冲击压实过程变得更加复杂。目前并没有有效的方法建立混合粉末冲击压缩方程,一般可通过0 K混合法则[9-10]或质量分数加权平均法则[11-12]先得到混合固体Hugoniot线,再经等容Mie-Grüneisen方法计算得到粉体的冲击压缩方程,但这种方法的计算过程繁琐,需要专门程序计算,而且其有效性缺乏实验验证。

    本研究测量了W-Cu混合粉末的Hugoniot数据,系统分析了测量误差;然后分别应用p-${\alpha} $p-$\lambda $等模型对测量数据进行拟合;利用Mie-Grüneisen方法[13]得到同孔隙度的W、Cu单质粉末的Hugoniot线,按照Barry等压混合法则[14]预测同孔隙度的W-Cu混合粉末的Hugoniot线,并与测量得到的数据进行对比,探讨p-${\alpha} $p-$\lambda $等模型在描述混合粉末冲击压缩行为时的适用性和经单质粉末Hugoniot线预测混合粉末Hugoniot线的可行性。

    • 实验使用的W和Cu单质粉末均为市购,纯度均为99.9%,粉末粒度均为 $1 \;{\text{μ}}{\rm{m}}$。W、Cu粉按照质量比76∶24混合,经V型混料机混合均匀得到W-Cu混合粉末。图2为实验测量所用装置的结构示意图,将无氧铜基板和约束环依次安装于中空的304不锈钢套筒后,形成一个容器,再依次装入W-Cu混合粉末和蓝宝石窗口,并通过背板确保粉末后表面与窗口材料紧密接触。实验使用二级轻气炮进行加载,飞片使用与基板同材质的无氧铜,应用光纤位移干涉仪(Displacement interferometer system for any reflector,DISAR)[15],通过基板背面和窗口背面设置的4个光纤探头,分别记录冲击波抵达粉末样品上下表面的时间,探头布置的空间位置如图3所示。

      图  2  零件拆分及装配图

      Figure 2.  Parts split diagram and assembly diagram

      图  3  探针位置

      Figure 3.  Probe location

    • 粉末初始密度${\rho _{00}}$主要通过调整粉体厚度来控制,粉末厚度由组装件整体厚度减去基板中心、蓝宝石窗口及背板厚度计算得到。厚度测量精度为10–4 mm,粉末质量称量精度为10–3 g。每组取多次测量的平均值,并在初始密度的不确定度分析中加入了每个位置测量的方差。粉末中存在局部的不均匀性,但无法确定不均匀部位的位置和范围,因此初始孔隙密度${\rho _{00}}$的测量误差取最大值,即直接取粉体厚度测量的最大值和最小值来确定。初始密度的测量误差会在阻抗匹配计算中传递,对Hugoniot参数(uP, pV)的测量误差影响较大。

    • 飞片速度vI由炮口的电磁线圈测量得到,该系统的测量误差为0.5%。由于飞片在炮膛中会发生偏转,形成的倾角可由图3中1、2号探针确定后扣除。

    • 冲击波速度uS根据冲击波分别到达图3中1、2号探针和3、4号探针记录的时间和粉末厚度计算得出,即${u_{\rm{S}}} = {\text{∆}} {x_i}/{\text{∆}} {t_i}$,其中${\text{∆}} {x_i}$${\text{∆}} {t_i}$分别是粉末厚度和冲击波传播时间。粉末厚度的标准差为${\text{δ}} \left( {{\text{∆}} {x_i}} \right) = {\left( {{\text{δ}} x_{\rm{B}}^2 + {\text{δ}} x_{\rm{W}}^2} \right)^{1/2}}$,其中xB表示基板位移,xW表示粉末厚度[16]${\text{∆}} {t_i}$即冲击波分别抵达1、2号探针和3、4号探针的时间差,其标准差为${\text{δ}} \left( {{\text{∆}} {t_i}} \right)$图4所示为粉体上下表面处探针的DISAR信号。由于飞片飞行过程中会发生偏转,1、2号探针信号存在时间差,也会导致冲击波速度测量误差,这一误差也体现在${\text{δ}} \left( {{\text{∆}} {t_i}} \right)$中。在两个标准差范围内一个实验的平均冲击速度总误差为

      图  4  典型的入射波(a)和传播波(b)的DISAR轮廓曲线(显示了结构化波形,以及50%的均衡到达时间和10%、90%的不确定到达时间,用于计算冲击波速度)

      Figure 4.  Typical extracted profile for input (a) and propagated waves (b) that illustrating structured waveforms (Locations of 50% equilibrium time of arrival and 10% and 90% uncertainty arrival times used for calculation of shock velocity are also marked.)

      $ \frac{{{\text{δ}}{u_{\rm{S}}}}}{{{u_{\rm{S}}}}} = {\left[ {{\sum_i}2{{\left[ {\frac{{{\text{δ}}\left( {{\text{∆}} {x_i}} \right)}}{{{\text{∆}} {x_i}}}} \right]}^2} + {\sum_i}2{{\left( {\frac{{{\text{δ}}\left( {{\text{∆}} {t_i}} \right)}}{{{\text{∆}} {t_i}}}} \right)}^2}} \right]^{\frac{1}{2}}} $

    • 粉末的粒子速度uP使用阻抗匹配法计算,计算原理和误差分布如图5所示。根据文献[17],粒子速度的总误差是系统误差和实验误差之和,即${\text{δ}}{u_{\rm{P}}} = {\text{δ}} u_{\rm{P}}^{\rm{E}} + {\text{δ}} u_{\rm{P}}^{\rm{S}}$。粉末粒子速度的系统误差是

      图  5  计算原理和误差分布

      Figure 5.  Calculation principle and error distribution

      $ {\text{δ}}{u_{\rm{P}}} = {\left[ {{{\left( {{\text{δ}}u_{{\rm{P}},{\rm{I}} - {\rm{B}}}^{\rm{S}}} \right)}^2} + {{\left( {{\text{δ}}u_{\rm{P}}^{\rm{S}}} \right)}^2}} \right]^{1/2}} $

      式中:$u_{{\rm{P}},{\rm{I}} - {\rm{B}}}^{\rm{S}}$表示靶板Hugoniot误差,$u_{\rm P}^{\rm S} $表示靶板等熵线上的误差。实验误差为

      $ {\text{δ}}u_{\rm{P}}^{\rm{E}} = {\left\{{{{\left[{{\text{δ}}{u_{\rm{P}}}\left( {{u_{\rm{S}}}} \right)} \right]}^2} + {{\left[{{\text{δ}}{u_{\rm{P}}}\left( {{v_{\rm I}}} \right)} \right]}^2} + {{\left[{{\text{δ}}{u_{\rm{P}}}\left( {{\rho _{00}}} \right)} \right]}^2} + {{\left[{{\text{δ}}{u_{\rm{P}}}\left( {{\rho _{0{\rm{I}}}}} \right)} \right]}^2}{{\left[{{\text{δ}}{u_{\rm{P}}}\left( {{\rho _{0{\rm{B}}}}} \right)} \right]}^2}} \right\}^{1/2}} $

      式中:${\rho _{0{\rm{I}}}}$${\rho _{0{\rm{B}}}}$分别是飞片和基板的初始密度。

    • 压力的初始值p0和粉末初始速度u0等于零,则粉末中的压力p可直接由下式计算

      $ p = {\rho _0}\left( {{u_{\rm{S}}} - {u_0}} \right){u_{\rm{P}}} $

      压力的误差${\text{δ}} p/p$也包括系统误差和实验误差。由于p${\rho _{00}}$uSuP的函数,因此,压力p的系统误差为

      $ \frac{{{{\text{δ}}\sigma }}}{\sigma } = \frac{{{\text{δ}} {\rho _{00}}}}{{{\rho _{00}}}} + \frac{{{\text{δ}} {u_{\rm{S}}}}}{{{u_{\rm{S}}}}} + \frac{{{\text{δ}} {u_{\rm{P}}}}}{{{u_{\rm{P}}}}} $

      实验误差为[17]

      $ {\left( {\frac{{{\text{δ}} \sigma }}{\sigma }} \right)_{\rm{E}}} = {\left\{{\begin{aligned} & {{{\left( {1 + \frac{{{u_{\rm{S}}}}}{{{u_{\rm{P}}}}}\frac{{\partial {u_{\rm{P}}}}}{{\partial {u_{\rm{S}}}}}} \right)}^2}{{\left( {\frac{{\partial {u_{\rm{S}}}}}{{{u_{\rm{S}}}}}} \right)}^2} + {{\left( {1 + \frac{{{\rho _{00}}}}{{{u_{\rm{P}}}}}\frac{{\partial {u_{\rm{P}}}}}{{\partial {\rho _{00}}}}} \right)}^2}{{\left( {\frac{{\partial {\rho _{00}}}}{{{\rho _{00}}}}} \right)}^2}}+ {{\left[ {\frac{{\partial {u_{\rm{P}}}\left( {{v_1}} \right)}}{{{u_{\rm{P}}}}}} \right]}^2} +\\ & {\left( {\frac{{\partial {u_{\rm{P}}}\left( {{\rho _{\rm {0I}}}} \right)}}{{{u_{\rm{P}}}}}} \right){{\left[ {\frac{{\partial {u_{\rm{P}}}\left( {{\rho _{\rm {0I}}}} \right)}}{{{u_{\rm{P}}}}}} \right]}^2} + \left( {\frac{{\partial {u_{\rm{P}}}\left( {{\rho _{\rm {0B}}}} \right)}}{{{u_{\rm{P}}}}}} \right){{\left[ {\frac{{\partial {u_{\rm{P}}}\left( {{\rho _{\rm {0B}}}} \right)}}{{{u_{\rm{P}}}}}} \right]}^2}} \end{aligned}} \right\}^{\textstyle\frac{1}{2}}} $

    • 粉末中的冲击压缩比容V$V = {V_0}\left[ {1 - \left( {{u_{\rm{P}}} - {u_0}} \right)/\left( {{u_{\rm{S}}} - {u_0}} \right)} \right]$计算得到,其中$V = 1/\rho $V的系统误差为${\left( {{\text{δ}}V/V} \right)_{\rm{S}}} = \left( {\eta - 1} \right)\left| {{\rm{\text{δ}}}u_{\rm{P}}^{\rm{S}}} \right|/{u_{\rm{P}}}$,其中$\eta $是粉末压缩度,$\eta = \rho /{\rho _{00}} $。实验误差是

      $ {\left( {\frac{{{\text{δ}}V}}{V}} \right)_{\rm E}} = \left\{ {{{\left( {\eta - 1} \right)}^2}{{\left( {1 - \frac{{{u_{\rm{S}}}}}{{{u_{\rm{P}}}}}\frac{{\partial {u_{\rm{P}}}}}{{\partial {u_{\rm{S}}}}}} \right)}^2}{{\left( {\frac{{\partial {u_{\rm{S}}}}}{{{u_{\rm{S}}}}}} \right)}^2} + {{\left[ {1 + \left( {\eta - 1} \right)\frac{{{\rho _{00}}}}{{{u_{\rm{P}}}}}\frac{{\partial {u_{\rm{P}}}}}{{\partial {\rho _{00}}}}} \right]}^2}{{\left( {\frac{{\partial {\rho _{00}}}}{{{\rho _{00}}}}} \right)}^2} + {{\left( {\eta - 1} \right)}^2}{F^2}} \right\} $

      式中:${F^2} = {\left[ {\dfrac{{\partial {u_{\rm P}}\left( {{V_1}} \right)}}{{{u_{\rm P}}}}} \right]^2} + {\left[ {\dfrac{{\partial {u_{\rm P}}\left( {{\rho _{\rm 0I}}} \right)}}{{{u_{\rm P}}}}} \right]^2} + {\left[ {\dfrac{{\partial {u_{\rm P}}\left( {{\rho _{\rm 0B}}} \right)}}{{{u_{\rm P}}}}} \right]^2}。$

    • 实验测量的所有Hugoniot参数及其误差范围如表1所示。可以看到uS的误差较大,主要是粉体厚度测量误差较大所致,也导致粉末初始密度和压缩密度误差均接近10%。压力超过4 GPa后,粉体已接近致密。

      No.${\rho _{00}}$/(g·cm–3)vI/ (km·s–1)uS/(km·s–1)uP/(km·s–1)p/GPa$\rho $/(g·cm–3)
      110.696±0.9100.622±0.0031.254±0.0410.457±0.0196.126±0.43416.826±1.122
      210.484±0.7310.270±0.0010.913±0.0550.213±0.0132.038±0.22313.668±0.536
      310.707±1.2290.784±0.0041.549±0.0890.546±0.0289.056±0.91316.639±1.324
      410.237±1.3700.508±0.0031.107±0.0300.388±0.0164.392±0.21115.752±1.370

      表 1  实验相关参数计算结果及误差范围

      Table 1.  Calculated results and errors of the related experimental parameters

    • p-${\alpha} $模型被广泛用于描述均质和非均质粉末的弹塑性压实响应[18]。其中${\alpha} =V/{V_{\rm S}}$是一个描述粉末膨胀程度的参数,V是粉末的比容,VS是相应固体材料在相同温度T和压力p下的比容。p-${\alpha} $模型具有不同的表达形式,本研究仅关注其中最常见的3种。

    • p-${\alpha} $模型常用于描述动态冲击加载下塑性多孔材料的冲击特性。当应力范围较大时,要分解为弹性和塑性两部分。此模型中忽略了剪切强度的存在,并且假设在等温等压下,多孔材料和固体材料的内能相等。在低压力段,弹性压实响应决定压实行为和从部分压实状态下的释放过程。通常初始的弹性状态${\alpha _{\rm E}} $很小,在计算中经常会忽略,或设为${\alpha _{\rm E}} ={\alpha _0} $直到达到阈值压力。

      塑性变形中的$\alpha$是对Butcher和Karnes方法[6]的二次改进。一般公式如下

      $ \alpha = 1 + \left( {{\alpha _{\rm E}} - 1} \right){\left( {\frac{{{p_{\rm S}} - p}}{{{p_{\rm S}} - {p_{\rm E}}}}} \right)^N} $

      式中:N是经验拟合参数,可以视实验具体情况改变。

    • Carroll和Holt[7]通过将多孔材料中孔洞的坍塌用一个由相同密度的空心球的压实相平衡来补充p-${\alpha} $模型。假设应力张量为球形,只考虑压力p。静态球形公式可对压实响应进行预测,因此压实响应由材料性质和孔洞的几何形状决定。假设固体材料不可压缩,将孔洞坍塌分为3个区域。在最初的弹塑性区,观察到在压力达到临界值前孔隙度基本没有发生变化。这使压实过程可以得到简化:在压实过程的前两个区域,压力小于pE时孔隙率的变化可以忽略不计。pE的大小使用下式计算

      $ {p_{\rm E}} = \frac{2}{3}Y\ln \frac{{{\alpha _0}}}{{{\alpha _0} - 1}} $

      式中:Y是固体材料的屈服强度,${\alpha _0} $是初始孔隙度。压力大于pE时为第3个阶段,材料塑性变形,压实由孔洞坍塌造成,此时的孔洞坍塌关系可以用下式简化

      $\alpha=\left(1 - {{\rm e}^{ - 3p/2Y}}\right)^{-1} $

      静态球面的p-${\alpha} $模型只考虑了与速率无关的屈服特性,因此可能不适用于高速率敏感粉体材料的压实响应,而且此模型忽略了颗粒形态可能造成的影响。

    • Brown等[19]p-${\alpha} $模型与前两种不同,他们用同一种机制描述弹塑性压实,直接描述了脆性微粒系统的动态压实,并没有将弹塑性机制分开研究。幂律公式的表达式如下

      $ \alpha = {\left( {P_{\rm S}}/p \right)^{\textstyle\frac{1}{n}}} $

      式中:PS是压碎强度,n是经验参数。

    • p-$\lambda $模型明确地将单个成分的物质属性包含在一个非均匀混合物中[20]。与p-${\alpha} $模型类似,该模型使用单内部状态变量λ描述压实过程并忽略粉末混合物的剪切强度,只考察了一般的p-V响应。该模型分为3个不同的区域:一端为不平衡的弹性响应,另一端为单个分量的压力平衡状态,两者之间为过渡区。假设压力较小时变形是弹性的。使用两种不同的附加模型定义混合物的体积弹性模量K和弹性响应:等压和等应变。等压响应是混合模量的下限,等应变响应是上限。在选择合适的混合模量后,通过应变的定义确定了非平衡弹性区混合物的p-V响应;在原方程的基础上,提出了一个简单的质量分数加性方程。在过渡区,假定材料以非平衡弹性状态和压力平衡状态的某种组合形式存在。材料在给定压力值下的相对量通过变量k来分配,k定义为完全压实状态下的物质质量分数。λ取值0~1。当材料完全处于不平衡的弹性状态下λ = 0;当材料处在压力平衡状态时,λ = 1。λ的具体函数关系为

      $\lambda = 1 - {{\rm e}^{ -{{\left( {p/Y} \right)}^{\scriptstyle{n}}}}} $

      式中:p是组件的局部压力差异;Y是固体的屈服强度,当混合物屈服时,会同时在混合物中引发屈服和压实;n无法从单个组件属性中进行定义。

    • 为了检验测量数据与模型拟合的精度,必须给出混合粉末的Hugoniot线,当压力超过粉体的致密压力后,测量数据点应该落在计算的粉末Hugoniot线上。本研究粉末的Hugoniot线使用等容Mie-Grüneisen方法计算得到。首先,基于W、Cu固体Hugoniot数据计算得到相同孔隙度的W、Cu粉体Hugoniot线,再经由Barry提出的两相混合Hugoniot原理计算得到W-Cu混合粉末的Hugoniot线。该原理主要存在两个基本假设:

      (1)各组分压力相等;

      (2)$V\left( p \right) = \mathop \sum \nolimits_{i = 1}^n {m_i}{V_i}\left( p \right)$${E_{\rm C}}\left( V \right) = \mathop \sum \nolimits_{i = 1}^n {m_i}{E_{{\rm C}i}}\left( V \right)$$\mathop \sum \nolimits_{i = 1}^n {m_i} = 1$EC是材料比内能。

      对于p-${\alpha} $模型,孔隙度$\alpha$可由Nn的值经拟合得到,pE是从弹性压缩转为塑性压缩的临界值,${\alpha _{\rm E}} $是相应的孔隙度,pS由实验数据和Mie-Grüneisen粉末状态方程的交点确定。p-$\lambda $模型中的曲线拟合与上述方法类似。混合物的屈服强度Y采用体积分数平均法计算。压力计算基于飞片和靶板材料的状态方程(EOS)。所用参数在表2中给出,其中C0S是线性冲击状态方程uS=C0+SuP参数 , γ0是Grüneisen常数。符合最佳拟合结果的数据在表3中给出。

      Material${\rho _0}$/(g·cm–3)C0/(km·s–1)S${\gamma _0}$V0/(cm3·g–1)Y/GPa
      W19.35 4.0641.2041.780.052 01.700
      Cu8.933.9101.5101.970.112 10.065
      W-Cu15.1174.0271.2771.830.066 41.030

      表 2  组分和混合物的相关冲击和材料特性

      Table 2.  Relevant shock and material properties of constituents and mixtures

      ModelParameters
      pS/GPapE/GPa${ \alpha _{\rm{E} } }$N
      p-${\alpha} $(MQ)4.600.057 71.3591.952
      p-${\alpha} $(PL)pS/GPan
      4.607.80
      p-${\alpha} $(SS)Y/GPa
      1.273 4
      p-$\lambda $n
      1.51

      表 3  模型拟合参数

      Table 3.  Model fitting parameters

      图6为多个模型的拟合结果。可以看到粉末在4~5 GPa压力下达到致密,压力继续升高后,数据点落到计算得到的粉体Hugoniot线上。除p-${\alpha} $ PL模型外,其他模型均较好地描述了W-Cu混合粉末的冲击压缩响应。但在孔隙度接近零时,预测结果均与计算的粉末Hugoniot线有一定偏差,这是由于Mie-Grüneisen方法忽略了材料强度和实验初始孔隙度测量误差所致。p-${\alpha} $ MQ和p-${\alpha} $ PL模型中都存在压溃强度pS,经拟合得到pS= 4.6 GPa;参数N影响曲线的曲率,N取值大则曲率增加,p-${\alpha} $ PL模型受N取值影响较大,拟合效果较差。p-${\alpha} $ SS模型仅受粉末的屈服强度Y的影响,完全忽略了粉末的微观尺度效应(颗粒形状和大小),但得到了较好的拟合结果,预测的压溃强度pS为4~5 GPa,与p-${\alpha} $ MQ和p-${\alpha} $ PL模型一致。p-$\lambda $模型也可以较好地描述低压段的冲击压实响应,但孔隙度无限接近于零时,在屈服强度Y一定的前提下,参数n的取值决定了曲线的曲率,当n减小时,曲率减小,pS急剧增大,而Y的取值影响pE的值。由于低于4 GPa的数据点仅有一个,所以很难确定其余模型哪个更适合描述W-Cu混合粉末的冲击压缩响应。拟合过程中发现,所有模型通过参数调整都可以获得较好的拟合效果,是典型的描述型模型,预测性较差。如p-${\alpha} $ MQ和p-${\alpha} $ PL模型中的参数pS无法通过实验测量,p-${\alpha} $ SS和p-$\lambda $模型中的Y受混合颗粒特性(如颗粒的尺寸、颗粒级配和颗粒形状等)的影响显著,无法准确预测。此外,除p-${\alpha} $ PL模型外,其余模型均能很好地描述粉末压缩中由弹性转塑性的行为,而p-${\alpha} $ PL模型经参数调整后也可以描述这一转变,但其曲线曲率变小,如图6所示,与实验数据差别较大。

      图  6  4种模型拟合结果

      Figure 6.  Fitting results of four models

    • 测量了W-Cu混合粉末的Hugoniot数据,系统分析了测量误差,发现粉末致密度测量精度对结果影响较大。利用Mie-Grüneisen方法和Barry等压混合法,可以利用单质W、Cu粉末的Hugoniot关系很好地预测同孔隙度W-Cu混合粉末的Hugoniot线,但由于忽略了粉体强度和粉体初始密度测量误差较大的原因,预测结果在低压段与实验偏差较大。分别应用3种p-${\alpha} $模型和p-$\lambda $模型对实验结果进行拟合,发现除p-${\alpha} $ PL模型外,其他模型均能较好地描述W-Cu混合粉末的冲击压缩响应,表明这些原本用于描述单质粉末的冲击压实响应的模型对于更为复杂的混合粉末也是适用的,但所有模型的拟合精度均受经验参数选择的影响,预测能力较差。

参考文献 (20)

目录

    /

    返回文章
    返回