李 帥,楊寶玉,陸 燕
一種基于空間光機熱模型的自適應(yīng)PID控制方法
李 帥1,2,楊寶玉1,陸 燕1
(1. 中國科學(xué)院上海技術(shù)物理研究所,上海 200083;2. 中國科學(xué)院大學(xué),北京 100049)
為提高空間相機的控溫穩(wěn)定度以保證成像質(zhì)量,本文提出一種基于空間光機熱模型的自適應(yīng)比例積分微分(proportional-integral-derivative, PID)控制方法。該控制器的設(shè)計從空間光機的熱平衡方程出發(fā),能夠?qū)崟r根據(jù)光機及與其輻射換熱對象的溫度修正光機的熱模型,繼而采用極點配置的方法實時校正PID控制器參數(shù),最終確定本控溫周期的加熱占空比。本文通過建立抽象的空間光機熱模型,分別施加上述自適應(yīng)PID控制方法與固定參數(shù)PID控制方法,對控溫效果進行了仿真及實驗對比。結(jié)果表明,對環(huán)境擾動引起的溫度波動,該自適應(yīng)PID控制器始終保持最佳動態(tài)響應(yīng),控溫穩(wěn)定度優(yōu)于±0.1K,具有更好的控溫穩(wěn)定性和環(huán)境適應(yīng)性。
空間光機;熱模型;極點配置;自適應(yīng);PID控制
空間相機是實現(xiàn)對地觀測、大氣和海洋探測及宇宙探測等應(yīng)用的主要有效載荷,根據(jù)不同的探測要求,相機精度和成像質(zhì)量也有不同要求。對高精度相機其光機溫度對成像質(zhì)量的影響受到廣泛關(guān)注,當(dāng)溫度發(fā)生變化時,一方面鏡片表面的曲率半徑、透鏡的厚度、光學(xué)材料的折射率均會發(fā)生改變[1],另一方面由于鏡框的熱膨脹系數(shù)較大,鏡框會相對于鏡片發(fā)生差分膨脹或收縮[2],這均會造成焦距偏離和系統(tǒng)失準(zhǔn),導(dǎo)致相機成像模糊,嚴(yán)重時甚至出現(xiàn)成像失敗的情況。
空間光學(xué)遙感相機在軌工作時,由于太陽輻射、地球紅外輻射、地球反照、空間背景輻射的交替變化,造成相機的光機周圍結(jié)構(gòu)溫度發(fā)生劇烈變化,周圍結(jié)構(gòu)通過輻射換熱影響光機溫度也發(fā)生波動。為使相機的光機維持在穩(wěn)定的溫度水平以保障成像質(zhì)量,除采用被動熱控措施外,一般光機還采用直接或間接的主動控溫加熱回路,通過控制加熱器的加熱功率,精確補償損失的熱量,實現(xiàn)光機溫度的穩(wěn)定控制。相機在具有寬視場及高分辨率的要求場合,光機對控溫穩(wěn)定度提出了更高的要求。例如“高分二號”衛(wèi)星,其上兩個相機的主鏡、次鏡、次鏡支撐件、前鏡筒等光機結(jié)構(gòu)都要求在控溫精度為20℃±2℃時,控溫穩(wěn)定度為±0.3℃/h[3];空間太陽望遠(yuǎn)鏡(space solar telescope,SST)的主鏡鏡面在等直徑的環(huán)上控溫穩(wěn)定度為±0.1℃/h[4];Herschel衛(wèi)星的遠(yuǎn)紅外光學(xué)系統(tǒng)要求10s內(nèi)控溫穩(wěn)定度≤3×10-4K[5]。
為提高控溫穩(wěn)定度,除了需要對光機及其周圍環(huán)境進行合理的熱設(shè)計,設(shè)計多個主動控溫回路以及采用多級控溫策略[6-8]外,改進溫度控制算法對提高控溫穩(wěn)定度的效果更為突出。目前,空間光學(xué)遙感相機上大多采用開關(guān)比例控制的控溫方式,一些重要結(jié)構(gòu)部件則采用控溫穩(wěn)定度更高的固定參數(shù)PID控制。固定參數(shù)PID控制具有控制算法簡單、技術(shù)成熟、可靠性高,為大多數(shù)工程師所熟悉的優(yōu)點,能夠滿足大多數(shù)被控對象的控溫要求[9]。但其僅依據(jù)被控對象控溫點的溫度做控制決策,不論與被控對象存在熱關(guān)聯(lián)的周圍熱環(huán)境如何改變都采用不變的PID控制參數(shù),這可能給被控對象的溫度帶來不必要的波動。雖然有些文獻(xiàn)中指出控制的動態(tài)和穩(wěn)態(tài)特性取決于所設(shè)置的PID參數(shù)[10-12],并由此發(fā)展出模糊自適應(yīng)PID控制、專家PID控制、基于模型辨識的自適應(yīng)PID控制等智能控制方法。模糊自適應(yīng)PID控制、專家PID控制等基于規(guī)則的PID參數(shù)自整定方法簡單,在線控制容易實現(xiàn),但整定出的PID參數(shù)并不是該控溫系統(tǒng)在該時刻的最優(yōu)參數(shù);基于模型辨識的自適應(yīng)PID控制的在線辨識較為復(fù)雜,計算機計算量大、工作時間長,而且還需要解決閉環(huán)可辨識性問題。因此,這些智能控溫方法尚未在空間相機的工程實踐得到應(yīng)用。目前,國外在空間相機的溫度控制中采取了基于被控對象熱特性的自適應(yīng)PID控溫算法,提高了系統(tǒng)控溫穩(wěn)定度。例如,“普朗克”衛(wèi)星光學(xué)系統(tǒng)的控溫算法是最優(yōu)PI(proportional-integral)控溫算法,該算法從被控對象的解析模型出發(fā),建立了一組可供探索的參數(shù)列表,以確定空間環(huán)境下儀器運行的最優(yōu)傳遞函數(shù),然后利用該溫控系統(tǒng)的傳遞函數(shù)模型實現(xiàn)PI控制參數(shù)的最優(yōu)整定[13],實現(xiàn)控溫穩(wěn)定度優(yōu)于±0.1K;“Herschel”衛(wèi)星的遠(yuǎn)紅外光學(xué)系統(tǒng)采用自適應(yīng)PI控溫算法,該算法通過熱平衡試驗獲取的溫度數(shù)據(jù)修正熱分析模型,根據(jù)修正后的熱分析模型,確定各個部件的傳遞函數(shù),從滿足其溫度穩(wěn)定性要求出發(fā),重新整定PI控制器參數(shù)[5],實現(xiàn)控溫穩(wěn)定度優(yōu)于±0.5K。我國在空間相機的溫度控制領(lǐng)域,暫未出現(xiàn)基于被控對象熱特性的自適應(yīng)PID控制方法應(yīng)用案例,PID控制器參數(shù)未能結(jié)合具體控溫對象進行設(shè)計,導(dǎo)致系統(tǒng)難以實現(xiàn)高穩(wěn)定度的溫度控制。
本文從空間光機抽象模型的熱平衡方程出發(fā),提出了一種能根據(jù)光機及與其輻射換熱對象的溫度實時修正光機數(shù)學(xué)模型,繼而采用極點配置的方法實時校正控制器參數(shù)的自適應(yīng)PID控制方法,利用建立的空間光機抽象模型對該控制方法的控溫效果進行了仿真及實驗驗證。
空間光機的主動加熱控溫回路一般包括溫度控制器、測溫用熱敏電阻(或熱電偶、測溫二極管等)、加熱器等。圖1為本文提出的基于空間光機熱模型的自適應(yīng)PID控溫原理圖。該溫度控制器需要同時采集被控對象及與其存在輻射換熱的周圍環(huán)境的溫度,通過計算出等效輻射換熱系數(shù)來修正被控對象模型,然后以二階最佳動態(tài)響應(yīng)模型為目標(biāo),利用極點配置設(shè)計PID參數(shù),輸出加熱占空比,最終通過加熱器實現(xiàn)被控對象的高穩(wěn)定度溫度控制。該自適應(yīng)PID控制方法具體設(shè)計過程詳細(xì)介紹如下。
圖1 基于空間光機熱模型的自適應(yīng)PID控溫原理圖
對于某空間光機的控溫系統(tǒng),設(shè)其控溫點的溫度為m,與該控溫區(qū)域存在傳熱的周圍環(huán)境的溫度分別為1,2,…,T,空間相機處于真空環(huán)境,只存在熱傳導(dǎo)和熱輻射兩種傳熱方式,忽略單個測溫區(qū)域的溫度分布不均勻性,建立光機的熱平衡方程:
式中:p為加熱器功率;為光機的總熱容;(其中,1,2,···,)為光機與換熱對象之間的等效換熱系數(shù)。
對于導(dǎo)熱換熱,一般視為不隨溫度變化的常數(shù),甚至由于光機的安裝一般使用接觸面積小、導(dǎo)熱率低的鈦合金螺釘和隔熱墊安裝,導(dǎo)熱換熱的等效換熱系數(shù)i可以忽略不計,而對于輻射換熱,可由(2)式進行計算:
式中:為斯忒藩-玻爾茲曼常量,其值為5.67×10-8W/(m2·K);i為光機與輻射換熱對象i之間的系統(tǒng)發(fā)射率;i為光機與輻射換熱對象i之間的輻射換熱面積。
對式(1)兩邊做拉普拉斯變換,可得光學(xué)對象溫度m與加熱器功率p和其他換熱對象溫度i之間的關(guān)系如下:
由于加熱器功率p的變化反映到光機溫度的變化存在一定的滯后時間,所以加熱器功率p和光機溫度m之間的傳遞函數(shù)p()為:
由式(4)可得p()帶零階保持器的廣義對象脈沖傳遞函數(shù)為:
由式(2)可知,當(dāng)光機及與其輻射換熱對象的溫度變化時,隨之發(fā)生變化,加熱器功率p和光機溫度m之間的傳遞函數(shù)p(-1)即發(fā)生改變。即能夠通過跟蹤光機及與其輻射換熱對象的溫度對光機的熱模型進行實時修正。
1)極點配置控制器算法
設(shè)被控對象受控自回歸滑動平均(controlled auto regressive moving average, CARMA)模型為:
式中:()、()、()分別為被控對象的輸出、輸入和不可檢測的零均值白噪聲序列;-1為后移算子;為滯后步數(shù);(-1)、(-1)、(-1)均為-1的多項式。
根據(jù)式(5)可得出:
常規(guī)的計算機閉環(huán)控制系統(tǒng)的結(jié)構(gòu)可表示為圖2。
圖2 計算機閉環(huán)控制系統(tǒng)結(jié)構(gòu)圖
由圖2可得,該系統(tǒng)的閉環(huán)方程為:
從式(10)中可知,根據(jù)增廣型自校正閉環(huán)極點配置[14]的要求,閉環(huán)特征多項式為:
式中:(-1)為期望特征多項式。
進行極點配置時,一般以典型的二階系統(tǒng)閉環(huán)傳遞函數(shù)的標(biāo)準(zhǔn)形式作為目標(biāo)。
與式(12)相對應(yīng)的離散特征多項式為:
式中:為阻尼比;n為無阻尼自然振蕩角頻率。
對式(13),當(dāng)二階系統(tǒng)最佳阻尼比=0.707時,為二階最佳動態(tài)模型,在單位階躍作用下的超調(diào)量%=4.3%,相角穩(wěn)定裕量(c)=65.5°。
采樣周期和n、的關(guān)系可按下式計算:
圖2系統(tǒng)可直接寫出極點配置自校正控制器輸出()的表達(dá)式為:
2)增量型PID控制算法
在數(shù)字系統(tǒng)中,通常采用帶有數(shù)字濾波器的PID控制器算法[15]:
式中:P為比例系數(shù);I為積分系數(shù);D為微分系數(shù)。
方程組(14)對P、I、D有唯一解,即:
3)極點配置自適應(yīng)PID控制算法
為將極點配置控制器變換為增量型PID控制器,對照式(16),式(15)選用下述形式:
由式(15)可知,為保證系統(tǒng)穩(wěn)態(tài)誤差最終為零,還需選擇:
式(11)中選擇多項式(-1)時應(yīng)防止在系統(tǒng)輸出端出現(xiàn)突變而導(dǎo)致振蕩[16],因此考慮選擇:
采用極點配置法整定PID控制器參數(shù),即是選取PID參數(shù),使系統(tǒng)具有期望的閉環(huán)特征方程。所以將式(8)、(9)、(13)、(19)、(20)、(22)代入式(11),令方程式(11)兩邊系數(shù)相等,解代數(shù)方程即可求得0、1、2,最后根據(jù)式(18)即可求得以二階最佳動態(tài)模型為目標(biāo)的PID參數(shù)。
該控溫方法與固定參數(shù)的PID控溫方法相比,當(dāng)光機及與其輻射換熱對象的溫度變化時,該控溫方法能夠?qū)鈾C的熱模型進行實時修正,然后采用極點配置的方法以二階最佳動態(tài)響應(yīng)模型為目標(biāo)實時校正PID控制器參數(shù),使PID控制器的動態(tài)控溫性能始終保持最佳狀態(tài),因而有助于減小光機由于熱慣性造成的不必要的溫度波動,具有更好控溫的穩(wěn)定性。
為了對該控溫方法進行仿真及實驗驗證,本文建立光機的熱控系統(tǒng)模型如圖3所示。
圖3 光機熱控系統(tǒng)模型
建立被控對象(鋁塊)的數(shù)學(xué)模型,取滯后時間=2s,由上文中分析可得:
加熱器功率p與鋁塊溫度m之間的傳遞函數(shù)p()為:
環(huán)境溫度s與鋁塊溫度m之間的傳遞函數(shù)s()為:
式中:
取采樣時間==2s,由上述1.1節(jié)的分析可得p()帶零階保持器的廣義對象脈沖傳遞函數(shù)為:
由上述1.2節(jié)中分析可得PID參數(shù),
上述熱控系統(tǒng)的仿真模型借助MATLAB軟件編程建立,該自適應(yīng)PID控溫方法與固定參數(shù)PID控溫方法的仿真結(jié)果如圖4所示。對于不同的環(huán)境溫度,PID控制器的最佳動態(tài)響應(yīng)參數(shù)不同,圖中兩個固定參數(shù)PID控制器的PID參數(shù)為控溫過程中某一溫度下被控系統(tǒng)模型對應(yīng)的最優(yōu)PID參數(shù),PID控制1的控制器參數(shù)為:P=0.412,I=0.050,D=0.007,PID控制2的控制器參數(shù)為:P=0.447,I=0.045,D=0.006。
圖4 仿真結(jié)果對比
圖4中0時刻為環(huán)境溫度保持298K時,控制器的階躍響應(yīng),結(jié)果顯示在鋁板升溫并穩(wěn)定在(333.15±0.1)K過程中,PID控制1所需時間為135s,超調(diào)量為0;PID控制2所需時間為132s,超調(diào)量為2.15K;自適應(yīng)PID控制所需時間為108s,超調(diào)量為0.55K。結(jié)果表明隨著被控對象溫度的改變,該PID控制器始終保持最佳動態(tài)響應(yīng),收斂更快。
圖4中在300s時刻,使環(huán)境溫度快速上升至323K,模擬空間環(huán)境外熱流的大幅值增大,結(jié)果顯示鋁塊溫度再次恢復(fù)穩(wěn)定到(333.15±0.1)K過程中,PID控制1超調(diào)量為3.55K;PID控制2超調(diào)量為3.15K;自適應(yīng)PID控制超調(diào)量為2.85K。結(jié)果表明隨著空間環(huán)境外熱流大幅值增大導(dǎo)致環(huán)境溫度快速升高,控溫系統(tǒng)的數(shù)學(xué)模型發(fā)生改變時,該PID控制器始終保持最佳動態(tài)響應(yīng),超調(diào)量小,調(diào)節(jié)時間短,抗干擾性更強,控溫穩(wěn)定性更好。
將上述光機熱控系統(tǒng)模型放置于真空罐內(nèi)進行實驗,該自適應(yīng)PID控溫方法與固定參數(shù)PID控溫方法的實驗結(jié)果如圖5所示,圖中固定PID控制器參數(shù)的選取與上文中仿真的控制器參數(shù)選取一致。
圖5 實驗結(jié)果對比
圖5中在0時刻,保持正方體殼溫度為298K不變,開啟鋁板控溫加熱,使其控溫穩(wěn)定在(333.15±0.1)K過程中,PID控制1所需時間為1684s,超調(diào)量為0.42K;PID控制2所需時間為1642s,超調(diào)量為1.06K;自適應(yīng)PID控制所需時間小于1404s,超調(diào)量小于0.54K,控溫穩(wěn)定度優(yōu)于±0.1K。結(jié)果表明隨著被控對象溫度的改變,該自適應(yīng)PID控制具有更好的收斂性。
圖5中在2100s時刻,鋁板穩(wěn)定控溫在333.15K時,對正方體殼六個面分別同時加熱18W,模擬空間環(huán)境外熱流的大幅值增大,使正方體殼六個面的溫度從298K升高至323K。鋁塊溫度再次恢復(fù)穩(wěn)定到(333.15±0.1)K過程中,PID控制1超調(diào)量為1.34K;PID控制2超調(diào)量為1.12K;自適應(yīng)PID控制超調(diào)量為0.94K,控溫穩(wěn)定度優(yōu)于±0.1K。結(jié)果表明隨著空間環(huán)境外熱流大幅值增大導(dǎo)致環(huán)境溫度快速升高,控溫系統(tǒng)的數(shù)學(xué)模型發(fā)生改變時,該PID控制器始終保持最佳動態(tài)響應(yīng),超調(diào)量小,抗干擾性更強,控溫穩(wěn)定性更好。
本文從空間光機的熱平衡方程出發(fā),提出一種基于空間光機熱模型的極點配置自適應(yīng)PID控制方法。該方法揭示了空間光機的熱模型與周圍輻射換熱對象的溫度之間的聯(lián)系,在對空間光機進行控溫時,能夠?qū)崟r根據(jù)采集的光機及與其輻射換熱對象的溫度修正光機的熱模型,繼而采用極點配置的方法實時校正PID控制器參數(shù),最終通過加熱器實現(xiàn)被控對象的高穩(wěn)定度溫度控制,文中通過仿真及實驗驗證,在空間環(huán)境外熱流大幅值增大導(dǎo)致環(huán)境溫度快速升高,控溫系統(tǒng)的數(shù)學(xué)模型發(fā)生改變時,該自適應(yīng)PID控制器始終保持最佳動態(tài)響應(yīng),超調(diào)量小,調(diào)節(jié)時間短,控溫穩(wěn)定度優(yōu)于±0.1K,具有更好的控溫穩(wěn)定性和環(huán)境適應(yīng)性。
[1] 付家鑫. 多波段共孔徑光學(xué)系統(tǒng)結(jié)構(gòu)設(shè)計與分析[D]. 長春: 長春理工大學(xué), 2014.
FU Jiaxin. Structure Design and Analysis of Multiband and Total Aperture Optical System[D]. Changchun: Changchun University of Science and Technology, 2014.
[2] 鄭建麗. 低溫光學(xué)溫度測控技術(shù)研究[D]. 北京: 中國科學(xué)院大學(xué), 2011.
ZHENG Jianli. Research on Low Temperature Optical Temperature Measurement and Control Technology[D]. Beijing: University of Chinese Academy of Sciences, 2011.
[3] 趙振明, 魯盼, 宋欣陽. “高分二號”衛(wèi)星相機熱控系統(tǒng)的設(shè)計與驗證[J]. 航天返回與遙感, 2015, 36(4): 34-40.
ZHAO Zhenming, LU Pan, SONG Xinyang. Thermal design and test for high resolution space camera on GF-2 satellite[J]., 2015, 36(4): 34-40.
[4] 李蓉. 空間太陽望遠(yuǎn)鏡主光學(xué)望遠(yuǎn)鏡熱效應(yīng)分析[J]. 航天返回與遙感, 2001, 31(12): 32-39.
LI Rong. Thermal effect analysis of main structure in space solar telescope[J]., 2001, 31(12): 32-39.
[5] Savino De Palo, Michele Cairola, Marco Compassi, et al. Herschel heaters control modeling and correlation[C]//, 2011, 4(1): 29-39.
[6] 胡旭曉, 楊克己, 臺憲青. 逐級遞推式超精密溫度控制策略研究[J]. 機械工程學(xué)報, 2004, 40(2): 64-67.
HU Xuxiao, YANG Keji, TAI Xianqing. Research on step by step recursive ultra-precision temperature control strategy[J]., 2004, 40(2): 64-67.
[7] Taeho Keem, Seung-Woo Kim. Thermal stabilization for accurate dimensional measurement using Gallium[J]., 2004, 44: 701-706.
[8] 李國強, 耿利寅, 童葉龍. 航天器鉚鐘的一種精密控溫系統(tǒng)[J]. 航天器工程, 2011, 20(4): 93-98.
LI Guoqiang, GENG Liyin, Tong Yelong. A precise temperature control system for spacecraft rubidium atomic clock[J]., 2011, 20(4): 93-98.
[9] 孟曉倩, 楊寶玉, 張家昆, 等. 機械制冷機溫控的系統(tǒng)仿真與分析[J]. 低溫與超導(dǎo), 2019, 47(5): 5-10.
MENG Xiaoqian, YANG Baoyu, ZHANG Jiakun, et al. System simulation and analysis of temperature control of mechanical cryocooler [J]., 2019, 47(5): 5-10.
[10] 王偉, 張晶濤, 柴天佑. PID參數(shù)先進整定方法綜述[J]. 自動化學(xué)報, 2000, 26(3): 347-355.
WANG Wei, ZHANG Jingtao, CHAI Tianyou. A survey of advanced PID parameter tuning methods[J]., 2000, 26(3): 347-355.
[11] 陶永華, 尹怡欣, 葛蘆生. 新型PID控制及其應(yīng)用[M]. 北京: 機械工業(yè)出版社, 1998.
TAO Yonghua, YIN Yixin, GE Lusheng.[M]. Beijing: China Machine Press, 1998.
[12] 李小亭, 韓冰, 李正坤. 關(guān)于精密控溫方法的比較研究[J]. 河北大學(xué)學(xué)報: 自然科學(xué)版, 2004, 24(1): 107-111.
LI Xiaoting, HAN Bing, LI Zhengkun. Comparative study on precision temperature control methods[J].: Natural Science Edition, 2004, 24(1): 107-111.
[13] Leroy C, Maisonneuve M, Piat M, et al. Simulation of the planck-HFI thermal control system[C]//, 2008, 7017: 701713.
[14] Allidina A Y, Hughes F M. Generalized self-tuning controller with pole-assignment[J].,, 1980, 127(1): 13-18(DOI:10.1049/ip-d:19800003).
[15] 劉瑤瑤, 楊寶玉, 吳亦農(nóng). 制冷機PID溫度控制參數(shù)設(shè)計[J]. 紅外, 2016, 37(4): 7-10.
LIU Yaoyao, YANG Baoyu, WU Yinong. Design of PID temperature control parameters for cryocooler [J]., 2016, 37(4): 7-10.
[16] 劉伯春. 離散時間自適應(yīng)控制算法的一種新統(tǒng)一格式[J]. 控制理論與應(yīng)用, 1989(6): 59-64.
LIU Bochun. A new unification of discrete-time adaptive control algorithms[J]., 1989(6): 59-64.
Adaptive PID Control Method Based on Space Optical Mechanical Thermal Model
LI Shuai1,2,YANG Baoyu1,LU Yan1
(1.,,200083,;2.,100049,)
To improve the temperature control stability of space cameras to ensure imaging quality, an adaptive proportional-integral-derivative(PID) control method based on a space optical mechanical thermal model is proposed. The design of the controller starts from the thermal balance equation of the space optical machinery and can correct the thermal model of the optical machinery in real time according to the temperature of the optical machinery and its radiating heat exchange object. Then, the parameters of the PID controller are corrected in real time using the pole assignment method, and the heating duty cycle of the temperature control period is finally determined. In this study, by establishing an abstract thermal model of space optical machinery and applying the above self-adaptive PID control method and PID control method with fixed parameters, the effect of temperature control is compared by simulation and experiment. The experimental results show that the adaptive PID controller always maintains the best dynamic response to the temperature fluctuation caused by environmental disturbance, and the temperature control stability is better than ±0.1K; thus, the controller has better temperature control stability and environmental adaptability.
space optical machinery, thermal model, pole assignment, self-adaption, PID control
TK323
A
1001-8891(2021)10-0934-06
2020-01-16;
2020-03-16.
李帥(1992-),男,河南省南樂縣人,碩士研究生,主要從事熱控技術(shù)研究。
楊寶玉(1979-),副研究員,研究生導(dǎo)師,主要從事航天有效載荷項目中的控制電路及相關(guān)信號處理研究。E-mail:byyang@mail.sitp.ac.cn。