劉愛元，劉持勝，戴洪德

（海軍航空工程學(xué)院a.控制工程系；b.外訓(xùn)系，山東煙臺264001）

飛機上實現(xiàn)交流電源系統(tǒng)輸出頻率恒定通常采用的方法是加裝頻率調(diào)節(jié)環(huán)路和恒速傳動裝置。然而，頻率調(diào)節(jié)環(huán)路和恒速傳動裝置在工作過程中因參數(shù)的變化，將直接影響交流電源系統(tǒng)的動靜態(tài)特性，甚至造成系統(tǒng)的不穩(wěn)定，且這種影響會隨著飛機發(fā)動機工作狀態(tài)的變化而變化[1-2]。

本文基于以上考慮，按照文獻[3-7]的研究方法，結(jié)合某型恒頻交流電源頻率調(diào)節(jié)環(huán)路的工作原理，在建立系統(tǒng)數(shù)學(xué)模型的基礎(chǔ)上，用MATLAB語言[8]對某型飛機發(fā)動機在中轉(zhuǎn)速負載狀態(tài)下，通過改變系統(tǒng)濾波整形放大電路第一級濾波器放大倍數(shù)k1和頻率敏感電壓變換電路傳導(dǎo)系數(shù)k5等參數(shù)進行仿真，得出其參數(shù)變化對系統(tǒng)動靜態(tài)特性所造成的影響，其仿真結(jié)果及所得結(jié)論既可為進行該系統(tǒng)研究的人員提供參考，又為某型飛機該系統(tǒng)的技術(shù)保障提供了有效的技術(shù)支持。

1 系統(tǒng)組成及工作原理

頻率調(diào)節(jié)環(huán)路主要由頻率敏感和頻率電壓變換電路、基準(zhǔn)信號電路、兩級一階有源濾波整形放大電路、PID 校正電路、電壓比較電路、脈寬調(diào)制與功率放大電路、電磁恒裝等構(gòu)成，系統(tǒng)的組成框圖如圖1 所示。[9-10]

圖1 頻率調(diào)節(jié)環(huán)路組成框圖

系統(tǒng)工作原理是：當(dāng)變速恒頻交流發(fā)電系統(tǒng)因受到某種干擾而導(dǎo)致系統(tǒng)頻率偏離穩(wěn)定值時，頻率敏感電路即敏感此頻率的變化，并經(jīng)頻率電壓變換電路轉(zhuǎn)換成相應(yīng)的電壓信號U2，該電壓信號U2與來自基準(zhǔn)信號電路的電壓U基進行比較，產(chǎn)生的偏差信號ΔU經(jīng)兩級一階有源濾波和整形放大電路，濾除其高次諧波，獲得幅值為0.5 V、頻率為400 Hz 的三角波，此三角波經(jīng)PID校正電路和電壓比較電路送入脈寬調(diào)制和功率放大電路，獲得脈沖方波信號。

此脈沖方波在不同的狀態(tài)下具有不同的導(dǎo)通截止比，這種導(dǎo)通截止比控制功率放大電路，以獲得相應(yīng)的導(dǎo)通截止比，進而控制電磁恒裝中激磁電流的大小，使其同步轉(zhuǎn)速n穩(wěn)定在一定精度范圍內(nèi)。例如，當(dāng)干擾使交流發(fā)電機的輸出頻率降低時，兩級一階有源濾波的反相輸入電位降低，PID 校正電路的輸出信號升高，電壓比較器的輸出降低，電磁恒裝中激磁電流增大，使其同步轉(zhuǎn)速n上升，輸出頻率升高，將其頻率調(diào)整到規(guī)定的范圍內(nèi)。

2 頻率調(diào)節(jié)環(huán)路數(shù)學(xué)模型的建立

1）頻率敏感和電壓變換電路。頻率敏感和電壓變換電路相當(dāng)于從頻率到電壓的一個變換器，在數(shù)學(xué)上可用一個比例系數(shù)來描述。當(dāng)電網(wǎng)頻率f=400 Hz，電壓U=0.96 V 時，該變換電路的比例系數(shù)為：

2）濾波整形放大電路。濾波整形放大電路采用了兩級一階有源濾波器，可使濾波時間常數(shù)下降到幾乎不影響系統(tǒng)的動態(tài)品質(zhì)，又可產(chǎn)生一個在形狀和幅值上能滿足后級使用的三角波。本級參數(shù)選擇應(yīng)遵循：輸出幅值為0.5 V而頻率為400 Hz的調(diào)制三角波，其時間常數(shù)盡可能小，盡可能提高其放大倍數(shù)。為此選定：

第一級濾波器：時間常數(shù)T1=6.8 ms，放大倍數(shù)k1=8.5；第二級濾波器：時間常數(shù)T2=4.7 ms，放大倍數(shù)k2=3；因一階濾波器相當(dāng)于一階慣性環(huán)節(jié)，其數(shù)學(xué)模型為K(1+Ts)。

故本級的數(shù)學(xué)模型是：

3）PID校正電路。根據(jù)系統(tǒng)固有的開環(huán)幅相頻特性，在作出期望的理想頻率特性的基礎(chǔ)上，減去其固有特性，即得出校正模塊的對數(shù)相頻特性，據(jù)此可得本級的數(shù)學(xué)模型為：

τ1=0.034 s，τ2=0.31 s，Tr1=1.2 s，Tr2=0.004 85 s。

4）脈寬調(diào)整與功率放大電路。脈寬調(diào)整與功率放大電路可近似看作一個比例環(huán)節(jié)，其比例系數(shù)為：

5）電磁恒裝。設(shè)電磁恒裝在最大力矩時對應(yīng)的滑速差為nsm，當(dāng)滑速差ns＜nsm時，系統(tǒng)穩(wěn)定，其對數(shù)幅相頻特性一一對應(yīng)；當(dāng)ns＞nsm時，系統(tǒng)不穩(wěn)定，其對數(shù)幅相頻特性不對應(yīng)。但考慮ns＞nsm恰是高空常出現(xiàn)的一種正常狀態(tài)，故本級數(shù)學(xué)模型為

式中：Kbr、Td、TM均隨發(fā)動機轉(zhuǎn)速和負載的變化而變化。

6）轉(zhuǎn)速頻率比例系數(shù)。在誤差允許范圍內(nèi)電磁恒裝輸出的轉(zhuǎn)速經(jīng)交流發(fā)電機轉(zhuǎn)換為頻率，其轉(zhuǎn)速頻率比例系數(shù)為

分析式（1）～（6）并根據(jù)系統(tǒng)工作原理及各個環(huán)節(jié)之間的信號傳遞關(guān)系，得出頻率調(diào)節(jié)環(huán)路總的數(shù)學(xué)模型如圖2所示。

圖2 頻率調(diào)節(jié)環(huán)路結(jié)構(gòu)圖

3 變參數(shù)仿真與動靜態(tài)特性分析

3.1 變參數(shù)仿真

在建立系統(tǒng)數(shù)學(xué)模型的基礎(chǔ)上，本文采用MATLAB 語言的矩陣運算、圖形繪制、信號處理和仿真等功能，對某型發(fā)動機的中轉(zhuǎn)速中負載狀態(tài)下頻率調(diào)節(jié)環(huán)路的參數(shù)k1、k2、k3、k4、k5的變化進行仿真研究。某型發(fā)動機在中轉(zhuǎn)速、中負載狀態(tài)下頻率調(diào)節(jié)環(huán)路的結(jié)構(gòu)圖如圖3所示。

圖3 中轉(zhuǎn)速中負載系統(tǒng)結(jié)構(gòu)圖

圖（3）中G(s)為：

由控制理論可知，系統(tǒng)的開環(huán)增益k1、k2、k3、k4對系統(tǒng)的影響實質(zhì)是一致的，只因數(shù)值的不同而相差一個比例系數(shù)，影響范圍有所不同而已。而系統(tǒng)的反饋系數(shù)k5的變化對系統(tǒng)動態(tài)特性的影響與k、k2、k3、k4對系統(tǒng)的影響完全不同。本文按兼顧全面、突出重點的原則，選擇具有代表性的濾波整形放大電路第一級濾波器放大倍數(shù)k1和頻率敏感電壓變換傳導(dǎo)系數(shù)k5進行變參數(shù)仿真。

1）變?yōu)V波整形放大電路第1 級濾波器放大倍數(shù)k1。取k2=3、k3=33.5、k4=0.1、k5=0.002 5，選擇不同的k1值，通過MATLAB 中Simulink 工具，得出系統(tǒng)階躍響應(yīng)曲線如圖4所示[5-7]。

圖4 k1 變化時系統(tǒng)階躍響應(yīng)曲線

2）變頻率敏感電壓變換傳導(dǎo)系數(shù)k5。取k1=8.5、k2=3、k3=33.5、k4=0.1，選擇不同的k5值，通過MATLAB中Simulink工具，得出系統(tǒng)階躍響應(yīng)曲線如圖5所示。

圖5 k5 變化時系統(tǒng)階躍響應(yīng)曲線

3.2 動靜態(tài)特性分析

對該系統(tǒng)動靜態(tài)特性的要求是（只討論時域部分，頻域部分略）：超調(diào)量σ≤17.5%、調(diào)節(jié)時間為ts≤15 s、峰值時間tp≤2 s。

1）變k1的動靜態(tài)特性分析。根據(jù)圖4 中變k1仿真的各響應(yīng)曲線，結(jié)合系統(tǒng)動靜態(tài)指標(biāo)要求，得出其仿真結(jié)果如表1所示。

表1 中轉(zhuǎn)速中負載狀態(tài)下變k1 仿真結(jié)果

由表1 可知，在0.434 ＜k1＜37.89 范圍內(nèi)，系統(tǒng)是穩(wěn)定的。當(dāng)k1由0.434向37.89逐漸增大時，峰值時間tp逐漸減小而調(diào)節(jié)時間ts先減小后稍微增大；超調(diào)量先減小后又急劇增大，k1在3.5～13.8之間時σ較小，當(dāng)k1≥13.8 時，系統(tǒng)開始出現(xiàn)震蕩，并逐步加劇，直至不穩(wěn)定。

2）變k5的動靜態(tài)特性分析。根據(jù)圖5，變k5仿真的各響應(yīng)曲線，結(jié)合系統(tǒng)動靜態(tài)指標(biāo)要求，得出其仿真結(jié)果如表2所示。

表2 中轉(zhuǎn)速中負載狀態(tài)下變k5 仿真結(jié)果

由表2 可知，在0.000 122 5 ＜k5＜0.010 7 范圍內(nèi)，系統(tǒng)是穩(wěn)定的。當(dāng)k5由0.000 122 5向0.010 7逐漸增大時，峰值時間tp逐漸減小而調(diào)節(jié)時間ts先減小后逐漸增大；超調(diào)量σ先減小后又急劇增大。

4 結(jié)論

k1適當(dāng)增大可提高系統(tǒng)的反應(yīng)速度，但k1過大則系統(tǒng)不穩(wěn)定。k1在7附近取值時，系統(tǒng)具有最佳的動靜態(tài)特性指標(biāo)。

系統(tǒng)對k5的變化異常敏感，且k5的微小變化對系統(tǒng)的動靜態(tài)特性的影響十分明顯。因系統(tǒng)參數(shù)稍有變化，經(jīng)過數(shù)級整形放大后其影響非常大，甚至導(dǎo)致系統(tǒng)不穩(wěn)定，因而在設(shè)計和使用過程中，應(yīng)對該部分電路加以注意。設(shè)計時應(yīng)盡量采取措施使k5穩(wěn)定在0.002附近，而在使用過程中應(yīng)注意檢測該部分電路的完好情況，以防器件損壞等影響系統(tǒng)的動靜態(tài)特性。

k5在0.002附近時具有最佳的動靜態(tài)特性。

[1] 劉愛元，盧建華，楊秀珍，等.飛機電源與電氣控制[M].北京：海潮出版社，2008：106-141.

LIU AIYUAN，LU JIANHUA，YANG XIUZHEN，et al.Aircraft power supply and electric control[M]. Beijing：Haichao Press，2008：106-141.（in Chinese）

[2] 劉慧英，姚遠，張國春.飛機電源系統(tǒng)保護控制的仿真研究[J].控制工程，2010，17(2)：135-138.

LIU HUIYING，YAO YUAN，ZHANG GUOCHUN.Simulation of protection control of aircraft power supply system[J]. Scientific Journal of Control Engineering，2010，17(2)：135-138.（in Chinese）

[3] 馬印才，王春平，賈占強. 二階PID 電路狀態(tài)觀測器MATLAB 實現(xiàn)[J].科學(xué)技術(shù)與工程，2009，7(19)：5089-5090.

MA YINCAI，WANG CHUNPING，JIA ZHANQIANG.Accomplish of state observer for one of two order pid circuit by MATLAB[J]. Science Technology and Engineering，2009，7(19)：5089-5090.（in Chinese）

[4] 洪陸合，吳順祥.基于Matlab/GUI 的逆變電路故障仿真研究[J]. 廈門大學(xué)學(xué)報：自然科學(xué)版，2010，49(3)：324-327.

HONG LUHE，WU SHUNXIANG. Research on simulation for inverter circuit under fault conditions based on Matlab/GUI[J]. Journal of Xiamen University：Natural Science，2010，49(3)：324-327.（in Chinese）

[5] 韋莉，張逸成，沈玉琢，等.基于MATLAB的大功率電源系統(tǒng)的建模與仿真[J].系統(tǒng)仿真技術(shù)，2009，5(3)：202-207.

WEI LI，ZHANG YICHENG，SHEN YUZHUO，et al.Modelling and simulation of high power system based on MATLAB[J].System Simulation Technology，2009，5(3)：202-207.（in Chinese）

[6] 劉建英.基于MATLAB的飛機電源系統(tǒng)仿真[J].中國民航學(xué)院學(xué)報，2004，22(6)：91-92.

LIU JIANYING.Simulation for aircraft power supply system based MATLAB[J]. Journal of Civil Aviation university of China，2004，22(6)：91-92.（in Chinese）

[7] 高朝暉，張曉斌，吳小華.基于MATLAB 的飛機高壓直流電源系統(tǒng)數(shù)字仿真[J]. 系統(tǒng)仿真學(xué)報，2003，15(4)：564-567.

GAO ZHAOHUI，ZHANG XIAOBIN，WU XIAOHUA.Simulation of high voltage direct current aircraft electrical power generating system based on MATLAB[J].Journal of System Simulation，2003，15(4)：564-567.（in Chinese）

[8] 薛定宇. 基于Matlab/Simulink 的系統(tǒng)仿真技術(shù)與應(yīng)用[M].北京：清華大學(xué)出版社，2002：192-220.

XUE DINGYU. Technology and application for system simulation based Matlab/Simulink[M]. Beijing：Tsinghua University Press，2002：192-220.（in Chinese）

[9] 劉愛元.恒頻交流電源頻率調(diào)節(jié)系統(tǒng)動靜態(tài)特性研究[J].計算機仿真，2003，20(6)：108-110.

LIU AIYUAN. Research on dynamic and static characteristic of frequency modulation system for constant frequency AC power supply[J].Computer Simulation，2003，20(6)：108-110.（in Chinese）

[10]沈頌華.航空航天器供電系統(tǒng)[M].北京：北京航空航天大學(xué)出版社，2005：2-20.

SHEN SONGHUA. Power supply system for aeronautics and astronautics[M]. Beijing：Beihang University Press，2005：2-20.（in Chinese）