[发明专利]固体火箭发动机装药设计方法

摘要

本发明提供一种固体火箭发动机装药设计方法，以代理模型为基础，直接对推力曲线进行近似，通过构造推力曲线而不是推力曲线与设计指标之间的最小二乘偏差的代理模型，有效刻画出推力随时间的变化规律，从而明显减少了高精度燃面与内弹道仿真的次数。代理模型能够更好地刻画推力曲线的变化规律，在后续的设计中更有效指导后续搜索，能显著提高固体火箭发动机装药设计效率。

主权项

1.一种固体火箭发动机装药设计方法，其特征在于，包括以下步骤：1)给定发动机装药几何构型设计指标要求的推力需求曲线F0(t)；2)建立发动机装药几何构型的参数化模型，根据所处理装药几何构型的类型确定所需处理的设计变量X及其范围；3)根据设计变量X的个数m及其范围建立设计空间，在设计空间内采用最优拉丁超立方采样法采集2m个采样点，在每个采样点处建立性能仿真模型，并运行性能仿真模型，得到每个采样点对应的2m条推力曲线，所得设计变量X在第i个采样点处的值Xi和与第i个采样点处的仿真推力曲线fi(t)之间的对应关系，如式(1)所示；其中，Xi为设计变量X在第i个采样点处的推力值，fi(t)为对第i个采样点处的各设计变量分别进行设计仿真得到的仿真推力曲线，将2m条仿真推力曲线在每个工作时间点t上分别均匀地离散为N个点，得到如式(2)所示的样本集，该样本集表示各离散时刻ti的推力值f2m(tN)与设计变量Xi之间的对应关系：迭代初始时，令迭代次数k＝0，指定搜索精度eps和最大搜索步数Kmax，4)根据式(2)构造每个离散时刻推力的代理模型si(X)，得到N个代理模型，同时si(X)满足式(3)：F(ti)＝si(X)    (3)其中，F(ti)为离散时刻ti对应的推力，其中ti中的i满足1在所得设计变量最优解Xk+2m处运行性能仿真模型，得到对应的最优推力曲线，并将最优推力曲线在每个工作时间点上离散为N个点，添加到矩阵(2)中，则公式(2)中的样本点个数从2m个变为2m+k个；6)收敛判定：按照公式(5)进行迭代计算得到任意两不同的离散时刻分别对应的推力曲线的均方偏差：其中，F2m+k‑1(ti)为在ti离散时刻第2m+k‑1个离散点处的推力曲线，F2m+k(ti)为在ti离散时刻第2m+k个离散点处的推力曲线，F0(ti)为设计指标要求的推力曲线，N为离散点个数；若error(k)＜eps或k＝Kmax，则停止搜索，输出设计变量最优解Xk+2m及其对应的最优推力曲线f2m+k(t)，[X2m+k,f2m+k(t)]，否则，转步骤4)直至满足该条件时迭代停止。