Constitutive Model of Aluminum Foams Considering Temperature Effect under Quasi-Static Compression
-
摘要: 采用MTS材料试验机研究了不同密度(0.322~0.726g/cm3)的闭孔泡沫铝在温度范围25~500℃下的准静态压缩力学性能,得到了泡沫铝在不同温度下的单轴压缩应力-应变曲线,分析了密度以及温度对其力学行为的影响。利用Liu和Subhash提出的本构模型对不同密度泡沫铝的应力-应变曲线进行拟合,分析并确定了模型中各参数随密度变化的函数,再代入Liu-Subhash模型,得到了泡沫铝的准静态压缩本构模型。通过引入温度软化项对准静态压缩本构模型进行修正,建立了考虑温度效应的泡沫铝准静态压缩本构模型,对闭孔泡沫铝的工程应用具有指导意义。Abstract: In this study, we employed the MTS machine to apply a compressive load to the closed cell aluminum foams with different densities (0.322-0.726g/cm3) at different temperatures (25-500℃) under quasi-static state, and obtained the stress-strain curves of the aluminum foam under uniaxial compression and analyzed the influence of density and temperature on its mechanical behavior.We fitted the stress-strain curves at different densities using the Liu and Subhash constitutive model, and the fitted results accorded well with the experimental curves.Then we analyzed and determined the function between the revised 5 parameters and the varying density, and established a constitutive model of the aluminum foams with density effect taken into consideration.With the temperature softening effect added to the modification of this constitutive model, we eventually established a constitutive model of aluminum foams considering both the temperature effect and the density effect under quasi-static compression.
-
泡沫金属质轻,具有较高的比刚度和比强度,以及隔热、电磁屏蔽等物理性能[1],在汽车交通、铁路、航空航天等领域广泛应用,如用于飞机外壳夹层、汽车防冲档。除此之外,泡沫金属在受到压缩时,由于其应变滞后于应力,压缩应力-应变曲线中有一个很长的低应力平台,可承受较大的塑性变形,因此泡沫金属具有良好的吸能特性,可用于缓和冲击的工程构件、能量吸收和防振构件[2]。
目前,常温下泡沫金属的静/动态本构关系已经得到广泛的研究。Chen和Lu[3]提出了一个依赖于特征应力和总应变的应力势,在此基础上建立了一个唯象的可压缩弹塑性本构模型的框架。该模型避免了人为区分应力-应变的弹塑性区带来的影响。王二恒等[4]利用Chen和Lu[3]提出的唯象本构模型框架,建立了一个泡沫金属准静态本构模型,得到了泡沫金属在三维等比例压缩和侧向受约束轴向压缩时的宏观应力-应变曲线。王志华等[5]提出了一个多参数的非线性弹塑性唯象本构模型,可以全面地描述泡沫金属材料线弹性段、应力平台段和密实段的典型三阶段变形特征。
泡沫金属是一种典型的低成本轻质材料,其复合结构有望作为近空间飞行器中的重要结构部件。但是,近空间飞行器苛刻的应用环境不但要求其结构轻质化,而且要求泡沫金属在高温下有较好的承载、隔热和冲击吸能能力。目前,考虑温度效应的泡沫金属力学行为的研究还比较少。Hakamada等[6]开展了ALPORAS闭孔泡沫铝及其基体材料在温度范围573~ 773 K内的准静态压缩试验,发现闭孔泡沫铝在高温下的变形机制与其基体材料的变形机制本质上是相同的。Aly[7]开展了ALPORAS闭孔泡沫铝在常温和高温下的压缩实验,研究了相对密度和实验温度的影响,研究发现,胞壁屈曲是闭孔泡沫铝的主要变形机制,增大密度与升高实验温度对泡沫金属力学性能的影响刚好相反。Cady等[8]研究了ALPORAS泡沫铝在不同应变率(0.001~1 800 s-1)和不同温度条件(77~295 K)下的力学性能,结果表明,闭孔泡沫铝的力学性能对温度具有很强的依赖性。
考虑泡沫金属温度效应的本构关系相关研究也比较匮乏。习会峰等[9]开展了-50~300 ℃范围内不同温度下泡沫铝的静态压缩实验,建立了考虑温度效应的泡沫铝静态压缩本构模型。王鹏飞等[10]基于Sherwood和Frost提出的本构关系框架,分析了泡沫铝本构方程中温度效应与应变率效应的耦合关系,对本构方程中的应变率敏感系数进行了适当修正,修正后的本构关系与实验结果的重合性较好,在此基础上得到了泡沫铝在一定密度范围内包含温度、应变率影响的较为完备的本构方程。
本研究通过开展不同密度(0.322~0.726 g/cm3)的泡沫铝在不同温度(25~500 ℃)下的准静态压缩实验,分析泡沫铝在不同温度下的力学性能,测得不同温度下泡沫铝的单轴压缩应力-应变曲线。基于Liu-Subhash模型[11],对不同温度下的实验数据进行拟合,建立考虑温度效应和密度效应的泡沫铝准静态压缩本构模型。
1. 实验方法及方案
实验材料为闭孔泡沫铝,由上海奥深特金属复合材料科技有限公司提供。本研究采用的闭孔泡沫铝材料原尺寸为400 mm×400 mm×600 mm,孔径介于2~5 mm之间。准静态压缩实验采用线切割技术加工圆柱形试件。为了保证均匀性,排除胞孔尺寸影响,试件高度和直径不小于5个胞孔,试件尺寸为Ø32 mm×15 mm。实验前,对每个试样进行称量,得到其真实的相对密度,测得试样密度范围为0.322~0.796 g/cm3。
准静态实验在MTS810材料试验机上进行,加载速率为0.001 s-1。实验温度分别取25、200、275、350、425和500 ℃,高温实验在MTS810自带的高温箱中进行,温度误差为±5 ℃。
2. Liu-Subhash本构模型
为了准确描述泡沫材料应力-应变曲线的三阶段特征,Liu与Subhash[11]提出了一维六参数本构模型
σ=p1ep2ε−1p6+ep3ε+ep4(ep5ε−1) (1) 式中:σ为压缩应力,ε为压缩应变,p1、p2、p3、p4、p5和p6为用于描述泡沫材料性质的参数。p1主要描述压缩时屈服应力的变化;p2和p3为无量纲参数,主要描述应力平台阶段的硬化性能;p4和p5主要描述密实阶段的起点和密实阶段斜率的大小;p6主要描述弹性段斜率的大小。从模型中可以看出,参数p6的作用可以用p1和p3代替。习会峰等[9]提出简化模型,用常数1取代参数p6,将其简化为五参数模型
σ=p1ep2ε−11+ep3ε+ep4(ep5ε−1) (2) 3. 实验结果与分析
3.1 包含密度效应的泡沫铝本构模型
3.1.1 密度对应力-应变曲线的影响
图 1为不同密度ρ下泡沫铝的准静态压缩应力-应变实验曲线。从实验结果来看,泡沫铝的应力-应变曲线的形状及发展趋势相似,呈典型的三阶段特性,即:较小应变下的线弹性阶段、对应于胞壁塑性坍塌而缓缓上升的应力平台阶段和密实后的应力迅速上升阶段。从图 1中可以看出,泡沫铝材料的初始密度对其应力-应变曲线影响巨大,随着密度的增大,曲线依次抬高,屈服强度也相应增大,但压缩至致密段的最大应变量减小。因此,在基体材料相同的情况下,在一定的密度范围内,泡沫铝材料的力学性能主要由材料的初始密度决定,初始密度越大,材料承载能力越强。
3.1.2 基于Liu-Subhash模型的实验数据拟合
基于简化的Liu-Subhash模型,对实验得到的常温下不同密度泡沫铝的准静态应力-应变曲线进行拟合,得到不同密度泡沫铝材料的模型参数值,如表 1所示。实验曲线和拟合曲线如图 2所示,可见拟合效果较好,可采用该模型描述泡沫铝的单轴准静态压缩应力-应变关系。
表 1 不同密度下的模型参数值Table 1. Parameter values for different densitiesDensity/(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 3.1.3 考虑密度效应的修正模型
由表 1可知,参数p1、p2、p3、p4、p5随密度的变化而变化,显然这5个参数都是密度ρ的函数。因此,(2)式可以写成
σ=p1(ρ)ep2(ρ)ε−11+ep3(ρ)ε+ep4(ρ)[ep5(ρ)ε−1] (3) (3) 式即为考虑了密度影响的泡沫铝的常温准静态压缩本构模型。图 3为模型中5个参数随密度的变化情况,对实验数据进行拟合,得到5个参数与密度的关系为
{p1(ρ)=20.34ρ1.69p2(ρ)=94.58ρ0.08p3(ρ)=94.41ρ0.08p4(ρ)=12.34ρ−10.55p5(ρ)=−8.6ρ+14.68 (4) 由(3)式和(4)式得到考虑密度效应的闭孔泡沫铝静态压缩本构模型,该模型可用于描述不同密度下泡沫铝的应力-应变曲线,具体形式如下
σ=20.34ρ1.69e94.58ρ0.08ε−11+e94.41ρ0.08ε+e12.34ρ−10.55[e(−8.6ρ+14.68)ε−1] (5) 3.2 包含温度效应的泡沫铝本构模型
3.2.1 温度对应力-应变曲线的影响
图 4是密度为0.322 g/cm3的泡沫铝在不同温度下的准静态应力-应变曲线。由图 4可知,随着温度的升高,材料由硬变软,由脆变韧。材料呈现明显的温度软化效应,当温度从25 ℃上升到500 ℃时,泡沫铝屈服强度从4.70 MPa下降到0.88 MPa。
3.2.2 考虑温度效应的修正模型
为了得到包含温度效应的泡沫铝本构模型,引入温度软化项
H(T)=1−T∗m (6) 式中:T*为无量纲温度项,m为指数。
T∗=T−TroomTmelt−Troom (7) 式中:Troom为室温,Troom=298 K;Tmelt为铝合金的熔点,Tmelt=933 K。
因此,(3)式可进一步写成
σ={p1(ρ)ep2(ρ)ε−11+ep3(ρ)ε+ep4(ρ)[ep5(ρ)ε−1]}⋅[1−(T−TroomTmelt−Troom)m] (8) (8) 式即为考虑温度效应和密度效应的泡沫铝准静态本构模型。基于该模型,对不同温度下得到的准静态压缩应力-应变曲线进行拟合,得到参数m=1.31。不同温度下泡沫铝的应力-应变实验及拟合曲线如图 5所示。从图 5中可以看到,应力平台阶段和压实段拟合曲线与实验曲线吻合较好,说明该模型可用于描述不同温度下泡沫铝的应力-应变曲线。但在弹性段,拟合曲线与实验结果吻合得不太好,主要原因是Liu-Subhash模型不能很好地描述应力-应变曲线中的应力降现象,另一方面是由于实验时温度控制不精确产生的误差。总体上拟合曲线可以反映实验曲线的特征,拟合得到的参数值是可信的。
4. 结论
利用MTS万能材料试验机研究了不同密度(0.322~0.726 g/cm3)的闭孔泡沫铝在不同温度(25~500 ℃)下的静态压缩力学性能,实验结果表明:泡沫铝材料的初始密度对其应力-应变曲线影响巨大,随着密度的增大,屈服强度也相应增大,材料承载能力增强;随着温度的升高,泡沫铝材料力学特性由硬变软,呈现明显的温度软化效应。利用Liu-Subhash模型对不同密度下的实验数据进行拟合,拟合效果很好,分析并确定了模型中5个参数随密度变化的函数,并代入Liu-Subhash模型,得到了考虑密度效应的本构模型;又引入温度软化项对本构模型进行修正,建立了综合考虑温度效应和密度效应的泡沫铝准静态压缩本构模型。
-
表 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 -
[1] 刘培生, 李铁藩, 傅超.多孔金属材料的应用[J].功能材料, 2001, 32(1):12-15. http://www.cqvip.com/QK/72003x/201607/epub1000000248429.htmlLIU 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.htmlSHANG 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.htmWANG 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=hebgcdxxb201308010XI 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.004WANG 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 -