核磁共振陀螺中原子氣室溫度場的研究

易 鑫，汪之國，夏 濤，徐 迪，楊開勇

(國防科技大學 光電科學與工程學院，湖南 長沙 410073)

?

核磁共振陀螺中原子氣室溫度場的研究

易 鑫，汪之國*，夏 濤，徐 迪，楊開勇

(國防科技大學 光電科學與工程學院，湖南 長沙 410073)

為研究核磁共振陀螺中加熱機構(gòu)對原子氣室性能的影響，設計了5種典型加熱方式。利用有限元分析軟件ANSYS建立了原子氣室的溫度場模型，給出了原子氣室表面的穩(wěn)態(tài)溫度場分布情況。同時設計了探測精度為0.01 ℃的測溫電路，對原子氣室表面不同位置的溫度進行監(jiān)控，獲得了不同加熱方式下原子氣室表面的溫度變化情況。將仿真和實驗結(jié)果進行比較，發(fā)現(xiàn)誤差在5%之內(nèi)，驗證了仿真模型的正確性。綜合仿真和實驗結(jié)果比較了不同加熱方式下原子氣室表面溫度分布情況，獲得了能夠使原子氣室表面溫度分布最均勻的加熱方式。

核磁共振陀螺；原子氣室；加熱機構(gòu)；溫度場

1 引 言

核磁共振陀螺是一種體積小、功耗低、精度高的新型陀螺儀，是下一代高精度微型陀螺儀的首選方案，目前已成為慣性導航領域中的研究熱點[1-3]。核磁共振陀螺中，原子氣室是核心部件。原子氣室內(nèi)充有一定量的銣蒸汽、氙氣等氣體，內(nèi)壁上鍍有抗弛豫膜層。為了使原子氣室具備優(yōu)異的性能，通常需要對其加熱到一定溫度并保持溫度均勻。不同加熱方式下，原子氣室的溫度場將不相同。同時，在實際情況中，原子氣室溫度場分布無法達到絕對均勻，溫度梯度會始終存在，并在一定程度上影響原子氣室內(nèi)物質(zhì)的分布均勻性，這是研究核磁共振陀螺所必須考慮的問題。因此，如何通過有效的方式獲得原子氣室溫度場分布信息并使得其均勻性良好，是核磁共振陀螺研究中的重要問題。

為了測量物體的溫度，通常采用熱敏電阻或者熱成像儀等。但是在核磁共振陀螺儀中，單純使用熱敏電阻只能探測到原子氣室部分位置的溫度值，無法獲得整個原子氣室的溫度信息；熱成像儀對于尺寸小的原子氣室測量精度不高，且在核磁共振陀螺中使用不便。這兩種方法都無法有效地獲得原子氣室溫度場分布信息。文獻[4-5]利用有限元軟件模擬并基于熱敏電阻檢測的手段，獲得了光纖陀螺的溫度場信息，文獻[6-7]利用同樣的方法獲得了激光陀螺的溫度場信息。本文將這種測溫方法應用到核磁共振陀螺中，獲得了核磁共振陀螺中原子氣室的溫度場信息。

在核磁共振陀螺的加熱方式上，文獻[8]采用加熱片加熱銅塊，通過銅塊傳遞熱量給中間的原子氣室；文獻[9]在小型原子氣室的上側(cè)貼有加熱片，利用加熱片直接傳遞熱量給原子氣室。但上述文獻中卻并未涉及原子氣室溫度場的分布信息，未考慮溫度不均勻的影響。本文根據(jù)實驗室現(xiàn)有的加熱機構(gòu)，設置了5種典型的加熱方式對原子氣室進行加熱并對溫度場進行分析。利用有限元分析軟件ANSYS建立了核磁共振陀螺加熱機構(gòu)的溫度場模型，介紹了熱分析的基本原理，給出了核磁共振陀螺中加熱機構(gòu)的結(jié)構(gòu)和有限元模型，同時給出了熱導率、熱生成率、熱輻射率等物理參數(shù)，得到了穩(wěn)態(tài)情況下原子氣室表面的溫度場分布信息。同時利用高精度的鉑電阻和測溫電路，測量了各加熱方式下原子氣室表面3個典型測溫點的溫度變化情況。綜合仿真和實驗結(jié)果，給出了不同加熱方式下原子氣室表面溫度場分布信息，并得到了現(xiàn)有加熱機構(gòu)中使原子氣室溫度最均勻的加熱方式，為核磁共振陀螺中加熱機構(gòu)的設計以及進一步優(yōu)化提供了良好的參考依據(jù)。

2 核磁共振陀螺加熱機構(gòu)原理與有限元模型

2.1 熱分析的基本原理

在熱量傳遞過程中，有3種基本傳遞方式：熱傳導、熱對流和熱輻射。

式中，q″是熱流密度，λ是導熱系數(shù)，gradt是空間某點的溫度梯度。

熱對流[11]是指溫度不同的流體之間發(fā)生相對運動所導致的熱量傳遞方式，包括自然熱對流和強迫熱對流兩種方式。熱對流滿足牛頓冷卻方程[12]：

式中，h是對流換熱系數(shù)，Ts是固體表面溫度，Tb是周圍流體溫度。

熱輻射[12]是指物體發(fā)射電磁能并被其它物體吸收的過程。與熱傳導和熱對流不同，熱輻射不需要傳熱介質(zhì)。熱輻射滿足斯蒂芬玻爾茲曼方程[13]：

式中,q為熱流率，ε為輻射率，σ為黑體輻射常數(shù)，A1為輻射面1的面積，F(xiàn)12為輻射面1到輻射面2的形狀系數(shù)，T1和T2分別為輻射面1和2的絕對溫度。

上述3式是進行熱分析的數(shù)學基礎。

ANSYS軟件進行穩(wěn)態(tài)熱分析的一般方程[14]為：

在磁浮列車的研究過程中，車輛和軌道之間的共振問題一直是人們研究的重點問題之一。在不同軌道條件、不同速度和不同車輛工況下，如何解決車軌耦合振動問題，一直是磁浮研究的難題。在長沙磁浮快線懸浮控制系統(tǒng)調(diào)試中，同樣面臨車軌共振問題。由于長沙磁浮快線 18.55 km的線路條件復雜，車軌共振問題更加復雜。在調(diào)試的初始階段，車輛在維修軌道、車站軌道以及部分線路上都存在車軌共振現(xiàn)象。為了對這個問題進行系統(tǒng)研究，建立了中低速磁浮列車車軌共振試驗臺(見圖2)。該試驗臺包括1個單懸浮架結(jié)構(gòu)、1套軌道及振動激勵系統(tǒng)和1套懸浮控制系統(tǒng)。

式中，[K]是傳導矩陣，{T}是節(jié)點溫度向量，{Q}是節(jié)點熱流向量。

2.2 核磁共振陀螺加熱機構(gòu)結(jié)構(gòu)及有限元模型

本文使用的核磁共振陀螺加熱機構(gòu)如圖1所示，包括外部的保溫框和內(nèi)部的底座、加熱片、原子氣室4個主要部件。保溫框和底座是由聚四氟乙烯材料制成，保溫框的邊長是50 mm，中間設置有通光孔；加熱片包括上、下、前、后、左、右加熱片，分別對稱地放置在原子氣室周圍，由細銅絲纏繞在云母片上制成，貼在保溫框內(nèi)壁上，加熱片尺寸相同，為22 mm×22 mm×2 mm，中間打有12 mm直徑的圓孔，且電阻值均為23 Ω；原子氣室位于保溫框內(nèi)部，由pyrex材料制成，其內(nèi)部充了一定量的銣蒸汽、氙氣、氮氣和氫氣，原子氣室尺寸為20 mm×20 mm×20 mm，內(nèi)壁厚度為1 mm。

圖1 核磁共振陀螺中加熱機構(gòu)示意圖 Fig.1 Schematic of the heating mechanism in NMRG

實驗中分別使用5種加熱方式對原子氣室進行加熱。其中加熱方式1僅上加熱片加熱；加熱方式2前、后加熱片同時加熱；加熱方式3前、后、左、右加熱片同時加熱；加熱方式4上、下、前、后加熱片同時加熱；加熱方式5使用全部加熱片進行加熱。并且各加熱方式總加熱功率相同。

在實際加熱機構(gòu)中，為了探測溫度，在加熱機構(gòu)內(nèi)放置了鉑電阻測溫電路。由于鉑電阻和導線對溫度的影響比較小，且不便于建模，故對核磁共振陀螺加熱機構(gòu)溫度模型進行簡化，得到了如圖2所示的有限元模型。有限元模型共劃分節(jié)點115 500個、單元41 434個。

圖2 核磁共振陀螺加熱機構(gòu)有限元模型 Fig.2 Finite element model of the heating mechanism in NMRG

2.3 參數(shù)計算

為了對有限元模型進行求解，需要獲得模型中各材料的熱導率、熱源的熱生成率、各輻射面的輻射率以及各接觸面的對流換熱系數(shù)。

在加熱機構(gòu)中，保溫框和底座是由聚四氟乙烯制成；加熱片是由導線繞制在云母片上，主要考慮云母片的熱導率；原子氣室是由pyrex材料制成，氣室內(nèi)有多種氣體存在，由于氮氣的量占多，因此只考慮氮氣的熱導率。材料熱參數(shù)如表1所示。

表1 材料熱參數(shù)

加熱機構(gòu)中的各加熱片，體積為7.4×10-7m3，電阻為23 Ω。在5種不同的加熱方式下，由于電流的不同，單個加熱片的熱生成率也不同。表2是不同加熱方式下單個加熱片的熱生成率。

表2 不同加熱方式下單個加熱片的熱生成率

再考慮熱對流對溫度場的影響。在本加熱機構(gòu)中，換熱面較多，內(nèi)部換熱比較復雜，主要考慮以下的換熱情況：(1)加熱框外表面與外界環(huán)境的無限空間對流換熱；(2)原子氣室前、后側(cè)外表面和前后加熱片外表面的有限空間對流換熱；(3)原子氣室左、右側(cè)外表面和左、右加熱片外表面的有限空間對流換熱；(4)原子氣室上側(cè)外表面和上加熱片的有限空間對流換熱；(5)原子氣室下側(cè)外表面和下加熱片的有限空間對流換熱。由于不同加熱方式下溫度的分布不同，因此原子氣室的上、下、前、后、左、右側(cè)各外表面的對流換熱不能夠直接等效，而應該分別進行考慮，并計算各對流換熱系數(shù)[15]。

熱輻射也是影響溫度場分布的重要因素。在本加熱機構(gòu)中，主要考慮兩種熱輻射情況：(1)加熱框外表面和外界環(huán)境的熱輻射；(2)原子氣室外表面和各加熱片外表面之間的熱輻射。其中，加熱框的熱輻射率為0.91，加熱片的熱輻射率為0.91，原子氣室的熱輻射率為0.9。

3 穩(wěn)態(tài)溫度場分析

圖3 加熱方式1下原子氣室表面穩(wěn)態(tài)溫度場 Fig.3 Temperature field of the surface of atom vapor cell under heating method 1

圖4 加熱方式2下原子氣室表面穩(wěn)態(tài)溫度場 Fig.4 Temperature field of the surface of atom vapor cell under heating method 2

圖5 加熱方式3下原子氣室表面穩(wěn)態(tài)溫度場 Fig.5 Temperature field of the surface of atom vapor cell under heating method 3

圖6 加熱方式4下原子氣室表面穩(wěn)態(tài)溫度場 Fig.6 Temperature field of the surface of atom vapor cell under heating method 4

圖7 加熱方式5下原子氣室表面穩(wěn)態(tài)溫度場 Fig.7 Temperature field of the surface of atom vapor cell under heating method 5

在ANSYS中設置完各參數(shù)并對模型加載之后，設置初始溫度為25 ℃，并對穩(wěn)態(tài)溫度場進行求解，如圖3～圖7所示，得到各加熱方式下原子氣室表面的穩(wěn)態(tài)溫度場分布。

在圖3～圖7中，可以很明顯地看到原子氣室表面不同的顏色，通過比照左邊的溫度標尺可以得到相應位置的大致溫度，再通過ANSYS軟件內(nèi)置的溫度探針，得到精確的溫度值。同時，通過各圖也可以得到原子氣室溫度的最大值和最小值。

4 溫度測量與實驗結(jié)果

為了獲得實際加熱機構(gòu)的溫度數(shù)據(jù)，并驗證仿真結(jié)果的正確性，在原子氣室左側(cè)面的3個測溫點放置了高精度的PT1000鉑電阻，利用測溫電路來探測這3個位置的溫度，測溫點位置如圖8所示。

圖8 測溫點位置示意圖 Fig.8 Position of temperature detecting point

當溫度變化時，鉑電阻的阻值發(fā)生改變，通過測溫電路的放大作用，最后輸出電壓信號到采集卡上，利用LABVIEW程序采集電壓信息并得到探測位置的溫度值，并記錄下來。程序的采樣速率為10 kHz，溫度測量精度約為0.01 ℃。不同加熱方式下各測溫點仿真和實驗數(shù)據(jù)如表3所示。

表3 各加熱方式下各測溫點仿真和實驗結(jié)果

通過實驗結(jié)果和仿真結(jié)果相比較可以看出，最大相對誤差在5%以內(nèi)，驗證了核磁共振陀螺加熱機構(gòu)仿真模型的正確性。

5 分析與討論

在核磁共振陀螺中，原子氣室的溫度場分布受多種因素的共同影響。不同加熱方式下，熱源分布會不一樣，加熱機構(gòu)中的熱傳導、熱對流和熱輻射效應也不相同，進而使得原子氣室的溫度場分布產(chǎn)生差異。實驗結(jié)果驗證了仿真模型的正確性，因此可以根據(jù)仿真模型中原子氣室表面溫度場分布情況來比較各加熱方式的優(yōu)劣。不同加熱方式下原子氣室的最大溫差仿真值如表4所示。

表4 各加熱方式下原子氣室表面溫度(單位：/℃)

通過表4可以發(fā)現(xiàn)，加熱方式3對應的原子氣室表面溫度均勻性最佳，其次是加熱方式5，而加熱方式2和4對應的溫度均勻性較差，加熱方式1對應的均勻性最差。

由于加熱方式1中的加熱片不具有對稱性，而另外4種方式中加熱片都是對稱放置，因此方式1中的溫度場分布均勻性最差，并在原子氣室中產(chǎn)生了一個從上到下逐漸降低的溫度梯度。加熱方式2加熱片雖然是對稱放置，但是由于只有兩個加熱片，對于立方體結(jié)構(gòu)的原子氣室來說，溫度均勻性會較差。

再比較加熱方式3、4、5。原子氣室的溫度會受到熱輻射、熱傳導和熱對流的共同影響，由于各加熱片的位置對稱，因此這3種方式下的熱輻射效應在各個方向上都是對稱的；再來看熱對流和熱傳導效應，根據(jù)加熱機構(gòu)的結(jié)構(gòu)可以發(fā)現(xiàn)，在前、后、左、右4個方向上這兩種效應是對稱的，但是在上、下這兩個方向上，由于原子氣室和下側(cè)底座直接接觸，同時原子氣室上表面和下表面接觸的氣體空間和另外4個表面接觸的氣體空間不一樣，導致熱對流和熱傳導效應在這兩個方向上和其它方向并不相同。因此，使用前、后、左、右加熱片的方式4相比于僅使用前、后加熱片的方式2并沒有改善溫度均勻性，同樣的，使用全部加熱片的方式5相比于使用前、后、左、右加熱片的方式3也沒有改善溫度均勻性。方式3由于在熱輻射、熱傳導和熱對流3個效應上具有更加良好的各向同性，因此該方式下的原子氣室表面溫度均勻性達到最佳。

6 結(jié) 論

本文在加熱機構(gòu)中設計了5種不同的加熱方式對原子氣室進行加熱，并建立了加熱機構(gòu)的有限元模型，對穩(wěn)態(tài)溫度場進行了仿真，得到了不同加熱方式下的原子氣室表面溫度場分布。同時利用鉑電阻和測溫電路，對測溫點溫度進行測量，得到的溫度實驗值和仿真值吻合度比較高，相對誤差在5%之內(nèi)，驗證了仿真模型的正確性。綜合仿真和實驗結(jié)果比較了這5種加熱方式下原子氣室表面溫度場均勻性的差異，得到了現(xiàn)有加熱機構(gòu)中溫度均勻性最佳的加熱方式。本文的測溫方法和結(jié)果為核磁共振陀螺中加熱機構(gòu)的設計與優(yōu)化提供了良好的參考。

[1] 劉院省,王巍,王學鋒.微型核磁共振陀螺儀的關鍵技術及發(fā)展趨勢[J].導航與控制,2014,13(4):1-6. LIU Y SH,WANG W,WANG X F. Key technology and development tendency of micro nuclear magnetic resonance gyroscope[J].NavigationandControl,2014,13(4):1-6.(in Chinese)

[2] 秦杰,汪世林,高溥澤,等.核磁共振陀螺技術研究進展[J].導航定位與授時,2014,1(2):64-69. QIN J,WANG SH L,GAO F Z,etal.. Advances in nuclear magnetic resonance gyroscope[J].NavigationPositioning&Timing,2014,1(2):64-69.(in Chinese)

[3] 周鵬,李夜蘭.國外核磁共振陀螺技術專利研究[J].戰(zhàn)術導彈技術,2015(4):5-9. ZHOU P,LI Y L. Patent information research on the nuclear magnetic resonance gyroscope technology[J].TacticalMissileTechnology,2015(4):5-9.(in Chinese)

[4] 張耀麒,曾輝,牟炳富,等.光纖陀螺光纖環(huán)溫度瞬態(tài)效應的模擬仿真[J].半導體光電,2014,35(2):197-200. ZHANG Y Q,ZENG H,MOU B F,etal.. Simulation of temperature transient effects in FOG fiber coil[J].SemiconductorOptoelectronics,2014,35(2):197-200.(in Chinese)

[5] 陳靜,董鈞.光纖陀螺非線性誤差的溫度特性研究[J].導航與控制,2014,13(1):33-38. CHEN J,DONG J. Research about nonlinear error of fiber optic gyroscope under different temperature[J].NavigationandControl,2014,13(1):33-38.(in Chinese)

[6] 于旭東,張鵬飛,湯建勛,等.機抖激光陀螺溫度場的有限元模擬與實驗[J].光學 精密工程,2010,18(4):913-920. YU X D,ZHANG P F,TANG J X,etal.. Finite element analysis and experiments of temperature fields of mechanically dithered ring laser gyroscopes[J].Opt.PrecisionEng.,2010,18(4):913-920.(in Chinese)

[7] 楊慧哲.溫度誤差對激光陀螺零偏溫度補償精度的影響分析[J].光學與光電技術,2014,12(2):98-100. YANG H Z. Effect of temperature error on the temperature compensation accuracy of RLG′s bias[J].Optics&OptoelectronicTechnology,2014,12(2):98-100.(in Chinese)

[8] WYLLIE IV R. The development of a multichannel atomic magnetometer array for fetal magnetocardiography[D]. Wisconsin:University of Wisconsin-Madison,2012.

[9] BULATOWICZ M,LARSEN M. Compact atomic magnetometer for global navigation[C]. 2012 IEEE on Position Location and Nabigation Symposium,IEEE,2012:1088-1093.

[10] 許光映,王晉寶,韓志.基于分數(shù)階熱傳導方程激光加熱瞬態(tài)溫度場研究[J].應用數(shù)學和力學,2015,36(8):844-854. XU G Y,WANG J B,HAN ZH. Study on the transient temperature field based on the fractional heat conduction equation for laser heating[J].AppliedMathematicsandMechanics,2015,36(8):844-854.(in Chinese)

[11] 彭嵐,毛娜,王飛,等.雙向溫度梯度作用下Cz池內(nèi)熱對流的數(shù)值模擬[J].工程熱物理學報,2015(6):1325-1328. PENG L,MAO N,WANG F,etal.. Numerical simulation of thermocapillary convection in czochralski configuration with bidirectional temperature gradients[J].J.EngineeringThermophysics,2015(6):1325-1328.(in Chinese)

[12] 李偉斌,熊永紅.熱輻射基礎實驗[J].物理實驗,2003,23(1):9-12. LI W B,XIONG Y H. Foundational experiments of thermal radiation[J].PhysicalExperiment，2003,23(1)：9-12.(in Chinese)

[13] 劉穎,范永青,徐金濤.光纖陀螺組合的熱設計及分析[J].紅外與激光工程,2015,44(12):3750-3754. LIU Y,FAN Y Q,XU J T. Thermal design and analysis for fiber optic gyroscope combination[J].InfraredandLaserEngineering,2015,44(12):3750-3754.(in Chinese)

[14] 吳國勇,顧啟泰,鄭辛,等.環(huán)形激光陀螺溫度模型[J].清華大學學報:自然科學版,2003,43(2):180-183. WU G Y,GU Q T,ZHENG X,etal.. Thermal model of ring laser gyro[J].JTsinghuaUniv(Sci & Tech),2003,43(2):180-183.(in Chinese)

[15] 楊世銘，陶文銓.傳熱學[M].北京：高等教育出版社，2006. YANG SH M,TAO W Q.HeatTransfer[M]. Beijing:Advanced Educational Press,2006.(in Chinese)

*Correspondingauthor,E-mail:eng.rad.laham@gmail.com

Abstract: High performance optical systems are characterized by high sensitivity to assembly procedure of the system and to any sensible change in environmental conditions. This sensitivity issue is handled in this paper through a methodology allowing a computational prediction of optical performance when the opto-mechanical system is subject to some external factors. This paper explains the methodology through an example of an optical objective of excellent performance undergoing static mechanical stresses which degrade the performance expressed by an MTF diagram. Then the objective is manufactured, assembled, and an optical interferometer is used to test the objective when stress is retained; and the experimental results of degraded MTF are compared to the analytical MTF. The excellent matching between the two sets of results confirms the validity of the proposed methodology.

Key words: predicting optical performance;assembly of optical system;tensile torque;MTF;objective;Ansys;Matlab

Research on temperature field in the vapor cell of nuclear magnetic resonance gyroscope

YI Xin, WANG Zhi-guo*, XIA Tao, XU Di, YANG Kai-yong

(CollegeofOptoelectricScienceandEngineering,NationalUniversityofDefenseTechnology,Changsha410073,China)

*Correspondingauthor,E-mail:maxborn@163.com

In order to study the influences of the heating mechanism on the vapor cell of the nuclear magnetic resonance gyroscope(NMRG), several different heating methods are designed to heat the vapor cell. The finite element analysis software ANSYS is used to establish a temperature field model of the vapor cell, afterwards the model is simplified and finally the steady-state temperature field of the surface is given. Meanwhile, a temperature measurement circuit with the precision of 0.01 ℃ is made to detect the temperature of different places on the surface of the vapor cell, thus getting the changes of the surface temperature of different methods. The experimental results are compared with the simulation results, and the calculation error turns out to be less than 5%, which verifies the validity of the temperature model. Based on both of the experimental and simulation results, and after comparison of temperature field of the vapor cell surface through different heating methods, the best heating method is given, which can makes the most uniform temperature field of the vapor cell surface.

nuclear magnetic resonance gyroscope(NMRG);atom vapor cell;heating mechanism;temperature field

Methodology for predicting optical system performance when subjected to static stresses

AL-LAHAM Radwan*, MOUSSELLY Mhd.Fawaz, NAIM Mamoun

(HigherInstituteforAppliedSciencesandTechnology,Damascus,Syria)

2016-06-13；

2016-07-29

國防科技大學校預研資助項目(No.JC140702) Supported by Advance Research Project of National University of Defense Technology(No.JC140702)

2095-1531(2016)06-0671-07

V241.557; TH703

A

10.3788/CO.20160906.0671

易 鑫(1992—)，男，江西宜春人，碩士研究生，2014年于國防科技大學獲得學士學位，主要從事激光與光電子技術方面的研究。E-mail:1091184065@qq.com

汪之國(1982—)，男，山東臨城人，博士，講師，2003年于吉林大學獲得學士學位，2010年于國防科技大學獲得博士學位，主要從事激光與光電子技術方面的研究。E-mail: maxborn@163.com