付少青，劉凱，熊偉，周偉, 3

交變氣動(dòng)載荷作用下高速列車風(fēng)機(jī)底板連接鉚釘疲勞特性研究

付少青1, 2，劉凱1, 2，熊偉1, 2，周偉1, 2, 3

(1. 中南大學(xué) 交通運(yùn)輸工程學(xué)院，湖南 長沙 410075；2. 軌道交通安全教育部重點(diǎn)實(shí)驗(yàn)室，湖南 長沙 410075；3. 軌道交通安全關(guān)鍵技術(shù)國際合作聯(lián)合實(shí)驗(yàn)室，湖南 長沙 410075)

針對(duì)復(fù)雜服役環(huán)境交變氣動(dòng)載荷下高速列車底板-鉚釘結(jié)構(gòu)的疲勞損傷問題，提出基于服役氣動(dòng)載荷-鉚釘結(jié)構(gòu)應(yīng)力分布的疲勞載荷譜編制及損傷評(píng)估方法。首先，建立列車底板/鉚釘/梁板有限元模型，分析不同方向氣動(dòng)載荷下底板彎曲效應(yīng)對(duì)鉚釘結(jié)構(gòu)強(qiáng)度的影響機(jī)理，獲取氣動(dòng)載荷-鉚釘應(yīng)力的非線性傳遞因子；然后在列車底板布置氣動(dòng)差壓測(cè)點(diǎn)捕捉底板服役氣動(dòng)載荷譜；最后，結(jié)合服役氣動(dòng)載荷譜和非線性傳遞因子編制鉚釘結(jié)構(gòu)的動(dòng)應(yīng)力載荷譜，基于雨流計(jì)數(shù)和Miners累積損傷準(zhǔn)則建立鉚釘結(jié)構(gòu)二維八級(jí)載荷譜，并采用第四強(qiáng)度理論應(yīng)力準(zhǔn)則進(jìn)行疲勞損傷評(píng)估。研究結(jié)果表明：在交變載荷下，高速動(dòng)車組達(dá)到服役壽命960萬km時(shí)，鉚釘結(jié)構(gòu)2個(gè)危險(xiǎn)點(diǎn)在3個(gè)分區(qū)氣動(dòng)載荷及平均氣動(dòng)載荷作用下的疲勞總損傷分別為0.64，0.53，0.56，0.51和0.19，0.52，0.36，0.13。疲勞總損傷達(dá)到1.0，即認(rèn)為鉚釘結(jié)構(gòu)氣動(dòng)致疲勞失效，由計(jì)算結(jié)果可知，高速動(dòng)車組服役壽命內(nèi)，鉚釘結(jié)構(gòu)能夠滿足使用要求。

高速列車；氣動(dòng)載荷；鉚釘結(jié)構(gòu)；載荷譜；疲勞損傷評(píng)估

在工程連接件應(yīng)用中，鉚接連接件以其良好的機(jī)械性能、連接可靠、質(zhì)量輕、低成本、較強(qiáng)的環(huán)境適用性等特點(diǎn)而廣泛應(yīng)用于航空航天、機(jī)械、汽車及船舶等領(lǐng)域[1]。在高鐵領(lǐng)域，我國高速列車的設(shè)備艙底板也廣泛采用了鉚接。近年來隨著列車運(yùn)行速度的不斷提高，顯著增加的空氣動(dòng)力學(xué)問題，列車所面臨的力學(xué)環(huán)境變得更加復(fù)雜[2]，因此高速動(dòng)車組的服役安全性能也逐漸受到國內(nèi)外學(xué)者的研究和關(guān)注。設(shè)備艙底板鉚釘結(jié)構(gòu)由于長期處在交變載荷的惡劣環(huán)境下工作，特別是開行在武廣、鄭西等存在大量隧道和隧道群的高速動(dòng)車組[3]，列車風(fēng)機(jī)底板要受到較劇烈的氣動(dòng)載荷和沖擊振動(dòng)載荷的共同作用，而設(shè)備艙風(fēng)機(jī)安裝底板的高應(yīng)力主要是由氣動(dòng)載荷引起的[4]，因此研究氣動(dòng)載荷對(duì)鉚釘疲勞的影響尤為必要。本文為研究交變氣動(dòng)載荷對(duì)鉚疲勞脫落的影響，從載荷角度分析高速動(dòng)車組在設(shè)計(jì)服役壽命960萬km內(nèi)，鉚釘是否可以正常使用，而對(duì)鉚釘疲勞損傷進(jìn)行了預(yù)測(cè)。針對(duì)鉚釘連接疲勞失效現(xiàn)象，國內(nèi)外學(xué)者從不同角度出發(fā)，對(duì)鉚接的疲勞壽命研究取得了一定成果。然而由于影響結(jié)構(gòu)疲勞的因素過多、疲勞破壞沒有先期預(yù)兆，研究困難較大，主要通過試驗(yàn)和仿真來分析鉚釘產(chǎn)生疲勞的主要因素及對(duì)比其他連接方式來研究鉚釘鉚接的疲勞性能[5]。李風(fēng)等[6]試驗(yàn)研究了在不同環(huán)境下的鉚接疲勞強(qiáng)度曲線，并對(duì)鉚釘搭接頭疲勞壽命表達(dá)式進(jìn)行了推導(dǎo)。Urban[7]對(duì)機(jī)身鉚接結(jié)構(gòu)進(jìn)行了疲勞試驗(yàn)。Blanchot等[8]對(duì)鉚釘進(jìn)行有限元模擬，分析了鉚釘結(jié)構(gòu)應(yīng)力分布等。Szymczyk等[9]對(duì)鉚釘接頭的位移和應(yīng)力進(jìn)行了分析，提出了預(yù)測(cè)鉚釘連接件疲勞壽命的方法。目前針對(duì)鉚釘疲勞的研究多局限于靜強(qiáng)度連接結(jié)構(gòu)的疲勞分析，對(duì)交變氣動(dòng)載荷因素對(duì)鉚釘結(jié)構(gòu)疲勞壽命的影響尚未有全面的研究，也缺乏全面的數(shù)據(jù)，本文主要考慮交變氣動(dòng)載荷作用下列車底板鉚釘結(jié)構(gòu)的疲勞問題。論文通過仿真方法得到氣動(dòng)載荷?鉚釘應(yīng)力的非線性傳遞因子，然后基于高速列車服役運(yùn)行條件下的真實(shí)交變氣動(dòng)載荷譜作為依據(jù)，并根據(jù)服役載荷譜和非線性傳遞因子編制鉚釘結(jié)構(gòu)的動(dòng)應(yīng)力載荷譜，基于雨流計(jì)數(shù)和Miners累積損傷準(zhǔn)則建立鉚釘結(jié)構(gòu)的二維八級(jí)載荷譜，最后采用第四強(qiáng)度理論應(yīng)力(VON MISES)準(zhǔn)則對(duì)列車鉚釘?shù)钠趽p傷及壽命進(jìn)行評(píng)估，為后續(xù)連接件更換和維護(hù)提供理論 依據(jù)。

1 仿真計(jì)算模型

1.1 幾何模型及材料參數(shù)

模型由底板，梁板以及用于連接底板和梁板的鉚釘組成。其中梁板(底架橫梁)是車體的一部分，底板通過鉚釘鉚接在梁板上。模型實(shí)際位于CRH380高速列車設(shè)備艙主風(fēng)機(jī)安裝底板上，底板作為設(shè)備艙的一部分，承受著來自設(shè)備艙上部的交變氣動(dòng)載荷作用。該模型底板長1 018 mm，寬363 mm，厚1.5 mm，梁板厚度5 mm，沿底板長度方向兩側(cè)各3個(gè)鉚釘孔，用于安裝鉚釘，孔徑5.5 mm，鉚釘直徑4.8 mm，模型尺寸位置如圖2所示。同時(shí)，在不影響計(jì)算結(jié)果的前提下，為了便于網(wǎng)格劃分，本模型刪除底板上的細(xì)小圓孔及倒角。模型幾何尺寸和實(shí)物圖如圖1(a)～1(b)所示。

安裝底板，梁板以及鉚釘，其各部分的材料參數(shù)如表1所示。

單位：mm

表1 模型各部分材料參數(shù)

1.2 網(wǎng)格劃分

本模型利用ABAQUS軟件，采用三維八節(jié)點(diǎn)C3D8R實(shí)體單元并采用縮減積分計(jì)算[10]，這種類型的單元可以使計(jì)算結(jié)果更容易收斂。為保證結(jié)果的準(zhǔn)確性，整體結(jié)構(gòu)采用六面體進(jìn)行網(wǎng)格劃分，網(wǎng)格數(shù)量對(duì)鉚釘連接處應(yīng)力有影響，因此對(duì)重點(diǎn)關(guān)注的鉚釘與底板和梁板的接觸區(qū)域進(jìn)行網(wǎng)格加密處理，因?yàn)榇颂幙赡軙?huì)產(chǎn)生較大變形或應(yīng)力集中，其他區(qū)域網(wǎng)格較為稀疏，從而建立起結(jié)構(gòu)有限元模型。有限元模型共產(chǎn)生C3D8R實(shí)體單元86 580個(gè), 節(jié)點(diǎn)112 103個(gè)，網(wǎng)格數(shù)量滿足工程要求。底板?鉚釘?梁板有限元模型如圖2所示。

1.3 載荷、約束與接觸

實(shí)驗(yàn)所得鉚釘?shù)睦鞓O限強(qiáng)度為1 kN, 即釘芯的拉斷力為1 kN，依據(jù)力的相互作用原理，得出拉鉚釘?shù)念A(yù)緊力：r=1 kN。氣動(dòng)載荷由武廣線上CRH380AL列車服役條件下的交變氣動(dòng)載荷獲得，在底板沿車長方向的前、中及尾部位置布置差壓傳感器測(cè)點(diǎn)3個(gè)，如圖3所示。

3個(gè)測(cè)點(diǎn)的氣動(dòng)載荷時(shí)程曲線如圖4所示。

圖2 底板?鉚釘?梁板有限元模型

圖3 風(fēng)機(jī)底板結(jié)構(gòu)氣動(dòng)載荷測(cè)點(diǎn)布置(沿車長方向)

根據(jù)圖4時(shí)程曲線可知，氣動(dòng)載荷在?450 Pa至450 Pa之間，因此在底板整個(gè)面上施加從上往下的梯度均布載荷壓力，大小450，350，…，50和從下往上的均布載荷?450，?350，…，?50，單位均為Pa。

約束：梁板作為車體的一部分，其在，，3個(gè)方向的位移是受到約束的，故約束平移自由度PINNED：1=2=3。

圖4 測(cè)點(diǎn)氣動(dòng)載荷時(shí)程曲線

接觸問題在工程中是一種高度的非線性行 為[11]，該實(shí)體模型，考慮到了接觸。在abquas 中鉚釘分別與梁板和底板的接觸1和3，梁板與底板之間接觸2，底板與鉚釘體接觸4以及鉚釘體與梁板之間的接觸5均設(shè)置了接觸，摩擦因數(shù)0.17，該模型共定義接觸30對(duì)。接觸如圖5示。

圖5 接觸設(shè)置位置示意圖

2 氣動(dòng)載荷?結(jié)構(gòu)應(yīng)力解算

2.1 結(jié)構(gòu)應(yīng)力分布

當(dāng)在底板施加不同方向的氣動(dòng)載荷時(shí)，MISES最大應(yīng)力為A與B 2點(diǎn)。下面就在不同方向施加氣動(dòng)載荷時(shí)鉚釘MISES最大值點(diǎn)應(yīng)力變化規(guī)律進(jìn)行分析。

在1 000 N預(yù)緊力下，自上而下分別施加450 Pa和50 Pa氣動(dòng)壓力時(shí)A點(diǎn)應(yīng)力云圖如圖6。

在1 000 N預(yù)緊力下，自下而上分別施加?450 Pa和?50 Pa氣動(dòng)壓力時(shí)B點(diǎn)應(yīng)力云圖如圖7所示。

(a) 450 Pa鉚釘應(yīng)力云圖；(b) 50 Pa鉚釘應(yīng)力云圖

(a) 450 Pa鉚釘應(yīng)力云圖；(b) 50 Pa鉚釘應(yīng)力云圖

在1 000 N預(yù)緊力下，無氣動(dòng)載荷施加時(shí)應(yīng)力云圖如圖8所示。

圖8 氣動(dòng)載荷為0，預(yù)緊力1 000 N時(shí)應(yīng)力云圖

通過靜力學(xué)分析得出其最大應(yīng)力分布理論上存在于鉚釘和被連接件的接觸部位以及桿部分。從仿真結(jié)果來看，鉚釘最大應(yīng)力出現(xiàn)在了鉚釘頭和被連接件的接觸部位，與理論分析相一致。從應(yīng)力云圖可以看出，氣動(dòng)載荷變化對(duì)鉚釘應(yīng)力的影響較預(yù)緊力相對(duì)較小。

2.2 應(yīng)力傳遞因子計(jì)算

模型計(jì)算了鉚釘在預(yù)緊力下，對(duì)底板施加從上向下及從下往上的梯度氣動(dòng)壓力值，得到鉚釘結(jié)構(gòu)的應(yīng)力大小及分布規(guī)律，并得到MISES應(yīng)力各分量與氣動(dòng)載荷之間的數(shù)學(xué)關(guān)系式。

根據(jù)仿真得到A點(diǎn)MISES應(yīng)力分量與氣動(dòng)載荷間的關(guān)系曲線如圖9(a)～9(f)所示。

(a) Sx與氣動(dòng)載荷曲線擬合式；(b) Sy與氣動(dòng)載荷曲線擬合式；(c) Sz與氣動(dòng)載荷曲線擬合式；(d) Sxy與氣動(dòng)載荷曲線擬合式；(e) Syz與氣動(dòng)載荷曲線擬合式；(f)Sxz與氣動(dòng)載荷曲線擬合式

根據(jù)結(jié)構(gòu)仿真結(jié)果繪制B點(diǎn)六應(yīng)力分量與氣動(dòng)載荷的關(guān)系曲線如圖10(a)～10(f)所示。

因此，根據(jù)圖9和圖10可以得到釘結(jié)構(gòu)MISES各應(yīng)力分量與氣動(dòng)載荷之間的一一映射關(guān)系?即可得應(yīng)力傳遞因子。

(a) Sx與氣動(dòng)載荷曲線擬合式；(b) Sy與氣動(dòng)載荷曲線擬合式；(c) Sz與氣動(dòng)載荷曲線擬合式；(d) Sxy與氣動(dòng)載荷曲線擬合式；(e) Syz與氣動(dòng)載荷曲線擬合式；(f) Sxz與氣動(dòng)載荷曲線擬合式

3 鉚釘結(jié)構(gòu)疲勞分析及壽命預(yù)測(cè)

3.1 鉚釘疲勞載荷譜的編制

3.1.1 風(fēng)載時(shí)程?應(yīng)力時(shí)程轉(zhuǎn)換

要對(duì)鉚釘結(jié)構(gòu)進(jìn)行疲勞評(píng)估，應(yīng)首先將得到的氣動(dòng)載荷譜轉(zhuǎn)換為應(yīng)力載荷譜。本文采用第4強(qiáng)度理論應(yīng)力大小EQV進(jìn)行分析。第4強(qiáng)度理論應(yīng)力大小EQV公式如式(1)所示

式中：S，S，S，S，S和S為仿真六應(yīng)力分量。

根據(jù)仿真得到的六應(yīng)力分量，代入式(1)便可得到在氣動(dòng)載荷下不同時(shí)刻EQV的大小。鉚釘A點(diǎn)和B點(diǎn)第4強(qiáng)度理論應(yīng)力分布如圖11(a)和11(b)所示(以3測(cè)點(diǎn)平均載荷為例)。

3.1.2 八級(jí)載荷譜編制

本文采用雨流計(jì)數(shù)法對(duì)載荷均幅值進(jìn)行統(tǒng)計(jì)后，得到應(yīng)力時(shí)間歷程中各次載荷循環(huán)中應(yīng)力均幅值。隨后通過編制載荷譜對(duì)載荷均幅值進(jìn)行簡化，從而得到能夠表征鉚釘結(jié)構(gòu)真實(shí)氣動(dòng)載荷作用下的應(yīng)力狀況以及可用于疲勞損傷計(jì)算的程序譜[12]。編制的疲勞載荷譜應(yīng)能真實(shí)模擬列車底板在實(shí)際運(yùn)行過程中所受的交變載荷和發(fā)現(xiàn)疲勞薄弱位置并能夠準(zhǔn)確預(yù)測(cè)鉚釘結(jié)構(gòu)的疲勞壽命[13]，而八級(jí)載荷譜可以滿足工程精度要求，故采用二維八級(jí)載荷譜計(jì)算。

鉚釘危險(xiǎn)點(diǎn)A在3個(gè)分區(qū)及其平均二維應(yīng)力均幅值載荷譜三維柱狀云圖如圖12所示。

(a) A點(diǎn)第4強(qiáng)度理論應(yīng)力分布；(b) B點(diǎn)第4強(qiáng)度理論應(yīng)力分布

(a) 3點(diǎn)平均值；(b) 測(cè)點(diǎn)1；(c) 測(cè)點(diǎn)2；(d) 測(cè)點(diǎn)3

危險(xiǎn)點(diǎn)A在底板分區(qū)3個(gè)測(cè)點(diǎn)應(yīng)力均值分布有所不同，測(cè)點(diǎn)1、測(cè)點(diǎn)2和測(cè)點(diǎn)3應(yīng)力均值分別主要分布在123～141 MPa，99～131 MPa和105～140 MPa。3個(gè)測(cè)點(diǎn)波動(dòng)中心也各有差異，分別依次在114～176 MPa之間、67～179 MPa之間和93～176 MPa之間。

鉚釘危險(xiǎn)點(diǎn)B在3個(gè)分區(qū)及其平均二維應(yīng)力均幅值載荷譜三維柱狀云圖如圖13所示。

(a) 3點(diǎn)平均值；(b) 測(cè)點(diǎn)1；(c) 測(cè)點(diǎn)2；(d) 測(cè)點(diǎn)3

危險(xiǎn)點(diǎn)B在底板分區(qū)3個(gè)測(cè)點(diǎn)應(yīng)力均值分布有所不同，測(cè)點(diǎn)1，測(cè)點(diǎn)2和測(cè)點(diǎn)3應(yīng)力均值分別主要分布在104～123 MPa，113～146 MPa和115～140 MPa。3個(gè)測(cè)點(diǎn)波動(dòng)中心也不相同，分別依次在66～133 MPa之間、64～179 MPa之間和67～152 MPa之間。

3.2 鉚釘疲勞壽命的評(píng)估

3.2.1 材料疲勞性能和疲勞參數(shù)

1) 參考相關(guān)文獻(xiàn)[14]得到A6N01S鋁合金材料的曲線如圖14所示。

2) 鋁合金材料疲勞參數(shù)，定義如表2所示。

圖14 A6N01S鋁合金材料的S-N曲線

表2 A6N01S鋁合金材料參數(shù)[14?15]

3.2.2 二維損傷評(píng)估

材料的曲線方程可表示為：

式(2)中：sn，為材料常數(shù)。

在第i級(jí)載荷幅ai作用下，載荷循環(huán)次數(shù)達(dá)i次時(shí)產(chǎn)生的疲勞損傷為：

將式(3)代入式(2)，于是單個(gè)幅值下的疲勞 損傷：

式中：N為該級(jí)載荷幅對(duì)應(yīng)的疲勞極限循環(huán)數(shù)；S為各級(jí)載荷水平的幅值；C為各級(jí)載荷幅對(duì)應(yīng)的循環(huán)次數(shù)。

二維載荷譜綜合考慮了載荷均幅值的影響，根據(jù)Goodman等壽命轉(zhuǎn)換方程將各級(jí)載荷轉(zhuǎn)換為對(duì)稱循環(huán)載荷，通過曲線計(jì)算疲勞損傷。根據(jù)Goodman等壽命轉(zhuǎn)換方程可得[16]：

式中：?1為對(duì)稱循環(huán)載荷幅值；a為載荷幅值；m為均值；u為材料極限強(qiáng)度。

由曲線可計(jì)算對(duì)應(yīng)的載荷循環(huán)次數(shù)N，對(duì)應(yīng)疲勞損傷D，則疲勞總損傷為

3.2.3 鉚釘結(jié)構(gòu)疲勞損傷和疲勞壽命計(jì)算結(jié)果

采用上述3.2.2節(jié)中二維載荷均幅值譜的疲勞評(píng)估方法，得到鉚釘危險(xiǎn)點(diǎn)A，B處底板分區(qū)3個(gè)測(cè)點(diǎn)單日疲勞損傷和服役960萬km時(shí)總疲勞損傷分別如表3～4所示。

表3 鉚釘A，B處3個(gè)測(cè)點(diǎn)疲勞損傷計(jì)算結(jié)果(單日運(yùn)行里程：2 100 km)

基于二維疲勞應(yīng)力譜的疲勞損傷計(jì)算結(jié)果可知：單日運(yùn)行里程2 100 km工況下，鉚釘結(jié)構(gòu)危險(xiǎn)點(diǎn)A 3個(gè)分區(qū)氣動(dòng)載荷及其平均氣動(dòng)載荷疲勞損傷分別為：01號(hào)1.41×10?4，02號(hào)1.16×10?4，03號(hào)1.22×10?4，平均1.11×10?4，對(duì)應(yīng)的達(dá)到動(dòng)車組設(shè)計(jì)服役壽命960萬km時(shí)鉚釘結(jié)構(gòu)疲勞總損傷分別為0.64，0.53，0.56，0.51；鉚釘結(jié)構(gòu)危險(xiǎn)點(diǎn)B 3個(gè)分區(qū)氣動(dòng)載荷及其平均氣動(dòng)載荷疲勞損傷分別為：01號(hào)為4.12×10?5，02號(hào)為1.15×10?4，03號(hào)為7.85×10?5，平均2.71×10?5，對(duì)應(yīng)的達(dá)到設(shè)計(jì)服役壽命960萬km時(shí)鉚釘結(jié)構(gòu)疲勞總損傷分別為0.19，0.52，0.36，0.13。疲勞總損傷達(dá)到1.0，即可以認(rèn)為鉚釘發(fā)生了疲勞破壞，從計(jì)算結(jié)果可以看出，可以認(rèn)為在高速動(dòng)車組服役960萬km里程壽命內(nèi)，鉚釘能夠滿足使用要求。

表4 對(duì)應(yīng)的服役960萬公里鉚釘結(jié)構(gòu)疲勞總損傷計(jì)算結(jié)果

4 結(jié)論

1) 有限元仿真是分析高速列車底板?鉚釘結(jié)構(gòu)疲勞損傷與壽命的重要方法。通過對(duì)底板?鉚釘結(jié)構(gòu)不同方向施加梯度壓力加載仿真，得到鉚釘結(jié)構(gòu)MISES應(yīng)力六分量與氣動(dòng)載荷之間數(shù)學(xué)關(guān)系式，為氣動(dòng)載荷譜→應(yīng)力載荷譜的關(guān)系轉(zhuǎn)換做好準(zhǔn)備。

2) 對(duì)于非對(duì)稱載荷時(shí)間歷程，載荷均值變化較大，由于二維載荷譜同時(shí)考慮載荷幅值和均值2個(gè)參數(shù)的影響，故采用二維八級(jí)載荷譜對(duì)鉚釘結(jié)構(gòu)進(jìn)行疲勞壽命估算可以準(zhǔn)確反映鉚釘?底板結(jié)構(gòu)的實(shí)際受載情況，具備較高的工程參考價(jià)值。

3) 由以上計(jì)算可知，單日運(yùn)行2 100 kM工況下，鉚釘結(jié)構(gòu)A和B危險(xiǎn)點(diǎn)3個(gè)分區(qū)測(cè)點(diǎn)氣動(dòng)載荷及其平均氣動(dòng)載荷疲勞損傷分別為：01號(hào)1.41×10?4，02號(hào)1.16×10?4，03號(hào)1.22×10?4，平均1.11×10?4和01號(hào)為4.12×10?5，02號(hào)為1.15×10?4，03號(hào)為7.85×10?5、平均2.71×10?5，對(duì)應(yīng)的A和B危險(xiǎn)點(diǎn)達(dá)到高速動(dòng)車組設(shè)計(jì)服役壽命960萬公里里程時(shí)鉚釘結(jié)構(gòu)疲勞總損傷分別為0.64，0.53，0.56，0.51和 0.19，0.52，0.36，0.13。疲勞總損傷達(dá)到1.0，即認(rèn)為鉚釘結(jié)構(gòu)氣動(dòng)致疲勞失效，由計(jì)算結(jié)果可知，高速動(dòng)車組服役壽命內(nèi)，鉚釘結(jié)構(gòu)能夠滿足使用要求。

[1] 陳福玉. 航空鉚接連接件疲勞壽命研究[D]. 南京: 南京航空航天大學(xué). 2011. CHEN Fuyu. Study on fatigue life of aeronautical riveted joints[D]. Nanjing: Nanjing University of Aeronautics & Astronzutics, 2011.

[2] 陳曉麗. 高速列車出隧道空氣動(dòng)力效應(yīng)數(shù)值模擬研究[D]. 蘭州: 蘭州交通大學(xué), 2016. CHEN Xiaoli. Numerical simulation of aerodynamic effects of high-speed train exiting tunnel[D]. Lanzhou: Lanzhou Jiaotong University, 2016.

[3] 罔旭東. 高速列車車下設(shè)備艙底板的可靠性研究[D].北京: 北京交通大學(xué), 2014. WANG Xudong. Reliability study of high speed train equipment bay bottom plates[D]. Beijing: Beijing Jiaotong University, 2014.

[4] 余思均. 高速列車氣動(dòng)載荷車體疲勞強(qiáng)度研究[D]. 成都: 西南交通大學(xué), 2011. YU Sijun. Study on the fatigue strength of high speed train aerodynamic load car body[D]. Chengdu: Southwest Jiaotong University, 2011.

[5] 許樂樂. 基于疲勞的鋁合金板件無鉚釘鉚接接頭布置規(guī)律研究[D]. 長春: 吉林大學(xué), 2016.XU Lele. Research on arrangement regularity of rivet less riveted joints of aluminum alloy plates based on fatigue[D]. Changchun: Jilin University, 2016.

[6] 李風(fēng), 陳海燕, 張毅. 鉚釘搭接頭疲勞強(qiáng)度曲線的研究[J]. 機(jī)械強(qiáng)度, 2002, 24(3): 453?454. LI Feng, CHEN Haiyan, ZHANG Yi. Study on the fatigue strength curves for riveted lap joints[J]. Journal of Mechanical Strength, 2002, 24(3): 453?454.

[7] Urban M R. Analysis of the fatigue life of riveted sheet metal helicopter airframe joints[J]. International Journal of Fatigue, 2003(25): 1013?1026.

[8] Blanchot V, Alain Daidie A. Riveted assembly modeling: Study and numerical characterization of a riveting process[J]. Materials Processing Technology, 2006(180): 201?209.

[9] Szymczyk E, Derewoko A, Jachimowicz J. Analysis of displacement and stress distributions in riveted joints[C]// III European Conference on Computational Mechanics Solids, Structures and Coupled Problems in Engineering, 2006: 434.

[10] 程祖鋒, 師歡歡, 徐光兵. 基于ABAQUS的深基坑變形和內(nèi)力三維有限元分析[J]. 煤炭工程, 2013, 45(9): 119?121. CHENG Zufeng, SHI Huanhuan, XU Guangbing. 3D finite element analysis of deformation and internal force of deep foundation pit based on ABAQUS[J]. Coal Engineering, 2013, 45(9): 119?121.

[11] 李妍. 基于ANSYS軟件的接觸問題分析及在工程中的應(yīng)用[D]. 長春: 吉林大學(xué), 2004. LI Yan. The contact element analysis and application in engineering based on the ANSYS software[D]. Changchun: Jilin University, 2004.

[12] 毛賀. 高速列車載荷譜編制方法的研究[D]. 北京: 北京交通大學(xué), 2009. MAO He. Research on the method of compiling load spectrum of high-speed train[D]. Beijing: Beijing Jiaotong University, 2009.

[13] 宋燁. 動(dòng)車組鋁合金車體疲勞壽命評(píng)估理論與試驗(yàn)研究[D]. 成都: 西南交通大學(xué), 2016. SONG Yi. Evaluation theory and experimental study of emu carbody fatigue life[D]. Chengdu: Southwest Jiaotong University, 2016.

[14] 陳書翔, 王飛, 寇福俊, 等. 服役后的A6N01S鋁合金型材焊接接頭組織與性能研究[J]. 輕合金加工技術(shù), 2015, 43(8): 64?68. CHEN Shuxiang, WANG Fei, KOU Fujun, et al. Study on microstructure and properties of welded joints of a6n01s aluminum alloy profiles after service[J]. Light Alloy Processing Technology, 2015, 43 (8): 64?68.

[15] 陳輝剛. 高速列車用A6N01S-T5鋁合金及其焊接接頭高寒條件下的性能[J]. 電焊機(jī), 2016, 46(9): 77?82 CHEN Huigang. Performance of the A6N01S-T5 aluminum alloy and the welding joint for high speed train at low temperature[J]. Electric Welder, 2016, 46(9): 77? 82.

[16] 宋振偉, 曹俊偉, 劉玲, 等. 船用鈦、鋼法蘭連接螺栓的疲勞強(qiáng)度性能分析[J]. 中國艦船研究, 2018, 13(增1): 106?110. SONG Zhenwei, CAO Junwei, LIU Ling, et al. Analysis of fatigue strength of titanium-steel flange connection bolts in ship[J]. China Ship Research, 2018, 13(Suppl 1): 106?110.

Research on fatigue characteristics of high speed train fan bottom plate rivets under alternating aerodynamic load

FU Shaoqing1, 2, LIU Kai1, 2, XIONG Wei1, 2, ZHOU Wei1, 2, 3

(1. School of Traffic and Transportation Engineering, Central South University, Changsha 410075, China; 2. Key Laboratory of Traffic Safety on Track of Ministry of Education, Changsha 410075, China; 3. Joint International Research Laboratory of Key Technology for Rail Traffic Safety, Changsha 410075, China)

Aiming at the fatigue damage of high-speed train floor-rivet structure under alternating aerodynamic load in a complex service environment, a fatigue load spectrum compilation and damage assessment method based on the service aerodynamic load-rivet structure stress distribution was proposed. First, a finite element model of the train floor/rivets/beam-slab was established, and the mechanism of the influence of the bending effect of the bottom plate on the strength of the rivet structure under aerodynamic loads in different directions was analyzed, and the nonlinear transfer factor of the aerodynamic load-rivet stress was obtained; The pressure measurement points capture the aerodynamic load spectrum of the bottom plate service. Finally, the dynamic stress load spectrum of the rivet structure was compiled based on the service aerodynamic load spectrum andthen on-linear transfer factor, and the two-dimensional eight-level load spectrum of the rivet structure was established based on the rain flow count and Miners cumulative damage criterion. The fourth strength theory stress criterion was used for fatigue damage assessment. The conclusions show that when the high-speed EMU reaches the service life of 9.6 million km under alternating loads, the total fatigue damage of the two dangerous points of the rivet structure under the effects of aerodynamic loads and average aerodynamic loads in the three zones are 0.64, 0.53, 0.56, 0.51 and 0.19, 0.52, 0.36, 0.13. The total fatigue damage reaches 1.0, which means the fatigue failure of the rivet structure under aerodynamic load is considered. From the calculation results, it can be known that the rivet structure can meet the service requirements during the service life of the high-speed EMU.

high-speed train; aerodynamic load; rivet structure; load spectrum; fatigue damage assessment

U266

A

1672 ? 7029(2020)12 ? 3003 ? 10

10.19713/j.cnki.43?1423/u.T20200108

2020?02?14

國家重點(diǎn)研發(fā)計(jì)劃資助項(xiàng)目(2017YFB1201201)

周偉(1982?)，男，湖南長沙人，副教授，博士，從事軌道交通裝備疲勞載荷譜在途監(jiān)測(cè)與壽命評(píng)估；E?mail：zhou_wei000@126.com

(編輯 涂鵬)