汽车发动机加速寿命试验数据处理方法的计算机模拟.pdfVIP

第九届全国地面机器系统学会学术年会论文集 140 U牟占中 汽车发动机加速寿命试验 数据处理方法的计算机模拟 姜华①曹正清 中国农业大学车辆工程学院 摘要本文提出了～种汽车发动机加速寿命试验相应失效分布的参数估计公式，并用蒙特卡 洛法进行模拟验证．根据模拟结果．表明参数估计公式正确，采用蒙特卡洛法进行模拟验证可 行． 关键词汽车》动机加速寿命试验数耱处理计算机模拟 目前在我国汽车发动机加速寿命试验中，发动机发生故障较多。分析故障发生的规律对 提高汽车发动机产品的可靠性很有价值。本文根据汽车发动机加速寿命试验数据的特点提 出相应失效分布的参数估计公式，并用蒙特卡洛法进行模拟验证。 1 汽车发动机加速寿命试验数据的特点 对多台汽车发动机进行加速寿命试验、可靠性试验或在正常使用过程中，可获得如下形 式的汽车发动机试验数据。 1ll，E。…………………………………………T1。T㈣ ……………………，，J，…………………………Z。z。 1 TⅢL2．．……………………………………L．^丁㈨ 1 中的数据，后者减去前者即得故障间隔值： X¨rXl2…………………………………………X1_Xl_lo ………………………X∥……………………”X…XⅢ 2 Xn，X船…………………………………………X^。‘X^。^o 式中： X。 为第i台第j个故障间的工作时间； 丁。， 为第i台第J个累积工作时间； 也 为第i台总故障数 ^ 为统计的总台数； x。。 第i台无故障截止时间； E。 为第i台试验结束的累积工作时间，若试验以故障结束丑4有T。。一7’。 ‘。。。。。。。。。。。。-。。。_。。____，。。●●--。。__-。。。。一 ④姜华，博士生，北京市海淀区清华东路17号中国农业大学东校区213．100083 141 姜华等：汽车发动机加速毒命试验数据处理方法的计算机模拟 当数据 1 遵从更新过程，数据 2 可用威布尔分布、正态分布、对数正态分布、指数分布 等模型加以描述。下文给出这种数据的相应失效分布公式的主要推导过程和结果。 2相应失效分布的参数估计公式“1 本文推导出多个可修系统情况下分布参数的似然函数通式，不仅适用于汽车发动机加 速寿命试验相应的失效分布形式，也适用于其它分布． 为欲估计的分布参数。将数据 2 中各系统故障间隔值，按从小列大的顺序重新排列，先从第 一个系统着手分析．按顺序统计量理论，这些子样同时出现的概率为： f* ’ ’r一 得第一个系统似然函数L， 啦+1 j1上， z∽O [1--F xⅫ。，口 ] ，。】 对于^个系统，由于它们之间是相互独立的，得似然函数 ^ 厶一Ⅱ 啦+1 ! o∞ — Fo∞ i一1 。Ⅱ㈦ 。．Ⅱ州 厂 。Ⅱ埘 n 眦 ] 为计算方便，引入对数 lnLx一∑∑l。， ％，目 +∑ln[1--F ‰。，∞] 函数．所以要使厶达到最大，等价于使．1nL。达到最大值． 值，分别对以，f一1，8-求偏导 蠢 警 3 就可求得未知参数F 即吼，以… 。对于我们掌握的汽车发动机加速寿命试验数据，经检 exp[ 1---X 。]代入 3 ，得威布尔分布的参数估计公式 ，∑∑搦1nx一∑“如嚣枷∑∑ln勺 扣簟五惹一专∑∑‰+∑梳 ∑m “’ f∑∑x：+∑‰。16 _ 立L型—r—蔓一 5 【 ∑m J 姜华等：汽车发动机加速寿命试验数据处理方法的计算机模拟 142 3计算机模拟验证 以上从理论上证明了笔者提出的威布尔分布的参数估计公式，下面用蒙特卡洛法来验 证该公式是否正确。 3．1多组均匀独立的随机数的产生 在这里采用线性同余法，给出鹭信度0．99，来产生每组400个，共4组的伪随机数，并 进行均匀性、独立性检验m”，结果如下：均匀性检验：柯尔莫哥洛夫一斯米尔洛夫检验D 0．01949，D n D，0．08150 。独立性检验：取每个伪随机数闯的间隔为I～12，结果如 下： 084≤Z¨9 Z”9 2．325 #k。 2．185 3．2符合某种分布的随机变量的产生 将随机数代入 6 式 X。 F一1 肚 6 对威布尔分布将随机数代入 7 式 x。 砸一ln 1一U， ]： 7 其中：m，7为事先给定的形状、尺度参数m为随机数，X，为随机变量，F1为分布函数的反 函数． 3．3参数估计公式的检验 根据所产生的符合某分布的随机变量和事先给定的某分布的参数0，由参数估计公式 算出估计值；．若0与吾比较接近，则说明参数估计公式是正确的，用蒙特卡洛法进行模拟 是可行的。对于威布尔分布的具体模拟结果如下： 1 全子样情况；其模拟结果如表1所示。 表1空子样模拟结果 m 7 m 日 m—tH ／m‘口一口 ／ 288 一