任彥松 馬志賽 王曉鵬 王加攀 符杰

摘 要：間隙在飛行器折疊舵面中普遍存在，會顯著影響結(jié)構(gòu)系統(tǒng)的動態(tài)行為，給折疊舵面的動力學建模與特性預(yù)示帶來巨大挑戰(zhàn)。針對含間隙飛行器折疊舵面的結(jié)構(gòu)設(shè)計問題，從控制間隙特征、改善工作性能的角度出發(fā)，以間隙為目標尺寸，綜合考慮間隙的產(chǎn)生機理，提出基于矢量環(huán)法的折疊舵面尺寸鏈分析方法，來獲取間隙對相關(guān)尺寸的敏感度以及相關(guān)尺寸對間隙的貢獻度指標。以典型含間隙折疊舵面為例，結(jié)合其實際工作原理構(gòu)建尺寸鏈，分別采用所提出的矢量環(huán)法和計算機輔助的統(tǒng)計分析法進行尺寸鏈分析，兩者所得的各相關(guān)尺寸的敏感度和貢獻度指標相互吻合，驗證了基于矢量環(huán)法進行折疊舵面尺寸鏈分析的可行性，可為折疊舵面的間隙特征設(shè)計優(yōu)化提供技術(shù)支撐。

關(guān)鍵詞：折疊舵面； 間隙非線性； 尺寸鏈分析； 矢量環(huán)法； 設(shè)計優(yōu)化

中圖分類號：TH122 文獻標識碼：A DOI：10.19452/j.issn1007-5453.2023.03.013

基金項目： 航空科學基金 （2020Z009048001）；國家自然科學基金 （12272258）；天津市自然科學基金 （21JCQNJC0360）；天津市青年人才托舉工程（TJSQNTJ-2020-01）

為適裝機載、艦載等平臺，滿足存儲、運輸和發(fā)射過程中嚴苛的空間約束條件，很多飛行器都選用折疊舵面以節(jié)省空間。但由于包含運動副以及加工、裝配誤差和使用磨損等因素，舵面的折疊機構(gòu)處不可避免地存在間隙[1-3]。作為典型的集中式非線性形式，間隙的存在會顯著影響結(jié)構(gòu)系統(tǒng)的動態(tài)特性，使傳統(tǒng)的基于線性理論的氣動彈性計算方法不再適用，給飛行器折疊舵面的動力學建模與特性預(yù)示帶來巨大挑戰(zhàn)，從而嚴重影響了伺服彈性、氣動彈性以及力學環(huán)境的設(shè)計品質(zhì)，甚至引發(fā)飛行故障[4-6]。因此，面向高性能飛行器的研制需求，既需要大力發(fā)展適用于折疊舵面的非線性動力學分析理論及方法，也需要從間隙的產(chǎn)生機理入手，對其進行尺寸鏈分析及設(shè)計優(yōu)化，以有效控制間隙特征，進而從根本上改善折疊舵面的實際工作性能。

在加工和裝配過程中，零部件尺寸及裝配形位均存在偏差，偏差的共同作用會影響最終產(chǎn)品的實際工作性能。所謂尺寸鏈分析，就是在已知零部件的尺寸分布和公差的情況下，通過考慮偏差的傳播和累積，計算得出最終裝配體的尺寸分布和公差的過程[7-9]。通過尺寸鏈分析可得出加工和裝配中各尺寸對目標尺寸的影響程度，進而實現(xiàn)尺寸鏈的設(shè)計優(yōu)化。目前可用于尺寸鏈分析的方法主要有極值法、統(tǒng)計分析法和模擬法（試驗法）等[10-13]。極值法在所有尺寸取公差上偏差或下偏差的情況下進行分析，多用于驗證簡單結(jié)構(gòu)是否滿足裝配要求。統(tǒng)計分析法和模擬法均以計算機作為輔助工具自動生成裝配尺寸鏈[14-16]，并利用相應(yīng)軟件計算得出最終的尺寸鏈分析結(jié)果。綜合來看，現(xiàn)有的尺寸鏈分析方法已能夠?qū)崿F(xiàn)對簡單結(jié)構(gòu)的定性分析，也能夠利用計算機輔助工具對相對復(fù)雜的結(jié)構(gòu)進行定量分析。

現(xiàn)階段，國內(nèi)外學者對于飛行器折疊舵面的間隙非線性動力學問題已開展了大量研究，在動力學模型降階[3，17-18]、氣動彈性分析[19-20]、參數(shù)識別[21-22]，以及振動測試[2，23-24]等方面均取得了顯著進展。然而，針對折疊舵面開展尺寸鏈分析及設(shè)計優(yōu)化的研究還比較少，尚難以滿足航空航天工程實際中對折疊舵面的間隙特征進行有效控制的迫切需求。

綜上所述，本文針對典型含間隙飛行器折疊舵面的結(jié)構(gòu)設(shè)計問題，以間隙為目標尺寸，提出基于矢量環(huán)法的折疊舵面尺寸鏈分析方法，得到各相關(guān)尺寸的敏感度和貢獻度指標，并通過統(tǒng)計分析法對分析結(jié)果進行驗證，可為飛行器折疊舵面的間隙特征設(shè)計優(yōu)化提供技術(shù)支撐。

1 折疊舵面的模型

典型折疊舵面的整體裝配圖如圖1所示，其各部件組成關(guān)系如圖2所示。其中定件與飛行器主體結(jié)構(gòu)相連，動件可繞定件進行轉(zhuǎn)動和折疊，二者通過轉(zhuǎn)軸和鎖銷連接，其中轉(zhuǎn)軸保證舵面動件可折疊，鎖銷保證舵面動件在展開狀態(tài)下可固定。在折疊狀態(tài)下，鎖銷位于動件內(nèi)部，待動件展開后，鎖銷可伸入定件內(nèi)部進行固定。

折疊舵面的工作原理如圖3所示，在舵面展開狀態(tài)下，定動件結(jié)構(gòu)配合面可限制動件逆時針轉(zhuǎn)動，鎖銷限位面可限制動件順時針轉(zhuǎn)動。對于該折疊舵面來說，間隙主要出現(xiàn)在各零部件的裝配部分，因此需要綜合考慮不同位置的間隙，以構(gòu)建合適的尺寸鏈進行分析。

2 尺寸鏈分析方法

2.1 矢量環(huán)法簡介

在加工和裝配過程中，存在若干個必要尺寸來確定整個裝配體的尺寸分布和公差，可將每個尺寸看成矢量，矢量的方向和大小對應(yīng)尺寸的位置和長度。相連的零部件尺寸可視為首尾相接的矢量，由此封閉的尺寸鏈可形成一個矢量環(huán)[25]，如圖4所示，該矢量環(huán)即可作為尺寸鏈分析的計算依據(jù)。

對于圖4中的矢量環(huán)，可建立坐標系，將其中每個矢量用坐標表示。在二維空間中確定矢量方向和大小的坐標有兩個，分別是相位角θ和長度l。根據(jù)各矢量的位置關(guān)系，可基于坐標將封閉環(huán)用多個自變量表示出來，通過進一步計算分析可得出各組成環(huán)對封閉環(huán)的影響程度。上述過程可歸納為基于矢量環(huán)法的尺寸鏈分析。

2.2 矢量環(huán)法示例

下面通過一個簡單示例對基于矢量環(huán)法的尺寸鏈分析流程進行說明。對如圖5所示的二維矢量環(huán)，建立直角坐標系，對其進行尺寸鏈分析。

3 折疊舵面的尺寸鏈分析

3.1 尺寸鏈構(gòu)建

針對圖1所示的含間隙折疊舵面，構(gòu)建尺寸鏈。對于該模型，引起結(jié)構(gòu)非線性現(xiàn)象的主要間隙出現(xiàn)在定動件配合面以及鎖銷限位面，由此引起的動件繞轉(zhuǎn)軸的平面轉(zhuǎn)動問題在實際結(jié)構(gòu)中最為明顯，在這種情況下可采用二維尺寸鏈模型進行分析。

選定最終的目標尺寸為動件繞轉(zhuǎn)軸相對于定件的可轉(zhuǎn)動角度，進一步結(jié)合折疊舵面模型的具體尺寸可構(gòu)建其尺寸鏈。轉(zhuǎn)動角度分為逆時針和順時針兩個方向，各自受不同尺寸的影響，可先對每個方向進行單獨分析，再組合構(gòu)建形成完整的尺寸鏈。

影響動件相對于定件逆時針轉(zhuǎn)動的主要尺寸如圖6所示，各尺寸的符號含義及尺寸公差見表1。其中旋轉(zhuǎn)角度θ1為封閉環(huán)，其公差受到各組成環(huán)偏差的傳播和累積影響，可通過矢量環(huán)法計算得出。

以上是基于矢量環(huán)法進行折疊舵面尺寸鏈分析的全過程，下面采用計算機輔助的分析方法對分析結(jié)果進行驗證。

3.3 分析結(jié)果對比驗證

本節(jié)運用統(tǒng)計分析法對折疊舵面模型進行尺寸鏈分析，借助計算機的運算能力，可得出尺寸鏈中封閉環(huán)對各組成環(huán)的敏感度和組成環(huán)對封閉環(huán)的貢獻度，并以此來驗證基于矢量環(huán)法進行折疊舵面尺寸鏈分析的正確性。

在計算機輔助的尺寸鏈分析過程中，首先基于圖論技術(shù)并結(jié)合連接矩陣，將模型中各零部件尺寸按實際裝配連接起來，構(gòu)建尺寸鏈[28]。在尺寸鏈生成后，進行偏差統(tǒng)計分析，將各零部件的偏差看作獨立的統(tǒng)計分布圖，綜合這些分布圖可預(yù)測出目標尺寸的分布圖[29]。本節(jié)對折疊舵面模型尺寸的統(tǒng)計分析選用系統(tǒng)矩方法，該方法通過計算統(tǒng)計系統(tǒng)矩對數(shù)據(jù)進行描述，可直接采用每個零部件變量的統(tǒng)計分布圖來計算目標尺寸的分布圖，在計算測量偏差時可考慮非線性幾何行為和非正態(tài)變量分布圖。

將模型中各尺寸特征及公差信息導入計算機中，并選用正態(tài)分布的公差概率對尺寸鏈進行統(tǒng)計分析，可得出敏感度和貢獻度指標，見表3。將統(tǒng)計分析法所得結(jié)果與矢量環(huán)法的計算結(jié)果進行對比可知，矢量環(huán)法分析得出的敏感度與貢獻度和統(tǒng)計分析法所得結(jié)果基本吻合，證明了基于矢量環(huán)法進行折疊舵面尺寸鏈分析的可行性。

4 結(jié)論

本文從控制飛行器折疊舵面間隙特征的角度出發(fā)，首先分析了典型折疊舵面中的間隙產(chǎn)生機理，并將間隙特征作為目標尺寸構(gòu)建尺寸鏈，進而基于矢量環(huán)法開展了折疊舵面的尺寸鏈分析，得到了各相關(guān)尺寸的敏感度和貢獻度指標，總結(jié)歸納出各組成環(huán)對目標尺寸的影響規(guī)律，最終通過計算機輔助的統(tǒng)計分析法驗證了上述尺寸鏈分析過程的正確性。上述研究工作是對折疊舵面間隙特征進行尺寸鏈分析的有益嘗試，可為飛行器折疊舵面的間隙特征設(shè)計優(yōu)化提供技術(shù)支撐。

參考文獻

[1]黃銳，胡海巖. 飛行器非線性氣動伺服彈性力學 [J]. 力學進展，2021，51（3）： 428-466. Huang Rui， Hu Haiyan. Nonlinear aeroservoelasticity of aircraft [J]. Advances in Mechanics， 2021， 51（3）： 428-466.（in Chinese）

[2]王博，馬志賽，丁千，等. 基礎(chǔ)激勵下含間隙折疊舵面非線性系統(tǒng)辨識 [J]. 振動與沖擊，2020，39（4）： 122-128. Wang Bo， Ma Zhisai， Ding Qian， et al. System identification of folding rudders with freeplay nonlinearity under base excitation [J]. Journal of Vibration and Shock， 2020， 39（4）： 122-128.（in Chinese）

[3]王強，馬志賽，張欣，等. 基于模態(tài)綜合法的含間隙折疊舵面動態(tài)特性分析 [J]. 航空學報，2020，41（5）： 202-210. Wang Qiang， Ma Zhisai， Zhang Xin， et al. Dynamic character analysis for a folding fin with freeplay nonlinearities based on mode synthesis method [J]. Acta Aeronautica et Astronautica Sinica， 2020， 41（5）： 202-210.（in Chinese）

[4]Yang Ning， Wang Nan， Zhang Xin， et al. Nonlinear flutter wind tunnel test and numerical analysis of folding fins with freeplay nonlinearities [J]. Chinese Journal of Aeronautics， 2016， 29（1）： 144-159.

[5]王斐，冉玉國，譚光輝，等. 全動翼面間隙非線性氣動彈性響應(yīng)分析 [J]. 航空工程進展，2022，13（2）： 18-26. Wang Fei， Ran Yuguo， Tan Guanghui， et al. Nonlinear aeroelastic response analysis of fully moving wing with free play [J]. Advances in Aeronautical Science and Engineering， 2022， 13（2）： 18-26.（in Chinese）

[6]Tang Deman， Dowell E H. Aeroelastic airfoil with free play at angle of attack with gust excitation [J]. Aiaa Journal， 2010， 48： 427-442.

[7]顏世鐺，劉忠明，張志宏，等. 大型起重機減速器尺寸鏈分析[J]. 機械傳動，2015，39（11）： 118-120. Yan Shidang， Liu Zhongming， Zhang Zhihong， et al. Analysis of the dimension chain of large crane reducer [J]. Journal of Mechanical Transmission， 2015， 39（11）： 118-120.（in Chinese）

[8]Chen Hua， Li Xin， Jin Sun. A statistical method of distinguishing and quantifying tolerances in assemblies [J]. Computers & Industrial Engineering， 2021， 156： 107259.

[9]朱彬，于乃江，孟祥海，等. 航空發(fā)動機裝配尺寸鏈公差設(shè)計方法研究[J]. 測試技術(shù)學報，2015，29（2）： 177-184. Zhu Bin， Yu Naijiang， Meng Xianghai， et al. Study on the tolerance analysis of aircraft engine assembly dimension chain[J]. Journal of Test and Measurement Technology， 2015， 29（2）： 177-184.（in Chinese）

[10]孟巧鳳，張林鍹，尹琳崢，等. 基于3DCS的三維尺寸公差的分析與優(yōu)化[J]. 系統(tǒng)仿真學報，2018，30（5）： 1730-1738. Meng Qiaofeng， Zhang Linxuan， Yin Linzheng， et al. Three dimension tolerance analysis and optimization based on 3DCS[J]. Journal of System Simulation， 2018， 30（5）： 1730-1738.（in Chinese）

[11]Chen Hua， Jin Sun， Li Zhimin， et al. A comprehensive study of three dimensional tolerance analysis methods [J]. ComputerAided Design， 2014， 53： 1-13.

[12]Shen Zhengshu， Ameta G， Shah J J， et al. A comparative study of tolerance analysis methods [J]. Journal of Computing and Information Science in Engineering， 2005， 5（3）： 247-256.

[13]容芷君，陳少云，陳奎生，等. 車身產(chǎn)品設(shè)計中的多層次尺寸鏈構(gòu)建 [J]. 計算機輔助設(shè)計與圖形學學報，2012，24（11）： 1394-1398. Rong Zhijun， Chen Shaoyun， Chen Kuisheng， et al.The devel‐opment of multiple hierarchy dimension chain in auto-body de‐sign [J]. Journal of Computer-Aided Design & Computer Graphics， 2012， 24（11）： 1394-1398.（in Chinese）

[14]Gao Zhenbo， Wang Zhexing， Wu Zijian， et al. Study on generation of 3D assembly dimension chain[C]//Proceedings of the 13th CIRP Conference on Computer Aided Tolerancing，2014： 163-168.

[15]郭崇穎，劉檢華，唐承統(tǒng)，等. 基于圖論的裝配尺寸鏈自動生成技術(shù) [J]. 計算機集成制造系統(tǒng)，2014，20（12）： 2980-2990. Guo Chongying， Liu Jianhua， Tang Chengtong， et al. Automatic generation technology of assembly dimension chain based on graph theory [J]. Computer Integrated Manufacturing Systems， 2014， 20（12）： 2980-2990.（in Chinese）

[16]Muholzoev A V， Masyagin V B. Probabilistic calculation of tolerances of the dimension chain based on the Floyd-Warshall algorithm[C]// Proceedings of the 2nd International Conference on Industrial Engineering （ICIE），2016： 959-962.

[17]Wu Zhigang， Yang Ning， Yang Chao. Identification of nonlinear structures by the conditioned reverse path method[J]. Journal of Aircraft， 2015， 52（2）： 373-386.

[18]He Haonan， Tang Hong， Yu Kaiping， et al. Nonlinear aeroelastic analysis of the folding fin with freeplay under thermal environment [J]. Chinese Journal of Aeronautics， 2020， 33（9）： 2357-2371.

[19]Chung K W， Chan C L， Lee B H K. Bifurcation analysis of a two-degree-of-freedom aeroelastic system with freeplay structural nonlinearity by a perturbation-incremental method[J]. Journal of Sound and Vibration， 2007， 299（3）： 520-539.

[20]Ding Qian， Wang Dongli. The flutter of an airfoil with cubic structural and aerodynamic non-linearities [J]. Aerospace Science and Technology， 2006， 10（5）： 427-434.

[21]Abdelkefi A， Vasconcellos R， Marques F D， et al. Modeling and identification of freeplay nonlinearity [J]. Journal of Sound and Vibration， 2012， 331（8）： 1898-1907.

[22]Sun Yukai， Yang Chao， Wu Zhigang， et al. Nonlinear system identification of an all movable fin with rotational freeplay by subspace-based method [J]. Applied Sciences， 2020， 10（4）： 1205.

[23]何昊南，于開平，唐宏，等. 有間隙折疊舵面的振動實驗與非線性建模研究 [J]. 力學學報，2019，51（5）： 1476-1488. He Haonan， Yu Kaiping， Tang Hong， et al. Vibration experiment and nonlinear modelling research on the folding fin with freeplay [J]. Chinese Journal of Theoretical and Applied Mechanics， 2019， 51（5）： 1476-1488.（in Chinese）

[24]Ma Zhisai， Wang Bo， Zhang Xin， et al. Nonlinear system identification of folding fins with freeplay using direct parameter estimation [J]. International Journal of Aerospace Engineering， 2019，30： 3978260.

[25]張蕾. 三維尺寸鏈的偏差分析 [J]. 機械工程師，2015（7）： 105-107. Zhang Lei. Deviation analysis of 3D dimension chain [J]. Mechanical Engineer， 2015（7）： 105-107.（in Chinese）

[26]侯天陽，宋冰. 基于振動水平的離心泵關(guān)鍵尺寸三維公差分析與優(yōu)化[J]. 計量與測試技術(shù)，2021，48（4）： 59-61. Hou Tianyang， Song Bing. 3D tolerance analysis and optimiza‐tion of centrifugal pump based on vibration level[J]. Metrology & Measurement Technique， 2021， 48（4）： 59-61.（in Chinese）

[27]Jiang Qing， Liu Sumei， Sun Jinxin， et al. 3D tolerance analysis of the 1/32 CFETR vacuum vessel [J]. Fusion Engineering and Design， 2020， 160： 111859.

[28]鄭素娟，黃美發(fā)，張奎奎，等. 基于裝配約束的尺寸鏈自動生成方法[J]. 機械設(shè)計與研究，2017，33（2）： 122-128. Zheng Sujuan， Huang Meifa， Zhang Kuikui， et al. Automatic generation of dimension chian based on assembly constraints[J]. Machine Design and Research， 2017， 33（2）： 122-128.（in Chinese）

[29]張志永. 統(tǒng)計公差設(shè)計及檢驗關(guān)鍵技術(shù)研究[D]. 鄭州：鄭州大學，2014. Zhang Zhiyong. Research on key techniques of statistical tolerance design and inspection [D]. Zhengzhou： Zhengzhou University， 2014.（in Chinese）

Dimensional Tolerance Analysis on Folding Fin with Freeplay for Flight Vehicles

Ren Yansong1， Ma Zhisai1， Wang Xiaopeng2， Wang Jiapan2， Fu Jie1

1. Tianjin Key Laboratory of Nonlinear Dynamics and Control， Tianjin University， Tianjin 300350， China

2. China Academy of Launch Vehicle Technology， Beijing 100076， China

Abstract： Freeplay exists widely in folding fins of flight vehicles， which can significantly affect their structural dynamic behaviors and bring great challenges to the corresponding dynamic modeling and characteristic prediction. This paper aims to investigate the structural design problem of a folding fin with freeplay to improve its operating performance by controlling the freeplay characteristics. Firstly， a dimensional tolerance analysis scheme of the folding fin with freeplay is proposed based on the vector loop method to obtain the sensitivity and contribution index of the relevant dimensions regarding to the freeplay. The freeplay is defined as the target dimension and its generation mechanism is comprehensively considered in the proposed method. Thereafter， the dimension chain of a typical folding fin with freeplay is constructed based on the proposed method by considering its actual operating principle. Finally， the proposed method and computer-aided statistical analysis method are subsequently used to analyze the dimension chain， and the sensitivity and contribution index of the relevant dimensions calculated by the two methods are consistent with each other， which validates the accuracy and feasibility of the proposed dimensional tolerance analysis method and provides technical support for the further design optimization of the folding fin with freeplay.

Key Words： folding fin； freeplay nonlinearity； dimensional tolerance analysis； vector loop method； design optimization