CN114595960A - 一种基于阈值双区间的火工品安全性评估方法 - Google Patents
一种基于阈值双区间的火工品安全性评估方法 Download PDFInfo
- 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
Links
- 239000002360 explosive Substances 0.000 title claims abstract description 80
- 230000000977 initiatory effect Effects 0.000 title claims abstract description 80
- 238000011156 evaluation Methods 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 claims abstract description 30
- 230000035945 sensitivity Effects 0.000 claims abstract description 24
- 230000005284 excitation Effects 0.000 claims abstract description 17
- 230000035939 shock Effects 0.000 claims description 6
- 101100213915 Bacillus subtilis (strain 168) yorX gene Proteins 0.000 claims description 3
- 150000001875 compounds Chemical class 0.000 claims description 3
- 230000006378 damage Effects 0.000 claims 1
- 238000012360 testing method Methods 0.000 abstract description 12
- 238000011076 safety test Methods 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 229940037201 oris Drugs 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION 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/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex 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.假设激源强度和火工品感度均值不确定性参数均在某区间范围内变化,则有:
S102.利用区间表示法,将式(1)表示成如下的区间形式:
S103.设激源的强度区间为yl,火工品的感度均值区间为xl,发生事故的概率函数M为:
M=yl-xl。
进一步地,所述步骤S2包括:
当意外激源作用于火工品时,火工品发生事故的概率取决于意外激源的强度和火工品感度均值的大小;假设在某一时刻,只有一种激源作用于火工品;设激源强度的阈值区间为火工品感度均值的阈值区间为则火工品发生事故的概率P为:
式中,w(yl)和w(xl)分别为区间yl和xl的宽度,即
根据这两个区间的位置关系,式(3)等价地改写为:
根据式(4),第1种情况case1发生事故的概率为P=1,区间yl完全处于区间xl的右边,即区间yl总是大于区间xl;第3情况case3发生事故的概率为P=0,区间yl完全处于区间xl的左边,即区间yl总是小于区间xl;第2种情况case2发生事故概率P的取值范围为[0,1],两区间发生相互干涉的情况,当两区间的宽度不变时,其概率值由区间yl的上界和区间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.假设激源强度和火工品感度均值不确定性参数均在某区间范围内变化,则有:
S102.利用区间表示法,将式(1)表示成如下的区间形式:
S103.设激源的强度区间为yl,火工品的感度均值区间为xl,发生事故的概率函数M为:
M=yl-xl。
进一步地,所述步骤S2包括:
当意外激源作用于火工品时,火工品发生事故的概率取决于意外激源的强度和火工品感度均值的大小;假设在某一时刻,只有一种激源作用于火工品;设激源强度的阈值区间为火工品感度均值的阈值区间为则火工品发生事故的概率P为:
式中,w(yl)和w(xl)分别为区间yl和xl的宽度,即
根据这两个区间的位置关系,式(3)等价地改写为:
根据式(4),第1种情况case1发生事故的概率为P=1,区间yl完全处于区间xl的右边,即区间yl总是大于区间xl;第3情况case3发生事故的概率为P=0,区间yl完全处于区间xl的左边,即区间yl总是小于区间xl;第2种情况case2发生事故概率P的取值范围为[0,1],两区间发生相互干涉的情况,当两区间的宽度不变时,其概率值由区间yl的上界和区间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组升降法试验数据
3号电点火头感度分布可假定为逻辑斯谛分布。利用极大似然估计原理求出参数估计值,见表2所示。
表2 3组升降法参数估计结果
根据表2,可得综合3组升降法试验的评估结果,可得均值为36.74,刻度参数为1.14。则在不同条件下,事故发生概率的计算结果见表3。
表3不同条件下事故发生概率
(2)采用本发明提出的方法确定事故发生概率
采用本发明提出的方法对激源强度在[22,34]范围内发生事故的概率进行评估,包含任意1组试验的感度均值的阈值区间为[32,41],根据公式(4)可得事故发生概率为:
可得通过本文提出方法获得的概率值与传统方法获得的概率值比较可得本文提出方法获得的值比较保守,且不依赖于分布模型的假设,对现场试验数据量要求较低,样本量少且评估精度高。
上述说明示出并描述了本发明的一个优选实施例,但如前所述,应当理解本发明并非局限于本文所披露的形式,不应看作是对其他实施例的排除,而可用于各种其他组合、修改和环境,并能够在本文所述发明构想范围内,通过上述教导或相关领域的技术或知识进行改动。而本领域人员所进行的改动和变化不脱离本发明的精神和范围,则都应在本发明所附权利要求的保护范围内。
Claims (5)
1.一种基于阈值双区间的火工品安全性评估方法,其特征在于:包括以下步骤:
S1.给出激源强度和火工品感度均值的阈值区间;
S2.根据区间模型计算火工品发生事故的概率;
S3.根据火工品发生事故的严重程度,确定事故严酷度等级;
S4.根据火工品发生事故的概率和严酷度等级,采用概率风险评价法对其进行安全性评估。
3.根据权利要求1所述的一种基于阈值双区间的火工品安全性评估方法,其特征在于:所述步骤S2包括:
当意外激源作用于火工品时,火工品发生事故的概率取决于意外激源的强度和火工品感度均值的大小;假设在某一时刻,只有一种激源作用于火工品;设激源强度的阈值区间为火工品感度均值的阈值区间为则火工品发生事故的概率P为:
式中,w(yl)和w(xl)分别为区间yl和xl的宽度,即
根据这两个区间的位置关系,式(3)等价地改写为:
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。
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)
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 | 北京理工大学 | 一种火工品可靠性指标的小样本验证方法 |
-
2022
- 2022-03-07 CN CN202210214960.2A patent/CN114595960A/zh active Pending
Patent Citations (3)
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)
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 |