呂江山，任宗金，張嘯風(fēng)，焦林虎，孟慶增

（大連理工大學(xué)機(jī)械工程學(xué)院，遼寧大連 116024）

0 引言

姿控火箭發(fā)動機(jī)作為一種重要的飛行器軌道控制裝置得到了越來越廣泛的應(yīng)用[1-3]，通過姿控固體火箭發(fā)動機(jī)產(chǎn)生的推力實(shí)現(xiàn)飛行器的姿/軌控制。為實(shí)現(xiàn)對飛行器的精確控制，需要對姿控發(fā)動機(jī)的推力進(jìn)行準(zhǔn)確測量，并且進(jìn)行地面試驗(yàn)[4-6]，從而對飛行器的不同飛行狀態(tài)進(jìn)行研究[7-8]。壓電式傳感器具有動態(tài)特性好等優(yōu)點(diǎn)，適合于進(jìn)行動態(tài)推力測試。但是壓電傳感器的靈敏度會隨著溫度的變化而發(fā)生改變，因此在進(jìn)行測試前需要進(jìn)行溫度加載實(shí)驗(yàn)，探究溫度對測力儀輸出性能的影響，以提高測力儀的測試精度。

針對溫度對測力儀測試精度的影響問題，相關(guān)領(lǐng)域的研究人員開展了大量工作。大連理工大學(xué)的盧江躍等[9-11]研制了一種新型壓電車削測力儀。該測力儀針對溫度對測力儀預(yù)緊力的影響進(jìn)行了熱補(bǔ)償設(shè)計(jì)，使得測力儀可以補(bǔ)償因切削熱引起的預(yù)緊力減小。S Poussiera等[12]采用可編程芯片對溫度應(yīng)變式傳感器進(jìn)行了溫度補(bǔ)償，在-30°~150℃之間能取得較好的補(bǔ)償效果。C Pra?manik等[13]提出了基于人工神經(jīng)網(wǎng)絡(luò)的智能補(bǔ)償方案，通過對網(wǎng)絡(luò)的訓(xùn)練，在25~80℃溫度范圍下漂移誤差與補(bǔ)償前相比減小了80%。綜上，目前少有研究報(bào)道對火箭發(fā)動機(jī)推力測試系統(tǒng)進(jìn)行溫度影響的研究。

針對上述問題，結(jié)合國內(nèi)外研究現(xiàn)狀，本文針對溫度對姿控火箭發(fā)動機(jī)推力測試系統(tǒng)的影響進(jìn)行了探究。首先通過ANSYS仿真軟件建立了姿控火箭推力測試系統(tǒng)的有限元模型，通過溫度仿真得到了測力儀中傳感器的溫度分布；進(jìn)行了測試系統(tǒng)的多物理場耦合仿真，得到測試系統(tǒng)在熱力耦合影響下的最大變形量不超過0.1 mm，符合設(shè)計(jì)要求；探究了測試系統(tǒng)在溫度下影響下的輸出性能，其最大輸出誤差不超過0.5 N，其性能符合設(shè)計(jì)要求。

1 姿控力測試系統(tǒng)結(jié)構(gòu)設(shè)計(jì)與仿真

1.1 壓電測力儀結(jié)構(gòu)

本文采用的壓電測力儀由轉(zhuǎn)接法蘭、上法蘭、三向力壓電傳感器和底座組成，如圖1所示，測力儀中傳感器的布置形式為4點(diǎn)支撐式，具有良好的穩(wěn)定性，可以實(shí)現(xiàn)對三向力/力矩的準(zhǔn)確測量。

圖1 壓電測力儀三維模型

1.2 三向力壓電傳感器設(shè)計(jì)

石英(SiO2)是壓電傳感器中常用的力-電轉(zhuǎn)換材料。石英材料具有晶體結(jié)構(gòu)，無色透明，硬度高，剛度大，莫氏硬度可達(dá)7級。當(dāng)晶體在某個(gè)方向受到力作用時(shí)，電平衡被破壞，在晶體表面積累電荷，從而產(chǎn)生壓電效應(yīng)。通過石英的壓電效應(yīng)可以對石英進(jìn)行不同的切割制造，從而制作成為壓電傳感器的核心力敏元件。

根據(jù)石英晶體在笛卡爾坐標(biāo)系中的不同位置可以將其劃分為X切族或Y切族。X切族的石英晶體其厚度方向平行于X軸，晶片表面垂直于Y軸，且沿坐標(biāo)軸沒有轉(zhuǎn)動，稱之為X切族的原始切割，也稱為x0°切型。同理，Y切族的石英晶體其厚度方向平行于Y軸，晶體表面垂直于X軸，同樣沿坐標(biāo)軸沒有轉(zhuǎn)動，稱之為Y切族的原始切割，也成為y0°切型。

由于本文中使用傳感器均為三向力傳感器，所以晶組中應(yīng)包含一個(gè)x0°單元晶組，用來測量法向力FZ；還應(yīng)包含兩個(gè)y0°單元晶組，用來測量切向力FX和F Y，關(guān)鍵在于兩個(gè)y0°單元晶組的最大靈敏軸互成90°夾角，同時(shí)晶組應(yīng)該具有同樣的輸出極性。

將設(shè)計(jì)好的晶組封裝在傳感器的殼體中并且導(dǎo)入電極用來引出產(chǎn)生的電荷，傳感器的中間設(shè)計(jì)了螺栓孔，通過預(yù)緊螺栓實(shí)現(xiàn)對傳感器上蓋和殼體的連接以及自身預(yù)緊。

1.3 傳感器布局方式分析

三向力傳感器作為測力儀的核心組成部分，其組成個(gè)數(shù)以及不同的布局形式將直接影響測試系統(tǒng)的整體性能?？臻g內(nèi)3個(gè)不共線的店可以確定唯一的一個(gè)平面，所以只要測力儀中的傳感器數(shù)量不小于3個(gè)便可以實(shí)現(xiàn)對力/力矩的測量。同時(shí)，為了降低測力儀的制造和裝配難度，降低測試數(shù)據(jù)解耦的難度，壓電測力儀內(nèi)傳感器的數(shù)量通常為4個(gè)，根據(jù)傳感器個(gè)數(shù)量以及布置方式的不同，測力儀的輸出性能也會存在差異。

4點(diǎn)支撐式測力儀的布局形式如圖2所示，根據(jù)傳感器安裝位置的不同可以將測力儀分為菱形布局形式和矩形布局形式。菱形布局中4個(gè)傳感器分布在坐標(biāo)軸上，每個(gè)傳感器與原點(diǎn)等距離對稱分布。矩形布局中4個(gè)傳感器分布在坐標(biāo)系的4個(gè)象限內(nèi)，與坐標(biāo)軸呈對稱分布。4點(diǎn)支撐式的測力儀整體剛度較高，在受外力作用時(shí)內(nèi)部支反力分布均勻，向間干擾較小，整體輸出性能好，適合于力值較大，力作用形式復(fù)雜的場合。但4支點(diǎn)布局形式屬于超靜定結(jié)構(gòu)，這要求傳感器零件需要有較高的制造和裝配精度。同時(shí)在菱形布置的條件下，當(dāng)作用力沿坐標(biāo)軸時(shí)會出現(xiàn)兩個(gè)傳感器輸出方向的重合，在測試過程中會導(dǎo)致測力儀的測試精度降低，而矩形布局則不會出現(xiàn)此問題，且矩形布局內(nèi)部傳感器應(yīng)力分布更為均勻，更適宜于多維力的測量。綜上所述，本文采用矩形布局形式的測力儀。

圖2 4點(diǎn)支撐式測力儀

1.4 測力儀溫度分布仿真

利用Creo軟件建立測力儀模型，尺寸嚴(yán)格按照1∶1繪制。將Creo中建立的測力儀模型導(dǎo)入ANSYS Work?bench中，其中測力儀轉(zhuǎn)接法蘭、上法蘭、底座和傳感器上、下蓋及殼體材料均為Q235結(jié)構(gòu)鋼，傳感器中力敏元件材料為石英，表1所示為各種材料的性能參數(shù)(包含楊氏模量、泊松比、密度及屈服極限)。

表1 材料參數(shù)屬性

對導(dǎo)入后的模型進(jìn)行參數(shù)設(shè)置，其中測力儀轉(zhuǎn)接法蘭、上法蘭、底座和傳感器上、下蓋及殼體的導(dǎo)熱系數(shù)為60.5 W/(m·K)，傳感器中石英晶片的導(dǎo)熱系數(shù)為1.5 W/(m·K)，根據(jù)現(xiàn)場的實(shí)際環(huán)境，設(shè)定環(huán)境溫度為22℃，考慮到姿控火箭發(fā)動機(jī)噴管的實(shí)際工作特點(diǎn)，熱傳導(dǎo)和輻射換熱對測力儀溫升的影響極小，僅考慮對流換熱的影響，對流換熱系數(shù)為20 W/(m2·℃)。根據(jù)發(fā)動機(jī)設(shè)計(jì)測試時(shí)長，仿真時(shí)間定為300s，作用在測力儀上法蘭的溫度由22℃逐步增加到120℃，仿真結(jié)果如圖3所示。由圖可以看出測力儀的最高溫度出現(xiàn)在與發(fā)動機(jī)進(jìn)行對流換熱的轉(zhuǎn)接法蘭上，最低溫度出現(xiàn)在底座上。

圖3 壓電測力儀溫度分布

由于火箭發(fā)動機(jī)工作產(chǎn)生的溫升對測力儀中傳感器的靈敏度造成影響，因此需要對測力儀中傳感器的溫度進(jìn)行探究，從而得知測力儀溫升對測試精度的影響，在傳感器晶片上施加探針，如圖4所示，得到傳感器的最大溫度為38.07℃，最大溫升不到20℃，符合要求。

圖4 壓電測力儀溫度分布

1.5 測力儀熱變形仿真

由于測力儀不僅需要測量被測對象產(chǎn)生的多維力，還要承擔(dān)整個(gè)發(fā)動機(jī)的重力，因此測力儀需要由較高的剛度。若測力儀的剛度不足會導(dǎo)致局部變形，導(dǎo)致測力儀輸出精度降低，同時(shí)由于在測試過程中發(fā)動機(jī)產(chǎn)生的溫度會對測力儀的力學(xué)性能產(chǎn)生影響，在進(jìn)行分析時(shí)測力儀上板的溫度為由室溫22℃逐步增加到120℃，測力儀的上、下板及連接件材料為結(jié)構(gòu)鋼，導(dǎo)熱系數(shù)為60.5 W/(m·K)，傳感器內(nèi)部為石英晶片，導(dǎo)熱系數(shù)為1.5 W/(m·K)，根據(jù)測試現(xiàn)場的真實(shí)環(huán)境，設(shè)定環(huán)境溫度為22℃，僅考慮對流換熱的影響，對流換熱系數(shù)為20 W/(m2·℃)。測試時(shí)間為300 s，同時(shí)在測力儀的上板中心處X向和Y向施加800 N載荷，其仿真結(jié)果如圖5所示。

圖5 測力儀應(yīng)力分布

在300 s時(shí)測力儀的應(yīng)力達(dá)到最大值，出現(xiàn)在測力儀的上板，此處為測力儀的連接部分，同時(shí)也是主要熱傳導(dǎo)區(qū)域。

在溫度仿真參數(shù)設(shè)置的基礎(chǔ)上約束測力儀的底座，限制其X、Y、Z和MX、MY、MZ六個(gè)自由度。在平行于測力儀轉(zhuǎn)接法蘭的X向和Y向同時(shí)施加800 N的力，以探究測力儀在溫度影響下的受力最大變形量，其仿真結(jié)果如圖6所示。從圖中可以看出測力儀的最大變形量為0.08 mm，小于0.1 mm，符合要求。

圖6 120°C下測力儀總變形量

2 測試系統(tǒng)溫度實(shí)驗(yàn)

測力儀的溫度實(shí)驗(yàn)系統(tǒng)主要由電加熱板、壓電測力單元、電荷放大器、USB數(shù)據(jù)采集卡、計(jì)算機(jī)及軟件等組成，溫度實(shí)驗(yàn)采用電加熱板對測力儀進(jìn)行加熱，采用220 V電源進(jìn)行加熱，如圖7所示，通過K型溫度傳感器實(shí)現(xiàn)溫度的實(shí)時(shí)反饋調(diào)節(jié)，通過給定預(yù)設(shè)溫度，使加熱板保持在特定溫度附近范圍。其中，X向與Y向連接的電荷放大器靈敏度均設(shè)置為8.00 pC/N，Z向連接的電荷放大器靈敏度設(shè)置均為4.00 pC/N。

圖7 測力儀溫度實(shí)驗(yàn)

測力儀在自身溫升時(shí)會產(chǎn)生熱變形，這種熱變形會導(dǎo)致傳感器產(chǎn)生無加載下的輸出，測力儀輸出的數(shù)據(jù)通過Dewesoft軟件進(jìn)行實(shí)時(shí)數(shù)據(jù)記錄，采樣頻率設(shè)置500 Hz，采樣時(shí)間為180 s。加熱板預(yù)設(shè)溫度為120℃，觀察熱傳遞過程中測力儀輸出情況。實(shí)驗(yàn)數(shù)據(jù)曲線圖如圖8所示。由圖可知，在測力儀溫升過程中，測力儀結(jié)構(gòu)產(chǎn)生了變形，且變形為非對稱變形，導(dǎo)致傳感器輸出之間并未相互抵消，從而導(dǎo)致測力儀在X、Y和Z三個(gè)方向上產(chǎn)生輸出，單輸出均較小。從測力儀的三向輸出曲線可以看出測力儀在Y方向的輸出誤差最大，最大誤差不到0.5 N，X方向誤差不到0.3 N，Z方向的誤差非常小，可以忽略不記，綜上所述，測力儀具有良好的熱穩(wěn)定性，可以在溫度變化的環(huán)境中進(jìn)行姿控力測試。

圖8 加溫空載標(biāo)定曲線

3 結(jié)束語

本文對小幅值姿控力測試系統(tǒng)熱可靠性進(jìn)行了仿真和實(shí)驗(yàn)研究，通過ANSYS有限元軟件建立了測力儀的參數(shù)化有限元模型，仿真求解了測力儀中傳感器的最高溫度和最大變形量，傳感器最高溫升不到20℃，最大變形量不到0.1 mm，符合設(shè)計(jì)要求；同時(shí)對測力儀進(jìn)行溫度實(shí)驗(yàn)，通過電加熱板對測力儀進(jìn)行加熱，在設(shè)定預(yù)定溫度后對測力儀進(jìn)行溫度加載，對測力儀進(jìn)行溫度加載后測力儀最大輸出誤差出現(xiàn)在Y方向，不到0.5 N，其他方向誤差均小于0.3 N，通過實(shí)驗(yàn)驗(yàn)證了測力儀具有良好的熱穩(wěn)定性。