固体火箭发动机内弹道多参数辨识及精准预示方法.pdf

第 卷第 期 固 体 火 箭 技 术 固体火箭发动机内弹道多参数辨识及精准预示方法张 凌，王德友，李世鹏，王宁飞（北京理工大学 宇航学院，北京）摘要：提高固体火箭发动机内弹道性能的预示精度对提升火箭导弹工作可靠性有着重要意义。实践中常根据发动机实测压强数据反算燃面，并结合实测燃速进行内弹道预示，其在预示精度上仍有较大优化空间。针对该问题，提出了一种发动机内弹道多参数辨识及精准预示方法，基于有限次试车实测数据，通过遗传算法对燃速系数、压强指数、推进剂密度和喷管喉衬烧蚀速率系数开展初步辨识，基于辨识结果逆向反算出燃面数据并更新模型，再开展二次辨识对关键参数进行修正。以美国 型发动机和国内某型长通道发动机为案例进行计算分析，对发动机内弹道性能进行预示。计算结果表明，该方法所得预示结果与实测压强数据高度吻合，验证了该方法对提高内弹道性能预示精度的有效性。关键词：固体火箭发动机；内弹道；参数辨识；性能预示；遗传算法中图分类号：文献标识码：文章编号：（）：，（，）：（），（），：；引言固体火箭发动机内弹道学的基本任务是在发动机各种工作条件下计算燃烧室压强随时间和空间的变化规律。零维内弹道计算方法假定燃烧室内参数均匀一致，与一维及多维内弹道计算方法相比更加简便、直观，是工程上应用较多的求解方法。发动机的内弹道特性直接决定着整机性能，进而影响着火箭导弹的飞行特性。因此，进一步提高内弹道预示的精度，对保障飞行器工作可靠性有着重要意义。工程上通常根据单台或多台发动机的实测数据和收稿日期：；修回日期：。作者简介：张凌，男，硕士，研究方向为固体火箭发动机优化设计。通讯作者：李世鹏，男，博士，研究方向为固体火箭发动机燃烧与流动特性。燃面几何退移规律，对燃速特性等参数进行辨识，以获得更精确的内弹道性能预示结果。目前，业内学者针对发动机内弹道性能精确预示问题，从理论算法修正、参数辨识方法改进等角度开展了研究。如刘杨等基于对发动机实测数据处理误差的理论分析，提出了修正的燃速系数处理方法，该方法可提高具有较高压强比的发动机内弹道预示精度。李春艳等针对短时间工作固体发动机特点，从理论角度对现有内弹道计算方法进行了适应性修改，所得预示结果与实测数据相比具有较高精度。在辨识参数选择方面，等通过研究 固体运载火箭发动机实际飞行性能影响因素，发现喷管喉部烧蚀率对推力和燃烧室压强的影响较大。等采取数值仿真方法研究喷管喉部烧蚀率和燃速变化对燃烧室压强的影响，发现这两个参数的改变对压强时间曲线影响十分明显。李晓斌等基于发动机工作过程的内弹道数学模型及实测燃烧室压强时变数据，通过参数辨识方法可对燃速、喉衬烧蚀率等参数进行计算。在参数辨识方法改进方面，由于最小二乘法、最大似然法等传统辨识方法常存在初始点敏感等问题，因而如遗传算法等具有初值稳健性和全局搜索性的方法在内弹道参数辨识领域逐渐得到深入应用。樊超等将遗传算法用于固体火箭发动机参数辨识，获得了推进剂燃速模型和喉径变化模型的全局最优辨识值，所得内弹道性能预示结果与试验测量值吻合良好。庄建华等对发动机贮存老化内弹道性能的影响因素进行参数建模，基于拉丁超立方设计方法进行了瞬态参数辨识，结合非定常内弹道模型对发动机贮存老化内弹道进行了预示。蒲晓航等基于部分实测压强数据，通过遗传算法逆向反算获得预示参数，再结合正向计算实现了内弹道的预示自修正，对复杂环境剖面下的发动机内弹道性能进行了准确预示。孙瑞阳等建立了基于经典 卷积神经网络的内弹道参数辨识方法，对长期贮存发动机推进剂装药的燃速特性参数进行了辨识。张小平等建立了高能固体推进剂的配方表征方法，利用遗传神经网络（）对推进剂高压燃烧性能及其变化规律进行了预示。传统发动机的实测数据处理常通过燃烧时间内的平均压强和药柱平均燃速计算获得燃速系数等参数并反算燃面，而实际上平均压强与平均燃速间并不符合燃速模型对应关系，因而计算存在一定误差。而目前的研究工作大多采用遗传算法和神经网络等方法对关键内弹道参数进行辨识，并没有根据实测数据对药柱燃面以及内弹道参数进行修正与二次辨识，内弹道参数辨识以及性能预示的精度还可以进一步提高。因此，有必要对关键内弹道参数开展二次辨识，得到更加精确的辨识结果，进而提高发动机内弹道预示精度。本文提出了一种发动机内弹道多参数辨识及性能精准预示方法，根据发动机实测压强数据，通过遗传算法对燃速系数、压强指数、推进剂密度和喷管喉衬烧蚀速率系数开展初步辨识，基于辨识结果修正燃面肉厚关系，在此基础上开展二次参数辨识，最终确定关键内弹道参数。根据第二次参数辨识结果和修正后的燃面肉厚关系，通过正向预示可获得发动机的内弹道性能。以美国 型发动机和国内某型长通道发动机为案例进行计算分析，验证了该方法内弹道性能预示的精度和有效性，相关结果可为固体火箭发动机的优化设计提供技术参考。内弹道预示模型本文以零维内弹道模型为基础进行发动机性能正向预示，该模型有如下假设：（）燃气为理想气体，气体流动为零维定常流动；（）装药点燃过程瞬时完成；（）喷管喉衬烧蚀速率与燃烧室压强正相关。在考虑点火药和喷管喉衬烧蚀作用对内弹道计算影响的情况下，结合燃烧室内质量守恒定律、能量守恒定律、燃烧室内体积守恒方程、燃速方程以及气体状态方程，推导出的内弹道微分方程组如下：（）（）（）其中，（）（）（）年 月张凌，等：固体火箭发动机内弹道多参数辨识及精准预示方法第 期（）（）式中为燃烧室内燃气密度；为燃烧室自由容积；为装药燃气生成速率；为点火药燃气生成速率；为喷管燃气流率；为燃烧室压强；为燃气比热容比；为热损失系数；为燃气气体常数；为推进剂绝热燃烧温度；为点火药的火药力；为燃烧室内燃气温度；为装药燃烧面积；为装药燃烧速率；为点火药燃烧面积；为点火药燃烧速率；为喷管喉径；为喷管喉衬线烧蚀速率；为喷管喉衬烧蚀速率系数；为燃速系数；为燃速压强指数；为推进剂密度；为点火药密度；为点火药燃去肉厚；为点火药颗粒直径；为比热容比函数。多参数辨识及精准预示方法考虑到发动机试验成本问题，基于一次或几次试车数据开展发动机内弹道关键参数的精确辨识研究，对降低固体发动机设计成本和提高设计效率有着重要意义。基于工程中的内弹道参数辨识方法，本文提出一种发动机内弹道多参数辨识及预示方法，可提高内弹道预示精度。工程中内弹道参数辨识方法工程上对内弹道参数的辨识，大体是根据实测燃烧室压强曲线、试车后喉部面积等数据计算燃速参数并反算燃面。常通过式（）计算平均燃速，根据式（）计算燃速系数，通过燃速系数以及式（）来进行燃面反算，获得燃面肉厚数据，最后根据这些参数进行内弹道正向预示，得到发动机性能预示曲线。（）（）（）式中 为装药肉厚；为燃烧时间；为燃面面积；为喷管喉部面积；为特征速度。多参数辨识及精准预示方法流程上述平均燃速和平均压强并不具有实际对应关系，计算存在一定误差，在内弹道预示精度上仍有较大优化空间。本文提出了一种内弹道多参数辨识及精准预示方法，该方法是对工程办法的改进。以单台发动机的实测压强数据为基础，在发动机稳定工作段选取部分时间段进行计算，基于所建立辨识修正方法获得燃速模型参数和燃面肉厚数据等内弹道预示所需全部关键参数，进行内弹道性能预示，该方法能提高内弹道预示精度，具体流程如图 所示。performance图 多参数辨识及精准预示方法流程 从发动机实测压强数据出发，基于 遗传算法对燃速系数、压强指数、推进剂密度 及喷管喉衬烧蚀速率系数 进行第一次参数辨识。通过参数辨识结果和燃面压强公式对燃面肉厚曲线进行反算，进而更新多参数辨识的计算模型，再通过所得燃面数据和遗传算法开展第二次参数辨识。根据二次辨识结果对初次辨识参数进行修正，在第一章介绍的内弹道预示模型的基础上，根据所得参数及燃面数据对发动机内弹道性能开展正向预示。基于遗传算法的内弹道多参数辨识内弹道多参数辨识主要根据发动机工作时间内的实测压强数据进行逆向反算，该工作基于 遗传算法完成，对燃速系数、压强指数、推进剂密度 及喷管喉衬烧蚀速率系数 四个参数进行辨识，具体步骤如下：（）根据实测数据信息及设计经验，选定时间间隔，作为计算两条曲线拟合程度的最小时间微元；（）给定参数初值范围，其中燃速系数 和压强指数 还需满足燃速模型关系，故需施加额外限制条件；（）在给定初值范围内，随机给出 组数据（，）作为模型初始群体，每组数据为初始群体的个体；（）根据零维内弹道理论建立正向预示数学模型，基于任意一组数据可进行一次内弹道性能预示，即 组数据可获得 条压强曲线；（）由于发动机稳定工作段数据更具代表性，故在稳定工作段内以 为间隔总计选取 个节点，建立内弹道预示模型的适应度函数，其为 个节点处预示压强与实测压强的差值组成的行向量，该行向量的 年 月固体火箭技术第 卷范数作为衡量该组数据适应度高低的判据；（）淘汰适应度低的初始个体，保留适应度高的初始个体，在保留个体间进行数据的交换和变异，形成初始群体的子代，再次通过适应度函数对新群体进行选择；（）重复上步操作直至结果收敛，将适应度最高的数据代入内弹道预示模型中，即可获得与实测数据吻合度最高的压强曲线。上述为初次参数辨识的步骤。基于初次参数辨识结果进行燃面反算，根据反算燃面数据更新内弹道预示数学模型，再进行二次参数辨识修正。第二次辨识的参数初值范围可适当缩小，以加快收敛速度并提高效率。案例计算分析为验证内弹道性能预示方法的准确性和有效性，本文选取 型发动机和一种长通道发动机作为研究案例。型发动机案例 型发动机的装药组成为 和 的复合推进剂，该发动机结构参数和试车数据较为详尽，许多关于发动机燃烧稳定性和内流场特性的研究均选取该发动机作为验证案例。该型发动机的具体参数如表 所示，结构如图 所示。表 发动机参数 ，（），图 发动机结构图 第一次参数辨识多参数辨识过程基于 软件实现，通过遗传算法逆向反算辨识出最优的一组参数据（，）。基于零维内弹道预示模型得到预示压强曲线，将预示压强曲线与实测压强曲线的拟合度作为适应度函数，求解适应度函数的最小值。初次参数辨识的适应度函数建立步骤如下：（）取最小时间微元 进行曲线拟合度对比计算；（）在稳定工作段内，以 为时间间隔选取 个点作为曲线拟合度计算时间节点；（）计算所有节点上预测压强值与实测压强值的差向量，将该行向量的 范数作为适应度函数。发动机实测压强曲线存在初始尖峰，根据工程经验可知其为点火药与装药共同燃烧导致产气量剧增所致。点火过程影响因素较为复杂，但其对内弹道曲线整体趋势影响并不显著。辨识计算中取发动机稳定工作部分时间段作为计算区间，避开初始点火过程和压强下降拖尾段，以提高多参数辨识的精度。在第一次多参数辨识的过程中，设定 项待辨识参数的初始取值范围分别为 （）（），（）（）。基于遗传算法进行多参数辨识，在经过多次迭代计算后结果收敛，求解过程如图 所示。初次参数辨识结果为：燃速系数 （），压强指数 ，推进剂密度 ，喷管喉衬烧蚀速率系数 （）。图 第一次遗传算法优化求解过程 由于原始文献未给出发动机点火装置参数，且点火过程对宏观内弹道特性影响较小，故将点火药假定为球形颗粒黑火药，再基于遗传算法和考虑点火过程的内弹道预示模型对点火药质量和颗粒粒径进行辨识。该过程所需其余参数采用初次参数辨识的结果，最终获得点火药质量为 ，粒径为 。根据参 年 月张凌，等：固体火箭发动机内弹道多参数辨识及精准预示方法第 期数辨识结果进行内弹道预示，可获得初步预示的压强时间曲线，如图 所示。1.81.51.20.90.60.30.00 1 23 4t/sp/MPaMeasured pressure dataFirst multi-parameter identification图 第一次辨识计算结果与实测数据对比 燃面反算及第二次参数辨识基于初次参数辨识结果，根据燃面压强公式进行燃面反算。为避免点火过程对反算燃面数据产生影响，在计算时不考虑点火药燃烧的影响，可获得如图 所示的反算燃面肉厚曲线。图 反算燃面肉厚曲线 在初次辨识基础上开展二次辨识，基于反算燃面数据对已有数学模型进行更新。为加速收敛，第二次辨识的参数初值范围相较第一次辨识适当缩小。设定 项待辨识参数的初始取值范围分别为 （）（），（）（）。如图 所示，经过多次迭代计算后结果收敛，第二次参数辨识结果为：燃速系数 （），压强指数 ，推进剂密度 ，喷管喉衬烧蚀率系数 （）。根据反算燃面数据和第二次参数辨识结果再次进行内弹道预示，可获得图 所示的修正压强时间曲线。由图 可知，第二次性能预示精度较第一次更高，且预示结果与实测数据高度吻合。图 第二次遗传算法优化求解过程 图 第二次辨识计算结果与实测数据对比 某型长通道发动机案例测试案例二为某型长通道固体发动机，该发动机装药为端面包覆的内外燃管形装药，采用双基推进剂并以黑火药为点火药，其相关参数较为详尽，具体结构参数如表 所示，结构图如图 所示表 长通道发动机参数 ，年 月固体火箭技术第 卷图 长通道发动机结构图 应用上述介绍的内弹道多参数辨识方法对该长通道发动机进行第一次参数辨识，第一次参数辨识结果为：燃速系数 （），压强指数，推进剂密度 ，喷管喉衬烧蚀率系数 （）。用第一次多参数辨识结果进行燃面反算，根据内弹道预示模型进行性能预示，计算得到如图 所示的燃烧室压强预示曲线和如图 所示反算燃面曲线。由于燃面数据计算采用式（）所示的工程方法，故上升段和下降段与实际燃面变化的物理过程存在一定差异。10864200 1 2 3 45 6t/sp/MPaMeasured pressure dataFirst multi-parameter identification图 长通道发动机第一次辨识预示曲线 图 长通道发动机反算燃面肉厚曲线 根据反算燃面数据更新模型后，进行第二次多参数辨识，得到最终参数辨识结果为：燃速系数 （），压强指数 ，推进剂密度，喷管喉衬烧蚀率系数 （）。由于辨识计算过程的主要依据为发动机稳定工作段压强数据，故所得辨识参数结果具有较好的合理性。把反算的燃面数据以及修正后的参数代入内弹道正向预示模型进行计算，可获得如图 所示的燃烧室压强曲线。由图 可知，第二次性能预示精度较第一次更高，且预示结果与实测数据高度吻合。图 长通道发动机第二次辨识预示曲线 结论（）提出了固体火箭发动机内弹道多参数辨识及精准预示方法，并以 型发动机和某型长通道发动机为研究案例进行计算分析，基于二次参数辨识修正和燃面反算，获得了与实测压强数据高度吻合的预示结果，验证了该方法对提高内弹道性能预示精度的有效性。（）除本文案例参数外，还可针对发动机使用环境剖面下具有较强散布特性的多项预示参数开展辨识修正，相关方法也适用于具有复杂装药构型和燃面退移规律的固体发动机，其应用具有较强的普遍性。（）基于该内弹道参数辨识与预示方法，可充分利用有限次发动机试车实测数据开展内弹道性能精准预示，掌握复杂环境剖面下的内弹道性能散布情况，提高发动机设计效率并降低成本。参考文献：唐金兰，刘佩进固体火箭发动机原理北京：国防工业出版社，：鲍福廷，侯晓固体火箭发动机设计北京：中国宇航出版社，：年 月张凌，等：固体火箭发动机内弹道多参数辨识及精准预示方法第 期 刘杨，丁淼，白彦军，等实测燃速处理误差对内弹道性能预示的影响分析固体火箭技术，（）：，（）：李春艳，王立武，刘沛，等短时间工作发动机内弹道计算方法弹箭与制导学报，（）：，（）：，李晓斌，王中伟，张为华推进剂稳态燃速最优化辨识固体火箭技术，（）：，（）：樊超，张为华基于遗传算法的固体火箭发动机参数辨识固体火箭技术，（）：，（）：庄建华，张哲伟，樊超，等固体发动机贮存条件下内弹道性能研究弹箭与制导学报，（）：，（）：蒲晓航，李冬，李富贵，等固体火箭发动机内弹道精准预示自修正方法固体火箭技术，（）：，（）：孙瑞阳，姜毅，牛钰森，等基于卷积神经网络的固体火箭发动机内弹道参数辨识固体火箭技术，（）：，（）：张小平，代志龙基于遗传神经网络的高能固体推进剂高压燃烧性能计算固体火箭技术，（）：，（）：，（）：，：郜冶，刘平安，胡伟长通道固体火箭发动机点火瞬态数值分析哈尔滨工程大学学报，（）：，（）：刘赟，王浩，陶如意，等点火过程对小型固体火箭发动机内弹道影响含能材料，（）：，（）：张智慧，李军伟，梅开，等多根管型装药固体火箭发动机初始压强峰影响因素研究推进技术，（）：，（）：（编辑：吕耀辉）年 月固体火箭技术第 卷