考虑温度效应的泡沫铝准静态压缩本构模型

李雪艳 李志斌 张舵

李雪艳, 李志斌, 张舵. 考虑温度效应的泡沫铝准静态压缩本构模型[J]. 高压物理学报, 2018, 32(4): 044103. doi: 10.11858/gywlxb.20170642
引用本文: 李雪艳, 李志斌, 张舵. 考虑温度效应的泡沫铝准静态压缩本构模型[J]. 高压物理学报, 2018, 32(4): 044103. doi: 10.11858/gywlxb.20170642
ZHANG Bingbing, XUE Zhongqing, LEI Yingchun, ZHANG Xizhu, FAN Zhiqiang. Study on Anisotropic Crushing Behavior of the Functionally Gradient Aluminum Foam[J]. Chinese Journal of High Pressure Physics, 2021, 35(2): 024101. doi: 10.11858/gywlxb.20200618
Citation: LI Xueyan, LI Zhibin, ZHANG Duo. Constitutive Model of Aluminum Foams Considering Temperature Effect under Quasi-Static Compression[J]. Chinese Journal of High Pressure Physics, 2018, 32(4): 044103. doi: 10.11858/gywlxb.20170642

考虑温度效应的泡沫铝准静态压缩本构模型

doi: 10.11858/gywlxb.20170642
基金项目: 

湖南省自然科学基金 2017JJ3359

详细信息
    作者简介:

    李雪艳(1991-), 女, 硕士, 主要从事材料的动态力学性能研究.E-mail:15073146797@163.com

    通讯作者:

    李志斌(1985-), 男, 博士, 讲师, 主要从事材料的动态力学性能研究

  • 中图分类号: O347;TJ410

Constitutive Model of Aluminum Foams Considering Temperature Effect under Quasi-Static Compression

  • 摘要: 采用MTS材料试验机研究了不同密度(0.322~0.726g/cm3)的闭孔泡沫铝在温度范围25~500℃下的准静态压缩力学性能,得到了泡沫铝在不同温度下的单轴压缩应力-应变曲线,分析了密度以及温度对其力学行为的影响。利用Liu和Subhash提出的本构模型对不同密度泡沫铝的应力-应变曲线进行拟合,分析并确定了模型中各参数随密度变化的函数,再代入Liu-Subhash模型,得到了泡沫铝的准静态压缩本构模型。通过引入温度软化项对准静态压缩本构模型进行修正,建立了考虑温度效应的泡沫铝准静态压缩本构模型,对闭孔泡沫铝的工程应用具有指导意义。

     

  • 全预制破片类杀伤战斗部的轻质壳体常采用高强度铝合金加工而成,主要用于提高战斗部结构强度、控制战斗部外形等,同时也便于装药过程的实施,而预制破片可布置于壳体内侧或外侧[1]。由于壳体材料与破片材料(如低碳钢、钨合金等)的阻抗存在较大差异,在爆炸加载过程中,不同的阻抗顺序将造成应力波的传播过程出现较大差别[2-3],可能导致破片的初速及完整性发生改变,进而影响战斗部的威力。因此,合理设计预制破片与轻质壳体的阻抗匹配顺序对于该类杀伤战斗部的优化设计具有重要意义。

    采用数值模拟方法研究该问题时,轻质壳体在加载初期由于变形过大而被失效删除,常难以对比不同阻抗顺序造成的差异[4]。因此,最直接的方式是对两种工况下的破片速度进行实际测量,并回收破片进行完整性分析。由于常用的断通靶等电测方式难以直接获得破片初速,一些科研人员基于圆筒结构实验模型,采用更高精度的转镜式高速扫描法及激光干涉测速法获取了爆炸加载过程中不同金属破片的初速,并分析了金属材料与炸药的匹配性[5-6];但该方法不适用于干扰因素较多的预制破片工况,且无法观测破片加速过程中的形态。而脉冲X射线测量方法虽然可观测破片加速过程中的形态,但常由于破片飞散方向杂乱难以准确获得破片的初速。正是由于上述观测及实验方法的限制,导致目前很少见到预制破片与轻质壳体的阻抗匹配方面的报道。

    本研究拟基于预制破片与轻质壳体的典型阻抗排布顺序,依据炸药装药爆轰加载释能特点,设计与杀伤战斗部爆轰加载过程较为相似的滑移爆轰单元结构实验模型(即采用滑移爆轰方式对“单条”破片进行加载),从而获得不同阻抗排布顺序下的破片初速及形态,并结合冲击波在不同介质间的传播规律,分析破片及壳体中的冲击波压力及粒子速度变化过程,以期能为相关杀伤战斗部的设计提供参考。

    图 1显示了杀伤战斗部全预制破片和轻质壳体的两种典型放置工况,即破片置于轻质壳体外侧和内侧。在爆轰加载过程中,破片外置时,爆轰产物首先冲击轻质壳体,使壳体中形成较强的冲击波,到达预制破片界面时,其冲击波主要由低阻抗介质向高阻抗介质传播;而破片内置时,冲击波则由高阻抗介质向低阻抗介质传播。因此,虽然这两种工况的装填比基本相同,但冲击波传播过程存在明显差异,进而可能影响破片的初速及变形。

    图  1  预制破片与轻质壳体的典型放置工况
    Figure  1.  Typical arrangements of prefabricated fragments and light shell

    然而,炸药爆轰加载过程并不是持续的冲击波驱动作用,由于爆轰产物膨胀过程分为高、中、低压3个阶段,且对战斗部破片及壳体的有效加载阶段主要位于高压及中压阶段,因此,一些研究认为炸药装药对破片及壳体的爆轰加载主要依靠两种作用力的共同作用,即冲击波的驱动力和气体爆轰产物的膨胀力[7]。冲击波驱动力可以使破片及壳体在加载初期速度迅速提高,但持续时间较短;而气体爆轰产物膨胀推动力虽然峰值相对偏低,但持续时间较长。以B炸药滑移爆轰驱动∅50 mm圆筒(即标准圆筒试验)为例。图 2显示了圆筒壁的加速历程。为了将两种作用力的作用效果区分开,采用us表示冲击波驱动力单独加载下圆筒壁的速度,ug表示气体产物膨胀力单独加载下圆筒壁的速度,ut(ut=us+ug)为两种作用力共同加载下圆筒壁的速度。从图 2中可以看出,整个加载过程可以明显地划分为两个阶段,即冲击碰撞加速阶段和气体产物膨胀力驱动加速阶段[7-8]。由于气体产物的膨胀力相对偏低,可假定在破片外置和内置两种工况下该阶段的加速效果差异较小,其差异主要在冲击碰撞加速阶段,因此可通过相应的单元实验对该阶段的加载规律进行分析。

    图  2  圆筒试验中圆筒壁膨胀速度-时间历程
    Figure  2.  Time histories of expansionvelocity of cylinder wall

    由于常用的断通靶测量方法在爆炸近场所受干扰因素较多,因此拟采用脉冲X射线系统对破片的分布及加速过程进行观测。从光学成像原理方面分析,脉冲X射线测量过程实际是将三维的实验模型投影到二维平面,因此,为了避免底片中的图像过于杂乱,只能对局部单位模块进行研究。为了尽可能地模拟实际战斗部的爆轰加载方式,采用滑移爆轰方式对“单条”破片进行加载,其中,起爆点为圆柱形装药一端面的中心处(如图 3所示),在柱形装药的中段,其径向冲量接近常量,则沿传爆方向的破片具有相同的加速度,且能够规则排列。此外,考虑到稀疏波的影响,破片距装药首、尾位置均空余一段距离,以保证破片所受爆炸冲量一致。这种加载方式对爆轰产物没有有效的约束,对破片的加载主要依靠冲击波驱动力作用。

    图  3  滑移爆轰驱动破片示意
    Figure  3.  Sketch of sliding detonation-driven fragments

    在滑移爆轰加载过程中,相邻破片运动的初始时间间隔固定,破片的飞散轨迹如图 4所示。根据X射线照片的判读结果可计算出各个破片的速度,其计算公式为

    图  4  破片飞散轨迹示意
    Figure  4.  Flying trajectory of fragments
    vi=ζ(SiSi1)ζL/D=SiSi1L/D
    (1)

    式中:vi为第i个破片的速度,ζ为实物尺寸与图像尺寸的比例,L为底片中相邻破片质心间的距离,Si为底片中第i个破片距装药下表面的距离,D为炸药的爆速。由于(1)式中不显含ζ,所以在获取破片的飞行速度时,不需要考虑底片中ζ的具体值。但有时在实验底片中无法显现装药的下表面,故难以判读破片运动的真实距离,只能获得其相对距离,同时考虑到不同破片的速度有时会出现波动,因此,实际分析时可根据(2)式计算破片的初速,即

    v=ni=2(SiSi1)(n1)L/D=DSnS1(n1)L=Dk
    (2)

    式中:n表示破片的总数,Sn-S1和(n-1)L分别表示底片中第n个破片与第1个破片之间的纵向和横向的距离。因此,(2)式也可以表示为:破片速度等于底片中破片曲线的斜率k与炸药爆速的乘积。

    实验时,破片及轻质壳体的材料、尺寸均根据实际战斗部进行选择。破片采用93钨合金方形破片,尺寸为6.5 mm×6.5 mm×3.8 mm,质量约为2.81 g;轻质壳体材料选用LY12铝合金,厚度为1.8 mm,宽度略大于破片。将破片排布并粘接在条状铝合金板上,然后再与炸药侧面粘接,以模拟破片内置及外置两种工况。实验所用主装药均为B炸药,密度为1.65 g/cm3,爆速约为7.8 km/s。由于装药的长度大于底片盒的宽度,故实验时将装药竖立放置,破片向侧向运动(实验布局如图 5所示);同时,由于采用的数字成像板不能二次曝光,图像中难以同时出现装药的基准线,因此起爆前先拍一发静止相片以判断装药的轴线是否与竖直悬挂的底片平行,然后再将成像板重新装入底片盒进行实验。

    图  5  实验布局
    Figure  5.  Experiment layout

    图 6为爆轰加载实验所获得的X射线摄影底片。图 6(a)显示了破片外置时钨破片和铝壳加速过程中的状态,其中并没有独立显示出铝壳体的迹线,且钨破片的下边界不太整齐,表明铝壳体与钨破片粘接在一起,两者速度基本一致;钨破片上方有部分细小的碎渣,其图像灰度与钨破片近似,表明爆轰加载过程中钨破片发生了轻微的破损。图 6(b)显示了破片内置时钨破片与铝壳加速过程中的状态,可以看出,钨破片下表面较为完好,但上表面处既存在一些细小的碎渣,同时还有较完整的条状碎裂体,其图像灰度与钨破片近似,可能是由于钨破片直接与炸药接触,其压力较高致使破片发生了较明显的层裂;钨破片图像上方出现了一条灰度相对较浅的线条,这是铝壳体的迹线(X射线摄影中,图像灰度与物质的密度相关),表明加载过程中铝壳体的速度明显高于钨破片。

    图  6  实验所获底片
    Figure  6.  Experimental film

    底片图像的横向距离(设坐标为x)可用于计算时间,纵向距离(设坐标为y)可用于计算破片的位移。由于每张底片的放大比不同,所以这些数据仅表示在底片中所处的位置,其数值并没有实际的物理意义,仅用于计算曲线的斜率。采用最小二乘法获得曲线的斜率k后,再按照(2)式计算破片的初速v

    图 7为底片的判读结果。对数据进行拟合计算时,考虑到边侧少量破片在加载时受到稀疏波等因素的影响,其速度明显偏低,故将其删去后再拟合,所获结果列于表 1

    图  7  底片判读结果
    Figure  7.  Results from experimental film
    表  1  不同工况下的破片及壳体速度
    Table  1.  Velocity of fragments and shell under different conditions
    Working condition Shell Fragments
    k v/(m·s-1) k v/(m·s-1)
    Tungsten fragments situated outside aluminum shell 0.126 983 0.126 983
    Tungsten fragments situated inside aluminum shell 0.143 1 115 0.110 858
    下载: 导出CSV 
    | 显示表格

    表 1中可以看出,破片内置时,钨破片的速度明显偏低,铝壳体的速度相对较高,且破片与壳体的动能之和也低于破片外置时的状态。这可能是因为此时钨破片直接与炸药接触,其上表面发生了相对严重的层裂(如图 6(b)所示),其碎裂部分的质量较大,且速度相对偏高,因此从整体上看,破片外置时,钨破片的完整性及速度均具有明显优势。轻质壳体和破片的完整性及初速均涉及其受力过程,故对应力波在这两种介质中的传播过程进行分析。

    为了简化分析过程,不考虑应力波在同一介质中的衰减,主要分析在不同介质间传播时的应力变化。分析过程中,首先需要描述爆轰产物作用过程的p(u)曲线以及铝合金壳体、钨合金破片的p(u)曲线,即波阵面压力与波阵面后粒子速度的关系曲线。若炸药以滑移爆轰方式驱动壳体及破片,则爆轰产物的p(u)曲线可根据(3)式计算[9]

    u=2Dγγ21[1(ppCJ)γ12γ]
    (3)

    式中:pCJ=ρ0D2/(γ+1),ρ0D分别为炸药密度和爆速;B炸药的多方指数γ=2.706[9]。铝合金和钨合金的p(u)曲线根据(4)式计算

    P=ρu(c0+λu)
    (4)

    式中:ρc0分别为金属材料的密度及声速,λ为本构方程参数。其具体数值列于表 2

    表  2  金属材料参数[9]
    Table  2.  Parameters of metals[9]
    Material ρ/(g·cm-3) c0/(km·s-1) λ
    Tungsten alloy 17.5 3.832 1.497
    Aluminum alloy 2.785 5.238 1.338
    下载: 导出CSV 
    | 显示表格

    图 8为应力波在不同介质中传播过程示意图,其中:3条实线分别表示B炸药爆轰产物的等熵线以及铝合金壳体、钨合金破片的p(u)曲线,爆轰产物的等熵线与金属材料p(u)曲线的交点坐标即为金属材料中的入射冲击波压力和波阵面后的粒子速度;金属介质间应力波的传播路径主要沿冲击材料在该状态点的反演线(如图 8中虚线所示),其计算表达式为[10]

    图  8  不同介质间应力波传播分析(①~⑥为界面编号)
    Figure  8.  Analysis of stress wave transmission betweendifferent materials (①-⑥:Interface No.)
    p=ρ(2u0u)[c0+λ(2u0u)]
    (5)

    式中:u0表示碰撞前的粒子速度。反演曲线与被冲击材料p(u)曲线的交点即表示碰撞后的状态。此外,由于空气的密度及声速均远低于金属材料,因此图 8中采用横坐标轴近似其p(u)曲线。

    图 8表 3可知, 破片置于轻质壳体外侧时:钨破片一端受到从壳体传入的冲击波作用,另一端受到空气反射的拉伸波作用,使其速度容易快速提升;而轻质壳体一端受到爆轰产物传入的冲击波作用,另一端受到钨破片反射的冲击波作用,两个作用方向相反,限制其速度过快提升,不会过多地占用炸药所转换的动能,同时没有拉伸波作用,轻质壳体也不容易快速破裂,在实际战斗部结构中可以增加对爆轰产物的约束,有利于对后期爆轰产物膨胀驱动能量的利用。

    表  3  材料界面处的状态参数
    Table  3.  State parameters of material interfaces
    Interface No. u/(km·s-1) p/GPa
    0.931 16.82
    0.385 29.70
    0.770
    0.310 23.30
    0.502 8.26
    1.004
    下载: 导出CSV 
    | 显示表格

    破片置于轻质壳体内侧时:虽然钨破片受到爆轰产物的冲击波及铝合金壳体反射的拉伸波作用,但从整体上看,其界面的平均粒子速度低于破片外置时的状态;而铝合金壳体由于外侧没有高阻抗介质阻挡,所受到的拉伸波作用可使其速度快速提升,但也易造成过早破裂。该工况在一些具有特定要求的战斗部中具有某些方面的优势,设计时可考虑在破片与装药之间增添一些性能较好的缓冲体,以提升破片的完整性;在轻质壳体与破片之间,也可适当增设缓冲材料,以减小轻质壳体与破片之间的速度差,不仅可以避免壳体过多地占用能量,也可减小壳体所受到的冲击力,延缓其破裂时间,从而增加对爆轰产物的约束时间,提升爆轰产物驱动能量的利用率。当然,本研究只针对冲击碰撞阶段的加速过程,未进行两个加速阶段的系统研究,也未涉及轻质壳体的动态破裂机理问题,还需要进一步深入研究。

    (1) 预制破片置于轻质壳体外侧时,初始应力波由低阻抗金属材料向高阻抗金属材料传播,使破片受到壳体传入的冲击波及空气传入的拉伸波作用,其速度易迅速提升,且破损轻微;而轻质壳体没有受到拉伸波作用,不易快速破裂,对后期爆轰产物膨胀驱动能量的利用较为有利。

    (2) 预制破片置于轻质壳体内侧时,初始应力波主要由高阻抗金属材料向低阻抗金属材料传播,虽然破片受到爆轰产物传入的冲击波及壳体反射的拉伸波作用,但其初速低于破片外置时的状态,且易发生破损,甚至有明显层裂现象;而铝合金壳体由于受到空气传入的拉伸波作用,其速度可快速提升,但也易造成其过早破裂,可能减少对爆轰产物的约束时间,降低炸药驱动能量的利用率。

  • 图  常温下不同密度泡沫铝的准静态应力-应变曲线

    Figure  1.  Quasi-static stress-strain curves of aluminum foam with different densities at room temperature

    图  常温下不同密度泡沫铝的应力-应变曲线拟合情况

    Figure  2.  Fitting of stress-strain curves of aluminum foam with different densities at room temperature

    图  常温下参数随密度的变化规律

    Figure  3.  Variation of parameters with density at room temperature

    图  不同温度下泡沫铝的准静态应力-应变曲线

    Figure  4.  Quasi-static stress-strain curves of aluminum foam at different temperatures

    图  不同温度下泡沫铝应力-应变曲线的拟合情况

    Figure  5.  Fitting of stress-strain curves of aluminum foam at different temperatures

    表  1  不同密度下的模型参数值

    Table  1.   Parameter values for different densities

    Density/(g·cm-3) Parameter
    p1 p2 p3 p4 p5
    0.322 3.08 86.09 86.22 -6.49 11.60
    0.481 5.98 90.83 90.74 -4.66 10.98
    0.639 9.31 91.20 90.89 -3.00 9.41
    0.726 12.01 92.04 92.14 -1.33 7.41
    下载: 导出CSV
  • [1] 刘培生, 李铁藩, 傅超.多孔金属材料的应用[J].功能材料, 2001, 32(1):12-15. http://www.cqvip.com/QK/72003x/201607/epub1000000248429.html

    LIU P S, LI T F, FU C.Application of porous metal materials[J].Functional Materials, 2001, 32(1):12-15. http://www.cqvip.com/QK/72003x/201607/epub1000000248429.html
    [2] 尚朝秋, 王应武.泡沫铝材料研究现状分析[J].云南冶金, 2016, 45(3):10-13. http://www.cqvip.com/QK/90321B/201603/669250357.html

    SHANG Z Q, WANG Y W.The research status analysis of aluminum foam material[J].Yunnan Metallurgy, 2016, 45(3):10-13. http://www.cqvip.com/QK/90321B/201603/669250357.html
    [3] CHEN C, LU T J.A phenomenological framework of constitutive modelling for incompressible and compressible elastic-plastic solids[J].International Journal of Solids and Structures, 2000, 37(52):7769-7786. doi: 10.1016/S0020-7683(00)00003-2
    [4] 王二恒, 虞吉林, 王飞, 等.泡沫铝材料准静态本构关系的理论和实验研究[J].力学学报, 2004, 36(6):673-679. http://mall.cnki.net/magazine/Article/LXXB200406004.htm

    WANG E H, YU J L, WANG F, et al.A theoretical and experimental study on the quasi-static constitutive model of aluminum foams[J].Acta Mechanica Sinica, 2004, 36(6):673-679. http://mall.cnki.net/magazine/Article/LXXB200406004.htm
    [5] WANG Z H, JING L, ZHAO L M.Elasto-plastic constitutive model of aluminum alloy foam subjected to impact loading[J].Transactions of Nonferrous Metals Society of China, 2011, 21(3):449-454. doi: 10.1016/S1003-6326(11)60735-8
    [6] HAKAMADA M, NOMURA T, YAMADA Y, et al.Compressive deformation behavior at elevated temperatures in a closed-cell aluminum foam[J].Materials Transactions, 2005, 46(7):1677-1680. doi: 10.2320/matertrans.46.1677
    [7] ALY M S.Behavior of closed cell aluminum foams upon compressive testing at elevated temperatures:experimental results[J].Materials Letters, 2007, 61(14):3138-3141. http://www.sciencedirect.com/science/article/pii/S0167577X06013206
    [8] CADY C M, Ⅲ G T G, LIU C, et al.Compressive properties of a closed-cell aluminum foam as a function of strain rate and temperature[J].Materials Science and Engineering A, 2001, 525(1):1-6. http://www.sciencedirect.com/science/article/pii/S0921509309007503
    [9] 习会峰, 刘逸平, 汤立群, 等.考虑温度效应的泡沫铝静态压缩本构模型[J].哈尔滨工程大学学报, 2013, 34(8):1000-1005. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hebgcdxxb201308010

    XI H F, LIU Y P, TANG L Q, et al.Constitutive model of aluminum foam with temperature effect under the quasi-static compression[J].Journal of Harbin Engineering University, 2013, 34(8):1000-1005. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hebgcdxxb201308010
    [10] 王鹏飞, 徐松林, 胡时胜.基于温度与应变率相互耦合的泡沫铝本构关系[J].高压物理学报, 2014, 28(1):23-28. doi: 10.11858/gywlxb.2014.01.004

    WANG P F, XU S L, HU S S.A constitutive relation of aluminum foam coupled with temperature and strain rate[J].Chinese Journal of High Pressure Physics, 2014, 28(1):23-28. doi: 10.11858/gywlxb.2014.01.004
    [11] LIU Q, SUBHASH G.A phenomenological constitutive model for foams under large deformations[J].Polymer Engineering and Science, 2004, 44(3):463-473. doi: 10.1002/(ISSN)1548-2634
  • 加载中
图(5) / 表(1)
计量
  • 文章访问数:  7095
  • HTML全文浏览量:  2993
  • PDF下载量:  277
出版历程
  • 收稿日期:  2017-09-18
  • 修回日期:  2017-11-03

目录

/

返回文章
返回