一種基于ABAQUS的火箭彈部段連接剛度模型修正方法研究

楊磊，郭治斌，尤春艷，楊元華，楊俊波

一種基于ABAQUS的火箭彈部段連接剛度模型修正方法研究

楊磊，郭治斌，尤春艷，楊元華，楊俊波

（四川航天系統工程研究所，四川 成都 610100）

火箭彈按照不同功用被設計成若干個部段，部段間通過各種結構形式連接，其連接部位的剛度下降會顯著影響全彈的固有特性。本文以某型火箭彈為例，首先，通過全彈剛度試驗獲取彈體受載后的變形數據；其次，基于ABAQUS軟件，采用剛度試驗數據結合有限元數值分析的剛度模型修正方法，對火箭彈部段連接剛度進行迭代修正，使得彈體變形量的計算值逐步逼近試驗值；最后，利用修正后的模型對火箭彈前六階固有頻率和振型進行分析，分析結果表明經剛度修正后的固有特性計算值更為接近試驗實測值，從而驗證了該方法在分析此類問題時的有效性。

連接剛度；ABAQUS；模型修正；固有特性

有控火箭彈通常由戰斗部、儀器艙、火箭發動機及尾段等部段組成，多部段的設計有利于日常維護和安裝生產，因此形成了若干部段間連接面。根據部段的載荷、結構形式、尺寸空間等因素，在連接處采用套接、盤式連接、螺紋連接等接頭形式將各部段組裝起來，形成完整的彈身。受接頭形式、配合間隙及預緊力的影響，連接處的剛度與附近部段的剛度相比較，必然會發生明顯的變化，圖1[1]給出了連接面對全彈剛度的局部擾動情況。

圖1 連接面對全彈剛度的擾動

大量試驗數據證明，根據接頭數量、類型和位置的不同，剛度損失可達10%～40%，由于連接面削弱了相鄰彈身的剛度，必將拉低大長細比火箭彈的固有頻率，表1[2]中列舉了一些導彈的一階頻率由于連接面的存在而引起的下降情況。

表1 連接面引起的剛度損失和頻率下降

除此以外，連接面的存在同時改變了全彈剛度分布，必然使得振型形狀、斜率、節點位置發生變化，這些模態參數是飛行姿態控制系統設計的依據。有控火箭彈的飛行控制原理[3]，簡單地說就是依靠陀螺和加速度傳感器感受彈體剛體姿態和加速度，進而控制其飛行。然而當彈體沿垂直于縱軸方向振動時，陀螺感受的姿態中含有彈體彈性振動的附加姿態角，這一附加姿態角將導致火箭彈不能按設計軌道正常飛行。姿態控制系統設計中的一項重要任務就是根據模態參數設計出合理的濾波器將飛行過程中的附加干擾濾掉，確保火箭彈飛行姿態的穩定。因此，產品設計人員必須重視連接剛度對火箭彈系統設計的內在影響，及早掌握其變化規律。

有限單元法（Finite Element Method，FEM）自1943年提出并發展至今，已成為工程數值分析的有力工具[4-5]，特別是隨著電子計算機的廣泛應用和發展，涌現了以ABAQUS、ANSYS為代表的一大批適用于工程領域的通用有限元軟件，它們基于有限元理論，采用數值計算方法，擁有完備的功能模塊和強大的前后處理能力[6]，可以模擬龐大復雜的模型，處理高度非線性問題。應用任意一款通用有限元軟件對部段本身的剛度進行分析均可得到非常準確的結果，但軟件算法原理上都沒有引入部段連接處的空隙、配合面表面粗糙度、接觸狀態等影響因素，使得很難直接對連接剛度進行有效預測。

隨著火箭彈重量和結構尺寸的不斷增加，彈體振動引起的姿態控制問題日益顯著，連接剛度的研究得到了廣泛關注。劉昌洪等[7]采用一種等效彈簧模型對導彈艙段連接接頭進行簡化，并應用這種模型分析了連接剛度對彈體固有特性的影響規律；鄭曉亞等[8]將各被連接件和接頭分開，建立不同粗細的網格模型進行計算，然后再將兩種計算模型通過連接單元結合到一起，降低計算規模；何國軍等[9]從盤式連接的接觸有限元方法出發，計算出部段連接結構的抗壓、抗剪、抗彎和抗扭剛度，采用三維模型描述導彈的整體特性，在部段分離面處用相應剛度的虛梁元來等效模擬連接剛度，分析全彈固有特性，取得了較好效果。

本文針對部段間連接剛度問題，首先通過全彈剛度試驗獲取彈體變形數據，然后基于ABAQUS建立部段及連接部位的三維有限元模型，并結合試驗數據對連接處剛度進行修正，最后利用修正后的連接剛度對彈體固有特性進行分析對比，以此驗證該修正方法的有效性。

1 剛度試驗

全彈剛度試驗的目的是通過測定彈體在外力作用下的變形或轉角，獲取彈體的剛度特性。本文對某型火箭彈進行剛度試驗時，考慮到彈體長細比的值較大，火箭彈如果水平固定，則會因為自重變形引入附加變形量，給剛度數據帶來偏差。為了消除重力的影響并保證所有部段間連接面均處于自由狀態，將火箭彈豎直放置，在發動機結構上選取兩個位置用于夾持固定全彈。將控制艙前端、尾段后端通過轉接工裝與多通道電液伺服加載系統連接，根據試驗件的尺寸規模，分別在彈身兩端加載2000 N單位載荷，試驗方案如圖2所示。

為了獲取彈體結構在加載過程中的變形情況，尤其是部段間的連接剛度，在各部段連接面兩側布置位移測點，如圖3所示。

試驗時，為保證載荷平穩地施加到試驗件上，以100 N為步長逐級加載到單位載荷，所有位移測量均待施加的載荷穩定后進行，彈身在單位載荷作用下的變形數據如圖4所示。

圖2 全彈剛度試驗方案

圖3 剛度試驗位移測點分布

圖4 彈身試驗位移曲線

2 剛度修正方法

為簡化計算，以往對全彈剛度進行分析時，通常忽略連接件，僅對各部段進行建模，然后將部段之間定義為固定連接。這種處理方式人為地提高了部段間的連接剛度，存在較大誤差。基于ABAQUS的連接剛度修正方法則是根據剛度試驗的位移結果來對剛度進行修正，以此獲得更為準確的結果。該方法可以分為兩個步驟：第一步，利用ABAQUS對全彈建立有限元模型，按試驗狀態進行約束和加載，得到全彈的位移數據；第二步，將上一步的計算結果與試驗測得的變形數據進行對比、修正連接剛度，以此擬合位移曲線。

2.1 建立火箭彈有限元模型

參加剛度試驗的某型火箭彈由控制艙、戰斗部、發動機和尾段四個部段組成，除尾段材料為鋁合金外，其余部段材料均為鋼，其材料的性能參數如表2所示。利用Creo對各部段進行幾何實體建模，建模過程中采取必要的簡化，原則是不改變結構的基本幾何特征、傳力路線、剛度分配等。部段建模完成后，將部段按實際連接關系進行裝配，然后采用ABAQUS與Creo之間專用接口程序導入幾何模型。

表2 材料性能參數

2.2 劃分網格

單元類型和網格的疏密程度直接影響運算時間及計算結果精確性，所以必須選擇合理的單元類型及網格劃分方法[10]。雖然采用線性完全積分單元的計算速度較快，但承受彎曲載荷時會出現剪切自鎖，計算精度較差，不利于獲得準確的變形結果。而非協調模式單元克服了剪切自鎖問題，且在厚度方向上只需很少的單元就可得到與二次單元相當的結果，計算成本卻明顯降低，非常適合用于剛度分析，但使用時需確保網格沒有大幅的扭曲，否則會大大降低精度。為確保網格規則，利用Partition Cell工具將各部段劃分成為若干個簡單區域，其中每個部段連接部位需單獨劃分為一個區域，以便對連接剛度進行修正，然后根據每個區域的特點生成結構網格或掃掠網格。經過對全彈進行網格劃分，有限元模型中包含413224個節點和292920個單元，如圖5所示。

2.3 定義作用關系

ABAQUS中常用的相互作用關系有接觸（Contact）、粘接（Tie）、耦合（Coupling）、殼－實體耦合（Shell-to-Solid Coupling）等。按照本文采取的方法思路，不考慮部段之間的連接件，通過修正連接部位的剛度來調整變形量，因此將部段之間的作用關系定義為粘接；固定夾具的剛度遠大于彈身，且與彈身之間為緊配合，將兩者之間的作用關系也定義為粘接。

2.4 施加載荷與邊界條件

在最終提交有限元靜力分析前，還需對彈身施加載荷與邊界條件，原則是力求與試驗狀態一致。在控制艙前端面和尾段后端面的所有節點上施加集中力（Concentrated force），即將單位載荷平均分配到各節點上，載荷方向垂直于彈身軸線；對固定夾具的底座設置固定約束（ENCASTRE）。完成了所有前處理步驟的全彈有限元模型如圖6所示。

2.5 靜力學求解

將全彈有限元模型提交求解器進行求解，全彈變形分布云圖如圖7所示。

圖5 全彈有限元網格

圖6 完成前處理的全彈有限元模型

從全彈變形結果中提取圖3中各測點位置在單位載荷作用下的位移結果，彈身位移曲線如圖8所示。通過對比計算與試驗的位移數據可以發現，沿著彈身固定端到自由端的路徑，在沒有經過部段連接的彈身位移計算值與試驗值較為一致（比如彈身后端從固定處到連接3之間），一旦經過三個部段間連接后位移計算值與試驗值則出現較大的差異，并且連接部位越靠近固定端造成的影響越大（連接2造成的影響大于連接3），這說明對于部段本體的剛度計算已經非常準確，但因部段間連接剛度的誤差導致彈身過于剛硬。

2.6 剛度修正

從圖8位移曲線的趨勢可以看出，因計算時部段間連接剛度偏大，造成彈身位移值偏小。按照剛度與變形成反比的思路，對各部段連接剛度進行修正，經過多輪反復迭代，位移計算曲線與試驗曲線比較吻合，修正前后的連接剛度如表3所示，位移曲線如圖9所示。

表3 修正前后連接剛度系數

3 固有特性分析

部段間連接面的存在使得連接剛度和阻尼發生非線性變化，從而導致結構固有特性的改變[11]。結構固有特性包括了結構的固有頻率和振型，固有特性分析是動力學分析的基礎，通過固有特性分析可以幫助設計人員確定火箭彈的固有頻率和振型，從而進行姿態控制設計和儀器部件安裝布局，避免發生頻率耦合。

圖7 全彈變形分布云圖

圖8 彈身位移曲線

圖9 彈身修正位移曲線

利用前面建立的有限元模型對全彈在自由狀態下的固有特性進行分析，由于低階模態對火箭彈的振動特性影響較大，所以選擇前六階的固有頻率和振型進行研究。為了便于比較，分別計算連接剛度修正前后的全彈固有頻率，結果如表4所示。

前六階振型均為彎曲振型，典型振型如圖10所示。

表4 固有頻率計算結果

圖10 典型振型

4 結論

本文根據全彈剛度試驗數據，采用有限元軟件ABAUQS對部段連接剛度進行了修正，進一步利用修正后的剛度參數對全彈固有特性進行了分析，其結果表明：

（1）部段之間的連接剛度與全彈剛度及固有頻率成正比，隨著連接剛度的降低，頻率也隨之降低；

（2）連接剛度對全彈固有特性的影響隨著頻率階數的升高，表現得愈加明顯；

（3）如果不考慮連接剛度的影響，就等于變相提高了全彈剛度，由此計算得到的固有頻率與模態試驗結果相比較，一彎、二彎、三彎振型對應頻率的誤差分別為5.9%、11.6%、20.1%，存在較大的偏差，因此在進行有限元計算時必須對連接剛度慎重考慮；

（4）在對連接剛度進行修正后，計算結果與試驗相比，一彎、二彎、三彎振型對應頻率的誤差分別為2.3%、4.1%、9.4%，說明計算結果可信；

（5）由計算結果來看，根據剛度試驗數據并采用ABAQUS對連接剛度進行修正的思路正確，對于連接面較多的火箭彈或者其他類似結構的仿真具有一定的工程實用價值。

[1]張延珍. 戰術導彈連接剛度的研究[J]. 強度與環境，1992（2）：28-34.

[2]余旭東，趙育善. 飛行器結構動力學[M]. 西安：西北工業大學出版社，1998.

[3]龍樂豪導彈與航天叢書——總體設計[M]. 北京：中國宇航出版社，2001.

[4]Courant R. Variational Method for Solutions of Problems of Equilibrium and Vibrations[J]. Bull. Am. Math. Soc.，1943（49）：1-23.

[5]王勖成，邵敏. 有限單元法基本原理和數值方法[M]. 北京：清華大學出版社，1997.

[6]Chen Yuhmin，TsaoTienheng. A structured methodology for implementing engineering data management: an information model for cooperative product development[J]. Robotics and Computer- Integrated Manufacturing，1998，14（4）：275-296.

[7]劉昌洪，成楚之. 戰術導彈艙段聯結剛度研究[J]. 現代防御技術，1997（3）：24-32.

[8]鄭曉亞，張鐸，姜晉慶. 連接剛度對導彈固有特性的影響[J]. 彈箭與制導學報，2005（4）：667- 669.

[9]何國軍，竺潤祥，尹云玉，等. 考慮連接剛度的導彈固有特性分析[J]. 強度與環境，2006（1）：14- 17.

[10]石亦平，周玉蓉. ABAQUS有限元分析實例詳解[M]. 北京：機械工業出版社，2006.

[11]王尚文. 航天器結構動力學[M]. 西安：西北工業大學出版社，1995.

Research on a Rocket Structure Connection Stiffness Model Updating Method Based on ABAQUS

YANG Lei，GUO Zhibin，YOU Chunyan，YANG Yuanhua，YANG Junbo

( System Engineering Institute of Sichuan Aerospace, Chengdu 610100, China )

The rocket consists of several structures according to different functions. And the structures have various forms of connections with different connection stiffness. The stiffness degradation of the connection significantly affects the inherent characteristics of the whole rocket. In this paper, a certain model of rocket is taken as an example. First of all, the deformation data of the rocket body after being loaded is obtained through the rocket’s structure connection stiffness test. Then, the stiffness model updating method based on the stiffness test data and the finite element analysis is adopted to iteratively correct the stiffness of the connection parts through the ABAQUS finite element software, so that the calculated value of the rocket body deformation gradually approximates the test value. Finally, the updated model is used to analyze the first six orders of the rocket’s natural frequencies and the mode shapes. The result shows that the calculated value of the inherent characteristics after the stiffness updating is closer to the test value, which verifies the validity of the method in analyzing such problems.

connection stiffness；ABAQUS；model updating；inherent property

TE415；TJ7

A

10.3969/j.issn.1006-0316.2021.12.006

1006-0316 (2021) 12-0043-07

2020－11－30

楊磊（1979－），男，重慶人，碩士，高級工程師，主要研究方向為力學仿真與試驗，E-mail：leiyangcqu@163.com。