文章快速检索  
  高级检索
民机配平空气系统小偏差化动态建模仿真方法
张絮涵, 吴成云     
中国商用飞机有限责任公司 上海飞机设计研究院, 上海 201210
摘要: 配平空气系统(TAS)作为民用飞机空调系统(ACS)中的子系统之一,其本质是舱室内温度调节的高温气体支路。该系统的动态特性较大地影响着下游乃至舱室内空气参数的变化,不利于下游参数的稳定。因此需要依据系统固有的动态特性,通过设计系统的控制规律,达到稳定下游空气参数的目的。随着系统方案及部件的更新换代,民机配平空气系统的压力控制方式逐步由机械气动式压力调节活门转变为电动式压力调节活门辅以压力传感器的方式。这种方式要求对配平空气系统的动态特性进行快速、准确的仿真,并能够依据仿真结果进行压力控制律设计。本文提出小偏差化的线性化处理方法作为动态特性建模仿真的理论推导依据,对采用电动式压力调节活门和压力传感器方案的典型民机配平空气系统建立高效的动态仿真模型。经校验,该方法结果准确。此外,还根据动态特性分析结果给出电动式压力调节活门的控制律设计实例,有效实现了良好的压力调节效果。根据设计流程提出了配平空气系统的压力控制需求,为民机空调系统的动态仿真建模方法提供指导思路。
关键词: 配平空气系统     民机     小偏差化     压力调节     控制律设计    
Small deviation linearization method for dynamic modeling and simulation of trim air system in civil aircraft
ZHANG Xuhan, WU Chengyun     
Shanghai Aircraft Design and Research Institute, Commercial Aircraft Corporation of China, Ltd., Shanghai 201210, China
Abstract: As one of the subsystems of Air Conditioning System (ACS) in civil aircraft, Trim Air System (TAS) is essentially the hot air branch to regulate cabin temperature. While the dynamic characteristics of TAS greatly affect the variation of downstream and cabin air parameters, it goes against with the stability of downstream parameters. Therefore, the system control law should be designed based on the inherent dynamic characteristics of the system, so as to stabilize the downstream parameters. With the updating of the system architecture and system components, the pressure is regulated by electrically operated Trim Air Pressure Regulating Valve (TAPRV) and Trim Air Pressure Sensor (TAPS) instead of the traditional mechanical pneumatic TAPRV. The dynamic characteristics of TAS is required to be simulated efficiently and accurately. In this paper, a small deviation linearization method is provided for theoretical derivation, and an efficient dynamic simulation model for typical electrically operated TAS is established. The method is validated to be accurate. Then the control law of electrically TAPRV is designed as an example, achieving good pressure regulating effect. Based on the design flow, this research provides TAS pressure control requirements, and provides guidance for the dynamic modeling and simulation of ACS.
Keywords: trim air system     civil aircraft     small deviation linearization     pressure regulation     control law design    

民用飞机空调系统的主要功能是从气源系统引气,经过空调组件将高温高压气体调节为具有合适流量、温度和压力的空调供气,并通过低压导管输送到舱内,满足座舱通风、增压和温度调节要求[1-4]。其中,舱室内温度控制方式是由来自制冷组件上游的高温气体与混合腔各送风总管内的冷空气混合实现。这里的高温气体支路被称为配平空气系统(Trim Air System, TAS),主要用于独立控制进入舱室的供气温度,同时控制通过制冷组件的流量,防止由于压力变化产生引射而影响系统性能[5-7]。系统一般由配平空气压力调节活门(Trim Air Pressure Regulating Valve, TAPRV)实现下游压力调节,由配平空气活门(Trim Air Valve, TAV)控制进入主分配管路的流量来控制进入舱室的气体温度[8]

传统配平空气系统在控制压力时多采用机械气动式TAPRV,其工作主要取决于控制气体在作动筒腔内的压力变化以及和作动筒腔内弹簧进行的压力比较结果[9]。随着民机的不断发展,新的压力控制方案开始采用电动式TAPRV活门阀瓣开度,并在TAPRV下游增设压力传感器来控制配平空气系统下游压力。该方案是具有压力信号反馈的闭环回路,其实施依靠控制器完成,即通过合理的控制律来控制驱动电机的输入信号,从而控制TAPRV开度。

该系统的动态特性较大地影响着下游乃至舱室内空气参数的变化,不利于下游参数的稳定[10-13]。因此需要依据系统固有的动态特性,通过设计系统的控制规律,达到稳定下游空气参数的目的[14-16]。这种方式要求对配平空气系统的动态特性进行快速、准确的仿真[17]。目前国内对该系统的动态特性仿真研究偏少。

本文提出小偏差化的线性化处理方法作为动态特性建模仿真的理论推导依据,搭建了典型配平空气系统至舱室区域的动态仿真模型。文中分析了系统压力参数和温度参数的固有动态特性,并提供了电动式TAPRV的控制律设计实例。

1 典型配平空气系统动态特性建模 1.1 模型简化

根据配平空气系统原理,可将系统简化为如图 1所示的流程。为聚焦压力值和温度值的动态特性,这里关注制冷组件上游的热空气来流和混合腔的冷空气来流汇集后输送至舱室的通路过程,暂不考虑空气分配及单向截止流动等问题。模型保留了与压力、温度相关的典型部件,即活门TAPRV和TAV、压力传感器(Trim Air Pressure Sensor, TAPS)、管路温度传感器(Duct Temperature Sensor, DTS)、舱室温度传感器(Cabin Temperature Sensor, CTS)、热空气与冷空气汇流点以及舱室空间。

图 1 典型TAS建模简化系统图 Fig. 1 Diagram of typical TAS simulating simplification

图 1中,W1为TAPRV入口空气流量,W2为TAV出口空气流量,W3为冷空气流量,W4为末端送风流量,流量单位均为kg/s;T1为TAPRV出口空气温度,T2为TAV入口空气温度,T3为TAV出口空气温度,T4为冷空气温度,Td为DTS处空气温度,T5为送风温度,Tz为CTS处空气温度,温度单位均为℃;P1为TAPRV入口空气压力,Pv为TAPS处空气压力,压力单位均为kPa。

建模时对典型配平空气系统的已知参数做一些常规化取值,具体如下:

1) 稳定状态参数

评估动态特性之初应选取某个稳定状态,观察参数阶跃变化时的动态过程。这里选取稳定状态为:①配平空气供气温度为400 ℉(204.44℃),供气压力为60 psi(413.69 kPa),供气压力变化范围为60~70 psi(413.69~482.63 kPa);②混合腔来流温度为35 ℉(1.67 ℃),流量为65 ppm(0.49 kg/s);③舱室内压力为14.7 psi(101.35 kPa),管路送至舱室的热流为4 220 W,舱室温度变化时间常数为120 s。

2) 温度及压力设定值

假定舱室温度设定值为75 ℉(23.89 ℃),温度变化范围为75~80 ℉(23.89~26.67 ℃);配平空气压力设定值为19.7 psi(135.83 kPa)。

3) 活门参数

假定TAPRV和TAV均采用步进电机驱动的3.5 in(5.08 cm)直径蝶阀,活门转角速率为9 (°)/s,转速时间常数为0.2 s。

4) 传感器参数

根据经验,压力传感器TAPS的时间常数取值0.05 s,管路温度传感器DTS的时间常数取值4 s,舱室温度传感器CTS的时间常数取值30 s。

5) 管路参数

TAPRV之后至TAV之前的管路为直径3 in(7.62 cm)、长度10 ft(3.05 m)、管壁厚度0.05 in(0.127 cm)的不锈钢管路;TAV之后至混合腔来流之前的管路为直径3 in(7.62 cm)、长度2 ft(0.61 m)、管壁厚度0.05 in(0.127 cm)的不锈钢管路;混合腔来流之后至送入舱室之前的管路为直径5 in(12.7 cm)、长度10 ft(3.05 m)、管壁厚度0.05 in(0.127 cm)的铝制管路。假定管路温度变化的时间常数为30~60 s。

1.2 理论公式推导

根据给定的参数能够推导出管路内的压力值及温度值等变量与其他参数之间的传递函数,完成物理模型向功能模型的转化。

1.2.1 稳态下的参数计算

稳定状态下管路中的参数有如下关系:

$ {T_{10}} = {T_{20}} = {T_{30}} $ (1)
 
$ {T_{40}} = {T_{{\rm{d}}0}} = {T_{50}} = \frac{{{W_{30}}{T_{60}} + {W_{20}}{T_{30}}}}{{{W_{30}} + {W_{20}}}} $ (2)
 
$ {W_{10}} = {W_{20}} $ (3)
 

式中:T为空气温度,℃;W为空气流量,kg/s;角标0指稳定状态。

对于输送至舱室的空气,有

$ \begin{array}{l} Q = {c_p}{W_{40}}\left( {{T_{z0}} - {T_{50}}} \right) = \\ \;\;\;\;\;\;{c_p}\left( {{W_{20}} + {W_{30}}} \right)\left( {\frac{{{W_{30}}{T_{60}} + {W_{20}}{T_{30}}}}{{{W_{30}} + {W_{20}}}}} \right) \end{array} $ (4)
 

式中:Q为热通量,W;cp为定压比热容,J/(kg·℃)。

代入已知条件可求得稳态下的配平管供气流量W20、混合空气流量W40及管路温度传感器温度Td0。根据活门流量与开度对应曲线可查得稳定状态下两个活门的开度值。

1.2.2 流量计算公式的线性化处理

经过活门的流量计算公式为[18]

$ W = K\frac{{A{P_{\rm{u}}}}}{{\sqrt {{T_{\rm{u}}}} }}f\left( {\frac{{{P_{\rm{u}}}}}{{{P_{\rm{d}}}}}} \right) $ (5)
 
$ f\left( {\frac{{{P_{\rm{u}}}}}{{{P_{\rm{d}}}}}} \right) = \sqrt {\frac{{2\gamma g}}{{\left( {\gamma - 1} \right)R}}} \cdot \sqrt {{{\left( {\frac{{{P_{\rm{u}}}}}{{{P_{\rm{d}}}}}} \right)}^{\frac{2}{\gamma }}} - {{\left( {\frac{{{P_{\rm{u}}}}}{{{P_{\rm{d}}}}}} \right)}^{\frac{{(1 + \gamma )}}{\gamma }}}} $ (6)
 

式中:A为流通面积,m2P为管路压力,角标u指元件上游,角标d指元件下游,Pa;γ为比热容比,对于双原子气体,有γ=1.4;g为重力加速度,m/s2R为理想气体常数,J/(kg·℃)

动态过程中,压力、温度和活门截面积均不是定值,因此需要对经过活门的流量做线性化处理[19]。采用的方法叫切线法或小偏差法,适合具有连续变化的非线性特性函数。对于有多个自变量的某工作点0点的非线性函数y=f(x1, x2, x3, …),有增量线性化方程:

$ \begin{array}{l} \Delta y = {\left( {\frac{{\partial f}}{{\partial {x_1}}}} \right)_0}\Delta {x_1} + {\left( {\frac{{\partial f}}{{\partial {x_2}}}} \right)_0}\Delta {x_2} + \\ \;\;\;\;\;{\left( {\frac{{\partial f}}{{\partial {x_3}}}} \right)_0}\Delta {x_3} + \cdots \end{array} $ (7)
 

故经过活门流量的增量线性化方程为

$ \begin{array}{l} \Delta W = {\left( {\frac{{\partial W}}{{\partial {P_{\rm{u}}}}}} \right)_0}\Delta {P_{\rm{u}}} + {\left( {\frac{{\partial W}}{{\partial {P_{\rm{d}}}}}} \right)_0}\Delta {P_d} + \\ \;\;\;\;\;\;{\left( {\frac{{\partial W}}{{\partial A}}} \right)_0}\Delta A + {\left( {\frac{{\partial W}}{{\partial {T_{\rm{u}}}}}} \right)_0}\Delta {T_{\rm{u}}} = \\ \;\;\;\;\;\;\left( {1 + {K_{{\rm{c}}0}}} \right){\left( {\frac{W}{{{P_{\rm{u}}}}}} \right)_0}\Delta {P_{\rm{u}}} - {K_{{\rm{c}}0}}{\left( {\frac{W}{{{P_{\rm{d}}}}}} \right)_0}\Delta {P_{\rm{d}}} + \\ \;\;\;\;\;\;{\left( {\frac{W}{A}} \right)_0}\Delta A - \frac{1}{2}{\left( {\frac{W}{{{T_{\rm{u}}}}}} \right)_0}\Delta {T_{\rm{u}}} \end{array} $ (8)
 
$ {K_{{\rm{c}}0}} = \frac{{(\gamma - 1)/2\gamma }}{{\left( {\frac{{{P_{\rm{u}}}}}{{{P_{\rm{d}}}}}} \right)_0^{1 - \frac{1}{\gamma }} - 1}} - \frac{1}{\gamma } $ (9)
 

计算Kc0时,应根据x=Pu/Pd的值确定是否超过超声速点,如果超过超声速点则Kc0=0。

1.2.3 活门流量增量传递函数

TAPRV的流通面积为

$ {A_1} = \frac{{{\rm{ \mathsf{ π} }}{d_1}}}{{{\lambda _1}}}\left( {1 - \cos {\theta _1}} \right) $ (10)
 

式中:d1为TAPRV之后至TAV之前的管路直径,m。

对该流通面积做线性化处理并求拉普拉斯变换,得到

$ \Delta {A_1}(s) = \frac{{{\rm{ \mathsf{ π} }}{d_1}}}{4}\sin {\theta _{10}}\Delta {\theta _1}(s) $ (11)
 

对线性化后的TAPRV增量做拉普拉斯变换,得到

$ \begin{array}{l} \Delta {W_1}(s) = \left( {1 + {K_{{\rm{c1}}}}} \right)\left( {\frac{{{W_{10}}}}{{{P_{10}}}}} \right)\Delta {P_1}(s) - \\ \;\;\;\;\;\;{K_{{\rm{cl}}}}\left( {\frac{{{W_{10}}}}{{{P_{{\rm{vo}}}}}}} \right)\Delta {P_{\rm{v}}}(s) + \left( {\frac{{{W_{10}}}}{{{A_{10}}}}} \right)\Delta {A_1}(s) - \\ \;\;\;\;\;\;0.5\left( {\frac{{{W_{10}}}}{{{T_{10}}}}} \right)\Delta {T_1}(s) \end{array} $ (12)
 

由于上游供气温度稳定,故ΔT1(s)=0。由已知条件可求得Kc1=0,即得到TAPRV流量增量ΔW1关于供气压力P1和TAPRV开度θ1的传递函数:

$ \begin{array}{l} \Delta {W_1}(s) = \left( {1 + {K_{{\rm{cl}}}}} \right)\left( {\frac{{{W_{10}}}}{{{P_{10}}}}} \right)\Delta {P_1}(s) + \\ \;\;\;\;\;\;\;\frac{{{\rm{ \mathsf{ π} }}{d_1}}}{4}\left( {\frac{{{W_{10}}}}{{{A_{10}}}}} \right)\sin {\theta _{10}}\Delta {\theta _1}(s) \end{array} $ (13)
 

同理可得到TAV流量增量ΔW2关于上游压力Pv和TAV开度θ2的传递函数:

$ \begin{array}{l} \Delta {W_2}(s) = \left( {1 + {K_{{\rm{c2}}}}} \right)\left( {\frac{{{W_{20}}}}{{{P_{{\rm{v}}0}}}}} \right)\Delta {P_{\rm{v}}}(s) + \\ \;\;\;\;\;\;\frac{{{\rm{ \mathsf{ π} }}{d_2}}}{4}\left( {\frac{{{W_{20}}}}{{{A_{20}}}}} \right)\sin {\theta _{20}}\Delta {\theta _2}(s) \end{array} $ (14)
 
1.2.4 管路压力增量的传递函数

由压力的动态特性计算公式,可得

$ \frac{{{\rm{d}}{P_{\rm{v}}}}}{{{\rm{d}}t}} = \frac{{\gamma R{T_1}}}{{{V_{1 - 2}}}}\left( {{W_1} - {W_2}} \right) $ (15)
 

式中:V为管段体积,m3

对其做线性化处理并求拉氏变换,得到

$ s\Delta {P_{\rm{v}}}(s) = \frac{{\gamma R{T_{10}}}}{{{V_{1 - 2}}}}\left[ {\Delta {W_1}(s) - \Delta {W_2}(s)} \right] $ (16)
 

代入传递函数ΔW1(s)和ΔW2(s)并写成传递函数形式,得到

$ \begin{array}{l} \Delta {P_{\rm{v}}}(s) = \frac{\varepsilon }{{\tau s + 1}}\left[ {\left( {1 + {K_{{\rm{cl}}}}} \right)\left( {\frac{{{W_{10}}}}{{{P_{10}}}}} \right)\Delta {P_1}(s) + } \right.\\ \;\;\;\;\;\;\frac{{{\rm{ \mathsf{ π} }}{d_1}}}{4}\left( {\frac{{{W_{10}}}}{{{A_{10}}}}} \right)\sin {\theta _{10}}\Delta {\theta _1}(s) - \\ \left. {\;\;\;\;\;\;\frac{{{\rm{ \mathsf{ π} }}{d_2}}}{4}\left( {\frac{{{W_{20}}}}{{{A_{20}}}}} \right)\sin {\theta _{20}}\Delta {\theta _2}(s)} \right] \end{array} $ (17)
 

式中:时间常数τ=$\varepsilon \left( \frac{{{V}_{1-2}}}{\gamma R{{T}_{10}}} \right)$ε=$\frac{{{P}_{\text{v}0}}}{{{K}_{\text{cl}}}{{W}_{10}}+{{K}_{\text{c}2}}{{W}_{20}}+{{W}_{20}}}$

代入已知条件,即得到压力关于输入压力、两个活门开度的传递函数。

1.2.5 管路温度传递函数

管路内空气温度传递函数为

$ F(s) = \frac{{\Delta {T_{\rm{d}}}(s)}}{{\Delta {T_{\rm{u}}}(s)}} = {k_{\rm{g}}}\frac{{{\tau _{\rm{d}}}s + 1}}{{{\tau _{\rm{u}}}s + 1}} $ (18)
 

式中:kg为根据稳态热损单独评估,若无损管路,kg=1;${{\tau }_{\text{u}}}=\frac{{{M}_{\text{m}}}{{c}_{\text{m}}}}{{{W}_{\text{f}}}{{c}_{\text{f}}}(1-U)}$,角标m为导管本体,角标f为流体,Mm为导管质量,kg;τd=uU=exp(-Hf/Wfcf), Hf=Nu(λ/dduct)AsurfaceAsurface为管路内壁表面积, Nu=0.023Ref0.8Prf0.3,当流体为空气时,Prf近似为0.7。

计算式(18)中的各参数值,可得到每个管段的传递函数。

1.2.6 混合位置温度增量的传递函数

配平空气与混合腔来流的混合点的温度为

$ {T_4} = \frac{{{W_3}{T_6} + {W_2}{T_3}}}{{{W_3} + {W_2}}} $ (19)
 

认为该式中W3T6恒定,做线性化处理后的拉氏变换,得到混合位置温度增量的传递函数:

$ \begin{array}{l} \Delta {T_4}(s) = \frac{{{W_{30}}\left( {{T_{30}} - {T_{60}}} \right)}}{{{{\left( {{W_{30}} + {W_{20}}} \right)}^2}}}\left[ {\left( {1 + {K_{c2}}} \right) \cdot } \right.\\ \;\;\;\;\;\;\;\left. {\left( {\frac{{{W_{20}}}}{{{P_{{\rm{v}}0}}}}} \right)\Delta {P_{\rm{v}}}(s) + \frac{{{\rm{ \mathsf{ π} }}{d_2}}}{4}\left( {\frac{{{W_{20}}}}{{{A_{20}}}}} \right)\sin {\theta _{20}}\Delta {\theta _2}(s)} \right] + \\ \;\;\;\;\;\;\;\frac{{{W_{20}}}}{{{W_{30}} + {W_{20}}}}\Delta {T_3}(s) \end{array} $ (20)
 
1.2.7 舱室温度增量的传递函数

舱室温度Tz的动态特性方程为

$ {M_{\rm{z}}}{c_{\rm{z}}}\left( {\frac{{{\rm{d}}{T_{\rm{z}}}}}{{{\rm{d}}t}}} \right) = \left( {{T_5} - {T_{\rm{z}}}} \right){W_4}{c_p} + Q $ (21)
 

对该式做线性化处理后拉氏变换,得到

$ \begin{array}{*{20}{l}} {\left( {{\tau _{\rm{z}}}s + 1} \right)\Delta {T_{\rm{z}}}(s) = \Delta {T_5}(s) + }\\ {\;\;\;\;\;\;\;\frac{{\left( {{T_{50}} - {T_{{\rm{z}}0}}} \right){c_p}}}{{{W_4}{c_p}}}\Delta {W_4} + \frac{1}{{{W_4}{c_p}}}\Delta Q(s)} \end{array} $ (22)
 

式中:${{\tau }_{\text{z}}}=\frac{{{M}_{\text{z}}}{{c}_{\text{z}}}}{{{W}_{4}}{{c}_{p}}}$,该值可预设。

对于W4,有W4=W2+W3,对其做线性化处理后拉氏变换,得到

$ \Delta {W_4}(s) = \Delta {W_2}(s) + \Delta {W_3}(s) $ (23)
 

由于W3为定值,故ΔW3(s)=0,因此:

$ \begin{array}{l} \Delta {T_{\rm{z}}}(s) = \frac{1}{{{\tau _{\rm{z}}}s + 1}}\left\{ {\Delta {T_5}(s) + \frac{{\left( {{T_{50}} - {T_{{\rm{z0}}}}{c_p}} \right)}}{{{W_4}{c_p}}}} \right.{\rm{ }} \cdot \\ \;\;\;\;\; \left[ {\left( {1 + {K_{{\rm{c}}2}}} \right)\left( {\frac{{{W_{20}}}}{{{P_{{\rm{v}}0}}}}} \right)\Delta {P_{\rm{v}}}(s) + \frac{{{\rm{ \mathsf{ π} }}{d_2}}}{4}\left( {\frac{{{W_{20}}}}{{{A_{20}}}}} \right)} \right. \cdot \\ \;\;\;\;\;\left. {\left. {\sin {\theta _{20}}\Delta {\theta _2}(s)} \right] + \frac{1}{{{W_4}{c_p}}}\Delta Q(s)} \right\} \end{array} $ (24)
 

若热载荷稳定,则ΔQ(s)=0,可得到舱室温度增量的传递函数。

1.2.8 电动式活门的活门开度传递函数

电动式活门与机械气动式活门的区别在于,机械气动式活门以控制活门开度来控制系统参数,而电动式活门以控制电机的输入信号来控制活门开度,从而控制系统参数,增加了输入信号对活门开度的控制。因此,上述以活门开度作为自变量的推导过程对于机械气动式活门同样适用,但电动式活门还应增加活门开度传递函数的推导。

假定电动式TAPRV的最大扭转速率为ω10,时间常数为τ1,则对于输入信号u1,有

$ {\tau _1}\left( {\frac{{{\rm{d}}{\omega _1}}}{{{\rm{d}}t}}} \right) + {\omega _1} = {\omega _{10}}{u_1} $ (25)
 
$ \frac{{{\rm{d}}{\theta _1}}}{{{\rm{d}}t}} = {\omega _1} $ (26)
 

对式(25)和式(26)做拉普拉斯变换并代入已知条件,做线性化处理得到增量传递函数:

$ \Delta {\theta _1}(s) = \frac{{{\omega _{10}}}}{{s\left( {{\tau _1}s + 1} \right)}}\Delta {u_1}(s) $ (27)
 

同理,假定电动式TAV的最大扭转速率为ω20,时间常数为τ2,输入信号u2,有增量传递函数:

$ \Delta {\theta _2}(s) = \frac{{{\omega _{20}}}}{{s\left( {{\tau _2}s + 1} \right)}}\Delta {u_2}(s) $ (28)
 
1.2.9 传感器预设值增量传递函数

根据传感器时间常数可快速得到预设值的传递函数。即压力传感器的压力增量传递函数为

$ \Delta {P_{{\rm{vs}}}}(s) = \frac{1}{{0.05s + 1}}\Delta {P_{\rm{v}}}(s) $ (29)
 

管路温度传感器的温度增量传递函数为

$ \Delta {T_{{\rm{ds}}}}(s) = \frac{1}{{3.6s + 1}}\Delta {T_{\rm{d}}}(s) $ (30)
 

式中:ΔTd(s)可以由管路模型得到与ΔT4(s)的关系。

舱室温度传感器的温度增量传递函数为

$ \Delta {T_{{\rm{zs}}}}(s) = \frac{1}{{30s + 1}}\Delta {T_{\rm{z}}}(s)$ (31)
 
1.3 线性化模型搭建

上述推导过程简而言之是分别建立两个活门处流量与压力、开度的传递函数以及各管段温度的传递函数,由此推导出传感器测点的输出,将其作为控制器的输入。模型的搭建也就是将这些传递函数关系连接起来。在Simulink软件中搭建线性化模型,如图 2所示。本文首先基于活门开度阶跃(可以认为是机械式活门调节过程)运行典型配平系统传递函数模型,输出各参数的动态特性;然后加入电动式活门开度传递函数,设计控制器中的压力控制律,完成控制架构的闭环运行。

图 2 小偏差化动态模型图 Fig. 2 Diagram of small deviation linearization dynamic model
2 典型配平空气系统仿真结果分析及校验 2.1 系统参数的时域动态特性

计算TAPRV开度阶跃变化(TAV开度固定)或TAV开度阶跃变化(TAPRV开度固定)时,900 s内管路压力值、管路温度值、舱室温度值的时域动态特性分别如图 3(a)图 3(b)所示。可以看到,TAPRV开度的阶跃变化会导致TAPRV下游压力增加,从而热空气流量增加,管路温度升高,带动舱室温度升高;900 s后各项参数基本稳定,其中压力响应非常迅速,管路内气体温度在150 s左右达到稳定,舱室温度在750 s左右达到稳定。然而,TAV开度的阶跃变化除了使TAPRV下游压力降低以外,管路温度和舱室温度几乎没有太大变化。经分析认为,由于上游存在TAPRV的节流作用,TAV开度的单位变化并不能大幅度地体现在对管路温度的调节上。

图 3 活门开度阶跃变化时900 s内时域动态特性 Fig. 3 Time domain dynamic characteristics in 900 s after valve angel step changing

为进一步观察,截取TAPRV或TAV开度阶跃变化时,前5 s内各参数值的时域动态响应如图 4所示。可以清晰地看到管路压力的响应迅速,在0.8 s左右基本达到稳定。此外,观察管路内的温度变化过程,TAPRV开度阶跃变化时管路内温度有较大提升,而TAV开度阶跃变化时管路内温度略有升高。当TAPRV或TAV开度阶跃变化时,舱室温度在前5 s内几乎无变化。

图 4 活门开度阶跃变化时5 s内时域动态特性 Fig. 4 Time domain dynamic characteristics in 5 s after valve angel step changing
2.2 结果校验

将时域动态响应结果与测试结果对比,以校验动态模型准确性。测试时分别实施两个活门的开度调节,即

1) 调节TAPRV开度由7.706°向8.706°阶跃变化1°(TAV开度固定在14.22°)。

2) 调节TAV开度由14.22°向15.22°阶跃变化1°(TAPRV开度固定在7.706°)。

图 5为测试结果。值得说明的是,测试结果表现的是真实的物理量,和时域幅值不同,因此图线的显示略有不同。模拟结果和测试结果的参数变化曲线吻合,证明模型有效。为了直观地对比结果,表 1中展示了仿真时域幅值与真实参数变化量的误差分析。

图 5 活门开度阶跃变化时测试结果 Fig. 5 Test results after valve angel step changing
表 1 误差分析 Table 1 Error analysis
参数 TAPRV开度阶跃 TAV开度阶跃
时域幅值 参数变化量 误差 时域幅值 参数变化量 误差
舱室温度/°F 5.92 6.1 3.0% 0.000 457 0
管路温度/°F 6.18 6.5 4.9% 0.000 976 0
下游压力/psi 2.65 3.0 11.7% -1.43 -1.2 19.2%
3 电动式压力调节活门控制律设计实例 3.1 控制架构

基于典型配平空气系统传递函数中TAPRV下游压力的时域动态响应,采用系统辨识方法求得TAPRV开度对管路压力传感器传递函数为

$ G(s) = \frac{{385.7}}{{{s^2} + 27.26s + 145.2}} $ (32)
 

由于TAPRV下游压力响应非常迅速,可以采用单输入单输出(Simple Input Simple Output,SISO)系统直接反馈的控制方式[20]。控制架构如图 6所示。

图 6 TAPRV下游压力控制架构 Fig. 6 TAPRV downstream pressure control configuration

图 6中,G1(s)和G2(s)为压力链路传递函数,K(s)为控制器内的压力调节控制律。由2.1节可知:

$ \begin{array}{l} {G_1}(s){G_2}(s) = \frac{{385.7}}{{{s^2} + 27.26s + 145.2}}{\theta _1}(s) = \\ \;\;\;\;\;\;\;\;\;\frac{{385.7}}{{{s^2} + 27.26s + 145.2}} \cdot \frac{9}{{s(0.2s + 1)}}{u_1}(s) \end{array} $ (33)
 

应注意的是式(33)中将自变量替换为控制电机输入信号,然后采用MATLAB里的控制器设计工具箱对式(33)中的传递函数设计控制律。

3.2 控制目标

压力控制目标的制定应根据系统原有特性来实施,主要提出TAPRV下游压力的时域响应时间常数及超调量。在设置控制目标时,基于舱室稳定时间逐步推导上游各点时间常数,最后确定压力控制时间常数。

由于舱室温度变化的时间常数为120 s,舱室温度传感器的时间常数为30 s,管路温度变化的时间常数为30~60 s;因此,总的温度变化时间常数为120 s+30 s+60 s=210 s。该时间常数略大,这里提出响应速度提升3倍的控制目标,即舱室温度的时间常数应在70 s左右。

民机配平空气系统的温度控制除座舱温度控制回路外,还有管道温度控制回路作为内环控制回路。为避免控制间的相互影响,提出管路温度的响应速度快于座舱温度响应速度5倍的控制目标,即管路温度的时间常数应在14 s左右。

同样地,为保证压力控制快于温度控制,提出管路压力的响应速度快于管路温度响应速度5倍的控制目标。制定TAPRV下游压力的控制目标如下:

1) 时间常数为2.8 s左右。

2) 超调量不超过10%。

3) 带宽频率为1/2.8 s=0.357 rad/s左右。

4) 无稳态误差。

3.3 控制律设计

在设计电动式TAPRV控制律时,应将式(33)作为被控对象(自变量为输入信号),设计K(s)的传递函数使输出结果满足控制目标。

被控对象的开环输出结果如图 7(a)所示。由时域响应可以看出被控对象输出发散,无法达到稳定。为使带宽频率向前调节至0.357 rad/s左右,幅相曲线应下降37 dB左右。经调试得到比例环节K(s)=1/67。

图 7 控制律设计过程 Fig. 7 Design process of control law

此时输出的闭环动态特性如图 7(b)所示。由时域响应可以看出此时输出稳定,时间常数为2.827 s,满足控制目标要求。

3.4 控制律合理性校验

将调试成功的控制律放入图 2中的控制器模型中做电动式TAPRV的闭环运行,分别对不同输入条件加载阶跃变化并观察输出结果,用以校验控制律设计的合理性。

图 8中分别是压力设定值、供气温度、供气压力及舱室热载荷输入阶跃变化时系统各参数的900 s输出结果。其中,图 8(a)中压力设定值阶跃变化的工况可以等同于机械气动式TAPRV活门开度阶跃变化的工况。可以看到图 8(a)图 3(a)中的曲线趋势吻合度较好,证明电动式TAPRV的控制器通过控制输入信号达到了控制活门开度的良好效果。另外,从图 8(b)~图 8(d)中可以看到,不同的系统外界条件突变时均可保证系统参数快速稳定,认为控制律设计合理。

图 8 控制律校验输出结果 Fig. 8 Output results of control law verification
4 结论

本文通过小偏差化的线性化处理仿真建模手段建立了民机配平空气系统中电动式压力调节活门辅以压力传感器的压力控制方案的动态仿真模型。证明方法可用之后,依据典型系统的动态特性完成了电动式压力调节活门的控制律设计。

1) 建立典型配平空气系统的线性化传递函数模型,以活门开度作为自变量。

2) 建立电动式压力调节活门开度的线性化传递函数,以控制电机的输出信号作为自变量。

3) 设计控制架构及控制目标,计算被控对象传递函数。

4) 调试控制律,使得系统输出稳定并满足控制目标。

5) 将控制律放入控制器模型,运行带有电动式压力调节活门的闭环模型,校验控制律合理性。

参考文献
[1] 朱春玲. 飞行器环境控制与安全救生[M]. 北京: 北京航空航天大学出版社, 2006: 322-326.
ZHU C L. Air vehicle environmental control and safety lifesaving[M]. Beijing: Beihang University Press, 2006: 322-326. (in Chinese)
[2] TERESA J L, ISABEL P G. A thermoeconomic analysis of a commercial aircraft environmental control system[J]. Applied Thermal Engineering, 2005, 25(2-3): 309-325.
Click to display the text
[3] ZHAO H L, HOU Y, ZHU Y F, et al. Experimental study on the performance of an aircraft environmental control system[J]. Applied Thermal Engineering, 2009, 29(16): 3284-3288.
Click to display the text
[4] 陈元先. 旅客机环境控制系统的发展[J]. 航空学报, 1999, 20(S1): 8-10.
CHEN Y X. Evolution of the environment control system for commercial aircraft[J]. Acta Aeronautica et Astronautica Sinica, 1999, 20(S1): 8-10. (in Chinese)
Cited By in Cnki (21) | Click to display the text
[5] 高铭. 民用飞机配平系统管路出口流场研究[J]. 科技创新导报, 2017, 11: 19-22.
GAO M. Study on flow field of pipeline outlet of civil aircraft trim air system[J]. Science and Technology Innovation Herald, 2017, 11: 19-22. (in Chinese)
Cited By in Cnki | Click to display the text
[6] 王荣, 唐敏杰. 空调系统热空气配平活门疑难故障分析[J]. 民航科技, 2010, 3(139): 106-108.
WANG R, TANG M J. Difficult fault analysis of trim air pressure regulation valve in air conditioning system[J]. Civil Aviation Science & Technology, 2010, 3(139): 106-108. (in Chinese)
Cited By in Cnki | Click to display the text
[7] DECHOW M, NURCOMBE C A H. Aircraft environmental control systems[M]. Berlin Heidelberg: Springer, 2005: 3-24.
[8] TU Y, LIN G P. Dynamic simulation of aircraft environmental control system based on flowmaster[J]. Journal of Aircraft, 2011, 48(6): 2031-2041.
Click to display the text
[9] 郭志刚.飞机气动活门类附件的特点[C]//2008中国航空维修论坛论文集, 2008: 153-159.
GUO Z G. Characteristics of pneumatic valves on aircraft[C]//2008 China Aviation Maintenance Conference Proceeding, 2008: 153-159(in Chinese).
Cited By in Cnki (1) | Click to display the text
[10] KNIGGE H, WORNER M, SCHMITZ G. Simulation with modelica for a dynamic representation of an aircraft cabin climate for comfort improved climate control[C]//25th International Congress of the Aeronautical Sciences, 2006: 1-7.
[11] LINARES D P. Modeling and simulation of an aircraft environmental control system[D]. Montréal: école Polytechnique De Montréal, 2016: 10-25.
Click to display the text
[12] 杨绍曾, 王浚. 飞机座舱温度控制系统动态特性研究[J]. 航空学报, 1983, 4(4): 78-85.
YANG S Z, WANG J. Study on dynamic characteristics of cockpit temperature control system[J]. Acta Aeronautica et Astronautica Sinica, 1983, 4(4): 78-85. (in Chinese)
Cited By in Cnki | Click to display the text
[13] PéREZ-GRANDE I, LEO T J. Optimization of a commercial aircraft environmental control system[J]. Applied Thermal Engineering, 2002, 22(17): 1885-1904.
Click to display the text
[14] 任明波, 王娟, 李荣军, 等. 大型飞机座舱温度控制系统控制律设计[J]. 航空学报, 2017, 38(S1): 721501.
REN M B, WANG J, LI R J, et al. Control law design for temperature control system of large-scale aircraft cabin[J]. Acta Aeronautica et Astronautica Sinica, 2017, 38(S1): 721501. (in Chinese)
Cited By in Cnki | Click to display the text
[15] 何君, 赵竞全, 袁修干. 飞机环境控制系统的模糊控制研究[J]. 北京航空航天大学学报, 2004, 30(12): 1511-1514.
HE J, ZHAO J Q, YUAN X G. Fuzzy control of aircraft environmental control system[J]. Journal of Beijing University of Aeronautics and Astronautics, 2004, 30(12): 1511-1514. (in Chinese)
Cited By in Cnki (6) | Click to display the text
[16] 屠毅, 林贵平. 大型飞机座舱温度控制系统仿真研究[J]. 航空学报, 2011, 32(1): 49-57.
TU Y, LIN G P. Simulation of the large scale aircraft cabin temperature control system[J]. Acta Aeronautica et Astronautica Sinica, 2011, 32(1): 49-57. (in Chinese)
Cited By in Cnki (10) | Click to display the text
[17] 姚洪伟, 王浚. 歼击机环境控制系统控制性能分析[J]. 中国工程科学, 2006, 8(6): 44-47.
YAO H W, WANG J. Analysis of control performance of ECS in fighter plane[J]. Engineering Science, 2006, 8(6): 44-47. (in Chinese)
Cited By in Cnki (2) | Click to display the text
[18] 沈维道, 童均耕. 工程热力学[M]. 北京: 高等教育出版社, 2007: 246-252.
SHEN W D, TONG J G. Engineering thermodynamic[M]. Beijing: High Education Press, 2007: 246-252. (in Chinese)
[19] 胡寿松. 自动控制原理[M]. 6版. 北京: 科学出版社, 2018: 21-29.
HU S S. Principle of automatic control[M]. 6th ed. Beijing: China Science Publishing Press, 2018: 21-29. (in Chinese)
[20] 王建辉, 顾树生. 自动控制原理[M]. 2版. 北京: 清华大学出版社, 2017: 41-57.
WANG J H, GU S S. Principle of automatic control[M]. 2nd ed. Beijing: Tsinghua University Press, 2017: 41-57. (in Chinese)
http://dx.doi.org/10.7527/S1000-6893.2019.22951
中国航空学会和北京航空航天大学主办。
0

文章信息

张絮涵, 吴成云
ZHANG Xuhan, WU Chengyun
民机配平空气系统小偏差化动态建模仿真方法
Small deviation linearization method for dynamic modeling and simulation of trim air system in civil aircraft
航空学报, 2019, 40(11): 122951.
Acta Aeronautica et Astronautica Sinica, 2019, 40(11): 122951.
http://dx.doi.org/10.7527/S1000-6893.2019.22951

文章历史

收稿日期: 2019-02-14
退修日期: 2019-03-25
录用日期: 2019-04-22
网络出版时间: 2019-05-09 08:50

相关文章

工作空间