课程名称：《随机过程》课程设计（论文）题目: 非齐次泊松过程在数控机床可靠性建模中的应用学院：理学院专业：数学与应用数学班级：数学12-1班***名：***学生学号： ********** ***师：***2015 年 1月 3 日随机过程课程设计目录任务书 (1)摘要 (1)前言 (2)１非齐次泊松过程理论 (2)1.1 非齐次泊松过程的基本理论简介 (2)1.2 基于试验总时间法的趋势检验 (2)2 数控机床的非齐次泊松过程可靠性建模 (3)2.1强度函数的建立................................................. . (3)2.2 Ｋ台数控机床强度函数的参数估计......................... (4)2.3 非齐次泊松过程下的可靠性指标............................... ....... (5)3实例分析 (5)4结束语 (7)5程序及结果 (8)6参考文献 (9)附录………………………………………………………………………………评阅书……………………………………………………………………………摘要基于试验总时间法对多样本随机截尾的数控机床现场数据进行趋势检验，在故障过程为浴盆曲线的趋势条件下构建了数控机床的非齐次泊松过程的可靠性模型。
本文使用极大似然估计法对非齐次泊松过程的强度函数进行参数估计得到了该模型的可靠性指标，以6台加工中心的现场数据为例建立了非齐次泊松过程的可靠性模型。
再通过matlab曲线拟合，绘制出故障时间的曲线，通过曲线的拟合程度，可以确定非齐次泊松过程能够更恰当地表现故障的趋势。
关键词：数控机床可靠性非齐次泊松过程浴盆曲线前言数控机床是由数目众多的零部件组成的复杂机电液可修系统。
在其可靠性研究中，需要考虑维修活动对其可靠度的影响【1】。
以往的数控机床可靠性建模方法，是将故障间隔时间视为独立同分布来分析其寿命分布，即假设维修活动是“修复如新”【2】而在实际生产中数控机床的维修活动是以调整或者更换一部分零部件和元器件为主的，对于复杂的系统来说这种维修活动只能使产品恢复到正常功能维修前后可靠度并没有很大改变，因此将数控机床的维修活动视为“修复如旧”更加合理。
非齐次泊松过程经常被用于建立“修复如旧”的维修策略且维修时间可忽略的可修系统可靠性模型用于模拟出现故障间隔时间的趋势【3-4】。
结合数控机床的维修特点，使用非齐次泊松过程建立的可靠性模型更能贴近于复杂系统的生产实际。
本文提出了非齐次泊松过程的数控机床可靠性建模方法，并结合数控机床的失效特点，建立故障率为浴盆曲线的非齐次泊松过程可靠性模型。
同时，结合具有随机截尾特点的多样本数控机床现场试验故障数据，对数控机床的可靠性进行了深入分析。
１．非齐次泊松过程理论1.1 非齐次泊松过程的基本理论简介：非齐次泊松过程是随机点过程的一种典型类型，当可修系统的相邻故障间隔呈现某种趋势时可以使用这种方法来描述。
非齐次泊松过程的重要参数【4】：)(t ω为强度函数，是非负函数；其 累 积 故 障 强 度 函 数⎰=tdu u t W 0)()(ω，表示在[0,t]中的平均故障数,即E[N(t)]=W（t ），N[t] 代表在[0,t]出现的故障次数ｔ表示机床从观测开始后的运行时间。
当强度函数为1)()(-=βααβωtt 时，成为威布尔过程。
其中，α、β>0，α为尺度参数，β为形状参数。
0<β<1, 表示不断改良的(好)系统; β>1, 表示不断恶化的(坏)系统；β=1，表示系统服从指数分布。
1.２基于试验总时间法的趋势检验：本文采用基于试验总时间的方法，对具有多样本随机截尾现场数据的故障过程进行趋势检验【5】。
将在观测期间采集到的所有故障数据按照从大到小时间进行排序，得到t (i)的时间序列。
根据试验总时间的建模思想【6】，得到该序列的第ｉ个故障发生时的试验总时间：⎰=)(0)()()(i t i du u n t T （1）=N ∑=ki i n 1（2）式中：n(u)表示在ｕ时刻观察到的数控机床数量，N 表示在观测期间的故障数；i n 表示第ｉ台机床的故障数（共有Ｋ台机床），当时间序列)(i t 的最后一个时间是故障数据时，)(i n =i n -1；当时间序列)(i t 的最后一个时间不是故障数据而是截尾数据时)(i n =i n 。
在实际检验时同时使用Ｕ－检验、Ｊ－检验和Ｖ－检验等检验方法综合确定有无趋势【7-8】。
其中Ｖ－检验如下：Ｈ0：齐次泊松过程； Ｈ1：具有非单调趋势；()11()()|()|24~(0,1)()48Ni i T S T S T t N V N N T S =--⨯=∑ (3)22()122()()|()|212()180Ni i T S T S T t N V NT S =--⨯=∑ ~(0,1)N (4)()1322log(()|2()()|)(2)Ni i T S T t T S V N χ=-=∑~(0,1)N (5)式中：Ｔ（Ｓ）表示总的观测时间。
当｜Ｖ１｜＜ｚα／２，｜Ｖ２｜＜ｚα／２，Ｖ３＜χ２（２ｎ）时，接受Ｈ０。
一般情况下，Ｕ检验和Ｊ检验是检验具有单调趋势或齐次泊松过程和更新过程的故障数据的，而当故障数据具有非单调趋势时，则可以考虑Ｖ检验中的统计量，如表１所示。
表1 故障率和故障强度函数变化特性２ 数控机床的非齐次泊松过程可靠性建模 2.1强度函数的建立对于浴盆曲线趋势的故障过程，假设其故障强度函数由早期故障期和偶然故障期两部分组成，并且每一个阶段都是一个威布尔过程［９］，参数为尺度参数αｍ和形状参数βｍ（ｍ＝１，２）。
结合以上假设和多重威布尔分布模型的性质，则该数控机床故障强度函数为：11111)()(-=βααβωtt （6）式中：α１，α２，β１，β２＞0在（0，ｔ］内的平均故障个数为累积故障强度函数，即由于该模型是由两重威布尔过程构成，其强度函数是具有非单调的浴盆曲线趋势，因此组成该模型的两个形状参数有（β１－１）（β２－１）＜0，则本文中假设β１＜１，β２＞１。
2.2 Ｋ台数控机床强度函数的参数估计本文使用极大似然估计法对k 台样本的强度函数进行参数估计，k 台数控机床的故障数据是随机截尾的，第ｉ台的故障观测时间为［0，T ｉ］,其中T i 为现场试验的截尾时间。
t 0≡0，因此,得到相应的似然函数为1212111211112212{[()()]exp[()()]}j n Kij ij i i K i j t t T T L ββββββαααααα--===+⨯--∏∏ (8)似然函数的对数函数以及此对数函数对模型参数的偏导数为1212111211112212{ln[()()][()()]}j n Kij ij i i K i j t t T T l ββββββαααααα--===+-+∑∑ (9)12211111121122()()1()()()()()jmmmij m n K m i Kij ij i j m m jt T l t t βββββααβββααααααα--==--∂=+∂+∑∑(10) 由（10）可以得到1211111121122[1ln()]()()ln()()()()()m jm ijijm m mn Kmi i Kij ij i j m m mt t T T l t t βββββαααββααααααα---==+∂=-∂+∑∑ (11)12111111121122()()()()m jmn KKm ijKij iji i j t Tt t βββββββαααα---====+∑∑∑(12)由累积故障函数可得1211112()[()(())]K K Ki ii i i i i T T W T n ββαα====+=∑∑∑ (13)以上式可导出：12111[()]Kii ii T T n βαα==-∑(14)将式（9）转换为三参数的函数，即121112111122{ln[()()()()]j n Kij ij Ki i j t tl n ββββαααα--=='=+-∑∑ (15) 最终，似然函数的参数估计转化成以下的求最大化问题：max K l ' 约束条件：一般情况，最大化问题都需要初始值。
根据经验，在没有合适的初始值选择下，可以假设2.3 非齐次泊松过程下的可靠性指标（１）首次故障间隔时间的可靠度函数【10-11】从t ＝0开始直到第一个故障发生的时间T １,T 1的可靠度函数为0()()1()Pr()tu duW t i R t T t ee ω--⎰=>== (16) 对于非齐次泊松过程模型的使用，如果能够估算出首次故障间隔时间的故障率函数，就能同时估计出产品整个寿命的强度函数。
（２）其他故障间隔时间的可靠度函数在t 0时刻后的可靠度函数00()000(|)Pr()t tt u duR t t t t t eω+-⎰=+>= (17)（３）平均故障间隔时间瞬时平均故障间隔时间———故障强度函数ω（t ）表示单位时间发生故障的次数，则ω(t)的倒数表示一次故障所经过的时间，定义瞬时故障平均间隔时间为1()IMTBF T t ω=(18) 累积平均故障时间间隔———表示一段时间内的平均故障间隔时间，即累积故障强度函数的倒数，21212112(,)()CMTBF t t t t t tT W t t u duω--==⎰(19)３实例分析以国内北京第一机床厂同一时期出厂的６台加工中心现场试验的故障数据为例，其发生故障的时间如表2所示。
首先，需要对这些数据进行趋势检验。
根据2.2节中多样本的趋势检验方法，在显著性水平α=0.05，得到这批加工中心的统计量值如表3所示。
编号故障时间/h1 2 3 4 5 650.99 423.72 753.06 760.65 795.63 1005.80 1209.40 2509.2 3350.16 3801.90 3915.62 4011.10 5109.03 5197.12 5353.92 5845.90 5942.81 6106.49 6323.63 6474.60 6526.03 6827.10 7059.69 7460.86 8240.27 8745.00 9142.65185.13 458.00 960.54 1005.87 3409.55 422.89 6061.44 6217.53 7479.45 7542.81 7775.96 7882.88 7994.25 8588.25 28.05 350.48 47.52 1560.23 1896.30 2541.10 3352.80 3915.12 4981.45 5112.97 5729.13 5812.46 5903.40 6109.13 6117.21 6275.28 6308.78 6348.21 6457.61 6620.46 6853.44 7005.85 7116.59 7249.74 7467.90 8088.96 8298.18 9509.28 131.09 785.61 287.51 870.56 2987.45 3500.75 4881.86 5136.51 5230.01 5376.53 5540.54 5746.57 6183.21 6505.13 6592.08 7125.03 7379.46 7703.03 7868.85 8275.74 8654.42 9032.10148.17 578.80 1014.14 1952.18 2893.01 3287.36 3747.55 4279.01 4714.12 4839.79 5558.09 5600.66 6694.61 6855.49 7120.48 7368.47 7496.84 7659.20 8451.26 8638.80551.98 359.4 956.72 1357.45 1549.56 2706.15 3417.46 4659.60 5150.64 5206.74 5483.61 5570.40 5651.25表２ 加工中心故障数据在V 检验中运行结果拒绝H 0可知故障数据具有非单调趋势，且由表１可知，故障发生过程呈浴盆曲线的趋势。