(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202210867418.7 (22)申请日 2022.07.22 (71)申请人 中国航天标准 化研究所 地址 100071 北京市丰台区小屯路89号 (72)发明人 李健　何鹏飞　刘金燕　王珏　 (74)专利代理 机构 北京理工大 学专利中心 11120 专利代理师 李微微 (51)Int.Cl. G06F 30/20(2020.01) G06F 17/18(2006.01) G06F 119/02(2020.01) (54)发明名称 一种液体火箭发动机多机同步可靠性评估 方法 (57)摘要 本发明公开了一种液体火箭发动机多机同 步可靠性评估 方法， 针对液体 火箭发动机多机同 步可靠性指标要求， 借助蒙特卡罗仿真算法， 在 保证计算精度的条件下， 对性能参数的分布类型 不再限制， 适用于各种分布类型参数的多机同步 可靠性评估， 能够扩充方法的适用范围； 同时能 够简化计算流程， 提高多机同步可靠性评估效 率。 从可靠性区间估计的原理出发， 整个计算过 程均未采用近似算法， 在足够抽样 仿真次数的情 况下， 能够保证 计算精度要求。 权利要求书2页 说明书7页 附图1页 CN 115391986 A 2022.11.25 CN 115391986 A 1.一种液体火箭发动机多机同步可靠性评估方法， 其特 征在于， 包括如下步骤： 步骤1)、 液体火箭发动机多机同步可靠性指标要求， 收集发动机的性能参数测量值数 据， 根据性能参数测量 值数据， 进行统计 推断， 判断数据的概 率分布类型； 步骤2)、 基于发动机的性能参数测量值数据， 计算概率分布的分布参数， 得到发动机关 键性能参数的概 率密度函数； 步骤3)、 确定循环次数M1； 步骤4)、 确定可靠性评估仿真次数M2； 步骤5)、 针对每次仿真， 根据概率密度函数生成与多机个数N数量一致的随机数rj； j＝ 1,2,…,M2； 步骤6)、 针对每次仿真， 计算 最大偏差: Dev(j)＝max(rj)‑min(ri) 其中， max(rj)和min(rj)分别表示rj中的最大值和最小值； 步骤7)、 判断最大偏差值Dev(j)是否大于多机同步要求的最大偏差上限值， 若大于， 则 第j次仿真的结果记为失效， 失效数fai l累积一次； 反 之则记为成功， 成功数suc c累积一次； 步骤8)、 重复执 行步骤5)至步骤7)M2次， 计算得到成功数suc c和失效数fai l； 步骤9)、 根据成功数suc c和失效数fai l， 开展多机同步可靠性估计， 计算可靠度： 式中， Ri为第i次循环计算的可靠度； 步骤10)、 重复执行步骤5)至步骤10)， 共计M1次， 对应得到M1个可靠度， 将计算得到的M1 个可靠度降序排序， 取序号K＝M1×γ排序的可靠度RK值为置信 下限为γ的可靠度RL； γ为 设定的置信下限； 根据可靠度RL对液体火箭发动机多机同步可靠性评估。 2.如权利要求1所述的一种液体火箭发动机多机同步可靠性评估方法， 其特征在于， 所 述步骤1)中， 根据性能参数测量值数据对常用的概率分布类型求取对数似然值， 则选择对 数似然值 最小的概 率分布类型， 作为 性能参数测量 值数据对应的概 率分布类型。 3.如权利要求1所述的一种液体火箭发动机多机同步可靠性评估方法， 其特征在于， 所 述步骤1)中， 根据性能参数测量值数据对常用的概率分布类型求取对数似然值， 则选择对 数似然值最小的概率分布类型进行假设检验， 如果接受原假设， 则该组数据服从之前确定 的概率分布类型； 若拒绝原假设， 选择对数似然值第二小的概率分布类型， 并进行假设检 验， 如果接受假设检验， 则属于该对数似然值第二小的概率分布类型； 以此类推， 直到找到 所属的分布类型。 4.如权利要求3所述的一种液体火箭发动机多机同步可靠性评估方法， 其特征在于， 所 述步骤1)中， 但概 率分布类型为 正态分布时， 采用Jarque ‑Bera方法进行假设检验。 5.如权利要求1所述的一种液体火箭发动机多机同步可靠性评估方法， 其特征在于， 所 述步骤3)中， 循环次数M1的确定方法为： 保证设置的置信下限γ与循环次数M1的乘积为整 数， 且为10的倍数， 满足该 条件的M1值再扩大10倍， 则为 最终的循环次数M1。 6.如权利要求1所述的一种液体火箭发动机多机同步可靠性评估方法， 其特征在于， 所 述步骤5)中， 采用反函数法生成随机数。权　利　要　求　书 1/2 页 2 CN 115391986 A 27.如权利要求1所述的一种液体火箭发动机多机同步可靠性评估方法， 其特征在于， 所 述设定的置信下限γ选取 70％。权　利　要　求　书 2/2 页 3 CN 115391986 A 3