一种齿轮接触疲劳可靠性分析方法技术

本发明专利技术公开了一种齿轮接触疲劳可靠性分析方法。该方法首先建立了齿轮的接触应力模型，并基于疲劳S‑N曲线获取了一定周期后的疲劳剩余强度，进而构建了齿轮接触疲劳失效的功能函数。然后，根据随机参数的分布特性，运用蒙特卡洛方法产生点集，并从中随机选出N个点作为实验设计点构建了初始Kriging模型，再通过主动学习策略，获得了最接近极限状态面或误差最大的样本点，并将其作为最佳样本点增加到原有实验设计点中，更新了Kriging模型，直至收敛。根据收敛的Kriging模型，获取了所有样本点集对应的功能函数近似值，并采用蒙特卡洛方法计算了其失效概率。最后，通过一具体算例验证了该方法的可行性。该发明专利技术的齿轮接触疲劳可靠性分析方法更具高效性和实用性。

【技术实现步骤摘要】
一种齿轮接触疲劳可靠性分析方法
本专利技术涉及一种可靠性分析方法，尤其涉及基于主动学习策略的齿轮接触疲劳可靠性分析方法。
技术介绍
齿轮往往作为装备的动力传输的主要载体，而齿面的接触疲劳引起的损伤破坏是最普遍的一种疲劳失效模式，约占报废齿轮的70％，针对齿轮进行齿面接触疲劳分析对装备的安全运行非常必要。传统的齿轮接触强度设计通常采用安全系数作为衡量标准。由于该安全系数的确定过程涉及人为主观因素，经常会造成结构粗笨，甚至存在安全性不足的隐患，已逐渐不能满足现代工程设计的需要。由于可以精确度量和跟踪结构服役过程中的不确定性因素，基于可靠性的结构强度分析和设计方法可保证结构具有足够的安全性，并且满足经济性需求。因此，实际工程中，非常有必要对齿轮结构进行接触疲劳可靠性分析，进而为齿轮疲劳强的设计提供更加科学而合理的依据。传统的齿轮接触可靠性分析方法通常依赖于试验，其成本高且周期长，当齿轮材料、模数及齿数等参数发生变化时，试验必须重新开展，需花费大量的成本和时间。为了有效避免试验、提高开发周期、节省研发成本，研究者们开始对齿轮结构进行有限元建模，发展基于物理模型的齿轮接触疲劳可靠性分析方法。目前，在基于物理模型的接触疲劳可靠性分析方法中，蒙特卡洛仿真方法(MCS)是最常见的方法，但该方法需调用大量的有限元分析模型，因计算成本过高而难以应用于实际工程。虽然，后续发展的最可能失效点(MPP)方法和矩方法等，在一定程度上解决了可靠性分析的效率问题，但其针对参数不确定性和非线性问题依然存在挑战。事实上，齿轮接触疲劳问题中不确定性参数的分散性通常较大，且随着疲劳循环次数的增加，疲劳寿命与参数之间的关系也变得极端非线性。因此，针对齿轮结构，迫切需要发展一种高效而稳健的接触疲劳可靠性分析方法。
技术实现思路
本专利技术提出了一种基于主动学习策略的齿轮疲劳接触可靠性分析方法，通过少量调用齿轮接触模型即可得到精确的可靠度，这对有效保障齿轮结构在服役过程中的可靠性、稳定性和经济性具有重要的工程意义。附图说明图1是齿轮接触疲劳疲劳可靠性分析流程图；图2是齿轮传动模型；图3是20CrMnTi的接触疲劳S-N曲线；图4是基于主动学习策略的实验设计图；图5是疲劳失效概率随循环周次的变化趋势。具体实施方式由于齿面接触力的循环作用，齿轮结构将发生接触疲劳破坏。因制造误差和环境因素等的影响，齿轮结构在几何特性、材料属性和载荷等方面均存在参数不确定性。本专利技术专利可考虑这些不确定性参数的影响，对齿轮结构进行高效的疲劳可靠性分析。实际工程中，影响齿面接触应力的因素众多，对不同的工况、系统效应、齿轮精度等而言，齿面接触应力往往相差较大。为了方面，本专利技术以渐开线直齿轮为例进行说明。根据美国ANSI/AGMA2001-D04标准，其齿面接触应力σH可计算如下：式中：ZE为弹性系数；Ft为传递的切向载荷；Ko为超载系数；Kv为动载系数；Ks为尺寸系数；KH为载荷分配系数；d为主动轮的工作节圆直径；b为齿宽；ZR为抗点蚀的表面状况系数；ZI为抗点蚀的几何系数。切向载荷Ft齿轮啮合时所承受的切向载荷取决于齿轮所传递的力矩或功率，该作用力作用于齿轮端面上并与分度圆相切，其计算公式如下：式中：T为力矩。对于任意复杂的齿轮结构，可利用有限元方法对齿轮结构进行分析，获取其齿面接触应力。由于不确定性因素的存在，齿面的接触应力不再为一确定性值，而是依赖于不确定性参数的变量。因此，针对齿轮的失效分析需考虑这些不确定性参数对接触应力的影响。齿轮结构接触疲劳可靠性表示在给定时间周期内结构不会因齿面循环接触和磨损而发生失效的概率，其可表示为：Pf＝P(g(x)＜0)\*MERGEFORMAT(3)式中：x为一q维向量，代表齿轮结构的几何尺寸、材料特征和载荷等不确定性因素。例如，若动载系数Kv和载荷分配系数KH为不确定性变量，则x＝[Kv,KH]；g(x)为功能函数：g(x)＝σf-σH(x)\*MERGEFORMAT(4)式中：σf表示齿轮运行一段时间后的疲劳剩余强度，σH为公式(1)中的齿面接触应力。为了获取式中的失效概率，下面将首先建立功能函数g(x)的Kriging模型，然后根据主动学习策略选择样本更新点，最终对齿轮结构进行接触疲劳可靠性分析。基于Kriging模型，齿轮疲劳失效的功能函数式(4)可近似表示为：式中：h(x)表示回归模型的基函数，代表了模型的趋势；β为趋势系数向量。对于通常情况下，h(x)可取为单位向量，β可简单取为当前训练样本点集处功能函数的均值向量。另外，Z(x)为一零均值的平稳高斯随机过程，表示偏离趋势的程度，其任意两点a和b之间的协方差函数可表示为：式中：为高斯过程的方差，R(a,b)为相关系数函数：式中：q代表变量的维数，θi为尺度参数，ai和bi分别表示a和b的第i个坐标。Kriging模型的均值和方差可被估计为：式中：R表示训练样本点之间的协方差函数；r表示x与训练样本点之间的协方差函数；g表示对应于训练样本点集的功能函数向量，F表示n×1的单位向量，n为训练样本点数。最终，通过以下最大似然估计可得到参数θi和σZ：式中：获得参数θi和σZ后，将其带入式(6-7)得到高斯随机过程Z(x)，进一步可建立如式(5)所示的齿轮接触疲劳失效Kriging模型，同时根据式(8)和(9)可继续得到Kriging模型在任意点处的均值和方差。主动学习策略首先运用蒙特卡洛方法产生点集S，并从S中随机选出N个点作为实验设计点构建初始Kriging模型。由于构建Kriging模型时，仅调用了少量齿轮接触疲劳功能函数，其精度不能得到保证。因此，如何从S中确定最好的下一个点来更新模型，进而构建精确的Kriging模型非常重要。下面本专利技术将引入一种主动学习策略来选取最佳样本点。该学习策略认为最影响失效概率的点应该被重点考虑，这些点包括极限状态面附近的点、模型中不确定度最大的点及两种情况都存在的点。基于此，该策略定义如下主动学习函数：式中，和分别为Kriging模型的均值和方差。根据公式和，可得到所有S中样本点对应的U值。然后，选取集合中对应U最小值的样本点作为新实验样本点，更新的实验设计点并计算新的Kriging模型。当U≥2时，Kriging模型被认为是足够准确的，停止学习。可靠性计算及其收敛准则当学习满足停止条件后，根据公式(5)中建立的齿轮接触疲劳Kriging模型预测点集S中所有样本点对应的功能函数并通过点集S中具有负的的数目与S的样本点总数之比来获得齿轮失效的概率的：上述失效概率Pf可根据下面的收敛准则判断是否收敛。根据MCS原理，Pf的变异系数可计算如下：通常，当变异系数低于5％时，Pf的结果精度是可以接受的。若变异系数过高，则需进行蒙特卡洛模拟新的点集合，将其加到原集合S中去，并重新进行上述分析。如果方差系数在允许的范围内，则结束，最后一次失效概率的计算作为最终的齿轮接触疲劳失效概率。如图1所示，本专利技术方法的流程可概括如下：(1)根据实际工况，确定齿轮接触疲劳模型中的不确定性参数及其分布信息，基于疲劳S-N曲线获取一定周期后的疲劳剩余强度，建立接触疲劳可靠性分析的功能函数；(2)根据齿轮接触疲劳模型中的不确定性参数分布特性，随机产生1e6个MCS样本集，并从本文档来自技高网...

【技术保护点】
1.一种齿轮接触疲劳可靠性分析方法，其特征在于，包括以下步骤：(1)根据实际工况，确定齿轮接触疲劳模型中的不确定性参数及其分布信息，基于疲劳S‑N曲线获取一定周期后的疲劳剩余强度，建立基于齿面接触应力的齿轮接触疲劳失效功能函数；(2)根据齿轮接触疲劳模型中的不确定性参数分布特性，随机产生1e6个MCS样本集，并从中选出n个样本点作为初始实验设计点，建立所述齿轮接触疲劳失效功能函数的Kriging模型；(3)根据构建的Kriging模型，获取所有样本点集对应的功能函数近似值，计算其失效概率；(4)根据主动学习策略评估S中的最佳点，将该样本点增加到原有实验设计点中，更新齿轮接触疲劳失效的Kriging模型，判断是否满足停止学习条件：若不满足停止条件，继续学习，返回步骤3；若满足停止条件，停止学习，跳到步骤5；(5)根据收敛准则判断齿轮失效概率是否收敛；若不收敛，在初始实验设计点中增加一个新的MCS样本点，返回步骤2；若收敛，获得最终的齿轮疲劳接触失效概率，结束。

【技术特征摘要】
1.一种齿轮接触疲劳可靠性分析方法，其特征在于，包括以下步骤：(1)根据实际工况，确定齿轮接触疲劳模型中的不确定性参数及其分布信息，基于疲劳S-N曲线获取一定周期后的疲劳剩余强度，建立基于齿面接触应力的齿轮接触疲劳失效功能函数；(2)根据齿轮接触疲劳模型中的不确定性参数分布特性，随机产生1e6个MCS样本集，并从中选出n个样本点作为初始实验设计点，建立所述齿轮接触疲劳失效功能函数的Kriging模型；(3)根据构建的Kriging模型，获取所有样本点集对应的功能函数近似值，计算其失效概率；(4)根据主动学习策略评估S中的最佳点，将该样本点增加到原有实验设计点中，更新齿轮接触疲劳失效的Kriging模型，判断是否满足停止学习条件：若不满足停止条件，继续学习，返回步骤3；若满足停止条件，停止学习，跳到步骤5；(5)根据收敛准则判断齿轮失效概率是否收敛；若不收敛，在初始实验设计点中增加一个新的MCS样本点，返回步骤2；若收敛，获得最终的齿轮疲劳接触失效概率，结束。2.如权利要求1所述的一种齿轮接触疲劳可靠性分析方法，其特征在于，基于Kriging模型，齿轮疲劳失效的功能函数可表示为：式中：x为一q维向量，代表齿轮结构的几何尺寸、材料特征和载荷等不确定性因素；动载系数Kv和载荷分配系数KH为不确定性变量，则x＝[Kv,KH]；h(x)表示回归模型的基函数，代表了齿轮疲劳失效模型的趋势；β为趋势系数向量；对于通常情况下，h(x)可取为单位向量，β可简单取为当前训练样本点集处齿轮失效功能函数的均值向量。3.如权利要求2所述的一种...

【专利技术属性】
技术研发人员：龙湘云，姜潮，田亚鑫，李博川，刘志成，
申请(专利权)人：湖南大学，
类型：发明
国别省市：湖南,43