表面粗糙度对表面应力集中系数和疲劳寿命影响分析.docx

表面粗糙度对表面应力集中系数和疲劳寿命影响分析.docx

《表面粗糙度对表面应力集中系数和疲劳寿命影响分析.docx》由会员分享，可在线阅读，更多相关《表面粗糙度对表面应力集中系数和疲劳寿命影响分析.docx（8页珍藏版）》请在冰豆网上搜索。

表面粗糙度对表面应力集中系数和疲劳寿命影响分析

黄河交通学院

学课论文（论文）

题目：

表面粗糙度对表面应力集中系数和疲劳寿命影响分析

完成人：

谷玉乐

班级：

14级车辆工程一班

学制：

4年

专业：

汽车工程学院车辆工程

指导教师：

XXX

完成日期：

2016年5月29日

表面粗糙度对表面应力集中系数和疲劳寿命影响分析

摘要把平板的表面形貌简化为半椭圆形微缺口，采用有限元法对不同表面粗

糙度下应力场进行分析，以得到不同表面粗糙度Rz时平板的表面应力集中系数，并预测出不同表面粗糙度的平板在不同存活率下的疲劳寿命，讨论表面粗糙度对

疲劳寿命的影响。

采用回归方法建立了表面应力集中系数与微缺口中心间距、表

面粗糙度之间的经验公式。

结果表明，在相同存活率下，表面粗糙度与平板的对数疲劳寿命呈二次曲线关系。

关键词表面粗糙度微缺口应力集中系数疲劳寿命

1、引言

表面粗糙度是反映零件表面微观几何形状误差的一个重要指标⑴7-8。

由断裂力学可知，表面粗糙度值愈大，表面的沟痕愈深，纹底半径愈小，应力集中越严重，抗疲劳破坏的能力就愈差。

因此表面粗糙度增大，会降低零件的疲劳强度［2］。

这一影响引起了广泛的关注。

在实验方面，Noll比较了低碳合金钢的表面分别在研磨、机加、热滚压和锻造这四种表面加工条件下，其应力-寿命关系曲线。

发现在高应力水平下，所有试件的寿命相近，在低应力区，寿命差异十分显着。

在应力幅值接近屈服强度的80%寸，研磨试件的寿命是锻造试件寿命的20倍［3］。

在数值计算方面，Andrews将表面粗糙度作为微观缺口，研究缺口对疲劳寿命的影响⑷。

文献［5］采用试验测试和有限元的方法，建立实际表面形貌的有限元网格，将表面犁沟视为微观缺口，研究裂纹萌生部位与缺口应力集中系和应力场的关系。

由此可见，

要确定表面粗糙度对疲劳寿命的影响，首先要建立表面粗糙度与表面应力集中系数的定量关系。

本文对不同表面粗糙度下的模型进行唯象统计分析，旨在建立微

缺口中心间距、粗糙度RZ与应力集中系数Kt之间的定量关系，进而

分析

存活

图1模型尺寸及受载示意图

Fig.1Dimensionandloadofmodels

疲劳

与表糙度系。

缺口

寿命

面粗的关

2、计算模型

分析计算用商用有限元分析软件ABAQU来完成。

表面形貌被简化为半椭圆

形微缺口。

如下图1所示：

在边长L为1mm平板的上边缘中间处有一组长短半径分别为a和b的半椭圆缺口；短半径b即为缺口深度，长半径a变化范围为10至50?

mb/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?

m10?

m的缺口部分Mises应力分布。

在此用Mises应力作为局部最大应力即可较为准确的求出缺口处的应力集中系数值⑹。

图4单微缺口Mises应力云图分布（局部）

Fig.4Misesstresscounterofmicro-notch

将图4所得的？

max?

?

201.1MPa和载荷P=100MP与净截面名义应力误差足够小，

可以不考虑）代入应力集中系数计算式中即得到此时K=2.011[11]。

该结果与式

（1）

中取入=1,n=2,?

=a2/b（缺口底部顶点曲率半径），Rz=b的单缺口情况得到K=2相比，可看到两者相差不大，说明模型单元网格划分是合适的。

同理通过长半径

a取值从10?

m到50?

m,R=b=10?

m且b/a逐渐减小时的应力集中系数计算得到表1所示结果；其中在b/a=0.5时增加计算b=5?

m和20?

m两组数据，确保以b=10?

ma可变来实现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

10

20

40

25

30

35

40

45

50

Kt

3.04

2.36

2.03

2.01

2.03

1.81

1.68

1.58

1.51

1.45

1.41

图5Kt与b/a关系曲线图

Fig.5b/avsstressconcentrationfactor

从图5中可以看到缺口短长半径比与K的关系用最小二乘法拟合得相关系数达0.998以上的线性关系式

Kt=12.044b/a（4）

将式（4）与式

（1）进行对比，式（4）的结果可为

（1）式中n加修正项所得。

在式

（1）的基础上，建立修正公式：

Kt=12.044,:

（5）

将b/a=0.5时b=5?

m和20?

m代入式（5）中，得到结果与表1中有限元解比较，发现两组误差均不超过1%此外将b=2?

m,a=10?

mffib=50?

m,a=50?

mB勺两种极限情况代入（5）式得K=1.4088和3.044,与采用相同单元尺寸的有限元算法得到的1.395和3.147误差也仅为1箱口3%这些说明式（5）准确，但该公式只适用于一定尺寸下单微缺口的情况，对多微缺口的状态需要进一步修正。

3.2多微缺口存在时的Kt及其受表面粗糙度影响

根据单变量分析法，只改变缺口的中心间距d，研究反映表面粗糙度的参数即单微缺口应力集中系数（K单）和缺口间距与缺口宽度比值（d/2a）对多微缺口状况Kt影响。

首先分析缺口个数n对Kt的影响。

取缺口中心间距d为40至200微米范围内某值不变，在实体上边缘中心附近有a=20?

m,b=10?

m勺2至4个连续微缺口，缺口分布为靠中心对称分布，见图1右所示。

建立模型计算得到多组K，数据见表2所示，并绘制相应曲线图6。

表2不同数目缺口时的K

Tab.2Kwithdifferentnumberofnotches

缺口数目

Numbern

应力集中系数StressConcentrationFactorKt

d=40

d=80

d=120

d=160

d=200

1

2.011

2.011

2.011

2.011

2.011

2

1.845

1.920

1.948

1.988

1.998

3

1.803

1.902

1.934

1.985

1.998

4

1.789

1.895

1.932

1.986

2.009

图6缺口数目n和Kt关系图示

Fig.6Numberofnotchesvsstressconcentrationfactor

从表2中和得到的图6结果可以看到多个缺口的存在可以降低单缺口存在时

的应力集中系数；一定范围内（从表中看是d取200?

m即5倍以内时），缺口中

心距增大，K逐渐变大，超出该范围Kt受间距影响就很小（8倍缺口宽度后几乎不受影响）；缺口数目增加时,Kt减小；且缺口中心间距足够大（比如图中d取120?

m为3倍缺口宽度时），缺口的数目的增加对Kt的影响很小，这些正是后面

选取合适间距做计算的依据。

此外，从曲线趋势上可以看到：

缺口数目增大到一定程度后，恒定间距多微缺口的K值趋于稳定，为此在定量分析时就应该保证

足够多的微缺口数以较好模拟表面的粗糙特征，以得到准确数据。

故选取10个相同的等间距微缺口，分布如图1右所示；采用对称分析左端轴线约束，右侧五个缺口均匀分布的模型，其他条件与前面相同。

缺口宽度为2a。

a=10?

mb=5?

m不变。

计算得到下面不同间距时的K，见表3所示。

表3不同缺口间距下的Kt

Tab.3Stressconcentrationfactorwithdifferentspacebetweennotches

d/?

m

30

40

50

60

70

80

90

100

1.5

2

2.5

3

3.5

4

4.5

5

1.823

1.863

1.914

1.931

1.950

1.965

1.977

2.017

考虑多微缺口时应力集中系数Kt与同条件单微缺口的Kt单有下面关系

Kt=：

Kt单（6）

式中K为多微缺口时的应力集中系数；？

为比例系数；Kt单为单微缺口时应力集中系数。

因此将所求得Kt与式（5）中a=10?

mb=5?

m的单微缺口K单=2.022进行合适拟合,得到如图7所示的线性关系，即得到与d/2a有关的系数参量?

，进而最终得到多微缺口的应力集中系数经验公式。

图7K与缺口间距关系拟合曲线

Fig.7Spacebetweennotchesvsstressconcentrationfactor

结合式（6）得到相关系数达0.987的连续相邻多微缺口与单微缺口应力集中系数关系

Kt=[0.8690.025（3）]心单（7）将公式（5）代入就2a

得到表面粗糙度对应力集中系数影响的最终修正公式

Kt=[0.8690.025（#）]（12.044Rz）（8）

选取间距d为90?

m长短半径a,b分别为15?

m,5?

m半椭圆缺口（d/2a=3）

来检验。

用有限元方法得到Kt0=1.671，与公式（8）得到的K=1.587的误差不超

过5.0%，说明该经验公式是适用的。

4、疲劳寿命预测

在分析多微缺口疲劳寿命之前应该首先确定危险缺口。

对a=10?

mb=5?

m的

左右五个缺口均匀分布的对称模型进行分析，计算得到三组不同间距下由外向内

五个缺口处的应力集中系数Kt，按图8中编号给出顺序缺口的K变化数值如表4

所示，变化曲线如图9所示。

图8缺口编号及Mises应力云图（d=50?

m）

Fig.8NumberandMisesstresscounterofmicro-notches

缺口编号

Number

d=30?

m

d=50?

m

d=80?

m

1

1.832

1.914

1.965

2

1.673

1.815

1.946

3

1.638

1.834

1.928

4

1.613

1.762

1.933

5

1.623

1.812

1.925

表4多缺口不同间距下Kt

Tab.4Kfornotcheswithdifferentspace

图9不同缺口的Kt变化

Fig.9Thechangesofdifferentnotches

从表4和图9可以看到对于表面多微缺口状况，最外边1号缺口处的应力集中程度最高，显然为危险缺口，分析缺口疲劳寿命应该重点分析这个最可能破坏处。

根据上面的分析，将加工后

LY12CZ铝合金的表面纹理组织简化

为宽度等于50?

m,中心间距与宽度的比d/2a为3的连续多微缺口模型。

研究不同粗糙度艮时缺口的疲劳寿命。

在式

（2）和式（3）基础上，文献［12］根据实验给出存活率（Sv）分别为

95%99.9%的LY12CZ铝合金切口件疲劳起始寿命表达式

N=3.301013Iveqv178-1371.78严95%S/（9）

N=3.0210131：

匚「8-1031.78严99.9%Sv（10）

由公式（3）所给的当量应力幅厶二eqv可得在疲劳载荷为R=-1，100MPa

时有式

*eqv/OOKt（11）

再根据式（8）及模型在不同粗糙度RZ下的应力场，结合式（9）、式（10）和式（11）就可以得到微缺口在存活率为95%99.9%寸的疲劳寿命如表5所示。

处理后得到粗糙度RZ与缺口疲劳寿命对数关系曲线如图10所示。

表5各个模型的疲劳寿命

Tab.5Fatiguelifeofmodels

R/?

m

Kt

AtJeqv/MPa

N（95%）/cycle

N（99.9%）/cycle

10

1.33

133

6211419

15

1.52

152

2065861

20

1.72

172

3274377

927460

25

1.91

191

1254313

514568—

30

2.10

210

629256

316164

图10粗糙度FZ与缺口的疲劳寿命关系

Fig.10Roughnessvsfatiguelifeofnotch

从图10可以看到，粗糙度艮增大，缺口的疲劳寿命减少。

存活率一定时，表面粗糙Rz与缺口的对数疲劳寿命之间有二次曲线拟合关系式

2

logN=11.625-0.378Rz0.00615Rz95%S（12）

2

logN=7.959-0.135Rz0.00179Rz99.9%S/（13）

式（12）和（13）的相关系数均在0.98以上。

对连续椭圆多微缺口化的表

面纹理，将Rz=23um（弋入式（12）和式（13），可得到缺口在不同存活率下的疲

劳寿命分别为1528797（95%S），632280（99.9%S）。

与Rz=23?

m时采用上面相同算法所得到的不同存活率下缺口的疲劳寿命1643383（95%S）、649575（99.9%S）进行对比，其相对误差分别为7%3%相对误差的平均值为5%

可以看到，用二次曲线来拟合缺口的表面粗糙度与缺口的对数疲劳寿命关系，其相对误差很小，通过验证平均值在5%左右。

因此，用二次曲线来拟合表面粗

糙度与缺口的对数疲劳寿命关系比较合适。

即可以认为，表面粗糙度与缺口的对

数疲劳寿命呈二次曲线关系。

5、结论

通过把平板表面纹理形貌特点连续相邻多微缺口化，采用有限单元方法计算得到不同表面粗糙度下的表面应力集中系数Kt，进而对缺口疲劳寿命进行了具有一定存活率的估算。

最后模拟分析有以下结果：

1）•多微缺口可以降低单缺口存在时的应力集中系数，且外沿缺口应力集中

程度高。

一定范围内，缺口半径保持不变时，缺口中心距的增大，K逐渐变大,

超过该范围即有Kt达到单缺口水平而不再受影响。

2）.缺口数目增加时，Kt减小，当数目增大到一定程度后K停止变化；缺口中心间距越大，缺口的数目增加对K的影响将减小。

3）•建立了表面粗糙度与表面应力集中系数之间的经验公式（8）。

4）•疲劳寿命随表面粗糙度增大而降低。

在相同存活率下，表面粗糙度与表面微缺口的对数疲劳寿命呈二次曲线关系。

6小结与展望

通过本次制作的论文，我深刻的发现自己的能力水平相对比较低，还有很多

不足之处，在以后的学习中，我应该更加的努力学习，打好基础，做好自己应该做的。

参考文献

1.吴松青.表面粗糙度应用指南•北京：

机械工业出版社，1990

2.张东初，裴旭明.加工工艺对表面粗糙度及疲劳寿命的影响.中国机械工

程,2003

3.郁大照，陈跃良，多缺口应力集中系数有限元研究.强度与环境,2002,

4.庄茁，张帆，岑松，等.ABAQUS非线性有限元分析与实例.北京：

科学出版

社.2003

5.魏建锋,郑修麟,吕宝桐.LY12CZ铝合金切口件的P—S—N曲线表达式及寿命估算.机械工程材料，1997

14级车辆工程一班

谷玉乐

2016年5月29日