無人機結冰探測器仿真研究

熊海霞，袁冬莉，張正明

（西北工業大學 自動化學院，陜西 西安 710129）

結冰探測器為無人機環境控制系統的重要組成部分，其主要是用來為無人機提供結冰告警信號。特定氣象條件下，在飛機的升力表面（如機翼、螺旋槳和旋翼、發動機進氣道、擋風玻璃、外露傳感器等）迎風表面會產生水分凝結成冰的現象。當發生輕度結冰時，會使飛機升阻比大幅下降，進而造成飛行姿態控制困難，嚴重結冰時可能造成飛機操縱翼面發生失效[1]。所以在飛機上安裝結冰探測器為飛機結冰情況提供預警信號，具有重要意義。

在無人機的開發過程中，往往是硬件和軟件設計并向進行，這樣可以大大縮減研制周期。而在軟件的開發過程中又需要硬件的接口數據，這就要求在硬件成形之前就要有數據輸出，建模仿真為此提供了很好的解決方法。文中建立了無人機結冰探測器的數學模型，為無人機控制軟件的設計與調試提供了重要的數據，縮短了飛機的開發周期。

1 結冰探測器介紹

結冰探測器探測飛機的結冰情況，對結冰速率進行實時解算，輸出結冰告警信號、結冰速率信號和結冰探測器故障信號。

結冰探測器由敏感部件、解算部件組成。解算部件給敏感部件提供激勵信號，敏感部件的探頭在激勵信號和位于敏感部件上的電磁系統作用下，產生穩定的振動即工作頻率。當有冰沉積到敏感部件的探頭上時，其質量變化，從而使工作頻率降低。解算部件根據頻率的變化，判定探頭上已有冰層的厚度，解算結冰速率。

結冰探測器具有自身加熱除冰的功能，輸出結冰速率信號后，自身電加熱除冰。圖1給出了結冰探測器的結構框圖。

圖1 結冰探測器結構框圖Fig.1 Structure diagram of the icing probe

為了建模方便需要對結冰探測器的物理框圖進行簡化，如圖2所示，將結冰探測器的數學模型簡化為4部分：結冰模擬、探頭頻率模擬、加熱模擬、結冰厚度解算。模型的輸出為結冰速率信號、結冰告警信號、結冰探測器故障信號。

圖2 結冰探測器簡化框圖Fig.2 Simplified structure diagram of the icing prob

2 結冰探測器模型

2.1 結冰模擬模塊模型

該模塊模擬結冰過程，其中rate模擬當前環境中探頭的結冰速率，heating模擬加熱除冰過程，thick模擬結冰厚度。當結冰厚度達到1 mm或結冰厚度達到0.5 mm超過60 s時，heating為1，模擬給探頭加熱除冰。探頭結冰厚度初始值為0，結冰厚度后面的飽和環節設置了結冰厚度的限定值，當前設置的上限值為2.5 mm （由于在探頭厚度達到1 mm后，探頭即加熱除冰，且探頭加熱功率達到300 W，因此探頭處的結冰厚度并不會達到2.5 mm），下限值為0 mm。如圖3所示為結冰模擬模塊結構。

圖3 結冰模擬模塊Fig.3 Structure diagram of icing

2.2 加熱模塊模型

該模塊在符合要求的條件下模擬給探頭加熱除冰，其中iceThick為輸入模擬結冰厚度，heat為輸出模擬加熱功能。加熱10 s后，若結冰厚度仍超過0.5 mm且持續60 s或結冰厚度達到1 mm，則對探頭繼續加熱，否則，結束加熱。如圖4所示。

圖4 加熱模塊Fig.4 Structure diagram of heating

2.3 探頭頻率模擬模塊模型

該模塊模擬結冰探測器探頭對應不同結冰厚度下的輸出頻率[2]，其中frequency為輸出模擬探頭頻率。如圖5所示。

圖5 探頭頻率模擬模塊Fig.5 Structure diagram of the probe

2.4 解算結冰速率模塊模型

該模塊模擬根據不同輸出頻率解算結冰厚度，同時該模塊實現以下功能：

1）探頭處結冰厚度達到0.5 mm時，輸出結冰告警信號；

2）達到1 mm后輸出結冰速率信號；

3）加熱結束后60 s內結冰厚度未達到0.5 mm，告警信號結束，告警信號結束的同時不再輸出結冰速率信號。

圖6為解算結冰速率模塊結構，其中frequency和heat為輸入，thickness、icewarn、iceRate為輸出分別代表解算出的結冰厚度、結冰告警和解算出的結冰速率。iceRateTest模塊模擬解算結冰速率后的輸出裝置。

3 模型仿真實現

當初始結冰速率為0.8 mm/s并在400 s時結冰速率轉換為2 mm/s，仿真曲線如圖7所示。

圖7中，當結冰速率為0.8 mm/s時，結冰63 s后，結冰厚度達到0.5 mm，此時結冰探測器告警信號變為1，發出結冰告警信號。但由于此時結冰厚度未達到1 mm，所以未輸出結冰速率信號。在123 s時，結冰厚度達到0.5 mm以上已經持續60 s，此時結冰探測器開始加熱除冰。在126 s時，結冰厚度下降到0.5 mm以下，結冰告警信號并未立即消失，而是在186 s即60 s后消失。在400 s時，結冰速率由0.8 mm/s變為2 mm/s。由于結冰厚度達到1 mm，所以當結冰速率為2 mm/s時，輸出結冰速率信號為2 mm/s。在421.5 s時，結冰厚度在加熱后降為0.5 mm以下，但當478 s時，結冰厚度又達到了0.5 mm以上，由于時間間隔未到60 s，所以結冰告警信號一直未消失。同時由于結冰速率增大，結冰加快，所以結冰探測器加熱頻率增快。

4 實時仿真

4.1 實時仿真平臺組成

結冰探測器仿真系統包括：兩臺PC機、兩臺工控機[3-4]。一臺PC機運行LabVIEW[5]控制界面，另一臺PC機運行LabV IEW顯示界面。一臺工控機運行結冰探測器模型，另一臺工控機運行實現控制功能的軟件，這兩臺工控機的操作系統均為VxWorks。實時仿真系統的具體構成如圖8所示。

圖6 解算結冰速率模塊結構Fig.6 Structure diagram of icing rate

圖7 結冰探測器仿真曲線Fig.7 Simulation results of the icing probe

圖8 實時仿真系統Fig.8 Real-time simulation system of the icing probe

兩臺PC機在XP系統下分別運行LabVIEW編寫的控制界面和顯示界面，控制界面用來輸入控制指令，這些控制指令通過以太網傳輸到兩臺工控機中；顯示界面用來顯示模型的運算結果，顯示界面可以查看模型結果是否符合預期要求。兩臺工控機均搭載VxWorks操作系統，分別為模型機和控制軟件機，模型機上運行的是在Simulink環境下建立的并經過RTW下載技術生成的結冰探測器數學模型可執行代碼；控制軟件機運行的是對模型進行控制和故障處理的軟件。

4.2 文件配置

由于軟件需要與多個外部系統進行收發數據的交互，故該仿真平臺采用UDP通信實現兩臺PC機與兩臺工控機間的數據通信。

為了能夠將建立好的仿真模型順利轉化為可加載的目標文件，需要對RTW以及Tornado做相應的設置，不同MATLAB版本的配置過程可能會略有區別，本文使用的是MATLAB R2012a版本，配置過程見參考文獻[6]。

5 結 論

本文通過對無人機結冰探測器功能和部分實驗數據的分析，建立了結冰探測器的Simulink模型能夠較好的模擬出無人機結冰探測器的工作特性。以搭載由VxWorks嵌入式實時操作系統的ARK-5260工控機和XP操作系統的PC機為仿真平臺，并利用RTW下載技術將建立的Simulink模型下載到工控機中實現了對無人機結冰探測器的實時仿真，其仿真結果在LabVIEW顯示界面中進行實時顯示，此仿真平臺很好的模擬了結冰探測器的功能，為類似仿真提供了一定的參考。

[1]蔣天俊.結冰對飛機飛行性能影響的研究[M].南京:南京航空航天大學，2008.

[2]王穎.壓電諧振式結冰傳感器數學模型研究 [M].武漢:華中科技大學，2006.

[3]鮑泳林，袁冬莉，張宇坤.某型飛機燃油系統交輸供油功能實時仿真[J].科學技術與工程，2012，20（25）:6390－6394.BAO Yong-lin，YUAN Dong-li，ZHANG Yu-kun.Real-time simulation of the cross-feed function of one aircraft fuel system[J].Science Technology and Engineering，2012，20（25）:6390－6394.

[4]張正明，袁冬莉，呂鵬.某型無人機液冷裝置的實時仿真[J].電子設計工程，2013，21（15）:1－3.ZHANG Zheng-ming，YUAN Dong-li，LV Peng.Real-time simulation of the liquid cooling system of one UAV[J].Electronic Design Engineering，2013，21（15）:1－3.

[5]楊樂平.LabVIEW高級程序設計[M].北京:清華大學出版社，2003.

[6]張正明.無人機燃油、環控系統仿真與機電軟件設計[M].西安:西北工業大學，2014.