- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
黄河交通学院学课论文(论文)题目:表面粗糙度对表面应力集中系数和疲劳寿命影响分析完成人:谷玉乐班级:14级车辆工程一班学制:4 年专业:汽车工程学院车辆工程指导教师:xxx完成日期:2016年5月29日表面粗糙度对表面应力集中系数和疲劳寿命影响分析摘要 把平板的表面形貌简化为半椭圆形微缺口,采用有限元法对不同表面粗糙度下应力场进行分析,以得到不同表面粗糙度Rz时平板的表面应力集中系数, 并预测出不同表面粗糙度的平板在不同存活率下的疲劳寿命,讨论表面粗糙度对疲劳寿命的影响。采用回归方法建立了表面应力集中系数与微缺口中心间距、表面粗糙度之间的经验公式。结果表明,在相同存活率下,表面粗糙度与平板的对数疲劳寿命呈二次曲线关系。关键词 表面粗糙度 微缺口 应力集中系数 疲劳寿命1 、引言表面粗糙度是反映零件表面微观几何形状误差的一个重要指标[1] 7-8。由断裂力学可知,表面粗糙度值愈大,表面的沟痕愈深,纹底半径愈小,应力集中越严重,抗疲劳破坏的能力就愈差。因此表面粗糙度增大,会降低零件的疲劳强度[2]。这一影响引起了广泛的关注。在实验方面,Noll比较了低碳合金钢的表面分别在研磨、机加、热滚压和锻造这四种表面加工条件下,其应力-寿命关系曲线。发现在高应力水平下,所有试件的寿命相近,在低应力区,寿命差异十分显著。在应力幅值接近屈服强度的80%时,研磨试件的寿命是锻造试件寿命的20倍[3]。在数值计算方面,Andrews将表面粗糙度作为微观缺口,研究缺口对疲劳寿命的影响[4]。文献[5]采用试验测试和有限元的方法,建立实际表面形貌的有限元网格,将表面犁沟视为微观缺口,研究裂纹萌生部位与缺口应力集中系和应力场的关系。由此可见,要确定表面粗糙度对疲劳寿命的影响,首先要建立表面粗糙度与表面应力集中系数的定量关系。本文对不同表面粗糙度下的模型进行唯象统计分析,旨在建立微缺口中心间距、粗糙度Rz与应力集中系数Kt之间的定量关系,进而分析一定存活率下,缺口疲劳寿命与表面粗糙度的关系。2、计算模型图1模型尺寸及受载示意图Fig.1Dimension and load of models分析计算用商用有限元分析软件ABAQUS来完成。表面形貌被简化为半椭圆形微缺口。如下图1所示:在边长L为1mm平板的上边缘中间处有一组长短半径分别为a和b的半椭圆缺口;短半径b即为缺口深度,长半径a变化范围为10至50?m,b/a在0.2~1范围内变化。多缺口中心间距d(下文简称间距)取1~5倍缺口宽度2a。结构两边受均布拉伸载荷P,大小为100MPa。材料为LY12CZ铝合金,其弹性模量E=68 GPa,泊松比μ= 0.3。对模型进行有限元网格划分时,先分区,再用不同划分方法和单元大小来刻画不同部分,以实现计算速度和质量的统一。划分时远离缺口部分选择二次减缩积分单元CPS8R,缺口附近的部分采用尺寸为1微米的二次完全积分单元CPS8来进行计算。在缺口附近部分采用自由划分法,其他位置采用结构划分法;同时选择四边形单元[10]。a为20?m,b为10?m单微缺口的模型网格划分最后结果如图2所示,缺口部分细节如图3所示。图2 单微缺口网格划分图示图3单微缺口的局部细化网格Fig.2Finite element model Fig.3Local mesh ofmicro-notch3、应力集中系数受影响分析3.1半径对单微缺口应力集中系数的影响建立单微缺口长短半径变化的多个模型,得到缺口处应力集中系数。如图4长短半径为20?m,10?m的缺口部分Mises应力分布。在此用Mises应力作为局部最大应力即可较为准确的求出缺口处的应力集中系数值[8]。图4 单微缺口Mises应力云图分布(局部)Fig.4Mises stress counter of micro-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不同缺口长短半径比时的KtTab.1Kt with different ratio between depth and breadthb/a12/30.52/51/32/71/42/91/5a/?m1015102040253035404550Kt3.042.36
文档评论(0)