邱蕾蕾,王雪,許光群,謝金標(biāo)
(國(guó)營(yíng)蕪湖機(jī)械廠,安徽 蕪湖 241007)
某型飛機(jī)活動(dòng)座艙蓋主要由前中后弧框、側(cè)骨架、通條壓板和有機(jī)玻璃組成[1],主要作用是為飛行員進(jìn)出座艙蓋提供通道及良好的飛行視野,同時(shí)在飛行中保護(hù)飛行員免受外界氣流的影響[2-3]。飛機(jī)在服役過(guò)程中不僅要承受各種循環(huán)加載載荷,同時(shí)還需抵擋外來(lái)物(如鳥(niǎo)撞等)帶來(lái)的沖擊力[4],這些性能要求座艙蓋結(jié)構(gòu)有足夠的強(qiáng)度來(lái)滿足飛行安全。
經(jīng)外場(chǎng)使用統(tǒng)計(jì),飛機(jī)在飛行過(guò)程中座艙蓋骨架側(cè)型材、前弧支臂及中弧拉桿處等部位裂紋故障的發(fā)生概率較高,屬于危險(xiǎn)薄弱部位,裂紋的擴(kuò)展對(duì)飛行安全產(chǎn)生較大的影響。針對(duì)裂紋擴(kuò)展規(guī)律,Sih[5]和Hartranft[6]研究了復(fù)變量方法,結(jié)果表明裂紋的強(qiáng)度因子與型材的厚度有著極大的關(guān)聯(lián);對(duì)于板材在受到拉伸和彎曲載荷時(shí),裂紋失穩(wěn)后延續(xù)擴(kuò)展準(zhǔn)則,Wynn等[7]進(jìn)行深入研究,結(jié)果表明“能量型”準(zhǔn)則與實(shí)際情況更加吻合;Lanciotti等[8]對(duì)鋁合金材質(zhì)的貫穿性裂紋在拉伸和彎曲應(yīng)力共同作用下的擴(kuò)展速率進(jìn)行了研究,結(jié)果表明裂紋擴(kuò)展速率與拉伸與彎曲應(yīng)力的占比有關(guān)。
本文針對(duì)座艙蓋骨架側(cè)型材危險(xiǎn)薄弱部位進(jìn)行了損傷容限分析,主要包括未修理下側(cè)型材損傷容限、采用1 mm厚鋼板加強(qiáng)后側(cè)型材損傷容限,以及1.5 mm厚鋼板加強(qiáng)后側(cè)型材損傷容限。研究結(jié)果表明,采用鋼板加強(qiáng)后,座艙蓋骨架側(cè)型材使用可靠性得到了極大的提升。
某座艙蓋側(cè)型材骨架通過(guò)多個(gè)鎖環(huán)與機(jī)身上的鎖鉤相連將座艙蓋固定到前機(jī)身上,經(jīng)過(guò)座艙蓋靜力疲勞試驗(yàn)和外場(chǎng)服役期間故障數(shù)據(jù)統(tǒng)計(jì),發(fā)現(xiàn)座艙蓋鎖環(huán)處側(cè)型材出現(xiàn)凹坑、裂紋(如圖1)。其主要原因?yàn)殒i鉤受載狀態(tài)受鎖鉤與鎖環(huán)的間隙影響,同時(shí)第4把鎖鉤結(jié)構(gòu)不同于其它鎖鉤,而左右兩側(cè)鎖環(huán)結(jié)構(gòu)形式均相同,靜力試驗(yàn)顯示第4把鎖環(huán)在飛行過(guò)程中所承受的載荷最大。

圖1 某型機(jī)座艙蓋鎖環(huán)處凹坑
基于ABAQUS有限元軟件建立圖2所示的有限元模型,分別建立以下3種表示不同工況的模型:修理前不含鋼板,修理后鋼板厚度分別為1 mm和1.5 mm。

圖2 有限元建立的側(cè)型材與鎖環(huán)結(jié)構(gòu)示意圖
其中側(cè)型材的材料為鋁合金,力學(xué)性能參數(shù)為:密度為2700 kg/m3,彈性模量為70 GPa,泊松比為0.33。鎖環(huán)、螺栓及鋼板的材料為鋼,力學(xué)性能參數(shù)為:密度為7800 kg/m3,彈性模量為202 GPa,泊松比為0.3。
網(wǎng)格劃分如圖3所示,對(duì)于鎖環(huán)采用C3D10網(wǎng)格類型,對(duì)于側(cè)型材、螺栓、鋼板采用C3D8網(wǎng)格類型。

圖3 鎖環(huán)、螺栓、鋼板和側(cè)型材的劃分網(wǎng)格示意圖
載荷條件和邊界條件如圖4所示,載荷施加在鎖環(huán)內(nèi)側(cè),沿離面方向(圖中-y方向)大小為20 549 N,沿指向艙內(nèi)方向(圖中-x方向)大小為3358 N,對(duì)兩側(cè)如圖所示位置進(jìn)行固支約束。

圖4 ABAQUS中施加載荷與邊界條件
考慮鎖環(huán)、鎖鉤由于接觸間隙不均勻而造成的載荷分散性,參考靜力試驗(yàn)中所測(cè)鎖環(huán)和鎖鉤上的試驗(yàn)數(shù)據(jù),將用于側(cè)型材與鎖環(huán)連接孔損傷容限分析的載荷譜按一定比例加重,設(shè)定初始裂紋設(shè)置為穿透裂紋3 mm。
修理前,整體結(jié)構(gòu)的Mises應(yīng)力云圖如圖5所示,最大應(yīng)力位于側(cè)型材孔邊。為方便按照?qǐng)D6所示的相鄰孔邊裂紋模型計(jì)算,即側(cè)型材與鎖環(huán)連接孔為有限大板相鄰孔邊裂紋,受釘傳載荷、遠(yuǎn)端(參考界面)彎矩和均勻拉應(yīng)力作用。結(jié)構(gòu)尺寸分別為:兩孔之間距離為26 mm;板厚t=3 mm;孔徑D=6.2 mm。選取S11方向觀察其載荷分布情況。

圖5 結(jié)構(gòu)Mises應(yīng)力云圖

圖6 相鄰孔邊裂紋
如圖7所示,在板的遠(yuǎn)端應(yīng)力較小,應(yīng)力主要集中于孔邊,并且在孔邊上下表面受力明顯不同,上表面最大應(yīng)力為拉應(yīng)力,下表面為壓應(yīng)力,上下表面間存在彎矩,故在計(jì)算時(shí)應(yīng)考慮彎矩的影響。

圖7 側(cè)型材及上下孔邊S11應(yīng)力云圖
結(jié)構(gòu)Mises應(yīng)力云圖如圖8所示。側(cè)型材及上下孔邊S11應(yīng)力云圖如圖9所示。

圖8 結(jié)構(gòu)Mises應(yīng)力云圖

圖9 側(cè)型材及上下孔邊S11應(yīng)力云圖
結(jié)構(gòu)Mises應(yīng)力云圖如圖10所示。側(cè)型材及上下孔邊S11應(yīng)力云圖如圖11所示。

圖10 結(jié)構(gòu)Mises應(yīng)力云圖

圖11 側(cè)型材及上下孔邊S11應(yīng)力云圖
如圖12所示,首先計(jì)算了Py和Pz載荷對(duì)螺栓處的彎矩(如圖12(a)),然后簡(jiǎn)化側(cè)型材與鎖環(huán)為一對(duì)邊簡(jiǎn)支矩形薄板在集中載荷作用下的彎矩模型(如圖12(b)),最后分別計(jì)算側(cè)型材兩個(gè)螺栓孔處的彎矩(如圖12(c))。
1)計(jì)算鎖環(huán)上Py和Pz載荷對(duì)螺栓的彎矩(如圖12(a))。Pz載荷與螺栓間的彎矩可簡(jiǎn)化為懸臂梁,此處Pz=3358 N,L=12.5 mm。Pz載荷與螺栓間的彎矩MPz=Pz×L=41975 N·mm;同理,可以計(jì)算Py載荷與螺栓間的彎矩,其中離心Py載荷與銷釘之間距離為1.5 mm。此處分別計(jì)算兩種不同Py值的情況:若Py=20549 N,彎矩M-Py=20549×1.5=30823.5 N·mm,其方向與M-Pz相反,故載荷對(duì)于螺栓連線的彎矩M螺栓=41975-30823.5=11151.5 N·mm;若Py=13000 N,彎矩M-Py=13000×1.5=19500 N·mm,M螺栓=41975-19500=22475 N·mm。

圖12 彎矩計(jì)算過(guò)程示意圖
2)計(jì)算鎖環(huán)當(dāng)量集中力,計(jì)算鎖環(huán)處彎矩,如圖12(b)所示,四邊簡(jiǎn)支板,與板的尺寸相比,兩螺栓孔間距可忽略不計(jì),認(rèn)為B點(diǎn)為載荷Py作用點(diǎn),若載荷Py=20549 N,根據(jù)彈性力學(xué)計(jì)算得到B點(diǎn)處的彎矩Mx=483114.68 N·mm;若載荷Py=13000 N,根據(jù)彈性力學(xué)計(jì)算得到B點(diǎn)處的彎矩Mx=305634.86 N·mm。
3)分別計(jì)算側(cè)型材兩個(gè)螺栓孔處的彎矩。如圖12(c)所示,將載荷等效為在兩個(gè)螺栓處,即將B點(diǎn)載荷等效到A和C。Py=20549 N時(shí),A點(diǎn)彎矩MA=483114.68×137÷150=441244.743 N·mm,C 點(diǎn)彎矩MC=483114.68×247÷260=458958.948 N·mm;Py=13000 N時(shí),A點(diǎn)彎矩MA=305634.86×137÷150=279146.50 N·mm,C點(diǎn)彎矩MC=305634.86×247÷260=290353.117 N·mm。
最后,分別計(jì)算兩個(gè)螺栓處的最終彎矩,即平板內(nèi)A、C各個(gè)點(diǎn)的彎矩加上螺栓處彎矩,即MA總=MA+M螺栓。Py=20549 N時(shí),A點(diǎn)彎矩MA總=452396.243 N·mm,C點(diǎn)彎矩MC總=483114.68×247÷260=470110.448 N·mm;Py=13000 N時(shí),A點(diǎn)彎矩MA總=301621.5 N·mm,C點(diǎn)彎矩MC總=312828.117 N·mm??梢钥闯觯槍?duì)不同的Py值,均顯示C點(diǎn)彎矩略大于A點(diǎn)。
此外,由于存在2個(gè)螺栓,故單個(gè)螺栓上的釘傳載荷P=Pz/2=1679 N。
裂紋類型為有限大板相鄰孔邊裂紋,受釘傳載荷、遠(yuǎn)端(參考界面)彎矩和均勻拉應(yīng)力作用。此處采用其1/2模型進(jìn)行計(jì)算分析,如圖13所示。
如圖13所示,尺寸分別為:W=260 mm(取鎖環(huán)至最后端面距離,圖12中C點(diǎn)距離);W=150 mm(圖12中A點(diǎn)距離);B=13 mm;D=6.2 mm。

圖13 計(jì)算模型示意圖(1/2模型):穿透裂紋
此處計(jì)算所采用材料為2系鋁合金包鋁板材,屈服強(qiáng)度σy=340 MPa,抗拉強(qiáng)度σb=452 MPa,其在高溫(150 ℃)下L-T 方向的裂紋擴(kuò)展速率曲線如圖14所示。

圖14 2系鋁合金包鋁板材L-T方向高溫da/dNΔK數(shù)據(jù)及Walker公式擬合線
采用Walker公式擬合:

相關(guān)參數(shù)具體擬合值為:C=2.209×10-8;n=3.60;M1=0.56(R≥0)。
結(jié)合圖13所示的計(jì)算模型,對(duì)于不同的工況,等效載荷計(jì)算公式為:

式中,寬度W按照完整模型寬度計(jì)算,W=410 mm。
有限元計(jì)算結(jié)果顯示在距離孔的遠(yuǎn)端,拉伸應(yīng)力非常小,近似等于0,這是由于邊界約束的存在,所以此處計(jì)算時(shí)取S0=0 MPa,對(duì)修理后的情形亦同。
1)臨界K準(zhǔn)則。Kmax超過(guò)斷裂韌度Kc,t=3 mm,Kc=67.64。
2)凈截面屈服準(zhǔn)則(NSY準(zhǔn)則)。
3)兩側(cè)裂紋尖端塑性區(qū)重疊連通,此處采用基于Irwin公式的循環(huán)塑性公式來(lái)計(jì)算:

4.5.1 修理前計(jì)算結(jié)果
此處分為4種情況,分別計(jì)算Py=20549 N、Py=13000 N時(shí)在W=260 mm(C點(diǎn))和W=150 mm(A點(diǎn))處對(duì)應(yīng)的裂紋擴(kuò)展壽命,然后取最危險(xiǎn)情況進(jìn)行后續(xù)分析。
1)第一種情況:Py=20549 N,W=260 mm。
S1=764.407 MPa,S3=90.268 MPa,失效方式為兩側(cè)裂尖塑性區(qū)重疊連通,臨界裂紋長(zhǎng)度為9.07 mm,最大應(yīng)力強(qiáng)度因子Kmax=52.07,對(duì)應(yīng)兩側(cè)裂紋擴(kuò)展剩余長(zhǎng)度為0.93 mm,計(jì)算臨界裂紋長(zhǎng)度循環(huán)塑性區(qū)長(zhǎng)度rp=1.86 mm,一半即為0.93 mm??梢员WC剩余裂紋可擴(kuò)展長(zhǎng)度即為兩側(cè)裂紋尖端塑性區(qū)重疊時(shí)的長(zhǎng)度。計(jì)算結(jié)果為:經(jīng)歷總循環(huán)數(shù)為857,等效飛行小時(shí)為1045.80 h。
2)第二種情況:Py=13000 N,W=260 mm。
S1=508.664 MPa,S3=90.268 MPa,失效方式為兩側(cè)裂尖塑性區(qū)重疊連通,臨界裂紋長(zhǎng)度為9.40 mm,最大應(yīng)力強(qiáng)度因子Kmax=41.77,對(duì)應(yīng)兩側(cè)裂紋擴(kuò)展剩余長(zhǎng)度為0.60 mm,計(jì)算臨界裂紋長(zhǎng)度循環(huán)塑性區(qū)長(zhǎng)度rp=1.20 mm,一半即為0.60 mm。可以保證剩余裂紋可擴(kuò)展長(zhǎng)度即為兩側(cè)裂紋尖端塑性區(qū)重疊時(shí)的長(zhǎng)度。計(jì)算結(jié)果:經(jīng)歷總循環(huán)數(shù)為3625;等效飛行小時(shí)為4462.52 h。
3)第三種情況:Py=20549 N,W=150 mm。
S1=735.604 MPa,S3=90.268 MPa;失效方式為兩側(cè)裂尖塑性區(qū)重疊連通,臨界裂紋長(zhǎng)度為9.105 mm,最大應(yīng)力強(qiáng)度因子Kmax=50.90,對(duì)應(yīng)兩側(cè)裂紋擴(kuò)展剩余長(zhǎng)度為0.895 mm,計(jì)算臨界裂紋長(zhǎng)度循環(huán)塑性區(qū)長(zhǎng)度rp=1.784 mm,一半即為0.892 mm。可以保證剩余裂紋可擴(kuò)展長(zhǎng)度即為兩側(cè)裂紋尖端塑性區(qū)重疊時(shí)的長(zhǎng)度。計(jì)算結(jié)果:經(jīng)歷總循環(huán)數(shù)為956,等效飛行小時(shí)為1170.70 h。
4)第四種情況:Py=13000 N,W=150 mm。
S1=490.441 MPa,S3=90.268 MPa,失效方式為兩側(cè)裂尖塑性區(qū)重疊連通,臨界裂紋長(zhǎng)度為9.42 mm,最大應(yīng)力強(qiáng)度因子Kmax=40.87,對(duì)應(yīng)兩側(cè)裂紋擴(kuò)展剩余長(zhǎng)度為0.58 mm,計(jì)算臨界裂紋長(zhǎng)度循環(huán)塑性區(qū)長(zhǎng)度rp=1.15 mm,一半即為0.58 mm。可以保證剩余裂紋可擴(kuò)展長(zhǎng)度即為兩側(cè)裂紋尖端塑性區(qū)重疊時(shí)的長(zhǎng)度。計(jì)算結(jié)果:經(jīng)歷總循環(huán)數(shù)為3950,等效飛行小時(shí)為4871.32 h。
可以看出在Py=20549 N、W=260 mm時(shí),裂紋擴(kuò)展壽命是最短的,故修理后均基于此固定參數(shù)進(jìn)行分析。
給出修理前Py=20549 N、W=260 mm時(shí)的裂紋長(zhǎng)度與飛行小時(shí)關(guān)系,如圖15所示。
剩余強(qiáng)度與裂紋長(zhǎng)度關(guān)系如圖16所示。
4.5.2 修理后加1 mm厚鋼板計(jì)算結(jié)果
失效方式為兩側(cè)裂尖塑性區(qū)重疊連通,臨界裂紋長(zhǎng)度為9.50 mm;最大應(yīng)力強(qiáng)度因子Kmax=38.47,對(duì)應(yīng)兩側(cè)裂紋擴(kuò)展剩余長(zhǎng)度為0.50 mm,計(jì)算臨界裂紋長(zhǎng)度循環(huán)塑性區(qū)長(zhǎng)度rp=1.01 mm,一半即為0.50 mm。可以保證剩余裂紋可擴(kuò)展長(zhǎng)度即為兩側(cè)裂紋尖端塑性區(qū)重疊時(shí)的長(zhǎng)度。經(jīng)歷總循環(huán)數(shù)為6838;等效飛行小時(shí)為8419.44 h。裂紋長(zhǎng)度與飛行小時(shí)關(guān)系如圖17所示。剩余強(qiáng)度與裂紋長(zhǎng)度關(guān)系如圖18所示。
4.5.3 修理后加1.5 mm厚鋼板
失效方式為兩側(cè)裂尖塑性區(qū)重疊連通;臨界裂紋長(zhǎng)為9.59 mm,最大應(yīng)力強(qiáng)度因子Kmax=34.08,對(duì)應(yīng)兩側(cè)裂紋擴(kuò)展剩余長(zhǎng)度為0.41 mm,計(jì)算臨界裂紋長(zhǎng)度循環(huán)塑性區(qū)長(zhǎng)度rp=0.80 mm,一半即為0.40 mm??梢员WC剩余裂紋可擴(kuò)展長(zhǎng)度即為兩側(cè)裂紋尖端塑性區(qū)重疊時(shí)的長(zhǎng)度。經(jīng)歷總循環(huán)數(shù)為15 720;等效飛行小時(shí)為19 356.20 h。裂紋長(zhǎng)度與飛行小時(shí)關(guān)系如圖19所示。剩余強(qiáng)度與裂紋長(zhǎng)度關(guān)系如圖20所示。
初始計(jì)算的裂紋擴(kuò)展壽命按照分散系數(shù)2進(jìn)行折算,計(jì)算結(jié)果如表1所示。

表1 側(cè)型材與鎖環(huán)連接孔裂紋擴(kuò)展壽命(飛行小時(shí)) h
座艙蓋骨架側(cè)型材修理后加1 mm和1.5 mm厚鋼板后,孔邊應(yīng)力明顯減小,采用1 mm鋼板加強(qiáng)后壽命比原始提高了約8倍,1.5 mm鋼板提高了約17倍,大大提高了座艙蓋危險(xiǎn)薄弱部位使用的可靠性,降低了座艙蓋外場(chǎng)使用的風(fēng)險(xiǎn),進(jìn)一步保障了飛行員的生命安全。

圖15 裂紋長(zhǎng)度-飛行小時(shí)曲線圖

圖16 剩余強(qiáng)度-裂紋長(zhǎng)度關(guān)系圖

圖17 裂紋長(zhǎng)度-飛行小時(shí)曲線圖

圖18 剩余強(qiáng)度-裂紋長(zhǎng)度關(guān)系圖

圖19 裂紋長(zhǎng)度-飛行小時(shí)曲線圖

圖20 剩余強(qiáng)度-裂紋長(zhǎng)度關(guān)系圖