Loading [MathJax]/jax/element/mml/optable/BasicLatin.js

应变率相关的橡胶本构模型研究

魏家威 石霄鹏 冯振宇

白晨, 杨昆, 吴艳青, 高洪泉, 薛海蛟. 不同类型装药侵彻安全性数值模拟[J]. 高压物理学报, 2021, 35(6): 065101. doi: 10.11858/gywlxb.20210754
引用本文: 魏家威, 石霄鹏, 冯振宇. 应变率相关的橡胶本构模型研究[J]. 高压物理学报, 2022, 36(2): 024205. doi: 10.11858/gywlxb.20210815
BAI Chen, YANG Kun, WU Yanqing, GAO Hongquan, XUE Haijiao. Numerical Simulation for PBX Charges Safety of Different Types During Penetration[J]. Chinese Journal of High Pressure Physics, 2021, 35(6): 065101. doi: 10.11858/gywlxb.20210754
Citation: WEI Jiawei, SHI Xiaopeng, FENG Zhenyu. Strain Rate Dependent Constitutive Model of Rubber[J]. Chinese Journal of High Pressure Physics, 2022, 36(2): 024205. doi: 10.11858/gywlxb.20210815

应变率相关的橡胶本构模型研究

doi: 10.11858/gywlxb.20210815
详细信息
    作者简介:

    魏家威(1996-),女,硕士研究生,主要从事材料动态力学行为研究. E-mail:13821971667@163.com

    通讯作者:

    冯振宇(1966-),男,博士,教授,主要从事飞机结构强度研究. E-mail:mhfzy@163.com

  • 中图分类号: O345

Strain Rate Dependent Constitutive Model of Rubber

  • 摘要: 为研究橡胶在不同应变率下的响应特性,建立应变率相关的橡胶黏超弹性本构模型,分别采用超弹性本构模型和黏弹性本构模型表征其非线性弹性行为和应变率相关的弹性行为。首先,对于超弹性模型,基于最小二乘法,对比了Mooney-Rivlin模型、修正的Mooney-Rivlin模型、Yeoh模型、修正的Yeoh模型、Ogden模型和Arruda-Boyce模型等超弹性本构模型的拟合能力。结果表明,经修正的Mooney-Rivlin模型和Yeoh模型的拟合优度与Ogden模型和Arruda-Boyce模型接近。在此基础上,基于一种参数较少且拟合效果良好的修正Mooney-Rivlin模型和应变率相关的Maxwell模型,建立了橡胶黏超弹性本构模型,考察了该黏超弹性本构模型在单轴拉伸和单轴压缩情况下中高应变率时的拟合能力。结果表明,对于这两种受力情况下的应变率相关的实验数据,该黏超弹性本构模型的拟合优度均在0.95以上。研究结果为大应变率范围内单轴拉伸和单轴压缩下橡胶的本构模型选择提供了参考。

     

  • 侵彻战斗部主装药广泛采用PBX炸药,由于服役环境可能同时存在多种外界刺激,炸药自身的力学响应会影响侵彻过程中装药的安全性和战斗部的毁伤能力[1],因此PBX装药对弹药的整体性能有着决定性意义。而在目标侵彻过程中,侵彻战斗部的内部装药面临复杂的力学环境[2],同时复杂的装药结构也会对内部装药产生影响。在这些外部刺激下,炸药中的局部损伤处可能产生热点,导致意外点火的发生,进而可能转化成为更剧烈的爆燃或者爆轰,引发多种事故。因此研究各种复杂侵彻环境下装药的力学性能与损伤响应,对保障弹药在整个服役过程中的安全性与可靠性具有重要意义。

    对于侵彻过程中内部装药的安全性而言,建立能够准确描述炸药力学响应和损伤行为的本构模型,是通过数值模拟方法研究侵彻装药安全性问题的关键。PBX炸药动态损伤本构模型可以划分为两类:第一类是基于黏弹性、黏塑性等连续介质损伤力学基础理论,考虑应变率、温度对PBX炸药力学行为的影响而建立的宏观唯象经验本构模型[3-9];第二类则考虑微缺陷是PBX炸药中一类重要的细观结构,是基于材料力学损伤-点火行为和细观缺陷演化与热点形成之间的联系而建立的PBX炸药细观损伤力学本构模型[10-13]。近年来,国内外学者基于以上模型对侵彻装药安全性进行了广泛的研究。张馨予等[14]将孔隙压塌损伤、炸药晶体破碎损伤、黏结剂脱黏等多种细观损伤形式耦合到炸药宏观本构模型中,研究了侵彻环境下弹体装药的损伤分布情况。石啸海等[15-16]基于内聚裂纹模型,模拟了战斗部侵彻半无限大混凝土过程中PBX装药的动态力学响应及损伤演化,并且对某种缩比弹侵彻混凝土靶板进行了数值模拟,从过载、裂纹宽度、裂纹含量等角度比较了弹头形状对装药损伤的影响,结果表明控制弹头曲径比有利于减小装药损伤。成丽蓉等[17]基于裂纹摩擦、孔洞塌缩两种热点生成机制细观模型,开展了侵彻单层和多层典型靶板时战斗部装药的动态响应、损伤演化及热点生成对比研究。Li等[18]基于微裂纹动态损伤模型,通过数值模拟方法预测了实弹侵彻实验中PBX1314的损伤分布、点火位置与点火区域的反应进程,结果表明,侵彻过程中经多次撞击载荷作用,PBX1314尾端面局部区域达到点火临界条件而发生点火,反复撞击过程中PBX1314内部微裂纹摩擦引起的能量局部化是引发点火的重要原因。

    上述模型往往仅考虑了剪切裂纹热点等单一缺陷机制的细观力热响应过程,在研究复杂载荷条件下不同类型PBX炸药损伤-热点主导机制的自适应能力方面,模型还需要进一步改进。此外,上述研究工作多围绕压装型炸药PBX展开,针对浇注型PBX炸药的研究较少,由于浇注类PBX炸药和压装类PBX炸药的材料组分、物理状态和成型工艺不同,力学性能与点火特性也存在差异,因此对比其在侵彻环境下的响应特性,对于战斗部装药材料的选取和装药结构设计具有重要意义。基于此,本研究应用前期发展的PBX炸药微裂纹-微孔洞力热化学耦合细观模型[18],考虑微裂纹-微孔洞两种细观缺陷演化对炸药损伤-热点形成的影响,通过动态分离式霍普金森压杆(SHPB)实验校核模型参数,并利用数值模拟方法分析两类典型装药(压装PBX04和浇注GOFL-5)在弹体侵彻混凝土薄板过程中的应力波传播、损伤演化和温升响应情况,为深入理解侵彻过程中装药的力学-损伤-点火响应提供参考。

    PBX炸药微裂纹-微孔洞力热化学耦合细观模型(CMM)中考虑了拉伸张开、剪切张开、纯剪切、剪切摩擦、摩擦自锁5种微裂纹演化模式,以及微孔洞坍塌与扭曲变形两种演化模式,如图1所示。同时,CMM模型还囊括了剪切裂纹热点与孔洞坍塌热点子模型,具有复杂应力状态-微缺陷演化模式自判断能力,以及两种微缺陷热点自启动能力。CMM模型中总体应力、应变分别分解为偏量部分与体量部分,依次建立微裂纹相关偏量本构关系与微孔洞相关体量本构关系,二者通过Gurson屈服准则进行耦合,通过微裂纹、微孔洞演化方程更新相关变量,建立剪切裂纹摩擦热点与孔洞坍塌热点子模型,模型细节与算法实现详见文献[1920]。

    图  1  CMM模型所考虑的微缺陷演化机制与变形机制示意图
    Figure  1.  Conceptual diagram of all kinds of considered mechanisms in the current model

    将总体偏应变(ε)分解为黏弹性应变(εve)、塑性应变(εp)以及由于微裂纹扩展所引起的微裂纹应变(εcr) 3部分,分别描述PBX炸药内黏结剂基体材料所表现出的黏弹性变形、塑性变形以及由于微裂纹扩展对材料力学行为的影响。

    PBX炸药黏弹性变形由广义Maxwell模型描述,微裂纹张开/剪切扩展引起的裂纹应变由SCRAM模型描述,材料的偏量本构关系可表示为

    ˙S=2GA0(˙ε˙εp)B0(S+C0) (1)

    式中:S为应力偏量,G为剪切模量,A0B0C0表示为

    A0=11+αe3(ˉc/a)3,B0=αe(ˉc/a)2˙ˉc/a1+αe3(ˉc/a)3,C0=Nn=1Sn/τnαe(ˉc/a)2˙ˉc/a (2)

    式中:ˉc为微裂纹平均尺寸;a为微裂纹初始特征尺寸,a3=6GβSnτn为广义黏弹性体元中第n个Maxwell黏弹性单元中的偏应力和松弛时间。

    αe={3p05νp<0,β=64π(1ν)N015(2ν)G (3)

    式中:p为压力,N0为初始裂纹密度,ν为泊松比。

    为描述材料内微孔洞演化所发生的不可逆损伤对PBX炸药体积变形的影响,采用孔隙率相关状态方程

    p(ρ,e,f)=(1f)[ρsc20ηs(1sηs)2(1Γsηs2)+Γsρses] (4)

    式中:f为孔隙率;ρ、e分别为孔隙材料的密度与比内能,密实材料对应的ρs = ρ/(1 − f),es = eηs = 1 − ρs0/ρsρs0为密实材料的初始密度;Γs为Grüneisen系数;c0s为材料参数。

    含孔隙PBX材料的塑性变形采用经典Gurson模型描述,模型中材料屈服面与von-Mises等效应力σe和压力p相关

    F(σe,p,f)=(σeYM)2+2fcosh(3p2YM)f21=0 (5)

    考虑动态加载下材料硬化效应与应变率效应,密实材料(f = 0)的屈服强度可表示为

    YM=[σ0+h(ˉεpM)n][1+Cln(1+˙ε)] (6)

    式中:σ0˙ε0=10 s−1加载下所对应的初始屈服应力,ˉεpM为等效塑性应变,˙ε=˙ε/˙ε0为等效应变率,h为硬化模量,C为应变率相关系数。

    基于Griffith能量释放率裂纹扩展准则,微裂纹扩展方程可表示为

    ˙ˉc=˙cmax (7)

    式中:{\dot c_{\max }}为最大扩展速度,\bar \gamma 为材料比表面能,{g_{\rm{dom}}}\left( {{{\sigma}},\bar c} \right)为主裂纹对应的能量释放率。

    由于材料内微裂纹方向分布具有随机性,存在临界微裂纹方向,其对应的能量释放率最大,即该方向的微裂纹在最小施加应力下最先发生失稳扩展,因此定义为主裂纹。主裂纹方向的确定与当前应力状态相关(\sigma 1\sigma 3平面),微裂纹扩展对应的材料整体损伤度定义为{d_{\rm{cr}}} = {{{{\bar c}^3}} / {\left( {{a^3} + {{\bar c}^3}} \right)}}

    考虑孔洞坍塌引起的孔隙率减小,以及孔洞扭曲引起的孔隙率增加两种变形机制,孔隙率演化方程可表示为

    \dot f = \left( {1 - f} \right)\dot \varepsilon _V^{\rm{p}} + f{k_{\rm{w}}}\omega \left( \sigma \right)\frac{{{s_{ij}}\dot e_{ij}^{\rm{p}}}}{{{\sigma _{\rm{e}}}}} (8)

    式中:\dot \varepsilon _V^{\rm{p}} = {{\dot \varepsilon _{\rm{kk}}^{\rm{p}}} / 3},为塑性体积应变率;kw为剪切相关材料参数;\omega 与应力状态相关,取值范围为0 ≤ \omega ≤ 1。

    \omega \left( {{\sigma}} \right) = 1 - {\left( {\frac{{27{{{J}}_3}}}{{2{{\sigma}} _{\rm{e}}^3}}} \right)^2} (9)

    式中:J3为表示应力张量的第三不变量。

    采用一维热传导方程来描述剪切裂纹表面及其周围区域摩擦生热、熔化、点火与传热等热力学过程

    \rho {c_V}{\dot T_{\rm{hs}}} = \frac{\partial }{{\partial x}}\left( {\kappa \frac{{\partial {T_{\rm{hs}}}}}{{\partial x}}} \right) + \rho {Q_{\rm{r}}}Z{{\rm{e}}^{ - E/(R{T_{\rm{hs}}})}} + \varphi {\mu _v}{\dot \varepsilon_{\rm{m}} ^2} (10)

    式中:x为沿微裂纹法向的坐标轴;等式右端3项分别表示热传导项、化学反应释放热量以及熔化区域(ThsTm)液相黏性流动生热;ThscVκQr分别为微裂纹热点温度、比定容热容、热传导系数以及单位质量化学反应放热;ZE为Arrhenius反应速率方程参数;\varphi (0 ≤ \varphi ≤ 1)为熔化百分数;{\;\mu _v}为熔化区域黏性系数;\dot \varepsilon_{\rm{m}}= \varphi {v_{\rm{c}}}/l为熔化区域的剪切应变率,vc为微裂纹滑动速度,熔化区域宽度l由熔化区域的传播速度确定,\dot l(t) = \Delta \left( {k{T_x}} \right)/\rho L,其中\Delta \left( {k{T_x}} \right)表示固液分界面热流量梯度。

    微孔洞周围材料温度(Tvo)分布情况由一维球形热传导方程进行描述

    \rho {c_V}{\dot T_{\rm{vo}}} = k\left[ {\frac{{{\partial ^2}{T_{\rm{vo}}}}}{{\partial {r^2}}} + \frac{2}{r}\frac{{\partial {T_{\rm{vo}}}}}{{\partial r}}} \right] + \dot w_{\rm{vp}}^* + \rho {Q_{\rm{r}}}Z{{\rm{e}}^{ - E/(R{T_{\rm{vo}}})}} (11)

    其中,孔洞周围单位体黏塑性功生成速率可表示为

    \dot w_{\rm{vp}}^*\left( r \right) = s_{ij}^*\dot \varepsilon _{ij}^* = \frac{{ - 2{Y_{\rm{M}}}\dot \varepsilon _V^{\rm{p}}}}{{{r^3}}}\frac{{1 - {f_0}}}{{1 - {f_{\rm{vc}}}}}b_0^3 + \frac{{4\eta {{\left( {\dot \varepsilon _V^{\rm{p}}} \right)}^2}}}{{{r^6}}}{\left( {\frac{{1 - {f_0}}}{{1 - {f_{\rm{vc}}}}}} \right)^2}b_0^6 (12)

    式中:等号右边两项分别表示基体材料的塑性效应与黏性效应所产生的功率,b0为微孔洞初始外径,\eta 为微孔洞周围基体材料相关的黏性系数。

    应用SHPB方法对中国工程物理研究院化工材料研究所提供的压装和浇注两类PBX炸药进行动态力学性能测试,浇注类复合炸药选取GOFL-5炸药(HMX、FOX-7、黏结剂的质量分数分别为50%、35%、15%),样品尺寸为20 mm × 20 mm,密度为1.75 g/cm3,初始孔隙率为1.01%。压装类复合炸药选取PBX04炸药(HMX、黏结剂的质量分数分别为95%、5%),密度为1.82 g/cm3,初始孔隙率为1.08%。在中国工程物理研究院化工材料研究所物理与力学性能实验室开展SHPB实验,实验中采用直径为20 mm的铝杆,子弹长度为300 mm,入射杆、透射杆长度分别为2 000、1 500 mm。Yang等[19]在前期研究了CMM模型中微缺陷参数({\bar c_0}, N0, \bar \gamma , f0, kw)对材料力学性能的影响规律。基于此,依次对CMM模型中的广义Maxwell黏弹性参数(G1G5, τ1τ5)、塑性参数({\sigma _0}, h, C, n)、微裂纹模型参数({\bar c_0}, N0, \bar \gamma, m)、微孔洞模型参数(f0, kw)以及Mie-Grüneisen状态方程参数(c0, s, Г)进行优化调整,使得标定的参数能够最大程度地符合计算与实验应力-应变曲线。压装炸药PBX04与浇注炸药GOFL-5标定曲线如图2所示,相关材料参数如表1所示,其中 {\mu _{\rm{s}}} 为静摩擦系数。由表1可知:两种炸药材料的黏弹性模量G1G5存在数量级差异,GOFL-5炸药的塑性参数({\sigma _0}, h, C)均远小于PBX-5炸药,GOFL-5炸药的初始微裂纹密度与微裂纹尺寸均小于PBX-5炸药。

    图  2  计算与实验得到的两类炸药在不同应变率下的应力-应变曲线
    Figure  2.  Calculated and experimental stress-strain curves of two kinds of explosives at different strain rates
    表  1  GOFL-5与PBX04材料参数
    Table  1.  Material parameters for GOFL-5 and PBX04
    Material\;\rho 0/(kg·m−3)G/GPa\nuG1/MPaG2/MPaG3/MPaG4/MPaG5/MPa{\tau{_1^{-1} } } /{\rm{s} }{^{-1} }
    GOFL-51 7500.550.3 167 30.45 90.03 185.6 120.00
    PBX041 8208.250.31 9401 175.001 521.001 909.01 688.00
    Material{\tau{_2^{-1} } } /{\rm{s} }{^{-1} }{\tau{_3^{-1} } } /{\rm{s} }{^{-1} }{\tau{_4^{-1} } } /{\rm{s} }{^{-1} }{\tau{_5^{-1} } } /{\rm{s} }{^{-1} }\sigma{{_0}}/MPaCh/MPan{\bar c}0/μm
    GOFL-57.32 × 1037.32 × 1047.32 × 1057.32 × 106 2.20.76 4.50.4530
    PBX049.00 × 1039.00 × 1049.00 × 1052.00 × 10640.00.101500.01.0030
    MaterialN0/cm−3\bar \gamma /(J·m−2){\dot c_{\max }}/(m·s–1)\;\mu_{\rm{s}}mf0/(m·s–1)kwc0/(m·s–1)s \varGamma
    GOFL-5 30.53000.35.00.012.01 0000.46 0.89
    PBX043001.43000.55.00.012.02 5002.26 1.50
    下载: 导出CSV 
    | 显示表格

    对比压装炸药与浇注炸药的应力-应变曲线,可以看到PBX04炸药峰值应力远高于GOFL-5炸药,但其破坏应变远小于GOFL-5炸药。两种材料的力学行为差异主要与其内部微缺陷数量和黏结剂含量相关。压装类PBX04炸药内微裂纹含量较多,且黏结剂含量较少,因此材料表现为准脆性材料破坏特征。浇注类GOFL-5炸药内部微缺陷数量较少,且对材料起到增韧效应的黏结剂含量较多,表现为韧性材料破坏特征。

    建立二维侵彻混凝土靶板的计算模型,设置轴对称边界条件和加载条件,建立1/2模型进行计算,选取侵彻速度为800 m/s,弹体尺寸如图3所示,侵彻的混凝土靶板厚度为0.5 m。计算过程中,由于侵彻过程中靶体内温度升高的区域是有限的,靶体外边界的温度不会受到侵彻过程中靶体温度变化区域的影响,因此在靶体边界处设置温度边界条件,边界温度设置成294 K。

    图  3  有限元计算模型
    Figure  3.  Finite element calculation model

    弹体材料本构模型使用Johnson-Cook模型。该模型是一个经验型的黏塑性本构模型,由Johnson和Cook在1983年首先提出,多用于描述金属材料在高载荷、高应变率和高温下的应力-应变关系,模型形式简单,本构参数少,能较好地描述材料的加工硬化效应、应变率效应和温度软化效应。在侵彻速度不高时弹体变形较小,因此本研究运用Johnson-Cook模型来描述壳体材料4340钢的力学响应,其主要材料参数见表2,其中:\rho为密度,c为比热容,α为热膨胀系数,Tm为熔化温度,\kappa 为导热系数,M为温度软化系数,N为应变强化参数,C1C4D1D5为材料常数。混凝土靶板的材料模型使用HJC本构模型[21]。该模型是针对混凝土材料提出的一种率相关损伤型本构模型,适用于大应变、高应变率和高压情况,其主要材料参数见表3,其中:ABN为材料硬化参数,FC为材料准静态单轴抗压强度,K1K3为材料体积模量参数,pCpLUCUL为压力参数,D1D2为损伤参数,\varepsilon_{\rm{min}}为最小断裂应变。

    表  2  弹体材料参数
    Table  2.  Parameters of projectile material
    Physical properties Johnson-Cook model
    \;\rho /(g·cm−3)c/(J·kg−1·K−1)\kappa /(kW·m−1·K−1)\alpha /(m·K−1)Tm/℃G/GPa NM
    7.8247838.113.24 × 10−51 793.15 774.970.26 1.03
    下载: 导出CSV 
    | 显示表格
    Johnson-Cook model Damage model
    C1/MPaC2/MPa C3 C4 D1 D2 D3 D4 D5
    792.21509.521.4 0 –0.8 2.1 –0.5 20.0 0.61
    下载: 导出CSV 
    | 显示表格
    表  3  混凝土靶板材料参数
    Table  3.  Material parameters of concrete plate target
    Physical properties HJC model
    \;\rho /(g·cm−3)c/(J·kg−1·K−1)κ/(kW·m−1·K−1)α/(m·K−1)G/GPaFC/MPaABNC
    2.286541.764.32 × 10−5 16.4040.680.751.650.767.0 × 10−3
    下载: 导出CSV 
    | 显示表格
    HJC model Damage model
    pC/MPaUC pL/GPa UL K1/GPa K2/GPa K3/GPa D1 D2 \varepsilon min
    13.565.80 × 10−4 1.05 0.10 17.40 38.80 29.80 0.03 1.0 0.01
    下载: 导出CSV 
    | 显示表格

    为分析药柱内压缩波传播情况,初始时刻沿药柱中心线由装药头部至尾部等间隔依次选取1#~5#位置,如图4(b)所示,图4(a)图5(a)分别为压装(PBX04)与浇注炸药(GOFL-5)内1#~5#位置处的压力变化曲线。由图4(a)可知,弹体与混凝土靶板初始接触时,弹体受到强烈的压缩波作用,装药头部压力最大,入射压缩波S1分别于104和179 μs到达装药中部(3#位置)和装药尾部(5#位置),当应力波传播至尾部自由端面时反射拉伸波,两种应力波相互抵消,装药尾部压力值最低。930 μs时,装药尾部在弹性恢复与拉伸波的作用下与弹体内壁面发生撞击,尾部压力迅速升高。由图5(a)可知,入射压缩波S1到达装药中部和装药尾部的时刻分别为123和213 μs,对比图4(a)图5(a)可知,相同侵彻速度下,浇注药内压缩波S1的波速(约2.11 km/s)小于压装药(约2.57 km/s)。

    图  4  压装药PBX04的压力演化情况
    Figure  4.  Pressure evolution of PBX04
    图  5  浇注药GOFL-5的压力演化情况
    Figure  5.  Pressure evolution of GOFL-5

    不同时刻压装、浇注药内压力云图演化如图4(b)图5(b)所示。随着压缩波的传播,药柱整体受压区域逐渐扩大,装药尾部与弹体内壁产生相对间隙,浇注GOFL-5炸药由于材料强度较低,抵抗变形的能力较弱,侵彻过程中装药尾部的最大间隙(dmax)较压装药更大(GOFL-5,dmax = 1.61 cm;PBX04,dmax = 0.68 cm)。结合图4(a)图5(a)可知,浇注GOFL-5炸药在侵彻过程中产生了较大的变形,装药尾部与弹体内壁发生撞击时尾部压力峰值(pmax)也更高(GOFL-5,pmax = 234 MPa;PBX04,pmax = 128 MPa)。

    图6为压装、浇注药柱内微裂纹损伤(dcr)演化情况。由图6(a)可知,加载初期压装药微裂纹损伤主要集中在头部,随着压缩波的传播,微裂纹损伤逐渐扩展至装药中部,装药尾部与弹体内壁发生撞击后尾部微裂纹损伤较严重。由图6(b)可知,浇注药的流动性较好,在侵彻过程中产生了较大的变形,尾部装药受到向药柱中心的拉应力作用产生裂纹扩展损伤。整个侵彻过程中,两类炸药微裂纹损伤较严重区域均为装药头部和尾部,由于开始侵彻时压装药内入射压缩波强度高于浇注药,同时压装药初始微裂纹密度较高,加载初期压装药头部微裂纹损伤高于浇注药。

    图  6  不同时刻两类炸药微裂纹损伤演化情况
    Figure  6.  Evolution of microcrack damage of two kinds of explosives

    图7为不同时刻压装、浇注药柱内不同位置处微孔洞损伤(dvo = f0ft)随时间演化曲线。对比图6图7可知,微孔洞坍塌时间尺度相对于微裂纹扩展时间尺度明显更小,微孔洞随着入射压缩波传播而发生坍塌,孔洞坍塌损伤发生的时刻与入射压缩波到达的时刻基本一致,而由于裂纹扩展损伤主要受拉伸和剪切状态影响,入射压缩波刚到达装药对应位置时并未发生明显的微裂纹扩展损伤。

    图  7  不同时刻两类炸药的微孔洞损伤演化
    Figure  7.  Evolution of microvoid damage of two kinds of explosives

    PBX04和GOFL-5内不同位置微裂纹与微孔洞热点温度随时间演化曲线分别如图8图9所示。由图8(a)可知,压装药微裂纹热点温升主要集中在装药头部和装药前端,1#和2#处的微裂纹热点峰值分别为362和487 K。根据CMM模型中由应力二轴度可确定5种微裂纹状态(拉伸张开、剪切张开、纯剪切、剪切摩擦、摩擦自锁),5种微裂纹演化模式如图10所示,由此可进一步分析温升较高位置处的微裂纹状态。整个侵彻过程中,1#和2#处5种微裂纹状态的频率分布如图11所示。侵彻过程中弹体内装药处于压力主导的应力状态,微裂纹大多处于摩擦自锁状态而非剪切裂纹扩展状态,同时热点密度较小可能引起热点湮灭现象,因此微裂纹热点温度达到峰值后逐渐降低。对比图8(a)图9(a)可知,由于GOFL-5材料的初始微裂纹密度、微裂纹尺寸以及微裂纹扩展速率较低,因此GOFL-5材料中剪切裂纹引起的热点温升较PBX04更低。结合图8(b)图9(b),随着入射压缩波的传播孔洞发生坍塌,与GOFL-5不同的是,PBX04装药尾部微孔洞温升峰值主要来自弹体撞击造成的孔洞坍塌,而坍塌孔洞周围黏塑性功引起的温升较小(约20 K),不足以引起点火。对比图8图9可知:对于压装药而言,裂纹摩擦相比孔洞坍塌引起的温升更高,剪切裂纹热点为压装药主导的温升机制;对于浇注药而言,两种热点机制引起的温升差别不大。

    图  8  不同时刻PBX04中微裂纹与微孔洞相关热点温度曲线
    Figure  8.  Microcrack and microvoid related hotspot temperature evolution curves for PBX04
    图  9  不同时刻GOFL-5中微裂纹与微孔洞相关热点温度曲线
    Figure  9.  Microcrack and microvoid related hotspot temperature evolution curves for GOFL-5
    图  10  微裂纹不同扩展形态示意图
    Figure  10.  Schematic diagram of different propagation patterns of microcracks
    图  11  PBX04装药中不同位置的微裂纹状态频率分布
    Figure  11.  Frequency distribution of microcracks state in different positions of PBX04 charge

    应用PBX炸药微裂纹-微孔洞力热化学耦合细观模型,研究了压装PBX04和浇注GOFL-5两类典型装药在弹体侵彻混凝土薄板过程中应力波的传播、损伤演化和温升响应情况,对比分析了相同侵彻条件下两类炸药力学-损伤-温升响应的差异性,得到以下主要结论。

    (1) 根据实验曲线标定了PBX04和GOFL-5微裂纹-微孔洞本构模型参数,两种炸药材料的弹性模量存在数量级差异,GOFL-5炸药的屈服强度、硬化模量、初始微裂纹密度和微裂纹尺寸均小于PBX04炸药。

    (2) 加载初期压装药头部微裂纹损伤高于浇注药,而浇注药的流动性较好,侵彻过程中产生了较大的变形,当装药尾部和壳体内表面发生撞击时形成高压区,整个侵彻过程中两类炸药微裂纹损伤较严重的区域均为装药头部和尾部,装药损伤较严重的区域往往容易引起能量聚集,进而在这些局部高温区形成热点,在侵彻弹体设计时应作为重点防护区域。

    (3) 通过计算可知,800 m/s速度侵彻混凝土薄板条件下,两种装药材料均未发生点火,裂纹摩擦热点为压装药主导的温升机制,而对于浇注药而言,两种热点机制引起的温升差别不大,且浇注药GOFL-5在侵彻过程中的温升较压装药PBX04更低。

  • 图  基于M-R模型拟合的工程应力-伸长比的结果

    Figure  1.  Fitting results of the principle stress versus the principle stretch using M-R model

    图  M-R模型和修正后的M-R模型拟合结果对比

    Figure  2.  Comparison of fitting results between M-R model and the modified M-R model

    图  基于Yeoh模型拟合的工程应力-伸长比的结果

    Figure  3.  Fitting results of the principle stress versus the principle stretch using Yeoh model

    图  Yeoh模型和修正后的Yeoh模型拟合结果对比

    Figure  4.  Comparison of fitting results between Yeoh model and the modified Yeoh model

    图  不同模型拟合Treloar3种实验数据的结果

    Figure  5.  Fitting results of different models of Treloar’s three kinds of experimental data

    图  不同实验类型下不同模型的拟合结果对比

    Figure  6.  Comparison of the fitting results of different models for different experiment types

    图  黏超弹性模型示意图

    Figure  7.  Schematic diagram of visco-hyperelastic model

    图  M-R黏超弹性本构模型拟合硅橡胶单轴拉伸实验数据[29]的结果

    Figure  8.  Fitting results of M-R visco-hyperelastic constitutive model for the uni-axial tensile experimental data of silicone rubber[29]

    图  M-R黏超弹性本构模型拟合硫化橡胶单轴压缩实验数据[28]的结果

    Figure  9.  Fitting results of M-R visco-hyperelasticconstitutive model for the uni-axial compressionexperimental data of vulcanized rubber[28]

    表  1  不同超弹性模型对 ST、PT 和 ET 实验数据的拟合效果比较

    Table  1.   Comparison of fitting results of different hyperelastic models on ST, PT and ET experimental data

    ModelEquationParametersR2
    M-R modelEq.(2)C10, C010.8043
    Modified M-R modelEq.(15)C10, C01, C200.9704
    Yeoh modelEq.(6)C10, C20, C300.9897
    Modified Yeoh modelEq.(21)C10, C20, C30, C010.9961
    Ogden model (N=2)Eq.(8) {\,\mu } 1, {\alpha } 1, {\,\mu } 2, {\alpha } 20.9769
    Ogden model (N=3)Eq.(8){\,\mu }1, {\alpha } 1, {\,\mu } 2, {\alpha } 2, {\,\mu }3, {\alpha } 30.9924
    A-B modelEq.(10)\,\mu, {\lambda } m0.9891
    下载: 导出CSV

    表  2  单轴拉伸和单轴压缩实验的拟合参数值

    Table  2.   Fitting parameter values of the uni-axial tensile and the uni-axial compression experiment

    ExperimentC10C01C20{E}1 {\theta } 0 \,\beta R2
    Uni-axial tensile experiment[29]0.731.001.0×10−147.9821.900.820.9811
    Uni-axial compression experiment[28]−0.86−0.100.30−4.409.300.700.9585
    下载: 导出CSV
  • [1] 胡小玲, 刘秀, 李明, 等. 炭黑填充橡胶超弹性本构模型的选取策略 [J]. 工程力学, 2014, 31(5): 34–42, 48.

    HU X L, LIU X, LI M, et al. Selection strategies of hyperelastic constitutive models for carbon black filled rubber [J]. Engineering Mechanics, 2014, 31(5): 34–42, 48.
    [2] 龚科家, 危银涛, 叶进雄. 填充橡胶超弹性本构参数试验与应用 [J]. 工程力学, 2009, 26(6): 193–198.

    GONG K J, WEI Y T, YE J X. Constitutive parametric experiment of tire rubber hyperelastic laws with application [J]. Engineering Mechanics, 2009, 26(6): 193–198.
    [3] ELIAS H G. Macromolecules, volume 1: structure and properties [M]. Boston, MA: Springer, 2012.
    [4] MOONEY M. A theory of large elastic deformation [J]. Journal of Applied Physics, 1940, 11(9): 582–592. doi: 10.1063/1.1712836
    [5] RIVLIN R S, SAUNDERS D W. Large elastic deformations of isotropic materials. Ⅶ. Experiments on the deformation of rubber [J]. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 1951, 243(865): 251–288.
    [6] YEOH O H. Some forms of the strain energy function for rubber [J]. Rubber Chemistry and Technology, 1993, 66(5): 754–771. doi: 10.5254/1.3538343
    [7] OGDEN R W. Non-linear elastic deformations [M]. New York: Dover Publications, 1997.
    [8] OGDEN R W. Large deformation isotropic elasticity: on the correlation of theory and experiment for incompressible rubberlike solids [J]. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 1972, 326(1567): 565–584.
    [9] ARRUDA E M, BOYCE M C. A three-dimensional constitutive model for the large stretch behavior of rubber elastic materials [J]. Journal of the Mechanics and Physics of Solids, 1993, 41(2): 389–412. doi: 10.1016/0022-5096(93)90013-6
    [10] TRELOAR L R G. Stress-strain data for vulcanized rubber under various types of deformation [J]. Rubber Chemistry and Technology, 1944, 17(4): 813–825. doi: 10.5254/1.3546701
    [11] BOYCE M C, ARRUDA E M. Constitutive models of rubber elasticity: a review [J]. Rubber Chemistry and Technology, 2000, 73(3): 504–523. doi: 10.5254/1.3547602
    [12] STEINMANN P, HOSSAIN M, POSSART G. Hyperelastic models for rubber-like materials: consistent tangent operators and suitability for Treloar’s data [J]. Archive of Applied Mechanics, 2012, 82(9): 1183–1217. doi: 10.1007/s00419-012-0610-z
    [13] MARCKMANN G, VERRON E. Comparison of hyperelastic models for rubber-like materials [J]. Rubber Chemistry and Technology, 2006, 79(5): 835–858. doi: 10.5254/1.3547969
    [14] 肖锐, 向玉海, 钟旦明, 等. 考虑缠结效应的超弹性本构模型 [J]. 力学学报, 2021, 53(4): 1028–1037. doi: 10.6052/0459-1879-21-008

    XIAO R, XIANG Y H, ZHONG D M, et al. Hyperelastic model with entanglement effect [J]. Chinese Journal of Theoretical and Applied Mechanics, 2021, 53(4): 1028–1037. doi: 10.6052/0459-1879-21-008
    [15] 施成, 周恒为, 丁明明, 等. 一种基于分子链统计理论的橡胶超弹性混合本构模型 [J]. 应用化学, 2021, 38(2): 228–235.

    SHI C, ZHOU H W, DING M M, et al. A hyperelastic mixed constitutive model for rubber based on molecular chain statistical theory [J]. Chinese Journal of Applied Chemistry, 2021, 38(2): 228–235.
    [16] 付宾, 杨晓翔, 李庆. 炭黑填充橡胶材料改进Mooney模型 [J]. 固体力学学报, 2017, 38(5): 408–415.

    FU B, YANG X X, LI Q. A revised Mooney model of carbon black filled-rubber materials [J]. Chinese Journal of Solid Mechanics, 2017, 38(5): 408–415.
    [17] 李雪冰, 危银涛. 一种改进的Yeoh超弹性材料本构模型 [J]. 工程力学, 2016, 33(12): 38–43. doi: 10.6052/j.issn.1000-4750.2015.05.0388

    LI X B, WEI Y T. An improved Yeoh constitutive model for hyperelastic material [J]. Engineering Mechanics, 2016, 33(12): 38–43. doi: 10.6052/j.issn.1000-4750.2015.05.0388
    [18] 魏志刚, 陈海波. 一种新的橡胶材料弹性本构模型 [J]. 力学学报, 2019, 51(2): 473–483. doi: 10.6052/0459-1879-18-303

    WEI Z G, CHEN H B. A new elastic model for rubber-like materials [J]. Chinese Journal of Theoretical and Applied Mechanics, 2019, 51(2): 473–483. doi: 10.6052/0459-1879-18-303
    [19] 段宇星, 杨强, 赵苗苗, 等. 弹性体材料应变率相关力学行为模型 [J]. 橡胶工业, 2020, 67(12): 899–903.

    DUAN Y X, YANG Q, ZHAO M M, et al. Strain rate-related mechanical behavior model of elastomer material [J]. China Rubber Industry, 2020, 67(12): 899–903.
    [20] CANDAU N, OGUZ O, PEUVREL-DISDIER E, et al. Effect of the strain rate on damage in filled EPDM during single and cyclic loadings [J]. Polymers, 2020, 12(12): 3021. doi: 10.3390/polym12123021
    [21] WANG Y L, LI Z, LI X, et al. Effect of the temperature and strain rate on the tension response of uncured rubber: experiments and modeling [J]. Mechanics of Materials, 2020, 148: 103480. doi: 10.1016/j.mechmat.2020.103480
    [22] 周相荣, 王强, 涂耿伟. 弯曲型橡胶缓冲器冲击试验与数值仿真 [J]. 振动与冲击, 2007, 26(4): 97–100. doi: 10.3969/j.issn.1000-3835.2007.04.023

    ZHOU X R, WANG Q, TU G W. Impact test and simulation for rubber shock absorbers of bending structures [J]. Journal of Vibration and Shock, 2007, 26(4): 97–100. doi: 10.3969/j.issn.1000-3835.2007.04.023
    [23] 林玉亮, 卢芳云, 卢力. 高应变率下硅橡胶的本构行为研究 [J]. 高压物理学报, 2007, 21(3): 289–294. doi: 10.3969/j.issn.1000-5773.2007.03.012

    LIN Y L, LU F Y, LU L. Constitutive behaviors of a silicone rubber at high strain rates [J]. Chinese Journal of High Pressure Physics, 2007, 21(3): 289–294. doi: 10.3969/j.issn.1000-5773.2007.03.012
    [24] 杨建兴, 张江涛, 乔炎亮, 等. 硅橡胶动态压缩性能SHPB测试方法及其本构模型研究 [J]. 高分子通报, 2021(4): 27–34.

    YANG J X, ZHANG J T, QIAO Y L, et al. SHPB test method for the dynamic compressive properties and dynamic constitutive model of silicon rubber [J]. Polymer Bulletin, 2021(4): 27–34.
    [25] TSCHOEGL N W. Constitutive equations for elastomers [J]. Journal of Polymer Science Part A-1: Polymer Chemistry, 1971, 9(7): 1959–1970. doi: 10.1002/pol.1971.150090714
    [26] 卢强, 王占江, 王礼立, 等. 基于ZWT方程的线黏弹性球面波分析 [J]. 爆炸与冲击, 2013, 33(5): 463–470. doi: 10.11883/1001-1455(2013)05-0463-08

    LU Q, WANG Z J, WANG L L, et al. Analysis of linear visco-elastic spherical waves based on ZWT constitutive equation [J]. Explosion and Shock Waves, 2013, 33(5): 463–470. doi: 10.11883/1001-1455(2013)05-0463-08
    [27] 曹侃. 聚碳酸酯动态拉伸力学行为的测试与表征 [D]. 合肥: 中国科学技术大学, 2012.

    CAO K. Experimental investigation and modeling of dynamic tension behavior of polycarbonate [D]. Hefei: University of Science and Technology of China, 2012.
    [28] 吴长河, 冯晓伟, 叶培, 等. 应变率对硫化橡胶压缩力学性能的影响 [J]. 功能材料, 2013, 44(8): 1098–1101. doi: 10.3969/j.issn.1001-9731.2013.08.009

    WU C H, FENG X W, YE P, et al. Effect of strain rate on mechanical properties of vulcanized rubber [J]. Journal of Functional Materials, 2013, 44(8): 1098–1101. doi: 10.3969/j.issn.1001-9731.2013.08.009
    [29] 郭玲梅, 汪洋, 徐伟芳. 硅橡胶拉伸行为的应变率相关性测试和表征 [J]. 高压物理学报, 2019, 33(5): 054101. doi: 10.11858/gywlxb.20180664

    GUO L M, WANG Y, XU W F. Experimental investigation and modeling of strain-rate dependence on tensile behavior of silicone rubbers [J]. Chinese Journal of High Pressure Physics, 2019, 33(5): 054101. doi: 10.11858/gywlxb.20180664
  • 加载中
图(9) / 表(2)
计量
  • 文章访问数:  804
  • HTML全文浏览量:  1134
  • PDF下载量:  103
出版历程
  • 收稿日期:  2021-06-18
  • 修回日期:  2021-07-07
  • 录用日期:  2021-07-16

目录

/

返回文章
返回