基于珠光体与贝氏体体积比的轮轨接触力学性能评估方法

本发明属于轮轨接触力学与材料性能评估，特别涉及一种轮轨接触力学性能评估方法，具体是一种基于珠光体与贝氏体体积比的轮轨接触力学性能评估方法。

背景技术：

1、珠光体和贝氏体是钢轨用钢中最常见的两类微观组织。抗拉强度一般在800~1200mpa范围内，塑性较好。贝氏体由过饱和铁素体和碳化物构成，硬度通常高于珠光体，但断后伸长率和冲击韧性相对偏低。通过调整热处理过程中的冷却速度和等温温度，可以改变这两相在钢中的相对含量，从而获得不同力学性能的材料。珠光体与贝氏体的体积配比不仅决定了材料的宏观弹性模量和屈服强度，还直接影响其在高应力循环加载条件下的变形行为和损伤演化机制。珠光体的层片状结构赋予材料较高的抗拉强度和良好的抗磨损能力，而贝氏体的细小板条结构则提供了较高的硬度和一定的韧性储备。两相的力学性能互补性使得珠光体与贝氏体双相钢在硬度、强度和韧性之间具有良好的匹配优势，但这种优势的发挥高度依赖于两相的相对含量及其空间分布状态。在实际工程中，钢轨材料的微观组织往往随着热处理工艺参数（如冷却速率、等温温度）的波动而发生变化，这种变化会直接反映在两相体积比上，进而影响轮轨接触面的力学响应。因此，弄清不同体积比下材料力学性能的变化规律，对于稳定控制钢轨质量具有重要意义。

2、轮轨接触面承受极高的法向载荷和切向摩擦力，接触应力可达数gpa，足以引发塑性变形和疲劳损伤。材料的微观组织构成直接影响其在高应力接触条件下的表现，因此弄清珠光体和贝氏体的体积配比如何影响轮轨接触力学行为，对于钢轨成分设计和热处理工艺制定有实际意义。特别是在重载铁路和高速铁路运营条件下，轮轨接触应力水平进一步提高，微观组织对接触力学行为的调控作用更加显著，系统掌握二者之间的定量关系具有重要的工程应用价值。重载铁路的轴重已从25吨提升至30吨甚至更高，接触应力随之增大，对钢轨材料的微观组织调控提出了更高要求。不同轴重条件下，接触应力的分布状态和幅值存在差异，材料微观组织对接触响应的影响程度也可能随之变化。

3、对于轮轨的材料微观组织参数与轮轨宏观接触力学性能不能快速关联而进行评估的问题，现有技术中存在以下不足：

4、第一，对于珠光体与贝氏体两相体系在0%到100%全成分范围内的系统研究尚不充分。多数工作集中在某几个特定成分点的实验测试，缺少覆盖整个成分区间的连续变化数据，不同体积比下接触力、应力分布和滑移的演变规律尚不完整。不同体积比下材料的等效弹性模量、泊松比等基础力学参数的变化规律尚缺乏系统表征，这直接制约了轮轨接触有限元模拟中材料本构参数的合理选取。此外，等效弹性模量作为有限元模拟中最重要的输入参数之一，其取值直接影响接触应力、应变和变形场的计算精度。然而，现有研究中对珠光体与贝氏体双相材料等效弹性模量随体积比变化的定量关系尚无系统报道，voigt混合法则在该体系中的适用性也缺乏充分验证。仅依靠经验取值或单点实验数据外推，难以保证不同体积比下接触计算结果的可靠性。

5、第二，在轮轨接触有限元模拟中，钢轨材料通常被处理成均质体，直接取一种固定成分的力学参数，未充分考虑实际钢轨中两相比例变化对接触计算结果的定量影响。这导致仿真选材时缺乏依据。对于不同牌号或不同热处理批次的钢轨，其珠光体含量可能存在较大差异，但现有仿真模型对此差异缺乏有效的参数化描述手段，通常只能采用同一组材料参数进行计算，无法反映实际组织状态对接触力学行为的影响。

6、第三，两相的力学性能差异会在相界面处产生应变不协调，这种局部效应如何映射到轮轨接触界面的宏观响应如接触力大小、滑移方向和幅值等，尚无明确结论。两相界面处的应力应变集中程度与体积配比密切相关，但现有研究多局限于定性描述，缺少能够定量描述这一关系的计算模型。在接触载荷作用下，软硬两相的变形差异可能导致界面脱粘或微裂纹萌生，这种微观损伤如何影响宏观接触性能，目前尚缺乏系统的定量研究。

7、第四，现有研究缺乏从微观本构参数到宏观接触响应的系统评估方法。材料微观组织参数获取后，如何快速、准确地评估其对轮轨接触力学性能的影响，缺少一套完整的仿真评估流程，限制了轮轨材料微观组织设计与宏观服役性能之间的有效衔接。由于缺少标准化的评估流程，不同研究者采用的本构模型、边界条件和评价指标各不相同，结果之间缺乏可比性。建立一套标准化的评估流程，可使不同材料体系、不同研究条件下的结果具有可比性，对于推动轮轨接触力学领域的发展具有重要意义。

8、相关专利文献：cn121682157a公开了一种基于逻辑回归的轮轨接触状态评估方法，基于轮轨接触理论，选取接触点分布百分比数、集中度复合指数及3mm等效锥度以多维度参数更全面、精准地表征轮轨接触状态，避免因参数片面性导致的评估偏差，使等级评估结果更贴合实际轮轨匹配性能，且采用逻辑回归模型作为评估核心算法，模型输入与输出的映射关系可通过回归系数量化解释，模型可解释性强，适配工程实践需求，且模型阈值可根据不同线路运维需求灵活调整，规避潜在安全风险，适配多样化的铁路运营场景。

9、以上技术并未解决已有技术中，在获取轮轨的材料微观组织参数后，轮轨的材料微观组织参数与轮轨宏观接触力学性能不能快速关联而进行评估的问题。

技术实现思路

1、本发明的目的在于提供一种基于珠光体与贝氏体体积比的轮轨接触力学性能评估方法，在获取轮轨的材料微观组织参数后，能快速、准确地评估所述材料微观组织对轮轨接触力学性能的影响，以解决轮轨的材料微观组织参数与轮轨宏观接触力学性能不能快速关联而进行评估的问题。

2、为解决上述技术问题，本发明采用的技术方案如下：

3、一种基于珠光体与贝氏体体积比的轮轨接触力学性能评估方法，其技术方案在于，它包括以下步骤：

4、s1：针对轮轨材料中珠光体与贝氏体，设置多组珠光体与贝氏体体积比。

5、s2：确定材料参数。

6、s3：根据voigt等应变混合法则，计算各组所述体积比下的等效弹性参数，等效弹性参数包括等效弹性模量和等效泊松比。上述混合法则基于等应变假设，适用于两相界面结合良好、变形协调性较好的双相材料体系，是计算双相材料等效弹性性能的方法。

7、s4：建立哑铃型单轴拉伸试件有限元第一模型，进行拉伸仿真验证：将步骤s3得到的各组所述体积比下的等效弹性参数分别赋予所述第一模型，获取各组所述体积比下的应力-应变曲线，验证所述混合法则在弹性阶段的预测精度；提取所述应力-应变曲线中的数据进行线性拟合，得到各组所述体积比下的等效弹性模量仿真值，将每个所述等效弹性模量仿真值与步骤s3中相对应的等效弹性模量理论值比对，相对偏差≤x时进入下一步；否则，返回检查参数与模型设置，重新验证。

8、s5：建立轮轨接触三维有限元模型；将步骤s4中比对后的各组所述体积比下的等效弹性模量仿真值（各组等效弹性参数）及按混合法则线性插值确定的塑性参数赋给所述模型中的钢轨材料；进行轮轨接触仿真。

9、s6：提取轮轨接触界面上的力学响应参数，得到各力学响应参数随珠光体体积分数的变化特征。

10、s7：根据各力学响应参数随珠光体体积分数的变化特征，判定轮轨接触力学性能演变阶段：当≤50%时性能稳定；当60%≤≤80%时性能发生突变；当≥85%时性能恢复稳定；为珠光体体积分数。

11、上述技术方案中，优选的技术方案可以是，步骤s1中，设置11组珠光体与贝氏体体积比，这11组珠光体与贝氏体体积比为：0∶10，1∶9，2∶8，3∶7，4∶6，5∶5，6∶4，7∶3，8∶2，9∶1，10∶0。上述11组珠光体与贝氏体体积比覆盖了从全贝氏体到全珠光体的完整成分区间，步长为10%，足以反映两相比例连续变化对材料力学性能的影响趋势。11组珠光体与贝氏体体积比数据构成了全成分扫描，可系统识别性能突变区间，避免了因步长过大而遗漏临界体积比。

12、步骤s2中，所述材料参数包括珠光体弹性模量、珠光体泊松比、贝氏体弹性模量、贝氏体泊松比。

13、所述材料参数的数值可以是：珠光体弹性模量=206.3gpa，珠光体泊松比=0.2885，贝氏体弹性模量=184.5gpa、贝氏体泊松比=0.314。材料参数还可以包括珠光体的屈服应力与塑性应变关系曲线中的数值、贝氏体的屈服应力与塑性应变关系曲线中的数值。

14、步骤s3中，等效弹性模量和等效泊松比的计算公式为：

15、 （1），

16、 （2），

17、其中，为珠光体体积分数，为贝氏体体积分数，且。

18、步骤s4中，哑铃型单轴拉伸试件（即哑铃型试件）的总长度为300mm，标距段的长度为80mm，夹持端的直径为60mm，标距段的直径为35mm，圆角过渡区域处的圆角半径为106mm，采用c3d8r八节点六面体线性减缩积分单元进行网格划分，标距段中心区域局部细化，单元总数为10000-12000个，单元总数可选择12000个。提取所述应力-应变曲线中的数据进行线性拟合、得到各组所述体积比下的等效弹性模量仿真值的方法是，采用位移加载方式，将所述试件的下端面固定，在所述试件上端面沿轴向施加总位移为12mm、加载速率为2mm/min的拉伸载荷，提取所述应力-应变曲线中标距段中心节点的应力-应变数据进行线性拟合，得到各组所述体积比下的等效弹性模量仿真值；所述标距段中心节点应力峰值为37.9~42.4mpa，仅占屈服强度的6.5%~10.8%，所有工况处于线弹性变形阶段。哑铃型试件是材料单轴拉伸试验的标准形式，其标距段为均匀应力区，中心节点应力-应变数据能够准确反映材料的本构响应。

19、步骤s4中，采用位移加载方式加载时，加载过程中动能与内能之比始终≤5%，满足准静态加载条件。步骤s4中，x=0.01%。

20、上述技术方案中，优选的技术方案还可以是，步骤s5中，所述轮轨接触三维有限元模型包括轮对、构架、外轨、内轨、浮置板、扣件、一系悬挂、接地弹簧和道床，轮对由外侧车轮、内侧车轮和车轴组成，外轨、内轨即为两根钢轨；轮轨是由车轮（外侧车轮、内侧车轮）和钢轨（外轨、内轨）组成的一个系统。进行轮轨接触仿真的方法是，以外侧车轮的表面、内侧车轮的表面为主面，以外轨的表面、内轨的表面为从面，采用有限滑移公式，法向接触行为为贴紧，即采用“硬”接触，切向行为采用罚函数摩擦公式，摩擦系数按典型钢-钢接触特性取0.35，进行轮轨接触仿真。

21、所述轮轨接触三维有限元模型完整反映了轮轨系统的实际传力路径：轮轨接触力通过一系悬挂传递至转向架的构架，再经二系悬挂传递至车体，扣件和浮置板则将轮轨力传递至道床。将各组等效弹塑性参数赋给钢轨材料，可模拟不同微观组织下轮轨接触响应的差异，轮轨接触仿真中，钢轨材料的本构参数直接影响接触斑内的应力分布和塑性应变积累。将不同体积比对应的等效弹塑性参数分别赋给钢轨材料，可定量评估微观组织差异对轮轨接触宏观响应的影响程度。摩擦系数的取值参考典型钢-钢接触的实验数据，确保仿真结果与工程实际具有可比性。

22、步骤s5中，所述一系悬挂的参数为：垂向刚度是1.30×106n/m，横向刚度是1.10×106n/m，纵向刚度是1.10×106n/m，垂向阻尼是5.00×103n·s/m；扣件本体的参数为：垂向刚度是4.07×107n/m，横向刚度是8.79×106n/m，纵向刚度是8.79×106n/m，垂向阻尼是9.898×103n·s/m，横向阻尼是1.928×103n·s/m，纵向阻尼是1.928×103n·s/m，间距是0.60m；浮置板中的钢弹簧隔振器的参数为：垂向刚度是5.33×106n/m，横向刚度是4.56×106n/m，纵向刚度是4.56×106n/m，垂向阻尼是1.50×104n·s/m，横向阻尼是1.00×104n·s/m，纵向阻尼是1.00×104n·s/m，间距是1.20m；转向架的轴距是2.50m。

23、步骤s5中，所述模型内钢轨的材料为弹塑性本构，各向同性硬化。

24、步骤s6中，所述力学响应参数包括接触法向力cnormf、接触法向应力cpress、横向切应力cshear1、纵向切应力cshear2、合成切向力cshearf、横向滑移cslip1和纵向滑移cslip2； 各组所述体积比下的等效摩擦系数=cshearfmax/cnormfmax，均处于0.350~0.351范围内。上述参数涵盖了轮轨接触界面的法向作用、切向摩擦和相对滑移三个方面的力学行为，构成了完整的接触力学评价指标体系。

25、上述技术方案中，优选的技术方案还可以是，步骤s7中，所述性能发生突变的微观机制为：珠光体体积分数=70%时，硬质相珠光体趋于形成连续承力骨架的临界阈值，珠光体与贝氏体两相间变形，不协调加剧，导致接触力学性能突变。

26、步骤s7中，所述性能发生突变的具体表现为：在珠光体体积分数=70%处，接触法向力cnormf降低，由1.178×104n降至6.524×103n，降幅44.6%；接触法向应力cpress降低，由15.44 mpa降至8.551mpa，降幅44.6%；横向切应力cshear1增加，由±0.11mpa增至±0.17mpa，增幅50%；纵向切应力cshear2的绝对值降低，由5.40mpa降至2.99mpa，降幅44.7%；合成切向力cshearf降低，由4.128×103n降至2.287×103n，降幅44.6%；横向滑移cslip1的幅值缩小，由0.075mm缩至0.050mm，缩减33%；纵向滑移cslip2由正向主导（+0.295~0.000mm）反转为负向主导（+0.061~-0.273 mm）。

27、所述的基于珠光体与贝氏体体积比的轮轨接触力学性能评估方法，它最好还包括以下步骤：

28、s8：在步骤s7之后，确定设计准则输出步骤，即当工程应用要求稳定的接触力学性能时，选择≤50%或≥85%，避开60%≤≤80%的临界敏感区间，为珠光体体积分数。

29、本发明通过abaqus有限元软件实现，采用python脚本进行参数化建模与批量仿真计算。

30、本发明针对轮轨材料中珠光体与贝氏体的体积配比与接触力学性能之间缺乏定量评估手段的问题，以多组（最好是11组）不同珠光体与贝氏体体积比（从0∶10至10∶0）为研究对象，首先采用voigt等应变混合法则计算各组等效弹性模量与等效泊松比，建立哑铃型单轴拉伸试件有限元模型进行拉伸仿真，获取各体积比下应力-应变曲线，验证混合法则在弹性阶段的预测精度；然后将验证后的等效弹塑性参数分别赋给轮轨接触三维有限元模型，进行轮轨接触仿真，提取接触面上的法向力、切向力及滑移量作为评价指标。分析结果表明：各体积比下应力-应变曲线呈高度线性，等效弹性模量随珠光体含量增加由184.5gpa线性升至206.3gpa；接触法向力与法向应力随珠光体含量呈非单调演变，当珠光体体积分数达到70%时出现临界突变，接触法向力骤降44.6%，横向切应力增加50%、纵向切应力绝对值下降44.7%、横向滑移幅值减小约33%、纵向滑移方向发生反转；各体积比下等效摩擦系数稳定在0.350~0.351之间，不受微观组织比例影响。本发明明确了珠光体体积分数60%~80%为轮轨接触力学性能临界敏感区间，材料设计应使珠光体含量即珠光体体积分数≤50%或≥85%，以避开所述的不稳定区间。本发明无需大量实验测试，仅通过有限元仿真即可实现轮轨材料微观组织参数与宏观接触力学性能的快速关联评估，可广泛应用于轮轨材料研发阶段的成分筛选与热处理工艺优化。该方法同样适用于不同钢种、不同热处理状态下的轮轨材料接触力学性能对比评估，为钢轨材料设计提供了一套通用的仿真分析工具。

31、本发明提出了一套从单相参数输入到接触响应输出的完整评估流程，通过有限元手段建立起体积比与接触性能之间的对应关系；能识别接触力学性能发生突变的临界体积比区间；阐明了微观组织体积比变化导致宏观接触力学性能突变的微观机制。采用本发明的上述技术方案，所产生的有益效果在于：

32、①建立了珠光体与贝氏体体积比与轮轨接触力学性能之间的定量评估方法。本发明设置了多组（优选11组）珠光体与贝氏体体积比，覆盖从0%到100%的全成分区间。在此基础上，系统获得了接触法向力、法向应力、切向力及滑移等力学响应参数随体积比的完整演变规律。该结果弥补了双相材料微观组织与宏观接触性能之间定量关系研究不确定的缺陷。

33、②识别了珠光体体积分数为60%~80%是轮轨接触力学性能的临界敏感区间。本发明提出珠光体含量在70%附近时，接触法向力骤降44.6%、横向切应力增大50%、纵向滑移方向发生反转等一系列临界现象，明确了该区间为材料性能不稳定的敏感区域，为轮轨材料的微观组织目标设计提供了明确的避开准则。

34、③建立了从单轴拉伸本构表征到轮轨接触响应的完整仿真评估流程。本发明通过哑铃型单轴拉伸试件拉伸仿真，获取了各组所述体积比下的弹性本构参数，并将弹塑性参数赋予轮轨接触模型进行后续分析，实现了从材料本构到宏观接触响应的全链条仿真评估，无需依赖大量实验测试，计算效率高、可重复性强。该流程中的每一个环节都有明确的输入输出和判定标准，可操作性强，不同研究人员采用该方法得出的结果具有良好的一致性。该评估流程采用有限元仿真手段，可大幅降低实验成本，缩短材料评估周期。

35、④揭示了微观组织体积比变化导致宏观接触力学性能突变的微观机制。本发明从两相变形不协调的角度，阐明了在珠光体体积分数=70%处力学性能突变的本质原因：硬质相珠光体形成贯穿连续网络所引发的软硬相界面应力集中与变形不协调。

36、⑤验证了等效摩擦系数与微观组织成分无关的特性。本发明提出等效摩擦系数在所有体积比下均稳定在0.350~0.351之间，确认摩擦系数由接触界面特性决定。

37、⑥验证了voigt等应变混合法则在珠光体/贝氏体双相材料弹性模量预测中的高精度（相对偏差≤0.01%），为有限元模拟中双相材料本构参数的选取提供了科学依据。

38、综上所述，本发明提供了一种基于珠光体与贝氏体体积比的轮轨接触力学性能评估方法，在获取轮轨的材料微观组织参数后，能快速、准确地评估所述材料微观组织对轮轨接触力学性能的影响，解决了轮轨的材料微观组织参数与轮轨宏观接触力学性能不能快速关联而进行评估的问题。本发明适用于轮轨材料研发阶段的成分筛选和有限元建模中的本构参数选取，还适用于不同热处理工艺下轮轨材料接触力学性能的对比评估，为钢轨材料的微观组织设计与服役性能优化之间的关联研究提供了一种方法，该方法还可用于分析不同珠光体与贝氏体体积比下轮轨接触疲劳寿命的差异，为延长钢轨服役周期提供依据。