(19)国家知识产权局 (12)发明 专利 (10)授权公告 号 (45)授权公告日 (21)申请 号 202110330699.8 (22)申请日 2021.03.26 (65)同一申请的已公布的文献号 申请公布号 CN 112949209 A (43)申请公布日 2021.06.11 (73)专利权人 北京航空航天大 学 地址 100191 北京市海淀区学院路37号 专利权人 中国兵器 工业第五九研究所 (72)发明人 马小兵　闫冰心　王晗　周堃　 吴建国　 (74)专利代理 机构 北京慧泉知识产权代理有限 公司 11232 专利代理师 王顺荣　唐爱华 (51)Int.Cl. G06F 30/27(2020.01) G06N 7/00(2006.01) G06F 111/08(2020.01) G06F 119/04(2020.01) G06F 119/08(2020.01) G06F 119/14(2020.01)(56)对比文件 CN 111859658 A,2020.10.3 0 CN 111832151 A,2020.10.27 CN 10262 2473 A,2012.08.01 US 20191389 26 A1,2019.0 5.09 Daojiang He 等.Objective Bayesian analysis for the ac celerated. 《STATISTICAL THEORY AND RELATED FIELD S》 .2018, HAN WANG 等.Mec hanism Equivalence i n Designing Optimum. 《IE EE》 .2018, Xiaobing Ma.A N ovel Parameter-Related Wiener Proces s Model W ith Its Estimati on. 《IEEE》 .2020, HONGYU WANG 等.Remai ning Useful L ife Predicti on Using. 《IEEE》 .2018, (续) 审查员 吴瑶 (54)发明名称 一种退化速率-波动联合更新的弹用密封橡 胶贮存寿命评估方法 (57)摘要 本发明提供一种退化速率 ‑波动联合更新的 弹用密封橡胶贮存寿命评估方法， 其步骤如下： 步骤一： 建立随机退化过程模型： 步骤二： 确定模 型参数的先验分布； 步骤三： 确定模型参数的后 验分布； 步骤四： 贮存寿命融合评估方法对比分 析； 通过以上步骤， 本发 明针对退化型高可靠、 长 寿命产品自然与加速退化试验数据， 根据随机过 程理论与加速模 型建立产品加速退化模型， 运用 极大似然估计方法对模型参数进行估计， 进而将 加速数据模 型参数等效到自然模 型参数， 进而确 定模型参数先验分布， 运用贝叶斯理论基于自然 退化数据更新模 型参数的后验分布， 最终保证了贮存寿命估计的准确性。 [转续页] 权利要求书6页 说明书11页 附图3页 CN 112949209 B 2022.05.17 CN 112949209 B (56)对比文件 周源等.基 于Copula函数耦合 性建模的二元 加速退化数据统计分析方法. 《兵器装备工程学 报》 .2018,(第0 5期),魏高乐等.基于多阶段-随机维纳退化过程 的产品剩余寿 命预测方法. 《科 学技术与工程》 .2015,(第26期),2/2 页 2[接上页] CN 112949209 B1.一种退化速率 ‑波动联合更新的弹用密封 橡胶贮存寿命评估方法， 其特 征在于： 设置1： 产品先验信息为加速退化数据， 利用维纳过程描述其退化， 且模型参数具有非 共轭先验分布； 设置2： 加速退化试验的应力为温度， 加载方式为恒定应力， 加速模型为阿伦尼斯模型， 阿伦尼斯加速模型 是指描述产品寿命和温度应力之间关系的物理模型； 具体步骤如下： 步骤一： 建立随机退化过程模型： 产品在时刻t的性能退化表示 为X(t)， 若随机过程X(t)为维纳过程， 则 X(t)＝ μt+σ B(t)                                (1) 其中， 漂移参数 μ表示退化速率， 波动参数σ 表示退化波动， B(t)表示标准布朗运动； 则产品退化 性能首次达 到失效阈值的时间为 其寿命L， 即 L＝inf{t|X(t)≥ω}                               (2) 其中， ω表示给定 的失效阈值， L表示产品寿命， X(t)表示产品退化过程， inf{ ·}表示 下确界； 产品寿命的概 率密度函数f(l)如下： 其中， l表示产品寿命， ω表示给定的失效阈值， μ表示退化速率， σ 表示退化波动； 产品 寿命的累计分布函数 F(l)如下： 其中， Φ(·)表示标准正态分布的分布函数， l表示产品寿命， ω表示给定的失效阈值， μ表示退化速率， σ 表示退化波动； 产品寿命的期望如下： 其中， 表示期望算子， 表示产品寿命的期望， ω表示给定的失效阈值， μ表示退化 速率； 步骤二： 确定模型参数的先验分布； 2.1加速退化试验数据参数估计； 针对恒定加速应力下退化试验数据， T0表示正常应力水平， Ti表示加速应力水平， 其中i ＝1,2,…,n， xi,j,k表示应力Ti下第j个样本在时刻k的退化观测值， 其中j＝1,2, …,m， k＝1, 2,…,l； 根据维纳过程的性质， 样本退化增量 服从正态分布 得到第j个样本的似然函数为： 权　利　要　求　书 1/6 页 2 CN 112949209 B 3