Processing math: 100%
  • ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

聊聊动态塑性和黏塑性

王礼立 董新龙

张延炜, 徐景德, 胡洋, 田思雨, 冯若尘, 秦汉圣. 柔性障碍物对甲烷空气爆炸波激励作用的实验研究[J]. 爆炸与冲击, 2021, 41(5): 055402. doi: 10.11883/bzycj-2020-0144
引用本文: 王礼立, 董新龙. 聊聊动态塑性和黏塑性[J]. 爆炸与冲击, 2020, 40(3): 031101. doi: 10.11883/bzycj-2020-0024
ZHANG Yanwei, XU Jingde, HU Yang, TIAN Siyu, FENG Ruochen, QIN Hansheng. Experimental study on incentive effect of flexible obstacle on methane-air explosion wave[J]. Explosion And Shock Waves, 2021, 41(5): 055402. doi: 10.11883/bzycj-2020-0144
Citation: WANG Lili, DONG Xinlong. Talk about dynamic plasticity and viscoplasticity[J]. Explosion And Shock Waves, 2020, 40(3): 031101. doi: 10.11883/bzycj-2020-0024

聊聊动态塑性和黏塑性

doi: 10.11883/bzycj-2020-0024
详细信息
    作者简介:

    王礼立(1934- ),男,教授,博士生导师,wanglili@nbu.edu.cn本文根据作者在2019年全国冲击动力学前沿论坛(2019年12月13~15日,海南万宁)的大会报告《聊聊动态塑性、黏塑性和损伤演化》删减整理而成

  • 中图分类号: O347.1

Talk about dynamic plasticity and viscoplasticity

  • 摘要: 固体力学研究者致力于具有应力-应变本构关系(以下简称为形变型本构关系)的变形体的力学响应研究,而流体力学研究者致力于具有应力-应变率本构关系(以下简称为流动型本构关系)的流动体的力学响应研究。当涉及结构和材料的动态塑性时,到底应该用“塑性变形”还是“塑性流动”来表示?本文从宏观塑性本构理论和微观位错动力学机理两个角度,分别讨论并指出塑性本构关系属于流动型黏塑性率相关本构关系,且同时适用于加载和卸载;因而不应该用应力-应变图来描述塑性加-卸载过程。弹塑性本构关系则是一种形变型和流动型本构关系的耦合。
  • 反应装甲作为对抗聚能射流侵彻的有效装置之一,广泛应用于现代装甲车辆的防护,根据内层材料及其产生的效应不同,可分为爆炸反应装甲(explosive reactive armor, ERA)和被动反应装甲(passive armor)两大类[1-2]。其中爆炸反应装甲(也称平板装药)的典型结构为两层钢板夹层炸药组成的三明治结构,夹层炸药在射流的高速撞击下被引爆,驱动包覆板反向运动切割射流,使其产生断裂、偏折而失去侵彻能力。研究结果[3]显示:射流高速侵彻反应装甲过程中,在撞击点处形成一个高速扩张的孔,由于飞板运动和孔扩张的耦合,背板仅能与穿透反应装甲后射流头部后某处开始作用,未受到飞板干扰作用的射流部分称为逃逸射流。逃逸射流越长,后效穿深越大。因此,如何降低逃逸射流的长度是爆炸反应装甲设计的一个重要方向。研究人员试图采用新的结构设计来减少逃逸射流的长度,赵慧英等[4]通过反应装甲后附加装陶瓷复合装甲结构来提高其对逃逸射流的防护,H.J.Lee等[5]通过在反应装甲后加衬板减小逃逸射流的长度,采用AUTODYN软件模拟了衬板材料对其防护性能的影响,S.Friling等[6]通过在反应装甲后附加被动反应装甲以提高其防护能力,将射流简化为长杆对其作用过程进行了2D模拟,研究结果显示背板具有更高的速度,但是不能反映倾斜条件下与射流的作用过程。

    被动反应装甲典型结构为两层钢板夹层惰性材料,由于具有射流防护效能好、附带损伤效应较小等优点而受到研究人员重视,目前已经有大量针对其防护性能及机理的研究[7-9]。本文中,通过将平板装药与橡胶复合板集成设计,采用实验和数值模拟方法研究其作为面板或背板时对射流的防护性能和机理,并与钢面板的反应装甲进行对比,以期为新型反应装甲的结构设计提供参考。

    为研究橡胶复合板位置对爆炸反应装甲防护性能的影响,设计了以橡胶复合板及钢面板作为面板或背板的3种反应装甲结构,如图 1所示。其中橡胶复合板由两层Q235钢板(厚度为1.2 mm)和硫化橡胶夹层(密度为1.01 g/cm3,厚度为1.5 mm)组成,其长度为150 mm,宽度为50 mm。将钢板表面清洗干净并进行喷砂处理,用调配好的环氧树脂溶液将钢板和橡胶粘结,室温下固化。橡胶复合板的等效钢厚为(0.12×2×7.85+1.5×1.01)/7.85 = 2.55 mm;而爆炸反应装甲包覆板材为Q235钢,厚度为2.5 mm,由此可知,3种结构爆炸反应装甲面密度基本相同。夹层炸药采用厚度为3 mm、装药密度为1.71 g/cm3的B炸药。此外,图 1中结构记号意义如下:“St”表示钢板,“Rubber”表示硫化橡胶层,“E”表示夹层炸药。

    图  1  3种反应装甲结构
    Figure  1.  Three ERA structures

    采用装药口径36 mm聚四氟乙烯塑料壳体的聚能装药对反应装甲作引爆实验,其中铜药型罩壁厚为1 mm,锥角为60°,装药为JH-2。该聚能装药射流头部速度约为6.2 km/s,直径为1.5 mm,炸高为85 mm时对均质装甲钢的平均穿深为150 mm。

    图 2为聚能装药对反应装甲作用的实验布置示意图。实验时聚能装药呈水平放置,口部距测试装甲表面、后效靶板分别为85、210 mm。反应装甲倾角30°,后效靶材为603均质装甲钢,厚度为50 mm,实验后通过测量残余穿深(depth of penetration, DOP)来比较3种装甲结构的防护性能,采用Scandish Flash-XT450脉冲X射线摄影系统对3种结构反应装甲与射流作用情况进行了观测。聚能装药采用电雷管起爆。

    图  2  聚能装药对反应装甲作用的实验布置示意图
    Figure  2.  Experimental scheme of shaped charge and ERA

    图 3为脉冲X射线拍摄的射流与结构(b)和结构(c)作用时的典型时刻的X射线照片,图 4为逃逸射流对靶板表面的损伤情况。表 1为逃逸射流对靶板表面的损伤测量结果。从图 3可以看出,射流发生了偏转,由于稀疏波的影响,飞板边缘速度略低于其它部分;结构(b)的飞板速度约为860~880 m/s,背板与射流作用部位凸起现象不明显,而结构(c)背板与射流作用后呈花瓣形破裂,总体厚度明显大于钢飞板。在结构(c)背板前出现了逃逸射流颗粒,其长度约为6 mm,速度约为3 km/s。

    图  3  典型时刻的X射线实验照片
    Figure  3.  X-ray photographs at t = 46 μs
    图  4  后效靶穿深实验结果照片
    Figure  4.  Experimenal pictures of DOP results
    表  1  实验结果
    Table  1.  Experimental results of penetration
    装甲结构 开坑尺寸/(mm×mm) 开坑深度/mm
    结构(a) 7×11 9
    结构(b) 6×11 11
    结构(c) 6×7 6
    下载: 导出CSV 
    | 显示表格

    图 4可以看出,射流与爆炸反应装甲作用后,在后效靶表面的损伤形成了多个开坑,由射流碎片高速撞击而成,大致可分为两个区域,一个是逃逸射流作用区(如图 4中箭头所示),另外一个区域是背板飞离射流轴线后,射流后部碎片侵彻后效靶形成。橡胶复合板无论作为面板和背板,都可以减小逃逸射流的穿深,作为背板时效果更优,与钢反应装甲相比,穿深降低了46%。

    从实验结果可以看出,橡胶复合板作为爆炸反应装甲面、背板时其防护性能优于钢反应装甲,特别是作为背板时后效穿深相比于钢反应装甲减小了5 mm。

    利用非线性动态有限元ANSYS/LS-DYNA软件ALE算法对聚能装药与反应装甲的作用过程进行了数值模拟,其中空气和聚能装药采用欧拉算法,反应装甲采用拉格朗日算法。根据结构的对称性,建立了1/2计算模型,建模过程中施加对称约束和无反射边界条件。

    JH-2装药采用JWL状态方程和高能材料燃烧模型,夹层炸药Comp.B采用JWL状态方程和Lee-Tarver反应模型[10]描述,其表达式为:

    p=A(1ωR1V)eR1V+B(1ωR2V)eR2V+ωE/V
    (1)
    dλdt=I(1λ)b(ρρ01a)x+G1(1λ)cλdpy+G2(1λ)eλgpz
    (2)

    式中:V = ρ0/ρρ为爆轰产物密度,ρ0为炸药初始装药密度;E = ρ0ee为内能;ABR1R2ω为输入参数;λ为反应速率分数;t为时间,p为压力,IbaxG1cdyG2egz为常数。炸药的主要参数和Lee-Tarver反应模型参数分别如表 2表 3所示。

    表  2  炸药计算参数
    Table  2.  Computational parameters for JH-2 and Comp.B
    炸药 ρ/(g·cm-3) D/(m·s-1) A/GPa B/GPa R1 R2 ω
    JH-2 1.685 8 130 625.3 23.29 5.25 1.6 0.28
    Comp. B 1.715 7 980 524.2 7.77 4.2 1.1 0.50
    下载: 导出CSV 
    | 显示表格
    表  3  Lee-Tarver反应模型参数
    Table  3.  Computational parameters for Lee-Tarver model
    I/s-1 b a x G1/GPa c d y G2/GPa e g z
    4.4×1017 0.667 0 20 310 0.667 0.111 1.0 400 0.333 1.0 2.0
    下载: 导出CSV 
    | 显示表格

    紫铜药型罩和包覆板材料Q235钢板的力学行为分别采用Johnson-Cook模型和Grüneison状态方程进行描述,材料的本构参数见表 4,其中A1B1C1mn为Johnson-Cook模型参数,c0为体积声速,Γ0为Grüneisen系数,s为常数。橡胶夹层和聚能壳体材料采用Grüneison状态方程和Hydro(Pmin)模型描述[1],材料参数取值见表 5,其中σb为抗拉强度,ε为延伸率。

    表  4  紫铜和Q235钢材料的本构方程计算参数
    Table  4.  Computational parameters for copper and Q235 steel
    材料 ρ/(g·cm-3) A1/GPa B1/GPa n C1 m c0/(km·s-1) s Γ0
    Q235 7.85 0.792 0.51 0.26 0.014 1.03 4.57 1.33 1.67
    Cu 8.96 0.090 0.29 0.31 0.025 1.09 3.94 1.49 1.99
    下载: 导出CSV 
    | 显示表格
    表  5  橡胶夹层和聚能壳体材料参数
    Table  5.  Computational parameters for rubber interlayer and polymer shell
    材料 ρ/(g·cm-3) c0/(m·s-1) s Γ0 σb/MPa ε/%
    橡胶 1.01 852 1.865 1.5 20 400
    Teflon 2.15 1 680 1.82 0.59 30 450
    下载: 导出CSV 
    | 显示表格

    图 5给出了射流与3种结构的反应装甲作用形态的数值模拟结果,其中t = 46 μs的全对称模型的数值模拟结果与X射线照片结果符合较好,验证了模型和参数的正确性。当t = 25 μs时,射流刚好穿透结构(a)橡胶复合装甲的背板,由于夹层的存在,射流头部在背板处发生了反射,此刻还未能引爆炸药;对于结构(b)和结构(c),射流则已穿透面板,并引爆了夹层装药。钢板在冲击波和爆炸产物驱动作用下反向运动,后续射流在爆轰波的作用下局部产生了向上的弯曲。当t = 37μs时,射流头部穿透了背板,形成了逃逸射流,结构(b)最长,结构(a)次之,结构(c)最短,这是由于射流在侵彻结构(c)的背板时头部在复合层产生了反射,同时由于结构(c)的背板在运动过程中存在间隙,逃逸射流的后部与背板作用后会“挤入”间隙,因而逃逸射流长度最短。逃逸射流后部与结构(a)和结构(b)的背板作用后,会沿飞板孔壁接触面发生“滑移”,仍然可以逃逸。逃逸射流断裂后形成的颗粒在运动过程中长度发生改变、速度降低,最终3种结构的最长逃逸射流颗粒在触靶前运动速度分别为2.65、2.71和3.10 km/s,长度分别为8.5、12.0和6.5 mm。当t = 58 μs时,橡胶复合装甲两板之间的距离更加明显,射流后部经面板干扰后的偏折角无明显差别,当背板进一步飞离射流运动轴线后,就会失去对射流的作用,在后效靶表面的形成另一个损伤区域。

    图  5  射流与3种结构的反应装甲在典型时刻的作用形态
    Figure  5.  Results of interaction between jet and three structures of ERA at different times

    图 6为飞板的运动速度(vp)的计算结果,其中“F”和“B”分别表示面板和背板,“F-B”表示结构面板的背板(此时结构的面板为复合板),前一个字母是相对于结构整体而言,后一个字母则是相对于复合板而言。由图 6可知,结构(b)钢面板和背板速度曲线基本相同,经过约4 μs的加速过程速度趋于平稳,终了速度为900 m/s。结构(a)和结构(c)钢板的加速历程基本相同,只是方向相反,钢飞板的终了速度约为920 m/s,而橡胶复合板的运动加速过程较复杂,与炸药相邻的钢板(内层板)存在着一个剧烈震荡过程,平均速度约为880 m/s;而外层钢板加速过程比较平稳,终了速度为1 050 m/s。橡胶复合板内、外层具有较大的速度差,其原因可能是炸药爆炸后在内层钢板中产生了较强的冲击波向橡胶层以及外层钢板传播,由于橡胶可压缩性较小,可作为良好的传压介质将冲击波传递给外层板,冲击波经外层板表面反射后产生拉伸波,拉伸波到达外层钢板与橡胶界面时由于不能承受拉应力而产生了“层裂”效应,使外层钢板获得了更高的速度,使逃逸射流长度减少,增加了其防护性能。

    图  6  飞板速度计算结果
    Figure  6.  Calculated plate velocities

    (1) 面密度基本相同条件下,复合板作为面板或背板的反应装甲防护性能优于钢反应装甲,其中橡胶复合板作为反应装甲背板时,防护性能最优。

    (2) 爆炸驱动下橡胶复合板的外层钢板具有更高的速度,相比于钢反应装甲飞板提高约16%。

    (3) 橡胶复合板界面效应和橡胶复合飞板的间隙可有效减小逃逸射流的长度。

  • 图  1  (a) 在τ-γ坐标中表示的Hooke弹性定律(τ=);(b)在τ-˙γ坐标中表示的Newton黏性定律(τ=η˙γ);(c) 在τ-γ坐标中表示的Newton黏性定律(τ=η˙γ

    Figure  1.  (a) The Hooke’s elastic law (τ=) described in τ-γ coordinates; (b) The Newton’s viscous law (τ=η˙γ) described in τ-˙γcoordinates; (c) The Newton’s viscous law (τ=η˙γ) described in τ-γ coordinates

    图  2  理想晶体的剪切滑移

    Figure  2.  Shear slip in a perfect crystal

    图  3  由位错运动形成的滑移。

    Figure  3.  Slip formations due to dislocation movement.

    图  4  一列平行位错的运动造成的宏观塑性切应变

    Figure  4.  The macroscopic plastic shear strain γp(=tanθ) caused by the motion of a row of parallel dislocations

    图  5  位错势垒示意图

    Figure  5.  Schematics of dislocation barrier

    图  6  应力空间中的Mises屈服圆柱和刘氏断裂钟面[14]

    Figure  6.  Mises yielding cylinder and Liu’s bell-like fracture surface in principal stress space[14]

    (1) Yield cylinder (Hencky-Mises); (2) (3) Liu’s bell-like rupture surface; (7) Brittle fracture cone; (8) Plane of pure shear; (9) Liu’s non-fracturing cone

  • [1] 王礼立, 胡时胜, 杨黎明, 等. 材料动力学 [M]. 合肥: 中国科学技术大学出版社, 2017.
    [2] WANG L L, YANG L M, DONG X L, et al. Dynamics of materials: experiments, models and applications [M]. London: Elsevier Science and Technology Press, 2019.
    [3] 李敏华, 王仁. 塑性应力应变关系理论的文献总结 [J]. 力学学报, 1958, 2(2): 167–180.
    [4] 王仁, 黄文彬, 黄筑平. 塑性力学引论: 修订本[M]. 北京: 北京大学出版社, 1992.
    [5] OROWAN E. Plasticity of crystals [J]. Zeitschrift für Physik A, 1934, 89(9-10): 605–659. DOI: 10.1007/BF01341478.
    [6] POLANYI M. Lattice distortion which originates plastic flow [J]. Zeitschrift für Physik A, 1934, 89(9-10): 660–662. DOI: 10.1007/BF01341481.
    [7] TAYLOR G I. The mechanism of plastic deformation of crystals: Part I : theoretical [J]. Proceedings of the Royal Society of London, 1934, A145: 362–387. DOI: 10.1098/rspa.1934.0106.
    [8] OROWAN E. Problems of plastic gliding [J]. Proceedings of the Royal Society of London, 1940, 52(1): 8–22. DOI: 10.1088/0959-5309/52/1/303.
    [9] DAVIDSON D L, LINDHOLM U S. The effect of barrier shape in the rate theory of metal plasticity based on crystal dislocations [C] // Proceedings of the Conference on Mechanical Properties at High Rates of Strain. Oxford: The Institute of Physics, 1974: 124−137.
    [10] KOCKS U F, ARGON A S, ASHBY M F. Thermodynamics and kinetics of slip [J]. Progress in Materials Science, 1975, 19: 1–5. DOI: 10.1016/0079-6425(75)90005-5.
    [11] WANG L L. A thermo-viscoplastic constitutive equation based on hyperbolic shape thermo-activated barriers [J]. ASME Journal of Engineering Materials and Technology, 1984, 106(4): 331–336. DOI: 10.1115/1.3225726.
    [12] ZERILLI F J, ARMSTRONG R W. Dislocation-mechanics-based constitutive relations for material dynamics calculations [J]. Journals of Applied Physics, 1987, 61(5): 1816–1825. DOI: 10.1063/1.338024.
    [13] FOLLANSBEE P S, KOCKS U F. A constitutive description of the deformation of copper based on the use of the mechanical threshold stress as an internal state variable [J]. Acta Metallurgica, 1988, 36(1): 81–93. DOI: 10.1016/0001-6160(88)90030-2.
    [14] 刘叔仪. 关于固体的现实应力空间 [J]. 物理学报, 1954, 10(1): 13–34. DOI: 10.7498/aps.10.13.
  • 期刊类型引用(3)

    1. 王东坡,何思明,李新坡,吴永,孙新坡. 滚石冲击闭孔泡沫铝夹芯板耗能缓冲机理研究. 四川大学学报(工程科学版). 2016(01): 43-49 . 百度学术
    2. 任新见,张庆明,刘瑞朝. 成层式结构泡沫空心球分配层抗爆性能试验研究. 振动与冲击. 2015(21): 100-104 . 百度学术
    3. 田镇华,石少卿,崔廉明,罗伟铭. 钢管在冲击荷载作用下的变形与失效破坏. 后勤工程学院学报. 2014(04): 12-18 . 百度学术

    其他类型引用(5)

  • 加载中
图(6)
计量
  • 文章访问数:  6161
  • HTML全文浏览量:  2310
  • PDF下载量:  280
  • 被引次数: 8
出版历程
  • 收稿日期:  2020-01-15
  • 修回日期:  2020-02-20
  • 刊出日期:  2020-03-01

目录

/

返回文章
返回