35CrMo 钢螺栓疲劳寿命预测方法研究 Fatigue life prediction method of 35CrMo alloy steel bolt 1 first-author 2 3 2 corresp 庞建超，博士，副研究员，从事材料疲劳与断裂研究；E-mail：jcpang@imr.ac.cn 庞建超，博士，副研究员，从事材料疲劳与断裂研究；E-mail：jcpang@imr.ac.cn jcpang@imr.ac.cn 2 1 2 南方电网调峰调频发电有限公司 检修试验分公司，广东 广州，511400 南方电网调峰调频发电有限公司 检修试验分公司，广东 广州，511400 Maintenance and Test Branch, China Southern Power Grid Power Generation Co. Ltd., Guangzhou 511400, China Maintenance and Test Branch, China Southern Power Grid Power Generation Co. Ltd., Guangzhou 511400, China 中国科学院金属研究所 师昌绪先进材料创新中心，辽宁 沈阳，110016 中国科学院金属研究所 师昌绪先进材料创新中心，辽宁 沈阳，110016 SHI Changxu Innovation Center for Advanced Materials, Institute of Metal Research, Chinese Academy of Sciences, Shenyang 110016, China SHI Changxu Innovation Center for Advanced Materials, Institute of Metal Research, Chinese Academy of Sciences, Shenyang 110016, China 沈阳航空航天大学 机电工程学院，辽宁 沈阳，110136 沈阳航空航天大学 机电工程学院，辽宁 沈阳，110136 School of Mechatronics Engineering, Shenyang Aerospace University, Shenyang 110136, China School of Mechatronics Engineering, Shenyang Aerospace University, Shenyang 110136, China 赵俊 针对大型螺栓的疲劳寿命预测问题，以35CrMo钢螺栓为研究对象，提出以应力集中系数为等效损伤参量的螺栓模拟构件模型，并进一步考虑疲劳综合修正系数，预测服役螺栓的疲劳寿命。研究结果表明：随着公称直径、螺纹圈数和牙底圆角的增加，螺纹应力集中系数呈现下降趋势，而螺距的增加使得应力集中系数逐渐增大；模拟构件的螺纹缺口模型与圆环缺口模型的应力集中系数差异很小，计算分析时可用圆环缺口模型代替螺纹缺口模型；综合修正系数方法预测的螺栓疲劳寿命曲线与有限元模拟的结果基本一致，验证了本研究所提出的基于模拟构件的疲劳寿命预测方法具有较高准确性及较好适用性。 For solving the problem of fatigue life prediction on large-scale bolt, 35CrMo steel bolt was selected as the research object, and a simulated bolt component was proposed by the stress concentration factor as equivalent damage parameter. Then fatigue life of service bolt was estimated based on the fatigue comprehensive correction factor. The results show the factor of stress concentration takes a drop trend with the increase of nominal diameter, screw number and bottom corner, while the increase of screw pitch makes stress concentration factor gradually rise. For equivalent models of screw thread and circular notch, the stress concentration factors are basically the same, and so the model of circular notch can be applied to replace the screw thread model in calculation and analysis. The curves of fatigue life obtained by comprehensive coefficient method and finite element simulation are in good agreement. This verifies the fatigue life prediction method based on simulated component proposed in this research has the higher accuracy and applicability. 35CrMo钢 螺栓模拟构件 应力集中系数 综合修正系数 疲劳寿命预测 35CrMo steel simulated bolt component stress concentration factor comprehensive correction factor fatigue life prediction 螺栓作为重要的紧固件之一，其连接的可靠性直接关系到设备的长期稳定运行。35CrMo钢作为典型的低合金高强钢，具有较高的静强度、冲击韧性以及疲劳极限，常被用于作为高强度螺栓的材料 [ 1 1 ] 。为了保证服役设备的安全可靠，装配时通常需要控制螺栓拧紧力矩来施加适当的轴向预紧力。由于服役过程中工况周期的往复，使得连接螺栓承受的拉伸载荷不断循环变化，导致其发生疲劳破坏 [ 2 2 - 4 4 2-4 ] 。 紧固螺栓，特别是大型构件(例如抽水蓄能机电站、船舶舰艇、重型汽车等)，一旦发生疲劳断裂，将严重影响设备的正常运行，并造成很大的安全隐患，甚至人员伤亡及恶劣社会影响 [ 5 5 ] ，故对预紧螺栓开展疲劳寿命进行评估十分重要。现阶段相关研究人员多利用数值仿真及理论公式，根据服役环境对螺栓的强度进行校核计算。ZHAO等 [ 6 6 ] 建立了构件螺栓的有限元模型，并利用获得的力学特性预测了其疲劳寿命。欧阳卿 [ 7 7 ] 模拟了螺栓的预紧力，综合考虑应力幅和应力比的影响估算了其疲劳寿命。谢里阳等 [ 8 8 - 9 9 8-9 ] 基于弹塑性有限元分析的局部应力法，研究了高强度螺栓构件的疲劳寿命。杜静等 [ 10 10 ] 通过引入“疲劳强度缩减因子”对高强度螺栓的疲劳寿命曲线进行修正，得出了塔筒环形法兰连接螺栓的疲劳寿命。练继建等 [ 11 11 ] 运用传统雨流计数法和疲劳损伤累积理论，确定了水轮机顶盖联接螺栓的疲劳寿命。TARAKANOV [ 12 12 ] 研究了腐蚀环境及循环载荷下螺栓的疲劳损伤，并提出腐蚀环境下螺栓耐久度的估算方法。MARCELO等 [ 13 13 ] 对不同冶金条件下螺栓疲劳特性进行了研究，发现热处理对螺栓疲劳强度有很大影响。刘胜祥等 [ 14 14 - 16 16 14-16 ] 对风力发电机的螺栓开展了疲劳强度特性分析，为整机的可靠性设计提供了参考方法。此外，特定条件下的疲劳实验是研究螺栓损伤行为的重要方法 [ 17 17 - 18 18 17-18 ] 。YILMAZ等 [ 19 19 ] 分析了高强度钢螺栓连接在典型疲劳加载模式下的疲劳性能。POOVAKAUD等 [ 20 20 ] 通过开发的简单模型，研究了高强度连接螺栓中接触表面的疲劳损伤机制。YANG等 [ 21 21 - 22 22 21-22 ] 开展了恒定应力幅疲劳试验，讨论了在不同应力范围下M22高强度螺栓的疲劳破坏特性。 上述研究为工程中螺栓的设计、计算及寿命评估提供了一定参考。然而，某些大型设备的螺栓规格巨大，公称直径可达100 mm以上，导致开展螺栓构件的疲劳寿命测试试验及数值仿真计算非常困难。为解决此类问题，在通过简单性能预测复杂性能及小样品预测大构件的研究思路 [ 23 23 - 26 26 23-26 ] 基础上，本文以35CrMo钢预紧螺栓为研究对象，提出了以应力集中系数为等效损伤参量的螺栓模拟构件模型，进一步基于疲劳综合修正系数开展螺栓疲劳寿命的等效评估，并利用有限元仿真疲劳寿命计算方法验证有效性。 1 35CrMo 显微组织与力学性能分析 1 抽水蓄能机组频繁起停及工况转换时，水轮机过流部件受到较大的水流作用力，尤其是顶盖作为一个大型的承压部件，其紧固螺栓承受着往复交变载荷作用，极易诱发疲劳损伤破坏。本文以工程中35CrMo钢螺栓为研究对象，研究服役工况下紧固螺栓的疲劳寿命，其剖切示意图如 图1 图1 所示。 (a) 正面视图；(b) 截面视图 35CrMo钢的微观结构特征SEM照片如 图2 图2 所示。从 图2 图2 可以看出，样品组织为回火索氏体，渗碳体颗粒极为细小，且较为均匀地分布在等轴 α 相晶粒内，原本马氏体板条特征已完全消失。 放大倍数：(a) 500；(b) 1 000 利用Instron 5982和GPS-100试验机，分别对35CrMo钢的拉伸性能和疲劳性能开展测试 [ 23 23 ] 。其中，依据GB/T 228—2010，拉伸试样尺寸规格为平行段30 mm、标距25 mm，拉伸速率为1×10 -3 s -1 ；疲劳实验采用正弦波形，实验频率为112~118 Hz，依据GB/T 24176—2009拟合疲劳寿命曲线，35CrMo钢拉伸性能和疲劳性能如 图3 图3 所示。疲劳寿命预测采用Basquin公式： <math display="block" id="M1"><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">a</mi></mrow></msub><mo>=</mo><msubsup><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">f</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msubsup><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">f</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mi>b</mi></mrow></msup><mo>=</mo><mn mathvariant="normal">663.07</mn><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">f</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mo>-</mo><mn mathvariant="normal">0.035</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705424&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705436&type= 46.82066727 4.82600021 (a) 应力-应变曲线；(b) 疲劳寿命曲线 式中： σ a 为应力幅； <math id="M2"><msubsup><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">f</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705438&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705448&type= 2.79399991 4.06400013 为疲劳强度系数； N f 为循环次数； b 为疲劳强度指数。 由 图3 图3 和 式(1) 式(1) 可知：35CrMo的抗拉强度 σ b 为745 MPa，屈服强度 σ s 为561 MPa，断后伸长率为20.81%，疲劳强度 σ w 为388 MPa。 2 服役螺栓受载特性计算 1 服役螺栓的螺纹牙型为普通螺纹，其参数如 图4 图4 (a)所示。其中， P 为螺距； H 为原始三角形高度； D 为内螺纹大径； D 1 内螺纹小径； D 2 为内螺纹中径； d 为外螺纹大径； d 1 为外螺纹小径； d 2 为外螺纹中经，公称直径为110 mm，螺距为4 mm，牙底圆弧半径为0.5 mm。根据预紧螺栓的受载情况，螺栓底部设置固定约束，顶部施加拉伸载荷300 MPa，如 图4 图4 (b)所示。 (a) 螺纹参数；(b) 载荷条件 依据螺纹几何尺寸建立服役螺栓的三维几何模型和有限元模型，并利用Abaqus 6.14-5软件对其进行弹塑性有限元分析，获得了拉伸预紧载荷作用下服役螺栓的受载特性，如 图5 图5 所示。通过对计算结果的分析可知，螺栓的应力、应变分布主要位于螺纹牙底圆角处，且最大应力、最大应变均位于第一圈牙底圆角，这与工程实际中螺栓的常规断裂失效位置一致。可见，螺纹牙底圆角直接影响着服役螺栓的受载特性。 (a) 应力 σ 分布；(b) 应变 ε 分布 3 螺栓疲劳寿命预测与结果讨论 1 3.1 　不同参数螺纹的应力集中分析 2 螺栓是典型的缺口构件，应力集中系数为其重要的结构特征之一。螺纹参数的改变直接影响螺栓的应力集中系数分布，这决定着螺栓工程服役的疲劳性能。对于普通螺纹而言，几何参数主要包括公称直径 M ，螺纹圈数 N ，螺纹螺距 P 以及螺纹底角 R 。为建立等效螺栓模拟构件，本研究首先研究了不同螺纹参数对应力集中系数 K 的影响，具体如 图6 图6 所示。 (a) 公称直径；(b) 螺纹圈数；(c) 螺纹螺距；(d) 螺纹底角 从 图6 图6 可以发现，随着公称直径 M 、圈数 N 的增加，螺纹应力集中系数 K 呈现下降趋势，且变化逐渐趋于平缓；螺距 P 、底角半径 R 对螺纹应力集中系数的影响呈线性关系，螺距 P 增加使得应力集中系数增大，而螺纹底角 R 增加使得应力集中系数减小。 3.2 　螺栓模拟构件的建立 2 为建立服役螺栓的模拟构件，本研究首先考虑螺栓预紧受载条件，利用有限元仿真计算得到了服役螺栓的应力集中系数： <math display="block" id="M3"><msub><mrow><mi>K</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>=</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>z</mi></mrow></msub><mo>/</mo><mi>F</mi><mo>≈</mo><mn mathvariant="normal">4.88</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705460&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705441&type= 23.45266724 3.72533321 其中： <math id="M4"><msub><mrow><mi>σ</mi></mrow><mrow><mi>z</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705452&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705461&type= 2.87866688 3.72533321 为线弹性条件下轴向拉伸应力； <math id="M5"><mi>F</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705462&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705442&type= 2.62466669 2.62466669 为施加的名义应力。 以应力集中系数作为等效参量，结合服役螺栓的应力集中系数以及螺纹参数对应力集中影响分析结果，初步选定服役螺栓模拟构件的参数为 N =1， M =10 mm， P =1 mm， R =0.2 mm，即如 图 7 图 7 (a)所示的单圈螺纹模型。为考虑加工制作及试验操作等要求，进一步对服役螺栓模拟构件作简化，将其等效为与螺纹断面一致的单圈圆环缺口模型，如 图7 图7 (b)所示。 (a) 单圈螺纹模型；(b) 单圈圆环状缺口模型 为验证服役螺栓模拟构件模型的有效性，对单圈螺纹和圆环缺口模型的应力集中系数进行了计算，其计算结果分别为5.16和5.18。 经比较分析可知，螺栓模拟构件与工程实际服役螺栓的应力集中系数近似相等，相对误差在5%左右，而模拟构件的螺纹缺口模型的应力集中系数与圆环缺口模型的基本一致，这验证了基于圆环缺口模型的服役螺栓模拟构件的有效性。 3.3 　考虑综合修正系数的疲劳寿命评估 2 通过对比现有的疲劳寿命预测理论和方法，以及对疲劳缺口系数、表面粗糙度系数和尺寸系数的定量分析，PAN等 [ 25 25 - 27 27 25-27 ] 提出了针对模拟构件疲劳寿命的综合修正系数，其定义如下： <math display="block" id="M6"><msub><mrow><mi>K</mi></mrow><mrow><mi>σ</mi></mrow></msub><mo>=</mo><mfrac><mrow><msub><mrow><mi>C</mi></mrow><mrow><mi mathvariant="normal">S</mi></mrow></msub><mo>⋅</mo><msub><mrow><mi>C</mi></mrow><mrow><mi mathvariant="normal">D</mi></mrow></msub></mrow><mrow><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">f</mi></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705455&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705454&type= 18.96533203 8.21266651 式中： K σ 为疲劳寿命的综合修正系数； C S 为表面粗糙度系数； C D 为尺寸系数； K f 为疲劳缺口系数。 若疲劳寿命曲线Basquin公式中光滑试样的疲劳强度系数为 <math id="M7"><msubsup><mrow><mi>σ</mi></mrow><mrow><mn mathvariant="normal">0</mn><mi mathvariant="normal">f</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705475&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705456&type= 3.89466691 4.06400013 ，则模拟构件的疲劳强度系数 <math id="M8"><msubsup><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">f</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705466&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705464&type= 3.64066648 4.06400013 可修正为： <math display="block" id="M9"><msubsup><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">f</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msubsup><mo>=</mo><msub><mrow><mi>K</mi></mrow><mrow><mi>σ</mi></mrow></msub><mo>⋅</mo><msubsup><mrow><mi>σ</mi></mrow><mrow><mn mathvariant="normal">0</mn><mi mathvariant="normal">f</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705488&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705476&type= 16.67933464 4.06400013 根据 式(3) 式(3) 可知，对于一定强度的钢，表面粗糙度越大，粗糙度系数 C S 越小，疲劳寿命将越低；疲劳缺口系数 K f 越大，缺口对疲劳强度影响越大，寿命相对越低，从而使 K σ 变小，导致螺栓构件的疲劳寿命越低。对于尺寸系数 C D ，随着构件尺寸越大，疲劳强度越低而变小，而 K σ 也随之变小，即也导致螺栓构件的疲劳寿命越低。 各系数具体定义如下。 1) 疲劳缺口系数 K f 。疲劳缺口系数 K f 定义 [ 28 28 ] 为： <math display="block" id="M10"><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">f</mi></mrow></msub><mo>=</mo><mfrac><mrow><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">m</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow></msub></mrow><mrow><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">n</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">h</mi></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705469&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705468&type= 16.84866714 8.21266651 式中： <math id="M11"><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">m</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705470&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705478&type= 8.55133343 3.72533321 和 <math id="M12"><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">n</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">h</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705480&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705471&type= 6.94266701 3.72533321 分别为光滑试样和缺口试样的疲劳强度。 K f 越大，缺口试样疲劳强度降低越明显。 利用Peterson缺口敏感系数计算方法，可得： <math display="block" id="M13"><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">f</mi></mrow></msub><mo>=</mo><mn mathvariant="normal">1</mn><mo>+</mo><mfrac><mrow><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">t</mi></mrow></msub><mo>-</mo><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mfrac><mrow><msub><mrow><mi>a</mi></mrow><mrow><mi mathvariant="normal">p</mi></mrow></msub></mrow><mrow><mi>r</mi></mrow></mfrac></mrow></mfrac><mo>≈</mo><mn mathvariant="normal">3.70</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705504&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705481&type= 33.69733429 12.53066730 其中： K t 为理论应力集中系数，这里取单圈圆环缺口模型的应力集中系数5.18； r 为缺口根部半径， r =0.2 mm； a p 为与晶粒尺寸和载荷有关的材料常数。 对于抗拉强度较高的钢材，在轴向拉压或弯曲载荷作用时， a p 可用下式进行估算 [ 27 27 ] ： <math display="block" id="M14"><msub><mrow><mi>a</mi></mrow><mrow><mi mathvariant="normal">p</mi></mrow></msub><mo>=</mo><mn mathvariant="normal">0.025</mn><mtext> </mtext><mn mathvariant="normal">4</mn><mo>×</mo><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">2</mn><mtext> </mtext><mn mathvariant="normal">079</mn></mrow><mrow><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">u</mi></mrow></msub></mrow></mfrac></mrow></mfenced></mrow><mrow><mn mathvariant="normal">1.8</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705505&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705492&type= 36.74533463 10.24466610 式中： σ u 为材料的抗拉强度，MPa。 2) 粗糙度系数 C S 。螺栓模拟构件的粗糙度系数 C S 依据钢表面粗糙度系数经验图表来确定 [ 29 29 ] ，对于粗糙度为1.6 μm的预紧螺栓，查表可得 C S =0.895。 3) 尺寸系数 C D 。利用有限元仿真确定螺栓模拟构件的应力场，应力提取路径为箭头所指的等效缺口平分线，并提取应力路径上每个节点的等效Mises应力，如 图8 图8 (a)所示。通过参数拟合(如 图8 图8 (b)所示)最终确定应力公式 [ 30 30 ] 为 <math display="block" id="M15"><mi>σ</mi><mo>=</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">x</mi></mrow></msub><mfenced open="[" close="" separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mfenced separators="|"><mrow><mfrac><mrow><mi>x</mi></mrow><mrow><mi>r</mi></mrow></mfrac></mrow></mfenced><mo>+</mo><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><mi>x</mi></mrow><mrow><mi>r</mi></mrow></mfrac></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><mfenced open="" close="]" separators="|"><mrow><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><mi>x</mi></mrow><mrow><mi>r</mi></mrow></mfrac></mrow></mfenced></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup><mo>+</mo><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msub><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><mi>x</mi></mrow><mrow><mi>r</mi></mrow></mfrac></mrow></mfenced></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup></mrow></mfenced></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705495&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705483&type= 70.95066071 11.43000031 其中：拟合系数 a 1 =-0.593 1， a 2 =0.202 4， a 3 =-0.028 7， a 4 =0.001 5。 (a) 应力提取路径；(b) 应力公式参数拟合 对于几何相似试样，令 r 0 与 r n 分别表示基本尺寸与缩放后缺口试样的缺口根部半径，则二者比值就是试样整体尺寸的比值，故可得尺寸系数为： <math display="block" id="M16"><msub><mrow><mi>C</mi></mrow><mrow><mi mathvariant="normal">D</mi></mrow></msub><mo>=</mo><mfrac><mrow><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">n</mi></mrow></msub></mrow><mrow><msub><mrow><mi>σ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><msub><mrow><mi>r</mi></mrow><mrow><mi mathvariant="normal">n</mi></mrow></msub></mrow></mfrac><mfrac><mrow><mi>f</mi><mfenced separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">2</mn><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><msub><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac></mrow></mfenced></mrow><mrow><mi>f</mi><mfenced separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">2</mn><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><msub><mrow><mi>r</mi></mrow><mrow><mi mathvariant="normal">n</mi></mrow></msub></mrow></mfrac></mrow></mfenced></mrow></mfrac><mo>=</mo><mi>s</mi><mfrac><mrow><mi>f</mi><mfenced separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">2</mn><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><msub><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac></mrow></mfenced></mrow><mrow><mi>f</mi><mfenced separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">2</mn><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><msub><mrow><mi>r</mi></mrow><mrow><mi mathvariant="normal">n</mi></mrow></msub></mrow></mfrac></mrow></mfenced></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705487&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705486&type= 51.30799866 19.98133469 式中， σ n 为大尺寸缺口试样的疲劳强度， σ 0 为基本尺寸缺口试样的疲劳强度， s = r 0 / r n 为缩放系数， a 0 为临界裂纹尺寸。 一般情况下，试样尺寸越大，疲劳强度 σ n 越小，因此， C D 随构件尺寸增大而变小。通过拟合得到尺寸系数与缺口半径和缩放系数的关系，如 图9 图9 所示。从 图9 图9 可见：随着缺口半径的增大、缩放系数 s 的减小，尺寸系数均减小。 (a) 尺寸系数 C D 与缺口半径的关系；(b) 尺寸系数 C D 与缩放系数的关系 根据 式(3) 式(3) 可知，尺寸系数 C D 减小会使疲劳寿命综合修正系数 K σ 降低，最终导致试样的疲劳强度降低。尺寸系数与缩放系数之间具有较好的线性关系，拟合公式为： <math display="block" id="M17"><msub><mrow><mi>C</mi></mrow><mrow><mi mathvariant="normal">D</mi></mrow></msub><mo>=</mo><mn mathvariant="normal">0.75</mn><mo>+</mo><mn mathvariant="normal">0.25</mn><mi>s</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705515&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705514&type= 25.23066711 3.72533321 利用综合修正系数法，根据螺栓35CrMo钢光滑试样的疲劳寿命曲线( 图3 图3 )，对服役螺栓的模拟构件进行疲劳寿命预测，结果如 图10 图10 所示。图中，红线贯穿的数据点代表光滑试样的载荷和疲劳寿命的关系。由于疲劳缺口系数 K f 的影响，应力与寿命关系降到黑线附近，进一步，由于粗糙度 C S 和尺寸效应 C D 的影响，应力与寿命关系降到绿线附近。 由于疲劳试验中光滑试样和螺栓实际受载的应力比不同，因此，采用Goodman公式对其 S - N 曲线加以修正，得到服役螺栓模拟构件的疲劳寿命公式为： <math display="block" id="M18"><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">a</mi></mrow></msub><mo>=</mo><mn mathvariant="normal">131</mn><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">f</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mo>-</mo><mn mathvariant="normal">0.028</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705525&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705517&type= 25.90800095 4.82600021 同时，为验证本研究所提出方法的有效性，基于商业疲劳分析软件FE-safe ® 2016计算了服役螺栓模拟构件的疲劳寿命，并将计算结果与综合影响系数法预测的疲劳寿命进行对比分析，如 图11 图11 所示。 (a) 强度安全系数云图；(b) 与综合系数法对比 从 图11 图11 可以看出，模拟结果获得的最小疲劳寿命位于螺纹牙底位置，这与工程实际疲劳破坏位置一致；综合影响系数法预测的疲劳寿命曲线与有限元模拟的结果也基本一致。因此，验证了本研究提出的基于综合影响系数法预测螺栓模拟构件疲劳寿命准确性较高、使用性较强，可为服役构件的安全可靠性评估提供一定的指导。 4 结论 1 1) 根据35CrMo钢的力学性能，分析了服役螺栓的受载特性，并讨论了不同螺纹参数对应力集中系数的影响。最大应力、应变均位于第一圈牙底圆角，这与工程中螺栓的常规断裂失效位置一致。随着公称直径、螺纹圈数和底角半径的增加，螺纹应力集中系数呈现下降趋势，而螺距的增加使得应力集中系数逐渐增大。 2) 以应力集中系数作为等效参量，结合实际服役螺栓应力集中系数及螺纹参数对应力集中影响分析，建立了预紧螺栓的模拟构件模型。分析可知，模拟构件的螺纹缺口模型的应力集中系数与圆环缺口模型的基本一致，计算分析时可用圆环缺口模型代替螺纹缺口模型。 3) 基于服役螺栓的模拟构件模型，利用课题组提出的综合修正系数法实现了对模拟构件的疲劳寿命预测。综合影响系数法预测的疲劳寿命曲线与有限元模拟预测的结果基本一致，验证了本研究所提出方法具有较高准确性及较好适用性。 参考文献 GAO C , YANG M Q , PANG J C , et al . Abnormal relation between tensile and fatigue strengths for a high-strength low-alloy steel [J]. Materials Science and Engineering: A , 2022 , 832 : 142418 . Abnormal relation between tensile and fatigue strengths for a high-strength low-alloy steel 孟亚惠 . 高强度螺栓断裂失效分析 [J]. 热加工工艺 , 2017 , 46 ( 6 ): 254 - 255 . 高强度螺栓断裂失效分析 MENG Yahui . Analysis of failure reason of high strength bolt [J]. Hot Working Technology , 2017 , 46 ( 6 ): 254 - 255 . Analysis of failure reason of high strength bolt 安会芬 . 35CrMo钢缸盖螺栓断裂失效分析 [J]. 热加工工艺 , 2012 , 41 ( 8 ): 223 - 224 . 35CrMo钢缸盖螺栓断裂失效分析 AN Huifen . Failure analysis of 35CrMo cylinder cover bolt [J]. Hot Working Technology , 2012 , 41 ( 8 ): 223 - 224 . Failure analysis of 35CrMo cylinder cover bolt 赵海丽 , 陈琦峰 , 王继良 . 35CrMo钢螺栓断裂失效分析 [J]. 金属热处理 , 2007 , 32 ( 6 ): 114 - 115 . 35CrMo钢螺栓断裂失效分析 ZHAO Haili , CHEN Qifeng , WANG Jiliang . Failure analysis of the fracture bolts made of 35CrMo steel [J]. Heat Treatment of Metals , 2007 , 32 ( 6 ): 114 - 115 . Failure analysis of the fracture bolts made of 35CrMo steel 杨建东 , 赵琨 , 李玲 , 等 . 浅析俄罗斯萨扬-舒申斯克水电站7号和9号机组事故原因 [J]. 水力发电学报 , 2011 , 30 ( 4 ): 226 - 234 . 浅析俄罗斯萨扬-舒申斯克水电站7号和9号机组事故原因 YANG Jiandong , ZHAO Kun , LI Ling , et al . Analysis on the causes of units 7 and 9 accidents at Sayano-Shushenskaya hydropower station [J]. Journal of Hydroelectric Engineering , 2011 , 30 ( 4 ): 226 - 234 . Analysis on the causes of units 7 and 9 accidents at Sayano-Shushenskaya hydropower station ZHAO Libin , LIU Fengrui , ZHANG Jianyu . 3D numerical simulation and fatigue life prediction of high strength threaded bolt [J]. Key Engineering Materials , 2009 , 417/418 : 885 - 888 . 3D numerical simulation and fatigue life prediction of high strength threaded bolt 欧阳卿 . 高强螺栓受力及疲劳性能研究 [D]. 长沙 : 湖南大学 , 2013 . 高强螺栓受力及疲劳性能研究 OUYANG Qing . Study on stress and fatigue performance of high-strength bolts [D]. Changsha : Hunan University , 2013 . Study on stress and fatigue performance of high-strength bolts 谢里阳 , 徐灏 . 螺栓受力分析的组合结构法及其疲劳寿命估算 [J]. 机械设计 , 1986 , 3 ( 4 ): 15 - 20 . 螺栓受力分析的组合结构法及其疲劳寿命估算 XIE Liyang , XU Hao . Constitutional structure method for the stress analysls of a bolt and fatigue life prediction [J]. Journal of Machine Design , 1986 , 3 ( 4 ): 15 - 20 . Constitutional structure method for the stress analysls of a bolt and fatigue life prediction 林晓龙 . 高强度螺栓的应力分析及结构疲劳强度优化 [D]. 沈阳 : 东北大学 , 2012 . 高强度螺栓的应力分析及结构疲劳强度优化 LIN Xiaolong . Stress analysis and structural fatigue strength optimization of high strength bolts [D]. Shenyang : Northeastern University , 2012 . Stress analysis and structural fatigue strength optimization of high strength bolts 杜静 , 丁帅铭 , 王秀文 , 等 . MW级风力发电机塔筒环形法兰连接高强度螺栓疲劳评估 [J]. 机械设计 , 2014 , 31 ( 1 ): 75 - 79 . MW级风力发电机塔筒环形法兰连接高强度螺栓疲劳评估 DU Jing , DING Shuaiming , WANG Xiuwen , et al . Fatigue assessment of bolted ring flange connections in MW-class wind turbine tower [J]. Journal of Machine Design , 2014 , 31 ( 1 ): 75 - 79 . Fatigue assessment of bolted ring flange connections in MW-class wind turbine tower 练继建 , 王婷杰 , 王海军 , 等 . 基于原型振动测试的水轮机顶盖螺栓疲劳分析 [J]. 排灌机械工程学报 , 2020 , 38 ( 3 ): 248 - 253 . 基于原型振动测试的水轮机顶盖螺栓疲劳分析 LIAN Jijian , WANG Tingjie , WANG Haijun , et al . Fatigue analysis of hydro-turbine head cover bolts based on prototype vibration measurement [J]. Journal of Drainage and Irrigation Machinery Engineering , 2020 , 38 ( 3 ): 248 - 253 . Fatigue analysis of hydro-turbine head cover bolts based on prototype vibration measurement TARAKANOV P . Cracked high-strength bolt under cycling and hydrogen environment durability estimation using accurate and approximate models [J]. Procedia Engineering , 2015 , 109 : 403 - 409 . Cracked high-strength bolt under cycling and hydrogen environment durability estimation using accurate and approximate models MARCELO A L , UEHARA A Y , UTIYAMA R M , et al . Fatigue properties of high strength bolts [J]. Procedia Engineering , 2011 , 10 : 1297 - 1302 . Fatigue properties of high strength bolts 刘胜祥 , 宋晓萍 , 曾毅 , 等 . 随机载荷作用下的风力机高强度螺栓疲劳特性分析 [J]. 太阳能学报 , 2013 , 34 ( 2 ): 201 - 206 . 随机载荷作用下的风力机高强度螺栓疲劳特性分析 LIU Shengxiang , SONG Xiaoping , ZENG Yi , et al . Fatigue properties analysis of high strength bolts of wind turbine under random loads [J]. Acta Energiae Solaris Sinica , 2013 , 34 ( 2 ): 201 - 206 . Fatigue properties analysis of high strength bolts of wind turbine under random loads 鄂加强 , 李振强 , 陈燕 , 等 . 兆瓦级风力机组偏航齿圈螺栓联接强度分析 [J]. 中南大学学报(自然科学版) , 2015 , 46 ( 1 ): 310 - 316 . 兆瓦级风力机组偏航齿圈螺栓联接强度分析 E Jiaqiang , LI Zhenqiang , CHEN Yan , et al . Strength analysis of bolt connection on yaw gear of MW wind turbine generator system [J]. Journal of Central South University(Science and Technology) , 2015 , 46 ( 1 ): 310 - 316 . Strength analysis of bolt connection on yaw gear of MW wind turbine generator system 李秀珍 , 邓华 , 吕杏梅 , 等 . 兆瓦级风力机偏航齿圈高强度联接螺栓螺纹表面裂纹扩展分析 [J]. 中南大学学报(自然科学版) , 2014 , 45 ( 1 ): 91 - 98 . 兆瓦级风力机偏航齿圈高强度联接螺栓螺纹表面裂纹扩展分析 LI Xiuzhen , DENG Hua , LÜ Xingmei , et al . Growth analysis of surface crack on high-strength connecting bolt thread of yaw gear ring of MW wind turbine generator system [J]. Journal of Central South University(Science and Technology) , 2014 , 45 ( 1 ): 91 - 98 . Growth analysis of surface crack on high-strength connecting bolt thread of yaw gear ring of MW wind turbine generator system 王源晨 , 宋竹满 , 李瑞 , 等 . 8.8级螺栓用ML40Cr钢的疲劳性能和裂纹扩展行为 [J]. 材料研究学报 , 2019 , 33 ( 10 ): 771 - 775 . 8.8级螺栓用ML40Cr钢的疲劳性能和裂纹扩展行为 WANG Yuanchen , SONG Zhuman , LI Rui , et al . Fatigue properties and crack growth behavior of ML40Cr steel for 8.8-grade bolts [J]. Chinese Journal of Materials Research , 2019 , 33 ( 10 ): 771 - 775 . Fatigue properties and crack growth behavior of ML40Cr steel for 8.8-grade bolts 黄宁 . 大型结构件的疲劳寿命预测方法研究 [D]. 长沙 : 中南大学 , 2013 . 大型结构件的疲劳寿命预测方法研究 HUANG Ning . Study on fatigue life prediction method of large structural parts [D]. Changsha : Central South University , 2013 . Study on fatigue life prediction method of large structural parts YILMAZ O , JIMÉNEZ-PEÑA C , DEBRUYNE D . Fatigue response of high-strength steel bolted connections under in-plane shear and out-of-plane bending loading modes [J]. Procedia Structural Integrity , 2019 , 19 : 302 - 311 . Fatigue response of high-strength steel bolted connections under in-plane shear and out-of-plane bending loading modes POOVAKAUD V V , JIMÉNEZ-PEÑA C , TALEMI R , et al . Assessment of fretting fatigue in high strength steel bolted connections with simplified Fe modelling techniques [J]. Tribology International , 2020 , 143 : 106083 . Assessment of fretting fatigue in high strength steel bolted connections with simplified Fe modelling techniques YANG X , LEI H , CHEN Y . Constant amplitude fatigue test research on M20 high-strength bolts in grid structure with bolt-sphere joints [J]. Advances in Structural Engineering , 2017 , 20 ( 10 ): 1466 - 1475 . Constant amplitude fatigue test research on M20 high-strength bolts in grid structure with bolt-sphere joints YANG Xu , LEI Honggang . Constant amplitude fatigue test of high strength bolts in grid structures with bolt- sphere joints [J]. Steel & Composite Structures , 2017 , 25 ( 5 ): 571 - 579 . Constant amplitude fatigue test of high strength bolts in grid structures with bolt- sphere joints PANG J C , LI S X , WANG Z G , et al . General relation between tensile strength and fatigue strength of metallic materials [J]. Materials Science and Engineering: A , 2013 , 564 : 331 - 341 . General relation between tensile strength and fatigue strength of metallic materials PANG J C , LI S X , WANG Z G , et al . Relations between fatigue strength and other mechanical properties of metallic materials [J]. Fatigue & Fracture of Engineering Materials & Structures , 2014 , 37 : 958 - 976 . Relations between fatigue strength and other mechanical properties of metallic materials PAN L , PANG J C , ZHANG M X , et al . A prediction method for fatigue life of large moving components as function of scale factor: a case of motor-generator rotor [J]. Fatigue & Fracture of Engineering Materials & Structures , 2018 , 41 : 1089 - 1101 . A prediction method for fatigue life of large moving components as function of scale factor: a case of motor-generator rotor PAN L , PANG J C , XIE Y M , et al . A fatigue life prediction method from magnetic pole material to simulation part [J]. International Journal of Damage Mechanics , 2019 , 28 ( 9 ): 1438 - 1456 . A fatigue life prediction method from magnetic pole material to simulation part 潘龙 . 转子磁极关键构件疲劳寿命预测方法研究 [D]. 抚顺 : 辽宁石油化工大学 , 2018 . 转子磁极关键构件疲劳寿命预测方法研究 PAN Long . Study on fatigue life prediction method for key component of magnetic pole of a generator rotor [D]. Fushun : Liaoning Shihua University , 2018 . Study on fatigue life prediction method for key component of magnetic pole of a generator rotor YU M , DUQUESNAY D , TOPPER T . Notch fatigue behaviour of SAE1045 steel [J]. International Journal of Fatigue , 1988 , 10 ( 2 ): 109 - 116 . Notch fatigue behaviour of SAE1045 steel LEE Y L . Fatigue testing and analysis: theory and practice [M]. Amsterdam : Elsevier Butterworth-Heinemann , 2005 . Fatigue testing and analysis: theory and practice 孙迪 , 胡振东 . 基于TCD理论的缺口试样疲劳强度尺寸系数计算 [J]. 力学季刊 , 2015 , 36 ( 2 ): 288 - 295 . 基于TCD理论的缺口试样疲劳强度尺寸系数计算 SUN Di , HU Zhendong . Research of the size factor of fatigue strength base on TCD theory [J]. Chinese Quarterly of Mechanics , 2015 , 36 ( 2 ): 288 - 295 . Research of the size factor of fatigue strength base on TCD theory 10.11817/j.issn.1672-7207.2023.05.010 TH114 A 1672-7207（2023）05-1748-10 杨梦起, 杨文军, 庞建超, 等. 35CrMo钢螺栓疲劳寿命预测方法研究[J]. 中南大学学报(自然科学版), 2023, 54(5): 1748-1757. 杨梦起, 杨文军, 庞建超, 等. 35CrMo钢螺栓疲劳寿命预测方法研究[J]. 中南大学学报(自然科学版), 2023, 54(5): 1748-1757. YANG Mengqi, YANG Wenjun, PANG Jianchao, et al. Fatigue life prediction method of 35CrMo alloy steel bolt[J]. Journal of Central South University(Science and Technology), 2023, 54(5): 1748-1757. YANG Mengqi, YANG Wenjun, PANG Jianchao, et al. Fatigue life prediction method of 35CrMo alloy steel bolt[J]. Journal of Central South University(Science and Technology), 2023, 54(5): 1748-1757. 南方电网调峰调频发电有限公司检修试验分公司科技项目(022200KK52180006)；国家自然科学基金资助项目(51871224) Project(022200KK52180006) supported by the Science and Technology Project of Maintenance and Testing Branch in China Southern Power Grid Power Generation Co. Ltd.; Project(51871224) supported by the National Natural Science Foundation of China 2022-01-19 2022-04-01 2023-05-26 2023-06-20 38910187 中南大学学报(自然科学版) 5 54