疲劳分析流程fatigue

1、摘要：疲劳破坏是结构的主要失效形式，疲劳失效研究在结构安全分析中扮演着举足轻重的角色。因此结构的疲劳强度和疲劳寿命是其强度和可靠性研究的主要内容之一。机车车辆结构的疲劳设计必须服从一定的疲劳机理,并在系统结构的可靠性安全设计中考虑复合的疲劳设计技术的应用。国内的机车车辆主要结构部件的疲劳寿命评估和分析采用复合的疲劳设计技术，国外从疲劳寿命的理论计算和疲劳试验两个方面在疲劳研究和应用领域有很多新发展的理论方法和技术手段。不论国内国外，一批人几十年如一日致力于疲劳的研究，对疲劳问题研究贡献颇多。关键词：疲劳  UIC标准  疲劳载荷  IIW标准&

2、#160;  S-N曲线   机车车辆一、国内外轨道车辆的疲劳研究现状6月30日15时，备受关注的京沪高铁正式开通运营。作为新中国成立以来一次建设里程最长、投资最大、标准最高的高速铁路，京沪高铁贯通“三市四省”，串起京沪“经济走廊”。京沪高铁的开通，不仅乘客可以享受到便捷与实惠，沿线城市也需面对高铁带来的机遇和挑战。在享受这些待遇的同时，专家指出，各省市要想从中分得一杯羹，配套设施建设以及机车车辆的安全性绝对不容忽略。根据机车车辆的现代设计方法，对结构在要求做到尽可能轻量化的同时，也要求具备高度可靠性和足够的安全性。这两者之间常常出现矛盾，因此

3、，如何准确研究其关键结构部件在运行中的使用寿命以及如何进行结构的抗疲劳设计是结构强度寿命预测领域研究中的前沿课题。在随机动载作用下的结构疲劳设计更是成为当前机车车辆结构疲劳设计的研究重点，而如何预测关键结构和部件的疲劳寿命又是未来机车车辆结构疲劳设计的重要发展方向之一。机车车辆承受的外部载荷大部分是随时间而变化的循环随机载荷。在这种随机动载荷的作用下，机车车辆的许多构件都产生动态应力，引起疲劳损伤，而损伤累积后的结构破坏的形式经常是疲劳裂纹的萌生和最终结构的断裂破坏。随着国内铁路运行速度的不断提高，一些关键结构部件，如转向架的构架、牵引拉杆等都出现了一些断裂事故。因此，机车车辆的结构疲劳设计已

4、经逐渐成为机车车辆新产品开发前期的必要过程之一，而通过有效的计算方法预测结构的疲劳寿命是结构设计的重要目标。1.1国外早在十九世纪后期德国工程师Wohler系统论述了疲劳寿命和循环应力的关系并提出了S-N曲线和疲劳极限的概念以来，国内外疲劳领域的研究已经产生了大量新的研究方法和研究成果。结构疲劳设计中主要有两方面的问题：一是用一定材料制成的构件的疲劳寿命曲线；二是结构件的工作应力谱，也就是载荷谱。载荷谱包括外部的载荷及动态特性对结构的影响。根据疲劳寿命曲线和工作应力谱的关系，有3种设计概念：静态设计(仅考虑静强度)；工作应力须低于疲劳寿命曲线的疲劳耐久限设计；根据工作强度设计，即运用实际使用条

5、件下的载荷谱。实际载荷因为受到车辆等诸多因素的影响而有相当大的离散性，它严重地影响了载荷谱的最大应力幅值、分布函数及全部循环数。为了对疲劳寿命进行准确的评价，必须知道设计谱的存在概率，并且考虑实际载荷离散性，才可以确定结构可靠的疲劳寿命。20世纪60年代，世界上第一条高速铁路建成，自那时起，一些国外高速铁路发达国家已经深入研究机车车辆结构轻量化带来的关键结构部件的疲劳强度和疲劳寿命预测问题。其中，包括日本对车轴和焊接构架疲劳问题的研究；法国和德国采用试验台仿真和实际线路相结合的技术开发出试验用的机车车辆疲劳分析方法；英国和美国对转向架累计损伤疲劳方面的研究等等。在这些研究中提出了大量有效的疲劳

6、寿命的预测研究方法。1.2、国内1.2.1 国内疲劳研究现状与方法国内铁路相关的科研院所对结构的疲劳寿命也展开了大量的研究和分析，并且得到了很多研究寿命预测方法的有效成果。疲劳寿命的估算主要分为裂纹形成阶段寿命估算和裂纹扩展阶段寿命估算两部分。常用的疲劳寿命计算方法包括名义应力寿命法、局部应变寿命法和裂纹扩展计算法。. SN方法名义应力法,又称SN方法,强调了长疲劳寿命和耐久性的限制,或者是假定疲劳失效不会发生时的安全应力。它以材料或零部件的疲劳寿命曲线为基础,利用名义应力或局部应力和寿命之间的关系,主要应用在线性弹性应力的分析中。SN曲线也称作应力循环曲线,是说明结构疲劳过程的重要

7、方法。这个方法经常被用在变化的环境里,包括那些几乎没有塑性的长寿命疲劳问题,以及部件裂纹的萌生或裂纹的生长等问题。同时,该法主要应用在循环载荷是弹性的环境里,也就是意味着结构限制在寿命轴的疲劳循环数量必须大于10000次,这也是为了确保结构没有大的塑性发生。实际上在小于10000次循环的环境,要使用SN方程必须更加小心。通常要使用最终的强度或者是真实的断裂应力进行分析。应力法以循环应力作为造成疲劳的原因,寿命的预测主要是应力循环曲线(SN曲线)。它是由常幅加载试验获得的,尽管大部分情况中实际工程结构承受的多是变幅载荷,但在疲劳寿命的估算,依然有一定的应用价值。如果是变幅载荷,需要做出SN曲线进

8、行修正。变化的应力时间历程一般要进行计数,用累积频次分布的形式进行寿命的预测计算。经常使用如下的公式表示整条SN曲线:(-A)m N= D  或= A(1+C/N)式中:1/m为SN曲线的负斜率;A为材料的疲劳极限;D为材料常数;为SN曲线的负斜率;C为材料常数。. EN方法局部应变法(也称EN法)以局部应变为基准, 局部应变或应变寿命法,通常参考裂纹萌生方法,其关注的仅仅是裂纹的萌生;结构的疲劳寿命, 通常是所有的寿命(SN或者名义应变),以及裂纹的产生或损伤误差分析之和,关注的是直到结构断裂产生的循环数。EN分析法使用循环应力应变建模和Neuber弹塑性修正。

9、可以选择分析参数的范围,包括平均应力修正模型、置信度参数、制造的细节(表面的粗糙度)以及材料的热处理。. 裂纹扩展方法Paris最早提出的裂纹前沿应力强度因子范围K和裂纹扩展速率da/dN之间的经验关系是计算疲劳裂纹扩展寿命的基础8。基于断裂力学的裂纹扩展寿命预测法主要用于较长裂纹的损伤容限设计,最近这种方法有被拓展到结构耐久性分析领域的趋势。利用不同技术的这几种方法有不同的精度。实际上预测总的寿命、裂纹萌生和裂纹扩展3种方法很少在相同的问题上同时使用。这是因为不同的设计理念决定了不同的分析方法。利用这3种基本方法可以分析结构的点焊、缝焊等多种焊接方式,以及旋转结构、振动疲劳、复合多轴等疲劳寿

10、命的预测问题。1.2.2国内近期疲劳研究领域的发展与突破近年来，随着疲劳理论研究的进步和计算机软硬件技术的迅猛发展，在发达国家的铁路行业，产品的强度设计由原来的主要依据静强度指标和无限寿命设计发展到定量寿命设计，大大提高了产品的使用可靠性，并且降低了产品的生产成本。现代的抗疲劳设计中使用了“一体化耐久管理”(IDMIntegrated Durablity Management)的方法，该方法包括虚拟分析和疲劳试验两个方面，下面将从这两个方面分别叙述疲劳研究和应用的新发展。1、理论分析方法的发展虚拟计算所谓的“虚拟”计算就是在产品的设计阶段，使用软件建立产品有关承载构件的有限元模型，使用软件在虚

11、拟的“新产品构件”上施加载荷进行寿命的计算并可反复调整设计方案，这种计算所需时间短，费用低，降低了产品开发周期和成本。由于疲劳理论的发展结合计算机技术的进步虚拟计算作为一种分析手段，也有很多新而且方便实用的方法。主要包括3个方面，它们分别是用SN方法估算全寿命、用局部应力应变估算裂纹起始寿命、用Paris公式估算裂纹扩展寿命。使用PSD功率谱密度信号计算构件的全寿命传统上人们都是使用时域信号进行寿命或损伤的计算，时域信号一般是应力、应变或载荷。在时域中可以很方便地表示周期性的信号，若用来完全地描述一个随机过程，需要非常长的信号记录，特别是对于有上个节点的有限元模型施加几十个通道的载荷的情况下，

12、使用时域信号计算模型的寿命将非常困难，在这种情况下，疲劳寿命的计算可以在频域里完成。其主要过程是：首先对有限元模型进行频响分析，求得载荷与结构上应力的传递函数；然后，将传递函数乘以载荷的功率谱密度函数得到应力功率谱密度函数最后使用有关的方法由应力功率谱推算结构的寿命。9用应力PSD信号计算构件的寿命的方法有2种：其一是直接使用PSD信号；其二是先根据PSD信号计算出应力范围的循环数，再用它估算寿命。目前，第2种做法使用较为流行。预测多轴应力状态下构件的起始寿命多数情况下，工程构件所受的载荷是多载荷的复杂组合。这就使得构件上的某些点所受的3个主应力处于非比例状态主应力的方向往往也随时变化，在这种

13、应力条件下的疲劳就称为多轴疲劳。多轴疲劳的研究最近也有了新的进展，发展了一种新方法。该方法不仅能够估算多轴应力情况下的寿命，还能够以损伤极坐标图的形式给出各个临界面上的损伤量。用这一方法来估算始裂寿命分以下几步：第一步，在承载构件上的关键点贴上应变花测得3个应变输出，或者在有限元模型上贴虚拟应变花，求得3个虚拟的应变输出；第二步，根据3个应变输出用MrozGarud循环塑性模型计算该关键点的其余的应力应变分量；第三步，根据材料的疲劳寿命曲线，用有关的多轴疲劳损伤寿命模型计算寿命。10预测形状不规则裂纹的扩展寿命最近，有一种计算形状不规则裂纹的扩展寿命的新方法得到了发展，该方法是线弹性断裂力学理

14、论与有限元方法的结合。其主要内容是首先对含裂纹的构件进行三维有限元网格的划分，网格包括裂纹块和非裂纹块，裂纹块由于包含裂纹需要较细致的划分，而非裂纹块可以划分得较粗一些。裂纹块和非裂纹块的组合可以使用多点位移约束法，使得二者的位移在其接触面上保持一致，从而使组合后的网格保持变形协调。疲劳裂纹扩展后新网格的生成可以只对裂纹块进行，非裂纹块阿格保持不变。网格划分完毕后，应用四分之一点位移法或者三维J积分法计算裂纹扩展前沿上的各个节点的应力强度因子，J积分的计算方法有回路法和虚拟节点位移法。裂纹前沿的曲线形状可以由一些离散点(节点)来定义，那么，就可以应用Paris公式分别求这些节点上的裂纹扩展量(

15、也就是这些节点的位移量)a，据此可以确定新的裂纹前沿。该方法能够模拟裂纹扩展过程，跟踪裂纹形状的发展变化，大大提高了裂纹扩展寿命的预测精度。2、疲劳试验使用软件来仿真计算虚拟“新产品”作为抗疲劳设计的一个手段，并不能代替疲劳试验，虚拟计算的目的是为了尽可能地减少疲劳试验的次数。新产品的疲劳可靠性最终要通过疲劳试验来检验，目前国内对铁道车辆零部件的疲劳试验，只能够做程序加载的疲劳试验。如TB1959和TB1960所规定的试验标准，用程序载荷谱加载的疲劳试验一般只能够对产品进行合格检验或产品的性能对比，不能够真正用于研究产品的寿命。在国外，主要采用随机载荷谱加载进行疲劳试验来研究产品的寿命这样就彻

16、底消除了程序载荷谱时高低载荷的加载次序和各加载位置载荷相位差对构件寿命的影响。这种加载方法能够较为真实地模拟构件所承受的载荷，使所得的试验寿命能够真实地反映运用寿命在各个构件之间布置载荷传感器十分困难，而在构件上贴上应变片则很方便，所以一般不能直接采集到构件的随机载荷谱，只能够采集到构件上某些点的应变谱。从运用当中的构件上采集的应变谱若直接用于疲劳试验，虽然试验寿命最为接近实际的运用寿命，但试验的时间将与构件的实际运用时间相当，疲劳试验就失去了其意义。所以，要对采集的应变谱进行一些编辑处理。合理的疲劳编辑应该满足2个方面，其一是经过编辑后的信号与原始信号对构件的总损伤基本一致；其二是编辑前后的

17、信号所产主的损伤分布基本一致。过去，用线性放大载荷信号或者用“等效”的常幅载荷信号来加速疲劳试验的做法，不符合上面所提的第2个方面疲劳试验有单通道和多通道之分，疲劳编辑也有所区别单通道编辑只编辑响应信号(应变谱)，多通道编辑需要同步编辑响应信号(应变谱)和驱动信号(随机载荷谱)。二、CW-200转向架疲劳寿命预测实例图1是转向架焊接构架待评估位置:  图1 CW-200转向架焊接构架待评估位置2.1 有限元模型图2是CW-200转向架有限元模型 图2 CW-200转向架有限元模型有限元模型采用平面四节点单元，单元大小20mm。疲劳试验载荷包括垂向载荷与横

18、向载荷，垂向载荷作用在空气弹簧座上，横向载荷作用在横向减振器座上。2.2 疲劳试验载荷根据UIC标准：垂直力（每个侧架）(取=0.2，=0.1)静态部分：FZS1=FZS2=FZ准静态部分（模拟曲线中的滚动）：FZQ1=FZQ2=+-FZ动态部分（跳动模拟）：FZD1=FZD2=+-FZ横向力（每个转向架）准静态部分：FYQ=+-0.25*(FZ+0.5*m+*g)动态部分：FYD=+-0.25*(FZ+0.5*m+*g)其中，FZ=135KN   m+*g=60KN 图13UIC标准中垂向动载荷模式 图14UIC标准中横向动载荷模式

19、由于产生疲劳破坏的主要为动态部分力，所以在疲劳评估时忽略静态部分和准静态部分对疲劳寿命的影响，所以需加的疲劳载荷为：垂直力（每个侧架）:FZ=2*(0.2*135)=54KN横向力（每转向架）: FY=2*0.25*(135+0.5*60)=82.5KN在UIC标准中,加载分三个阶段，分别以1倍、1.2倍和1.4倍的比例关系逐段强化，每个阶段加载的循环次数分别为6e106、2e106和 e106次。 图15UIC标准中三段加载模式垂向力动态部分与横向力动态部分作用频率与相位均相同，因而视其为一个载荷事件。2.3 有限元计算按照UIC515-4的实验室疲劳实验加载要求，首先将

20、其静载到要有限元模中,然后基于动载与静载的比例关系进行应力换算，即得动载应力（任一点的值都可由下图读出）。图3是部分有限元计算结果:  图3疲劳载荷作用下应力云图 2.4 S-N曲线选取将待评估点所处的结构细部与IIW标准钢结构细部疲劳强度值表中提供的结构细部对比，选择合适的类型FAT级别。CW-200转向架构架各点的疲劳强度级别FAT值列入下表： 表1 待评估位置FAT等级列表编号123456789101112131415接头类型管接头管接头T型接头管接头T型接头端部接头管接头端部接头对接接头对接接头T型接头非承载T接头T型接头管接头T型接头

21、FAT5050805080455045717180100805080 根据各点的疲劳强度级别FAT值，即可在IIW标准中（表2）查找相应的S-N曲线参数,即m、C值。表2  IIW标准中各级别S-N曲线参数（MPa）Palmgren-Miner累积特征抗疲劳-S-N曲线的修正值最初斜率m1=3.0，循环数5e6时的常幅疲劳极限l,x二次斜率m2=5.0,在106循环时截止级别循环数N<5e6，m1=3时S-N曲线的常数C常幅疲劳极限循环数N>5E6，m2=5时S-N曲线的常数C截止极限2552001801601401251121009080716356

22、5045403632282522201816142.27e131.600e131.166e138.192e125.488e123.906e122.810e122.000e121.458e121.012e127.158e115.001e113.512e112.500e111.823e111.280e119.331e106.554e104.390e103.125e102.130e101.600e101.166e108.192e95.488e916614713311810392.182.573.766.358.952.346.441.336.833.229.526.523.620.618.416.2

23、14.713.311.810.36.261e173.474e172.052e171.139e175.840e163.313e161.913e161.086e166.411e153.558e151.959e151.078e155.980e143.393e142.004e141.112e146.565e133.643e131.869e131.060e135.596e123.474e122.052e121.139e125.840e1191.180.972.864.856.750.645.340.536.432.428.725.522.720.218.216.214.613.011.310.18.98

24、.17.36.55.7 然后就可以根据迈纳尔损伤累积理论计算疲劳寿命，这里为疲劳损伤比之和。疲劳寿命极限 ；（对应i）累计损伤率D=计算结果列表如下： 表3  CW-200焊接构架待评估点疲劳损伤比之和标号FAT 垂向横向1倍(加载6e6次)1.2倍(加载2e6次)1.4倍(加载2e6次)150应力范围2.533.5S-N曲线m无C无m无C无m无C无疲劳寿命无限无限无限损伤比000累积损伤比0250应力范围11.213.4415.68S-N曲线m无C无m无C无m=无C无疲劳寿命无限无限无限损伤比000累积损伤比0380应力范围9.911.881

25、3.86S-N曲线m无C无m无C无m无C无疲劳寿命无限无限无限损伤比000累积损伤比0450应力范围16.119.3222.54S-N曲线m无C无m无C无m无C无疲劳寿命无限无限无限损伤比000累积损伤比0580应力范围41.649.9258.24S-N曲线m=5C=3.558e15m=5C=3.558e15m=5C=3.558e15疲劳寿命2.86e71.15e75.31e6损伤比0.210.1740.377累积损伤比0.761645应力范围18.822.5626.32S-N曲线m=5C=2.004e14m=5C=2.004e14m=5C=2.004e14疲劳寿命8.53e73.43e71.

26、59e7损伤比0.070.0580.126累积损伤比0.254750应力范围3.13.724.34S-N曲线m无C无m无C无m无C无疲劳寿命1.58e3无限无限损伤比000累积损伤比0845应力范围6.67.929.24S-N曲线m无C无m无C无m无C无疲劳寿命无限无限无限损伤比000累积损伤比0971应力范围7.38.7610.22S-N曲线m无C无m无C无m无C无疲劳寿命无限无限无限损伤比000累积损伤比01071应力范围11.013.215.4S-N曲线m无C无m无C无m无C无疲劳寿命无限无限无限损伤比000累积损伤比01180应力范围26.531.837.1S-N曲线m无C无m无C无m

27、=5C=3.558e15疲劳寿命无限无限5.06e7损伤比000.04累积损伤比0.0412100应力范围37.244.6452.08S-N曲线m无C无m=5C=1.086e16m=5C=1.086e16疲劳寿命无限6.13e72.83e7损伤比00.0330.071累积损伤比0.1041380应力范围7.278.72410.178S-N曲线m无C无m无C无m无C无疲劳寿命无限无限无限损伤比000累积损伤比01450应力范围8.279.92411.578S-N曲线m无C无m无C无m无C无疲劳寿命无限无限无限损伤比000累积损伤比01580应力范围4.205.045.88S-N曲线m无C无m无C

28、无m无C无疲劳寿命无限无限无限损伤比000累积损伤比0       二、根据这一点的接头形式，查出对应的S-N曲线，方法：1. UIC标准-FAT值-（判断和常幅疲劳极限、截止极限的大小关系）-m、c值-(由公式 )-疲劳寿命极限值N;2. 同时（疲劳载荷引起的应力范围）对应一个循环次数n；3. 最后求得损伤比4. 将疲劳载荷逐步加强到1.2倍和1.4倍，重复以上过程，又会得到两个损伤比；5. 将三个损伤比相加，得到累计损伤比，并判断是否小于1（合格） 一、对于其它每一点重复以上过程，完成对每一评估点的合格与否检验

29、 编号123456789101112131415接头类型管接头管接头T型接头管接头T型接头端部接头管接头端部接头对接接头对接接头T型接头非承载T接头T型接头管接头T型接头FAT5050805080455045717180100805080 Palmgren-Miner累积特征抗疲劳-S-N曲线的修正值最初斜率m1=3.0，循环数5e6时的常幅疲劳极限l,x二次斜率m2=5.0,在106循环时截止级别循环数N<5e6，m1=3时S-N曲线的常数C常幅疲劳极限循环数N>5E6，m2=5时S-N曲线的常数C截止极限255200180160140125112100908

30、07163565045403632282522201816142.27e131.600e131.166e138.192e125.488e123.906e122.810e122.000e121.458e121.012e127.158e115.001e113.512e112.500e111.823e111.280e119.331e106.554e104.390e103.125e102.130e101.600e101.166e108.192e95.488e916614713311810392.182.573.766.358.952.346.441.336.833.229.526.523.620.61

31、8.416.214.713.311.810.36.261e173.474e172.052e171.139e175.840e163.313e161.913e161.086e166.411e153.558e151.959e151.078e155.980e143.393e142.004e141.112e146.565e133.643e131.869e131.060e135.596e123.474e122.052e121.139e125.840e1191.180.972.864.856.750.645.340.536.432.428.725.522.720.218.216.214.613.011.31

32、0.18.98.17.36.55.7  标号FAT 垂向横向1倍(加载6e6次)1.2倍(加载2e6次)1.4倍(加载2e6次)150应力范围2.533.5S-N曲线m无C无m无C无m无C无疲劳寿命无限无限无限损伤比000累积损伤比0250应力范围11.213.4415.68S-N曲线m无C无m无C无m=无C无疲劳寿命无限无限无限损伤比000累积损伤比0380应力范围9.911.8813.86S-N曲线m无C无m无C无m无C无疲劳寿命无限无限无限损伤比000累积损伤比0450应力范围16.119.3222.54S-N曲线m无C无m无C无m无C无疲劳寿命无限无限无限

33、损伤比000累积损伤比0580应力范围41.649.9258.24S-N曲线m=5C=3.558e15m=5C=3.558e15m=5C=3.558e15疲劳寿命2.86e71.15e75.31e6损伤比0.210.1740.377累积损伤比0.761645应力范围18.822.5626.32S-N曲线m=5C=2.004e14m=5C=2.004e14m=5C=2.004e14疲劳寿命8.53e73.43e71.59e7损伤比0.070.0580.126累积损伤比0.254750应力范围3.13.724.34S-N曲线m无C无m无C无m无C无疲劳寿命1.58e3无限无限损伤比000累积损伤比

34、0845应力范围6.67.929.24S-N曲线m无C无m无C无m无C无疲劳寿命无限无限无限损伤比000累积损伤比0971应力范围7.38.7610.22S-N曲线m无C无m无C无m无C无疲劳寿命无限无限无限损伤比000累积损伤比01071应力范围11.013.215.4S-N曲线m无C无m无C无m无C无疲劳寿命无限无限无限损伤比000累积损伤比01180应力范围26.531.837.1S-N曲线m无C无m无C无m=5C=3.558e15疲劳寿命无限无限5.06e7损伤比000.04累积损伤比0.0412100应力范围37.244.6452.08S-N曲线m无C无m=5C=1.086e16m=

35、5C=1.086e16疲劳寿命无限6.13e72.83e7损伤比00.0330.071累积损伤比0.1041380应力范围7.278.72410.178S-N曲线m无C无m无C无m无C无疲劳寿命无限无限无限损伤比000累积损伤比01450应力范围8.279.92411.578S-N曲线m无C无m无C无m无C无疲劳寿命无限无限无限损伤比000累积损伤比01580应力范围4.205.045.88S-N曲线m无C无m无C无m无C无疲劳寿命无限无限无限损伤比000累积损伤比0当前，提速/重载在迎合社会需求的同时，随之也出现了一些值得关注的问题，其中，机车车辆许多构件，尤其是焊接转向架，服役不久频繁发生疲劳断裂，并导致了重大经济损失。客观上，由于提速/重载，机车车辆动负荷加剧，然而，这不过是问题的外因，而内因则是我们对这些机车车辆装备没有注入与提速/重载相匹配的且有相当储备的抗疲劳能力。许多疲劳断裂案例分析表明：在设计阶段，我们仍停留在传统的静强度设计层次上，对焊接转向架尤其缺少抗疲劳设计的理念，缺少相应的技术，缺少必要的基础数据。因此，许多疲劳断裂隐患始于设计阶段。事实上，自从人们认定“预防疲劳断裂”是一个公认的难题以来