兰利法方差与可靠度估计影响因素研究.pdfVIP

兰利法方差与可靠度估计影响因素研究 花成1张盛国2 1．中国工程物理研究院化工材料研究所，四川绵阳1919信箱311分箱621900 2．二炮驻中国工程物理研究院型号办，四川i绵阳919信箱621900 摘要本文通过Monte-Carlo模拟的方法研究了主要用于估计感度分布参数的兰利法试验，对方差 参数估计较真值系统偏小的问题进行了若干定量研究，同时研充了单组试验发数、门宽等圉素对 方差参数及可靠性估计的影响． 关键词兰刺法模拟参数估计可靠,g／q-估 1引言 对于炸药、火工品等测试样品，其临界刺激量t无法直接进行测试，只能对其施加一外界刺 激算，当z≥‘时，样品出现响应；z‘时，样品不出现响应。通过各个样品的试验结果，只能 知晓其I临界刺激量是否大于或小于外界刺激量一试验刺激量。由对各样品施加的试验刺激量及大 于或小于关系，构成了特有的试验数据一敏艨性试验数据。 敏感性试验主要就是临界刺激量五的统计问题，每一个样品都具有各自特有的临界刺激量， 所有样品的临界刺激量组成的分布，对炸药及火工品称为感度分布。感度试验就是通过对临界刺 激量的试验研究，利用统计方法推导感度分布参数。分布参数的获取和分析，是炸药、火工品等 爆炸产品可靠性评估的基础。 通常，人们使用非序贯试验方法和序贯试验方法获取感度分布参数。非序贯试验方法最典型 的为爆炸百分率法和计数法；序贯试验方法有传统的Bruceton升降法“1和变步长的兰利法 (Langlie)[n、08lR法、Wu”法及Neyer系列0lmmal方法。 2兰利法试验与模拟方法 Bruceton方法为固定步长的升降法，不能使刺激量较精确地向均值逼近，直接影响估值的精 度和效率。1963年Lal捌e KJ提出了一种变步长的升降法一兰利法tLanglie)曲。它的试验程序中 规定了按当前的响应与否的试验结果，依据从试验开始至当前所获得的数据，确定下一试验所用 刺激量，所用方法与Bruceton方法有较大差异．根据最大似然函数及牛顿法可算出均值与方差， MII，STD．33lA埘给出了由数值模拟方法得到的经验结果：均值是无偏的，而方差的估值较真值系 统偏小。 兰利法要根据试验前的经验，选定参数：最大全不响应刺激量札和最小全响应刺激量五，使 试验区间能覆盏临界刺激量的?全部分布区间。使用兰利法在测试炸药的安全性能时。规定一旦出 现响应则响应情况记为D(Oown)，取降步(降低试验刺激量)继续试验；没出现响应时，响应情况即 记为u(up)．取升步继续试验。试验时，除了记录其u、D结果外，还要记录台阶数k及在刺激量 而下的响应数川。第一个刺激量试验XI=(x．+xD／2，在用试验刺激量如做完第i个单元试验后，如果 228 还需要继续试验。取刺激量 而+，=(玉十膏)，2 数至』时第一次出现U、D个数相等，则i=- 而若无相等位置出现，则 ， f毛，当第龃结果为U时 而51却当第龃结果为D日j 试验的终止准则，可以用预定试验数，或者预定玑D结果转换次数等选择．并且同样要求 数据应含有混合结果区，即试验结果显示最小试样爆炸刺激量札最大试样不爆炸刺激量劫口Ⅳ1。 对兰剥法进行Monte-Carlo模拟时，对产生的伪随机数集合进行检验与分析，使其具有较好的 随机性及均匀性和较长的周期，通过，检验法对随机数进行检验。 采用Monte．Carlo方法对感度试验兰利方法进行了模拟，先使用模拟程序产生符合要求的随机 数集合?，以代替该组试验试样临界刺激量的集合。利用线性摸余方法产生随机数，线性模余方法 其中M为模数，C为调整系数，)【0为初始数，本文使用的随机数周期为107量级。 兰利法试验的主要模拟步骤如下： (1)产生【0，1】均匀分布的随机数ri，并使随机数满足均匀性和随机性，按正态分布反函数产 生始，作为临界刺激量． (2)产生随机刺激量卸．第一试验点取初始参数瓤=‰z+xov)／2，下一个试验点的估计量公式 为而t】=(而+《)，2。 (3)比较两种变量大小，当Xc,Xi，记Y。=l，当拍≥而，记Yi=o。 (4)产生下一对鼬和而，比较变量大小并记录结果Yj，直至完成一组试验，求肛、a