攔阻著艦環境下的攔阻索索力識別方法研究

王陶 黃丹羽 陳國平 何歡

摘要：攔阻索的安全性是艦載機在航母上安全著艦的關鍵因素之一。針對艦載機攔阻著艦過程，提出了一種通過艦載機著艦攔阻過程中機體加速度求解攔阻索張力的方法。首先，根據艦載機攔阻鉤掛索后的受載情況，確定攔阻索載荷的傳力路徑，實現攔阻索索力的間接識別。然后，通過考慮攔阻索的彎折波特性，對索力識別方法進行了改進和優化。最后，采用多體動力學軟件MSC.ADAMS建立了離散攔阻索模型，對基于加速度的索力識別方法進行了分析和驗證。分析結果顯示，攔阻索索力識別的平均誤差可以控制在5%以內。

關鍵詞：艦載機；攔阻索；多體動力學；索力識別；彎折波

中圖分類號：V222文獻標識碼：ADOI：10.19452/j.issn1007-5453.2022.01.005

航空母艦是現代遠洋艦隊的核心，而艦載機則是航空母艦形成有效戰斗力的關鍵因素[1-3]。由于航空母艦甲板長度有限，艦載機通常需要借助攔阻系統進行著艦回收，而攔阻索是攔阻系統的關鍵部件。攔阻索的安全性不僅直接決定艦載機是否能安全著艦，也是保證航空母艦正常運行的核心要素之一。為保障艦載機的著艦安全性，需要及時發現攔阻索損傷特征并及時更換。因此，準確地確定攔阻索損傷情況對保證航母編隊的戰斗力具有非常重要的意義。

在橋梁領域，通過監測拉索索力來判斷索結構的安全性已有大量成功應用案例，取得了良好的社會和經濟效益[4-7]。在橋梁索力識別問題中，拉索長度通常被認為是固定的，所以橋梁拉索可認為是一種定常系統，可以根據實測拉索的動響應數據，再結合拉索動力學模型來識別拉索的索力，從而判斷拉索的力學狀態。

與橋梁拉索不同，在艦載機攔阻系統的攔阻過程中，阻攔索載荷是由于艦載機攔阻鉤的拖拽而產生的，所以需要監測索力的阻攔索并不是固定長度的繩索，而是會隨著艦載機攔阻減速，不斷從甲板下的攔阻系統中被抽出。這類時變系統模型與傳統的橋梁拉索力學模型之間存在顯著區別。此外，艦載機攔阻鉤在與攔阻索嚙合時，攔阻索會受到很大的沖擊作用，繩索在攔阻過程中并非始終保持軸線是平直的狀態，而是由直線預緊狀態變為三角形布局的應力傳播狀態，這會對索內張力造成較大影響[8]。針對攔阻索這一特點，張新禹[9]建立了考慮彎折波效應的攔阻索模型，開展了攔阻索應力波研究。Xu等[10]通過將繩索離散為等長圓柱剛體，建立了繩索的多體動力學模型，開展了由于撞擊產生的應力波的傳播特性研究。羅青等[11]建立了一種計及彎折波的攔阻系統模型。張萍等[12]詳細分析了彎折波的形成及傳播過程，建立了一種離散彎折波構型的數學模型，對攔阻過程最優軌跡進行了分析。沈文厚等[13-14]根據絕對節點坐標法建立了考慮彎折波效應的繩索動力學模型，研究了不同工況條件下的繩索應力變化規律。然而，目前在索力識別方面，國內外主要的研究對象為橋梁拉索，在艦載機攔阻過程中攔阻索索力識別方面，國內外還鮮有公開文獻涉及。

針對艦載機著艦攔阻過程，本文提出一種基于艦載機加速度的攔阻索索力識別方法[15]。首先，本文根據艦載機攔阻著艦過程的受載特點，建立了機體動力學方程，利用采集到的機體加速度數據并結合攔阻索傳力路徑完成攔阻索索力識別。然后，考慮繩索彎折波的影響，對基于加速度的索力識別方法進行了改進和優化。最后，通過一個二維平面攔阻過程數值算例對本文所提方法的有效性進行了說明。

1攔阻索索力識別

1.1艦載機著艦受力分析

1.2攔阻索索力

在艦載機攔阻著艦的過程中，攔阻索通過傳遞攔阻機構提供的攔阻力實現艦載機的攔停。在對攔阻索載荷進行研究時，若只關注攔阻索的傳力路徑，則可以將繩索簡化為直線。假設攔阻鉤撞擊攔阻索后，攔阻索在運動過程中總是保持直線狀態，則可以認為攔阻索內部張力的方向會一直由鉤、索嚙合點指向滑輪處。于是，艦載機在攔阻著艦過程中所受攔阻力FLZ與攔阻索內部張力T之間的幾何關系如圖2所示。

在圖2中，M為艦載機尾鉤與索之間的接觸點；TL、TR表示攔阻索在M點兩側的繩索內部張力；FLZ為攔阻索對機體的攔阻力；γ為FLZ與甲板之間的夾角；βL、βR表示攔阻索在M點兩側的繩索相對于掛索前位置的夾角。xd代表飛機飛行的方向。為了便于說明問題，提出如下假設：（1）假設攔阻索系統對稱，且攔阻鉤垂直鉤住攔阻索中心位置；（2）假設在攔阻過程中攔阻索各處張力相等；（3）假設整個攔阻索采用同一種材質。

3數值算例

3.1模型簡介

利用MSC.ADAMS進行攔阻索建模。首先，將繩索分割為多段圓柱體，每個圓柱體的長度為ΔL，攔阻索的幾何及材料參數見表1。然后，采用柔性連接單元（即六自由度彈簧-阻尼器）將相鄰圓柱體進行連接，如圖4所示，柔性連接單元的剛度和阻尼系數見表2。最后，建立定滑輪以及地面支撐架，并通過Contact定義攔阻索與它們之間的接觸關系。最終的攔阻系統模型如圖5所示。

為了說明問題方便，本算例假設艦載機在攔阻著艦運動過程中攔阻鉤與攔阻索掛點始終保持在攔阻索所在的平面內運動，并且在攔阻鉤掛索后攔阻鉤與攔阻索之間不存在相對滑動的現象。所以，在進行二維平面內的攔阻索動力學分析時，可以將艦載機與攔阻鉤作為一個整體部件進行剛性化建模。建模過程中需要保證其質量與初始條件與原艦載機模型相同，攔阻鉤部分保留其原有的幾何構型，通過接觸屬性定義其與攔阻索之間的碰撞作用，實現對攔阻索的抽出和拖拽。等效艦載機模型的相關參數見表3。

3.2攔阻索索力識別分析

通過仿真計算，艦載機可以在3s左右完成整個攔阻過程。提取圖4所示的離散攔阻索模型中的柔性連接單元軸向力等效替代索力值，圖6給出了攔阻索索力時域曲線，該索力曲線可以作為后續索力識別的參考值。

基于離散繩索模型進行攔阻仿真，圖7～圖9分別給出了艦載機速度、位移以及加速度變化曲線。從圖中可以看出，經2.6s左右，艦載機位移約為95m。

基于艦載機實時位置可以確定攔阻索張角隨時間變化的曲線，如圖10所示。利用本文提出的基于加速度的索力識別方法處理后，所得的攔阻索索力識別結果曲線與參考值之間的對比如圖11所示。從圖中可以看出，在攔阻的初期，尤其是攔阻鉤碰撞攔阻索瞬間，通過加速度識別方法得到的索力與實測索力有較大誤差。但在隨后的攔阻過程中，索力識別結果都可以與實測結果較好地符合。對索力識別結果分掛索早期與掛索中后期兩個階段進行誤差分析，這兩段時間歷程上的最大誤差與平均誤差見表4。

圖12為仿真分析得到的攔阻索不同階段的變形情況。可見，在實際攔阻過程中，攔阻索無法保證掛點與滑輪間繩索滿足直線構型的假設。由于攔阻鉤的撞擊，繩索產生彎折波效應對攔阻索掛點處的張角有較為明顯的影響。從圖12可以看出，這種彎折效應在攔阻初期的影響較大，兩端繩索還保持原位置的直線狀態時，鉤與索的掛點已向前運動了一段距離，此種狀態下，機體所受攔阻力與繩索內部張力之間的夾角與繩索直線構型假設下的角度有明顯不同。

因此，通過引入繩索彎折波特性重新對攔阻過程攔阻索索力進行識別。根據彎折波的波速與飛機位移可以重新構建考慮彎折波效應的繩索構型，從而得到更加精確的攔阻力與繩索內部張力之間的幾何關系。

圖13和圖14分別為采用直線假設的攔阻過程繩索構型和考慮折波效應的攔阻過程繩索構型。對比修正前后的攔阻索幾何構型可以看出，采用直線構型假設時，攔阻初期的繩索幾何構型與實際情況有較大的差異，從而導致在求解攔阻力與繩索內部張力之間的夾角時會引入較大的誤差。所以，考慮繩索彎折效應十分必要。在此基礎上，對掛索初期索力進行進一步的優化識別，其識別結果如圖15所示。

從圖15索力識別結果可以看出，攔阻鉤撞擊攔阻索初期的攔阻索張力識別的最大誤差約為80kN。將考慮彎折波效應前后攔阻過程中0～0.5s時間段的索力識別情況進行對比，見表5，相比直線構型假設下的約124kN的最大誤差，考慮彎折波效應后的識別方法的識別精度得到了顯著的提高，整個攔阻早期階段的平均誤差與誤差百分比均大幅降低。

4結論

針對攔阻工作狀態下攔阻索索力難以直接測試的問題，本文提出了通過艦載機著艦攔阻過程中機體加速度求解攔阻索張力的方法，然后引入繩索彎折波效應的影響，對識別方法進行了改進和優化。主要結論如下。

（1）通過監測機體加速度的方法可以有效識別攔阻索在艦載機攔阻著艦過程中沿軸向的索內張力。

（2）通過考慮攔阻索彎折波效應，可以有效降低攔阻初期階段的攔阻索索力識別誤差。而攔阻中后期仍可保留嚙合點與滑輪間繩索的直線假設。結果顯示，該階段的平均誤差仍可控制在3%左右。

參考文獻

[1]劉廣，顏世偉，李海旭，等.國外航母航空保障技術發展現狀[J].航空科學技術，2020，31（2）： 1-11. Liu Guang，Yan Shiwei，Li Haixu，et al. Aviation support technologydevelopmentforforeignaircraftcarrier[J]. Aeronautical Science & Technology，2020，31（2）： 1-11. （in Chinese）

[2]劉小川，王彬文，白春玉，等.航空結構沖擊動力學技術的發展與展望[J].航空科學技術，2020， 31（3）： 1-14. Liu Xiaochuan，Wang Binwen，Bai Chunyu，et al. Progress and prospect of aviation structure impact dynamics[J]. Aeronautical Science & Technology， 2020， 31（3）： 1-14. （in Chinese）

[3]董冬妮，危虹，蔣慶喜，等.美國海軍艦載機通用型輔助動力系統的PBL分析及啟示[J].航空科學技術，2020，31（12）： 74-79. Dong Dongni，Wei Hong，Jiang Qingxi，et al. Analysis and enlightenment of PBL for the APU installed in the US Navy’s carrier based aircraft[J]. Aeronautical Science & Technology， 2020， 31（12）： 74-79. （in Chinese）

[4]Yan B，Chen W，Yu J，et al. Mode shape-aided tension force estimation of cable with arbitrary boundary conditions[J]. Journal of Sound and Vibration，2018，440：315-331.

[5]Huang P，Wang X，Wen Q，et al. Active control of stay cable vibration using a giant magnetostrictive actuator[J]. Journal of Aerospace Engineering，2018，31（5）：04018074.

[6]Xiong W，KongB，TangP，etal. Vibration-based identification for the presence of scouring of cable-stayed bridges[J]. Journal of Aerospace Engineering，2018，31（2）：04018007.

[7]Liu M，Zheng L，Zhou P，et al. Stability and dynamics analysis of in-plane parametric vibration of stay cables in a cable-stayed bridge with superlong spans subjected to axial excitation[J]. Journal of Aerospace Engineering，2020，33（1）：04019106.

[8]宋旦鋒，盧永錦.彈性鋼索中縱波與彎折波的耦合作用研究[J].船舶工程， 2017， 39（1）：64-67. SongDanfeng，Lu Yongjin.Researchoncouplingof longitudinal wave and kink wave in the elastic cable[J]. Ship Engineering， 2017， 39（1）：64-67. （in Chinese）

[9]張新禹.阻攔索的動力學特性分析及仿真研究[D].哈爾濱：哈爾濱工程大學， 2011. Zhang Xinyu. Dynamic analysis and simulation of arresting cable[D]. Harbin： Harbin Engineering University， 2011. （in Chinese）

[10]Xu B，Dong P，Zhang J，et al. A three dimensional approach to model steel wire ropes used in high energy absorber apparatus[J]. International Journal of Solids & Structures，2014，51（25-26）：4280-4293.

[11]羅青，馮蘊雯，馮元生.基于彎折波的艦載機攔阻著艦動力學分析及仿真研究[J].機械強度， 2009， 31（4）：543-547. Luo Qing， Feng Yunwen， Feng Yuansheng. Dynamic analysis and simulation of carrier aircraft arrested deck-landing based on kink-wave[J]. Journal of Mechanical Strength， 2009， 31（4）： 543-547. （in Chinese）

[12]張萍，金棟平.計及彎折波的艦載機攔阻過程控制[J].航空學報， 2011， 32（11）： 2008-2015. Zhang Ping， Jin Dongping. Control of arresting process for carrier aircraft considering kink-wave[J]. Acta Aeronautica et Astronautica Sinica， 2011， 32（11）：2008-2015. （in Chinese）

[13]沈文厚，丁祥，郭叔偉.攔阻索應力的仿真研究[J].科學技術與工程， 2015（14）： 224-229. Shen Wenhou， Ding Xiang， Guo Shuwei. Stress simulation of the purchase cable[J]. Science Technology and Engineering， 2015（14）： 224-229. （in Chinese）

[14]沈文厚，趙治華，任革學，等.攔阻索沖擊的多體動力學仿真研究[J].振動與沖擊， 2015， 34（5）：73-77. Shen Wenhou， Zhao Zhihua， Ren Gexue， et al. Multi-body dynamic simulation of impact on cross deck pendant[J]. Journal of Vibration and Shock， 2015， 34（5）： 73-77. （in Chinese）

[15]黃丹羽.艦載機攔阻過程攔阻索分析與索力識別方法研究[D].南京：南京航空航天大學， 2020. Huang Danyu. Research on arresting cable and cable force identification under arresting condition[D]. Nanjing： Nanjing University ofAeronautics andAstronautics， 2020. （in Chinese）

Research on Identification Method of Arresting Cable Force for Carrier Based Aircraft Arresting Landing

Wang Tao，Huang Danyu，Chen Guoping，He Huan

State Key Laboratory of Mechanics and Control of Mechanical Structures，Nanjing University of Aeronautics and Astronautics，Nanjing 210016，China

Abstract： The safety of arresting cable is one of the key factors for carrier based aircraft’s safe landing on aircraft carrier. Aiming at the arresting and landing process of carrier based aircraft， a cable force identification method based on acceleration of carrier aircraft is proposed. Firstly， according to the dynamic relationship between the carrier aircraft and the arresting cable， the transmission path of arresting cable load is determined to realize the indirect identification of arresting cable force. Then， the cable force identification method is improved by considering the kink wave characteristics of the arresting cable. Finally， a discrete arresting cable model is established by using multi-body dynamics software MSC. ADAMS， and the acceleration based cable force identification method is analyzed and verified. The results show that the average error of arresting cable force identification can be controlled within 5%.

Key Words： carrier aircraft； arresting cable； multi-body dynamics； identification of the cable force； kink wave

3118500338200