CN107239634A - 一种航空发动机过渡过程建模方法 - Google Patents
一种航空发动机过渡过程建模方法 Download PDFInfo
- Publication number
- CN107239634A CN107239634A CN201710482860.7A CN201710482860A CN107239634A CN 107239634 A CN107239634 A CN 107239634A CN 201710482860 A CN201710482860 A CN 201710482860A CN 107239634 A CN107239634 A CN 107239634A
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtd
- msub
- mtr
- rho
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明提供一种航空发动机过渡过程建模方法,属于航空发动机过渡态建模与仿真领域。本发明通过辨识航空发动机不同稳定状态下的动态模型,与其相应稳态模型相结合,利用多项式拟合及最小二乘支持向量机的方法将航空发动机过渡过程表征为一种动态变参数模型。本发明对多离散稳态点近似过渡过程加以改进,动态特性更好,并改进了现有方法利用输出累加进行过渡态建模存在误差累加的情况;同时本方法针对单输入多输出进行建模,普适性更广,适当的调整即可推广到多输入多输出航空发动机过渡过程建模。
Description
技术领域
本发明涉及航空发动机过渡态建模与仿真领域,具体说,涉及一种航空发动机过渡态控制显示数学模型建立方法。
背景技术
在航空发动机控制领域,航空发动机过渡态建模一直是难以解决的问题,针对过渡态进行数学模型的建立,一方面可以使航空发动机过渡态使用的控制方法更具有针对性,另一方面可以有效的缩减航空发动机控制方面的研究时间,降低研究成本,因此本发明针对航空发动机进行过渡态建模就有很大的意义。
文献表明,现有的航空发动机建模多针对稳态工作点建模,如只能针对白箱模型进行设计的稳态点抽功法,该类方法建模以非线性计算程序为基础,通过抽取航空发动机高压转子和低压转子轴上功率和依次改变发动机的输入量,利用平衡计算获取发动机各部件的参数变化量,并根据获取的参数变化量的偏导数计算航空发动机稳态工作点模型。而对于航空发动机的过渡过程,目前通用的建模方法仍然基于稳态工作点模型,即对特征稳态工作点采用插值算法,近似过渡过程。由于特征稳态工作点较少,插值误差较大,导致模型动态响应误差较大。此外,该建模方法利用输出累加进行过渡态建模,即利用稳态点模型的小偏离输出量和上一时刻的输出的加和作为此稳态点处的模型输出量,但此种方法由于输出误差不断累加,模型精确度也在不断下降。因此,为了克服现有的控制用航空发动机过渡态建模方式的不足,即多点近似的效果不能很好的反映发动机过渡态动态性能,且输出累加进行过渡态建模降低模型精度的特点,本发明提出一种控制用航空发动机过渡态等效建模方法建模方法。此种方法对多离散稳态点近似过渡过程加以改进,动态特性更好;并针对现有方法利用输出累加进行过渡态建模存在误差累加的情况进行改进。同时本方法针对单输入多输出进行建模,普适性更广,适当的调整即可推广到多输入多输出航空发动机慢车到最大状态过程建模。
发明内容
针对现有技术中航空发动机控制器只能针对某稳态附近的数学模型进行设计,过渡过程无法用模型进行表述的问题,本发明提供一种基于模型辨识的航空发动机过渡态建模方法。
本发明的技术方案为:
一种航空发动机过渡过程建模方法,包括以下步骤:
步骤1:产生一组正弦激励信号序列{Δu}
正弦激励信号序列{Δu}由1-80Hz内不同幅值与相位相互叠加的正弦信号组成,幅值在[-1,1]之内。
步骤2:得到给定油量Wfsi下的高压涡轮相对换算转速序列{yi1}和低压涡轮相对换算转速序列{yi2},i=1,2,3…n。
2.1)向航空发动机输入给定的油量Wfs1,待航空发动机高压涡轮相对换算转速达到相应稳定状态ρ1后,将正弦激励信号序列{Δu}输入到航空发动机中,分别采集航空发动机输出高压涡轮相对换算转速响应序列{y11}和低压涡轮相对换算转速响应序列{y12}。
2.2)向航空发动机输入给定的油量Wfs2,待航空发动机高压涡轮相对换算转速达到相应稳定状态ρ2后,将正弦激励信号序列{Δu}输入到航空发动机中,分别采集航空发动机输出高压涡轮相对换算转速响应序列{y21}和低压涡轮相对换算转速响应序列{y22}。
2.3)重复上述过程n次,得到n个相应稳定状态ρi,分别采集得到给定油量Wfsi下的高压涡轮相对换算转速序列{yi1}和低压涡轮相对换算转速序列{yi2},i=1,2,3…n。
步骤3:将每个采集的高压涡轮相对换算转速序列{yi1}和低压涡轮相对换算转速序列{yi2}进行去均值化处理,得到相应序列{Δyi1}{Δyi2};将由正弦激励信号序列{Δu}作为输入数据,{Δyi1}、{Δyi2}作为输出数据组成的数据样本采用子空间迭代算法对系统的结构参数进行辨识,得到系统的多个动态模型Mi表达式为:
其中:x为状态变量不代表物理含义;ΔYi表示每个采集的高压涡轮相对换算转速序列{yi1}和低压涡轮相对换算转速序列{yi2}均值化处理后的相应序列{Δyi1}{Δyi2}的列组合,表示为:ΔYi=[Δyi1 Δyi2]T;Ai、Bi、Ci表示根据子空间迭代算法对系统的结构参数进行辨识后得到的系数矩阵,其中i=1,2,3…n。
步骤4:根据步骤3得到的n个动态模型Mi能够得到系数矩阵,其中i=1,2,3…n,系数矩阵包括{a111,a211,a311…an11},{a112,a212,a312…an12},{a121,a221,a321…an21},{a122,a222,a322…an22},{b11,b21,b31…bn1},{b12,b22,b32…bn2},{c111,c211,c311…cn11},{c112,c212,c312…cn12},{c121,c221,c321…cn21},{c122,c222,c322…cn22};
将系数矩阵拟合为高压涡轮相对换算转速稳定状态ρ的函数,分别表示为:a11(ρ),a12(ρ),a21(ρ),a22(ρ);b1(ρ),b2(ρ);c11(ρ),c12(ρ),c21(ρ),c22(ρ),拟合方法为:
以步骤2.3)得到的n个相应稳定状态ρi为横坐标,以系数矩阵为纵坐标,得到散点图,根据散点图确定高压涡轮相对换算转速稳定状态ρ为平滑函数或非平滑函数;
若系数矩阵为相应高压涡轮相对换算转速稳定状态ρ的平滑函数,则采用高次多项式拟合法或高次傅里叶函数拟合法进行拟合;
若系数矩阵为相应高压涡轮相对换算转速稳定状态ρ的非平滑函数,则采用最小二乘支持向量机法LS-SVM行拟合,其中核函数选为径向基函数宽度为σ。
步骤5:以步骤2.3)得到的n个相应稳定状态ρi即{ρ1,ρ2,ρ3…ρn},作为输入,以给定油量{Wfs1,Wfs2Wfs3…Wfsn}作为输出,采用高次多项式拟合法拟合函数φ(ρ)。
定义{u1,u2,u3…un}={Wfs1,Wfs2Wfs3…Wfsn},将{u1,u2,u3…un}作为输入,步骤2.3)得到的n个相应稳定状态ρi即{ρ1,ρ2,ρ3…ρn}作为输出,采用高次多项式拟合法拟合函数ψ(u)。
步骤6:根据步骤4、步骤5的结果,得到航空发动机过渡过程模型Z为:
其中,w为状态变量不具有实际的物理含义,u为航空发动机过渡过程输入的油量信号,y为航空发动机过渡过程相应输出的高压涡轮相对换算转速和低压涡轮相对换算转速。
本发明还可以将步骤2中输出的输出值调整为压气机出口总压P3、涡轮出口总压P5、涡轮出口总温T5中的任意两种,输入值不变,对航空发动机建模。此时,步骤4中模型Mi调整为其中,ΔYi=[Δyi1 Δyi2]T, Di=[d]i;Di表示系数矩阵,根据子空间迭代算法对系统的结构参数进行辨识后得到,最终得到的航空发动机过渡态模型L为:
其中,d(ρ)根据步骤4的方法拟合得到。
本发明的有益效果为:通过本发明的方法设计的航空发动机可以有效地对航空发动机的过渡过程进行描述,相比于以往的建模方法,精度更高,动态性能更好,实现了对现有方法利用输出累加进行过渡态建模存在误差累加情况的改进。同时本方法针对单输入多输出进行建模,普适性更广,适当的调整即可推广到多输入多输出航空发动机慢车到最大状态过程建模,为航空发动机过度过程控制器设计提供了更具有针对性的模型。
附图说明
图1为航空发动机过渡过程建模流程图;
图2为频率为1-80Hz由不同幅值与相位相互叠加的幅值在[-1,1]之内的正弦激励信号Δu;
图3为系数矩阵a11(ρ)最小二乘支持向量机拟合曲线;
图4为系数矩阵a11(ρ)最小二乘支持向量机拟合曲线的误差值;
图5为系数矩阵b1(ρ)多项式拟合曲线;
图6为测试油量输入信号;
图7为本发明建模仿真验证对比图。
具体实施方式
下面结合附图对本发明作进一步说明,本发明的研究对象为某小型无人机用涡扇发动机,其过渡过程模型建立如图1流程图所示,详细设计步骤如下:
步骤1:产生一组由1-80Hz内不同幅值与相位相互叠加的正弦信号组成,幅值在[-1,1]之内的正弦激励信号序列{Δu},其中样例中采用的信号如图2所示。
步骤2得到给定油量Wfsi下的高压涡轮相对换算转速序列{yi1}和低压涡轮相对换算转速序列{yi2},i=1,2,3…n
2.1)向航空发动机输入给定的油量Wfs1,待航空发动机高压涡轮相对换算转速达到相应稳定状态ρ1后,将正弦激励信号序列{Δu}输入到航空发动机中,分别采集航空发动机输出高压涡轮相对换算转速响应序列{y11}和低压涡轮相对换算转速响应序列{y12}。
2.2)向航空发动机输入给定的油量Wfs2,待航空发动机高压涡轮相对换算转速达到相应稳定状态ρ2后,将正弦激励信号序列{Δu}输入到航空发动机中,分别采集航空发动机输出高压涡轮相对换算转速响应序列{y21}和低压涡轮相对换算转速响应序列{y22}。
2.3)重复上述过程n次,得到n个相应稳定状态ρi,分别采集得到给定油量Wfsi下的高压涡轮相对换算转速序列{yi1}和低压涡轮相对换算转速序列{yi2},i=1,2,3…n。
步骤3:将每个采集的高压涡轮相对换算转速序列{yi1}和低压涡轮相对换算转速序列{yi2}进行去均值化处理,得到相应序列{Δyi1}{Δyi2}。将由正弦激励信号序列{Δu}作为输入数据,{Δyi1}、{Δyi2}作为输出数据组成的数据样本采用子空间迭代算法对系统的结构参数进行辨识,得到系统的动态模型Mi,表达式为:
其中:
x为状态变量不代表物理含义;ΔYi表示每个采集的高压涡轮相对换算转速序列{yi1}和低压涡轮相对换算转速序列{yi2}均值化处理后的相应序列{Δyi1}{Δyi2}的列组合,表示为:ΔYi=[Δyi1 Δyi2]T;Ai、Bi、Ci表示系数矩阵根据子空间迭代算法对系统的结构参数进行辨识后得到,其中i=1,2,3…n。
步骤4:根据步骤3得到的n个动态模型Mi能够得到系数矩阵,其中i=1,2,3…n,系数矩阵包括{a111,a211,a311…an11},{a112,a212,a312…an12},{a121,a221,a321…an21},{a122,a222,a322…an22},{b11,b21,b31…bn1},{b12,b22,b32…bn2},{c111,c211,c311…cn11},{c112,c212,c312…cn12},{c121,c221,c321…cn21},{c122,c222,c322…cn22};
将系数矩阵拟合为高压涡轮相对换算转速稳定状态ρ的函数,分别表示为:a11(ρ),a12(ρ),a21(ρ),a22(ρ);b1(ρ),b2(ρ);c11(ρ),c12(ρ),c21(ρ),c22(ρ),拟合方法为:
以步骤2.3)得到的n个相应稳定状态ρi为横坐标,以系数矩阵为纵坐标,得到散点图,根据散点图确定若系数矩阵相应高压涡轮相对换算转速稳定状态ρ的平滑函数或非平滑函数;
若系数矩阵为相应高压涡轮相对换算转速稳定状态ρ的平滑函数,则采用高次多项式拟合法或高次傅里叶函数拟合法进行拟合。如图5所示,在对b1(ρ)进行拟合时,由于{b11,b21,b31…bn1}是高压涡轮相对换算转速稳定状态ρ的平滑函数所以采用高次多项式的方式进行拟合。
若系数矩阵为相应高压涡轮相对换算转速稳定状态ρ的非平滑函数,则采用最小二乘支持向量机法LS-SVM进行拟合,核函数选用径向基函数宽度为σ。如图3所示:在对a11(ρ)进行拟合时,由于{a111,a211,a311…an11}序列为函数高压涡轮相对换算转速稳定状态ρ的非平滑函数,因此用最小二乘支持向量机法LS-SVM进行拟合,核函数选用径向基函数宽度为σ=0.6,得到a11(ρ)的表达式,在径向基函数中zi为支持向量机,通过数据训练得到,z即为a11(ρ)中的ρ,即mi,n通过LS-SVM方法计算得到。其拟合精度如图4所示,误差接近于0,符合误差要求。
步骤5:以步骤2.3)得到的n个相应稳定状态ρi即{ρ1,ρ2,ρ3…ρn},作为输入,以给定油量{Wfs1,Wfs2Wfs3…Wfsn}作为输出,采用高次多项式拟合法拟合函数φ(ρ)。
定义{u1,u2,u3…un}={Wfs1,Wfs2Wfs3…Wfsn},将{u1,u2,u3…un}作为输入,步骤2.3)得到的n个相应稳定状态ρi即{ρ1,ρ2,ρ3…ρn}作为输出,采用高次多项式拟合法拟合函数ψ(u)。
步骤6:航空发动机过渡态模型Z为:
步骤7:对模型Z进行验证,向模型输入任意斜坡叠加信号如图6所示,验证模型输出如图7所示,实线为Z模型输出曲线,虚线为航空发动机输出曲线,其动态特性近似相同。稳态特性基本一致,从数据分析瞬时模型输出误差为12.3%,平均动态模型输出误差为7%。
Claims (2)
1.一种航空发动机过渡过程建模方法,其特征在于以下步骤:
步骤1:产生一组正弦激励信号序列{Δu}
正弦激励信号序列{Δu}由1-80Hz内不同幅值与相位相互叠加的正弦信号组成,幅值在[-1,1]之内;
步骤2:得到给定油量Wfsi下的高压涡轮相对换算转速序列{yi1}和低压涡轮相对换算转速序列{yi2},i=1,2,3…n;
2.1)向航空发动机输入给定的油量Wfs1,待航空发动机高压涡轮相对换算转速达到相应稳定状态ρ1后,将正弦激励信号序列{Δu}输入到航空发动机中,分别采集航空发动机输出高压涡轮相对换算转速响应序列{y11}和低压涡轮相对换算转速响应序列{y12};
2.2)向航空发动机输入给定的油量Wfs2,待航空发动机高压涡轮相对换算转速达到相应稳定状态ρ2后,将正弦激励信号序列{Δu}输入到航空发动机中,分别采集航空发动机输出高压涡轮相对换算转速响应序列{y21}和低压涡轮相对换算转速响应序列{y22};
2.3)重复上述过程n次,得到n个相应稳定状态ρi,分别采集得到给定油量Wfsi下的高压涡轮相对换算转速序列{yi1}和低压涡轮相对换算转速序列{yi2},i=1,2,3…n;
步骤3:将每个采集的高压涡轮相对换算转速序列{yi1}和低压涡轮相对换算转速序列{yi2}进行去均值化处理,得到相应序列{Δyi1}{Δyi2};将由正弦激励信号序列{Δu}作为输入数据,{Δyi1}、{Δyi2}作为输出数据组成的数据样本采用子空间迭代算法对系统的结构参数进行辨识,得到系统的多个动态模型Mi表达式为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mi>x</mi>
<mo>+</mo>
<msub>
<mi>B</mi>
<mi>i</mi>
</msub>
<mi>&Delta;</mi>
<mi>u</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;Y</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>C</mi>
<mi>i</mi>
</msub>
<mi>x</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:x为状态变量不代表物理含义;ΔYi表示每个采集的高压涡轮相对换算转速序列{yi1}和低压涡轮相对换算转速序列{yi2}均值化处理后的相应序列{Δyi1}{Δyi2}的列组合,表示为:ΔYi=[Δyi1 Δyi2]T;Ai、Bi、Ci表示根据子空间迭代算法对系统的结构参数进行辨识后得到的系数矩阵,其中i=1,2,3…n;
步骤4:根据步骤3得到的n个动态模型Mi能够得到系数矩阵,其中i=1,2,3…n,系数矩阵包括{a111,a211,a311…an11},{a112,a212,a312…an12},{a121,a221,a321…an21},{a122,a222,a322…an22},{b11,b21,b31…bn1},{b12,b22,b32…bn2},{c111,c211,c311…cn11},{c112,c212,c312…cn12},{c121,c221,c321…cn21},{c122,c222,c322…cn22};
将系数矩阵拟合为高压涡轮相对换算转速稳定状态ρ的函数,分别表示为:a11(ρ),a12(ρ),a21(ρ),a22(ρ);b1(ρ),b2(ρ);c11(ρ),c12(ρ),c21(ρ),c22(ρ),拟合方法为:
以步骤2.3)得到的n个相应稳定状态ρi为横坐标,以系数矩阵为纵坐标,得到散点图,根据散点图确定高压涡轮相对换算转速稳定状态ρ为平滑函数或非平滑函数;
若系数矩阵为相应高压涡轮相对换算转速稳定状态ρ的平滑函数,则采用高次多项式拟合法或高次傅里叶函数拟合法进行拟合;
若系数矩阵为相应高压涡轮相对换算转速稳定状态ρ的非平滑函数,则采用最小二乘支持向量机法LS-SVM行拟合,其中核函数选为径向基函数宽度为σ;
步骤5:以步骤2.3)得到的n个相应稳定状态ρi即{ρ1,ρ2,ρ3…ρn},作为输入,以给定油量{Wfs1,Wfs2Wfs3…Wfsn}作为输出,采用高次多项式拟合法拟合函数φ(ρ);
定义{u1,u2,u3…un}={Wfs1,Wfs2Wfs3…Wfsn},将{u1,u2,u3…un}作为输入,步骤2.3)得到的n个相应稳定状态ρi即{ρ1,ρ2,ρ3…ρn}作为输出,采用高次多项式拟合法拟合函数ψ(u);
步骤6:根据步骤4、步骤5的结果,得到航空发动机过渡过程模型Z为:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>w</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>11</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>12</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>21</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>22</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>w</mi>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mo>(</mo>
<mrow>
<mi>u</mi>
<mo>-</mo>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>y</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>11</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>12</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>21</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>22</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>w</mi>
<mo>+</mo>
<mi>&psi;</mi>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,w为状态变量不具有实际的物理含义,u为航空发动机过渡过程输入的油量信号,y为航空发动机过渡过程相应输出的高压涡轮相对换算转速和低压涡轮相对换算转速。
2.根据权利要求1所述的一种航空发动机过渡过程建模方法,其特征在于,还可以将步骤2中输出的输出值调整为压气机出口总压P3、涡轮出口总压P5、涡轮出口总温T5中的任意两种,输入值不变,对航空发动机建模;此时,步骤4中模型Mi调整为其中,ΔYi=[Δyi1 Δyi2]T,Di=[d]i;Di表示系数矩阵,根据子空间迭代算法对系统的结构参数进行辨识后得到,最终得到的航空发动机过渡态模型L为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>w</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>11</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>12</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>21</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>22</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>w</mi>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mo>(</mo>
<mrow>
<mi>u</mi>
<mo>-</mo>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>y</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>11</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>12</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>21</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>c</mi>
<mn>22</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>w</mi>
<mo>+</mo>
<mi>&psi;</mi>
<mrow>
<mo>(</mo>
<mi>u</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>d</mi>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
<mo>(</mo>
<mrow>
<mi>u</mi>
<mo>-</mo>
<mi>&phi;</mi>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,d(ρ)根据步骤4的方法拟合得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710482860.7A CN107239634B (zh) | 2017-06-26 | 2017-06-26 | 一种航空发动机过渡过程建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710482860.7A CN107239634B (zh) | 2017-06-26 | 2017-06-26 | 一种航空发动机过渡过程建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107239634A true CN107239634A (zh) | 2017-10-10 |
CN107239634B CN107239634B (zh) | 2019-10-11 |
Family
ID=59986689
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710482860.7A Active CN107239634B (zh) | 2017-06-26 | 2017-06-26 | 一种航空发动机过渡过程建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107239634B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108663948A (zh) * | 2018-05-17 | 2018-10-16 | 西北工业大学 | 一种航空发动机控制系统数值仿真分析平台的设计方法 |
CN108762089A (zh) * | 2018-06-15 | 2018-11-06 | 大连理工大学 | 一种基于模型预测的航空发动机在线优化及多变量控制设计方法 |
CN109611217A (zh) * | 2018-11-07 | 2019-04-12 | 大连理工大学 | 一种航空发动机过渡态控制规律优化的设计方法 |
CN109829238A (zh) * | 2019-02-02 | 2019-05-31 | 哈尔滨工业大学 | 基于折合平衡流形展开模型的航空发动机系统辨识方法 |
WO2019144386A1 (zh) * | 2018-01-26 | 2019-08-01 | 大连理工大学 | 一种航空发动机过渡态关键性能参数预测方法 |
CN110348079A (zh) * | 2019-06-25 | 2019-10-18 | 南京航空航天大学 | 基于机器学习算法的涡轴发动机稳态开环控制设计方法 |
WO2020118513A1 (zh) * | 2018-12-11 | 2020-06-18 | 大连理工大学 | 一种基于自抗扰理论的涡扇发动机稳态过渡态多变量控制方法 |
CN111679574A (zh) * | 2020-05-13 | 2020-09-18 | 大连理工大学 | 一种基于大规模全局优化技术的变循环发动机过渡态优化方法 |
CN113065206A (zh) * | 2021-03-24 | 2021-07-02 | 北京航空航天大学 | 一种过渡态控制方法、装置、电子设备及存储介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020016665A1 (en) * | 1998-10-22 | 2002-02-07 | Ulyanov Sergei V. | System for intelligent control of an engine based on soft computing |
CN103778295A (zh) * | 2014-01-26 | 2014-05-07 | 南京航空航天大学 | 多失效模式下多模型集成航空发动机运行可靠性评估方法 |
CN103838920A (zh) * | 2014-02-13 | 2014-06-04 | 南京航空航天大学 | 一种涡轴发动机动态实时模型建模方法及系统 |
CN104317997A (zh) * | 2014-10-17 | 2015-01-28 | 北京航空航天大学 | 一种高负荷风扇/压气机端壁造型优化设计方法 |
CN105512429A (zh) * | 2015-12-31 | 2016-04-20 | 中国航空工业集团公司沈阳发动机设计研究所 | 一种三轴燃气轮机的总体方案计算方法 |
CN106647253A (zh) * | 2016-09-28 | 2017-05-10 | 南京航空航天大学 | 航空发动机分布式控制系统多性能鲁棒跟踪控制方法 |
-
2017
- 2017-06-26 CN CN201710482860.7A patent/CN107239634B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020016665A1 (en) * | 1998-10-22 | 2002-02-07 | Ulyanov Sergei V. | System for intelligent control of an engine based on soft computing |
CN103778295A (zh) * | 2014-01-26 | 2014-05-07 | 南京航空航天大学 | 多失效模式下多模型集成航空发动机运行可靠性评估方法 |
CN103838920A (zh) * | 2014-02-13 | 2014-06-04 | 南京航空航天大学 | 一种涡轴发动机动态实时模型建模方法及系统 |
CN104317997A (zh) * | 2014-10-17 | 2015-01-28 | 北京航空航天大学 | 一种高负荷风扇/压气机端壁造型优化设计方法 |
CN105512429A (zh) * | 2015-12-31 | 2016-04-20 | 中国航空工业集团公司沈阳发动机设计研究所 | 一种三轴燃气轮机的总体方案计算方法 |
CN106647253A (zh) * | 2016-09-28 | 2017-05-10 | 南京航空航天大学 | 航空发动机分布式控制系统多性能鲁棒跟踪控制方法 |
Non-Patent Citations (1)
Title |
---|
闫锋: "一种航空发动机燃油流量基线的建模方法", 《计算机测量与控制》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019144386A1 (zh) * | 2018-01-26 | 2019-08-01 | 大连理工大学 | 一种航空发动机过渡态关键性能参数预测方法 |
US11124317B2 (en) | 2018-01-26 | 2021-09-21 | Dalian University Of Technology | Method for prediction of key performance parameters of aero-engine in transition condition |
CN108663948A (zh) * | 2018-05-17 | 2018-10-16 | 西北工业大学 | 一种航空发动机控制系统数值仿真分析平台的设计方法 |
CN108762089A (zh) * | 2018-06-15 | 2018-11-06 | 大连理工大学 | 一种基于模型预测的航空发动机在线优化及多变量控制设计方法 |
CN109611217B (zh) * | 2018-11-07 | 2020-12-11 | 大连理工大学 | 一种航空发动机过渡态控制规律优化的设计方法 |
CN109611217A (zh) * | 2018-11-07 | 2019-04-12 | 大连理工大学 | 一种航空发动机过渡态控制规律优化的设计方法 |
US11392094B2 (en) | 2018-12-11 | 2022-07-19 | Dalian University Of Technology | Steady state and transition state multi-variable control method of a turbofan engine based on an active disturbance rejection control theory (ADRC) |
WO2020118513A1 (zh) * | 2018-12-11 | 2020-06-18 | 大连理工大学 | 一种基于自抗扰理论的涡扇发动机稳态过渡态多变量控制方法 |
CN109829238A (zh) * | 2019-02-02 | 2019-05-31 | 哈尔滨工业大学 | 基于折合平衡流形展开模型的航空发动机系统辨识方法 |
CN109829238B (zh) * | 2019-02-02 | 2022-06-17 | 哈尔滨工业大学 | 基于折合平衡流形展开模型的航空发动机系统辨识方法 |
CN110348079B (zh) * | 2019-06-25 | 2021-06-15 | 南京航空航天大学 | 基于机器学习算法的涡轴发动机稳态开环控制设计方法 |
CN110348079A (zh) * | 2019-06-25 | 2019-10-18 | 南京航空航天大学 | 基于机器学习算法的涡轴发动机稳态开环控制设计方法 |
CN111679574B (zh) * | 2020-05-13 | 2021-05-07 | 大连理工大学 | 一种基于大规模全局优化技术的变循环发动机过渡态优化方法 |
CN111679574A (zh) * | 2020-05-13 | 2020-09-18 | 大连理工大学 | 一种基于大规模全局优化技术的变循环发动机过渡态优化方法 |
CN113065206A (zh) * | 2021-03-24 | 2021-07-02 | 北京航空航天大学 | 一种过渡态控制方法、装置、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN107239634B (zh) | 2019-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107239634A (zh) | 一种航空发动机过渡过程建模方法 | |
CN103306822B (zh) | 一种基于喘振裕度估计模型的航空涡扇发动机控制方法 | |
CN106874569B (zh) | 一种变几何分轴式燃气轮机实时仿真建模方法 | |
CN105209984B (zh) | 用于确定技术系统输出值的模型的方法 | |
CN103838920B (zh) | 一种涡轴发动机动态实时模型建模方法及系统 | |
CN103942357B (zh) | 包线内全状态的涡扇发动机机载实时模型构建方法 | |
CN108108528B (zh) | 一种分轴式航改燃机动力涡轮一维匹配设计方法 | |
CN105700380B (zh) | 二次再热机组汽轮机调速系统仿真模型及其建模方法 | |
CN105275508B (zh) | 一种基于功率值计算的汽轮机流量曲线辨识及优化方法 | |
CN107391858A (zh) | 一种获取风洞模型静气动弹性变形影响量的方法 | |
KR101474102B1 (ko) | 임펠러 형상 최적설계방법 | |
CN110219736A (zh) | 基于非线性模型预测控制的航空发动机直接推力控制方法 | |
CN107002636A (zh) | 用于估计风速,包括计算针对叶片扭转调节的桨距角的方法 | |
CN105676640B (zh) | 基于贝塞尔曲线的涡扇发动机加速过程控制律设计方法 | |
CN110221537A (zh) | 控制方法、装置及航空发动机限制保护控制方法、装置 | |
Tsoutsanis et al. | Part-load performance of gas turbines: part I—a novel compressor map generation approach suitable for adaptive simulation | |
CN105389427A (zh) | 基于自适应粒子滤波的航空发动机气路部件故障检测方法 | |
CN108196443A (zh) | 变循环发动机的非线性预测控制设计方法 | |
CN101813920B (zh) | 电站汽轮机组温度传感器虚拟冗余方法 | |
CN110516395A (zh) | 一种基于非线性模型预测的航空发动机控制方法 | |
CN109871653A (zh) | 航空发动机数学模型部件特性修正方法 | |
CN111680823A (zh) | 一种风向信息预测方法及系统 | |
CN109829238B (zh) | 基于折合平衡流形展开模型的航空发动机系统辨识方法 | |
CN107301268A (zh) | 一种船舶燃气轮机可调静叶压气机偏转角度优化方法 | |
CN109635512A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CB03 | Change of inventor or designer information |
Inventor after: Sun Ximing Inventor after: Liu Tongjian Inventor after: Ma Yanhua Inventor after: Wang Rui Inventor before: Liu Tongjian Inventor before: Sun Ximing Inventor before: Ma Yanhua Inventor before: Wang Rui |
|
CB03 | Change of inventor or designer information |