基于實測模型的反射面天線裝配定位與補償技術研究

毛 喆 孫瑞峰 郭百森 李 源 朱士琦

基于實測模型的反射面天線裝配定位與補償技術研究

毛 喆1孫瑞峰1郭百森2李 源1朱士琦1

（1. 上海航天電子技術研究所，上海 201109；2. 中國航天科技集團有限公司，北京 100048）

反射面天線裝配質量直接影響著天線的電性能指標，傳統的天線裝配定位對工裝和操作人員均提出了較高的要求，使得天線的裝配質量和裝配效率得不到有效控制。為了提高天線的裝配精度，提出一種基于實測模型的天線裝配定位與補償方法，以實際掃描測量的點云數據為基礎，對點云分割、模型重構、裝配定位等方法展開研究，構建天線實測模型。在三維空間中以實測模型代替理論模型進行虛擬裝配，確定零件修搓量。

模型重構；裝配定位；軸線擬合；點云分割

1 引言

星載反射面天線是衛星有效載荷的重要組成部分，由主副反射面、支撐桿和緊固件等構成[1，2]。天線主反射面和副反射面裝配后的同軸度和安裝面間距離是影響天線性能的關鍵指標，也是裝配質量控制的核心要素。當前天線主副反射面裝配后的同軸度和距離采用專用工藝裝備保證，通過高精度的裝配工裝對主副反射面定位后確定支撐桿的修配或補償量。這種基于模擬量的裝配協調方式具有人工操作誤差累積大、協調路線長、工裝精度要求高等缺點[3]。

隨著探測頻段的上升，對天線主副反射面裝配精度提出了更高的要求。傳統的模擬量裝配協調體系越來越難以滿足現代反射面天線的裝配要求，隨之而來的是數字量裝配協調方法[4～7]。通過數字化設計、數字化測量設備等，將裝配協調部位的尺寸與形位數據直接傳遞到裝配對象，實現裝配對象配合尺寸和形位參數的協調一致。

采用數字化測量設備獲取的高精度點云數據往往需要復雜的分析和處理過程，主要包括點云分割[8，9]、模型重構[10，11]、精度評估[12]、裝配定位等幾個部分。當前這些操作大多數依賴于國外成熟的數據處理軟件完成，如PolyWorks、GeometrixDesign、CATIA、NHRE等[13]。這些商業點云處理軟件雖然具備點云分割和模型重構的能力，但多數需要將理論模型與點云數據對齊后，以理論模型為參照基礎開展分割工作，由于反射面天線理論模型中缺少與點云數據匹配的基準點，采用ICP[14]算法又無法約束理論模型與測量數據兩者之間的相對轉動，因此不適用于反射面天線點云的分割工作。且這些軟件缺乏一定的專用性，處理效率較低，操作過程繁瑣。

針對上述問題，開展反射面天線模型重構與裝配定位技術研究，探索反射面天線裝配的數字量裝配協調方法，構建主副反射面實測模型，定義主副反射面裝配約束條件，在三維空間中對主副反射面實測模型預裝配，根據裝配結果定制支撐桿長度，實現反射面天線的高精度裝配。

2 點云分割

天線主副反射面主要由拋物面、平面和圓柱面三種特征構成。從參數方程和曲率兩個角度出發，分割點云數據。

2.1 平面特征與柱面特征點云分割

平面特征和柱面特征屬于規則幾何特征，具有明確的參數表達式，可通過種子點確定平面、柱面特征初始參數方程。由種子點搜索其附近的鄰域點，計算鄰域點到參數方程的距離，距離滿足一定閾值認為鄰域點是該平面或柱面特征的點，并以鄰域點作為新的種子點，重復向外生長直到分割出滿足目標特征類型的測量點。

2.1.1 平面特征點云分割

平面特征的參數方程可由++=0表達，初始參數方程擬合可由式（1）表示，選取至少4個點作為種子點(x,y,z)，通過特征值法計算出矩陣AA特征值和特征向量，絕對值最小的特征值對應的特征向量即為平面初始參數方程的解。

=0 （1）

其中和分別由下式表示：

2.1.2 柱面特征點云分割

圓柱面方程需要7個參數唯一確定，分別為圓柱中心軸線的方向向量(,,)，圓柱軸線上的點坐標0(,,)，圓柱的半徑，圓柱的方程可由式（2）表示[15]：

函數(,)在初值0處一階泰勒展開可得式（4）：

為求得最小值，對上式求導，令導數值為0可得：

平面特征和柱面特征點云分割具體步驟如下：

a. 選取種子點，平面特征種子點數量不少于3個，柱面特征種子點不少于7個，將種子點存入集合Seed中；

d. 重復執行步驟c開始新一輪的區域增長，將篩選出的特征點存到幾何中；直至所有種子點生長完畢，完成特征分割。

2.2 自由曲面特征點云分割

反射面天線工作面為旋轉拋物面，雖然拋物面具有明確的參數方程，但方程由多個參數構成，擬合求解極為復雜。采用基于曲面平滑度和邊界曲率的區域生長分割法分割拋物面特征。如圖1所示，種子點鄰域內相鄰兩個生長點法向量變化緩慢，在曲面邊界，相鄰兩個生長點曲率之和變化明顯，邊界處的曲率絕對值之和明顯大于曲面內部相鄰兩點曲率絕對值之和。

圖1 自由曲面點云法向及曲率分布特點

由于采用二次曲面擬合的曲面特征與理論特征之間的偏差小，且二次曲面的曲率計算相對簡便，因此本文將局部鄰域點云擬合為二次曲面，根據二次曲面方程估算局部鄰域內任意一點的曲率值。二次曲面可由方程(,)=2++2++表述，構造最小二乘誤差函數為[16]：

為獲取的最小值，對,,,,求一階偏導數，令一階偏導數的值為零可得式（8）。

由式（8）可得出二次曲面參數,,,,的值，局部點云中任意一點的曲率可由下式估算：

反射面工作面點云分割具體步驟如下：

c. 對所有初始種子點執行步驟b后對分割點集匯總，對集合內重復數據點剔除。

3 特征重構

平面和柱面特征重構可由式（1）和式（2）求得，本節對式（1）和式（2）的運用與上一節的差別在于：上一節參與計算的點為選取的種子點，而本節參與計算的點為平面或柱面特征的所有分割點。

用于虛擬裝配的實測模型曲面特征應反映出曲面的尺寸形狀信息，基于拋物面參數方程擬合重構的曲面往往不能精準地反映主副反射面拋物面的變形情況，擬合誤差較大。由于天線工作面的拋物面是開展虛擬裝配的核心，是最需要精確重構的區域，因此采用基于擬合誤差控制的B樣條曲面擬合法對拋物面進行重構[17]。

式（12）中Pi（i=0～n）表示曲線的型值點，曲面控制點的計算需以曲線控制點計算為基礎，如圖2所示，首先對u方向的n+1組點按曲線控制點反算方法，得到各條曲線的多邊形頂點Qij,然后把上述方法計算得到的Qij看作為w方向的m+1組點列，再次按照B樣條反算方法即得到。

拋物面點云重構步驟如下：

b. 基于(+1)×(+1)個點采用式（12）計算出曲線多邊形頂點Q，再次調用式（12）基于Q計算出曲面控制點V,l；

e. 在偏差較大區域增加型值點數量重新執行步驟a，直到輸出重構曲面。

4 天線軸線擬合

主反射面和副反射面的同軸是主副反射面虛擬裝配的重要約束條件之一。由于重構后的反射面為自由曲面，無法通過參數方程計算其軸線方向，因此提出一種基于截面曲線的反射面天線軸線擬合方法。

以重構獲取的反射面的安裝平面作為基準平面P，將其沿平面法向方向生成一組等距面P，將等距面P與重構曲面相交生成一組截交線J。對截交線按照等弧長離散法離散后可獲得每組截交線的離散點P，將每組截交線離散點P擬合為圓C，并將所有圓心C擬合為空間直線作為反射面軸線。

將上式矩陣簡化書寫為=由最小二乘法可計算出擬合圓的圓心公式：

一條空間直線可由式（16）表示為：

由最小二乘原理，對式（16）變形整理后，可推導出式（17）。

上述方法擬合的軸線可能因為雜點的存在導致誤差較大，需剔除雜點。

5 天線裝配定位

固定主反射面，調整副反射面至適宜位置，實現主副反的虛擬裝配。以天線裝配后的要求和裝配連接形式為約束條件，確定副反射面在空間中的位置姿態。

定義主反射面軸線與副反射面軸線重合，可約束副反射面在空間中的2個平移自由度和2個旋轉自由度。定義主反射面和副反射面安裝面中心點之間的距離，可限制副反射面在空間中的一個平移自由度。至此僅剩副反射面繞自身的轉動自由度未得到限制，如圖4所示，主反射面安裝孔中心點應位于副反射面支撐桿安裝面的中垂面上。

圖4 主副反射面裝配定位方式

主副反射面之間的空間位置姿態被唯一確定后，即可在三維模型中確定每一根支撐桿的長度范圍、截面尺寸和形狀等，對這些支撐桿按照分析結果定制、裝配后完成主副反射面的裝配。

6 試驗驗證

為驗證方法的可行性與有效性，基于三維軟件平臺開發反射面天線模型重構與裝配定位模塊，并投產主副反射面試驗件以應用驗證，主反射面口徑450mm，副反射面口徑200mm。采用CWS100雙目視覺掃描系統測量試驗件，該系統測量精度0.02，測量范圍不小于1m，獲取反射面全尺寸點云數據，分割出與實測模型構建相關聯的點云特征，包括反射面工作面、安裝面、支撐桿安裝面等，部分點云特征分割如圖5所示。

圖5 試驗件測量點及點云分割

圖6 試驗件模型重構與偏差分析

對每一分割后的點云特征進行重構，以測量點與擬合面之間的距離作為偏差值，分析重構偏差，如圖6所示。其中主副反射面工作面的重構偏差在±0.1mm范圍內，重構精度較好，而主副反射面安裝孔附近的重構精度較差，主要原因在于該區域獲取的掃描數據缺乏連續性，孔洞數量較多，孔洞尺寸較大。由于該區域對后續裝配定位的影響較小，測量數據仍具有較高的分析價值。

擬合主副反射面軸線，如圖7所示。以天線安裝面和工作面作為輸入條件，輸入等距面間距及截交線離散點間距，將天線工作面解構為一組有序有限離散點集。處理離散點集后，擬合出主副反射面軸線。

圖7 反射面軸線擬合

在三維空間中對主副反射面虛擬裝配。以主副反射面軸線、距離和支撐桿連接形式作為主副反射面裝配的約束條件，實現主副反射面的裝配定位。根據裝配結果對主副反射面支撐桿的截面尺寸、長度等定制，實現主副反射面的高精度裝配。主副反射面試驗件虛擬裝配及實際裝配結果如圖8所示。

圖8 主副反射面虛擬裝配

7 結束語

以反射面天線裝配為研究對象，對基于實測模型的主反射面和副反射面裝配定位方法展開理論研究與實際應用。主要包括點云分割、模型重構、軸線擬合、裝配定位等幾個方面。在上述研究的基礎上基于三維軟件開發原型系統，并投產試驗件驗證本方法的可行性及有效性。結果表明，本方法可以有效地將主副反射面同軸度及距離偏差控制在設計要求的范圍內，有效避免了傳統反射面裝配依賴于工裝定位導致的裝配效率低、裝配精度差等問題。

1 柴艷紅，魏鵬鵬，胡佳，等. 雙反射面天線位姿分析方法及調整技術研究[J]. 航天制造技術，2021(3)：42～44，48

2 施錦文，周衛來，禹旭敏. 星載大口徑太赫茲反射面天線設計與實現[J]. 太赫茲科學與電子信息學報，2019(3)：373～378

3 何勝強. 飛機數字化裝配技術體系[J]. 航空制造技術，2010(23)：32～37

4 Chase K W, Gao J, Magleby S P. General 2-D tolerance analysis of mechanical assemblies with small kinematic adjustments[J]. Journal of Design and Manufacturing, 1995(5)：263～274

5 Sivakumar K, Balamurugan C, Ramabalan S. Concurrent multi-objective tolerance allocation of mechanical assemblies considering alternative manufacturing process selection[J]. The International Journal of Advanced Manufacturing Technology, 2011(5-8)：711～732

6 Ramesh K L, Padmanaban K, Balamurugan C. Optimal tolerance allocation in a complex assembly using evolutionary algorithms[J]. International Journal of imulation Modelling, ISBN, 2016(1)：121～132

7 蘭翔，李曉星，周國鋒，等. 激光跟蹤儀組網測量在大型緊縮場反射面裝調中的應用[J]. 航天制造技術，2015(2)：18～22

8 Wolf，J Prankl，M Vincze. Fast semantic segmentation of 3D point clouds using a dense CRF with learned parameters[C]// 2015 IEEE International Conference on Robotics and Automation (ICRA). 2015

9 鄭書富，李渭，劉莉. 基于超體素聚類的三維點云輪廓特征提取[J]. 蘭州文理學院學報（自然科學版），2018(3)：56～60

10 Ben M A, Louhichi B, Mahjoub M A, et al. Reconstruction of a CAD model from the deformed mesh using B-spline surfaces[J]. International Journal of Computer Integrated Manufacturing, 2019(7)：669～681

11 Peng Weilong, Xu Chao, Feng Zhiyong. 3D face modeling based on structure optimization and surface reconstruction with B-Spline[J]. Neurocomputing, 2016(179)：228～237

12 周石恩. 基于數字孿生的復雜產品裝配建模與精度分析方法[D]. 杭州：浙江大學，2019

13 汪俊. 飛機結構件反求建模關鍵技術研究[D]. 南京：南京航空航天大學，2007

14 Senin N, Colosimo B M, Pacella M. Point set augmentation through fitting for enhanced ICP registration of point clouds in multisensor coordinate metrology[J]. Robotics and Computer-Integrated Manufacturing, 2013(1)：39

15 袁建剛，潘軼. 一種圓柱面擬合方法[J]. 工程勘察，2017(12)：60～64

16 劉勝蘭. 逆向工程中自由曲面與規則曲面重建關鍵技術研究[D]. 南京：南京航空航天大學，2005

17 張麗艷，周來水，蔡煒斌，等. 基于截面測量數據的B樣條曲面重建[J].應用科學學報，2002(2)：173～175

Research on Assembly Location and Compensation Technology of Reflector Antenna Based on Measured Model

Mao Zhe1Sun Ruifeng1Guo Baisen2Li Yuan1Zhu Shiqi1

(1. Shanghai Aerospace Electronic Technology Institute, Shanghai 201109;2. China Aerospace Science and Technology Corporation, Beijing 100048)

The assembly quality of the reflector antenna has a direct influence on the electrical performance of the antenna. The traditional antenna assembly and positioning puts forward higher requirements on both the tooling and the operators, which makes the assembly quality and efficiency of the antenna hardly to be effectively controlled. In order to improve the assembly accuracy of the antenna, an assembly positioning and compensation method based on the measurement data model was proposed. Based on the actual scanning point cloud data, point cloud segmentation, model reconstruction, assembly positioning and other methods were studied, and accordingly the antenna measurement data model was constructed. In the three-dimension space, the actual model is used to replace the theoretical model for virtual assembly to determine the amount of parts rubbing or compensation.

model reconstruct；assembly positioning；axis fitting；point cloud segmentation

V465

A

上海市青年科技英才楊帆計劃（21YF1417300）。

毛喆（1995），助理工程師，航空宇航制造工程專業；研究方向：數字化裝配技術及裝備

2021-07-28