买专利,只认龙图腾
首页 专利交易 科技果 科技人才 科技服务 商标交易 会员权益 IP管家助手 需求市场 关于龙图腾
 /  免费注册
到顶部 到底部
清空 搜索

【发明授权】一种多元参数退化相依的细纱机可靠性评估方法_西安工程大学_202110784692.3 

申请/专利权人:西安工程大学

申请日:2021-07-12

公开(公告)日:2023-07-28

公开(公告)号:CN113361142B

主分类号:G06F30/20

分类号:G06F30/20;G06F119/02

优先权:

专利状态码:有效-授权

法律状态:2023.07.28#授权;2021.09.24#实质审查的生效;2021.09.07#公开

摘要:本发明公开了一种多元参数退化相依的细纱机可靠性评估方法,通过分析细纱机性能参数之间关系,引入主客观综合权重辨识方法对监测的性能参数进行了筛选;接着利用非线性Wiener过程对多元参数退化失效轨迹进行了拟合,同时利用Copula函数对多元参数退化失效过程进行了相关性分析,并利用Gibbs算法进行了参数估计,进而利用Weibull分布对细纱机性能退化过程进行了表征。对细纱机的可靠性实现了评估。本发明实现了对细纱机的可靠性评估,有利于提高提升细纱性能的可靠性水平。

主权项:1.一种多元参数退化相依的细纱机可靠性评估方法,其特征在于,包括以下步骤:S1、从细纱机性能退化失效的影响因素入手,对细纱机运行过程中设备参数、工艺数据、以及每个班次的细纱机运行数据进行收集,同时对细纱机各类停机故障进行统计分析,认为在成纱过程中,当实际监测的运行参数与设定的工艺参数之间数值不一致时,被认定为细纱机性能发生了退化;主要包括考虑捻度、转速、钢领板上升速度和下降速度关键参数,并将其作为细纱机性能退化的原始性能参数,记为X;设细纱机历史性能数据包括N个样本点,性能退化参数为M维,则细纱机在监测时间ti的监测值Xi由M维性能监测参数构成,而且每一维性能参数都是长度为N的序列;即在时间序列t=t1,t2,…tn监测得到的细纱机监测数据的特征序列为X=x1,x2,…xn;S2、通过引入主客观综合权重对细纱机的关键特征参数进行识别,计算参数影响细纱机性能退化的比重,从而选取出表征细纱机退化的关键参数; 其中,pj表示第j个特征参数的重要度,在上式中,引入熵权法进行综合权重计算,具体的计算公式为: 上式中,在样本数据yij中y′ij表示为第j个参数的第i个样本数据所占的比重,m为样本量,n为特征参数数量,ej为第j个特征参数的熵权,qj为利用熵权法得到的第j个特征参数的权重值;为了提高权重评价结果的准确性,利用最小二乘法,构建面向细纱机性能退化的指标权重w,具体过程详见如下所示: 上式中,wj为待求的综合权重;在此基础上,选择综合权重大于0.5的参数,作为细纱机性能退化的表征参数,并将其应用于细纱机自身性能退化的全过程;S3、性能退化模型的构建1构建基于非线性维纳过程的细纱机性能退化模型:Xt=X0+μt+σβt7上式中,X0为初始性能退化量,μ为飘移系数,βt为标准布朗运动,σ为扩散系数,设定X0=0;像细纱机这种复杂系统的退化过程往往具有非线性的特征,因此要将非线性转化为线性特征;进而,在式7的基础上,利用时间尺度转换函数τ将其进行转换,结果详见如下:τ=τt,γ=tγ8Yτ=Y0+μτ+σβτ9上式中,γ为形状参数,式9为线性Wiener过程;同时,令细纱机的寿命为T,失效阈值ω,则细纱机的性能退化轨迹由9式的Wiener过程进行描述;由此,细纱机的寿命T定义:T={t:Yt≥ω|Y0<ω}10由式10可知,细纱机寿命T达到失效阈值ω的时间分布,服从逆高斯分布,进而得到对应的概率密度函数为: 相应地,寿命分布函数为: 上式中,Φ·为标准正态分布的分布函数;由此,可靠度函数表示为: 上式中,ΔYi=Yit-Yi-1t为细纱机在时刻τi-1和τi之间的退化增量,Δτi=τi-τi-1为时刻τi-1和τi之间的时间间隔;由Wiener过程的性质可知,细纱机性能退化增量之间相互独立,而且服从正态分布ΔYi~NμΔτi,σ2Δτi;于是,通过细纱机的性能退化数据,得到如下所示的似然函数; 通过式14,对漂移系数和扩散系数求偏导,令偏导为零并进行方程求解,从而得到μ和σ2的估计值,详见如下: 将ΔYi、Δτi代入式15、式16,得未知参数的估计值,并代入式11、式13得细纱机一元性能退化过程中的寿命密度函数和可靠度;S4、多退化过程的细纱机可靠度建模假设细纱机的最大失效阈值为ωj,当第j个退化过程的退化量Yjτ超过其对应的最大失效阈值ωj时,细纱机发生退化生效;为式中,fyjτ为第j个退化过程对应的概率密度函数;在此基础上,对可靠度Rjτ进一步构建多元退化相关失效可靠性模型R'τ,具体见式17: 为便于细纱机多退化相关性问题的分析,利用Copula函数进行细纱机多元退化相关性的建模;设细纱机的性能退化量为Y1τ、Y2τ,与其对应的性能退化函数为则的联合分布函数表示为:HY1τ,Y2τ;θ=CFTY1τ,FTY2τ;θ18 的联合密度函数为: 在常见的Copula函数形式中,选择合适的函数进行相关性分析;由于细纱机性能退化复杂,其在相关性建模后变得更加复杂,而且未知参数较多;为此,利用基于贝叶斯理论的Gibbs算法进行参数估计;由此,将两个过程所有的未知参数设为ξ,则ξ为一个n维变量,即ξ=ξ1,ξ2,…ξn,对应的先验分布为pξ1,ξ2,…ξn,具体步骤如下:S41,给定一组符合马尔科夫链的初始值ξ0=ξ10,ξ20,…ξn0;S42,从条件概率密度pξn|ξ10,ξ20,…ξn-10抽取ξ11;S43,重复步骤S42,直到从pξn|ξ10,ξ20,…ξn-10抽取ξn1,完成一次迭代;S44,重复步骤S42和步骤S43,迭代m次,得到样本ξm=ξ1m,ξ2m,…ξnm;当m足够大时,ξm看作ξ的真值,从而求得样本ξ的后验分布qξ1,ξ2,…ξn,并实现未知参数值的估计;若给定失效阈值ω,通过上述算法进行参数估计,得到未知参数并将其代入式17,得到多元参数的细纱机可靠度;S5、突发退化失效模型的构建将细纱机t时刻退化量为x对应的突发失效故障概率为λct|x,则根据威布尔分布建立如下关系式: 通过式12,式13、式14为得到的突发失效时间Tc的概率密度函数和累积分布函数; 由此,细纱机在t时刻性能突发失效的可靠度Rct|x表示成如下关系式: 记录n台细纱机在性能突发失效时间Tc1,Tc2…,TcN对应的性能退化量为x1,x2,…xN,则根据式14得到似然函数为: 通过式24,求解得到突发失效参数和并将其代入式23,得到突发退化失效的可靠度;S6、利用主客观方法筛选出能够表征细纱机性能退化的参数,利用Copula函数考虑多元参数之间的相关性,进而对基于非线性维纳过程的性能退化模型中的未知参数进行估计,将估计得到的参数值带入推导得到的可靠性函数和剩余寿命密度函数,即实现细纱机性能退化过程中可靠性的评估;对性能退化发生突变的细纱机进行跟踪检测,在获取监测数据的基础上,将数据代入基于威布尔分布的突发退化失效模型,计算得到突发失效的可靠度,进而判断细纱机突发失效和退化失效之间独立或者相依两种情况下;计算得到细纱机性能突发退化失效的可靠度,对细纱机的可靠性实现了评估。

全文数据:

权利要求:

百度查询: 西安工程大学 一种多元参数退化相依的细纱机可靠性评估方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息,力求客观、公正,但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解,仅供参考使用,不能作为本公司承担任何法律责任的依据或者凭证。