劉俊增 曾玖海

摘 要：隨著新構(gòu)型直升機技術(shù)的發(fā)展，對直升機槳葉有限元仿真技術(shù)的要求也越來越高。參數(shù)化建模仿真方法能夠降低有限元法仿真技術(shù)的使用難度，對槳葉快速化、精細(xì)化設(shè)計具有重要意義。本文基于ABAQUS二次開發(fā)方法構(gòu)建了一種適用于直升機槳葉根部段的快速建模方法，在有限元建模過程中參數(shù)化建立材料屬性和接觸屬性，完成直升機復(fù)合材料槳葉根部段的快速化仿真分析。利用該方法對某直升機槳根進(jìn)行應(yīng)力計算，與常規(guī)建模方法相比，快速建模方法有效提高了建模效率，保證了仿真結(jié)果的一致性?？焖俳７椒ǚ抡娼Y(jié)果與常規(guī)建模方法仿真結(jié)果相比應(yīng)變誤差不超過1%，快速建模方法與試驗結(jié)果相比應(yīng)變誤差不超過10%，常規(guī)建模與試驗結(jié)果應(yīng)變誤差不超過91%，驗證了該方法的可靠性。

關(guān)鍵詞：直升機； 復(fù)合材料槳葉； 二次開發(fā)； 快速建模方法； 建模效率

中圖分類號：V215.5 文獻(xiàn)標(biāo)識碼：A DOI：10.19452/j.issn1007-5453.2023.05.010

旋翼槳葉是直升機的關(guān)鍵動部件，復(fù)合材料占比成為槳葉先進(jìn)性的主要指標(biāo)之一，先進(jìn)直升機槳葉基本采用復(fù)合材料設(shè)計[1]。Z9、貝爾-214、AH-64D等直升機均在槳葉上使用了玻璃纖維復(fù)合材料或碳纖維復(fù)合材料，以增強槳葉的力學(xué)性能[2]。

復(fù)合材料的應(yīng)用降低了槳葉重量（質(zhì)量），大幅度提升了槳葉飛行壽命，但也給槳葉強度設(shè)計帶來了新的困難，主要表現(xiàn)為其結(jié)構(gòu)復(fù)雜，易產(chǎn)生應(yīng)力集中、分層等問題，尤其是直升機槳根，它主要由蒙皮、大梁、填塊、加強肋、填充物等部件組成。其應(yīng)力分析方法主要包括工程梁算法和有限元法[3]。工程梁算法簡潔高效，能快速進(jìn)行模型重建和應(yīng)力分析，但由于工程梁理論具有局限性，其應(yīng)用范圍受到較大限制。有限元法因其應(yīng)用范圍廣，可對模型進(jìn)行精細(xì)的分析而獲得青睞。在直升機槳葉參數(shù)化建模方面，國外學(xué)者Ceil等[4]使用了參數(shù)化方法建立了單閉室和雙閉室的剖面結(jié)構(gòu)模型。S. Moffatt等[5]在BERP IV槳葉基礎(chǔ)上使用了參數(shù)化方法建立了滿足氣彈剪裁要求的槳葉參數(shù)化模型。國內(nèi)學(xué)者楊建靈等[6]為提高槳葉結(jié)構(gòu)設(shè)計效率，以C形梁為研究對象提出了基于參數(shù)化組件定義的槳葉結(jié)構(gòu)優(yōu)化設(shè)計方法。趙秋華[7]在楊建靈的研究的基礎(chǔ)上進(jìn)一步討論了槳葉接頭填塊、大梁和后緣條的參數(shù)化定義方法。這些方法推進(jìn)了槳葉設(shè)計技術(shù)的發(fā)展。在直升機槳葉參數(shù)化有限元仿真方面的相關(guān)研究較少，有許多學(xué)者通過ANSYS、ABAQUS、NASTRAN等軟件對復(fù)合材料槳葉進(jìn)行了分析，獲得了精確的槳葉強度仿真結(jié)果，但都采用了常規(guī)的有限元手動建模方法，建模效率低，不利于迭代設(shè)計，且分析結(jié)果會因人為因素產(chǎn)生差異[8-11]。

近年來，隨著共軸高速直升機、雙復(fù)合推力高速直升機等新構(gòu)型直升機技術(shù)的發(fā)展，使用的復(fù)合材料占比越來越高，對槳葉有限元仿真技術(shù)的要求也越來越高[12-14]。但是在直升機槳葉實際設(shè)計過程中，往往需要多次迭代進(jìn)行有限元建模計算，有限元法因繁瑣的建模過程而在實際工程應(yīng)用上受到較大限制。參數(shù)化建模仿真方法能夠降低有限元法仿真技術(shù)的使用難度，提高槳葉迭代計算效率，對槳葉快速化、精細(xì)化設(shè)計具有重要意義。

為解決常規(guī)有限元建模方法耗時長、效率低的問題，本文建立了一種適用于直升機槳葉根段的快速建模方法，該方法基于ABAQUS平臺，利用PYTHON開發(fā)接口，實現(xiàn)對槳根模型的快速建模，并對某直升機槳根進(jìn)行應(yīng)力分析，驗證了該方法的可靠性。

1 槳根結(jié)構(gòu)

槳根是槳葉結(jié)構(gòu)受力狀態(tài)最復(fù)雜的部位[3]，槳葉的動載荷和靜載荷通過根段傳遞到槳轂。復(fù)合材料槳根典型幾何模型如圖1所示，橫剖面如圖2所示，主要包括大梁、蒙皮、加強肋、泡沫等。

槳根組成部件中，蒙皮由0/90預(yù)浸碳布和預(yù)浸玻璃布鋪設(shè)成形，可以起到定型和傳力作用，對扭轉(zhuǎn)剛度和擺振剛度貢獻(xiàn)很大。大梁是槳根的主要承力部件，主要使用沿展向的碳纖維/環(huán)氧樹脂基與玻璃纖維/環(huán)氧樹脂基兩種復(fù)合材料。蒙皮和大梁均為有限元仿真過程中的重點關(guān)注部件。

填塊和泡沫使用質(zhì)輕的短切纖維和硬質(zhì)泡沫，主要起到支撐和固定作用。本文中采用的模型各部件材料見表1。

2 基于二次開發(fā)的槳根有限元建模方法

通過編寫內(nèi)核程序自動執(zhí)行人工操作過程中重復(fù)性任務(wù)，可快速實現(xiàn)用戶在前處理階段的仿真建模過程，便于對建模過程進(jìn)行修改和參數(shù)化分析[15-17]。

由于槳根結(jié)構(gòu)復(fù)雜，在有限元建模過程中，建立材料屬性和接觸屬性兩個部分操作耗時最長。本文通過總結(jié)槳根幾何模型的特性，構(gòu)建了一種適用于直升機槳根段的快速建模方法，可以自動建立材料屬性和接觸屬性，有效簡化建模步驟，同時可以避免復(fù)雜操作過程中產(chǎn)生的建模錯誤。

2.1 自動建立材料屬性

槳根組成部件中，蒙皮由多層復(fù)合材料鋪設(shè)而成，在有限元模型中多用殼單元來表示，建立其材料屬性最為繁瑣；而在實際工程設(shè)計中，蒙皮多用參數(shù)化方法來描述，見表2，剖面位置是指沿槳葉展向方向，剖面位置到槳轂中心的距離，扭轉(zhuǎn)角度是指槳葉全局坐標(biāo)系到剖面局部坐標(biāo)系的扭轉(zhuǎn)角度，扭轉(zhuǎn)中心是指在弦線上距剖面前緣的距離，一般為0。起點和終點坐標(biāo)表示鋪層起點與終點距前緣的距離， I、E分別代表上、下翼面，如圖3所示。

內(nèi)核程序首先對蒙皮參數(shù)化數(shù)據(jù)進(jìn)行讀取，通過判斷各層的起點及終點坐標(biāo)位置分析蒙皮是否存在加強布以及加強布與其他鋪層之間的相互關(guān)系。在內(nèi)核程序中定義的常見鋪層類型如圖4所示，其中，圖4（a）為不含加強布的常規(guī)鋪層類型，圖4（b）為不含加強布在相同位置鋪層大于等于兩層的類型，圖4（c）～圖4（f）為加強布在相對其他鋪層不同位置的4種類型，圖4（g）為既含有加強布又存在相同位置鋪層大于等于兩層的類型。

內(nèi)核程序確定蒙皮鋪層類型之后，讀取蒙皮模型的幾何特征，獲取其所有面上的節(jié)點坐標(biāo)，并將每個面沿槳葉展向坐標(biāo)與弦向坐標(biāo)值最大的節(jié)點作為特征點與該面一一對應(yīng)。

分析所有特征點的坐標(biāo)，根據(jù)槳葉展向坐標(biāo)篩選出在目標(biāo)剖面的特征點，通過槳葉厚度方向坐標(biāo)區(qū)分上下翼面，并根據(jù)槳葉弦向坐標(biāo)排序并編號，上翼面坐標(biāo)排序示意圖如圖5所示，其中①～⑤代表層1到層5起點坐標(biāo)，⑥代表前緣坐標(biāo)。

根據(jù)蒙皮鋪層類型確定每一層鋪層的起始點與終止點，并賦予特征點對應(yīng)幾何面相應(yīng)的材料屬性，完成蒙皮材料屬性的自動建立，如圖6所示。

槳根各部件除蒙皮外，均為實體單元。內(nèi)核程序通過讀取表1中各部件與材料的對應(yīng)關(guān)系可直接為其他部件建立材料屬性。

2.2 自動建立接觸屬性

槳根在正常使用時，各部件連接處只有微小的相對位移和變形，因此在有限元建模過程中，采用綁定約束連接各部件接觸面對模型進(jìn)行簡化。綁定約束要求每個表面只能做一次從面，但手動選擇綁定約束的主從表面操作較為煩瑣，且容易發(fā)生過約束錯誤。

為避免以上問題[18-19]，內(nèi)核程序首先通過用戶界面讀取用戶輸入的部件，為方便理解，各部件一般按照從外側(cè)指向內(nèi)側(cè)的順序，如先由剖面上翼面和下翼面指向內(nèi)側(cè)，再從前緣和后緣指向內(nèi)側(cè)。

由于槳根局部結(jié)構(gòu)微小特征較多，為搜索到所有面，內(nèi)核程序首先為蒙皮劃分尺寸較小的網(wǎng)格，并獲取蒙皮上所有節(jié)點[20]。如圖7所示，根據(jù)節(jié)點坐標(biāo)利用ABAQUS內(nèi)置的getClosest函數(shù)搜索剩余部件與蒙皮的接觸面，根據(jù)模型精細(xì)程度設(shè)置好容差，將搜索結(jié)果作為主面，蒙皮所有表面作為從面，建立綁定約束。

本文中內(nèi)核程序根據(jù)用戶輸入順序依次將蒙皮、大梁、肋、泡沫等部件表面作為從面，搜索結(jié)果作為主面，建立TIE約束，完成槳根接觸屬性的建立，其流程如圖8所示。

3 槳根強度分析

使用常規(guī)建模方法和本文的快速建模方法完成某型號槳葉根段的建模。將槳根連接處一側(cè)固定，在另一側(cè)建立參考點與槳根關(guān)注區(qū)域右側(cè)剖面耦合，并在參考點分別施加揮舞力和擺振力，如圖9所示。

圖10和圖11分別為采用快速建模方法和常規(guī)建模方法時揮舞應(yīng)變和擺振應(yīng)變仿真結(jié)果與試驗結(jié)果的對比圖，對比分析可知，快速建模方法仿真結(jié)果與常規(guī)建模方法仿真結(jié)果的揮舞和擺振應(yīng)變誤差不超過1%，快速建模方法仿真仿真結(jié)果與試驗結(jié)果應(yīng)變誤差在-10%～8.3%之間。常規(guī)建模方法仿真結(jié)果與試驗結(jié)果應(yīng)變誤差在-9.1%～9%之間，證明了快速建模方法的可靠性。

4 結(jié)論

針對直升機槳葉根段仿真建模繁瑣、耗時長、不利于迭代等問題，本文基于ABAQUS二次開發(fā)構(gòu)建了槳根快速建模方法，可以自動建立槳根有限元模型的材料屬性和接觸屬性，并通過對比某型號直升機槳根的快速建模方法仿真結(jié)果、常規(guī)方法仿真結(jié)果和試驗結(jié)果，驗證了該方法的可靠性。

在槳根有限元建模過程中，常規(guī)的建模方法在定義材料屬性和接觸屬性時雖然因模型復(fù)雜程度用時不同，但多數(shù)耗時超過1h，而快速建模方法耗時一般不超過5min。快速建模方法有效解決了在槳根迭代設(shè)計過程中有限元建模的效率問題，并通過固定的程序邏輯避免了復(fù)雜操作過程中產(chǎn)生的建模誤差，保證了仿真結(jié)果的一致性。

參考文獻(xiàn)

[1]黃文俊，李滿福.直升機旋翼設(shè)計技術(shù)應(yīng)用現(xiàn)狀及發(fā)展綜述[J].航空制造技術(shù)，2011（17）：32-35. Huang Wenjun，Li Manfu.Application and development of rotor design technology for helicopter[J].Aeronautical Manufactur‐ing Technology，2011（17）：32-35.（in Chinese）

[2]孟雷，程小全，胡仁偉，等.直升機旋翼復(fù)合材料槳葉結(jié)構(gòu)設(shè)計與選材分析[J].高科技纖維與應(yīng)用，2014，39（2）：16-23. Meng Lei， Cheng Xiaoquan， Hu Renwei， et al. Material selection in helicopter composite blade structure design[J]. Hitech-Fiber and Application，2014，39（2）：16-23.（in Chinese）

[3]程小全，高宇劍.含穿透損傷復(fù)合材料槳葉結(jié)構(gòu)靜強度分析[J].失效分析與預(yù)防，2011，6（1）：23-27. Cheng Xiaoquan，Gao Yujian.Static analysis of composite rotor blade with penetration damage[J].Failure Analysis and Preven‐tion， 2011，6（1）： 23-27.（in Chinese）

[4]Celi R. Aeroelastic modeling of swept tip rotor blades using finite elements[J]. Journal of The American Helicopter Society， 1988， 33（2） ： 43-52.

[5]Moffatt S，Griffiths N.Structural optimisat Ion and aeroelastic Tailoring of the berp IV demonstrator blade[C]. American Helicopter Society 65th Annual Forum， 2009：15-20.

[6]楊建靈.基于整體參數(shù)化定義的直升機槳葉結(jié)構(gòu)設(shè)計與優(yōu)化[D].南京：南京航空航天大學(xué)，2013. Yang Jianling.Structure design and optimization for helicopter blade based on global parametric definition[D]. Nanjing： Nanjing University of Aeronautics and Astronautics， 2013.（in Chinese）

[7]趙秋華.直升機槳葉結(jié)構(gòu)參數(shù)化幾何建模與剖面特性分析技術(shù)研究[D].南京：南京航空航天大學(xué)，2015. Zhao Qiuhua. Researches on parametric geometry modeling and profile characteristics analysis for helicopter blade structure[D]. Nanjing：Nanjing University of Aeronautics and Astronau‐tics， 2015.（in Chinese）

[8]楊文志，朱錫，陳悅，等.復(fù)合材料螺旋槳纖維鋪層對強度影響研究[J].艦船科學(xué)技術(shù)，2016，38（3）：85-90. Yang Wenzhi， Zhu Xi， Chen Yue， et al. Research for the ply intensity on the strength of composite propeller[J].Ship Science and Technology，2016，38 （3）：85-90.（in Chinese）

[9]José P B， Christian B， Poul A. Hydro-elastic analysis and optimization of a composite marine propeller[J]. Marine Structures， 2009，23（1）：22-38.

[10]尚麗娜.重型直升機及旋翼槳葉加筋薄壁結(jié)構(gòu)設(shè)計[D].南京：南京航空航天大學(xué)，2009. Shang Lina. Structural analysis of stiffened thin walled rotor blade for heavy-lift helicopter [D].Nanjing： Nanjing University of Aeronautics and Astronautics，2009.（in Chinese）

[11]李雪芹，陳科，劉剛.基于ANSYS的復(fù)合材料螺旋槳葉片有限元建模與分析[J].復(fù)合材料學(xué)報，2017，34（4）：591-598. Li Xueqin， Chen Ke， Liu Gang. Finite element modeling and analysis of composite propeller blade on ANSYS[J]. Acta Materiae Compositae Sinica，2017，34（4）：591-598.（in Chinese）

[12]鄧景輝.直升機技術(shù)發(fā)展與展望[J]. 航空科學(xué)技術(shù)，2021，32

（1）：10-16. Deng Jinghui. Development and prospect of helicopter technology[J].Aeronautical Science & Technology，2021，32（1）： 10-16.（in Chinese）

[13]吳希明，張廣林，牟曉偉. 中國直升機產(chǎn)業(yè)的現(xiàn)狀及發(fā)展建議[J].航空科學(xué)技術(shù)，2021，32（1）： 3-9. Wu Ximing， Zhang Guanglin， Mu Xiaowei. China helicopter industry status and develo- pment proposal[J]. Aeronautical Science & Technology，2021，32（1）：3-9.（in Chinese）

[14]李春華，樊楓，徐明.共軸剛性旋翼構(gòu)型高速直升機發(fā)展研究[J].航空科學(xué)技術(shù)，2021，32（1）： 47-52. Li Chunhua， Fan Feng， Xu Ming. Development of high-speed helicopter with coaxial rigid rotor configuration[J]. Aeronautical Science & Technology， 2021，32（1）：47-52. （in Chinese）

[15]石亦平，周玉蓉. ABAQUS 有限元分析實例詳解[M]. 北京：機械工業(yè)出版社，2006. Shi Yiping， Zhou Yurong. Detailed explanation of ABAQUS finite element analysis example [M]. Beijing： China Machine Press，2006.（in Chinese）

[16]曹金風(fēng)，王旭春.Python語言在Abaqus中的應(yīng)用[M].北京：機械工業(yè)出版社，2011. Cao Jinfeng， Wang Xuchun. The application of Python language in Abaqus[M].Beijing：China Machine Press，2011.（in Chinese）

[17]霍永強，李言，崔蒞沐，等.ABAQUS二次開發(fā)在冷滾打成形參數(shù)化建模與仿真研究中的應(yīng)用[J/OL].機械科學(xué)與技術(shù)：1-12.https：//doi.org/10.13433/j.cnki.1003-8728.20220100 Huo Yongqiang， Li Yan， Cui Limu， et al. Application of ABAQUS secondary development in parametric modeling and simulation of cold roll-beating forming[J/OL]. Mechanical Science and Technology for Aerospace Engineering：1-12.https：//doi.org/10.13433/j.cnki.1003-8728.20220100.（in Chinese）

[18]劉湘云，趙荃.開孔復(fù)合材料層壓板剩余強度分析建模ABAQUS二次開發(fā)[J].機械設(shè)計與制造工程，2019，48（12）： 29-32. Liu Xiangyun， Zhao Quan. The secondary development of ABAQUS for residual strength of open-hole composite laminates[J]. Machine Design and Manufacturing Engineering， 2019，48（12）：29-32.（in Chinese）

[19]吳航.ABAQUS二次開發(fā)在頂管開挖面分析模擬中的應(yīng)用[J].建設(shè)科技，2021（11）：101-104. Wu Hang.Application of ABAQUS secondary development in analysis and simulation of pipe jacking excavation[J]. Construction Science and Technology， 2021（11）： 101-104. （in Chinese）

[20]安康. ABAQUS二次開發(fā)在基于網(wǎng)格的變形幾何體重構(gòu)中的應(yīng)用[D].蘭州：蘭州理工大學(xué)，2020. An Kang.Application of ABAQUS secondary development in mesh based deformation geometry reconstruction[D].Lanzhou： Lanzhou University of Technology，2020.（in Chinese）

Rapid Modeling Method and Stress Analysis on Helicopter Propeller Root Based on ABAQUS Secondary Development

Liu Junzeng， Zeng Jiuhai

China Helicopter Research and Developent Institute，Tianjin 300300，China

Abstract： In order to improve the modeling efficiency of helicopter composite propeller root and ensure the consistency of simulation results， based on ABAQUS secondary development method， a rapid modeling method for helicopter propeller root segment is constructed in this paper， which can quickly establish material properties and contact properties in the process of finite element modeling. This method is used to calculate the stress of a helicopter propeller root. Compared with the conventional modeling method， the rapid modeling method can improve the modeling efficiency effectively. The strain error of the simulation results of the rapid modeling method is less than 1% compared with the simulation results of the conventional modeling method， and less than 10% compared with the experiment results， which verifies the reliability of the method.

Key Words： helicopter； composite blade； secondary development； rapid modeling method； modeling efficiency