基于sines准则的钢轨疲劳分析

《铁道机车车辆强度设计方法》基于sines准那么的钢轨疲劳分析说明书学院：机械工程学院专业：机械工程姓名：XX学号：2023210227指导老师：米彩盈日期：2023年5月11日钢轨疲劳概况在工程结构和机械设备中，疲劳的现象极为广泛，它普遍存在于每个承受反复荷载的结构中。疲劳潜藏于机械设备运转的过程中，它常常会造成加剧零部件的磨损和断裂，从而使机械零部件失效，使机械设备不能正常工作，甚至造成生命危险，因此疲劳特性是人们一直非常关注的问题。钢轨是轨道的主要部件，它是轮轨接触中的直接承重局部，引导机车车辆的车轮前进，承受车轮巨大的压力，并将所承受的荷载传布于轨枕、道床及路基。同时，钢轨必须为车轮提供连续、平顺和阻力最小的滚动外表。在电气化铁道或自动闭塞区段，钢轨还可兼作轨道电路之用。因此，钢轨的状态和使用寿命将直接影响铁路运输。所以，国内外纷纷投入大量科研力量进行理论分析和实验，研究在轮轨接触中钢轨的动力响应特征，并根据动力响应结果，预测钢轨寿命，指导工务部门对钢轨进行何时打磨、何时换轨等。基于abaqus6.10对钢轨疲劳应力的分析2.1钢轨的建模为使钢轨轧制冷却均匀，要求轨头、轨腰及轨底的面积分配，有一个较适宜的比例。根据我国对75kg/m钢轨断面尺寸在AutoCAD中绘制钢轨的截面如下列图1。根据已绘好的截面导入Pro/E中，根据0.65m的轨枕距装配如下列图2。图1钢轨截面图图2钢轨装配图2.2钢轨模型参数为了发挥钢轨的性能，它应具有足够的刚度，以抵抗由动荷载引起的弹性挠曲变形；具有一定的韧性，以防止在动荷载作用下，发生折断或损坏；具有足够的硬度，以防止车轮压陷或磨耗太快；其顶面应具有一定的粗糙度，以利实现机车的黏着牵引与制动力。我国用于轧制钢轨的材料有U71、U74、U71Cu、U71Mn、U71MnSi等，本次模拟疲劳应力选用钢轨材料为U71Mn，其弹性模量E＝2.145E5MPa，伸长率δ＝8%－10%，泊松比u＝0.3，密度ρ＝78kg/m3。2.3abaqus仿真分析钢轨模型使用一对轨枕进行边界模拟约束，因此轨枕距中心截面为钢轨受力的危险截面，模拟CRH2车型中24t轴重荷载情况下，对建立好的模型选用C3D8R六面体单元进行网格划分，并对钢轨进行四个工况的受力进行模拟分析：①轨枕距中心截面正中心节点受力；②轨枕距中心截面顶部边缘节点受力；③轨枕距中心截面线节点受力；④轨枕距中心截面面节点受力，分别如下列图3所示。将写好的文件导入工作站进行分析计算，导出计算结果如下列图4所示，并导出各个节点的MISESMAX、S.Mises、S.Max、S.Mid、S.Min应力文件，以便进行应力分析。②④③①②④③①图3钢轨四种工况下的受力④③①②④③①②图4abaqus分析结果图三．基于matlab2023a对abaqus6.10仿真分析结果进行分析3.1钢轨材料的力学性能钢轨材料U71Mn的力学性能反映了钢轨的工作能力，其根本力学性能有，强度极限δb＝964MPa，屈服极限δs＝552MPa，由材料力学性能经验公式可知，材料对称循环弯曲疲劳极限为材料强度极限的一半，因此可知钢轨材料的对称循环弯曲疲劳极限δ－1N=482MPa。塑性变形后碳钢扭转疲劳极限与弯曲疲劳极限的关系的相关知识，15号碳钢扭转与弯曲疲劳极限的比在0.53－0.57之间，45号碳钢扭转与弯曲疲劳极限的比为0.57－0.58，且有含碳量高的材料其疲劳极限比也较高。因此利用插值计算的方法，可以计算出钢轨材料U71Mn含碳量w(C)＝0.65－0.77的扭转疲劳极限ζ－1＝288MPa。3.2matlab程序编制及验证3.2.1修正的Goodman曲线由上述的参数可以对钢轨各节点的应力进行编程，并对其与Haigh图形式的修正的Goodman曲线比照，Matlab程序如附录一。其比照图形如图5所示。由图5可知，各个不同的受力工况下，钢轨各节点的应力均在Goodman曲线内部，因此钢轨在不同工况下，其危险截面的受力均在其材料力学性能内，其工作状态平安。图5钢轨应力与Haigh－Goodman的比拟3.2.2sines准那么sines准那么是最早提出并被普遍接受的高周多轴疲劳准那么，其内容如下：其中：①为第二应力偏张量不变量,＝;②;③;④,=。δb——材料的抗拉强度极限；δ－1——材料的旋转弯曲疲劳极限；Ζ－1——材料的扭转疲劳极限。由上述参数及计算公式，对钢轨各节点的应力进行matlab编程，并将其与sines准那么图线进行比照，其matlab程序如附录二。其比照图形如图6所示。由图6可知，各个不同的受力工况下，钢轨各节点的各种应力均在sines准那么图线的内部，因此钢轨在不同工况下，其危险截面的受力均在其材料力学性能内，其工作状态平安。图6钢轨应力与sines图线的比拟总结经过对CRH2型钢轨24t荷载情况下，模拟钢轨在不同工况受力的情况下，对其进行了abaqus仿真和基于Haigh－Goodman曲线和sines准那么下进行了matlab编程的理论分析。分析结果均说明，U71Mn的钢轨材料是满足理论的运行条件和要求，对实际铁路机车车辆的运行具有极大的参考和指导意义。同时，也说明了铁路机车车辆的开展必须要有足够的理论支撑，理论知识的不断完善和进步将毫无疑问地推动了铁路机车车辆向前快速开展。参考文献[1]杨宝峰,于春华.铁路轨枕现状及开展[J].铁道工程学报,2007,101(2):37-40.[2]肖生灵.我国铁路轨枕使用现状及轨枕材料开展趋势的研究[J].森林工程,2005,21(4):50-52.[3]SinesG.Behaviorofmetalsundercomplexstress.In:sinesG.WaismanJL,eds,MetalFatigue,NewYork:McGraw-Hill,1959,145-169.[4]HENRYOttenFuchs,RALPHIvanStephen.MetalFatigueinEngineering[M].NewYork:Wiley-Interscience,1980.[5]时新红,张建宇,鲍蕊等.材料多轴上下周疲劳失效准那么的研究进展[J].机械强度,2023,30(3):515－521.[6]苗德华.塑性变形后碳钢扭转疲劳极限与弯曲疲劳极限的关系[J].天津职业技术师范大学学报,2023,21(3):1-4.[7]韦树军.高速铁路钢轨的疲劳寿命研究[D].成都：西南交通大学,2006.[8]刘月明,李建勇,蔡永林等.钢轨打磨技术现状和开展趋势[J].中国铁道科学,2023,35(1):30－35.[9]刘启跃.钢轨的安定状态研究[J].西南交通大学学报,1995,30(4):466-471.[10]刘建锋,魏庆朝.高速铁路钢轨疲劳寿命[J].铁道工程学报,2000,66(2):30-33.[11]王元丰,梁亚平.高性能混凝土的弹性模量与泊松比[J].北方交通大学学报,2004,28(1):5-7.[12]赵雪芹.钢轨接触疲劳裂纹形成机理研究[D].成都：西南交通大学,2007.附录一〔1〕功能函数：求应力幅和平均应力function[y_a,y_m]=myfunction_pilao(x,y)y_a=zeros(1,12);y_m=zeros(1,12);fork=1:12y_a(k)=(x(k)-y(k))/2;enddisp('应力幅矩阵y_a＝');disp(y_a);fork=1:12y_m(k)=(x(k)+y(k))/2;enddisp('平均应力矩阵y_m＝');disp(y_m);〔2〕绘制Haigh-Goodman曲线并检验%绘制Haigh-Goodman曲线并检验:%其中y_b=964MPa,y_s=552MPa,y_n=482MPax=[0552];x1=-x;y=[5520];y_b=[-964964;9640];plot(x,y,'-k',...x1,y,'-k',...y_b(:,1),y_b(:,2),'-k');axis([-600,1000,0,600]);xlabel('y_m/MPa');ylabel('y_a/MPa');disp('工况一：');y_max=[187.64487.9213.6506137.842-36.3812...155.369155.369152.72155.36936.7922...63.763336.7924];y_min=[-87.8606-118.752-20.7267-223.32...-253.108-223.32-282.897-282.897...-282.897-50.1325-43.1023-50.1326];[y_a,y_m]=myfunction_pilao(y_max,y_min);holdonscatter(y_m,y_a,'s','markerfacecolor','b',...'markeredgecolor','b');disp('工况二：');y_max=[38.853349.797610.122961.7767...72.997961.777243.4437511.858E-03...43.4431-42.42-71.8975E-03-42.4193];y_min=[-57.8156-66.4293-14.6476-64.6593...-79.4268-64.6595-130.38-275.939...-130.38-421.497-265.026-421.498];[y_a,y_m]=myfunction_pilao(y_max,y_min);holdonscatter(y_m,y_a,'d','markerfacecolor','g',...'markeredgecolor','g');disp('工况三：');y_max=[7.65017-600.245e-037.6501416.6076...1.8125916.607739.768411.1223...4.7200517.524412.450317.5245];y_min=[-41.3266-51.317-41.3266-31.3219...-93.705-110.292-24.0324-105.924...-110.292-101.555-86.851-101.555];[y_a,y_m]=myfunction_pilao(y_max,y_min);holdonscatter(y_m,y_a,'p','markerfacecolor','r',...'markeredgecolor','r');disp('工况四：');y_max=[5.671185.624998.481277.36345...-3.69083-3.69456-10.6576-10.6372...8.335578.36779-9.91144-9.94554];y_min=[-34.5103-34.5401-31.4678-37.4927...-69.5824-69.4646-90.2238-90.4277...-37.2167-35.2162-77.3741-77.1295];[y_a,y_m]=myfunction_pilao(y_max,y_min);holdonscatter(y_m,y_a,'o','markerfacecolor','k',...'markeredgecolor','k');附录二sines准那么疲劳强度评定%疲劳强度评定%sines准那么:其中t=288MPa,y_n=482MPa,y_b=964MPaclear;clc;t=288;y_n=482;y_b=964;k=(3*t*(y_b+y_n))/(y_b*y_n)-sqrt(6);x=-350:2000;y=-k*x+t;plot(x,y);xlabel('y_H,_m/MPa');ylabel('sqr.J_2,_a/MPa');axis([-350,2000,-150,500]);disp('工况一：');y_max=[187.64487.9213.6506137.842-36.3812...155.369155.369152.72155.36936.7922...63.763336.7924];y_mid=[76.015376.0153-558.285e-0332.0417...-122.09164.7092-147.633-70.7018...-147.633-16.4688-14.8813-16.4687];y_min=[-87.8606-118.752-20.7267-223.32...-253.108-223.32-282.897-282.897...-282.897-50.1325-43.1023-50.1326];y_mises=[142.003179.39822.8006159.826...189.165159.826144.198181.445...144.19877.995696.025877.9957];y_h=(y_max+y_mid+y_min)/3;aq_j=y_mises/(sqrt(3));holdonscatter(y_h,aq_j,'s','markerfacecolor','b',...'markeredgecolor','b');disp('工况二：');y_max=[38.853349.797610.122961.7767...72.997961.777243.4437511.858E-03...43.4431-42.42-71.8975E-03-42.4193];y_mid=[-18.6161-24.3711-3.95117-8.65852...-3.41798-8.65873-30.3961-64.6681...-30.3959-98.9401-43.6302-98.9402];y_min=[-57.8156-66.4293-14.6476-64.6593...-79.4268-64.6595-130.38-275.939...-130.38-421.497-265.026-421.498];y_mises=[84.4484102.03122.0141110.551...132.686110.552206.045252.659...206.045200.232246.157200.233];y_h=(y_max+y_mid+y_min)/3;aq_j=y_mises/(sqrt(3));holdonscatter(y_h,aq_j,'d','markerfacecolor','g',...'markeredgecolor','g');disp('工况三：');y_max=[7.65017-600.245e-037.6501416.6076...1.8125916.607739.768411.1223...4.7200517.524412.450317.5245];y_mid=[-4.82279-10.9099-4.82283-7.37709...-20.7407-7.377088.12969-35.1937...-23.2061-47.1813-37.9153-47.1813];y_min=[-41.3266-51.317-41.3266-31.3219...-93.705-110.292-24.0324-105.924...-110.292-101.555-86.851-101.555];y_mises=[44.206346.466144.206368.9947...86.478868.994779.8837103.579...79.883769.291