2A12铝合金微动疲劳全寿命预测方法研究什么的问候
陈跃良;徐丽;张勇;刘旭
【摘 要】目的：对于2A12铝合金，提出基于成核寿命和扩展寿命的微动疲劳全寿命预测方法。方法基于损伤力学法计算裂纹成核寿命，利用扩展有限元计算裂纹尖端应力强度因子，应用断裂力学计算裂纹扩展寿命，并对预测者和试验值进行比较。结果损伤力学法能考虑接触面应力三维度的作用来反映多轴状态作用，能有效模拟微动疲劳多轴行为。基于损伤力学法的微动疲劳全寿命预测模型能有效预测微动疲劳全寿命。由于微动作用，裂纹成核非常早，扩展寿命从试件的近表面开始，占全寿命的主要部分。结论考虑成核寿命和扩展寿命的微动疲劳全寿命分析是完善的，预测值与试验值比较吻合。%Objective The holistic life prediction approach for fretting fatigue was propod for 2A12 aluminum alloy bad on initiation life and propagation life. Methods The initiation life was calculated using continuum damage mechanics approach, and the propagation life was calculated using fracture mechanics approach with SIF obtained from XFEM. The prediction life was compared with the test life. Results The results indicated that the damage mechanics appro
ach could consider the effect of the stress triaxiality which reflected the multiaxial stress state property. The model could effectively predict the fretting fatigue holistic life. The crack initiated quickly becau of fretting, and the crack propagation began from the near surface, and the propagation life accounted for the main part of holistic life. Conclusion The holistic life prediction model considering initiation life and propagation life was perfect, and the prediction life coincided well with the test life.
【期刊名称】人脉《装备环境工程》
【年(卷),期】2014(000)004
【总页数】6页(P1-6)
【关键词】2A12 铝合金;损伤力学法;微动疲劳;全寿命
【作 者】陈跃良;徐丽;张勇;刘旭
【作者单位】海军航空工程学院 青岛校区，山东 青岛 266041;海军航空工程学院 青岛校
雨韵区，山东 青岛 266041; 海军航空兵学院，辽宁 葫芦岛 125001;海军航空工程学院 青岛校区，山东 青岛 266041;海军航空工程学院 青岛校区，山东 青岛 266041
【正文语种】中 文
【中图分类】TG146.2+1
关于微动疲劳寿命的研究,目前应用较多的方法主要有:局部应力-应变法、名义应力法、断裂力学法、基于临界面的疲劳参数法以及基于损伤力学法的微动疲劳寿命预测方法[1—4]。
qq图片头像
微动疲劳寿命是裂纹成核寿命和扩展寿命的总和,但在寿命预测中往往只考虑成核寿命或只考虑扩展寿命。文献[5—7]基于断裂力学进行微动疲劳扩展寿命预测,研究发现裂纹扩展寿命占全寿命的主要部分,然而,Hill[8]和Szolwinski[9]等其他研究者支持微动疲劳是一个成核控制的过程这一论点,认为微动疲劳成核寿命占全寿命的主导部分。产生这一矛盾的主要原因与材料本身的特性和初始裂纹的选取有关。近年来,Houghton和Wavish[10]又提出了微动疲劳全寿命预测方法,其基本原理是用多轴疲劳参数计算裂纹成核寿命,用线弹性断裂力学计算裂纹扩展寿命,并考虑短裂纹作用进行修正。
对于2A12铝合金,文中提出了基于损伤力学和断裂力学的微动疲劳全寿命预测方法,利用扩展有限元计算应力强度因子,结合全寿命预测模型计算微动疲劳寿命,并分析了裂纹成核寿命、扩展寿命比例。
材料损伤主要是微裂纹和夹杂成核及扩展累积的过程。Kachanov[11]首先引入损伤变量,考虑应力应变建立本构方程来预测裂纹成核。
损耗势能函数φ为损伤弹性能释放率Y、塑性应变率˙P、累积微塑性应变˙π等3个主要变量的函数,根据热力学可表示为[12]:
式中:C和β为材料常数;Y为与损伤参数D有关的热力学变量,根据热力学势函数(ψ)可表示为:
式中:ρ为密度;E为弹性模量;σ*为等效多轴损伤应力。用Von Mis等效应力σeq和三维函数Rv表示为:
式中Rv可表示为:
式中:σH为静水应力。徐州加勒比>咸鸭蛋煮多久
在微动疲劳情况下,假定微观塑性为线性累积,则损伤演化方程为:
式中:A=1/(KmCβ(2E)β)。
损伤变量为载荷循环数的函数,当N=Ni→D= 1时,Ni为微动疲劳宏观裂纹成核寿命,表示为:
式中:σeq-max和σeq-min为一个循环内的最大Von Mis等效应力和最小Von Mis等效应力,Rv为三维函数平均值,σeq-max和σeq-min以及Rv的值都可以根据有限元模型计算得到;m为Ramberg-Osgood方程中的幂常数,通常m=2;A和β为损伤参数,可通过用方程(6)对根据试验计算得到的成核寿命的拟合得到。
用断裂力学方法进行微动疲劳扩展寿命预测是应用较多的一种方法。首先确定裂纹萌生的位置,然后计算裂纹尖端应力强度因子的范围ΔK。根据断裂力学方法,裂纹从假定初始裂纹长度ai扩展到临界裂纹长度af的寿命为:
式中:f(ΔK)为给定的裂纹扩展准则;C,m是材料常数。文中选用Paris裂纹扩展准则来计算扩展寿命,假定裂纹以模型Ⅰ模式扩展。在计算扩展寿命时,必须计算得到相应的应力强度因子(SIF)值。由于循环载荷裂纹的闭合作用(Kmin=0),则相应的SIF范围等于Kmax。
扩展有限元可以通过数值计算来体现不连续状态,而不需要重新划分网格。这个方法由 Moës, Belystschko等[13]提出,是基于单位分解法,该方法能在标准有限元结构中导入局部节点集。研究表明,该方法能准确地计算应力强度因子[14]。
利用扩展有限元法,应力强度因子SIF计算如下:
对于平面应力:E＇=E,对于平面应变:E＇=E/(1-υ2)。
4.1 计算方法
独行踽踽
基于损伤力学和断裂力学的全寿命预测为:用损伤模型来计算损伤演化至宏观裂纹成核,根据从有限元模型中计算的应力应变,基于断裂力学估算微动疲劳裂纹扩展寿命;然后根据试验结果反推裂纹成核寿命,结合成核寿命方程,根据得到的成核寿命拟合得到损伤参数A和β;再利用成核寿命方程进行成核寿命计算,成核寿命和扩展寿命相加即为预测的全寿命。在计算裂纹扩展寿命时,假定初始裂纹长度ai=100 μm,根据最大等效塑性应变位置确定裂纹成核位置,裂纹方向与接触面垂直。2A12铝合金的弹性模量E=69.6 GPa,泊松比为υ= 0.33,极限强度σu=447.86 MPa,屈服强度σy= 342.02 MPa[15],扩展直到试件断裂af=3.5 mm。文献[16]研
究表明,初始裂纹的选取对裂纹演化预测有一定的影响,但对最终寿命没有影响。国外通常用涡电流无损检测来检测搭接件间的微裂纹,最小可检尺寸为760 μm×380 μm的半椭圆裂纹[17],而且没有特定的准则来定义微动疲劳初始裂纹长度,介于上述原因,将初始裂纹长度选为100 μm,裂纹扩展增量Δa为10 μm。
选用2A12铝合金材料进行微动疲劳试验,微动垫半径r分别为115,180,225 mm。具体实验设备和实验方法见文献[18],试验结果和根据断裂力学计算出的扩展寿命见表1。
估算微动疲劳全寿命的具体步骤为:
1)用ABAQUS有限元软件模拟微动疲劳接触;
2)用El Haddad参数确定过程区尺寸lc;
3)根据有限元模型从过程区中得到σeq-max,σeq-min,σ*和Rv的平均值;
4)用 XFEM模拟裂纹扩展得到扩展寿命(Np-calculated);
5) 反推裂纹成核寿命Ni-experimental=Nf-experimental-Np-calculated;
6)选取其中4个试件,根据计算的裂纹成核寿命用方程(6)来拟合得到损伤参数A和β;
7)根据拟合得到的损伤参数值和从(3)中得到的各参数,利用方程(6)可求出其余8个试件的成核寿命(Ni-estimated);
8)相加计算成核寿命和扩展寿命,得到预测的微动疲劳全寿命为Nf-predicted=Ni-estimated+Np-calculated。
4.2 参数确定
对于圆柱/平面接触微动疲劳,在接触载荷和轴向载荷作用下,接触边缘产生微观尺度的局部塑性变形,从而产生较高的应力梯度,则应力分量σeq-max,σeq-min和Rv受网格划分尺度的影响很大。为了减少网格划分尺度的影响,应用基于过程区(面)法(如图1所示),以初始裂纹成核点为过程区中心,以临界距离lc为过程区半径r,计算半圆过程区内应力分量σeq-max,σeq-min和Rv的平均值。 临界距离lc根据El Haddad参数[19]来确定,其定义为:心得报告
式中:ΔKth为长裂纹门槛应力强度因子;Δσe为光滑试件疲劳极限。对于2A12铝合金,ΔKth=2.47 MPa(m)1/2,Δσe=137 MPa,则lc=100 μm。
通过对试件有限元模型进行计算,从有限元中得到的σeq-max,σeq-min和Rv的平均值见表2。
因为损伤参数A和β为未知材料常数,从12个试验件中选取4个试验件来确定材料的损伤参数。试件包含4种不同的轴向应力状态,法向接触载荷为450 N,微动垫半径为180 mm。根据表3所选试验件过程区内σeq-max,σeq-min和Rv的平均值,将这些值代入到寿命计算方程(6)中,即可用回归分析计算试验成核寿命 Ni-experimental与估算成核寿命Ni-estimated间的最小二乘拟合,从而确定损伤参数的值。试验值与估算值的误差为回归差(误差=Ni-experimental-Ni-estimated),误差总量用均方根差来估计,其公式为:
式中:n为分散数据量。非线性回归使RMSE取最小值,得到的损伤参数值见表3。损伤参数A和β与等效多轴损伤应力σ*之间的关系曲线分别如图2和图3所示,由此计算出剩余试件的损伤参数值。
