本发明属于飞机强度设计技术领域，具体涉及一种基于局部应力-应变法的疲劳寿命校准方法。
背景技术：
当前，对飞机疲劳关键部位进行随机谱下的寿命分析时，由于随机谱中存在复杂的载荷交互作用，因此目前尚未有一个普遍适用的高精度损伤累积模型能较好地预测其疲劳寿命。
但通过研究发现，系列化飞机往往在结构和载荷谱上具有很大的继承性，因此可以采用已有型号的疲劳试验结果，对后续型号的疲劳寿命进行预测。
技术实现要素：
本发明的目的：为了解决上述问题，本发明提出了一种基于局部应力-应变法的疲劳寿命校准方法，可在变换随机载荷谱的情况下，仍保证计算所得的疲劳寿命具有较高的精度。
本发明的技术方案：一种基于局部应力-应变法的疲劳寿命校准方法，所述局部应力-应变法的εa-2nf曲线的表达式：
式中，εa为应变幅值，σ'f为疲劳强度系数，b为疲劳强度指数，ε'f为疲劳延性系数，c为疲劳延性指数，2nf为到破坏的反复数；
包括以下步骤：
步骤一、确定局部应力-应变法中的最敏感参数；
所述εa-2nf曲线包括5个参数：εa、σ'f、b、ε'f、c；
a)对εa进行优化，相当于在εa-2nf曲线表达式左侧乘以一个系数k，相关敏感参数有kf、σ'f、b、ε'f、c；
b)由于局部应力-应变法适用于应力集中较大及应力集中区较大范围进入屈服的情况，弹性段性能参数不应作为主要考虑的因素，因此，在上述参数中取kf为最敏感的参数；
步骤二、对kf进行修正，以期实现对疲劳性能曲线参数的修正；
c)根据预先指定的载荷应力谱，确定其对应的中值寿命t50为目标进行优化；
d)设定初始值，并不断调整kf的取值，反复用局部应力-应变法计算中值寿命，直至计算寿命等于设定的目标中值t50，此时对应的kf值即为修正后的s-n曲线三参数式中的kf值；
步骤三、对确定的局部应力-应变法及其对应kf值进行验证；
利用优化得到的参数和局部应力-应变法进行载荷谱下的寿命预测并将预测计算的该载荷谱下的计算的中止寿命t'与所述载荷谱下的实测寿命中值比较误差；
若误差小于20％,则确定该校准方法。
本发明技术方案的技术效果：本发明一种局部应力-应变法的疲劳寿命校准方法，利用该校准方法优化得到的参数和名义应力法进行其他载荷谱下的寿命预测，使得其他载荷谱下的寿命计算结果仍具有较高的精度。
附图说明
图1为本发明一种局部应力-应变法的寿命校准方法的一优选实施例的流程示意图。
具体实施方式
为使本发明实施的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行更加详细的描述。在附图中，自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本发明一部分实施例，而不是全部的实施例。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。下面结合附图对本发明的实施例进行详细说明。
疲劳寿命校准方法的核心是，以已有的试验结果为基础，采用合理的疲劳分析方法，确定疲劳分析的敏感参数，通过调整参数使得计算寿命和试验寿命相同，然后采用该参数和相同的寿命分析方法进行其他载荷谱下的寿命预测，使得其他载荷谱下的寿命计算结果具有较高的精度。
局部应力-应变法将作用于结构细节的名义应力谱，通过缺口弹塑性分析，转换为结构细节危险点的局部应力-应变历程，然后通过当量循环的方法，把局部谱用雨流计数得到的循环等效为光滑试件对称循环，最后由光滑试件的应变-寿命曲线估算结构危险点的损伤，进而预测结构的疲劳寿命。而疲劳累积损伤理论通常采用miner线性累积损伤理论，该理论是理论界和工程界最普遍接受的，因此提高疲劳寿命预测结果精度的核心在于细节疲劳应变-寿命曲线(εa-2nf曲线)的确定。εa-2nf曲线的表达式如下：
式中，εa为应变幅值；σ'f为疲劳强度系数；b为疲劳强度指数；ε'f为疲劳延性系数；c为疲劳延性指数；2nf为到破坏的反复数。上式可知εa-2nf曲线包括5个参数：εa、σ'f、b、ε'f、c。
其中εa相当于疲劳强度，对εa进行优化相当于在εa-2nf曲线表达式左侧乘以一个系数k，而k相当于在对局部应力应变进行求解时，直接在求得的应变值上乘以一个系数k，因此，相当于对疲劳缺口系数kf进行优化。于是，相关敏感参数有kf、σ'f、b、ε'f、c。
局部应力-应变法中所使用的εa-2nf曲线为材料特性参数曲线，包含的参数较多，在优化参数时，参数越多，需要的样本数量越多，且优化结果的稳健性越差。考虑到局部应力-应变法适用于应力集中较大、应力集中区较大范围进入屈服的情况，弹性段性能参数不应作为主要考虑的因素。因此，在上述参数中取kf为最敏感的参数，在获得局部应力-应变场的过程中综合考虑细节几何、表面加工质量等因素的影响，对kf进行修正，以期实现对疲劳性能曲线参数的修正。
在优化kf值时，可采用某个谱下试验寿命反推确定。反推方法为：依据指定载荷(应力)谱，以其对应的中值寿命t50为目标，假定kf的初始值，并以一定步长不断调整kf值，反复用局部应力-应变法计算中值寿命，直至计算寿命等于t50，计算得到kf值。
利用优化得到的参数和局部应力-应变法进行载荷谱下的寿命预测并将预测计算的该载荷谱下的计算的中止寿命t'与所述载荷谱下的实测寿命中值比较误差，若误差小于20％,则确定该校准方法。
经验证利用该校准方法优化得到的参数和名义应力法进行其他载荷谱下的寿命预测，使得其他载荷谱下的寿命计算结果仍具有较高的精度。
如图1所示：具体实施如下：
1)原始数据的准备，包括载荷谱、材料的循环应力—应变曲线、应变—疲劳曲线和危险部位的疲劳缺口系数；
2)通过缺口弹塑性分析，将作用于结构细节的名义应力谱转换为局部应力—应变历程；
3)进行“雨流”计数统计处理，获得局部应变谱所有全循环和半循环；
4)计算当量应变，将所有全循环或半循环的真实应变范围转换为光滑试件对称循环的当量应变范围；
5)通过累积损伤理论计算中值寿命。
为使εa-2nf曲线反映实际细节的原始疲劳质量，拟以典型细节模拟试件在实际载荷谱下的成组疲劳试验结果为依据，从相关材料手册中查找与结构关键部位相同材料的εa-2nf曲线参数取为初始值，采用局部应力-应变疲劳分析方法进行试验谱下的疲劳寿命分析，以计算寿命与试验寿命相同为目标，不断优化εa-2nf曲线参数；
局部应力-应变法中所使用的εa-2nf曲线包括5个参数：εa、σ'f、b、ε'f、c。其中εa相当于疲劳强度，对εa进行优化相当于在εa-2nf曲线表达式左侧乘以一个系数k，而k相当于在对局部应力应变进行求解时，直接在求得的应变值上乘以一个系数k，因此，相当于对疲劳缺口系数kf进行优化。于是，相关敏感参数有kf、σ'f、b、ε'f、c。
考虑到局部应力-应变法适用于应力集中较大、应力集中区较大范围进入屈服的情况，弹性段性能参数不应作为主要考虑的因素。因此，在上述参数中取kf为最敏感的参数。
在优化kf值时，可采用某个谱下试验寿命反推确定。反推方法为：依据指定载荷(应力)谱，以其对应的中值寿命t50为目标，假定kf的初始值，并以一定步长不断调整kf值，反复用局部应力-应变法计算中值寿命，直至计算寿命等于t50，计算得到kf值。
利用优化得到的参数和局部应力-应变法进行载荷谱下的寿命预测并将预测计算的该载荷谱下的计算的中止寿命t'与所述载荷谱下的实测寿命中值比较误差，若误差小于20％,则确定该校准方法。
经验证利用该校准方法优化得到的参数和名义应力法进行其他载荷谱下的寿命预测，使得其他载荷谱下的寿命计算结果仍具有较高的精度。
最后需要指出的是：以上实施例仅用以说明本发明的技术方案，而非对其限制。尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。