表面粗糙度对表面应力集中系数和疲劳寿命影响分析文档格式.docx
- 文档编号:5116948
- 上传时间:2023-05-04
- 格式:DOCX
- 页数:10
- 大小:44.14KB
表面粗糙度对表面应力集中系数和疲劳寿命影响分析文档格式.docx
《表面粗糙度对表面应力集中系数和疲劳寿命影响分析文档格式.docx》由会员分享,可在线阅读,更多相关《表面粗糙度对表面应力集中系数和疲劳寿命影响分析文档格式.docx(10页珍藏版)》请在冰点文库上搜索。
本文对不同表面粗糙度下的模型进行唯象统计分析,旨在建立微缺口中心间距、粗糙度Rz与应力集中系数Kt之间的定量关系,进而分析一定存活率下,缺口疲劳寿命与表面粗糙度的关系。
2、计算模型
图1模型尺寸及受载示意图
Fig.1Dimensionandloadofmodels
分析计算用商用有限元分析软件ABAQUS来完成。
表面形貌被简化为半椭圆形微缺口。
如下图1所示:
在边长L为1mm平板的上边缘中间处有一组长短半径分别为a和b的半椭圆缺口;
短半径b即为缺口深度,长半径a变化范围为10至50?
m,b/a在0.2~1范围内变化。
多缺口中心间距d(下文简称间距)取1~5倍缺口宽度2a。
结构两边受均布拉伸载荷P,大小为100MPa。
材料为LY12CZ铝合金,其弹性模量E=68GPa,泊松比μ=0.3。
对模型进行有限元网格划分时,先分区,再用不同划分方法和单元大小来刻画不同部分,以实现计算速度和质量的统一。
划分时远离缺口部分选择二次减缩积分单元CPS8R,缺口附近的部分采用尺寸为1微米的二次完全积分单元CPS8来进行计算。
在缺口附近部分采用自由划分法,其他位置采用结构划分法;
同时选择四边形单元[10]。
a为20?
m,b为10?
m单微缺口的模型网格划分最后结果如图2所示,缺口部分细节如图3所示。
图2单微缺口网格划分图示图3单微缺口的局部细化网格
Fig.2FiniteelementmodelFig.3Localmeshofmicro-notch
3、应力集中系数受影响分析
3.1半径对单微缺口应力集中系数的影响
建立单微缺口长短半径变化的多个模型,得到缺口处应力集中系数。
如图4长短半径为20?
m,10?
m的缺口部分Mises应力分布。
在此用Mises应力作为局部最大应力即可较为准确的求出缺口处的应力集中系数值[8]。
图4单微缺口Mises应力云图分布(局部)
Fig.4Misesstresscounterofmicro-notch
将图4所得的?
max?
?
201.1MPa和载荷P=100MP(与净截面名义应力误差足够小,可以不考虑)代入应力集中系数计算式中即得到此时Kt=2.011[11]。
该结果与式
(1)中取λ=1,n=2,?
=a2/b(缺口底部顶点曲率半径),Rz=b的单缺口情况得到Kt=2相比,可看到两者相差不大,说明模型单元网格划分是合适的。
同理通过长半径a取值从10?
m到50?
m,Rz=b=10?
m,且b/a逐渐减小时的应力集中系数计算得到表1所示结果;
其中在b/a=0.5时增加计算b=5?
m和20?
m两组数据,确保以b=10?
m,a可变来实现b/a变化的正确性。
表1不同缺口长短半径比时的Kt
Tab.1Ktwithdifferentratiobetweendepthandbreadth
b/a
1
2/3
0.5
2/5
1/3
2/7
1/4
2/9
1/5
a/?
m
10
15
20
40
25
30
35
45
50
Kt
3.04
2.36
2.03
2.01
1.81
1.68
1.58
1.51
1.45
1.41
图5Kt与b/a关系曲线图
Fig.5b/avsstressconcentrationfactor
从图5中可以看到缺口短长半径比与Kt的关系用最小二乘法拟合得相关系数达0.998以上的线性关系式
(4)
将式(4)与式
(1)进行对比,式(4)的结果可为
(1)式中n加修正项所得。
在式
(1)的基础上,建立修正公式:
(5)
将b/a=0.5时b=5?
m代入式(5)中,得到结果与表1中有限元解比较,发现两组误差均不超过1%。
此外将b=2?
m,a=10?
m和b=50?
m,a=50?
m的两种极限情况代入(5)式得Kt=1.4088和3.044,与采用相同单元尺寸的有限元算法得到的1.395和3.147误差也仅为1%和3%。
这些说明式(5)准确,但该公式只适用于一定尺寸下单微缺口的情况,对多微缺口的状态需要进一步修正。
3.2多微缺口存在时的Kt及其受表面粗糙度影响
根据单变量分析法,只改变缺口的中心间距d,研究反映表面粗糙度的参数即单微缺口应力集中系数(Kt单)和缺口间距与缺口宽度比值(d/2a)对多微缺口状况Kt影响。
首先分析缺口个数n对Kt的影响。
取缺口中心间距d为40至200微米范围内某值不变,在实体上边缘中心附近有a=20?
m,b=10?
m的2至4个连续微缺口,缺口分布为靠中心对称分布,见图1右所示。
建立模型计算得到多组Kt,数据见表2所示,并绘制相应曲线图6。
表2不同数目缺口时的Kt
Tab.2Ktwithdifferentnumberofnotches
缺口数目
Numbern
应力集中系数StressConcentrationFactorKt
d=40
d=80
d=120
d=160
d=200
2.011
2
1.845
1.920
1.948
1.988
1.998
3
1.803
1.902
1.934
1.985
4
1.789
1.895
1.932
1.986
2.009
图6缺口数目n和Kt关系图示
Fig.6Numberofnotchesvsstressconcentrationfactor
从表2中和得到的图6结果可以看到多个缺口的存在可以降低单缺口存在时的应力集中系数;
一定范围内(从表中看是d取200?
m即5倍以内时),缺口中心距增大,Kt逐渐变大,超出该范围Kt受间距影响就很小(8倍缺口宽度后几乎不受影响);
缺口数目增加时,Kt减小;
且缺口中心间距足够大(比如图中d取120?
m为3倍缺口宽度时),缺口的数目的增加对Kt的影响很小,这些正是后面选取合适间距做计算的依据。
此外,从曲线趋势上可以看到:
缺口数目增大到一定程度后,恒定间距多微缺口的Kt值趋于稳定,为此在定量分析时就应该保证足够多的微缺口数以较好模拟表面的粗糙特征,以得到准确数据。
故选取10个相同的等间距微缺口,分布如图1右所示;
采用对称分析左端轴线约束,右侧五个缺口均匀分布的模型,其他条件与前面相同。
缺口宽度为2a。
a=10?
m,b=5?
m不变。
计算得到下面不同间距时的Kt,见表3所示。
表3不同缺口间距下的Kt
Tab.3Stressconcentrationfactorwithdifferentspacebetweennotches
d/?
60
70
80
90
100
1.5
2.5
3.5
4.5
5
1.823
1.863
1.914
1.931
1.950
1.965
1.977
2.017
考虑多微缺口时应力集中系数Kt与同条件单微缺口的Kt单有下面关系
(6)
式中Kt为多微缺口时的应力集中系数;
为比例系数;
Kt单为单微缺口时应力集中系数。
因此将所求得Kt与式(5)中a=10?
m的单微缺口Kt单=2.022进行合适拟合,得到如图7所示的线性关系,即得到与d/2a有关的系数参量?
,进而最终得到多微缺口的应力集中系数经验公式。
图7Kt与缺口间距关系拟合曲线
Fig.7Spacebetweennotchesvsstressconcentrationfactor
结合式(6)得到相关系数达0.987的连续相邻多微缺口与单微缺口应力集中系数关系
(7)将公式(5)代入就得到表面粗糙度对应力集中系数影响的最终修正公式
(8)
选取间距d为90?
m,长短半径a,b分别为15?
m,5?
m半椭圆缺口(d/2a=3)来检验。
用有限元方法得到Kt0=1.671,与公式(8)得到的Kt=1.587的误差不超过5.0%,说明该经验公式是适用的。
4、疲劳寿命预测
在分析多微缺口疲劳寿命之前应该首先确定危险缺口。
对a=10?
m的左右五个缺口均匀分布的对称模型进行分析,计算得到三组不同间距下由外向内五个缺口处的应力集中系数Kt,按图8中编号给出顺序缺口的Kt变化数值如表4所示,变化曲线如图9所示。
图8缺口编号及Mises应力云图(d=50?
m)
Fig.8NumberandMisesstresscounterofmicro-notches
表4多缺口不同间距下Kt
Tab.4Ktfornotcheswithdifferentspace
缺口编号Number
d=30?
d=50?
d=80?
1.832
1.673
1.815
1.946
1.638
1.834
1.928
1.613
1.762
1.933
1.623
1.812
1.925
图9不同缺口的Kt变化
Fig.9Thechangesofdifferentnotches
从表4和图9可以看到对于表面多微缺口状况,最外边1号缺口处的应力集中程度最高,显然为危险缺口,分析缺口疲劳寿命应该重点分析这个最可能破坏处。
根据上面的分析,将加工后LY12CZ铝合金的表面纹理组织简化为宽度等于50?
m,中心间距与宽度的比d/2a为3的连续多微缺口模型。
研究不同粗糙度RZ时缺口的疲劳寿命。
在式
(2)和式(3)基础上,文献[12]根据实验给出存活率(
)分别为95%、99.9%的LY12CZ铝合金切口件疲劳起始寿命表达式
95%
(9)
99.9%
(10)
由公式(3)所给的当量应力幅
可得在疲劳载荷为R=-1,Smax=100MPa时有式
(11)
再根据式(8)及模型在不同粗糙度RZ下的应力场,结合式(9)、式(10)和式(11)就可以得到微缺口在存活率为95%、99.9%时的疲劳寿命如表5所示。
处理后得到粗糙度RZ与缺口疲劳寿命对数关系曲线如图10所示。
表5各个模型的疲劳寿命
Tab.5Fatiguelifeofmodels
RZ/?
/MPa
N(95%)/cycle
N(99.9%)/cycle
1.33
133
6211419
1.52
152
2065861
1.72
172
3274377
927460
1.91
191
1254313
514568
2.10
210
629256
316164
图10粗糙度RZ与缺口的疲劳寿命关系
Fig.10Roughnessvsfatiguelifeofnotch
从图10可以看到,粗糙度RZ增大,缺口的疲劳寿命减少。
存活率一定时,表面粗糙Rz与缺口的对数疲劳寿命之间有二次曲线拟合关系式
(12)
(13)
式(12)和(13)的相关系数均在0.98以上。
对连续椭圆多微缺口化的表面纹理,将Rz=23um代入式(12)和式(13),可得到缺口在不同存活率下的疲劳寿命分别为1528797(95%Sv),632280(99.9%Sv)。
与Rz=23?
m时采用上面相同算法所得到的不同存活率下缺口的疲劳寿命1643383(95%Sv)、649575(99.9%Sv)进行对比,其相对误差分别为7%、3%。
相对误差的平均值为5%。
可以看到,用二次曲线来拟合缺口的表面粗糙度与缺口的对数疲劳寿命关系,其相对误差很小,通过验证平均值在5%左右。
因此,用二次曲线来拟合表面粗糙度与缺口的对数疲劳寿命关系比较合适。
即可以认为,表面粗糙度与缺口的对数疲劳寿命呈二次曲线关系。
5、结论
通过把平板表面纹理形貌特点连续相邻多微缺口化,采用有限单元方法计算得到不同表面粗糙度下的表面应力集中系数Kt,进而对缺口疲劳寿命进行了具有一定存活率的估算。
最后模拟分析有以下结果:
1).多微缺口可以降低单缺口存在时的应力集中系数,且外沿缺口应力集中程度高。
一定范围内,缺口半径保持不变时,缺口中心距的增大,Kt逐渐变大,超过该范围即有Kt达到单缺口水平而不再受影响。
2).缺口数目增加时,Kt减小,当数目增大到一定程度后Kt停止变化;
缺口中心间距越大,缺口的数目增加对Kt的影响将减小。
3).建立了表面粗糙度与表面应力集中系数之间的经验公式(8)。
4).疲劳寿命随表面粗糙度增大而降低。
在相同存活率下,表面粗糙度与表面微缺口的对数疲劳寿命呈二次曲线关系。
6、小结与展望
通过本次制作的论文,我深刻的发现自己的能力水平相对比较低,还有很多不足之处,在以后的学习中,我应该更加的努力学习,打好基础,做好自己应该做的。
参考文献
1.吴松青.表面粗糙度应用指南.北京:
机械工业出版社,1990
2.张东初,裴旭明.加工工艺对表面粗糙度及疲劳寿命的影响.中国机械工程,2003
3.郁大照,陈跃良,多缺口应力集中系数有限元研究.强度与环境,2002,
4.庄茁,张帆,岑松,等.ABAQUS非线性有限元分析与实例.北京:
科学出版社.2003
5.魏建锋,郑修麟,吕宝桐.LY12CZ铝合金切口件的P—S—N曲线表达式及寿命估算.机械工程材料,1997
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 表面 粗糙 应力 集中 系数 疲劳 寿命 影响 分析