CN114595960A - 一种基于阈值双区间的火工品安全性评估方法 - Google Patents

一种基于阈值双区间的火工品安全性评估方法 Download PDF

Info

Publication number
CN114595960A
CN114595960A CN202210214960.2A CN202210214960A CN114595960A CN 114595960 A CN114595960 A CN 114595960A CN 202210214960 A CN202210214960 A CN 202210214960A CN 114595960 A CN114595960 A CN 114595960A
Authority
CN
China
Prior art keywords
initiating explosive
explosive device
severity
interval
probability
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210214960.2A
Other languages
English (en)
Inventor
穆慧娜
刘炜
伊枭剑
温玉全
马文涛
崔宇航
李晓刚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN202210214960.2A priority Critical patent/CN114595960A/zh
Publication of CN114595960A publication Critical patent/CN114595960A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0635Risk analysis of enterprise or organisation activities
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Abstract

本发明公开了一种基于阈值双区间的火工品安全性评估方法,包括以下步骤:S1.给出激源强度和火工品感度均值的阈值区间;S2.根据区间模型计算火工品发生意外事故的概率;S3.根据火工品发生意外事故严重程度,确定火工品事故严酷度等级;S4.根据火工品发生事故的概率和严酷度等级,采用概率风险评价法对其进行安全性评估。本发明可以在评估火工品的安全性试验数据和信息较少的情况下,综合利用历史数据或同类产品的试验数据对其进行安全性评估,具有评估精度高、所需样本量少的优点。

Description

一种基于阈值双区间的火工品安全性评估方法
技术领域
本发明涉及火工品评估,特别是涉及一种基于阈值双区间的火工品安全性评估方法。
背景技术
火工品在各类武器装备中广泛应用,具有独立首发、极端敏感和一次性作用等特点,其安全性影响整个武器效能的发挥,是武器系统能够正常作用的首要前提,因此对其进行安全性评估就显得尤为重要。
现有风险评估中事故发生概率的计算均需要大量的统计数据来确定参数或变量的概率密度函数和隶属函数,建立符合实际情况且较为精确的计算模型,以便准确地计算其发生概率。然而,实际进行火工品安全性评估的时候常常无法被满足。造成这个问题的原因一方面是由于火工品价值比较高,火工品的感度试验往往是小样本,无法或不可能获得大量的统计数据;另一方面,由于感度试验中模型的假设和参数估计的影响,建立完全准确的感度模型是不太可能的,也不太现实,所以基于感度模型进行计算,必将影响事故发生概率值的准确性,从而影响评估结果的正确性。
鉴于上述原因,为了能更有效地处理因统计数据不足、感度模型本身导致的不确定性问题,提出基于非概率理论的双区间模型来获得火工品发生事故的概率进行安全性评估的方法,只需知道参数的阈值区间而不要求其具体的分布形式,且不涉及概率的概念,没有必要求其概率密度函数或者隶属函数,解决统计数据缺乏或者难以得到试验数据的问题,提高安全性评估的精度。
发明内容
本发明的目的在于克服现有火工品安全性评估方法中统计数据缺乏、需求样本量大、模型精度不高等技术的不足,提供一种基于阈值双区间的火工品安全性评估方法,在评估产品安全性试验数据和信息较少的情况下,利用少量试验数据进行安全性评估,具有评估精度高、所需样本量少的优点。
本发明的目的是通过以下技术方案来实现的:一种基于阈值双区间的火工品安全性评估方法,包括以下步骤:
S1.给出激源强度和火工品感度均值的阈值区间;
S2.根据区间模型计算火工品发生事故的概率;
S3.根据火工品发生意外事故的严重程度,确定事故严酷度等级;
S4.根据火工品发生事故的概率和严酷度等级,采用概率风险评价法对其进行安全性评估。
进一步地,所述步骤S1包括:
S101.假设激源强度和火工品感度均值不确定性参数均在某区间范围内变化,则有:
Figure BDA0003534082860000021
式中,
Figure BDA0003534082860000022
Figure BDA0003534082860000023
为不确定性参数的上界,YyXx为不确定性参数的下界。
S102.利用区间表示法,将式(1)表示成如下的区间形式:
Figure BDA0003534082860000024
S103.设激源的强度区间为yl,火工品的感度均值区间为xl,发生事故的概率函数M为:
M=yl-xl
进一步地,所述步骤S2包括:
当意外激源作用于火工品时,火工品发生事故的概率取决于意外激源的强度和火工品感度均值的大小;假设在某一时刻,只有一种激源作用于火工品;设激源强度的阈值区间为
Figure BDA0003534082860000025
火工品感度均值的阈值区间为
Figure BDA0003534082860000026
则火工品发生事故的概率P为:
Figure BDA0003534082860000027
式中,w(yl)和w(xl)分别为区间yl和xl的宽度,即
Figure BDA0003534082860000028
根据这两个区间的位置关系,式(3)等价地改写为:
Figure BDA0003534082860000029
根据式(4),第1种情况case1发生事故的概率为P=1,区间yl完全处于区间xl的右边,即区间yl总是大于区间xl;第3情况case3发生事故的概率为P=0,区间yl完全处于区间xl的左边,即区间yl总是小于区间xl;第2种情况case2发生事故概率P的取值范围为[0,1],两区间发生相互干涉的情况,当两区间的宽度不变时,其概率值由区间yl的上界
Figure BDA00035340828600000210
和区间xl的下界x的相对位置所决定。
进一步地,所述步骤S3包括:
根据下表所示严酷度等级的判别标准,在已知严酷度类别和后果严重程度的基础上,给出火工品事故严酷度等级S。
严酷度类别 后果严重程度描述 严酷度等级S
非常严重 灾难性伤亡和损失 ≥27
严重 重大伤亡和损失 [23,27)
中等严重 较大伤亡和损失 [19,23)
比较严重 一定伤亡和损失 [12,19)
轻度 较小伤亡和损失 <12
进一步地,所述步骤S4包括:
根据火工品发生事故的概率P和严酷度等级S,采用概率风险评价法对其进行安全性评估,则对应的风险值R为:
R=P×S。
本发明的有益效果是:本发明可以在评估产品的安全性试验数据和信息较少的情况下,综合利用少量试验数据进行安全性评估,具有评估精度高、所需样本量少的优点。
附图说明
图1为本发明的方法流程图;
图2为计算火工品事故发生概率的双区间关系示意图。
具体实施方式
下面结合附图进一步详细描述本发明的技术方案,但本发明的保护范围不局限于以下所述。
当我们所掌握的数据信息较少时,可采用双区间模型来描述输入不确定性因素,获得参数的不确定性变化范围,由此变化范围与允许的变化范围作比较,进而得到用于风险值计算的事故发生概率。在概率论与数理统计的研究体系,基于双区间的方法属于非概率方法的范畴。非概率方法与传统的概率方法相比,它并不是强调可接受行为的概率,而是强调可接受行为的范围。因此,利用双区间模型来对火工品进行安全性分析的优点在于对原始数据要求很低,只需知道不确定参量的界限而不要求其具体的分布形式,且不涉及概率的概念,没有必要求其概率密度函数或者隶属函数,当统计数据缺乏或者难以得到试验数据时可首选此方法。
双区间模型它具有概念清晰、模型简单、运算方便等优点,具体地:
如图1所示,一种基于阈值双区间的火工品安全性评估方法,包括以下步骤:
S1.给出激源强度和火工品感度均值的阈值区间;
S2.根据区间模型计算火工品发生事故的概率;
S3.根据火工品发生意外事故的严重程度,确定火工品事故严酷度等级;
S4.根据火工品发生事故的概率和严酷度等级,采用概率风险评价法对其进行安全性评估。
进一步地,所述步骤S1包括:
S101.假设激源强度和火工品感度均值不确定性参数均在某区间范围内变化,则有:
Figure BDA0003534082860000041
式中,
Figure BDA0003534082860000042
Figure BDA0003534082860000043
为不确定性参数的上界,YyXx为不确定性参数的下界。
S102.利用区间表示法,将式(1)表示成如下的区间形式:
Figure BDA0003534082860000044
S103.设激源的强度区间为yl,火工品的感度均值区间为xl,发生事故的概率函数M为:
M=yl-xl
进一步地,所述步骤S2包括:
当意外激源作用于火工品时,火工品发生事故的概率取决于意外激源的强度和火工品感度均值的大小;假设在某一时刻,只有一种激源作用于火工品;设激源强度的阈值区间为
Figure BDA0003534082860000045
火工品感度均值的阈值区间为
Figure BDA0003534082860000046
则火工品发生事故的概率P为:
Figure BDA0003534082860000047
式中,w(yl)和w(xl)分别为区间yl和xl的宽度,即
Figure BDA0003534082860000048
根据这两个区间的位置关系,式(3)等价地改写为:
Figure BDA0003534082860000049
根据式(4),第1种情况case1发生事故的概率为P=1,区间yl完全处于区间xl的右边,即区间yl总是大于区间xl;第3情况case3发生事故的概率为P=0,区间yl完全处于区间xl的左边,即区间yl总是小于区间xl;第2种情况case2发生事故概率P的取值范围为[0,1],两区间发生相互干涉的情况,当两区间的宽度不变时,其概率值由区间yl的上界
Figure BDA00035340828600000410
和区间xl的下界x的相对位置所决定。
进一步地,所述步骤S3包括:
根据下表所示严酷度等级的判别标准,在已知严酷度类别和后果严重程度的基础上,给出火工品事故严酷度等级S。
严酷度类别 后果严重程度描述 严酷度等级S
非常严重 灾难性伤亡和损失 ≥27
严重 重大伤亡和损失 [23,27)
中等严重 较大伤亡和损失 [19,23)
比较严重 一定伤亡和损失 [12,19)
轻度 较小伤亡和损失 <12
进一步地,所述步骤S4包括:
根据火工品发生事故的概率P和严酷度等级S,采用概率风险评价法对其进行安全性评估,则对应的风险值R为:
R=P×S。
下面结合具体实施例对本申请的方案作进一步说明:
电点火头用于DRS12电子时间等引信。该产品可靠度指标为:γ=0.90,R=0.999,分析其在电流输入强度为[22,32]内的安全性。
本发明提出的安全性评估方法的特点主要体现在火工品发生事故概率的确定方法上,所以主要从这方面进行对比分析。
(1)传统确定火工品事故发生概率的方法
按小样本方法的要求,先进行3组升降法试验,试验原始数据列于表1中。
表1 3组升降法试验数据
Figure BDA0003534082860000051
3号电点火头感度分布可假定为逻辑斯谛分布。利用极大似然估计原理求出参数估计值,见表2所示。
表2 3组升降法参数估计结果
Figure BDA0003534082860000052
根据表2,可得综合3组升降法试验的评估结果,可得均值为36.74,刻度参数为1.14。则在不同条件下,事故发生概率的计算结果见表3。
表3不同条件下事故发生概率
Figure BDA0003534082860000061
(2)采用本发明提出的方法确定事故发生概率
采用本发明提出的方法对激源强度在[22,34]范围内发生事故的概率进行评估,包含任意1组试验的感度均值的阈值区间为[32,41],根据公式(4)可得事故发生概率为:
Figure BDA0003534082860000062
可得通过本文提出方法获得的概率值与传统方法获得的概率值比较可得本文提出方法获得的值比较保守,且不依赖于分布模型的假设,对现场试验数据量要求较低,样本量少且评估精度高。
上述说明示出并描述了本发明的一个优选实施例,但如前所述,应当理解本发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述发明构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。

Claims (5)

1.一种基于阈值双区间的火工品安全性评估方法,其特征在于:包括以下步骤:
S1.给出激源强度和火工品感度均值的阈值区间;
S2.根据区间模型计算火工品发生事故的概率;
S3.根据火工品发生事故的严重程度,确定事故严酷度等级;
S4.根据火工品发生事故的概率和严酷度等级,采用概率风险评价法对其进行安全性评估。
2.根据权利要求1所述的一种基于阈值双区间的火工品安全性评估方法,其特征在于:所述步骤S1包括:
S101.假设激源强度和火工品感度均值不确定性参数均在某区间范围内变化,有:
Figure FDA0003534082850000011
Figure FDA0003534082850000012
Figure FDA0003534082850000013
Figure FDA0003534082850000014
式中,
Figure FDA0003534082850000015
Figure FDA0003534082850000016
为不确定性参数的上界,YyXx为不确定性参数的下界。
S102.利用区间表示法,将式(1)表示成如下的区间形式:
Figure FDA0003534082850000017
Figure FDA0003534082850000018
Figure FDA0003534082850000019
Figure FDA00035340828500000110
S103.设激源的强度区间为yl,火工品的感度均值区间为xl,发生事故的概率函数M为:
M=yl-xl
3.根据权利要求1所述的一种基于阈值双区间的火工品安全性评估方法,其特征在于:所述步骤S2包括:
当意外激源作用于火工品时,火工品发生事故的概率取决于意外激源的强度和火工品感度均值的大小;假设在某一时刻,只有一种激源作用于火工品;设激源强度的阈值区间为
Figure FDA00035340828500000111
火工品感度均值的阈值区间为
Figure FDA00035340828500000112
则火工品发生事故的概率P为:
Figure FDA00035340828500000113
式中,w(yl)和w(xl)分别为区间yl和xl的宽度,即
Figure FDA00035340828500000114
根据这两个区间的位置关系,式(3)等价地改写为:
Figure FDA00035340828500000115
根据式(4),第1种情况case1发生事故的概率为P=1,区间yl完全处于区间xl的右边,即区间yl总是大于区间xl;第3情况case3发生事故的概率为P=0,区间yl完全处于区间xl的左边,即区间yl总是小于区间xl;第2种情况case2发生事故概率P的取值范围为[0,1],两区间发生相互干涉的情况,当两区间的宽度不变时,其概率值由区间yl的上界
Figure FDA0003534082850000021
和区间xl的下界x的相对位置所决定。
4.根据权利要求1所述的一种基于阈值双区间的火工品安全性评估方法,其特征在于:所述步骤S3包括:
根据如下所示严酷度等级的判别标准,给出火工品事故严酷度等级S:
后果严重程度描述:灾难性伤亡和损失;严酷度类别:非常严重,严酷度等级S≥27;
后果严重程度描述:重大伤亡和损失;严酷度类别:严重,严酷度等级S在[23,27)之间;
后果严重程度描述:较大伤亡和损失;严酷度类别:中等严重,严酷度等级S在[19,23)之间;
后果严重程度描述:一定伤亡和损失;严酷度类别:比较严重,严酷度等级S在[12,19)之间;
后果严重程度描述:较小伤亡和损失;严酷度类别:轻度,严酷度等级S<12。
5.根据权利要求1所述的一种基于阈值双区间的火工品安全性评估方法,其特征在于:所述步骤S4包括:
根据火工品发生事故的概率P和严酷度等级S,采用概率风险评价法对其进行安全性评估,则对应的风险值R为:
R=P×S。
CN202210214960.2A 2022-03-07 2022-03-07 一种基于阈值双区间的火工品安全性评估方法 Pending CN114595960A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210214960.2A CN114595960A (zh) 2022-03-07 2022-03-07 一种基于阈值双区间的火工品安全性评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210214960.2A CN114595960A (zh) 2022-03-07 2022-03-07 一种基于阈值双区间的火工品安全性评估方法

Publications (1)

Publication Number Publication Date
CN114595960A true CN114595960A (zh) 2022-06-07

Family

ID=81815727

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210214960.2A Pending CN114595960A (zh) 2022-03-07 2022-03-07 一种基于阈值双区间的火工品安全性评估方法

Country Status (1)

Country Link
CN (1) CN114595960A (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008269354A (ja) * 2007-04-20 2008-11-06 Toshiba Corp インシデント・アクシデントレポート分析装置、方法、およびプログラム
CN107093022A (zh) * 2017-04-24 2017-08-25 中国工程物理研究院化工材料研究所 高置信度的燃爆产品可靠性的小样本评估方法
CN110458436A (zh) * 2019-08-02 2019-11-15 北京理工大学 一种火工品可靠性指标的小样本验证方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008269354A (ja) * 2007-04-20 2008-11-06 Toshiba Corp インシデント・アクシデントレポート分析装置、方法、およびプログラム
CN107093022A (zh) * 2017-04-24 2017-08-25 中国工程物理研究院化工材料研究所 高置信度的燃爆产品可靠性的小样本评估方法
CN110458436A (zh) * 2019-08-02 2019-11-15 北京理工大学 一种火工品可靠性指标的小样本验证方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张跃 等: "可靠度的区间算法", 化学工程与装备, no. 11, 30 November 2020 (2020-11-30), pages 198 - 199 *

Similar Documents

Publication Publication Date Title
CN111046564B (zh) 两阶段退化产品的剩余寿命预测方法
WO2020186718A1 (zh) 一种直流电缆绝缘状态判断方法及系统、存储介质
CN111523251A (zh) 一种随机环境应力下的产品寿命快速评估方法
CN107220907B (zh) 一种采用秩和比综合评价的谐波污染用户分级方法
KR20180041727A (ko) 신용점수 모델 훈련 방법,신용 점수 계산 방법,장치 및 서버
Ma et al. Optimum step‐stress accelerated life test plans for log‐location‐scale distributions
CN113779714A (zh) 一种焊接接头用p_s_n曲线测定方法、装置及系统
CN113360983B (zh) 一种边坡可靠度分析与风险评估方法
CN114595960A (zh) 一种基于阈值双区间的火工品安全性评估方法
CN106570616A (zh) 一种科技项目评估用定量评价方法
EP3842984A1 (en) Probabilistic simulation method of a microelectronic device
CN109670593B (zh) 一种评估、以及预测深度学习模型中层计算时间的方法
CN103852562A (zh) 判断样品检测数据可疑值的方法
CN114236332A (zh) 一种电力电缆绝缘状态判断方法及系统
EP2815246B1 (de) Verfahren zum bestimmen eines stroms, batteriemanagementeinheit, batterie und kraftfahrzeug
Dave et al. Predicting diabetic neuropathy risk level using artificial neural network based on clinical characteristics of subjects with diabetes
Mauro et al. Extreme values calculation of multi-modal peak distributions
CN112598259A (zh) 产能测算方法、装置及计算机可读存储介质
CN112132435A (zh) 一种基于半定量计算的化工风险评价方法
Knofczynski Sample sizes for predictive regression models and their relationship to correlation coefficients
CN104346534A (zh) 一种测量随机数据非指数性的熵方法
CN114386525B (zh) 一种评估标准的自适应调整方法、装置、设备和存储介质
CN112859793B (zh) 基于改进滑动时间窗口的工业生产过程动态时延辨识方法
CN113404742B (zh) 一种基于测试数据的电液伺服机构健康评估方法及系统
CN112577461B (zh) 一种基于挠度分离的大跨径梁桥状态的预测方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination