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

恒定高应变率拉伸条件下泡沫金属力学性能

张晓阳 谭仕锋 刘泽宇 赵飘

马彦, 袁福平, 武晓雷. 双相高强钢FeNiAlC的动态剪切行为及微结构机理[J]. 爆炸与冲击, 2021, 41(1): 011404. doi: 10.11883/bzycj-2020-0224
引用本文: 张晓阳, 谭仕锋, 刘泽宇, 赵飘. 恒定高应变率拉伸条件下泡沫金属力学性能[J]. 爆炸与冲击, 2024, 44(1): 013105. doi: 10.11883/bzycj-2023-0128
MA Yan, YUAN Fuping, WU Xiaolei. Dynamic shear behaviors and microstructural deformation mechanisms in FeNiAlC dual-phase high strength alloy[J]. Explosion And Shock Waves, 2021, 41(1): 011404. doi: 10.11883/bzycj-2020-0224
Citation: ZHANG Xiaoyang, TAN Shifeng, LIU Zeyu, ZHAO Piao. Mechanical property of metallic foams under dynamic tension with constant high strain rate[J]. Explosion And Shock Waves, 2024, 44(1): 013105. doi: 10.11883/bzycj-2023-0128

恒定高应变率拉伸条件下泡沫金属力学性能

doi: 10.11883/bzycj-2023-0128
基金项目: 国家自然科学基金项目(11902140);湖南省自然科学基金项目(2018JJ3425);
详细信息
    作者简介:

    张晓阳(1988- ),男,博士,讲师,xyzhang@usc.edu.cn

  • 中图分类号: O347.1

Mechanical property of metallic foams under dynamic tension with constant high strain rate

  • 摘要: 为了探究泡沫金属恒定应变率动态拉伸力学行为,基于3D Voronoi模型,采用双向拉伸加载方式和1.55倍等效胞孔直径高度的试件,实现了5000 s−1恒定高应变率动态拉伸条件下泡沫金属力学性能测试数值模拟实验,模拟结果显示:动态拉伸过程满足应力均匀性和变形均匀性要求,且试件破坏位置合理;在恒定应变率(0.5~5000 s−1)动态拉伸时,泡沫金属的破坏应变基本不受应变率的影响;当应变率不超过500 s−1 时,破坏应力受应变率影响很小,当应变率在 500~5000 s−1 时,破坏应力随着加载速率的增大而线性增大。
  • 高应变速率载荷下,金属材料往往因局部化形变而失效破坏[1]。在军事工业领域,子弹的穿甲侵彻与材料的局域化绝热剪切损伤紧密相关。对实心穿甲弹,弹头的绝热剪切有利于在侵彻过程中自锐化而增大穿刺深度[2]。在航空航天领域,飞行器的安全威胁主要来自空间碎片的高速撞击,材料的绝热剪切直接导致飞行器构件的灾难性破坏。所以,金属材料的动态剪切行为是材料在特殊使用环境下必须考虑的关键问题[1-3]

    高应变速率载荷诱发材料的局部化形变,有两个基本特征:加载速度快变形过程接近绝热状态;变形区域小易形成剪切带。由于绝热剪切带(adiabatic shear band, ASB)的形成和扩展在约105 s内完成且宽度一般在102 μm内,目前仍缺乏有效的实验技术实现原位直接观测[2]。现阶段,主要手段是通过表征材料变形失效后的微观结构来分析和推演ASB形成和扩展过程。Meyers等[4]在低层错能304L不锈钢动态剪切实验中,发现ASB内形成晶粒尺寸100~200 nm的等轴晶,通过形变织构分析,得出这些等轴晶是通过变形过程中初始晶粒的旋转、滑移和孪生机制形成的。Yang等[5]则通过应变冻结实验,系统地研究1Cr18Ni9Ti奥氏体型不锈钢形成ASB过程中微结构的演化,证实ASB内晶粒的动态旋转再结晶机制。Li等[6]在高锰钢的子弹侵彻试验中,发现ASB内产生明显的纳米晶区域和非晶区域,表明在ASB内由于极高的温升产生了金属熔化和快速冷却。Wei等[7-8]则分别研究了粗晶(coarse grained, CG)、超细晶(ultrafine grained, UFG)、纳米晶(nanocrystalline, NC)结构铁在高应变率下的压缩变形,结果表明,减小晶粒尺寸会导致材料形成ASB的概率提高,同时定量分析了应变硬化率和率敏感性对形成ASB的影响。

    目前,关于金属材料ASB失效的基本准则是由Zener等[9]提出的材料应变硬化与热软化相互竞争的机制,即当材料丧失应变硬化能力时发生ASB失效。Batra等[10]提出ASB形成判据(最大应力准则:ASB在流变应力达到最大值时形成并快速扩展,导致流变应力快速跌落至峰值的80%),但该准则往往仅适用于均质材料。在梯度材料[11]和多层材料[12]中,ASB形成于材料中硬化能力较弱的区域,但由于材料的非均匀性,当局部形成ASB后材料整体依然具有硬化能力,当ASB贯穿整个样品时材料才会丧失承载能力。在这些异构材料[11-12]中,局部微区的力学属性差异能够延缓ASB的形成和传播,从而延缓材料失效的发生。由此可见,材料的微观结构对形成ASB失效非常重要。本文中,以FeNiAlC双相钢为例,讨论不同微结构材料的剪切变形行为,进一步理解微结构对材料形成ASB的影响以及温升和热软化对材料发生局部化变形的影响。

    FeNiAlC合金由纯度为99.99%的铁、镍、铝金属单质进行熔炼而成。首先,在真空电弧炉中经过5次反复熔炼形成20 mm的铸锭,测得铸锭中铁以外化学成分的质量分数为w(Ni)=24.86%, w(Al)=5.8%, w(C)=0.38%,记为Fe-24.86Ni-5.8Al-0.38C。接着,将铸锭加热至1 100 ℃后热锻成10 mm×10 mm×50 mm的钢板,在950 ℃下进行热轧,终轧厚度为6 mm,冷却至室温。然后,将热轧板在真空马弗炉中进行均匀化固熔,固熔温度为1 150 ℃,固熔24 h后进行快速淬火处理。最后,将固熔后的金属板进行大应变冷轧处理,压下量(截面积减少量)分别为39%、53%(记为CR39、CR53),未进行冷轧处理的材料为粗晶态(记为CG)。

    动态冲击样品采用强迫剪切构形试样(帽形试样),强迫剪切区宽度为200 μm,剪切方向垂直于板材轧制方向,帽形试样的厚度即为冷轧后板材的厚度,如图1所示。实验时,将试样紧密放置于用马氏体时效钢制成的限位环中,限位环在冲击过程中变形量极小,保证样品主要受剪应力的作用。由有限元模拟结果[13],剪切区的高度(2 mm)远大于剪切区的宽度(0.2 mm),压应力相对剪应力是小量,因此帽形样品的变形过程可以近似等效为纯剪切变形。在13 mm分离式霍普金森压杆(split Hopkinson pressure bar, SHPB)装置上进行冲击实验,子弹速度约20 m/s,加载应变率范围为10 000~30 000 s−1,每个状态的样品重复5次实验,重复性相对较好,实验误差小。通过杆上的应变片和动态应变仪测量,得到入射杆和透射杆中的反射波应变εR和透射波应变εT,进行空载标定应变片敏感系数。利用两波法和一维弹性理论,计算样品在一维应力加载下的剪切应力τs、剪切应变γs、应变率˙γs与和位移us

    图  1  动态剪切实验装置和样品
    Figure  1.  Dynamic shear experimental device and its sample
    {τs=E(A/As)εTus=2c0t0εRdtγs=Us/ts˙γs=2c0tsεR
    (1)

    式中:EAc0分别为霍普金森杆的杨氏模量(210 GPa)、横截面积(132.7 mm2)和纵波波速(5 050 m/s),Asts分别为帽形样品的剪切区面积和宽度[14]

    实验前后样品的微观结构采用Zeiss Gemini 300扫描电镜(scanning electron microscope, SEM)和JEOL 2100F透射电镜(transmission electron microscope, TEM)进行表征。使用电子背散射衍射(electron backscatter diffraction, EBSD)技术获取材料的晶体学取向及位错密度等,选用Oxford Symmetry探头进行扫描,扫描步长为50~1 000 nm,使用CHANNEL 5软件进行数据后处理。观测前的样品准备,详见文献[15]。

    冲击实验前,对CG、CR53样品进行微结构表征。图2为CG样品晶粒的反极图(inverse pole figure, IPF)和相图(phase figure, PF)。固熔后样品由奥氏体相(austenite, γ)和马氏体相(martensite, α)构成,奥氏体平均晶粒尺寸较大,约为11.5 μm,马氏体相析出在奥氏体晶界处,平均晶粒尺寸约为3.55 μm,马氏体相约占15%,属于一种双相合金。由图2(a)可见,固熔后晶粒取向随机分布,无明显织构。图3为CR53样品的IPF和PF。轧制过程中发生少量马氏体相变,马氏体相比例略有增加,约占18%。奥氏体晶粒被拉长,统计晶粒的长轴和短轴比约为2.1,据晶粒被拉长计算等效应变约为0.49[16],而冷轧53%换算成等效应变约为0.50。这说明冷轧过程中,变形主要发生在奥氏体相中,马氏体相几乎不承担变形。由图3(a)可见,与粗晶态相比,冷轧后出现明显的〈111〉丝织构[17]。为了研究样品冷轧前后的位错形貌和位错密度,对CG、CR53样品进行TEM观察。图4(a)为CG样品奥氏体相的形貌和[011]晶向衍射谱。固熔后样品内部有少量的平形位错,长度较短,位错密度极低。图4(b)为CR53样品的晶粒衍衬像、α相[001]晶向的衍射谱和γ相[011]晶向的孪晶衍射谱,样品内部位错密度明显增大,位错大多缠结在晶界和相界附近,奥氏体晶粒内部产生少量变形孪晶。CG和CR53样品中晶粒形貌、尺寸和位错密度的差异,影响两种微结构样品的动态剪切行为,这将在以下作详细介绍。

    图  2  实验前CG样品的微观结构
    Figure  2.  Microstructures of CG sample before experiment
    图  3  实验前CR53样品的微观结构
    Figure  3.  Microstructures of CR53 sample before experiment
    图  4  实验前CG和CR53样品的TEM
    Figure  4.  TEM observations of CG and CR53 sample before experiment

    用动态剪切实验所得的数据,通过式(1)获得的剪切应力与变形位移、剪切应变的关系曲线如图5所示。圆点为名义动态剪切屈服强度,CG样品约为480 MPa,CR39、CR53样品分别为910、1 250 MPa,冷轧量越大,样品的剪切屈服强度越高。方点为样品的最大剪切应力,根据最大剪应力准则[10],此时样品剪切区内部形成ASB,此后剪切应力快速下降,样品发生失稳破坏。在形成剪切带前,变形主要集中在宽200 μm的剪切区内,一般认为这个过程中剪切区内的变形是均匀的[14]图5(b)中,曲线所包围的面积可以近似为材料的动态剪切韧性。

    图  5  不同微结构的动态剪切性能
    Figure  5.  Dynamic shear properties of various microstructures

    与其他合金一样,利用帽形试样获得FeNiAlC合金的动态剪切性能,并进行比较。图6为FeNiAlC合金和其他常用合金(301、316不锈钢等)的均匀剪切应变和剪切韧性与动态剪切强度的关系[13-14, 18-24]。在同等动态剪切强度情况下,FeNiAlC合金剪切韧性优于其他合金,拥有较优异的动态剪切性能。

    图  6  不同金属的动态剪切性能
    Figure  6.  Dynamic shear properties of various metals

    为了研究变形过程中微结构的变化,采用SEM中的EBSD技术对变形后的样品进行表征。由图7(a)(b)可见,CG样品中剪切区的中心有一条识别率较低的窄带(两条虚线之间)是剪切失稳后形成的ASB。ASB两侧约100 μm区域是主要的剪切变形区,剪切变形区内产生大量小角晶界,在奥氏体和马氏体的相界附近有大量的小晶粒产生。在等轴的奥氏体晶粒内形成大量的变形孪晶(deformation twin, DT),宽度为2~3 μm。由图7(b)可见,马氏体相所占比例约14%,与剪切变形前相比几乎没有变化,说明变形过程中没有发生马氏体相变。在图7(c)中,红线为该点与相邻点之间的取向差,蓝线为该点与初始点的取向差。每跨过一个(孪晶)界面时,取向角改变60°,这是典型的FCC晶体Σ3孪晶取向关系[25]。CG样品剪切变形后孪晶密度较高,发生孪生变形的区域占总体积的60%以上,说明在此应力状态下,奥氏体晶粒大部分发生孪生变形。

    图  7  实验后CG样品的微观结构
    Figure  7.  Microstructures of CG sample after experiment

    图8(a)(b)可见,CR53样品中靠近ASB区域产生大量的小晶粒,与变形前相比晶粒明显细化,远离ASB区域产生少量的孪晶。由图8(b)可见,α相所占比例约18%,与变形前相近,说明冷轧态样品也几乎没有发生马氏体相变。由图8(c)可见,当跨过孪晶界面时,CR53样品界面两侧取向角相差60°,同CG样品结果类似。

    图  8  实验后CR53样品的微观结构
    Figure  8.  Microstructures of CR53 sample after experiment

    KAM (kernel average misorientation)即mka,为晶粒内一点与周围相邻点的取向差的平均值,是一个与几何必需位错(geometrically necessary dislocation, GND)密度正相关的物理量[26]图910分别为CG和CR53样品剪切变形前后局部晶粒的KAM及分布。分布图中,统计分布时取向差步长为0.25°。对比CG样品剪切变形前后局部晶粒的KAM,可以得到两个有价值的结论。(1)变形前,由于经过高温长时间的固熔处理,晶粒内部取向差极小,说明此时晶粒内部位错密度极低,仅在晶界周围小于1 μm内存在明显的取向差分布(0.5°~1.1°)。这可能是因固熔后进行快速淬火处理,在界面处会产生残余应变,引起GND密度升高。(2)与变形前相比,剪切变形后奥氏体晶粒承担主要变形,因此在奥氏体晶粒内部出现了高密度的GND,KAM在0.8°~2.0°,尤其在孪晶界周围,KAM约为2°。说明变形产生大量的位错钉扎在孪晶界附近,位错与孪晶的相互作用为CG样品在剪切变形过程中提供了较大的应变硬化能力。对比CG态样品变形前后γ相KAM分布可见,由于样品中α相含量较少,硬度较大[15],几乎不参与变形,变形前后KAM几乎没有变化,所以这里只进行γ相KAM的比较。变形前KAM平均值约为0.96°,变形后为1.26°,根据KAM与GND密度的关系[26],可以估计变形后GND增加约30%。这也说明,CG样品均匀剪应变之所以能达到5,晶粒内部有充足的可动位错参与变形。

    图  9  实验前CG和CR53样品的KAM及分布
    Figure  9.  KAM values and distributions of CG and CR53 sample before experiment
    图  10  实验后CG和CR53样品的KAM及分布
    Figure  10.  KAM values and distributions of CG and CR53 sample after experiment

    对比CR53样品剪切变形前后局部晶粒的KAM可见,由于经过大应变冷轧过程(等效应变约0.50),样品表现出很高的KAM,尤其是在晶界附近KAM为1.8°~3.0°,远高于固熔后的CG样品。不同于CG样品,CR53样品变形后晶粒内部的KAM明显小于变形前的KAM。同样,对比CR53样品变形前后γ相KAM分布可见,变形前分布范围为0.8°~3.1°,变形后减小为0.8°~2.3°,变形前平均值为1.73°,变形后减小至1.33°,根据KAM与GND密度的线性关系,GND密度减小约23%。通常,变形会使位错密度增大,但由于CR53样品变形前的位错密度已经很高,且动态实验应变率高,变形时间短,剪切区内发生绝热剪切,有较大的温升,会使局部晶粒发生部分回复。同时,该样品均匀剪切应变较小,约为1,所以剪切变形后,回复的作用大于变形的效果,KAM不增反降,GND密度降低,致使CR样品应变硬化能力降低,过早的发生塑性失稳(形成ASB)。

    在动态加载过程中,由于加载速度极高,变形时间在100 μs内,热传导来不及进行,所以整个剪切过程可近似成绝热过程[1-3]。同时,认为在均匀剪切阶段,剪切区内的应变是均匀的。计算剪切过程中剪切区温升的塑性功转化温升公式为[14]

    ΔT=ηρcVγ0τdγp
    (2)

    式中:η为总的非弹性功转化为热量的系数(Taylor-Quinney系数),本文中取0.8;ρcV为材料的密度和比热容,τγp为剪切应力和剪切塑性应变。对FeNiAlC合金,测得ρ为9.8 g/cm3cV取540 J/(kg·K)[27]图11为用实验所测的剪切应力和应变、据式(2)计算的温升曲线。CR样品的温升斜率明显高于CG样品,这是因为CR样品由于冷轧过程中引入大量的可动位错及位错胞等亚结构,在塑性变形初期剪切流变应力大,导致温升的斜率高于CG样品。当剪应变达到0.5时,CR样品的温升约为100 K,CG样品只有50 K,而当剪应变达到最大均匀剪应变时,CR53样品最终温升达到190 K,而CG样品最终温升高达720 K。这表明,绝热温升可能不是ASB形成的主导或者唯一因素[27-28]

    图  11  最大应力点前剪切区塑性功引起的温度升高
    Figure  11.  Temperature rise due to plastic dissipation work in shear zone before maximum stress point

    一般认为,ASB是材料受到冲击载荷作用时的一种失效模式[1-2]。当材料受到高速率剪切力时,通常认为材料的失效是塑性应变硬化和绝热温升所致的热软化相互竞争的结果[9]。变形初期,温升较小,材料的应变硬化占据主导,发生均匀变形,剪切应力不断提高。当剪切区温升达到一定值时,材料的应变硬化不足以抵抗温升带来的热软化,材料发生局部化变形,剪切应力骤然下降,形成ASB。图78中未能解析的带状部分就是变形形成的ASB。因为带内发生了比带外高1~2个数量级的变形,所以当前的扫描步长不足以解析其微结构。对于ASB内的微结构采用40 nm的步长重新扫描,如图12所示。CG样品ASB内晶粒不同于剪切区的晶粒发生孪生变形或局部晶粒细化,而是形成了大小几乎相等的等轴晶。晶粒尺寸集中在200~500 nm,远远小于初始的晶粒尺寸(13 μm)。由此可得,在ASB形成时带内组织发生了动态再结晶(dynamic recrystallization, DRX),与大多数金属材料类似[14],在大应变和高温的作用下,ASB内形成取向随机、大小相等的再结晶晶粒。同CG样品相似,CR53样品同样发生了DRX。CG和CR53两种样品ASB内虽然应变极大,但均不发生相变,只是在原有马氏体相的基础上发生晶粒细化,没有形成新的马氏体相。

    图  12  绝热剪切带内微观结构
    Figure  12.  Microstructures within ASBs

    根据热塑模型理论[1],材料发生ASB失效时,剪切带的宽度可以通过材料所受最大剪应力和应变率等进行估算。Dodd等[1]提出剪切带宽度的计算公式:

    δ=2(λTτ˙γ)1/2
    (3)

    式中:τ˙γ分别为形成剪切带时的剪应力和剪应变率,T为试样温度,λ为热传导系数(不锈钢约为15 W/(m·K)[28]。基于前面的实验和计算结果,CG样品剪切带的理论宽度为29.7 μm(T=1 018 K,τ=1 100 MPa,˙γ=63 000 s−1),CR53样品的为14.6 μm (T=468 K,τ=1 340 MPa,˙γ=98 000 s−1)。图78中的剪切带宽度并不是均匀的,受到局部微结构和受力状态的影响较大。通过划线法统计剪切带的宽度,结果如图13所示。CG样品剪切带的宽度范围为12~19 μm,平均值为15.8 μm,而CR53样品剪切带宽度范围为9~15 μm,平均值为12.3 μm。对比来看,CR53样品的理论计算ASB宽度更接近于实验所测真实宽度,而CG样品的实际剪切带宽度范围小于理论计算值。这可能是根据塑性功计算样品温度带来的影响。热塑模型的基础是理想的绝热条件,而实际实验过程中:CR53样品从加载到ASB失效过程较短,更接近绝热条件;CG样品具有较好的应变硬化能力,因此从加载至产生ASB失效经历时间更长,考虑到热耗散因素,实际的剪切带宽度小于理论值。

    图  13  绝热剪切带宽度分布
    Figure  13.  ASBs width distributions

    在FeNiAlC合金准静态拉伸变形的研究[15]中,CG样品的主要变形方式为位错滑移和少量孪生变形,硬化能力较弱。而在本文中,由于应变率提高约108倍,面心立方结构一般具有正的率敏感系数,动态条件下会显著提升流变应力水平,极大地促进孪晶的形成[29-30],变形方式以孪晶与位错交互作用为主导,因此剪切韧性得到充分提升。在准静态拉伸下,CR样品变形方式以位错滑移和马氏体相变为主,相变诱导塑性(transformation induced plasticity,TRIP)效应显著,所以保持了较高的均匀塑性。但在动态剪切加载下,由于应变速率的提高,局部温升导致材料层错能升高,马氏体相变被抑制[31],因此应变硬化能力不足导致过早地发生局部失稳。同时,CR样品冷轧后,由于存在明显的织构不易于发生孪生变形[17],所以在动态剪切时,屈服后硬化能力不足,较早地发生失稳。

    ASB的萌生、绝热温升和应力骤降三者的因果关系,一直是材料绝热剪切变形研究的热点问题之一[3, 27]。在热塑模型理论中,认为绝热温升是形成ASB失效的主要原因,最大应力准则也是建立在温升导致的ASB形成后产生应力的骤然下降上。但是,由于动态冲击的变形特点(高应变率和极短时间),ASB范围极窄,形成速度极快与载荷下降时间相差无几,所以一直没有明确的实验证据证实三者的先后关系。Guo等[27-28]试图用实验证实材料剪切带失效与温升的关系,结果表明:在Ti合金中ASB的萌生先于塑性变形产生明显温升,而材料的最大应力也出现在ASB萌生前。该结果与传统的热塑模型所给的预测结论相悖。本文中关于绝热剪切带热塑模型分析计算的结果与实验结果对比显示,当剪切变形时间较长时(材料的硬化能力较强),变形过程与热塑模型的理论前提(理想绝热条件)有所差异,材料的热耗散和微观结构变形不均匀等原因会影响ASB的形成和扩展,使热塑模型的结论与真实实验结果相差较大。同时,本文中也证实了绝热温升不是双相钢材料形成ASB的主导因素或唯一因素,还应考虑共同影响材料ASB失效过程的局部化应变和材料微观结构演化等。

    对FeNiAlC合金进行帽形冲击实验,研究其动态剪切变形行为和微结构机理。实验主要对均匀化处理的粗晶结构样品和冷轧后的层状结构样品进行强迫剪切,研究不同结构样品均匀变形过程和发生ASB失效前后的微观结构变化,利用热塑模型对ASB宽度和塑性功的热转化温升进行分析,主要结论如下。

    (1)FeNiAlC双相钢拥有较优异的动态剪切性能,剪切强度为400~1 300 MPa,剪切应变为1.5~5.1。同等剪切强度的剪切韧性均高于传统的TWIP钢、304、316不锈钢等合金材料。

    (2)孪生与位错主导剪切变形过程中,CG样品在变形过程中发生大量的孪生变形,位错密度显著提高,位错与孪晶界的交互作用使材料获得良好的硬化能力,提高了剪切韧性。CR样品因在剪切变形前冷轧引入大量的位错,变形初期位错密度极高,屈服强度达1.3 GPa,同时因冷轧织构的存在抑制孪生变形;变形过程中,绝热温升使位错发生部分回复,材料的应变硬化能力不及CG样品,屈服后在较小的均匀应变下形成绝热剪切带,发生失稳。动态剪切下,材料的马氏体相变被抑制。

    (3)绝热剪切带内发生动态再结晶,形成晶粒取向随机的超细晶晶粒,同时因剪切带内极高的温度,带内不发生马氏体相变。

    (4)热塑模型在应变较小的情况能够与实验结果高度吻合;绝热温升不是ASB形成的主导因素或唯一因素,ASB的形成可能由局部化变形、材料的微观结构演化等影响。

  • 图  1  3D Voronoi模型及网格划分

    Figure  1.  3D Voronoi model and grid division

    图  2  模型A、B拉伸加载示意

    Figure  2.  Schematic under tensile loading of model A and B

    图  3  动态拉伸加载应变率-时间曲线

    Figure  3.  Stain rate-time curve of dynamic tensile simulations

    图  4  模型A和B拉伸向均分示意图

    Figure  4.  Schematic of drawing direction equalization of model A and B

    图  5  各组模型伪应变能与内能比值(η=Eas/U

    Figure  5.  Ratio of the artificial strain energy to internal energy (η=Eas/U) for dynamic tension simulations

    图  6  300 s−1拉伸时模型z=0截面变形过程

    Figure  6.  Deformations of foams at z=0 section under dynamic tension loading with strain rate of 300 s−1

    图  7  模型动态拉伸破坏示意图

    Figure  7.  Dynamic tensile failure diagram of each model

    图  8  拉伸模型应力-应变曲线及其应力不均匀性指标

    Figure  8.  Stress-strain curves and stress inhomogeneity of models

    图  9  各组模型应变不均匀性指标

    Figure  9.  Strain inhomogeneity of foam models

    图  10  模型B重复性验证

    Figure  10.  Repeatability verification of model B

    图  11  模型B高应变率拉伸破坏

    Figure  11.  Failure of the model B under high strain rate tensile loading

    图  12  模型B动态拉伸应力-应变曲线及不均匀性指标(5000 s−1

    Figure  12.  Dynamic tensile stress-strain curves and non-uniformity of the model B (5000 s−1)

    图  13  模型B不同应变率动态拉伸应力−应变曲线

    Figure  13.  Stress-strain curves of model B under dynamic tension loading with different strain rates

    图  14  应变率对泡沫铝动态拉伸破坏应力和应变的影响

    Figure  14.  Effects of strain rate on failure stress and strain of aluminum foam under dynamic tension loading

  • [1] GIBSON L J, ASHBY M F, ZHANG J, et al. Failure surfaces for cellular materials under multiaxial loads—I. Modelling [J]. International Journal of Mechanical Sciences, 1989, 31(9): 635–663. DOI: 10.1016/S0020-7403(89)80001-3.
    [2] YUAN H, ZHANG J X. Dynamic response of slender multilayer sandwich beams with fiber-metal laminate face-sheets subjected to low-velocity impact [J]. Thin-Walled Structures, 2022, 172: 108932. DOI: 10.1016/j.tws.2022.108932.
    [3] BABAKHANI A, GOLESTANIPOUR M, ZEBARJAD S M. Modelling of aluminium foam core sandwich panels under impact perforation [J]. Materials Science and Technology, 2016, 32(13): 1330–1337. DOI: 10.1080/02670836.2015.1122297.
    [4] 朱凌, 郭开岭, 余同希, 等. 泡沫金属夹芯梁在重复冲击下的动态响应 [J]. 爆炸与冲击, 2021, 41(7): 073101. DOI: 10.11883/bzycj-2020-0198.

    ZHU L, GUO K L, YU T X, et al. Dynamic responses of metal foam sandwich beams to repeated impacts [J]. Explosion and Shock Waves, 2021, 41(7): 073101. DOI: 10.11883/bzycj-2020-0198.
    [5] 张元瑞, 汪高飞, 张永亮, 等. 基于泡沫铝复合结构的汽车座椅横梁填充设计与优化 [J]. 机械强度, 2022, 44(1): 140–147. DOI: 10.16579/j.issn.1001.9669.2022.01.019.

    ZHANG Y R, WANG G F, ZHANG Y L, et al. Filling design and optimization of automobile seat crossbeam based on aluminum foam composite structure [J]. Journal of Mechanical Strength, 2022, 44(1): 140–147. DOI: 10.16579/j.issn.1001.9669.2022.01.019.
    [6] 刘伟明, 程和法, 黄笑梅, 等. 开孔泡沫铝填充圆管的准静态压缩行为 [J]. 爆炸与冲击, 2009, 29(6): 654–658. DOI: 10.11883/1001-1455(2009)06-0654-05.

    LIU W M, CHENG H F, HUANG X M, et al. Quasi-static compression behaviors of cylindrical tubes filled with open-cell aluminum foam [J]. Explosion and Shock Waves, 2009, 29(6): 654–658. DOI: 10.11883/1001-1455(2009)06-0654-05.
    [7] NAMMI S K, MYLER P, EDWARDS G. Finite element analysis of closed-cell aluminium foam under quasi-static loading [J]. Materials & Design, 2010, 31(2): 712–722. DOI: 10.1016/j.matdes.2009.08.010.
    [8] 庄蔚敏, 王恩铭. 胞孔孔径对泡沫铝压缩力学性能影响的仿真研究 [J]. 机械工程学报, 2022, 58(12): 75–82,92. DOI: 10.3901/JME.2022.12.075.

    ZHUANG W M, WANG E M. Simulation study on the influence of pore size on the compression mechanical properties of aluminum foam [J]. Chinese Journal of Mechanical Engineering, 2022, 58(12): 75–82,92. DOI: 10.3901/JME.2022.12.075.
    [9] 王二恒, 虞吉林, 王飞, 等. 泡沫铝材料准静态本构关系的理论和实验研究 [J]. 力学学报, 2004, 36(6): 673–679. DOI: 10.3321/j.issn:0459-1879.2004.06.005.

    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]. Chinese Journal of Theoretical and Applied Mechanics, 2004, 36(6): 673–679. DOI: 10.3321/j.issn:0459-1879.2004.06.005.
    [10] BALCH D K, O’DWYER J G, DAVIS G R, et al. Plasticity and damage in aluminum syntactic foams deformed under dynamic and quasi-static conditions [J]. Materials Science and Engineering: A, 2005, 391(1/2): 408–417. DOI: 10.1016/j.msea.2004.09.012.
    [11] AN Y K, YANG S Y, ZHAO E T, et al. Characterization of metal grid-structure reinforced aluminum foam under quasi-static bending loads [J]. Composite Structures, 2017, 178: 288–296. DOI: 10.1016/j.compstruct.2017.07.031.
    [12] YANG B, TANG L Q, LIU Y P, et al. Localized deformation in aluminium foam during middle speed Hopkinson bar impact tests [J]. Materials Science and Engineering:A, 2013, 560: 734–743. DOI: 10.1016/j.msea.2012.10.027.
    [13] 李忠献, 张茂轩, 师燕超. 闭孔泡沫铝的动态压缩性能试验研究 [J]. 振动与冲击, 2017, 36(5): 1–6. DOI: 10.13465/j.cnki.jvs.2017.05.001.

    LI Z X, ZHANG M X, SHI Y C. Tests for dynamic compressive performance of closed-cell aluminum foams [J]. Journal of Vibration and Shock, 2017, 36(5): 1–6. DOI: 10.13465/j.cnki.jvs.2017.05.001.
    [14] 常白雪, 郑志军, 赵凯, 等. 具有恒定冲击载荷的梯度泡沫金属材料设计 [J]. 爆炸与冲击, 2019, 39(4): 041101. DOI: 10.11883/bzycj-2018-0431.

    CHANG B X, ZHENG Z J, ZHAO K, et al. Design of gradient foam metal materials with a constant impact load [J]. Explosion and Shock Waves, 2019, 39(4): 041101. DOI: 10.11883/bzycj-2018-0431.
    [15] 习会峰, 姚一鸣, 刘逸平, 等. 泡沫金属中高速拉伸的试验数据处理与材料力学性能测量 [J]. 实验力学, 2020, 35(6): 978–984. DOI: 10.7520/1001-4888-19-043.

    XI H F, YAO Y M, LIU Y P, et al. Data processing and mechanical properties measurement of foamed metal in medium to high-speed tensile test [J]. Journal of Experimental Mechanics, 2020, 35(6): 978–984. DOI: 10.7520/1001-4888-19-043.
    [16] HE M C, LI C, GONG W L, et al. Dynamic tests for a constant-resistance-large-deformation bolt using a modified SHTB system [J]. Tunnelling and Underground Space Technology, 2017, 64: 103–116. DOI: 10.1016/j.tust.2016.12.007.
    [17] LEDFORD N, PAUL H, GANZENMÜLLER G, et al. Investigations on specimen design and mounting for Split Hopkinson Tension Bar (SHTB) experiments [J]. EPJ Web of Conferences, 2015, 94: 01049. DOI: 10.1051/epjconf/20159401049.
    [18] 邹广平, 谌赫, 唱忠良. 一种基于SHTB的Ⅱ型动态断裂实验技术 [J]. 力学学报, 2017, 49(1): 117–125. DOI: 10.6052/0459-1879-16-239.

    ZOU G P, CHEN H, CHANG Z L. A modified mode II dynamic fracture test technique based on SHTB [J]. Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(1): 117–125. DOI: 10.6052/0459-1879-16-239.
    [19] MIRONE G, BARBAGALLO R, CADONI E. Tensile test of a HSLA steel at high strain rates with two different SHTB Facilities [J]. Procedia Engineering, 2017, 197: 89–98. DOI: 10.1016/j.proeng.2017.08.085.
    [20] 李尚昆, 胡文军, 徐伟芳, 等. 高温霍普金森拉杆实验技术研究进展 [J]. 中国测试, 2018, 44(10): 35–42. DOI: 10.11857/j.issn.1674-5124.2018.10.006.

    LI S K, HU W J, XU W F, et al. Research progress on SHTB experiment technique at elevated temperature [J]. China Measurement & Test, 2018, 44(10): 35–42. DOI: 10.11857/j.issn.1674-5124.2018.10.006.
    [21] 康建功, 刘芳, 马惠香. 泡沫铝衰减一维应力波特性研究 [J]. 爆破, 2018, 35(1): 180–185. DOI: 10.3963/j.issn.1001-487X.2018.01.029.

    KANG J G, LIU F, MA H X. Study on attenuation characteristic of one-dimensional stress wave in foam aluminum [J]. Blasting, 2018, 35(1): 180–185. DOI: 10.3963/j.issn.1001-487X.2018.01.029.
    [22] XIE B X, TANG L Q, LIU Y P, et al. Numerical analysis on usability of SHPB to characterize dynamic stress–strain relation of metal foam [J]. International Journal of Applied Mechanics, 2017, 9(5): 1750075. DOI: 10.1142/s1758825117500752.
    [23] ZHANG X Y, WU Y D, TANG L Q, et al. Modeling and computing parameters of three-dimensional Voronoi models in nonlinear finite element simulation of closed-cell metallic foams [J]. Mechanics of Advanced Materials and Structures, 2018, 25(15/16): 1265–1275. DOI: 10.1080/15376494.2016.1190426.
    [24] WU Y D, QIAO D, TANG L Q, et al. Global topology of failure surfaces of metallic foams in principal-stress space and principal-strain space studied by numerical simulations [J]. International Journal of Mechanical Sciences, 2019, 151: 551–562. DOI: 10.1016/j.ijmecsci.2018.12.003.
    [25] ZHANG X Y, TANG L Q, LIU Z J, et al. Yield properties of closed-cell aluminum foam under triaxial loadings by a 3D Voronoi model [J]. Mechanics of Materials, 2017, 104: 73–84. DOI: 10.1016/j.mechmat.2016.10.007.
    [26] ZHANG X Y, TANG L Q, JIANG Z Y, et al. Effects of Meso shape irregularity of metal foam on yield features under triaxial loading [J]. International Journal of Structural Stability and Dynamics, 2015, 15(7): 1540014. DOI: 10.1142/S0219455415400143.
  • 加载中
图(14)
计量
  • 文章访问数:  213
  • HTML全文浏览量:  71
  • PDF下载量:  92
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-04-10
  • 修回日期:  2023-10-25
  • 网络出版日期:  2023-11-07
  • 刊出日期:  2024-01-11

目录

/

返回文章
返回