CN105715384A - 改善航空发动机加速性能的增量式pi控制器参数调整方法 - Google Patents
改善航空发动机加速性能的增量式pi控制器参数调整方法 Download PDFInfo
- Publication number
- CN105715384A CN105715384A CN201410736305.9A CN201410736305A CN105715384A CN 105715384 A CN105715384 A CN 105715384A CN 201410736305 A CN201410736305 A CN 201410736305A CN 105715384 A CN105715384 A CN 105715384A
- Authority
- CN
- China
- Prior art keywords
- need
- fuel flow
- controller
- increment
- final
- 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
Landscapes
- Electrical Control Of Air Or Fuel Supplied To Internal-Combustion Engine (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
本发明属于航空发动机数控系统领域,提供一种改善航空发动机加速性能的增量式PI控制器参数调整方法,用于调整PI控制器的Kp和Ti,使发动机获取优良的加速特性。技术方案包括:将设计值Kp_design、Ti_design缩放,得到Kp0、Ti0,进行加速试验,获取试验数据,提取满足WF=WFacc、dup<0的数据,计算dWF_need、dui_need、K_need、K_rsv,根据分离点A转速N2r0对K分段取值,然后得到Ti_final、Kp_final,并写入PI控制器,进行台架试车。
Description
技术领域:
本发明属于航空发动机数控系统领域。
背景技术:
加速性能是航空发动机一项非常重要的性能指标,减小发动机加速时间是提高飞机作战敏捷性的重要因素之一。航空发动机加速性能完全取决于加速过程中的供油情况。
发动机加速过程中的主燃油供油量计算方法如图1所示。当飞行员在1s内快速将油门杆从慢车状态推动至中间状态,主燃油流量控制算法使得高压转子转速给定值N2_dem(k)(k为系统采样时刻)从慢车状态控制计划突升到中间状态控制计划,使得高压转子转速N2的控制误差“eN2(k)=N2_dem(k)-N2(k)”突升,进而导致增量式PI(比例、积分)控制器计算出的当前周期主燃油流量增量“dWF(k)=Kp×[eN2(k)-eN2(k-1)]+Kp/Ti×eN2(k)”突升,主燃油流量总量“WF_PI(k)=WF_PI(k-1)+dWF(k)”也跟着突升。主燃油流量控制算法将WF_PI与加速燃油流量限制线WFacc进行低选,得到主燃油供油量WF。一般在加速过程的初期dWF较大,后期dWF较小,导致初期按照加速燃油流量限制线供油,即WF=WFacc,后期按照增量式PI控制器运算结果供油,即WF=WF_PI。加速过程中的主燃油供油变化规律如图2所示,其中A为分离点,即在分离点之前按照加速燃油流量限制线WFacc供油,在分离点之后按照增量式PI控制器运算结果WF_PI供油。
在发动机加速过程中,分离点A是个关键因素。如果WF与WFacc过早分离,则加速性能不会太好,加速时间偏长;如果WF与WFacc分离晚,虽然加速时间短,但是发动机转子转速会产生较大的超调,影响发动机的安全。
目前国内工程应用中加速过程主燃油供油量分离点A靠增量式PI控制器的运算结果WF_PI与加速燃油流量限制线WFacc的自适应低选确定,当主机单位一旦确定加速燃油流量限制线后,发动机的加速性能就完全取决于增量式PI控制器的运算结果,而增量式PI控制器的运算结果又唯一取决于其控制参数Kp、Ti。因此,PI控制参数的确定是影响发动机加速性能的关键所在。
主燃油增量式PI控制器控制参数Kp、Ti设计时通常只能考虑稳态时频域内的带宽特性和稳定裕度储备,不能直接反映加速特性,因此台架试车时发动机的加速特性通常不会令人满意,需要根据试验结果调整PI控制参数以改善加速特性。目前国内工程应用中调整PI控制参数以改善发动机加速特性的方法欠缺,通常采用增大比例增益Kp或者减小积分时间常数Ti的方法,但不能很好地把握,需要反复试凑。增大Kp能适当改善加速特性,但是系统稳定裕度储备会降低,而且Kp增大到一定程度后加速特性反而会降低,加速过程中由于加速燃油流量限制线WFacc的干扰,同时增量式PI控制器比例环节计算出的燃油流量增量“dup(k)=Kp×[eN2(k)-eN2(k-1)]”为负值(加速过程中随着转速的上升eN2(k)逐渐减小),因此会阻碍加速过程。减小积分时间常数Ti能改善发动机加速特性,但会降低系统稳定裕度储备,而且降低多少很难评估,同时分离点A也难以准确预估。
发明内容:
发明目的:提供一种改善航空发动机加速性能的增量式PI控制器控制参数调整方法,用于快速调整主燃油流量PI控制器的比例增益Kp和积分时间常数Ti,以使发动机获取优良的加速特性。
技术方案:
一种改善航空发动机加速性能的增量式PI控制器参数调整方法,包括:
(1)将主燃油流量增量式PI控制器的比例增益设计值Kp_design乘以比例因子kKp,作为加速性能摸底试验的比例增益初始值Kp0:
(2)将积分时间常数设计值Ti_design除以缩放因子kTi,结果与0.2比较,取二者中最大值作为加速性能摸底试验的积分时间常数初始值Ti0:
(3)将比例增益Kp0、积分时间常数Ti0写入增量式PI控制器作为发动机台架加速性能摸底试验的控制参数初始值,获取试验数据;
(4)从所述试验数据中提取同时满足WF=WFacc、dup<0条件的高压转子换算转速N2r、加速燃油流量限制线WFacc、主燃油供油量WF、增量式PI控制器比例环节计算出的燃油流量增量dup和增量式PI控制器积分环节计算出的燃油流量增量dui,作如下计算:
dWF_need=WFacc-WFacc_pre>0,其中dWF_need是保证主燃油供油走加速燃油流量限制线所需的最小单步燃油流量增量,WFacc_pre是WFacc上一周期的值;
dui_need=dWF_need–dup>0,其中dui_need是除去比例环节燃油流量增量dup的反作用之后,还能保证最小单步燃油流量增量dWF_need所需的积分环节燃油流量增量;
K_need=dui_need/dui<1,其中K_need是所需积分量比实际所用积分量的比例;
K_rsv=K_need×1.3,其中K_rsv是考虑鲁棒性后所需积分量比实际所用积分量的比例;
(5)设期望的加速过程燃油流量分离点A的转速为N2r0,则对积分时间常数缩放因子K作如下取值:
当N2r≤(N2r0-2%)时,K=K_rsv;
当(N2r0-2%)<N2r<N2r0时,K为K_rsv与K_need的线性插值;
当N2r=N2r0时,K=K_need;
当N2r0<N2r<(N2r0+3%)时,K从K_need线性平滑过渡到1;
当N2r0≥(N2r0+3%)时,K=1;
(6)计算得到增量式PI控制器最终所需的积分时间常数:Ti_final=Ti0/K;
(7)计算得到增量式PI控制器最终所需的比例增益Kp_final:
当N2r≤N2r0时,Kp_final=Kp0;
当N2r0<N2r<(N2r0+3%)时,Kp_final从Kp0线性平滑过渡到Kp_design;
当N2r≥(N2r0+3%)时,Kp_final=Kp_design;
(8)将Kp_final、Ti_final写入主燃油流量增量式PI控制器,进行后续台架试车。
其中,比例因子kKp的取值范围为0.5≤kKp≤0.8;缩放因子kTi的取值范围为3≤kTi≤5;加速过程燃油流量分离点A的转速N2r0的取值范围为92%≤N2r0≤95%。
有益效果:
本方法对主燃油流量增量式PI控制器的控制参数设计结果进行一次调整,便可最终确定控制参数,在后续台架试车中既能获取准确的加速过程燃油流量分离点A,保证发动机具有良好的加速性能,又能使增量式PI控制器在具备一定鲁棒性的基础上,将系统的稳定余度储备降低到最小,解决了当前工程应用中试凑法存在的弊端,即既需要反复试凑,又不能准确获取加速过程中的燃油流量分离点A,且很难最大限度得保证系统稳定裕度储备等问题。
附图说明
图1为发动机加速过程中主燃油供油量计算法示意图。
图2为发动机加速过程中主燃油供油变化规律示意图。
具体实施方式
(1)将主燃油流量增量式PI控制器的比例增益设计值Kp_design乘以比例因子kKp(0.5≤kKp≤0.8),作为发动机台架加速性能摸底试验的比例增益初始值,例如Kp0=Kp_design×0.8;
因为后续要减小增量式PI控制器的积分时间常数Ti,这样会降低系统的稳定裕度储备,而在此处将比例增益设计值Kp_design乘以比例因子kKp,能够提高系统的稳定裕度储备,从而一定程度弥补减小积分时间常数Ti对系统稳定裕度储备的影响;
(2)将主燃油流量增量式PI控制器的积分时间常数设计值Ti_design除以缩放因子kTi(3≤kTi≤5),但所得值不能小于0.2,作为发动机加速性能摸底试验的积分时间常数初始值,例如Ti0=max(Ti_design/5,0.2);
按照该方法得到的积分时间常数初始值Ti0会使增量式PI控制器产生一个较强的积分效果(dui=Kp/Ti×eN2(k)),从而导致加速过程中燃油流量分离点A较晚产生,产生转子转速超调(但是由于转子转速限制计划的存在,超调不会太大,因此对发动机台架试车不会有风险),同时系统稳定裕度储备下降较多。
(3)将根据上述步骤获得的比例增益初始值Kp0、积分时间常数初始值Ti0写入增量式PI控制器,作为发动机台架加速性能摸底试验的控制参数初始值,获取试验数据;然后根据试验数据重新调整积分时间常数,以保证获取准确的加速过程燃油流量分离点A,同时在保证一定鲁棒性的基础上积分时间常数尽可能减小得少,使得系统稳定裕度储备在满足加速特性和鲁棒性的基础上尽可能下降得少;
(4)获取试验数据中主燃油流量走加速燃油流量限制线的数据,即满足如下要求:
a)WF与WFacc重合,即WF=WFacc;
b)dup<0,
从满足上述要求的试验数据中抽离出如下参数作为后续计算的基础:
a)N2r:发动机高压转子换算转速;
b)WFacc:加速燃油流量限制线;
c)WF:增量式PI控制器运算结果与加速燃油流量限制线低选后的结果;
d)dup:增量式PI控制器比例环节计算出的燃油流量增量;
e)dui:增量式PI控制器积分环节计算出的燃油流量增量,即dui=Kp/Ti×eN2(k)。
(5)计算保证主燃油供油走加速燃油流量限制线所需的最小单步燃油流量增量:dWF_need=WFacc-WFacc_pre>0,其中WFacc_pre是WFacc上一周期的值;
(6)计算除去比例环节燃油流量增量dup的反作用之后,还能保证最小单步燃油流量增量dWF_need所需的积分环节燃油流量增量:dui_need=dWF_need–dup>0;
(7)计算K_need(所需积分量比实际所用积分量的比例)和K_rsv(考虑鲁棒性后所需积分量比实际所用积分量的比例,在K_need基础上预留1.3倍的余量):
a)K_need=dui_need/dui<1;
b)K_rsv=K_need×1.3;
(8)设期望的加速过程燃油流量分离点A的转速为N2r0(92%≤N2r0≤95%),则对积分时间常数缩放因子K作如下取值:
a)当N2r≤(N2r0-2%)时,K=K_rsv;
b)当(N2r0-2%)<N2r<N2r0时,K为K_rsv与K_need的线性插值;
c)当N2r=N2r0时,K=K_need;
d)当N2r0<N2r<(N2r0+3%)时,K从K_need线性平滑过渡到1;
e)当N2r0≥(N2r0+3%)时,K=1;
(9)计算得到增量式PI控制器最终所需的积分时间常数:Ti_final=Ti0/K;
(10)计算得到增量式PI控制器最终所需的比例增益Kp_final:
a)当N2r≤N2r0时,Kp_final=Kp0;
b)当N2r0<N2r<(N2r0+3%)时,Kp_final从Kp0线性平滑过渡到Kp_design;
c)当N2r≥(N2r0+3%)时,Kp_final=Kp_design;
(11)将最终所得的比例增益Kp_final、积分时间常数Ti_final写入主燃油流量增量式PI控制器,进行后续台架试车,这样在试验时既能获取准确的加速过程燃油流量分离点A,保证发动机具有良好的加速性能,又能使增量式PI控制器在具备一定鲁棒性的基础上(通过积分时间常数固定预留1.3倍的余量),将系统的稳定余度储备降低到最小。
Claims (2)
1.一种改善航空发动机加速性能的增量式PI控制器参数调整方法,其特征在于,包括:
(1)将主燃油流量增量式PI控制器的比例增益设计值Kp_design乘以比例因子kKp,作为加速性能摸底试验的比例增益初始值Kp0:
(2)将积分时间常数设计值Ti_design除以缩放因子kTi,结果与0.2比较,取二者中最大值作为加速性能摸底试验的积分时间常数初始值Ti0:
(3)将比例增益Kp0、积分时间常数Ti0写入增量式PI控制器作为发动机台架加速性能摸底试验的控制参数初始值,获取试验数据;
(4)从所述试验数据中提取同时满足WF=WFacc、dup<0条件的高压转子换算转速N2r、加速燃油流量限制线WFacc、主燃油供油量WF、增量式PI控制器比例环节计算出的燃油流量增量dup和增量式PI控制器积分环节计算出的燃油流量增量dui,作如下计算:
dWF_need=WFacc-WFacc_pre>0,其中dWF_need是保证主燃油供油走加速燃油流量限制线所需的最小单步燃油流量增量,WFacc_pre是WFacc上一周期的值;
dui_need=dWF_need–dup>0,其中dui_need是除去比例环节燃油流量增量dup的反作用之后,还能保证最小单步燃油流量增量dWF_need所需的积分环节燃油流量增量;
K_need=dui_need/dui<1,其中K_need是所需积分量比实际所用积分量的比例;
K_rsv=K_need×1.3,其中K_rsv是考虑鲁棒性后所需积分量比实际所用积分量的比例;
(5)设期望的加速过程燃油流量分离点A的转速为N2r0,则对积分时间常数缩放因子K作如下取值:
当N2r≤(N2r0-2%)时,K=K_rsv;
当(N2r0-2%)<N2r<N2r0时,K为K_rsv与K_need的线性插值;
当N2r=N2r0时,K=K_need;
当N2r0<N2r<(N2r0+3%)时,K从K_need线性平滑过渡到1;
当N2r0≥(N2r0+3%)时,K=1;
(6)计算得到增量式PI控制器最终所需的积分时间常数:Ti_final=Ti0/K;
(7)计算得到增量式PI控制器最终所需的比例增益Kp_final:
当N2r≤N2r0时,Kp_final=Kp0;
当N2r0<N2r<(N2r0+3%)时,Kp_final从Kp0线性平滑过渡到Kp_design;
当N2r≥(N2r0+3%)时,Kp_final=Kp_design;
(8)将Kp_final、Ti_final写入主燃油流量增量式PI控制器,进行后续台架试车。
2.如权利要求1所述的改善航空发动机加速性能的增量式PI控制器参数调整方法,其特征在于,
比例因子kKp的取值范围为0.5≤kKp≤0.8;缩放因子kTi的取值范围为3≤kTi≤5;加速过程燃油流量分离点A的转速N2r0的取值范围为92%≤N2r0≤95%。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410736305.9A CN105715384B (zh) | 2014-12-05 | 2014-12-05 | 改善航空发动机加速性能的增量式pi控制器参数调整方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410736305.9A CN105715384B (zh) | 2014-12-05 | 2014-12-05 | 改善航空发动机加速性能的增量式pi控制器参数调整方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105715384A true CN105715384A (zh) | 2016-06-29 |
CN105715384B CN105715384B (zh) | 2017-09-01 |
Family
ID=56144103
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410736305.9A Active CN105715384B (zh) | 2014-12-05 | 2014-12-05 | 改善航空发动机加速性能的增量式pi控制器参数调整方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105715384B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110850788A (zh) * | 2019-12-06 | 2020-02-28 | 中国航发控制系统研究所 | 航空发动机用多电分布式控制系统架构 |
CN111042929A (zh) * | 2019-12-09 | 2020-04-21 | 中国人民解放军总参谋部第六十研究所 | 采用油泵供油的中小型涡喷发动机转速回路控制方法 |
WO2020093264A1 (zh) * | 2018-11-07 | 2020-05-14 | 大连理工大学 | 一种航空发动机过渡态控制规律优化的设计方法 |
CN112412631A (zh) * | 2020-11-18 | 2021-02-26 | 中国航发沈阳发动机研究所 | 一种具有自刹车功能的加速过程控制方法 |
CN113357017A (zh) * | 2021-06-24 | 2021-09-07 | 中国航发沈阳发动机研究所 | 一种航空发动机加速过程转速控制方法 |
CN114545805A (zh) * | 2020-11-24 | 2022-05-27 | 中国航发商用航空发动机有限责任公司 | 用于在线调整航空发动机控制参数的方法、装置和系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5365446A (en) * | 1990-07-04 | 1994-11-15 | Aerospatiale Societe Nationale Industrielle | System for integrated pitch and thrust control of any aircraft |
RU2474712C2 (ru) * | 2010-12-29 | 2013-02-10 | Открытое акционерное общество "СТАР" | Способ управления расходом топлива в газотурбинный двигатель |
CN103592844A (zh) * | 2012-08-15 | 2014-02-19 | 王景芳 | 增量式pi参数时变智能优化控制 |
CN103967622A (zh) * | 2014-05-23 | 2014-08-06 | 中国航空动力机械研究所 | 微型燃气轮机的起动供油控制方法及系统 |
CN104022508A (zh) * | 2014-05-26 | 2014-09-03 | 中南大学 | 一种三相四开关型有源电力滤波器及其控制方法 |
-
2014
- 2014-12-05 CN CN201410736305.9A patent/CN105715384B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5365446A (en) * | 1990-07-04 | 1994-11-15 | Aerospatiale Societe Nationale Industrielle | System for integrated pitch and thrust control of any aircraft |
RU2474712C2 (ru) * | 2010-12-29 | 2013-02-10 | Открытое акционерное общество "СТАР" | Способ управления расходом топлива в газотурбинный двигатель |
CN103592844A (zh) * | 2012-08-15 | 2014-02-19 | 王景芳 | 增量式pi参数时变智能优化控制 |
CN103967622A (zh) * | 2014-05-23 | 2014-08-06 | 中国航空动力机械研究所 | 微型燃气轮机的起动供油控制方法及系统 |
CN104022508A (zh) * | 2014-05-26 | 2014-09-03 | 中南大学 | 一种三相四开关型有源电力滤波器及其控制方法 |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020093264A1 (zh) * | 2018-11-07 | 2020-05-14 | 大连理工大学 | 一种航空发动机过渡态控制规律优化的设计方法 |
US10934944B2 (en) | 2018-11-07 | 2021-03-02 | Dalian University Of Technology | Method for optimization of transient control law of aero-engine |
CN110850788A (zh) * | 2019-12-06 | 2020-02-28 | 中国航发控制系统研究所 | 航空发动机用多电分布式控制系统架构 |
CN111042929A (zh) * | 2019-12-09 | 2020-04-21 | 中国人民解放军总参谋部第六十研究所 | 采用油泵供油的中小型涡喷发动机转速回路控制方法 |
CN111042929B (zh) * | 2019-12-09 | 2022-06-24 | 中国人民解放军总参谋部第六十研究所 | 采用油泵供油的中小型涡喷发动机转速回路控制方法 |
CN112412631A (zh) * | 2020-11-18 | 2021-02-26 | 中国航发沈阳发动机研究所 | 一种具有自刹车功能的加速过程控制方法 |
CN112412631B (zh) * | 2020-11-18 | 2022-03-04 | 中国航发沈阳发动机研究所 | 一种具有自刹车功能的加速过程控制方法 |
CN114545805A (zh) * | 2020-11-24 | 2022-05-27 | 中国航发商用航空发动机有限责任公司 | 用于在线调整航空发动机控制参数的方法、装置和系统 |
CN114545805B (zh) * | 2020-11-24 | 2023-08-04 | 中国航发商用航空发动机有限责任公司 | 用于在线调整航空发动机控制参数的方法、装置和系统 |
CN113357017A (zh) * | 2021-06-24 | 2021-09-07 | 中国航发沈阳发动机研究所 | 一种航空发动机加速过程转速控制方法 |
CN113357017B (zh) * | 2021-06-24 | 2022-11-22 | 中国航发沈阳发动机研究所 | 一种航空发动机加速过程转速控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105715384B (zh) | 2017-09-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105715384A (zh) | 改善航空发动机加速性能的增量式pi控制器参数调整方法 | |
US8833085B2 (en) | System and method for gas turbine startup control | |
CN102193532B (zh) | 火电厂单元机组自启停优化控制系统 | |
CN110824926B (zh) | 一种火电机组深度调峰一次调频控制方法 | |
CN103711598B (zh) | 液压系统调节设备、方法、功率匹配控制系统和工程机械 | |
CN104865830A (zh) | 一种机组负荷双重智能优化控制方法 | |
CA2943023C (en) | Feedforward control with intermittent re-initialization based on estimated state information | |
EP3161310B1 (en) | Wind turbine controller with pitch feedback control loop in partial load | |
CN105934862B (zh) | 具有改进的额定值跃变特性的风电场调节装置 | |
US11279366B1 (en) | Feedforward mechanism with signal decay for torque adjustment in diesel engine operation | |
DE102016113660A1 (de) | Verfahren und Vorrichtung zur Optimierung einer Vermischung mittels Dampfsammler und Aufladung mittels Gasturbine in Kombikraftwerken | |
CN110513158B (zh) | 汽轮机前馈的多级速率调节方法 | |
CN104712378A (zh) | 火电机组主蒸汽压力闭环节能控制方法和系统 | |
CN104696077A (zh) | 燃料流量调节方法和pid控制器 | |
CN103236715B (zh) | 一种发电机组调速器机械开限的调整方法 | |
CN102828944B (zh) | 工程机械及其泵流量控制系统和方法 | |
CN110081410A (zh) | 一种汽动给水泵小汽轮机的控制方法 | |
CN115808876A (zh) | 发动机尾喷管执行机构自适应控制方法及装置 | |
CN113482733B (zh) | 汽轮机进汽调节阀的顺序阀直接配汽方法、装置和汽轮机 | |
US10647413B2 (en) | Enhanced engine load demand anticipation | |
CN103670539A (zh) | 补偿发电机组动态特性联合调频控制方法、系统与装置 | |
JP5091712B2 (ja) | ガバナフリー制御装置およびガバナフリー制御方法 | |
CN102817394A (zh) | 一种挖掘机液压泵控制系统及方法及挖掘机 | |
US20100298996A1 (en) | Method for operating a power station | |
JP4905835B2 (ja) | 火力発電プラントの負荷制御装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP03 | Change of name, title or address | ||
CP03 | Change of name, title or address |
Address after: No. 104, Liangxi Road, Binhu District, Wuxi City, Jiangsu Province Patentee after: AECC AERO ENGINE CONTROL SYSTEM INSTITUTE Address before: 214063 Jiangsu city of Wuxi province Liangxi Road No. 792 Patentee before: AVIC AVIATION MOTOR CONTROL SYSTEM INSTITUTE |