基于ANSYS的起落架關鍵部件有限元分析及后處理

孫樂萌 陳閔葉 楊思力

摘 ?要： 起落架是飛機結構中的重要部件之一，不僅承受著飛機與地面接觸時產生的靜、動載荷，還需要吸收飛機在著陸、高速滑跑時所產生的能量?；钊麠U作為起落架的關鍵受力部件，本文以活塞桿為研究對象，通過ANSYS對活塞桿進行應力分析。通過APDL存儲求解結果并導出，利用MATLAB搭建數據批量后處理系統。將活塞桿的有限元分析結果和數據后處理相結合，并對有限元的求解結果進行了深入研究，使得工程技術人員可以針對求解數據進行后續(xù)開發(fā)，為有限元分析提出了新的思路，改善后期數據的處理與交互。

關鍵詞： 起落架;有限元分析;APDL;力學性能;數據分析

中圖分類號： TP319 ? ?文獻標識碼： A ? ?DOI：10.3969/j.issn.1003-6970.2020.06.034

本文著錄格式：孫樂萌，陳閔葉，楊思力，等. 基于ANSYS的起落架關鍵部件有限元分析及后處理[J]. 軟件，2020，41（06）：164168

【Abstract】： Landing gear is one of the important parts of aircraft structure. It not only bears the static and dynamic loads generated when the aircraft contacts the ground， but also needs to absorb the energy generated when the aircraft lands and slides at high speed. The piston rod is the key part of the landing gear. This paper takes the piston rod as the research object， and analyzes the stress of the piston rod by ANSYS. The solution results are stored and exported by APDL， and the data batch post-processing system is built by MATLAB. The finite element analysis results of piston rod are combined with data post-processing， and the solution results of the finite element are studied in depth， so that the engineers and technicians can carry out the follow-up development for the solution data， and put forward new ideas for the finite element analysis， improve the treatment and interaction of the later data.

【Key words】： Landing gear; Finite element analysis; APDL; Mechanical properties; Data analysis

0 ?引言

飛機起落架作為飛機安全著陸時極其重要的承力部件，影響著整個飛機的著陸安全，一旦起落架結構發(fā)生損壞會造成不可估量的損失。利用ANSYS這一款有限元分析軟件準確地進行起落架力學有限元仿真分析是全面評估起落架結構力學性能的有效手段，能夠系統地研究起落架結構的承載能力。但由于ANSYS軟件本身的數據后處理、繪圖等功能有限，常常需要把ANSYS模擬后的結果數據導出以做進一步的處理分析。本文提出了一種新的方法，基于Matlab平臺編寫腳本，將ANSYS APDL模塊計算完成的數據表批量整齊的導出成csv文件，便于后期數據處理交互，提高工作效率。

1 ?起落架結構與關鍵零件建模

1.1 ?飛機起落架結構

飛機起落架主要用于對飛機的起飛、著陸、地面滑跑和地面停放。采用油氣式減震支柱進行減震，通過操縱液壓系統進行起落架正常收放。起落架在飛機著陸時承受巨大沖擊載荷，其振動與沖擊載荷主要靠油氣式減震支柱來吸收，并消耗著陸和滑行時的撞擊能量[1]。工作原理主要是在主活塞與次活塞之間形成的空腔形成緩沖支柱，其內部充滿液壓油與干燥空氣或高壓氮氣，當減震支柱壓縮時，氣體與油液受到壓縮，吸收能量，起到緩沖減震作用。通過機輪剎車裝置吸收水平方向的能量。通過轉彎操縱機構或者差動剎車控制飛機轉彎和地面運動。減緩飛機滑跑時由于跑道不平導致的振動[2]。

1.2 ?零件選擇與建模

減震支柱作為飛機起落架的關鍵零部件，直接影響到整個起落架落地時的安全性能。而減震支柱中又包括外筒與活塞桿兩個主要零部件，考慮到有限元分析的后續(xù)步驟順利進行（有限元網格的劃分、材料和物理特性的定義以及邊界條件的施加等）[3]，選擇活塞桿為本文的研究對象。

具體到活塞桿建模時，由于起落架活塞桿比較復雜，對其整體進行實體分析是一件比較困難的事情，而且計算量也很大。需要得到的是活塞桿的應力情況，同時考慮到計算機資源的限制，需要對活塞桿各個部位以及加載部位進行簡化分析，將對應力影響不大的部位略去。對應力分布只產生較小局部影響的特征，可以忽略[4]。對于一下加載附件和約束構件，在保證載荷正確傳遞的情況下，也可以進行較大簡化，這樣得到保持較高的計算精度但卻大大減少計算量和求解時間的簡化模型，如圖1所示。

2 ?基于ANSYS的活塞桿有限元分析

2.1 ?活塞桿模型導入及參數設計

將所建模型導入ANSYS軟件，以某型號民航飛機為例，其著陸重量大約為66 t。建立一個新的項目樹，進行材料設置。由于起落架結構為飛機的主要受力結構之一，飛機起降頻繁、載重量大，因此對用鋼的抗沖擊性、抗疲勞強度等要求苛刻，本文選用的起落架材料為30CrMnSiNi2A鋼，彈性模量為210 GPa，泊松比為0.28[5]。

2.2 ?活塞桿網格劃分與網格質量檢查

ANSYS有限元網格劃分是進行數值模擬分析至關重要的一步，直接影響著后續(xù)數值計算分析結果的精確性。網格數據的多少將影響計算結果的精度和計算規(guī)模的大小，網格劃分的質量也會影響計算分析結果的準確性[6]。一般來講，網格數量增加，計算精度會有所提高，但同時計算規(guī)模也會增加，所以在確定網格數量時應權衡兩個因素綜合考慮。由于ANSYS19.2具備人工設置網格大小控制疏密分布以及選擇分網算法等，本文采用自由網格劃分對起落架活塞桿進行網格劃分。選中Mesh模塊Details of mesh選項網格參數。其中sizing中的Relevance Center選項設置為Medium，其余采用默認設置。網格劃分圖像如圖2所示。

網格質量檢查是有限元分析過程中的一個重要步驟[7]。通過檢查單元形狀、局部網格細化、改進網格等操作完善網格劃分，使分析結果盡可能準確。ANSYS軟件提供了形狀檢查手段來幫助判斷網格質量，形狀不好的單元可能是計算結果非正常結束，需要清除或修改。

針對活塞桿網格分析的網格質量檢查如圖3所示。通過檢查單元發(fā)現，整個活塞桿共劃分節(jié)點16377個，單元8907個。形狀不好的單元并未出現在所關心的關鍵部位上，Element Quality的均值為0.53098。綜合看來，整體網格質量是良好的。

2.3 ?活塞桿約束及載荷定義

由于活塞桿結構經簡化后作為一個整體進行分析，按照實際結構和力學關系施加約束[8]。結合動力學分析的一般方法，活塞桿與其他構建之間在空間中有三個平動自由度：UX、UY、UZ，三個轉動自由度：RX、RY、RZ。選定Fixed support選項，選擇需要施加固定約束的面。本文需要施加的固定約束面有：活塞桿頂部中心點、輪轂接頭處、活塞桿與扭臂連接處。均施加全約束（即UX、UY、UZ、RX、RY、RZ六個方向全約束）。選擇Pressure選項，同時在Magnitude選項設置壓力，加載模型仿真。

2.4 ?ANSYS與仿真結果分析

對所建模型仿真，其中等效力云圖、等效彈性變形云圖與等效總變形云圖結果如圖4-6所示。

圖4顯示為當前工況下活塞桿的第一應力云圖。從圖中可以看出內筒接觸擠壓區(qū)域出現了較高的第一主應力。圖5顯示為當前工況下活塞桿的第三應力云圖?；钊麠U與外筒接觸區(qū)域出現了較大的應力值。除此之外，下部車架受擠壓作用，也出現

較高的應力值。圖6顯示為當前工況下活塞桿的VonMises應力云圖，在此特殊工況下，活塞桿下端要承受上剎車拉桿的剎車載荷，因此引起活塞桿底端用于安裝上剎車拉桿的耳片拉應力較高。

3 ?ANSYS數據后處理與二次開發(fā)

3.1 ?ANSYS APDL模塊下數據后處理

Mechanical APDL是ANSYS的經典界面，通常所說的ANSYS指的是這個經典界面，workbench ? 是ANSYS推出的一個CAE軟件整合平臺。由于Mechanical APDL與CAE軟件交互困難，所以前期的有限元分析本文使用了workbench進行求解[9]。為了后續(xù)方便批量的處理求解結果的信息表，我們在workbench插入Command調用APDL模塊。

調入APDL模塊后，通過Read Results導入最后一次求解后的結果。打開Sorted nodes模塊，針對我們所需要的節(jié)點坐標、應力分量、主應力和節(jié)點位移分別進行參數設置。由于ANSYS的數據文件保存為lis格式，不方便后續(xù)對數據的批量處理。因此本文通過Matlab軟件編寫腳本，將ANSYS APDL模塊計算完成的數據表批量整齊的導出成csv文件，便于后期數據處理交互，提高工作效率。

3.2 ?m/p腳本

MATLAB是美國MathWorks公司出品的商業(yè)數學軟件，用于算法開發(fā)、數據可視化、數據分析以及數值計算的高級技術計算語言和交互式環(huán)境。本文通過MATLAB軟件對文件格式實現轉化，對多個文件中的海量數據進行批量處理。

采集數據出現次數，并根據節(jié)點序號從小到大排列起來，方便后面獲取每個數的位置。依次識別節(jié)點位移、應力分量、總應力和節(jié)點位移中的數據。通過識別數據集中的空格將數據按照所對應的物理量進行分類，并對應相關的節(jié)點[10]。

運行腳本，依次導入四個lis文件點擊運行，分為坐標操作、應力分量操作、主應力操作、提取S1、S2、S3、位移操作及結果保存五個步驟，在提示音響后數據處理結束，csv文件保存在源數據文件夾下。

3.3 ?數據結果分析

打開生成的csv文件，由左往右依次生成節(jié)點序號、該節(jié)點的X、Y、Z軸坐標、該節(jié)點分別在X、Y、Z軸受力的大小、該節(jié)點所受第一拉應力矢量大小與坐標軸的方向余弦、該節(jié)點所受第二拉應力矢量大小與坐標軸的方向余弦、該節(jié)點所受第三拉應力矢量大小與坐標軸的方向余弦、該節(jié)點發(fā)生位移的適量大小與坐標軸的方向余弦，共26列數據截取部分節(jié)點如表1所示。

4 ?結語

起落架裝置是飛機重要的具有承力兼操縱性的部件，在飛機安全起降過程中擔負著極其重要的使命，是飛機的主要部件之一，其性能的優(yōu)劣直接關系到飛機的使用于安全。減震器作為起落架的重要受力部件，飛機在著陸接地瞬間或在不平的跑道上高速滑跑時，與地面發(fā)生劇烈的撞擊，除充氣輪胎可起小部分緩沖作用外，大部分撞擊能量要靠減震器吸收。本文針對減震器的主要受力部件活塞桿進行研究，并根據已有數據進行建模，由ANSYS對活塞桿進行網格劃分并檢查和優(yōu)化網格質量，分析起落架活塞桿的受力與變形，通過APDL模塊查看和導出分析結果。通過對結果的初步分析可知，內筒接觸擠壓區(qū)域出現了較高的切向應力，活塞桿與外筒接觸區(qū)域以及下部車架收到擠壓作用出現較大的拉應力?；钊麠U底端用于安裝上剎車拉桿的耳片受到較大的拉應力。再通過MATLAB平臺編寫程序，將ANSYS APDL模塊計算完成的lis文件批量整齊的導出成csv文件，便于工程技術人員后期對有限元分析結果的交互處理，簡化工作流程，提高工作效率，改善后期數據的處理與交互。

參考文獻

[1] 姚光生. 飛機起落架關鍵零件強度分析及結構優(yōu)化[D]. 中北大學， 2014.

[2] 姚光生， 秦可偉， 席園. 飛機起落架撐桿結構有限元分析及結構優(yōu)化[J]. 河北農機， 2014（01）： 45-46.

[3] 孫緣. 基于ADAMS和ANSYS的起落架落震仿真及強度優(yōu)化設計[D]. 南京航空航天大學， 2012.

[4] 蔡高存. 飛機主起落架力學建模與分析[J]. 機械工程師， 2018（04）： 56-58+62.

[5] 張成玉. 基于ANSYS的起落架減震支柱設計技術研究[D]. 南京航空航天大學， 2010.

[6] 嵇昆浦. 基于MATLAB-GUI的GNSS基準站網坐標序列數據處理軟件[J/OL]. 測繪地理信息： 1-7[2020-02-22]. http：// kns.cnki.net/kcms/detail/42.1840.P.20200103.1600.030.html.

[7] 高喆， 禹朝帥， 劉釗賓， 等. 基于Python-Matlab的Abaqus后處理技術在柴油機有限元分析中的應用[J]. 拖拉機與農用運輸車， 2017， 44（04）： 46-49.

[8] 張小輝， 王顯會， 張明， 等. 基于MATLAB的螺栓組連接優(yōu)化設計與有限元驗證[J]. 機械制造， 2017， 55（07）： 68-71.

[9] 李亮. 基于ANSYS的起落架扭力臂受力分析及結構優(yōu)化設計[D]. 南京航空航天大學， 2012.

[10] 陳玉振. 主起落架減震支柱的建模方法與強度分析[D]. 南京航空航天大學， 2007.