一種特殊衛星推力器矢量測量方法

王偉，郭潔瑛，任春珍，劉浩淼

（1.北京衛星環境工程研究所；2.北京市航天產品智能裝配技術與裝備工程技術研究中心；北京　100094）

?

一種特殊衛星推力器矢量測量方法

王偉1，2，郭潔瑛1，2，任春珍1，2，劉浩淼1，2

（1.北京衛星環境工程研究所；2.北京市航天產品智能裝配技術與裝備工程技術研究中心；北京100094）

摘要：常規衛星推力器的總裝精度測量是通過經緯儀系統準直測量立方鏡工裝實現的。某衛星上出現了一種帶有一個135°拐角的特殊推力器，受其結構及安裝位置的限制，使用常規的經緯儀，現有的精度測量軟件和立方鏡精測工裝無法直接測得此特殊推力器在衛星坐標系下的角度。設計了一種新型精測工裝，并提出了基于該精測工裝的特殊衛星推力器矢量的測量方法。該測量方法已經過驗證和認可，并應用到后續型號的此種特殊推力器矢量測量工作中。

關鍵詞：特殊衛星推力器；矢量測量；新型精測工裝

推力器是為衛星進行變化軌道和調整姿態提供動力的重要控制設備［1］，其在衛星上的總裝精度將直接影響自身控制的精度及燃耗，從而影響航天任務的成敗。因此，在進行推力器的總裝時必須準確測量出其在衛星上的安裝精度并將其調整到設計要求的范圍以內［2］。推力器總裝精度測量至關重要，是衛星的一項關鍵測試工作。推力器在衛星上的總裝精度包括角度和位置度，推力器的角度是指其噴口端面的法線與衛星坐標系所成的角度；推力器的位置度是指推力器噴口端面中心點在衛星坐標系下的坐標值。

目前，推力器在衛星上的總裝精度測量是通過電子經緯儀系統實現的［3］。由于經緯儀無法直接準直測量推力器，所以須借助精測工裝并以其作為測量基準，推力器的精測工裝一般為標準的立方鏡。在進行推力器總裝精度測量之前，首先須在推力器噴口的端面上安裝立方鏡工裝，然后使用經緯儀系統準直測量立方鏡，平行于推力器噴口端面的立方鏡鏡面的法線與衛星坐標系所成的角度就代表推力器在衛星坐標系下的角度；使用經緯儀測量系統布站定標測點可以測得推力器在衛星坐標系下的位置度。

常規推力器的結構都是直線型的，且安裝時指向衛星的下側或外側，可通過經緯儀系統準直測量立方鏡工裝，然后再使用精測軟件對測量數據進行計算就能夠直接獲得推力器在衛星坐標系下的角度。但在某平臺的衛星上出現了一種帶有一個135°拐角的特殊推力器，同時它的噴口又指向衛星內側，受該推力器結構及安裝位置的限制，常規的精測方法無法直接實現其在衛星坐標系下角度的測量。為了準確測量出型號中出現的這種特殊推力器在衛星上的安裝角度，保證其能夠正常工作，本文首次提出了帶有一個135°拐角的特殊衛星推力器矢量的測量方法，成功解決了此種特殊推力器在衛星上安裝角度無法測量的難題，并在型號中進行了推廣應用，效果良好。

1　測量難點分析

衛星上使用的常規推力器的結構均為直線型，其結構如圖1所示。

圖1　常規衛星推力器結構示意圖

常規衛星推力器的精測工裝是金屬立方鏡，其結構如圖2所示。金屬立方鏡的表面進行了鍍膜，以便反射經緯儀的準直光并保證其測量精度。

圖2　常規衛星推力器的立方鏡精測工裝結構示意圖

在進行常規推力器在衛星坐標系下的角度精度測量時，首先將立方鏡工裝通過螺釘安裝固定在推力器的噴口端面上，如圖3所示。

從圖3中可以看出，常規推力器噴口的法線豎直向下。在進行常規推力器在衛星上安裝角度的測量時，使用經緯儀測量系統分別準直立方鏡垂直于推力器噴口端面的兩個相鄰并相互垂直的鏡面，則可直接獲得這兩個被測鏡面的法線在衛星坐標系下的矢量，利用矢量叉乘計算則可間接獲得與這兩個鏡面都垂直的另一鏡面的矢量，再對該矢量進行反余弦計算，即可得到常規推力器在衛星上的安裝角度。

圖3　立方鏡精測工裝在常規衛星推力器上的安裝示意圖

但是，近來在某平臺的衛星上出現了一種帶有一個135°拐角的特殊推力器，其噴口與推力器的主體成135°夾角，如圖4所示。

圖4　特殊衛星推力器結構示意圖

特殊推力器在衛星上的總裝位置如圖5所示，從圖中可以看出該推力器的噴口指向衛星內側。

圖5　特殊推力器在衛星上的總裝位置示意圖

常規立方鏡精測工裝在特殊衛星推力器上的安裝示意圖如圖6所示。

圖6　常規精測工裝在特殊衛星推力器上的安裝示意圖

在對這種特殊衛星推力器進行安裝角度的測量時，衛星是放置在支架車上的，同時該推力器的噴口又指向衛星的內側。在此種工況下如使用常規的總裝精度測量方法，就會因操作空間的限制只能使用經緯儀準直測量到立方鏡工裝的P面（如圖6所示）及與其平行的鏡面，只能直接得到這兩個鏡面的法線在衛星坐標系下的矢量，無法準直測量到與立方鏡工裝P面相鄰并垂直的鏡面，也就無法通過矢量叉乘的方式來計算得到推力器的矢量H。這樣就造成了使用現有的經緯儀系統、精度測量軟件及立方鏡精測工裝的常規方法無法完成此種特殊衛星推力器安裝角度的測量工作。

2　測量方法

使用常規的精度測量方法無法實現此種特殊推力器在衛星上安裝角度的測量工作。為了解決該難題，本文提出了一種基于新型精測工裝的特殊推力器矢量的測量方法，準確地獲得了這種特殊衛星推力器的矢量，進而獲得了該特殊推力器在衛星上的安裝角度。此種測量方法得到了各方的一致認可，并已在型號上進行了推廣應用，取得了較好的效果。

結合特殊衛星推力器的結構及其在星體上的安裝位置，設計了適用于該推力器矢量測量的新型精測工裝，其結構如圖7所示。新型精測工裝由立方鏡和連接桿組成；采用不銹鋼材料一體化加工而成，連接桿與立方鏡上表面之間的夾角為135°；立方鏡表面進行鍍膜；連接桿的另一端通過螺釘安裝固定在推力器噴口端面上，連接桿的B面即為該工裝與特殊衛星推力器的安裝面。

圖7　新型精測工裝示意圖

新型精測工裝在特殊衛星推力器上的安裝示意圖如圖8所示。

圖8　新型精測工裝在特殊衛星推力器上的安裝示意圖

新型精測工裝之所以采用如此結構，就是為了利用工裝與特殊推力器的結構特點，將不可見的噴口法線轉化為方便測量且雙方向可見的立方鏡鏡面法線，建立特殊衛星推力器矢量與新型精測工裝中立方鏡鏡面矢量之間的空間幾何關系，通過對由經緯儀系統測量直接得到的立方鏡鏡面矢量的計算來間接獲得該特殊衛星推力器的矢量，再對該矢量進行反余弦計算，即可獲得特殊推力器在衛星坐標系下的安裝角度。

為了保證帶有一個135°拐角的特殊衛星推力器矢量的測量精度，在設計新型精測工裝時對其加工精度要求如下：連接桿與立方鏡上表面的夾角為135°±2"；連接桿的B面及立方鏡各個鏡面的面形精度均為λ/4（λ=633nm）；立方鏡6個鏡面之間的平行度/垂直度誤差均小于3"；工裝6個鏡面的反射率在可見光波段（λ=0.45～0.65μm）內均大于90%。

帶有一個135°拐角的特殊衛星推力器及新型精測工裝的安裝示意圖如圖9所示，其噴口端面的法線在衛星坐標系下的單位矢量用H來表示。

圖9　特殊推力器及新型精測工裝在衛星上的安裝示意圖

在對該特殊衛星推力器進行精度測量時，首先將新型精測工裝通過螺釘牢固地安裝在推力器的噴口端面上。如圖10所示，設新型工裝立方鏡靠近衛星內側的鏡面為A面，其單位矢量為L；新型工裝立方鏡上表面的單位矢量為M；新型工裝立方鏡下表面的單位矢量為N。

圖10　特殊衛星推力器矢量與

新型精測工裝的立方鏡各鏡面矢量關系圖

從圖10中可以知道，特殊衛星推力器的單位矢量H與新型精測工裝中立方鏡的矢量E相等并且分別與立方鏡鏡面的單位矢量N和單位矢量L成45°夾角，由這四個矢量的空間關系［4］，可知道單位矢量H等于單位矢量N與單位矢量L之和，同時單位矢量N與單位矢量L的模又均為1，所以：

圖11為特殊衛星推力器矢量測量光路圖。如圖11所示，在利用經緯儀系統進行特殊衛星推力器的矢量測量時，首先架設經緯儀T1與T2分別準直衛星基準立方鏡兩個相鄰并垂直的鏡面，然后再架設經緯儀T3與T4分別準直特殊推力器新型精測工裝立方鏡的兩個相鄰并垂直的鏡面，并且經緯儀T2、T3及T4分別與經緯儀T1互瞄，這樣就可以建立由新型精測工裝立方鏡各鏡面的法線構成的坐標系與由衛星基準立方鏡各鏡面的法線構成的坐標系之間的角度矩陣關系，再利用衛星基準立方鏡坐標系與衛星坐標系之間的矩陣關系（已知），就可以傳遞計算出新型精測工裝立方鏡各鏡面法線與衛星坐標系所成的角度，對角度進行余弦計算則可以得到新型精測工裝立方鏡各鏡面法線在衛星坐標系下的單位矢量，再通過式（1）就可以計算出特殊衛星推力器噴口端面法線在衛星坐標系下的單位矢量，對該矢量進行反余弦計算則可獲得特殊推力器噴口法線與衛星坐標系三個軸所成的角度，即其在衛星上的安裝角度。

圖11　特殊衛星推力器矢量測量光路圖

至此就得到了特殊推力器在衛星上的安裝角度，再依據其安裝角度的設計值及精度要求將其調整到合格的范圍內，從而完成了帶有一個135°拐角的特殊衛星推力器的總裝精度測量工作。

3　測量精度分析

帶有一個135°拐角的特殊衛星推力器矢量的測量誤差源自三個方面，第一個方面是由新型精測工裝的加工精度所引起的測量誤差；第二個方面是經緯儀系統及人眼的誤差所引起的測量誤差；第三個方面是建立衛星基準立方鏡與衛星坐標系之間矩陣關系時的測量誤差。

在測量衛星基準立方鏡與衛星坐標系的關系矩陣時，建立衛星基準立方鏡坐標系時的測量誤差約為7″。

以上對特殊推力器矢量的測量精度分析時，都是采用各種誤差因素的最大值，而其實際測量精度會優于12.8″。

4　計算實例

使用本文中提出的矢量測量方法，順利地得到了這種特殊衛星推力器在衛星坐標系下的安裝角度，圓滿完成了該衛星的總裝精度測量工作。

在對該特殊衛星推力器矢量進行測量時，通過經緯儀測量系統及精度測量軟件可以直接得到新型精測工裝立方鏡的上表面和A面的法線與衛星坐標系三個軸所成的角度分別為（0.3746°，90.0213°，90.3740°）、（89.9806°，0.2877°，90.2871°），使用文中提出的計算方法，得到了特殊衛星推力器的單位矢量為H（-0.706852246，0707360737，0.001072445）。再對該矢量進行反余弦計算，則可以得到該特殊推力器與衛星坐標系三個軸所成的角度為（134.9794°44.9794°89.9386°），而其安裝角度的設計值為（135°45°90°）。特殊衛星推力器的精測數據與設計值的符合性較好，也驗證了測量及計算結果的正確性。

目前該衛星已經發射，特殊衛星推力器在軌運行良好，進一步證明了此測量方法的正確性。

5　結論

本文提出了一種帶有一個135°拐角的特殊衛星推力器矢量的測量方法，包括設計了新型精測工裝、提出了推力器矢量及其總裝角度的計算方法，同時也對該測量方法的測量精度進行了分析，并對此測量方法進行了實例驗證。

本文提出的特殊衛星推力器矢量的測量方法，成功解決了因衛星推力器的特殊結構和安裝位置所導致的使用常規的精度測量方法無法進行其總裝角度測量的難題，得到了各方的一致認可并已經應用到后續型號的精度測量中，取得了良好的效果。本文只介紹了一種帶有一個135°拐角的特殊衛星推力器矢量的測量方法，其它類型衛星推力器矢量的測量方法也可采用本文的思路進行。

參考文獻

［1］王憑慧，范本堯，傅惠民.衛星推力器可靠性評估和壽命預測［J］.航空動力學報，2004，19（6）：745-748.

［2］楊再華.航天器總裝精度測量方法分析［J］.航天器環境工程，2007，24（6）：390-392.

［3］王偉，劉笑，郭潔瑛，等.航天器總裝精度測量中一種不規則棱鏡矢量計算方法［J］.航天器環境工程，2013，30 （1）：103-106.

［4］蔣大為.空間解析幾何及其應用［M］.北京：科學出版社，2004.

The Vector Measurement Method of a Special Satellite Thruster

WANG Wei1，2，GUO Jieying1，2，REN Chunzhen1，2，LIU Haomiao1，2
（1.Beijing Institute of Spacecraft Environment Engineering，Beijing 100094；2.Beijing Engineering Research Center of the Intelligent Assembly Technology and Equipment for Aerospace Product，Beijing 100094）

Abstract：Gerneral satellite thruster’s integration precision on the satellite is measured by the electronic theodolite which alignment the assistant tool of cubic prism. There is a special thruster which contains a 135 degree corner on the satellite. Because of the limitation of special thruster’s structure and its assembly position on the the satellite，the use of the electronic theodolite，the current precision measuring software and the cubic prism can’t measure directly the special thruster’s integration angles relative to the satellite coordinate system. In this paper，a new type precision measurement assistant tool was designed and the vector measurement method of the special satellite thruster was proposed，which based on the new type precision measurement assistant tool. The method is verified，approved，and applied in the vector measurement of special thruster in the follow-up satellite models.

Key words：special satellite thruster；the vector measurement；new type precision measurement assistant tool

中圖分類號：V448

文獻標識碼：A

文章編號：1672-9870（2016）02-0030-05

收稿日期：2015-07-13

作者簡介：王偉（1976-），男，工程師，E-mail：wangw0719@sina.com