基于疲劳驱动力能量损伤参数修正的非线性疲劳寿命预测

疲劳破坏是轨道车辆中焊接结构常见的一种失效形式，这些焊接结构通常的承载状态是承受循环载荷作用，且往往承受变幅载荷作用[1]。对于轨道车辆结构，急需更为合理且实用的疲劳寿命分析模型，以有效解决车辆焊接结构的疲劳问题，避免重大灾害或安全事故，具有十分重要的工程应用研究价值[2-3]。

针对焊接结构疲劳问题，最为常用的疲劳寿命分析模型是线性的Miner法则，计算模型简单且便于工程应用。但该模型的损伤累积规则简化了疲劳失效的过程机理，使疲劳寿命的估计具有不确定性，易出现较大的偏差。大量的工程实际数据也表明，对于复杂的载荷加载，Miner法则的分析结果不准确，与实际结果偏差较大，常不能满足工程要求。对于多级变幅载荷作用时，加载历史效应不能忽略，因载荷加载次序、载荷间的干涉以及实时的疲劳损伤程度等相关因素对损伤累积有显著影响[4-5]。对于实际中普遍存在的变幅载荷作用下疲劳问题，迫切需要建立有效的模型处理变幅载荷下损伤累积过程，从而更好地预测疲劳寿命。

疲劳损伤从本质上分析是一个不可逆、材料性能退化和能量耗散的过程，可通过评估机械性能或物理性能的变化进行判定。非线性损伤累积理论是当前进行疲劳分析的2大类方法之一，所涉及的理论模型有损伤曲线、连续损伤力学、能量法、物理性能退化等模型[6]。虽然这些模型在一定程度上能处理载荷加载历史效应，但是其数学形式多为多层嵌套，计算形式复杂，有的还需引入额外的参数，难以直接应用于工程实际[7-8]。如何合理地描述复杂的载荷加载历史效应，进而构建准确的疲劳寿命预测模型，一直是工程专家迫切需要解决的实际问题。

PCCPL、SP采用无溶性特加强级环氧煤沥青喷涂，干膜厚度600 μm；DIP管外喷涂金属锌 （含锌量不小于99.99%），重量不少于 130 g/m2，再涂敷沥青漆或环氧树脂漆，涂层干馍平均厚度≥70μm，局部最小厚度不小于50μm，最大干膜厚度<250 μm。

近年来，部分学者将应变能密度参数用于跟踪损伤累积行为，提出了一种有效的疲劳寿命预测方法，且得到较好的应用[9-10]。文献[11]结合应变能密度参数，提出了疲劳驱动力能量的概念，建立了适用于高周疲劳状态下疲劳损伤累积的非线性损伤模型，可用于变幅载荷下疲劳寿命的预测，同时还引入应力比函数，探讨了载荷相互作用的影响。但在变幅载荷作用下，疲劳损伤累积规律因载荷幅值的变化而变得非常复杂[12-13]。疲劳损伤累积规律与加载历程密切相关，载荷加载次序、载荷间的干涉、疲劳损伤程度等因素直接影响着疲劳损伤累积规律，相关研究还需进一步深入探讨。

本文基于疲劳驱动力能量损伤参数，利用变幅载荷应力比、载荷循环比和实时疲劳损伤程度构建指数函数，对相邻载荷间的交互因子进行修正，建立改进的非线性疲劳寿命预测模型。所构建的疲劳寿命预测模型有效性和可行性通过不同材料的焊接结构在多级变幅载荷作用下的疲劳试验数据结果给予验证，并与其他数值模型进行对比。

前段时间网上有段曝光问题肥料的视频，我们暂不讨论该肥料是否违规，是否不合格，但是在视频中提到这个向当地分销的农资商，竟然实现了零库存管理。这是多少企业梦寐以求的事情呀。如果这套管理办法再进行完善，形成一套管理软件，那对我们农资行业可真是一大幸事。

1 基于疲劳驱动力能量的疲劳损伤累积

对于高周疲劳问题，疲劳寿命与应力之间的S—N曲线关系常用Basquin公式描述，即

两口薄棺停放在凉棚下。灵床拆除了，门板上了回去，长明灯已吹灭；唯有俩人的牌位，被摆到客堂里墙的壁几上，依旧点着香烛，供着糕果。桃花和梨花因为连夜哭丧与守灵，俩人站起来都摇摇晃晃的，终于被人劝扶到隔壁叔叔家休息了。梨花躺在床上，整个人像躺在摇晃的船上，始终无法入睡。倒是桃花，一会儿就起了鼾声。

(1)

式中：C为材料固有的疲劳强度常数；b为对应的材料性能参数，为常数；σ为载荷作用下的应力,常用应力幅表示；Nf为相应应力幅对应的疲劳寿命。

由于往复作用的外部载荷是产生疲劳损伤的主要因素，可按式(2)定义疲劳驱动力参数[14]，该参数与载荷循环数、应力幅以及对应的疲劳寿命直接相关。

(2)

式中：σD为疲劳驱动力参数；n为循环载荷作用的次数。

首先，要加大对农业科技投入的力度，加强对农业科学研究的投入。注重农业技术的改造与创新，使各项技术能够因地制宜地发挥应有的作用，同时把自主开发和引进新技术相结合，提高安徽省的农业科技水平。

式(2)定义的参数是与实际载荷循环相关的增量函数,具有非线性特征。

关于高周疲劳问题，可忽略不计很小的塑性应变，依据式(2)定义的疲劳驱动力参数以及能量准则，可进一步定义疲劳驱动力能量参数，并将其简化为

(3)

式中：E为弹性模量；WDn为与n次循环载荷作用对应的疲劳驱动力能量。

由此可见，对意义的正确理解，是口译的核心，译员的知识储备、心理素质等多方面原因会影响其对意义的理解。按释意理论重意义的模式，从各方面研究对意义的理解，以辅助译员正确理解意义，确实是口译研究的一个重要方向。

式(3)定义的疲劳驱动力能量随着循环载荷的作用而不断变化，不断地循环累积，最终达到临界阈值。

总而言之，建筑节能设计十分重要和必要，符合国家战略导向，是绿色生态建设的重要举措，对人类社会生存与发展意义重大，建筑节能设计作为一项系统化工程，应从屋顶节能、外墙节能、门窗节能以及再生能源利用等诸多方面进行考量设计，未来关于此类课题的研究任重而道远。

损伤参数的确定是寻找疲劳损伤累积规律的一个关键问题。结合损伤累积理论以及基于能量耗散考虑，损伤参数D可采用疲劳驱动力能量的增量进行描述，为

(4)

式中：ΔWDn,D0为实际加载状态下耗散的能量；ΔWDc,D0为失效状态下耗散的总能量。

式(4)定义的损伤参数D直接与载荷循环次数相关，且两者之间为非线性关系。以损伤参数D定义的疲劳寿命预测模型则属于非线性疲劳寿命预测模型的范畴。预测模型中载荷与损伤是不可分离的量，在一定程度上可以体现载荷顺序效应。

依据式(4)定义的损伤参数，对于多级载荷作用，可通过迭代计算，得到第i级载荷作用下应力幅对应的疲劳寿命预测计算公式为[15]

会展英语教学要体现后现代理论中的可行性，在教学过程中，要注意培养学生的情感和能力素质，从始业教育、专业作业团队、综合实训、职业技能培训、企业合作项目、职业技能考试、专业比赛等方面，营造职业工作氛围，创造职业工作环境，让学生在不断进行的团队作业中，重视工作责任感的培养和创新能力的培育，提高人际交往能力和沟通能力，同时也能让专业教师及时了解学生职业素质培养的效果。

(5)

其中，

式中：ni为第i级载荷作用下应力幅对应的循环次数；Nfi为第i级载荷单独作用下应力幅对应的疲劳寿命；Nni为第i级载荷作用下应力幅对应的循环次数比；下标p表示对应级别载荷根据预测模型计算得到的预测值。

在利用上述给定的模型、即按式(5)进行疲劳寿命预测时，只需要确定b和Nf这2个基本参数即可，而这2个参数可由试验所得的S—N曲线拟合确定，不需要引入其他参数，计算模式相对明显且直接，具有很好的工程应用性。

境外单位或个人在境内发生应税行为而在境内未设有经营机构的，以购买方为增值税扣缴义务人，《规定》中新增“代扣代交增值税”明细科目核算扣缴义务人代扣代交增值税，实务中扣缴义务人除了代扣代缴增值税外，还有增值税金附加、企业所得税或个人所得税等，按简并原则可将该其改为“代扣代交税金”明细科目。

2 改进的非线性预测模型

对于恒幅载荷作用，疲劳损伤累积规律相对简单，累积规律比较容易探究。而在变幅载荷作用下，疲劳损伤累积规律因载荷幅值的变化而变得复杂。变幅载荷作用下的疲劳损伤累积规律与加载历程密切相关，载荷大小次序、载荷间的干涉、疲劳损伤程度等因素直接影响着损伤累积规律。

部分学者研究表明，载荷间的影响可用相邻两级载荷大小的比值进行描述，比值的大小直接决定交互影响的程度。在多级载荷作用下，利用相邻两级载荷应力比进行修正，对应的疲劳寿命预测模型为

(6)

其中，

式中：σi为第i级载荷对应的应力幅。

式(6)所示预测模型虽然有一定的改进，数值验证结果也有一定的提升，但是该模型仅考虑了载荷先后作用顺序的影响，而未能考虑高低载荷先后作用下实时疲劳损伤对疲劳寿命的影响，损伤累积的演化规律尚需进一步探讨。在变幅载荷作用下，破坏时对应的疲劳损伤累积往往不等于1，前期的低应力具有初始锻炼效应，而前期的高应力则使裂纹提前形成。如果直接仅用应力比表示载荷交互作用的影响，那么实际低应力的锻炼程度以及高应力裂纹的形成状态未能予以体现；可在交互因子中引入载荷循环比，体现载荷循环的作用程度。

按照疲劳损伤等效原理，损伤参数与循环比分数n/Nf的关系均可用幂函数的形式表示。对于多级载荷作用，疲劳寿命对应的循环载荷比分数可以经下式推导求得。

(F((Nni-1)e, σi-1))ωi-2,i-1=

(F(1-Nni,σi))ωi-1,i

(7)

式中：F为描述疲劳损伤状态的函数；ωi-1,i为第i-1级与第i级载荷间的交互因子；(Nni-1)e为前i-1级载荷作用下的等效循环次数比。

有机磷类农药在 DB-17色谱柱上的出峰时间及顺序依次为敌敌畏2.706分钟、氧化乐果8.461分钟、甲基对硫磷11.213分钟，这3种有机磷农药在DB-17色谱柱上未出现出峰时间相近和峰面积重叠的现象；有机氯类农药在 HP-5色谱柱上的出峰时间及顺序依次为乙烯菌核利11.035分钟、联苯菊酯19.157分钟、氯氰菊酯23.628分钟、氰戊菊酯25.637分钟，这4种有机氯类农药能较好地显示在 HP-5色谱柱上。

在式(7)对应的疲劳损伤模型中，体现载荷变幅影响的关键参数为交互因子ω。先高后低变幅加载时ω大于1，反之则ω小于1。在疲劳寿命等效模型中同时考虑载荷大小、加载次序以及载荷间的干涉影响，对应的疲劳寿命预测模型可统一写为

(8)

其中，

γ=ωi-1,i-2ωi-1,i

据上述分析，交互因子的定义是疲劳寿命预测模型的关键，如果将交互因子仅定义为相邻载荷的应力比，式(8)则与式(6)相同。在保证交互因子取值特征的同时，可进一步将相邻两级载荷作用的交互因子定义为应力幅比、载荷循环比或实时损伤相关的函数，即

(9)

式中：αi-1为与载荷循环比以及实时损伤相关的参数。

由上文推导可知，损伤参数可用驱动力能量的变化定义，且该参数与载荷循环比密切相关。据此分析，可以采用下面的2种形式对交互因子的指数参数进行定义，对交互因子进行修正，具体如下。

方法1：将交互因子的指数参数定义为与载荷循环次数比ni/Nfi相关的函数，为

(10)

将式(10)代入式(9),可得到按照方法1修正的疲劳寿命预测模型的交互因子为

(11)

(12)

方法2：将交互因子的指数参数定义为实时疲劳损伤Di的指数函数，为

还有一次去延安采访，那时我从未去过延安，只知道那是个革命圣地，有宝塔山、延河水，还有枣园的灯光、领袖住过的窑洞，但经济尚不发达。当时我住在当地最好的延安宾馆，虽然设施陈旧些，但比我想象中好多了，房间宽敞大气。风尘仆仆采访了一天，有点累，吃过晚饭，看了一会儿无聊的电视，就想早点休息，因为第二天还要早起采访。进了卫生间洗澡，先冲一遍水，然后关上水龙头，打浴液再次冲洗就齐活了。谁知打完浴液准备再次冲澡时，水龙头竟然“罢工”了，无论冷水还是热水踪影皆无。给服务台打电话，方知停水了，第二天上午8时恢复供水。没办法，只好把暖瓶里喝的水倒到脸盆里，晾凉了，蘸湿毛巾，把身上的浴液擦掉了事。

αi=expDi

(13)

式(13)中损伤参数含有载荷循环比值，可将其代入式(9),可得到按照方法2修正的疲劳寿命预测模型的交互因子为

(14)

(15)

在疲劳寿命预测分析时，结合式(8)、式(10)和式(13)进行计算。式(10)和(13)分别代表2种改进的疲劳寿命预测模型的指数参数，模型中既有前后2级载荷的应力幅比，也含有载荷循环比或实时损伤量，能够有效地体现载荷次序的影响、载荷间的干涉以及实时疲劳损伤影响程度。此外，模型中仅含有常用的2个基本参数，通过试验数据比较容易确定，不需要引入其它多余的参量，方便应用于实际结构的疲劳寿命分析。

3 数值算例

为验证所提模型预测结果的可行性和有效性，依据文献中2种常用材料的焊接结构多级载荷加载的试验数据，将所建疲劳寿命预测模型与文献中模型的结果进行对比。

算例1：Q235B焊接接头疲劳寿命预测

Q235B是车辆转向架焊接构架的常用焊接材料，依据文献[12]中不同焊接形式在多级变幅载荷作用下的试验数据验证所建模型的疲劳寿命预测能力。对于不同级数的变幅载荷作用，采用不同分析模型得到的疲劳寿命预测结果见表1—表3,疲劳寿命预测结果与试验结果的对比如图1所示。表中，加载模式是指多级递增应力幅载荷作用时所对应的应力幅值。

由表1—表3和图1可知：对于Q235B焊接接头，在多级递增应力幅载荷作用下，文中所建非线性模型均能有效地进行疲劳寿命预测；相比较于其他模型，所改进模型预测结果更符合于试验数据。依据不同方式交互因子修正的对比结果，可以看到实时损伤修正模式更为合理。

算例2：铝合金焊接接头疲劳寿命预测

铝合金材料是轨道交通领域常用的焊接接头材料，以文献[15]中铝合金试件在多级应力幅载荷作用下的疲劳试验数据验证本文预测模型的疲劳寿命预测能力。多级应力幅试验加载包括由高到低、由低到高以及随机加载3种情况，所对应载荷的应力幅σi分别为305，280，260和240 MPa这4个等级，各级应力幅值对应的单级应力幅疲劳寿命Nfi分别为38 000，87 612，180 660和394 765次。不同预测模型所对应最后一级载荷的疲劳损伤以及疲劳寿命预测结果见表4，疲劳损伤累积结果的对比如图2所示。

这样白净斯文的人做防赤化工作，而且是个头头，一定是个厉害角色。要知道，这类人肚子里永远有使不完的计策和阴谋。刘雁衡揣度自己这一回讨不了好去，打定主意，索性逗他一回。

国家对建筑材料的质量有着严格的规定和划分，个别企业也有自己的材料使用质量标准，但是个别施工单位为了追求额外效益，会有意无意地在工程项目建设过程中使用一些不规范的工程材料，造成工程项目的最终质量存在问题。

由表4和图2可知：对于铝合金焊接接头，在递增、递减以及随机3种类型的变幅载荷作用下，文中改进的非线性疲劳寿命预测模型均能对该焊接接头的疲劳寿命进行有效的预测；与传统的Miner法则以及文献[15]中模型的结果对比可知，改进的非线性预测模型预测精度均有一定的提升，预测结果与实际试验数据更为接近；在文中改进的2种方式中，直接采用实时疲劳损伤对交互因子进行修正的分析结果更为准确一些，实时损伤程度越大，交互因子的作用效果则更为明显。

表1 Q235B焊接接头在3级载荷作用下采用不同模型得到的疲劳寿命预测结果

加载模式应力幅值σi/MPa各级应力幅的循环次数ni/(104次)单级应力幅作用的疲劳寿命Nfi/(104次)试验寿命NE/(104次)采用不同模型时的预测寿命Nf/(104次)Miner模型文献[12]预测模型本文预测模型1本文预测模型238.0040.00107.72140.0040.0086.40104.7588.20107.66100.8492.9842.0024.7570.0540.0040.0086.40242.0040.0070.05107.6371.0185.9790.9876.5944.0027.6357.3520.6049.60179.77321.5050.30144.15110.58151.44181.77200.33166.0122.4010.68116.64

表2 Q235B焊接接头在4级载荷作用下采用不同模型得到的疲劳寿命预测结果

加载模式应力幅值σi/MPa各级应力幅的循环次数ni/(104次)单级应力幅作用的疲劳寿命Nfi/(104次)试验寿命NE/(104次)采用不同模型时的预测寿命Nf/(104次)Miner模型文献[12]预测模型本文预测模型1本文预测模型2120.6049.60179.7721.5050.30144.1522.4040.60116.6423.3039.9495.17180.44129.57152.90142.75142.37228.9540.00146.4730.4840.00119.0232.0040.0097.7035.058.2467.62128.24112.45184.71129.27127.10330.4840.00119.0233.5240.0080.9435.0540.0067.6238.1014.0548.26134.0578.26123.62103.71100.79

表3 Q235B焊接接头在5级载荷作用下采用不同模型得到的疲劳寿命预测结果

加载模式应力幅值σi/MPa各级应力幅的循环次数ni/(104次)单级应力幅作用的疲劳寿命Nfi/(104次)试验寿命NE/(104次)采用不同模型时的预测寿命Nf/(104次)Miner模型文献[12]预测模型本文预测模型1本文预测模型2122.4025.03116.6423.3053.0795.1725.1044.3564.8126.0021.5654.0330.9014.0422.15158.0563.49110.43117.94203.67232.7640.0088.8534.2940.0073.8735.8140.0061.9837.3340.0052.3938.8613.0344.53173.0364.2593.40121.1151.61332.0040.0097.7033.5240.0080.9835.0540.0067.6036.5740.0056.9338.1092.8748.24252.8761.3375.43168.24167.82

图1 多级载荷作用下Q235B焊接接头疲劳寿命结果对比

图2 4级载荷作用下铝合金6082T6接头疲劳损伤结果对比

表4 铝合金6082T6在4级载荷作用下采用不同模型得到的疲劳寿命预测结果

加载应力幅/MPa各级应力幅的循环数ni试验剩余寿命循环比(n/Nf)Miner模型文献[15]预测模型本文预测模型1本文预测模型2循环数ni循环比n/Nf循环数ni循环比n/Nf循环数ni循环比n/Nf循环数ni循环比n/Nf递减σ1=305n1=10 950σ2=280n2=19 427σ3=260n3=26 258σ4=240n4=52 5000.133136 1900.345103 0300.26191 5850.23288 1030.223递增σ1=240n1=103 000σ2=260n2=26 258σ3=280n3=19 427σ4=305n4=16 8000.44214 1360.37217 5940.46314 5160.38214 9340.393随机σ1=280n1=19 427σ2=305n2=10 950σ3=260n3=26 258σ4=240n4=43 4000.110136 1900.345106 9800.27191 9800.23386 3500.219

注：表中第3列各级应力幅的循环数所对应第4级应力幅的循环数n4为试验值。

4 结 论

(1)利用变幅载荷的应力比、循环比以及实时损伤状态参量对相邻载荷间交互因子进行修正，所改进的非线性疲劳寿命预测模型能较好地处理变幅载荷加载的历史效应，在不需要引入额外参数的前提下，可有效地处理变幅载荷加载历史对损伤累积的影响，能够较为准确地对焊接结构进行疲劳寿命预测。

1.首先应加强饲养管理，严格卫生消毒措施，注意通风换气，保持舍内空气清新。减少各种应激因素的影响，保持猪群足够均衡的营养水平。

(2)根据Q235B焊接接头的疲劳寿命对比结果，在不同级数的变幅载荷作用下，改进模型均具有可行性，在精度方面也有不同程度的提高，预测分析结果更接近于试验结果。

(3)根据铝合金焊接接头的疲劳寿命对比结果，在递增、递减以及随机的变幅载荷作用下，改进的非线性疲劳寿命预测模型均能对焊接接头疲劳寿命进行有效的预测。与文献模型以及传统的Miner模型相比，预测精度有明显的提升，与文献中所给出的试验结果更加接近。

(4)改进的非线性疲劳寿命预测模型涉及的参数仅为2个，且比较容易通过试验确定，计算递推格式也较为简单，不需要多层嵌套，非常便于应用于实际焊接结构的疲劳寿命预测。

参考文献

[1] 王红，刘万选，商跃进. 基于热点应力的国产化Y25型转向架焊接构架疲劳强度分析[J].中国铁道科学,2013,34(4):66-70.

(WANG Hong, LIU Wanxuan, SHANG Yuejin. Fatigue Strength Analysis on Welded Frame of Domestic Y25 Bogie Based on Hot Spot Stress [J].China Railway Science, 2013, 34(4):66-70. in Chinese)

[2] 薛齐文，杜秀云，王生武.基于载荷加载次序的疲劳寿命预测改进模型[J].中国铁道科学,2019,40(1):88-93.

(XUE Qiwen, DU Xiuyun, WANG Shengwu. An Improved Fatigue Life Prediction Model Based on Loading Sequence [J].China Railway Science, 2019,40(1):88-93. in Chinese)

[3] 刘德昆，李强，林浩博，等.高速列车转向架构架结构的疲劳可靠性模型[J].中国铁道科学,2017,38(5):100-106.

(LIU Dekun, LI Qiang, LIN Haobo, et al. Fatigue Reliability Model of Bogie Frame Structure for High-Speed Train [J].China Railway Science,2017,38(5):100-106. in Chinese)

[4] YUE P, LEI Q, ZHANG C L, et al. A Modified Cumulative Damage Model for Fatigue Life Prediction under Variable Amplitude Loadings [J].Applied Mechanics & Materials, 2016, 853:62-66.

[5] YUAN R, LI H, HUANG H, et al. A New Non-Linear Continuum Damage Mechanics Model for Fatigue Life Prediction under Variable Loading [J].Mechanics, 2013, 19(5): 506-511.

[6] DJEBLI A, AID A, BENDOUBA M, et al. A Non-Linear Energy Model of Fatigue Damage Accumulation and Its Verification for Al-2024 Aluminum Alloy [J].International Journal of Non-Linear Mechanics, 2013, 51:145-151.

[7] GAO H, HUANG H Z, LV Z, et al. An Improved Corten-Dolan’s Model Based on Damage and Stress State Effects [J].Journal of Mechanical Science & Technology, 2015, 29(8):3215-3223.

[8] BENKABOUCHE S, GUECHICHI H, AMROUCHE A, et al. A Modified Nonlinear Fatigue Damage Accumulation Model under Multiaxial Variable Amplitude Loading [J].International Journal of Mechanical Sciences, 2015, 100:180-194.

[9] LAGODA T. Energy Models for Fatigue Life Estimation under Uniaxial Random Loading. Part I: The Model Elaboration [J].International Journal of Fatigue, 2001, 23(6):467-480.

[10] KWOFIE S, RAHBAR N. A Fatigue Driving Stress Approach to Damage and Life Prediction under Variable Amplitude Loading [J].International Journal of Damage Mechanics, 2013, 22(3):393-404.

[11] 彭兆春. 基于疲劳损伤累积理论的结构寿命预测与时变可靠性分析方法研究[D]. 成都: 电子科技大学, 2017.

(PENG Zhaochun. Research on Methods for Structural Life Prediction and Time-Dependent Reliability Analysis Using Cumulative Fatigue Damage Theories [D].Chengdu: University of Electronic Science and Technology of China, 2017. in Chinese)

[12] 高会英.复杂应力状态下焊接接头的疲劳寿命预测方法研究[D]. 成都: 电子科技大学, 2016.

(GAO Huiying. Research on Fatigue Life Prediction Methods of Welded Joints under Complex Stress States [D].Chengdu: University of Electronic Science and Technology of China, 2016. in Chinese)

[13] AID A, AMROUCHE A, BOUIADJRA B B, et al. Fatigue Life Prediction under Variable Loading Based on a New Damage Model [J].Materials & Design, 2011, 32(1):183-191.

[14] El Aghoury I, GALAL K. A Fatigue Stress-life Damage Accumulation Model for Variable Amplitude Fatigue Loading Based on Virtual Target Life [J].Engineering Structures, 2013, 52:621-628.

[15] PENG Z, HUANG H Z, ZHU S P, et al. A Fatigue Driving Energy Approach to High-Cycle Fatigue Life Estimation under Variable Amplitude Loading [J].Fatigue & Fracture of Engineering Materials & Structures, 2016, 39(2):180-193.