張沛勇，孫 麗，武麗花

(1.北京控制工程研究所，北京 100190； 2.天津航海儀器研究所，天津 300131)

?

一種二浮陀螺儀熱仿真分析及試驗(yàn)驗(yàn)證

張沛勇1，孫 麗1，武麗花2

(1.北京控制工程研究所，北京 100190； 2.天津航海儀器研究所，天津 300131)

針對(duì)某衛(wèi)星姿態(tài)控制系統(tǒng)中的重要測(cè)量部件二浮陀螺儀的精度及可靠性不斷提高的要求，為了分析、掌握二浮陀螺儀的熱設(shè)計(jì)及內(nèi)部溫度場(chǎng)分布對(duì)其精度和可靠性的影響，利用Ansys Fluent軟件對(duì)二浮陀螺儀溫度場(chǎng)進(jìn)行仿真分析，得出其內(nèi)部溫度場(chǎng)分布情況，并與熱真空試驗(yàn)數(shù)據(jù)進(jìn)行比對(duì)，驗(yàn)證了二浮陀螺儀有限元模型及邊界條件的準(zhǔn)確性，根據(jù)仿真結(jié)果提出改善二浮陀螺儀熱設(shè)計(jì)的優(yōu)化措施.

二浮陀螺儀；熱仿真.

0 引 言

二浮陀螺儀是衛(wèi)星姿態(tài)控制系統(tǒng)中的主要測(cè)量部件之一，用于測(cè)量衛(wèi)星本體相對(duì)慣性空間的運(yùn)動(dòng)角速度.隨著衛(wèi)星對(duì)二浮陀螺儀精度及可靠性需求的不斷提高，二浮陀螺儀內(nèi)部溫度場(chǎng)對(duì)精度的影響凸顯.之前，在進(jìn)行熱學(xué)特性的研究時(shí)，基本上只通過(guò)試驗(yàn)的方式對(duì)設(shè)計(jì)結(jié)果進(jìn)行驗(yàn)證，對(duì)二浮陀螺儀自身的熱學(xué)特性的分析很少，既不能做到提前識(shí)別存在的問(wèn)題，又不利于產(chǎn)品的設(shè)計(jì)改進(jìn).當(dāng)前，隨著數(shù)值計(jì)算方法的不斷發(fā)展，通過(guò)有限元分析方法可獲得陀螺儀溫度場(chǎng)的分布情況[1-4].本文通過(guò)Ansys Fluent有限元分析軟件[5-6]對(duì)某衛(wèi)星用二浮陀螺儀進(jìn)行熱仿真及分析，得到二浮陀螺儀內(nèi)部的溫度場(chǎng)分布情況，并通過(guò)熱真空試驗(yàn)驗(yàn)證二浮陀螺儀熱仿真模型的正確性，為后續(xù)二浮陀螺儀熱設(shè)計(jì)提供參考數(shù)據(jù).

1 二浮陀螺儀熱模型的建立及仿真

1.1 二浮陀螺儀的熱設(shè)計(jì)

二浮陀螺儀三大組件：浮子組件，傳感器組件和殼體組件.浮子組件主要由陀螺儀馬達(dá)、浮子框架和浮筒組成，其內(nèi)部間隙充滿氦氣.二浮陀螺儀熱環(huán)境復(fù)雜，陀螺儀既要加熱保證浮液對(duì)浮子組件的浮力和重力完全相等，又要考慮陀螺儀內(nèi)部馬達(dá)的散熱，因此二浮陀螺儀的熱設(shè)計(jì)至關(guān)重要又有難度.

圖1 二浮陀螺儀溫控系統(tǒng)示意圖Fig.1 Diagram of temperature control system of floated gyro

圖1為二浮陀螺儀溫控系統(tǒng)示意圖.系統(tǒng)中加熱元件采用帶狀薄膜加熱片，纏繞在陀螺儀殼體上，通過(guò)傳導(dǎo)對(duì)浮液加熱.測(cè)溫元件采用熱敏電阻，布置在傳感器組件上.感溫元件采用兩個(gè)熱敏電阻串聯(lián)，一個(gè)內(nèi)熱敏電阻布置在傳感器組件上，另一個(gè)外熱敏電阻布置在殼體外，用硅橡膠固定.上述元件和外部溫控電路共同組成陀螺儀的溫控系統(tǒng)，對(duì)陀螺儀的內(nèi)部溫度進(jìn)行控制，進(jìn)而達(dá)到浮液對(duì)浮子組件的浮力和重力完全相等，消除摩擦，提高陀螺儀精度的目的.

1.2 有限元模型的建立

二浮陀螺儀的幾何結(jié)構(gòu)比較復(fù)雜，幾何模型的建立是前期的一項(xiàng)重要工作.建立的幾何模型，既要保證計(jì)算模型與真實(shí)模型在結(jié)構(gòu)、尺寸上的一致性，又要保證模型的計(jì)算效率.因而在建立模型前對(duì)二浮陀螺儀的結(jié)構(gòu)模型進(jìn)行了適當(dāng)?shù)暮?jiǎn)化.最終得到的有限元幾何模型及網(wǎng)格模型如圖2所示.

圖2 二浮陀螺儀有限元幾何模型及網(wǎng)格模型Fig.2 Geometric model and grid model of floated gyro

1.3 邊界條件的確定

1)內(nèi)熱源的確定

二浮陀螺儀的內(nèi)熱源主要是陀螺儀馬達(dá)功耗3.5 W.

2)主動(dòng)溫控?zé)嵩?/p>

二浮陀螺儀的外熱源為殼體外側(cè)的加熱片功率,根據(jù)測(cè)試數(shù)據(jù)陀螺儀在真空條件下各溫度點(diǎn)加熱片功率如表1所示.

表1 二浮陀螺儀穩(wěn)態(tài)加熱片功率

3)真空輻射邊界條件

二浮陀螺儀外罩在真空下主要是輻射換熱.陀螺儀外罩設(shè)為金屬的半球發(fā)射率，吸收率與發(fā)射率相等，輻射環(huán)境溫度為陀螺儀所處真空環(huán)境溫度.

1.4 仿真結(jié)果

圖3是真空環(huán)境為20 ℃時(shí)，陀螺儀溫度場(chǎng)切片分布云圖，切片位置分別是Y=0，X=0，Z=0.

圖3 陀螺溫度切片分布云圖Fig.3 Slices of temperature field distribution of floated gyro

圖4是真空環(huán)境溫度為20 ℃時(shí)，殼體與浮筒間陀螺儀浮液薄層的溫度場(chǎng)分布云圖.

圖4 陀螺儀浮液薄層溫度分布云圖Fig.4 Slices of temperature field distribution of floated gyro oil lamella

1.5 仿真結(jié)果分析

由圖3知：1)二浮陀螺儀溫度場(chǎng)分布范圍是51.1～57 ℃，溫差為5.9 ℃，最高溫部位是陀螺儀馬達(dá)處，最低溫位于陀螺儀法蘭處.2)整個(gè)二浮陀螺儀的溫度梯度主要集中在浮子內(nèi)氦氣處，溫度范圍53.5～55.5 ℃，溫差為2 ℃，原因是氦氣導(dǎo)熱率小，其向外傳熱主要靠氦氣的流動(dòng).

由圖4知：1)陀螺儀浮液薄層內(nèi)側(cè)溫度場(chǎng)分布范圍是51.6～52.1 ℃，溫度差為：0.5 ℃，浮液層四周是高溫區(qū)域，馬達(dá)軸正對(duì)的中間區(qū)域是低溫區(qū).原因是浮液層四周氦氣流動(dòng)速度快，傳熱快，而中間正對(duì)馬達(dá)軸處氦氣速度幾乎為零.2)陀螺儀浮液薄層外側(cè)溫度場(chǎng)分布范圍是51.5～51.9 ℃，溫度差為：0.4 ℃，溫度分布基本是兩端高，中間溫度最低.原因是中間正對(duì)陀螺儀的法蘭，法蘭通過(guò)陀螺儀支架導(dǎo)熱，熱量傳遞快導(dǎo)致中間溫度低.

2 二浮陀螺儀熱仿真的試驗(yàn)驗(yàn)證

2.1 熱分析試驗(yàn)陀螺儀的設(shè)計(jì)

為了驗(yàn)證仿真模型的正確性，設(shè)計(jì)了熱分析試驗(yàn)陀螺儀.熱試驗(yàn)陀螺儀主要是在目前陀螺儀零件的基礎(chǔ)上，在陀螺內(nèi)部設(shè)置測(cè)溫點(diǎn)[7]，通過(guò)測(cè)量測(cè)溫點(diǎn)的溫度，進(jìn)而了解陀螺儀內(nèi)部溫度場(chǎng)的分布情況.

本設(shè)計(jì)中，共布置5個(gè)測(cè)點(diǎn)(浮子組件內(nèi)部2個(gè)，浮子組件外部3個(gè))，分別為：在浮子組件內(nèi)部，安裝陀螺儀電機(jī)的框架內(nèi)壁(C1)、框架臂壓塊端面(C2)各1個(gè)測(cè)點(diǎn)；在浮子組件外部，浸于浮油內(nèi)的左(C3)、右端部(C4)、殼體安裝法蘭上(C5)各1個(gè)測(cè)點(diǎn).位置如圖5所示.

圖5 二浮陀螺儀內(nèi)部測(cè)溫點(diǎn)的位置圖Fig.5 The position of thermal resistance in floated gyro

2.2 試驗(yàn)內(nèi)容

在真空環(huán)境條件下，選取4個(gè)環(huán)境溫度點(diǎn)，分別為10 ℃、20 ℃、30 ℃、40 ℃.試驗(yàn)時(shí)，將試驗(yàn)陀螺儀置于專用的密封罐內(nèi)，利用真空泵實(shí)現(xiàn)真空狀態(tài)，再將密封罐置于恒溫箱內(nèi).試驗(yàn)陀螺儀的工作溫度設(shè)置在50 ℃，實(shí)時(shí)測(cè)試各測(cè)溫點(diǎn)阻值.

2.3 試驗(yàn)結(jié)果及分析

真空狀態(tài)下，試驗(yàn)陀螺儀的工作溫度設(shè)置為50 ℃ 保持不變，改變恒溫箱的溫度(即環(huán)境溫度)，穩(wěn)態(tài)時(shí)對(duì)應(yīng)的各測(cè)溫點(diǎn)的溫度數(shù)據(jù)如表2所示.

表2 真空狀態(tài)各環(huán)境溫度下對(duì)應(yīng)的各測(cè)點(diǎn)溫度

從表2數(shù)據(jù)可見(jiàn)：1)溫度最高點(diǎn)位于框架內(nèi)壁上，比框架臂壓塊端面溫度高0.1 ℃左右；2)浮子內(nèi)溫度比浮子外高1.5 ℃左右；3)環(huán)境溫度為10、20、30 ℃時(shí)，表體內(nèi)溫度較為接近，相差0.2 ℃左右，當(dāng)環(huán)境溫度升至40 ℃時(shí)，表體內(nèi)溫度比30 ℃高4.5 ℃左右.

2.4 試驗(yàn)結(jié)果與仿真結(jié)果比對(duì)

在第1章仿真模型的基礎(chǔ)上，改變環(huán)境溫度，分別得到10 ℃、30 ℃、40 ℃時(shí)陀螺儀內(nèi)部各測(cè)溫點(diǎn)的溫度值.

圖6是真空環(huán)境溫度為20 ℃時(shí)，陀螺儀內(nèi)部各測(cè)溫點(diǎn)的仿真溫度值.

圖6 真空環(huán)境溫度為20 ℃時(shí)，二浮陀螺儀內(nèi)部各測(cè)溫點(diǎn)的仿真溫度值Fig.6 The temperature values of measurement points inside floated gyro(vacuum ambient temperature is 20 ℃)

表3是真空環(huán)境溫度分別為10 ℃、20 ℃、30 ℃、40 ℃時(shí)，陀螺儀熱試驗(yàn)測(cè)溫點(diǎn)數(shù)據(jù)與有限元仿真的溫度數(shù)據(jù)對(duì)比情況.

表3 二浮陀螺儀熱試驗(yàn)測(cè)點(diǎn)數(shù)據(jù)與仿真數(shù)據(jù)對(duì)比表

由表3知：通過(guò)有限元仿真得到的二浮陀螺儀溫度場(chǎng)各測(cè)點(diǎn)的溫度值與熱試驗(yàn)得到的各測(cè)點(diǎn)的溫度值相差很小，充分證明二浮陀螺儀熱仿真模型準(zhǔn)確性.

3 二浮陀螺儀熱設(shè)計(jì)的優(yōu)化

根據(jù)仿真及試驗(yàn)結(jié)果對(duì)二浮陀螺儀的熱設(shè)計(jì)提出如下優(yōu)化建議：

1)浮液的導(dǎo)熱率低是溫度梯度集中在浮液處的主要原因，因此在不改變浮液其他特性的情況下，提高浮液的導(dǎo)熱率也是一種提高二浮陀螺儀內(nèi)部溫度場(chǎng)均勻性的有效措施;

2)根據(jù)仿真和試驗(yàn)結(jié)果，采取分區(qū)溫控設(shè)計(jì)，調(diào)整加熱功率分布，研制異型云圖式加熱片，也是進(jìn)一步減小陀螺儀內(nèi)部溫度梯度的方法.

4 結(jié) 論

本文通過(guò)Ansys Fluent軟件對(duì)二浮陀螺儀的溫度場(chǎng)進(jìn)行了仿真及分析，并通過(guò)熱真空試驗(yàn)驗(yàn)證了仿真模型的正確性，同時(shí)根據(jù)仿真及試驗(yàn)結(jié)果提出了針對(duì)當(dāng)前二浮陀螺儀熱設(shè)計(jì)的優(yōu)化措施，對(duì)后續(xù)二浮陀螺儀的熱設(shè)計(jì)具有重要的參考意義.

[1] 李德才.基于流固耦合理論的長(zhǎng)壽命陀螺儀浮子內(nèi)部流場(chǎng)及溫度場(chǎng)分析[J].中國(guó)慣性技術(shù)學(xué)報(bào)，2007，15(6)：721-723. LI D C. Analysis on temperature field and inner flow field of long-life gyro's floater with fluid-structure coupling theory[J].Journal of Chinese Inertial Technology，2007，15(6)：721-723.

[2] 卜石，劉娜.某型號(hào)液浮陀螺溫度場(chǎng)分析[J].彈箭與制導(dǎo)學(xué)報(bào)，2006，26(4)：95-98. BU S,LIU N.The analysis on temperature field of a type of floated gyro[J].Journal of Projectiles, Rockets, Missiles and Guidance，2006，26(4)：95-98.

[3] 液浮陀螺溫度場(chǎng)優(yōu)化方法研究及實(shí)現(xiàn)[J].計(jì)算機(jī)測(cè)量與控制.2010,18(9)：2145-2147.

[4] APPEL N,HOLLAND J.Three dimensional finite element analysis of Inertial gyros[C]//AIAA Guidance and Control Conference.Washington D.C.:AIAA,1977：242-246.

[5] 韓占忠.FLUENT—流體工程仿真計(jì)算實(shí)例與應(yīng)用[M].北京理工大學(xué)出版社.2004，185-265.

[6] 田杰，陳培譞，莊園.基于Fluent的ITER杜瓦的輻射換熱分析[J].機(jī)械工程與自動(dòng)化,2012,175(6):4-6 TIAN J, CHEN P X, ZHUANG Y. Thermal radiation analysis of ITER dewar by fluent[J].Mechanical Engineering & Automation,2012,175(6):4-6.

[7] 喻文煥 胡順菊.液浮陀螺儀浮筒表面瞬態(tài)溫度的確定[J].系統(tǒng)科學(xué)與數(shù)學(xué),1983,5(3):216-222 YU W H, HU S J.Identification of transient temperature on the float surface of a floated gyroscope[J].Journal of Systems Science and Mathematical Sciences，1983,5(3)：216-222.

Thermal Simulation Analysis and Experiment Verification for a Type of floated gyro

ZHANG Peiyong1, SUN Li1， WU Lihua2

(1.BeijingInstituteofControlEngineeringBeijing100190,China；2.TianjinNavigationInstrumentResearchInstitute,Tianjin300131,China)

As one of the important measuring instruments in satellite attitude control system, the requirements for accuracy and reliability of floated gyro are increasing. In order to analyze the influences of thermal design and internal temperature field distribution of floated gyro on its accuracy and reliability, the internal temperature field is simulated and analyzed via Ansys Fluent software. By compared with the thermal vacuum test data, the accuracy of the finite element model and boundary conditions of floated gyro are verified. According to the simulation results, some optimization measures are proposed to improve the thermal design of floated gyro.

floated gyro； thermal simulation

2015-03-27

V249

A

1674-1579(2015)05-0039-04

10.3969/j.issn.1674-1579.2015.05.008

張沛勇(1978—)，男，工程師，研究方向?yàn)閼T性姿態(tài)敏感器設(shè)計(jì)；孫 麗(1973—)，女，高級(jí)工程師，研究方向?yàn)閼T性姿態(tài)敏感器設(shè)計(jì).