一种考虑随机载荷和强度退化的疲劳可靠性分析方法技术

本发明专利技术公开一种考虑随机载荷和强度退化的疲劳可靠性分析方法，通过建立产品的安全余量方程，并对一定次数随机载荷进行仿真，计算安全余量最小值和对应的载荷作用次数，利用核函数估计建立安全余量向量和对应载荷作用次数向量的联合概率密度函数，最终建立考虑随机载荷和强度退化的产品疲劳可靠性分析模型，相比于现有疲劳可靠性分析方法中的剩余强度模型，本申请不需要考虑剩余强度和随机载荷概率密度函数类型以及二者之间的关系，直接以安全余量方程为指标进行求解，本申请对全面提高产品的设计水平和制定合适的维修策略具有较高的理论支撑和工程实践意义。

A fatigue reliability analysis method considering stochastic loading and strength degradation

The invention discloses an analysis method considering the fatigue reliability degradation of random load and strength, the safety margin equation of the product, and for a certain number of simulation times of random load, load calculation and the corresponding minimum safety margin, estimation of the joint probability density function of the establishment of safety margin vector and the corresponding frequency of loading vectors using nuclear function, finally establish the analysis model of fatigue reliability degradation products random loading and strength, compared to the existing reliability analysis of fatigue residual strength model method, the application does not need to consider the relationship between the residual strength and the random load probability density function type and two, directly to the safety margin equation as the index for this application. To improve the design level of the products and high appropriate maintenance strategy theory Supporting and practical significance of engineering.

【技术实现步骤摘要】
一种考虑随机载荷和强度退化的疲劳可靠性分析方法
本专利技术属于可靠性分析
，尤其针对机电产品疲劳可靠性分析方法。
技术介绍
在载荷的反复作用下，各种结构及其加工而成的产品都会产生疲劳。据统计，疲劳破坏在机械结构破坏中所占比例高达50％-90％，疲劳破坏往往十分突然，并且造成灾难性的事故，引起巨大的经济损失。特别是近40年来，在机械向高温、高速以及大型化方向发展，机械零件所承受的应力越来越高，服役工况越来越复杂，导致疲劳失效事故越来越多。因此，疲劳可靠性分析方法的研究，对预防疲劳破坏事故的发生和经济发展意义重大。疲劳可靠性的研究目的主要在于从经济性和维修性出发，在规定的工况下，规定的使用时间内，以及完成规定功能的条件下，使得产品因疲劳而失效的可能性降到最低程度。目前结构元件的疲劳分析方法主要有三类：一是疲劳累积损伤模型，该模型基于传统的机构疲劳寿命模型，结合疲劳累积损伤理论和应力统计分布，得到疲劳寿命的相关统计参数，同时该模型比较简单，只能反映疲劳破坏的基本趋势；二是疲劳寿命模型，该模型以疲劳寿命作为参数，结合疲劳累积损伤进行疲劳可靠性分析，但是在随机载荷作用下疲劳寿命的瞬时分布的解析表达式很难给出；三是剩余强度模型，该模型通过对元件剩余强度在疲劳载荷下随时间变化规律的研究，从而给出其寿命参数。目前的研究工作主要关注在以剩余强度为指标的疲劳可靠性分析方法，而直接以安全余量为指标的疲劳可靠性研究尚未开展。
技术实现思路
本专利技术为解决上述技术问题，提出了一种考虑随机载荷和强度退化的疲劳可靠性分析方法，以安全余量为指标，通过随机载荷仿真，构建安全余量向量和对应的载荷作用次数向量，从而建立疲劳可靠性分析模型。本专利技术采用的技术方案是：一种考虑随机载荷和强度退化的疲劳可靠性分析方法，包括：S1、根据建立的产品随机载荷数学模型和剩余强度数学模型，构建产品安全余量数学模型；S2、给定随机载荷作用次数，通过仿真得到安全余量向量和相对应的载荷作用次数向量；S3、结合步骤S2所得的安全余量向量和相对应的载荷作用次数向量，利用核函数估计建立安全余量向量和相对应的载荷作用次数向量的联合概率密度函数，最终建立考虑随机载荷和强度退化的产品疲劳可靠性分析模型。进一步地，所述步骤S1包括：S11、通过分析产品在实际工况中的随机载荷作用情况，建立随机载荷数学模型；S12、根据步骤S11得到的随机载荷数学模型，建立剩余强度模型；S13、根据步骤S11得到的随机载荷数学模型，以及步骤S12得到的剩余强度模型，建立产品安全余量模型。进一步地，所述步骤S2包括：S21、给定随机载荷作用次数，通过仿真得到安全余量最小值和相对应的载荷作用次数；S22、若步骤S21所得的安全余量最小值小于零，则用安全余量小于零的第一时间点对应的安全余量值和载荷作用次数作为最小值点对应的安全余量最小值和载荷作用次数；否则保持原来的安全余量最小值和对应的载荷作用次数不变；S23、重复N次步骤S21、步骤S22的操作，得到N维的安全余量向量，和N维的载荷作用次数向量。更进一步地，所述步骤S21具体包括以下分步骤：S211、设定随机载荷作用次数为L1，则仿真一次得到包含L1个元素的随机载荷向量；S212、根据步骤S211得到的随机载荷向量得到随机载荷谱，并根据该随机载荷谱计算得到剩余强度向量；S213、结合步骤S211得到的随机载荷向量和S212得到的剩余强度向量，根据安全余量数学模型的计算公式，得到安全余量向量；所述的安全余量数学模型的计算公式为：M(ω)＝R(ω)-s(ω)；其中，M(ω)表示安全余量，R(ω)表示剩余强度，s(ω)表示随机载荷，ω表示载荷作用次数；S214、对步骤S213所得的安全余量向量进行比较，选择安全余量的最小值和安全余量最小值对应的载荷作用次数。更进一步地，步骤S3所述核函数K(v)的表达式为：其中：v表示核函数参数。进一步地，所诉步骤S3还包括：A1、计算安全余量向量的均方差，得到核密度估计中的安全余量向量最佳带宽；A2、计算载荷作用次数向量的均方差，得到核密度估计中的载荷作用次数向量最佳带宽；A3、结合步骤A1得到的安全余量向量最佳带宽，步骤A2得到的载荷作用次数向量最佳带宽，以及核函数公式得到安全余量向量和载荷作用次数的联合概率密度函数；A4、根据步骤A3得到的联合概率密度函数，建立考虑随机载荷和强度退化的产品疲劳可靠性分析模型。更进一步地，所述步骤A4具体为：根据安全余量向量和载荷作用次数向量的联合概率密度函数，得到疲劳失效概率；根据疲劳失效概率，得到考虑随机载荷和强度退化的产品疲劳可靠性分析模型。更进一步地，所述疲劳失效概率表达式为：其中，w表示随机载荷作用次数，x为取最小值点时对应的载荷作用次数；y为最小值点对应的安全余量数值。本专利技术的有益效果：本专利技术通过建立产品的安全余量方程，并对一定次数随机载荷进行仿真，计算安全余量最小值和对应的载荷作用次数，利用核函数估计分别建立安全余量向量的概率密度函数、安全余量向量和对应载荷作用次数向量的联合概率密度函数，最终建立考虑随机载荷和强度退化的产品疲劳可靠性分析模型，相比于现有的剩余强度模型，本申请不需要考虑剩余强度和随机载荷概率密度函数类型以及二者之间的关系，直接以安全余量方程为指标进行求解，本申请对全面提高产品的设计水平和制定合适的维修策略具有较高的理论支撑和工程实践意义。附图说明图1是本专利技术疲劳可靠性分析方法的流程图。图2是本专利技术实施例中滤波减速器固定齿轮的三维模型图。图3是本专利技术实施例中滤波减速器固定齿轮安全余量向量的直方图。图4是本专利技术实施例中滤波减速器固定齿轮安全余量向量对应的载荷作用次数向量的直方图。图5是本专利技术实施例中可靠度随载荷作用次数的变化图。具体实施方式为便于本领域技术人员理解本专利技术的
技术实现思路
，下面结合附图对本
技术实现思路
进一步阐释。以滤波减速器固定齿轮为例进行说明，在以本专利技术技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本专利技术的保护范围不限于下述的实施例。如图1所示为本专利技术的方案流程图，本专利技术采用的技术方案是：一种考虑随机载荷和强度退化的疲劳可靠性分析方法，针对如图2所示滤波减速器固定齿轮，包括以下步骤：S1、建立产品随机载荷和剩余强度数学模型，结合随机载荷模型和剩余强度模型得到产品的安全余量数学模型；具体包括：S11、分析根据谐波减速器固定齿轮在实际工况中和双联齿轮的啮合情况，固定齿轮所受的随机载荷s(ω)的数学模型如公式(1)所示：公式(1)中，ZH为区域系数，本实施例中取值为2.5；ZE为材料的弹性影响系数，本实施例中取值为189.8；Zε为重合度系数，本实施例中取值为0.9642；K为载荷系数，本实施例中取值为0.3635；T(ω)为输出转矩，为随机过程，在给定ω时，假定T(ω)服从正态分布N(50933，4244.12)；b为固定齿轮齿宽，服从正态分布N(5，0.0332)；d为固定齿轮分度圆直径，服从正态分布N(79.5，0.062)；μ是滤波减速器固定齿轮与双联齿轮的齿数比，本实施例中取值为0.9056。S12、结合步骤S11得到的随机载荷作用数学模型，得到固定齿轮的剩余强度R(ω)模型如公式(2)所示：公式(2)中，δs表示固定齿本文档来自技高网...

【技术保护点】
一种考虑随机载荷和强度退化的疲劳可靠性分析方法，其特征在于，包括：S1、根据建立的产品随机载荷数学模型和剩余强度数学模型，构建产品安全余量数学模型；S2、给定随机载荷作用次数，通过仿真得到安全余量向量和相对应的载荷作用次数向量；S3、结合步骤S2所得的安全余量向量和相对应的载荷作用次数向量，利用核函数估计建立安全余量向量和相对应的载荷作用次数向量的联合概率密度函数，最终建立考虑随机载荷和强度退化的产品疲劳可靠性分析模型。

【技术特征摘要】
1.一种考虑随机载荷和强度退化的疲劳可靠性分析方法，其特征在于，包括：S1、根据建立的产品随机载荷数学模型和剩余强度数学模型，构建产品安全余量数学模型；S2、给定随机载荷作用次数，通过仿真得到安全余量向量和相对应的载荷作用次数向量；S3、结合步骤S2所得的安全余量向量和相对应的载荷作用次数向量，利用核函数估计建立安全余量向量和相对应的载荷作用次数向量的联合概率密度函数，最终建立考虑随机载荷和强度退化的产品疲劳可靠性分析模型。2.根据权利要求1所述的一种考虑随机载荷和强度退化的疲劳可靠性分析方法，其特征在于，所述步骤S1包括：S11、通过分析产品在实际工况中的随机载荷作用情况，建立随机载荷数学模型；S12、根据步骤S11得到的随机载荷数学模型，建立剩余强度模型；S13、根据步骤S11得到的随机载荷数学模型，以及步骤S12得到的剩余强度模型，建立产品安全余量模型。3.根据权利要求1所述的一种考虑随机载荷和强度退化的疲劳可靠性分析方法，其特征在于，所述步骤S2包括：S21、给定随机载荷作用次数，通过仿真得到安全余量最小值和相对应的载荷作用次数；S22、若步骤S21所得的安全余量最小值小于零，则用安全余量小于零的第一时间点对应的安全余量值和载荷作用次数作为最小值点对应的安全余量最小值和载荷作用次数；否则保持原来的安全余量最小值和对应的载荷作用次数不变；S23、重复N次步骤S21至步骤S22的操作，得到N维的安全余量向量，和N维的载荷作用次数向量。4.根据权利要求3所述的一种考虑随机载荷和强度退化的疲劳可靠性分析方法，其特征在于，所述步骤S21具体包括以下分步骤：S211、设定随机载荷作用次数为L1，则仿真一次得到包含L1个元素的随机载荷向量；S212、根据步骤S211得到的随机载荷向量得到随机载荷谱，并根据该随机载荷谱计算得到剩余强度向量；S213、结合步骤S211得到的随机载荷向量和S212得到的剩余强度向量，根据安全余量数...

【专利技术属性】
技术研发人员：汪忠来，俞水，程晓文，刘静，王磊，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川,51