基于汉密尔顿蒙特卡罗算法的数控刀架贝叶斯可靠性评估.pdfVIP

摘要

基于汉密尔顿蒙特卡罗算法的数控刀架贝叶斯可靠性评估

数控刀架是数控机床的核心功能部件，其可靠性水平与整机的可用性密切相

关。因此研究数控刀架的可靠性对于提升国产数控机床的可靠性水平，减少与国

外先进水平的差距具有重要意义。针对数控刀架在实际试验过程中所产生的故障

数据，需要对其进行可靠性建模评估。然而在新时代和机床专项的背景下，面对

高可靠性、小批量生产的数控刀架，在进行可靠性试验的过程中，呈现了故障数

据小样本的局面。处理大样本数据的经典统计方法不再适用解决此类问题，如何

解决小样本数据的评估问题变得愈发重要，对展开国产数控刀架的可靠性全面研

究具有重要学术意义以及工程应用价值。

本文以某国产型号为SLT**的数控刀架为研究对象，基于贝叶斯框架研究了

如何解决数控刀架小样本可靠性评估的问题，具体的研究工作如下：

1.确定数控刀架小样本可靠性评估模型。将两参数威布尔分布作为数控刀

架的小样本评估模型，基于贝叶斯的理论框架推导模型参数的后验分布，并指出

解决小样本评估的两大关键点：确定先验分布和计算后验分布，将平均故障间隔

时间MTBF作为可靠性评估指标，并结合参数模型进行了理论推导。

2.建立了模型参数的先验分布。选择与研究对象相近且具有大样本数据的

数控刀架作为参考系统，通过专家打分法定性给出可靠性水平对比；通过建立专

家等级排序评价指标体系，基于层次分析法计算出专家判断的权重；基于参考刀

架在选定分位数处的时间函数值和可靠性水平的高低对比，最终通过专家打分的

方法，建立了提取专家判断信息，并通过数学转换的技巧将专家判断结果转化为

参数先验分布的流程。

3.解决后验分布参数无解析解的问题，从两个角度出发：M-H算法和传统

贝叶斯软件WinBUGS。具体开发了MCMC算法族中的基本算法:Metropolis-Ha

stings算法，结合研究对象的样本模型研究了算法中的建议分布、接受概率、迭

代步骤等并在MATLAB中编程实现，运行计算了可靠性指标结果；基于传统的

贝叶斯分析软件WinBUGS，研究建立了威布尔参数的BUGS模型，细致给出了

I

参数的仿真流程，从六个角度出发，全面对比不同迭代次数下参数的收敛性，得

到处于收敛状态下的参数仿真估计值，并最终完成可靠性指标计算。

4.针对传统MCMC算法迭代效率低、收敛速度慢和WinBUGS软件停止更

新的问题，采用了贝叶斯分析新工具Stan执行HMC算法更高效率地解决参数无

解析解的问题。结合小样本可靠性模型研究建立参数的Stan概率模型，在Rstudio

平台下编写了执行HMC算法的抽样程序，从六个角度判断参数的收敛，以更少

的迭代次数得到参数的仿真估计值，并完成可靠性指标的计算。

5.基于马氏链收敛速度的角度，对比了传统MCMC、Slice和HMC算法的

结果；基于通用性的角度，将Stan的运算结果与网格近似法进行了对比；Stan

作为新型贝叶斯分析工具，具有迭代效率高和通用性广的优点，对于其他复杂的

贝叶斯统计模型，具有很好的通用延展性。

关键词：

数控刀架，贝叶斯理论，可靠性评估，汉密尔顿蒙特卡罗，Rstan

II

Abstract

BayesianReliabilityEvaluationofNCturretbasedonHamilton

MonteCarloalgorithm

AsthekeyfunctionalcomponentofCNCmachinetool,thereliabilitylevelofCNCturretis

closelyrelatedtotheavailabilityofthewholemachine.Therefore,itisofgreatsignificanceto

studythereliabilityofNCturretforimprovingthereliabilitylevelofdomesticCNCmachinetools

andreduci