一种基于多因素融合修正的结构件剩余寿命预测方法
【技术领域】
[0001] 本发明属于机械检测领域，尤其是一种结构件裂纹扩展预测方法。
【背景技术】
[0002] 随着社会的进步，大量设备都朝着大型化的方向发展。在这样的背景前提下，微小 的故障就会造成巨大的损失，因此对设备或构件的剩余寿命进行预测是必要且紧迫的。随 着科学技术的发展，人们对疲劳这一学科的认识更加深入，一系列疲劳寿命预测方法被得 到广泛应用：如名义应力法、局部应力-应变法、临界区域法等。临界区域法是一类具有工 程应用前景的疲劳寿命估算方法。应力场强法和临界距离法作为临界区域法中常见的两类 方法已得到一些应用。工程实际中则更倾向于使用名义应力法和局部应力-应变法。传统 名义应力法的最大优点是分析步骤简单、适用面广且实用性强。局部应力-应变法以应变 分析为基础，用应变寿命e-N曲线代替S-N曲线并配以循环加载下的应力-应变曲线来预 测疲劳寿命。这些现有的预测技术主要存在以下不足：（1)应力场强法和临界距离法大多 只针对平面结构，而以结构件作为分析对象的还很少（2)对于传统名义应力法来说，由于 结构件与材料试样间的疲劳性能参数存在差异，所以疲劳寿命的估算精度无法得到保证。
[3] 局部应力-应变法的估算精度较高，但其适用面较窄。
[0003] 故，需要一种新的技术方案以解决上述问题。
【发明内容】
[0004] 为了克服现有技术中存在的问题，本发明提供了一种基于多因素融合修正的结构 件剩余寿命预测方法。
[0005] 为达到上述目的，本发明可采用如下技术方案：
[0006] 一种基于多因素融合修正的结构件剩余寿命预测方法，其特征在于：该方法包括 以下步骤：
[0007] S1?结构件应力强度的计算：
[0008] 建立结构件的三维模型，并运用有限元分析软件对三维模型的应力强度进行分 析；
[0009] S2?残余应力的测量与修正：
[0010] 对实体的结构件的残余应力进行测量以获得其残余应力的分布状态，并将测量结 果与S1的计算结果一起代入修正公式完成修正；
[0011] S3.基于分段原理的应力比计算：
[0012] 根据采集的实际工况下的非规律性随机载荷历程，采用分段原则计算出应力比；
[0013] S4.尺寸修正参数的计算：
[0014] 根据S1的计算结果及材料试样的应力状态分布计算结构件所对应的尺寸修正参 数；
[0015] S5?应力集中修正参数的计算：
[0016] 根据S1的计算结果选取合适的应力场积分路径，在等效于应力缓释区域的最佳 积分路径下计算应力集中修正参数；
[0017] S6.表面制造加工质量修正参数的选取：
[0018] 确定工艺类型，并根据实测的表面硬度选取出与结构件所对应的表面制造加工质 量修正参数；
[0019] S7.基于多因素修正的结构件剩余寿命预测：
[0020]将S1-S6的计算结果代入基于多因素修正的结构件剩余寿命预测算法来完成结 构件剩余寿命的预测。
[0021] 与现有技术相比，本发明能够将所有能够对结构件疲劳寿命造成影响的众多因素 考虑在内，并通过上述操作步骤从而可以有效地避免现有技术存在的考虑因素少的弊端， 预测结果精度得到大幅度提高，从而实现结构件剩余寿命的准确有效预测。
【附图说明】
[0022] 图1是本发明基于多因素融合修正的结构件剩余寿命预测方法的流程框图。
【具体实施方式】
[0023] 请参阅图1所示，本发明公开了一种基于多因素融合修正的结构件裂纹扩展预测 方法。
[0024] 该方法包括以下步骤：
[0025] 一种基于多因素融合修正的结构件剩余寿命预测方法，其特征在于：该方法包括 以下步骤：
[0026] S1.构件应力强度的计算
[0027] 根据图纸建立分析构件的三维模型，将建立完毕后的三维模型导入到有限元软件 中，按网格划分、约束设置、载荷施加步骤完成构件有限元分析的前处理，待前处理完毕后 运用有限元分析软件对其应力强度进行分析，以获取其应力状态分布；
[0028] S2?残余应力的测量与修正
[0029] 对实体构件的残余应力进行测量以获得结构件所对应的残余应力分布状态，并将 测量结果与S1的计算结果一起代入修正公式1完成修正；
【主权项】
1. 一种基于多因素融合修正的结构件剩余寿命预测方法，其特征在于：该方法包括以 下步骤： S1?结构件应力强度的计算： 建立结构件的三维模型，并运用有限元分析软件对三维模型的应力强度进行分析； S2?残余应力的测量与修正： 对实体的结构件的残余应力进行测量以获得其残余应力的分布状态，并将测量结果与S1的计算结果一起代入修正公式完成修正；
53. 基于分段原理的应力比计算： 根据采集的实际工况下的非规律性随机载荷历程，采用分段原则计算出应力比；
54. 尺寸修正参数的计算： 根据S1的计算结果及材料试样的应力状态分布计算结构件所对应的尺寸修正参数； S5?应力集中修正参数的计算： 根据S1的计算结果选取合适的应力场积分路径，在等效于应力缓释区域的最佳积分 路径下计算应力集中修正参数；
56. 表面制造加工质量修正参数的选取： 确定工艺类型，并根据实测的表面硬度选取出与结构件所对应的表面制造加工质量修 正参数；
57. 基于多因素修正的结构件剩余寿命预测： 将S1-S6的计算结果代入基于多因素修正的结构件剩余寿命预测算法来完成结构件 剩余寿命的预测。
2. 根据权利要求1所述的基于多因素融合修正的结构件剩余寿命预测方法，其特征在 于：所述的步骤S2中，包括以下具体工序： S2. 1、据S1的计算结果，根据应力幅值的大小决定疲劳寿命数值的关系来确定结构件 中需要进行重点测量残余应力的区域范围； S2. 2、对2. 1中已确定的区域进行电解抛光； S2. 3、对S2. 2已处理完毕的区域进行残余应力的测量； S2. 4、对S2. 3的结果进行数据统计以合成三维残余应力分布图，并将分析所得到的三 维应力分布结果与S1的分析结果一起代入修正公式来完成修正过程。
3. 根据权利要求1所述的基于多因素融合修正的结构件剩余寿命预测方法，其特征在 于： S2中，所述修正公式为：
〇 "为基于残余应力修正后的应力幅值；〇 _为最大应力；r为应力比；〇 bb为抗拉强 度；〇K为对应的分析强度；w为富裕应力缓释放敏感系数。
4. 根据权利要求1所述的基于多因素融合修正的结构件剩余寿命预测方法，其特征在 于： S3中，采用分段原则计算出应力比的方法为： 对于一个存在（ti,oi)对应关系的载荷谱来说，其平均应力〇ni可表示为：
其中载荷谱表示为0iifUi)，令f(t) = 〇m，可得对应于平均应力〇m的时间分量tq,q= 1, 2, 3......m； 由于载荷的循环特性，按照载荷谱上的时间顺序在时间轴上依次选取两个点，则根据 载荷谱的循环特性对其进行分段后则形成（Uh)，（t2，t4)，（t5，t7)......;即在分段数据 中至少存在两点，使得。，(^且〇 待分段完毕后得应力比r
〇
5. 根据权利要求1所述的基于多因素融合修正的结构件剩余寿命预测方法，其特征在 于： S4中，尺寸修正参数为
结构件尺寸修正参数为结构件与参照标准试样间应力场函数积分的比值； 应力场函数也可通过应力场中应力积分路径下某点距离最大局部应力根部的距离L(i)来表示，即 f(x1；x2) =A!+A2L(i) +A3L2 (i) +A4L3 (i) 式中，ApA2,A3,A4及BpB2,B3,B4表示拟合函数表达式中的系数。
6. 根据权利要求1所述的基于多因素融合修正的结构件剩余寿命预测方法，其特征在 于： S5中，应力集中修正参数的计算具体求解过程包括： 1) 、根据S1的计算结果，沿最大应力的横截面选取反映应力场分布的积分路径； 2) 、在积分路径上取点，即可获得应力根部截面方向上点距L与应力s之间的对应值， 通过拟合上述数据即可获得对应的应力场函数表达式； 3) 、将应力场函数表达式代入下式，求解与应力场所对应的名义应力：
4) 将名义应力代入下式，获得结构件所对应的应力集中修正参数：
7. 根据权利要求1所述的基于多因素融合修正的结构件剩余寿命预测方法，其特征在 于： 结构件剩余寿命的预测通过下式得到
其中V为基于残余应力修正后的应力幅值；Kt为应力集中修正参数；e为尺寸修正 参数；0为表面质量修正参数；N为疲劳寿命；m，C为与材料、应力比有关的参数。
【专利摘要】本发明公开一种基于多因素融合修正的结构件剩余寿命预测方法。对结构件剩余寿命造成影响的众多因素考虑在内，包括结构件的应力强度、残余应力的测量与修正、应力比计算、尺寸修正参数的计算、应力集中修正参数、表面制造加工质量修正参数的选取等，使预测结果精度得到大幅度提高，从而实现结构件剩余寿命的准确有效预测。
【IPC分类】G06F17-50
【公开号】CN104809311
【申请号】CN201510247506
【发明人】贾民平, 朱林, 许飞云, 胡建中, 黄鹏
【申请人】东南大学
【公开日】2015年7月29日
【申请日】2015年5月14日