针对铝合金材料,确定不同缺陷特征参数对铝合金疲劳寿命的影响,通过Pearson系数选取与疲劳损伤影响较大的特征参数进行研究,基于Abaqus分析软件建立了缺陷特征的疲劳损伤模型,在特征权重公式的基础上,获得了不同缺陷特征对疲劳性能影响的相对权重,最后达到量化不同缺陷特征对疲劳性能的影响。 For aluminum alloy, the influence of different defect characteristic parameters on the fatigue life of aluminum alloy was determined. Pearson coefficient was selected and the characteristic parameters with great influence on fatigue damage were studied. The damage model of defect characteristics was established based on Abaqus analysis software. The relative weights of the effects of different defect characteristics on the fatigue properties were obtained. Finally, the mesoscopic damage factors were obtained through the multiple linear regression equation, which can effectively characterize the damage degree of aluminum alloy materials.
针对铝合金材料,确定不同缺陷特征参数对铝合金疲劳寿命的影响,通过Pearson系数选取与疲劳损伤影响较大的特征参数进行研究,基于Abaqus分析软件建立了缺陷特征的疲劳损伤模型,在特征权重公式的基础上,获得了不同缺陷特征对疲劳性能影响的相对权重,最后达到量化不同缺陷特征对疲劳性能的影响。
汽车轻量化,缺陷特征,疲劳损伤,损伤模型
Hejun Wang
Shanghai University of Engineering Science, Shanghai
Received: Jan. 16th, 2023; accepted: Feb. 17th, 2023; published: Feb. 23rd, 2023
For aluminum alloy, the influence of different defect characteristic parameters on the fatigue life of aluminum alloy was determined. Pearson coefficient was selected and the characteristic parameters with great influence on fatigue damage were studied. The damage model of defect characteristics was established based on Abaqus analysis software. The relative weights of the effects of different defect characteristics on the fatigue properties were obtained. Finally, the mesoscopic damage factors were obtained through the multiple linear regression equation, which can effectively characterize the damage degree of aluminum alloy materials.
Keywords:Automotive Lightweight, Defect Characteristics, Fatigue Damage, Damage Model
Copyright © 2023 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
随着汽车轻量化研究的发展,铝合金材料广泛应用于汽车行业,且铝合金有着密度小、加工性能好同时强度高等优点,可以给汽车提供更轻的车身质量 [
针对缺陷特征对疲劳损伤和寿命的影响,很多学者都进行了广泛研究。缺陷产生一个应力强度因子,它取决于缺陷的大小,对于表面的缺陷,应力强度因子更大。Williams等 [
本次疲劳拉伸实验采用的是基于ASTM E8标准制备的6061铝合金试件。6061铝合金是一种以硅和镁为主要合金元素的铝合金与其他材料相比,具有易加工和焊接,易电镀,韧性较好,不易变形,容易抛光和表面处理,极强的耐腐蚀性的优点,是6系列铝合金中非常主要的一种合金产品。6061铝合金试件材料相比于其他系的铝合金,具有更好的加工性能以及抗腐蚀性等优良特点,常用于装饰、包装、建筑、运输、电子、航空、航天、兵器等各行各业。铝合金材料的主要成分如表1所示,铝合金材料的具体形状及尺寸如图1所示。
Al | Mg | Fe | Si | Cr | Cu | Zn | Mn | Ti | 其他 |
---|---|---|---|---|---|---|---|---|---|
95.8 | 0.8~1.2 | 0.7 | 0.6 | 0.3 | 0.2~0.4 | 0.25 | 0.15 | 0.15 | 0.15 |
表1. 6061-T6铝合金的成分
图1. 实验试件的尺寸图
通常情况下,铝合金材料的表面粗糙度越高的地方,很容易产生较大的应力集中,在受到较大的疲劳载荷时,会严重损害铝合金材料的疲劳性能。于是,为了避免铝合金材料表面粗糙度的不利影响而导致的实验误差,需要对试件进行抛光打磨后的工艺处理,然后采用粗糙度测量仪(图2)对铝合金试件的表面粗糙度测量,处理过后的试件表面粗糙度结果如表2所示。
图2. 表面粗糙度测量仪
试件组别 | 试件编号 | 粗糙度Ra/μm |
---|---|---|
第一组 | 1 | 0.395 |
2 | 0.421 | |
3 | 0.412 | |
第二组 | 4 | 0.390 |
5 | 0.402 | |
6 | 0.386 | |
第三组 | 7 | 0.395 |
8 | 0.423 | |
9 | 0.411 | |
第四组 | 10 | 0.421 |
11 | 0.396 | |
12 | 0.405 | |
第五组 | 13 | 0.429 |
14 | 0.412 | |
15 | 0.410 |
表2. 试件的表面粗糙度
基于相关的实验标准需要对试件进行静力学拉伸,得到了铝合金材料的基本力学参数,本次静力学拉伸实验所使用的设备是MTS Landmark (如图3所示),该设备是具备动拉伸和静拉伸多种应用功能,能够在恶劣(腐蚀、极寒、炎热)的环境下工作,应用场景主要是测试各种金属、非金属、复合材料及其结构件的动静态力学性能。可以进行腐蚀、高低温环境下的力学性能测试。
图3. MTS Landmark电液伺服动态测试系统
试件的拉伸采用恒定拉伸速度的方式进行,拉伸的速度设为0.2 mm/s。为避免实验的偶然性,每组数据取三个试件拉伸参数的平均值作为该试件最终的力学参数,如表3所示。并以此得出铝合金材料拉伸的应力–应变关系图,如图4所示。
为获取不同损伤阶段的试件,使用MTS Landmark电液伺服动态测试系统对铝合金材料试件进行阶段性的疲劳拉伸实验。基于已获得的力学参数对试件进行多次疲劳预实验,最终获得试件疲劳寿命约为10万次的加载载荷值为3600 N,因此,实验采用3600 N作为拉伸载荷对试件进行阶段性疲劳,实验相关的控制参数如表4所示。
拉伸强度Rm (MPa) | 屈服强度R0.2(MPa) | 弹性模量E (MPa) | 伸长率A (%) |
---|---|---|---|
339.8 | 326.6 | 58368 | 3.96 |
表3. 6061-T6铝合金的力学性能
图4. 6061-T6铝合金试件的应力–应变关系曲线
参数名称 | 温度 | 加载方式 | 应力比 | 频率 | 载荷 |
---|---|---|---|---|---|
参数 | 24℃ | 正弦波形 | 0.1 | 52 Hz | 3600 N |
表4. 疲劳实验的实验参数
因此,按照表4中的加载条件,将实验的疲劳阶段分为5组,分别对应疲劳循环次数的2万次、4万次、6万次、8万次和10万次,每个疲劳阶段3个试件,共15个试件,各试件如图5。
图5. 疲劳拉伸的5组试件
本文采用皮尔逊相关系数(Pearson)分析与疲劳损伤相关性较大的特征,皮尔逊相关系数的定义为:
ρ X , Y = cov ( X , Y ) σ x σ y = E [ ( X − E X ) ( Y − E Y ) ] σ x σ y = E ( X Y ) − E ( X ) E ( Y ) E ( X 2 ) − E 2 ( X ) E ( Y 2 ) − E 2 ( Y ) (1)
其中 cov 为协方差, σ 为标准差。皮尔逊相关系数有一重要特点,不受不同参数单位的影响,使计算的值分布在−1到1中间。根据疲劳拉伸实验获取得到的特征参数如表5所示,引入中间变量疲劳循环周次来表征材料的疲劳损伤程度,与缺陷特征(数目、位置、形状和孔隙率)进行相关性分析,将表5中缺陷特征参量代入到各疲劳阶段的缺陷特征损伤量代入到皮尔逊相关系数求解公式中,从而得到各个特征与疲劳损伤过程之间的相关系数,如图6所示。
循环次数(万次) | 2 | 4 | 6 | 8 | 10 | |
---|---|---|---|---|---|---|
孔洞数量/个 | 第一组 | 35,106 | 215,398 | 125,703 | 35,684 | 20,982 |
第二组 | 36,061 | 233,821 | 104,895 | 30,569 | 16,456 | |
第三组 | 34,854 | 190,636 | 76,695 | 24,534 | 18,672 | |
孔隙率/% | 第一组 | 0.035 | 0.166 | 0.268 | 0.338 | 0.343 |
第二组 | 0.076 | 0.141 | 0.219 | 0.288 | 0.376 | |
第三组 | 0.06 | 0.176 | 0.242 | 0.343 | 0.439 | |
最大孔洞体积/μm3 | 第一组 | 996 | 1965 | 2145 | 3360 | 4806 |
第二组 | 1117 | 2048 | 2238 | 3045 | 4463 | |
第三组 | 1064 | 1896 | 2405 | 3251 | 4907 | |
最大损伤表面积/μm2 | 第一组 | 0.0124 | 0.0231 | 0.0614 | 0.0895 | 0.1541 |
第二组 | 0.0097 | 0.0291 | 0.0541 | 0.0932 | 0.1932 | |
第三组 | 0.0113 | 0.0340 | 0.0481 | 0.0747 | 0.2056 |
表5. 各试件的部分缺陷特征和特征量
图6. 特征参数与循环次数之间的相关系数
从图6的相关系数表中可以看出,循环次数与缺陷数量特征的Pearson相关系数小于0.5,呈现弱相关性趋势,而缺陷的孔隙率,形状和位置呈现强相关性。因此,本文选取缺陷的孔隙率、位置和形状三个特征量三个特征量与疲劳损伤的相关系数值都在0.9以上,呈现强相关性。因此,本文将量化缺陷的孔隙率、形状和位置参数对铝合金材料疲劳损伤的影响。
随着加载循环周次的提高,内部缺陷的特征参量的数值也在不断变化。相关研究表明 [
假设材料的细观损伤程度为评价目标G,三个缺陷特征参数(孔隙率、形状和位置)为评价指标 { S 1 , S 2 , S 3 } 。设xg为评价目标G的取值变量,xi为指标Si的取值变量 ( i = 1 , 2 , 3 ) 。同时,细观损伤参量根据缺陷特征量的变化而变化。基于ABAQUS有限元分析软件建立疲劳损伤模型,使其中一个缺陷特征参量xi增加10%,细观损伤参量xg会对应增加ki%。则ki为缺陷特征指标Si对细观损伤评价目标G的弹性系数。计算方法如公式2所示:
ω i = k i ∑ n = 1 3 k n (2)
其中, ω i 是为该缺陷特征对疲劳损伤影响的权重。
根据上述权重公式,求出孔隙率、形状和位置分别对疲劳损伤的影响程度。
在疲劳拉伸过程中,试件不可避免地会出现大小和位置随机的缺陷。利用仿真软件模拟这些缺陷的变化是研究缺陷对试件疲劳损伤影响的有效方法。在研究中,利用ABAQUS仿真软件对试件的缺陷进行了仿真。试件模型参照CT扫描中试件中心面积,设定为长、宽、高分别为5 mm、2 mm、6 mm的矩形六面体,分别对应坐标轴的Y、X、Z轴,如图7所示。通过ABAQUS子例程导入缺陷的特征信息,生成所需的损伤模型进行仿真。
图7. 试件CT扫描区等效损伤模型
根据缺陷特征将仿真实验分为3组,每组通过改变特征参量进行5次仿真实验。在仿真中,采用位移控制方法对损伤模型进行有限元拉伸分析。固定模型的一端,在另一端施加应力P,如图8(a)所示。从最后的拉伸模拟结果,可得到最大应力点的 P max 值,如图8(b)所示。为了表征缺陷模型的疲劳性能,引入应力集中因子K作为中间变量,应力集中系数的计算公式如公式3所示:
K = P max P (3)
其中, P max 最大应力,P施加应力。
图8. (a) 张力模型,(b) 仿真结果
在仿真实验中,缺陷的分布如下图所示,在Z轴方向有三层缺陷,分别在Z轴的2.8 mm,3 mm,3.2 mm处,每层均有160个缺陷,在XY平面上的排布方式如图9所示。
图9. 模型中缺陷在xy平面上的分布
根据权重计算方法,缺陷的特征量每增加10%。为了避免缺陷间距的相互影响,当缺陷的孔隙率增大时,在保持原有缺陷位置不变的情况下,按比例增加缺陷。同时保证了相同的参数设置和相同的网格划分方法。最后得到缺陷孔隙度特征、形状(纵横比)、位置变化的5个最大应力集中系数,如图10所示。
如图10(a)所示随着铝合金试件加载循环次数的增加,缺陷的总体积随着加载次数的增加而增加。由于缺陷的总体积(孔隙率)增大,材料应力集中系数急剧增大,大幅度降低材料的疲劳损伤;如图10(b)所示,缺陷在加载循环的初期,大部分缺陷的初始形状接近为圆形,缺陷的纵横比为1,但是由于加载循环持续的不断增加,导致缺陷在不断被拉长,拉长的过程中与其他缺陷连接一起,同时缺陷在拉长过程中尖端产生恶劣的应力集中,减少铝合金材料的疲劳性能;SERRANO [
图10. 不同缺陷特征对应的应力集中系数的变化,(a) 孔隙率,(b) 形状,(c) 位置
由上述的仿真分析输出了不同缺陷参数的应力集中系数,采用最小二乘法将曲线拟合成一条直线,拟合后的曲线如图10所示。三个曲线的线性函数如公式4所示。
{ y 1 = 39900 x 1 + 0.3080 y 2 = 0.6793 x 2 + 1.5998 y 3 = − 1.649 x 3 + 3.0851 (4)
其中,x1、x2、x3分别是缺陷的孔隙率、形状、位置的特征参量,同理,y1、y2、y3是三个特征参量对应的应力集中系数。通过公式4可得到孔隙度、形状和位置缺陷特征对细观损伤的弹性系数k1、k2和k3分别为8.38、2.98和5.98。基于弹性系数可以得到孔隙度、形状和位置对疲劳损伤的权重分别为0.483、0.172和0.345。
本文通过Pearson系数选取相关性最大的缺陷特征参数,得到缺陷的孔隙率、位置以及形状和疲劳性能的相关性较高。基于改变缺陷的特征参量建立仿真损伤模型,利用应力集中系数表征损伤的程度,最后达到量化不同缺陷特征对疲劳性能的影响。
王何俊. 缺陷特征对铝合金疲劳性能的研究Study on Fatigue Performance of Aluminum Alloy with Defect Characteristics[J]. 运筹与模糊学, 2023, 13(01): 361-370. https://doi.org/10.12677/ORF.2023.131038
https://doi.org/10.1016/j.scriptamat.2016.05.002
https://doi.org/10.1038/s41598-017-06504-5
https://doi.org/10.1016/j.prostr.2017.11.063
https://doi.org/10.1016/j.ijfatigue.2020.105584
https://doi.org/10.1038/srep45239