利用MATLAB计算AP1000可靠性数据初步探索.docVIP

利用MATLAB计算AP1000可靠性数据初步探索 　　摘 要：可靠性数据分析为核电厂设备的可靠性试验和可靠性维修等可靠性管理工作提供决策依据，在可靠性工程中占有重要地位。该文以AP1000的DTS泵为例，初步探索了AP1000电站设备可靠性数据的采集与处理、可靠性参数的分析计算这一过程。对于可靠性数据的计算，本文在行业通用数据及电站专有设备运行数据的基础上，利用Matlab平台实现了经典概率统计和贝叶斯估计两种方法的算法，摸索出一套可行的可靠性数据分析体系，引领AP1000电站设备可靠性数据分析工作的开展。 　　关键词：失效率 可靠性数据 MATLAB 贝叶斯估计 AP1000 　　中图分类号：TL353.13 文献标识码：A 文章编号：1674-098X（2015）03（a）-0033-04 　　可靠性数据是系统可靠性设计、研究、分析、评定和改进的基本依据，数据分析则是可靠性工作的基础。通过有计划、有目的地收集和分析系统各组成单元的可靠性数据，可定量评定系统的可靠性水平，发现其可靠性的薄弱环节和产生原因，经过改进，使系统的性能与可靠性水平不断提高。 　　核电厂使用的设备可靠性参数一般为以下四种：设备的运行失效率λ，设备的需求失效概率γ，设备的平均故障维修时间τ，设备由于维修及试验形成的设备总不可用度P。因篇幅和数据所限，该文的计算集中在前两个数据。 　　可靠性数据的计算和分析，涉及较多概率统计和高等数学计算，较为复杂和繁重，难以人工计算，可靠性数据分析普遍借助专用定制软件完成。MATLAB的统计工具箱是建立在MATLAB数值计算环境上的工具集，具有保证可靠性数据分析工作及时、准确、直观简便完成的能力。该文以全球首堆的AP1000工程的DTS的清水泵为例，利用MATLAB这一在工程上得到普遍应用的数学软件，实现了可靠性数据的核心计算过程。 　　1 设备基础信息 　　1.1 设备类及设备选取 　　设备可靠性数据的计算均基于设备类进行。一个设备类表示具有相似的工艺性能、相似的功能和相似的运行条件的一组设备。多数情况下，一个设备类相当于一组明确的设备组合，如该文的DTS清水泵类。下表列出了该设备类的基本信息（表1）。 　　1.2 设备边界 　　设备边界是可靠性数据库的一个重要属性。定义设备边界主要是为了提高设备失效历史数据的采集质量，同时，设备边界的划分也会影响到设备的失效模式。 　　该文所涉的DTS清水泵的边界为泵的进出口接管嘴，以法兰或焊缝为分界点。主要装置包括泵本体、润滑装置、冷却装置、以及相应的传动装置，不包括电动机。主要零部件包括泵壳、泵轴、泵轴承、叶片、冷却设备、润滑设备、联轴器。 　　2 设备失效模式 　　失效模式是设备失效的表现形式。失效模式可采用多种分类方法。为规范化和标准化，该文将设备的失效模式分为两大类：其一为运行失效，即设备在运行期间或备用期间发生的失效；其二为需求失效，即设备在状态变换时刻发生的失效。表2为基本失效模式定义描述详表。 　　2.1 设备运行和失效数据 　　2.1.1 数据源 　　表3列出了三门核电设备可靠性数据的来源和用途。 　　2.1.2 清水泵运行数据分析 　　表4对DTS清水泵的运行数据进行了详细的收集和整理。 　　综上，清水泵累计投用小时数为19680 h，净运行小时数为1010.35 h，累计启动次数为1354次。 　　2.1.3 清水泵失效数据分析 　　清水泵完全失效的标准包括：要求时不启动、非计划停泵、严重外漏、性能丧失、冷却丧失、卡住、卡轴、丧失润滑。 　　清水泵性能降级的标准包括振动、渗漏油、甩水、异音、油位低、辅助管件脱落、锈蚀等。 　　表5为清水泵失效记录 　　综上，清水泵累计运转失效次数为1次，平均维修时间为3.3 h，累计启动失效次数为0次。 　　3 设备可靠性参数的经典估计及Matlab实现 　　3.1 大样本空间下的经典估计 　　经典估计对大样本空间下的可靠性参数估计效果很好。当设备失效样本充分时，Matlab的经典估计按如下步骤进行。 　　绘制可靠性数据分布的频数直方图。失效概率密度的形状与频数直方图一致，因此可以直观假设可靠性数据服从的分布。Matlab工具箱提供hist（data，k）函数用于作频数直方图。 　　参数估计。根据直方图的形状，假设可靠性数据服从负指数分布。指数分布可用[muhat，muci]=expfit（data）进行参数估计，该命令用极大似然法完成指数分布的点估计和区间估计。 　　假设检验，观测样本数据的分布与选定的理论分布之间符合程度。通用的拟合优度检验有2检验和柯尔莫哥洛夫检验。后者可用[H，P，KSSTAT，CV]= kstest（X，cdf，alpha）命令完成。