一种液体冲压发动机风斗油路燃油流量闭环控制方法与流程

本发明属于发动机技术领域，具体涉及一种液体冲压发动机风斗油路燃油流量闭环控制方法。

背景技术：

液体冲压发动机风斗油路采用伺服阀进行燃油流量调节，传统的控制方法是试验标定伺服阀驱动电流与燃油流量和伺服阀前后压差之间的关系，通过实时采集伺服阀前后压差计算控制伺服阀电流来实现燃油流量的开环控制。虽然这种方法实现简单、易于工程实现。但是由于伺服阀的流量和其前后压差存在正反馈关系，开环控制容易出现不稳定现象，无法适应燃油流量的大范围调节。

技术实现要素：

(一)要解决的技术问题

本发明要解决的技术问题是：如何设计一种液体冲压发动机风斗油路燃油流量闭环控制方法，能够解决开环控制中出现的振荡问题，提高系统控制的稳定性。

(二)技术方案

为了解决上述技术问题，本发明提供了一种液体冲压发动机风斗油路燃油流量闭环控制方法，包括以下步骤：

(s1)根据当前周期输入的伺服阀5驱动电流is(k)和采集到的伺服阀5前油压信号pep、伺服阀5后油压信号pev，利用公式(1)计算伺服阀反馈燃油流量qmf_fk(k)：

qmf_fk(k)＝p00+p10*△p(k)+p01*is(k)+p11*△p(k)*is(k)+p02*is(k)²+p12*△p(k)*is(k)²+p03*is(k)³(1)

其中，△p(k)＝pep(k)-pev(k)(2)

p00、p10、p01、p11、p02、p12、p03为根据伺服阀5试验数据得到的拟合参数；

(s2)根据当前周期风斗燃油流量指令值qmf(k)计算当前时刻指令滤波后的风斗燃油流量指令值qmf_lv(k)，以及下一周期的滤波器状态量x(k+1)，具体计算方法如下：

qmf_lv(k)＝(0.4066*x(k)+0.09091*(qmf(k)-y_0))+y_0(3)

x(k+1)＝0.8182*x(k)+(qmf(k)-y_0))*0.4066(4)

其中，x(k)为滤波器状态，初始值x(0)＝0，y_0为燃油流量指令滤波基值，取值为第一个周期计算的伺服阀反馈燃油流量值qmf_fk(0)；

(s3)根据回路成型的方法求解流量闭环控制器2的参数，回路成型控制器参数确定方法即为通过控制指标要求确定系统开环对数频率特性模型参数，从而确定流量闭环控制器2参数的方法；

根据选定的相角裕度pm和系统调节时间要求值ts利用公式(5)、(6)、(7)计算系统谐振峰值mr、中频段宽度h、剪切频率ωc。

选取开环增益kv和高频道系数δ，利用公式(8)、(9)和(10)确定中频段系数γ、α、低频段系数β，确定开环对数频率特性模型中的低频段转折频率中频段衔接频率中频段转折频率γωc和高频道转折频率γδωc；然后根据公式(11)确定系统开环传递函数g(s)，根据公式(12)确定流量闭环闭环控制器2传递函数gc(s)；

其中，s为拉普拉斯算子，gg(s)为被控对象线性化传递函数模型，所述被控对象为伺服阀5；

根据总体提出系统调节时间ts，相角欲度pm、式(13)的被控对象线性化传递函数gg(s)、选取的相角裕度、开环增益、高频道系数，确定闭环控制器传递函数gc(s)：

(s4)根据步骤s1计算得到的伺服阀反馈燃油流量qmf_fk(k)、步骤s2计算的流量闭环控制器2传递函数及流量闭环控制器2的输入qmf_lv(k)计算流量闭环控制器2的输出；

首先根据qmf_lv(k)，qmf_fk(k)，求解归一化后的燃油流量误差量e(k)。

e(k)＝(qmf_lv(k)-qmf_fk(k))/qmf_fd0(15)

其中，qmf_fd0为设计点风斗油路燃油流量；

将公式(14)表示的流量闭环控制器2传递函数转换为离散状态空间形式，见公式(16)、(17)：

x(k+1)＝ax(k)+be(k)(16)

y(k)＝cx(k)+de(k)(17)

其中a、b、c、d为系数；

根据离散化后的流量闭环控制器2传递函数解算当前时刻流量闭环控制器2输出的伺服阀5驱动电流is(k),见公式(18)：

is(k)＝y(k)*i0+is(0)(18)

其中i0为流量闭环控制器2设计点伺服阀驱动电流，本实施例中，i0＝90。is(0)为流量闭环控制器2输出基值；

(s5)把驱动电流is转换为电压vs，并通过da电路3输出给伺服阀驱动电路4；

(s6)伺服阀驱动电路4输出响应的驱动电流i给伺服阀5。

优选地，利用阀前压力传感器6采集伺服阀5前油压信号pep。

优选地，利用阀后压力传感器7采集伺服阀5后油压信号pev。

优选地，拟合参数p00＝-0.002654，p10＝0.01263，p01＝0.0007698，p11＝0.0003604，p02＝3.08e-06，p12＝5.14e-07，p03_1＝-1.57e-08。

优选地，总体提出系统调节时间ts＝2s，相角欲度pm不小于40°，被控对象线性化传递函数为式(13)。

优选地，选取相角裕度为60°，开环增益kv＝10和高频道系数δ＝2。

优选地，qmf_fd0＝0.1。

优选地，

c＝[1.928-1.447-0.5678]，d＝0.2673。

优选地，伺服阀驱动电路与电压转换关系为vs＝20is。

(三)有益效果

本发明给出了一种液体冲压发动机风斗油路燃油流量闭环控制方法，能够解决开环控制中出现的振荡问题，提高了系统流量控制的稳定性；同时该方法可根据不同发动机需求，通过给定的控制指标直接计算出控制器参数，也可根据伺服阀特性更换伺服阀，节省了闭环控制算法参数整定的时间。

附图说明

图1为本发明的闭环控制方法原理图；

图2为系统开环对数频率特性模型曲线图；

图3是本发明的方法示意图。

具体实施方式

为使本发明的目的、内容、和优点更加清楚，下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。

图1为一种液体冲压发动机风斗油路燃油流量闭环控制方法原理图，其对应的控制系统组成包括指令滤波器1、控制器2、da电路3、驱动电路4、伺服阀5、阀前压力传感器6、阀后油压传感器7、燃油流量解算装置8。本发明的具体实施步骤如下：

(s1)燃油流量解算装置8根据当前周期输入的伺服阀5驱动电流is(k)和利用阀前压力传感器6采集到的伺服阀5前油压信号pep、利用阀后压力传感器7采集到的伺服阀5后油压信号pev，利用公式(1)计算伺服阀反馈燃油流量qmf_fk(k)：

qmf_fk(k)＝p00+p10*△p(k)+p01*is(k)+p11*△p(k)*is(k)+p02*is(k)²+p12*△p(k)*is(k)²+p03*is(k)³(1)

其中，△p(k)＝pep(k)-pev(k)(2)

p00、p10、p01、p11、p02、p12、p03为根据伺服阀5试验数据得到的拟合参数，例如某伺服阀p00＝-0.002654，p10＝0.01263，p01＝0.0007698，p11＝0.0003604，p02＝3.08e-06，p12＝5.14e-07，p03_1＝-1.57e-08；

(s2)指令滤波器1根据当前周期风斗燃油流量指令值qmf(k)计算当前时刻指令滤波后的风斗燃油流量指令值qmf_lv(k)，以及下一周期的滤波器状态量x(k+1)，具体计算方法如下：

qmf_lv(k)＝(0.4066*x(k)+0.09091*(qmf(k)-y_0))+y_0(3)

x(k+1)＝0.8182*x(k)+(qmf(k)-y_0))*0.4066(4)

其中，x(k)为滤波器状态，初始值x(0)＝0，y_0为燃油流量指令滤波基值，取值为第一个周期计算的伺服阀反馈燃油流量值qmf_fk(0)；

(s3)根据总体提出的控制指标要求，根据回路成型的方法求解流量闭环控制器2的参数。回路成型控制器参数确定方法即为通过控制指标要求确定系统开环对数频率特性模型(见图2)参数，从而确定流量闭环控制器2参数的方法。

图2为系统开环对数频率特性模型曲线图；

根据选定的相角裕度pm和系统调节时间要求值ts利用公式(5)、(6)、(7)计算系统谐振峰值mr、中频段宽度h、剪切频率ωc。

选取开环增益kv和高频道系数δ，利用公式(8)、(9)和(10)确定中频段系数γ、α、低频段系数β，确定图2中的低频段转折频率中频段衔接频率中频段转折频率γωc和高频道转折频率γδωc；然后根据公式(11)确定系统开环传递函数g(s)，根据公式(12)确定流量闭环闭环控制器2传递函数gc(s)。

其中，s为拉普拉斯算子，gg(s)为被控对象(伺服阀5)线性化传递函数模型。

例如总体提出系统调节时间ts＝2s，相角欲度pm不小于40°，被控对象线性化传递函数为式(13)；选取相角裕度为60°，开环增益kv＝10和高频道系数δ＝2，可确定闭环控制器传递函数，见公式(14)。

(s4)根据步骤s1计算得到的伺服阀反馈燃油流量qmf_fk(k)、步骤s2计算的流量闭环控制器2传递函数及流量闭环控制器2的输入qmf_lv(k)计算流量闭环控制器2的输出。

首先根据qmf_lv(k)，qmf_fk(k)，求解归一化后的燃油流量误差量e(k)。

e(k)＝(qmf_lv(k)-qmf_fk(k))/qmf_fd0(15)

其中，qmf_fd0为设计点风斗油路燃油流量，本实施例中qmf_fd0＝0.1。

将公式(14)表示的流量闭环控制器2传递函数转换为离散状态空间形式，见公式(16)、(17)：

x(k+1)＝ax(k)+be(k)(16)

y(k)＝cx(k)+de(k)(17)

c＝[1.928-1.447-0.5678]，d＝0.2673

根据离散化后的流量闭环控制器2传递函数解算当前时刻流量闭环控制器2输出的伺服阀5驱动电流is(k),见公式(18)。

is(k)＝y(k)*i0+is(0)(18)

其中i0为流量闭环控制器2设计点伺服阀驱动电流，本实施例中，i0＝90。is(0)为流量闭环控制器2输出基值(根据实际起控点选取)，本实施例中，取is(0)＝100。

(s5)把驱动电流is转换为电压vs，并通过da电路3输出给伺服阀驱动电路4。

例如，某型发动机伺服阀驱动电路与电压转换关系为vs＝20is。

(s6)伺服阀驱动电路4输出响应的驱动电流i给伺服阀5。

以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。