(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202210924553.0 (22)申请日 2022.08.02 (71)申请人 西安航天动力研究所 地址 710100 陕西省西安市航天基地飞 天 路289号 (72)发明人 张航　陈晖　高玉闪　龚南妮　 李程　张淼　苏展　 (74)专利代理 机构 西安智邦专利商标代理有限 公司 6121 1 专利代理师 王少文 (51)Int.Cl. G06F 30/20(2020.01) G06F 17/12(2006.01) (54)发明名称 一种基于大型方程组求解的液体火箭发动 机工况调整方法 (57)摘要 本发明提供的一种基于大型方程组求解的 液体火箭发动机工况调整 方法， 用以解决当初始 值与正确的求解结果间相差较大时， 大型方程组 出现无法正确求解， 导致无法对 液体火箭发动机 的工况进行准确及时调整的技术问题。 本方法包 括： 建立工 况调整计算模型； 将100％推力工况参 数作为初始值， 输入大型非线性方程组中； 设置 待求解工况并进行求解； 当求解顺利时， 将当前 结果作为发动机工况调整目标的对应参数； 否 则， 增加一个中间推力工况作为发动机工况调整 目标进行求解， 直至求解顺利， 将当前辅助待求 解工况参数作为初始值输入重新求解， 以此类推 直至求解顺利， 将当前结果作为发动机工况调整 目标的对应参数； 对液体 火箭发动机进行工况调 整。 权利要求书2页 说明书5页 附图5页 CN 115329557 A 2022.11.11 CN 115329557 A 1.一种基于大型方程组求解的液体火箭发动机工况调整方法， 其特征在于， 包括以下 步骤： 1】 建立液体火箭发动机 工况调整计算模型 根据发动机组件特性和物理连接关系， 建立液体火箭发动机工况调整计算模型； 所述 液体火箭发动机 工况调整计算模型为大 型非线性方程组； 2】 输入初始值 将100％推力工况中的参数作为初始值， 输入步骤1】 建模得到的大 型非线性方程组中； 3】 设置待求 解工况， 并对大 型非线性方程组进行求 解 将液体火箭发动机工况调 整目标作为待求解工况， 在输入了初始值的大型非线性方程 组中设置待求 解工况的参数， 并对大 型非线性方程组进行求 解； 4】 获得发动机 工况调整目标的对应参数 当大型非线性方程组 的求解顺利时， 将当前求解结果作为发动机工况调整目标的对应 参数； 否则， 在当前待求解工况与100％推力工况之间增加一个中间推力工况作为辅助的待 求解工况， 并返回步骤3】 ， 以辅助的待求解工况作为 发动机工况调整目标， 对大型非线性方 程组求解， 直至求解顺利， 得到 当前辅助的待求解工况的参数； 将当前辅助的待求解工况的 参数作为初始值输入大型非线性方程组， 将液体火箭发动机待求解工况作为工况调整目 标， 重新对大型非线性方程组求解， 以此类推， 直至求解顺利， 将当前求解结果作为发动机 工况调整目标的对应参数； 5】 根据步骤4】 获得的发动机工况调整目标的对应参数， 对液体火箭发动机进行工况调 整。 2.根据权利要求1所述的一种基于大型方程组求解的液体火箭发动机工况调整方法， 其特征在于： 步骤4】 中， 大型非线性方程组的求解顺利进行具体为： 大型非线性方程组的求解结果 正常， 未出现求解速度慢、 多解或发散的问题。 3.根据权利要求2所述的一种基于大型方程组求解的液体火箭发动机工况调整方法， 其特征在于： 步骤1】 中， 建立液体火箭发动 机工况调整计算模型具体为： 采用Modelica语言， 建立液 体火箭发动机 工况调整计算模型。 4.根据权利要求3所述的一种基于大型方程组求解的液体火箭发动机工况调整方法， 其特征在于： 步骤2】 具体为： 2.1】 对步骤1】 建模得到的大型非线性方程组进行解耦， 确定大型非线性方程组中的方 程分组以及各分组的求 解顺序； 2.2】 通过 数值求解算法确定待输入初始值的参数类别； 2.3】 根据步骤2.2】 确定的待输入初始值的参数类别， 将100％推力工况中对应的参数 作为初始值， 输入大 型非线性方程组中。 5.根据权利要求4所述的一种基于大型方程组求解的液体火箭发动机工况调整方法， 其特征在于： 步骤3】 中， 对大型非线性方程组进行求解具体为： 根据步骤2.1】 确定的方程分组以及权　利　要　求　书 1/2 页 2 CN 115329557 A 2各分组的求 解顺序， 对大 型非线性方程组进行求 解。权　利　要　求　书 2/2 页 3 CN 115329557 A 3