Loading [MathJax]/jax/output/HTML-CSS/jax.js
  • ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

泡沫铝率相关性能的有限元模拟

范志庚 陈常青 万强

张凤国, 周洪强, 胡晓棉, 王裴, 邵建立, 冯其京. 强动载作用下孔洞汇合对延性金属层裂损伤演化过程的影响[J]. 爆炸与冲击, 2016, 36(5): 596-602. doi: 10.11883/1001-1455(2016)05-0596-07
引用本文: 范志庚, 陈常青, 万强. 泡沫铝率相关性能的有限元模拟[J]. 爆炸与冲击, 2014, 34(6): 742-747. doi: 10.11883/1001-1455(2014)06-0742-06
Zhang Fengguo, Zhou Hongqiang, Hu Xiaomian, Wang Pei, Shao Jianli, Feng Qijing. Influence of void coalescence on spall evolution of ductile polycrystalline metal under dynamic loading[J]. Explosion And Shock Waves, 2016, 36(5): 596-602. doi: 10.11883/1001-1455(2016)05-0596-07
Citation: Fan Zhi-geng, Chen Chang-qing, Wan Qiang. Finite element simulation on the rate-dependent properties of aluminum foams[J]. Explosion And Shock Waves, 2014, 34(6): 742-747. doi: 10.11883/1001-1455(2014)06-0742-06

泡沫铝率相关性能的有限元模拟

doi: 10.11883/1001-1455(2014)06-0742-06
基金项目: 国家自然科学基金项目(11102196, 11372295);国家重点基础研究发展计划(973计划)项目(2010CB832700);中国工程物理研究院科学技术发展基金项目(2014A0203006)
详细信息
    作者简介:

    范志庚(1978—), 男, 博士, 高级工程师

  • 中图分类号: O344.7

Finite element simulation on the rate-dependent properties of aluminum foams

Funds: Supported by the National Natural Science Foundation of China (11102196, 11372295); the National Basic Research Program of China (973 Program) (2010CB832700)
More Information
  • 摘要: 构建了三维随机分布球形泡孔模型,模拟开、闭孔混合结构泡沫铝材料的微细观结构,并通过有限元方法计算了10~104 s-1应变率范围内、孔隙率35%~65%泡沫铝材料的率相关性以及应变率和相对密度变化对泡沫铝动态压缩力学性能的影响。研究表明:中、低应变率下,泡沫铝材料率相关性能主要取决于基体材料的应变率敏感性;高应变率下,泡沫铝材料率相关性能受基体材料的应变率敏感性以及微结构惯性联合作用,且相对密度较低泡沫铝材料的微结构惯性效应更显著。
  • 在不同厚度飞片相互撞击过程中,冲击波从材料自由表面反射,材料内部或近自由表面区域受到卸载稀疏波作用,如果延性金属材料内部的拉伸应力足够高,以孔洞形式的损伤将经历成核、增长和汇合演化过程,当微孔洞的数量和尺寸增长到一定阶段之后,它们的(强)相互作用就变得十分重要,这时,微孔洞之间的汇合以及细观尺度上的不稳定性将导致材料的最终破坏。实验结果分析表明:在强动态拉伸载荷作用下,孔洞汇合可能开始于孔洞成核后,并伴随损伤演化整个过程,此时,孔洞汇合对孔洞增长起到抑制作用[1];孔洞增长不仅以孤立孔洞形式增长,而且可能主要以孔洞汇合形式增长[2]。因此,在连续损伤模型中耦合孔洞汇合的影响,有助于精细描述延性材料损伤演化过程、分析不同因素对损伤演化过程的影响[3-5]

    由于微孔洞之间强相互作用在理论分析方面遇到的困难,至今还没有一个较好的模型对孔洞汇合进行合理的描述[6]。一个合理的孔洞汇合模型不仅应给出开始汇合的准则,而且还应该对其间的应力松弛或材料软化予以适当的描述。当然要作以上分析是非常困难的,因而许多学者通过数值模拟两个孔洞之间的相互关系给出一些半经验性的孔洞汇合临界条件,即以孔洞之间的相互距离或临界损伤作为孔洞汇合的临界判据:T.Pardoen等[7]认为孔洞汇合的临界损伤依赖于微孔洞初始体积分数、加载应力三轴度及微孔洞形状;M.F.Horstemeyer等[8]采用微力学有限元计算方法研究准静态加载条件下两个孔洞之间的贯通行为,指出孔洞开始汇合的临界韧带距离依赖于边界条件和材料韧性,且约为2~8倍的孔洞直径;E.T.Seppala等[9]采用分子动力学方法对动态条件下铜材料中孔洞的增长和聚集行为进行了数值模拟,结果显示临界韧带距离约为0.5倍孔洞直径时,孔洞周围塑性区相遇,孔洞增长速度突然增加,从而导致孔洞汇合。此外,在理论分析方面,P.F.Thomson[10]简化了塑性极限载荷的滑移线场解, 提出了一种孔洞汇合的判断准则;D.L.Tonks等[11]进一步用高应变率下的随机滤渗理论解释了断裂点, 提出了一种通过孔洞聚集形成小尺度孔洞簇增长的模型;T.Pardoen等[12]讨论了基于微结构机理的孔洞汇合模型的进展,并指出因微孔洞汇合引起金属断裂过程最后阶段的模拟需要反映微结构演化重要信息的延性损伤模型。在过去的十多年里,孔洞汇合研究的重要进步是基于简单的孔隙度准则或临界应变准则发展为基于微观机理认识的进展。不过,目前的研究主要是分析相同大小孔洞间的汇合,还需要将其扩展到更一般的情况,即不同大小的孔洞之间的汇合问题;此外,对于孔洞汇合最初发生于相同的大孔洞之间、小孔洞之间亦或不同大小孔洞之间,目前还没有确切的说明。同时,现有的孔洞汇合模型严格来说是描述孔洞变形对损伤的影响,且一般采用指数函数来唯像描述孔洞汇合后损伤的快速增长[13],而根据损伤度的定义,孔洞汇合时刻并没有引起损伤的改变,但实际却促进了损伤的快速发展,目前对于其物理机理至今仍处于探讨之中。

    本文中基于两个不同大小孔洞之间的几何关联,给出孔洞汇合的临界判定方法,并基于能量守恒原理,揭示孔洞汇合后引起损伤增长的物理机理,同时,耦合孔洞汇合的影响,采用数值方法讨论孔洞汇合对延性金属层裂损伤演化特性的影响。

    现阶段对于孔洞汇合机理的分析已经成为层裂损伤研究的主要关注点[14-16],不过,现有的汇合判据还存在不足之处:距离判据没有考虑损伤的影响;损伤(或应力、应变)判据没有考虑孔洞大小的影响。此外,在强动态加载情况下,孔洞增长过程不易发生塑性局域化,孔洞几乎仍然保持球形形状,孔洞汇合时,损伤度较高,且孔洞之间产生相互接触,这可能是因为惯性对孔洞间的颈缩起了阻碍作用[14],这种情况在一些层裂实验[17-18]以及裂纹扩展实验[19]中均可以观测到(见图 1)。因汇合前孔洞基本保持球形形状,且相邻孔洞间距较小,则因孔洞间基体材料颈缩引起孔洞变化对损伤的影响较小,因此相对于基于颈缩过程分析得到的孔洞汇合判据,采用孔洞间距离判据更加简单、适宜。

    图  1  孔洞汇合实验结果
    Figure  1.  Void coalescence by direct impingement in tantalum

    为此,我们考虑两个空心球壳之间的几何关联(见图 2),并定义孔隙度:

    α=b3/(b3a3)
    图  2  孔洞间的几何关系
    Figure  2.  Porous material model

    则孔洞间距离可以表示为:

    d=[(αα1)1/31](a1+a2) (1)

    本文中仍将采用孔洞间的距离作为孔洞开始汇合的判据,并根据实验结果确定临界孔洞距离dcr(根据文献[2]的实验结果,本文中采用dcr=4min(a1, a2)),即当两个孔洞间的距离ddcr时,认为相应的两个孔洞之间发生汇合。与现有的距离计算方法不同,公式(1)不仅包含了孔洞大小,同时也包含了损伤度的影响。众所周知,在相同加载条件下小孔洞周边的应力集中更明显,同时,在大孔洞不变的情况下,小孔洞越小,其对应的损伤度越小,因此,孔洞汇合应该首先发生于最大孔洞与最小孔洞之间,这也从侧面说明了实验最后的观测结果以及相关的数值模拟结果中小孔洞很少的可能原因。

    a1=ka2d=ma1=mka2,则公式(1)可以转化为损伤度D的表达式:

    D=α1α=[1+k1+(m+1)k]3 (2)

    两个特例:

    m=k=1时,D=0.296,这与分析层裂问题时常采用的临界断裂损伤度相近;

    m=1时,有D=α1α=limk[1+k1+2k]3=0.125

    这与L.M.Brown等[20]和D.L.Tonks等[11]分析的孔洞汇合初始临界损伤度完全一致,换句话说,文献[11, 20]中所采用的孔洞汇合临界损伤度只是我们所讨论的特例,图 3显示了孔洞汇合损伤度与汇合孔洞相对大小比值之间的对应关系,同时,计算结果也显示了孔洞汇合损伤度的取值范围在0.296~0.125之间。

    图  3  孔洞相对大小与孔洞汇合临界损伤度关系
    Figure  3.  Critical damage for void coalescence vs. relative difference in size between two voids

    k=10时,图 4显示了孔洞汇合损伤度与孔洞距离之间的对应关系。综合分析图 3~4和公式(2)可知,孔洞汇合是损伤、孔洞相对大小以及孔洞间距离综合影响的结果。对于材料性质、应力状态对孔洞汇合的影响:一方面,损伤与应力、应变耦合在一起[21],损伤对孔洞汇合的影响也间接地反映了材料性质和应力状态对孔洞汇合的影响;另一方面,针对强加载情况,直至汇合前孔洞仍基本保持球形形状,因此,材料的性质和应力状态的影响主要反映在汇合前的孔洞增长过程。

    图  4  孔洞间距离与孔洞汇合临界损伤度关系
    Figure  4.  Critical damage for void coalescence vs. distance between two voids

    现有的孔洞汇合判据(模型),包括距离判据或应力、应变判据(如Tonks模型和Thomason模型),都是基于孔洞间材料颈缩过程得到的,描述的是汇合前孔洞的变化情况,而没有涉及孔洞汇合后对损伤的影响。对于在强动态拉伸载荷作用下,因汇合前孔洞基本保持球形形状,且孔洞间距很小,则孔洞间材料颈缩过程对损伤的影响减弱。目前已经有大量的文献研究孔洞间材料颈缩过程对损伤的影响,因此本文中忽略颈缩的影响,主要讨论孔洞汇合后对损伤度增长的影响。

    孔洞汇合遵循能量守恒原则,即在孔洞汇合前后有:

    Ei+Ek=ˉEi+ˉEk (3)

    式中:EiEk分别表示孔洞周围基体材料的内能和动能,EiEk为孔洞汇合后对应的值。

    对于内能的变化,在绝热条件下有:

    ˙e=pvt (4)

    式中:epv分别表示单位体积的内能、压力和材料的相对比容。因基体材料不可压假设,并且假设孔洞汇合并没有改变材料内部的孔隙度,因此有:

    vt=0

    即孔洞汇合时内能守恒。

    此外,针对单一的空心球壳,设r0为基体材料内部的Lagrange坐标,相应的r为Euler坐标,a0α0为初始孔洞半径和初始孔隙度,则根据不可压缩假设有[21]

    r3=r30a30α0αα01 (5)
    ˙r=a303r2(α01)˙α (6)

    则其动能可以表示为:

    Ek=12ba4πρr2˙r2dr=2πρ9a60(α01)2[1(α1α)1/3]˙α21a (7)

    损伤材料内部含有大量不同大小的孔洞,不妨将其等效为不同大小的空心球壳,且每一个空心球壳的内外径之比相同,即αα0相同,这样不仅方便损伤的计算,而且可以间接地考虑孔洞之间的相互作用[22-23]。因孔洞汇合前后αα0和基体材料的密度ρ不变,则根据动能守恒原则,由公式(7)有:

    ˙ˉα=˙α(ni/ai)(Nj/Rj) (8)

    式中:˙αai以及ni为孔洞汇合前的孔隙度增长率、孔洞半径以及对应的孔洞数,˙ˉαRj以及Nj为孔洞汇合后的值。式(8)显示了因孔洞汇合引起的孔洞数减少、孔洞尺寸增加,从而造成孔隙度增长速度的提高,这也明确了因孔洞汇合引起损伤增长的物理机理。

    对于延性金属层裂损伤的研究,我们已经建立了反映初始损伤及孔洞大小、惯性和材料弹塑性效应的层裂损伤物理统计描述方法[22],同时,基于材料晶粒尺寸与潜在孔洞成核数之间的关系,构建了一个耦合晶粒尺寸影响的孔洞成核方程[24]。在前期工作的基础上,我们将孔洞汇合的影响引入到已有的层裂损伤模型中,即在计算损伤演化过程中,根据公式(1),当孔洞间距离ddcr时,采用公式(8)调整孔隙度的增长率,同时结合相关文献对实验结果的分析[2],探讨孔洞汇合对层裂损伤演化过程的影响。层裂实验靶材料分别选用30、60、100和200 μm等4种晶粒尺寸的高纯铜,靶厚均为4 mm;飞片材料为石英(Z-cut quartz),厚度约2 mm,飞片速度约131 m/s。对4种靶材料的实验结果进行了数值模拟,同时,在30 μm晶粒尺寸高纯铜的计算模型中考虑了孔洞汇合的影响。计算中所采用的模型参数和材料参数均与我们前期的工作[22, 24]所采用的参数相同。

    图 5显示了汇合对孔洞尺寸d增长的影响:因汇合以及汇合导致孔隙度增长率的提高,相对于没有考虑汇合的计算结果,孔洞尺寸迅速增长。图 6的计算结果显示:虽然孔洞汇合时损伤度没有增加,但因为汇合引起孔隙度增长率的提高,从而导致了损伤的快速增长。

    图  5  孔洞汇合对孔洞增长的影响
    Figure  5.  Influences of void coalescence on void size
    图  6  孔洞汇合对损伤发展的影响
    Figure  6.  Influences of void coalescence on spall damage

    图 7左图显示了晶粒尺寸、孔洞汇合对自由面速度vf的影响:随着材料平均晶粒尺寸的增加,自由面速度曲线的回跳点降低,即层裂强度增加,这与P.B.Trivedi等[25]对不同晶粒尺寸高纯铝层裂实验结果的分析在定性上一致;同时,对于平均晶粒尺寸为30 μm的材料,考虑孔洞汇合的影响,曲线回跳后的斜率显著增加,上升的幅度增高,图 7右图为对应的自由面速度曲线实验结果:平均晶粒尺寸dg=30、200 μm的损伤材料中存在大量的孔洞汇合现象,相对于没有发生孔洞汇合的自由面速度曲线(60、100 μm),自由面速度曲线回跳后的斜率和上升的幅度增加[2],关于孔洞汇合的计算结果和实验结果定性上符合也较好。

    图  7  晶粒尺寸、孔洞汇合对自由面速度曲线的影响
    Figure  7.  Influences of void coalescence and grain size on free surface velocities

    表 1列出了不同平均晶粒尺寸dg损伤材料内部层裂面处孔洞数N和平均孔洞直径dv的实验统计结果和数值计算结果(注:实验给出的是可观测到的孔洞数,而计算给出的是单位体积(cm3)内的孔洞数)。实验观测结果显示:随着晶粒尺寸的增加,损伤材料内部的孔洞数减少、平均孔洞尺寸增加。计算结果与实验定性符合。同时,对实验观测结果和计算结果的分析表明(30 μm样品):孔洞汇合引起孔洞数减少、平均孔洞尺寸增加。

    表  1  损伤材料内部孔洞数及孔洞大小的统计结果
    Table  1.  Damage statistics
    dg /μm N dv
    实验 计算 实验 计算
    30 236 11.460 38.1 12.17
    0.044(考虑汇合) 37.12(考虑汇合)
    60 363 3.236 22.7 22.36
    100 267 1.421 33.0 34.51
    200 111 0.566 55.1 42.60
    下载: 导出CSV 
    | 显示表格

    图 8显示了选取不同孔洞汇合临界损伤度D对自由面速度曲线的影响:孔洞汇合发生越早,材料内部损伤发展越快,即自由面曲线回跳后曲线上升的斜率和幅度越大。

    图  8  孔洞汇合临界损伤度对自由面速度曲线的影响
    Figure  8.  Influences of critical damage for void coalescence on free surface velocities

    耦合孔洞汇合的影响是精细化描述延性金属材料层裂损伤演化过程的关键。针对强动态拉伸载荷作用下延性金属层裂损伤问题,尝试建立了一个反映材料损伤和材料内部孔洞之间几何信息的孔洞汇合判定方法,从而弥补了现有判据只考虑单一影响因素的不足。同时,基于孔洞汇合前后能量守恒原理,给出了孔洞汇合对损伤增长影响的关系式,明确了孔洞汇合引起损伤增长的物理机理。

    将材料平均晶粒尺寸影响和孔洞汇合影响引入到层裂损伤模型中,结合相关文献的实验分析结果,数值计算分析结果显示:晶粒尺寸越小,损伤材料内部成核孔洞越多、平均孔洞尺寸越小,自由面速度回跳点增高(即层裂强度降低)、回跳后速度曲线上升的斜率降低;孔洞汇合引起回跳后速度曲线上升的斜率增加、损伤材料内部的孔洞数减少、平均孔洞尺寸增加。计算结果与实验结果定性上符合较好,从而在一定程度上推进了延性金属层裂损伤的微细观物理建模研究。

  • 图  1  泡沫材料三维随机有限元模型

    Figure  1.  Three-dimensional finite element models of aluminum foams

    图  2  泡沫模型初始弹性模量随相对密度的分布

    Figure  2.  Initial elastic modulus of foams versus foam relative densities

    图  3  基体材料率相关泡沫模型动态压缩应力应变曲线

    Figure  3.  Dynamic compressive stress-strain curves of foam models with rate-dependent cell material

    图  4  基体材料率无关泡沫模型动态压缩应力应变曲线

    Figure  4.  Dynamic compressive stress-strain curves of foam models with rate-independent cell material

    图  5  动态压缩下泡沫铝流动应力与相对密度关系

    Figure  5.  Flow stresses of aluminum foams versus foam relative densities under dynamic compressive loadings

    表  1  拟合结果

    Table  1.   Fitting results of parameters

    ˙ε/s-1an
    100.6981.72
    1000.7841.72
    1 0000.9421.72
    2 5001.031.68
    5 0001.131.66
    10 0001.271.64
    下载: 导出CSV

    表  2  泡沫模型及其基体材料强度的变化

    Table  2.   Variation of flow stress for foams and its cell material

    ˙ε2/s-1αf/%αs/%
    ρ=0.35ρ=0.50ρ=0.65
    5 00012.0411.2410.668.34
    10 00028.5926.7725.4418.25
    下载: 导出CSV
  • [1] 卢天健, 何德坪, 陈常青, 等.超轻多孔金属材料的多功能特性及应用[J].力学进展, 2006, 36(4): 517-535.

    Lu Tian-jian, He De-ping, Chen Chang-qing, et al. The multi-functionality of ultra-light porous metals and their applications[J]. Advances in Mechanics, 2006, 36(4): 517-535.
    [2] 王永胜, 左孝青, 尹志其, 等.应变率对泡沫铝压缩性能的影响[J].材料导报, 2009, 23(3): 47-52.

    Wang Yong-sheng, Zuo Xiao-qing, Yin Zhi-qi, et al. Effect of strain rate on compressive property of aluminum foam[J]. Materials Review, 2009, 23(3): 47-52.
    [3] Deshpande V S, Fleck N A. High strain rate compressive behavior of aluminium alloy foams[J]. International Journal of Impact Engineering, 2000, 24(3): 277-298.
    [4] 郭伟国, 李玉龙, 黄福增.不同应变率下泡沫铝的形变和力学性能[J].爆炸与冲击, 2008, 28(4): 289-292.

    Guo Wei-guo, Li Yu-long, Huang Fu-zeng. Deformation and mechanical property of aluminium foam at different strain rates[J]. Explosion and Shock Waves, 2008, 28(4): 289-292.
    [5] Dannemann K A, Lankford Jr J. High strain rate compression of closed-cell aluminum foams[J]. Materials Science and Engineering: A, 2000, 239(1/2): 157-164.
    [6] Mukai T, Kanahashi H, Miyoshi T, et al. Experimental study of energy absorption in closed-cell aluminum foam under dynamic loading[J]. Scripta Materialia, 1999, 40(8): 921-927.
    [7] Mukai T, Miyoshi T, Nakano S, et al. Compressive response of a closed-cell aluminum foam at high st rain rate[J]. Scripta Materialia, 2006, 54(4): 533-537.
    [8] 程和法.应变率对泡沫Al-Mg合金压缩性能的影响[J].特种铸造及有色合金, 2003, 5: 1-3.

    Cheng He-fa. The effect of strain rate on the compressive properties of foamed Al-Mg alloy[J]. Special Casting & Nonferrous Alloys, 2003, 5: 1-3.
    [9] 程和法, 黄笑梅, 王强, 等.通孔泡沫铝的动态压缩行为[J].爆炸与冲击, 2006, 26(2): 169-173.

    Cheng He-fa, Huang Xiao-mei, Wang Qiang, et al. The dynamic compressive behaviors of an open-cell aluminum foam[J]. Explosion and Shock Waves, 2006, 26(2): 169-173.
    [10] 田杰, 胡时胜.基体性能对泡沫铝力学行为的影响[J].工程力学, 2006, 23(8): 168-176.

    Tian Jie, Hu Shi-sheng. Effect of matrix properties of the mechanical behaviors of aluminum foams[J]. Engineering Mechanics, 2006, 23(8): 168-176.
    [11] 张健, 赵桂平, 卢天健.闭孔泡沫铝应变率效应的试验和有限元分析[J].西安交通大学学报, 2010, 44(5): 97-101.

    Zhang Jian, Zhao Gui-ping, Lu Tian-jian. Experimental and numerical study on strain rate effects of close-celled aluminum foams[J]. Jounal of Xi'an Jiaotong University, 2010, 44(5): 97-101.
    [12] 胡时胜, 王悟, 潘艺, 等.泡沫材料的应变率效应[J].爆炸与冲击, 2003, 23(1): 13-18.

    Hu Shi-sheng, Wang Wu, Pan Yi, et al. Strain rate effect on the properties foam materials[J]. Explosion and Shock Waves, 2003, 23(1): 13-18.
    [13] Gibson L J, Ashby M F. Cellular solids: Structure and properties[M]. 2nd ed. Cambridge: Cambridge University Press, 1997.
  • 期刊类型引用(1)

    1. 张凤国,刘军,何安民,赵福祺,王裴. 强冲击加载下延性金属卸载熔化损伤/破碎问题的物理建模及其应用. 物理学报. 2022(24): 284-292 . 百度学术

    其他类型引用(1)

  • 加载中
图(5) / 表(2)
计量
  • 文章访问数:  3462
  • HTML全文浏览量:  359
  • PDF下载量:  491
  • 被引次数: 2
出版历程
  • 收稿日期:  2013-04-11
  • 修回日期:  2013-06-26
  • 刊出日期:  2014-11-25

目录

/

返回文章
返回