復(fù)雜交變載荷條件下微動運行狀態(tài)

田帥，李欣，楊建偉，劉樂強(qiáng)，張育軒

復(fù)雜交變載荷條件下微動運行狀態(tài)

田帥，李欣，楊建偉，劉樂強(qiáng)，張育軒

（北京建筑大學(xué) 機(jī)電與車輛工程學(xué)院，北京 100044）

研究典型的加載參數(shù)對復(fù)雜交變載荷作用下微動運行狀態(tài)的影響。利用ABAQUS建立了二維（平面應(yīng)變）有限元模型，模擬微動墊和試件在循環(huán)法向載荷和軸向載荷下的接觸狀態(tài)，對不同加載條件下的微動運行狀態(tài)進(jìn)行分析，提出一種曲線的分析方法，再結(jié)合雙軸微動疲勞試驗驗證其適用性。比例加載條件下，在微動墊夾具剛度較小時的t的曲線為直線，而微動墊夾具剛度較大時的t曲線表現(xiàn)為不規(guī)則的四邊形，這些條件下的曲線均為線性函數(shù)，并且隨著微動墊夾具剛度增加，曲線的斜率增大；非比例加載條件下，t-曲線的形狀不再是四邊形，形狀較為復(fù)雜，此時曲線為橢圓函數(shù)。隨著微動墊夾具剛度的增大，橢圓的短軸增加。當(dāng)相位差不為90°時，橢圓兩半軸與坐標(biāo)系不平行，橢圓發(fā)生旋轉(zhuǎn)；特別地，相位差為90°時，曲線為標(biāo)準(zhǔn)橢圓函數(shù)。當(dāng)曲線與直線=±相交時，由于滑移的產(chǎn)生，橢圓曲線將發(fā)生變化。在復(fù)雜交變載荷作用下，不同的加載參數(shù)（法向載荷、軸向載荷、微動墊夾具剛度等）或其相互間的組合會影響曲線的大小和位置，提出的曲線方法可以為分析復(fù)雜交變載荷條件下的運行狀態(tài)研究提供手段，為進(jìn)一步地討論微動疲勞或微動磨損行為提供指導(dǎo)。

微動；多軸交變載荷；動力學(xué)；接觸狀態(tài)；有限元仿真

微動是指在交變載荷的作用下，兩個相互接觸的表面之間發(fā)生微小滑動而引起的磨損和疲勞現(xiàn)象[1]。微動常見于各種緊固連接的結(jié)合面上，由微動磨損導(dǎo)致接觸表面的損傷，在交變載荷作用下，結(jié)構(gòu)更容易發(fā)生疲勞破壞。相同試件的疲勞試驗，存在微動作用時的疲勞壽命比不存在微動作用時的疲勞壽命降低30%～50%[2-4]。

航空發(fā)動機(jī)葉片和輪盤一般采用榫連接結(jié)構(gòu)，這是一種接觸結(jié)構(gòu)的連接形式[5-6]。在發(fā)動機(jī)工作過程中，葉片受到離心力、氣動力、振動等復(fù)雜的載荷作用，使葉片根部的榫頭與輪盤的榫槽之間產(chǎn)生微小的運動，造成兩者接觸界面上的微動疲勞[7-9]。榫連接結(jié)構(gòu)是航空發(fā)動機(jī)故障的多發(fā)部位之一，微動疲勞會加速榫連接結(jié)構(gòu)接觸部位裂紋的萌生和擴(kuò)展，大大降低零件的疲勞壽命[10-11]。顯然，在飛行過程中，輪盤和葉片接觸界面之間的法向載荷不斷變化，這意味著微動通常發(fā)生在多軸可變載荷條件下。一些研究者注意到了這一局限性，并進(jìn)行了相關(guān)研究。

米雪等[12]通過690合金管/405不銹鋼塊（線接觸）的切向微動試驗，探究了690合金管在不同法向載荷作用下的切向微動磨損機(jī)制和損傷演變規(guī)律。景鵬飛等[13]在不同的載荷和位移幅值下，結(jié)合微動圖以及不同微動接觸狀態(tài)和滑移狀態(tài)下材料的損傷機(jī)理，得出切向摩擦力、微動振幅是影響微動磨損的重要因素。石煒等[14]在自行設(shè)計的微動疲勞試驗裝置上進(jìn)行了各向異性材料的微動疲勞試驗，分析了軸向和法向應(yīng)力對微動疲勞壽命以及微動接觸區(qū)最大等效應(yīng)力、滑移幅值、法向應(yīng)變幅值等微動疲勞參數(shù)的影響。

Peteghem等[15]研發(fā)了雙伺服液壓執(zhí)行器試驗系統(tǒng)利用微動磨損試驗，通過控制切向位移和法向載荷研究了恒載荷和變載荷條件下的磨損過程和磨損動力學(xué)。Fouvry等[16]討論了接觸區(qū)尺寸、頻率和循環(huán)法向載荷對微動裂紋形狀和磨損率的影響。Abbasi等[17-18]研究了接觸壓力以及軸向和法向接觸載荷之間的相位差對循環(huán)法向接觸載荷條件下Al7075-T6 微動疲勞行為的影響。Liao等[19]研究了合金690在285 ℃純水中不同法向載荷下的微動腐蝕疲勞行為，研究表明隨著法向載荷的增加，疲勞壽命降低，表面損傷和亞表面變形更加嚴(yán)重。Sun等[20]開發(fā)了一種新型的高溫摩擦疲勞試驗裝置，研究了鎳基單晶（NBSX）超合金在高溫下的摩擦疲勞機(jī)制，在600 ℃下進(jìn)行了5種加載條件的疲勞試驗，發(fā)現(xiàn)摩擦疲勞壽命隨著軸向載荷和法向載荷的增加而降低。

目前大部分的研究工作都集中在可變切向力與軸向載荷的組合上，然而，許多實際工況下，廣泛存在著可變法向載荷和軸向載荷的多軸加載條件，因此，本文在傳統(tǒng)微動疲勞模型分析的基礎(chǔ)上，討論了受循環(huán)法向載荷和軸向載荷作用的修正微動疲勞模型的接觸行為，重點研究了法向載荷和軸向載荷的相位差、微動墊夾具（“彈簧”）剛度等對微動動力學(xué)和接觸狀態(tài)的影響，提出了一種曲線分析方法，研究復(fù)雜多軸交變載荷條件下微動運行機(jī)制。

1 微動疲勞研究模型及動力學(xué)分析

1.1 研究模型概述

微動疲勞是燕尾接頭、螺栓接頭、花鍵接頭等機(jī)械接頭的主要失效形式之一[21-22]。幾十年來，許多研究都集中在這個問題上，并逐漸形成了一個統(tǒng)一的研究模式。如圖1a所示，2個圓柱形微動墊靠在一個扁平的試樣上，他們之間通過一個恒定的法向載荷形成接觸。同時，微動墊的固定裝置可以被當(dāng)作一個“彈簧”來抵抗微動墊的運動。然后，在“彈簧”和施加在試件下端部的軸向載荷的共同作用下，會產(chǎn)生一個可變的切向力，因此，在微動墊和試件的接觸區(qū)域中產(chǎn)生了微動效應(yīng)，該模型已被廣泛用于研究微動疲勞和磨損機(jī)理，基于該模型，研究者們進(jìn)行了一系列的試驗和仿真。但是，該模型仍有一定的局限性，不能完全反映實際工況下的微動條件。

圖1 當(dāng)前的微動疲勞模型（a），實際航空發(fā)動機(jī)燕尾榫工況的載荷模型（b）

一個典型的例子如航空發(fā)動機(jī)燕尾接頭的微動問題，如圖1b所示。在飛行中，發(fā)動機(jī)榫連接處的法向載荷不斷變化，與恒定法向加載條件下相比，實際微動系統(tǒng)的動力學(xué)和接觸應(yīng)力分布有很大的不同。

Barber等[23]研究了周期性法向載荷和切向載荷下的接觸情況，并給出了分析接觸狀態(tài)的解析解。Hojjati-Talemi等[24]建立了有限元模型，分析了變法向載荷和切向力作用下的應(yīng)力分布。

1.2 微動動力學(xué)分析

微動疲勞模型各部分的位移如圖2所示，其中“彈簧”（微動墊夾具）的剛度為f，試件的剛度為s。點和點是在接觸區(qū)中心的2個對應(yīng)點。這兩點的位移量分別為和。s為試樣的整體位移。將定義為與點在試樣中的位置相關(guān)的參數(shù)，α=δ/s，則=αδs，1。切向力可描述為：

式中：為微動墊與試件的摩擦因數(shù)；為施加在微動墊上的法向載荷。

圖2 微動疲勞模型示意圖

Fig.2 Schematic diagram of fretting fatigue model

運動過程如圖3所示。如果在試樣的末端施加周期性的軸向載荷B，試件和“彈簧”開始變形。當(dāng)加載循環(huán)開始軸向載荷B較小時，接觸區(qū)中心處于黏滯狀態(tài)，即點和點具有相同的位移和切向力，f=kf。在、||<||的條件下，隨著B的增加，接觸狀態(tài)由部分滑移變?yōu)橥耆?。在部分滑移情況下，點和點始終具有相同的位移；在完全滑移情況下，微動墊在第一個循環(huán)中調(diào)整其相對位置，并在點移動到一個新的平衡位置，然后隨著軸向載荷的變化，其平衡位置將在第二個循環(huán)中保持穩(wěn)定或再次改變。

在上述分析中，沒有考慮法向荷載的變化，如果考慮到法向載荷的變化，接觸狀態(tài)將會非常復(fù)雜。由式（1）可知，當(dāng)接觸狀態(tài)在整個周期內(nèi)保持部分滑移時，切向力不受法向荷載變化的影響，仍與軸向荷載呈線性關(guān)系，而在法向載荷和軸向載荷的作用下，微動墊會移動到一個新的平衡位置或不斷變化，這種情況下的運行狀態(tài)模式與正常的恒定完全滑移條件相似，但更復(fù)雜。

圖3 試件與微動墊運動過程示意圖

2 微動疲勞有限元仿真分析

2.1 微動問題的接觸有限元模型

接觸區(qū)應(yīng)力應(yīng)變場分析是進(jìn)行微動疲勞損傷研究的重要基礎(chǔ)。通過接觸分析得到的接觸區(qū)的應(yīng)力應(yīng)變分布是裂紋萌生位置、裂紋初始擴(kuò)展方向和微動疲勞損傷的主要控制參數(shù)[25-26]。

由于多軸加載條件下的微動過程非常復(fù)雜，可以建立有限元分析模型來輔助分析這一問題。本文利用有限元軟件ABAQUS建立了二維（平面應(yīng)變）模型（見圖4），模擬了微動墊和試件在循環(huán)法向載荷和軸向載荷下的接觸狀態(tài)。模型分為3個主要部分：微動墊、試件和橫向彈簧。圓柱體的半徑為100 mm，微動墊夾具模擬為橫向彈簧，通過為微動墊夾具材料定義小的彈性模量，橫向彈簧具有低剛度，橫向彈簧左邊緣節(jié)點的所有自由度都是固定的，這限制了微動墊的剛性運動。對于平面試件，左邊緣節(jié)點的方向自由度和下邊緣節(jié)點的方向自由度是固定的，以限制試件的剛性運動。然后，在微動墊的上邊緣施加法向載荷，在試件的右邊緣施加軸向載荷?？勺兎ㄏ蜉d荷和軸向載荷見式（2）。

圖4 多軸加載有限元分析模型

式中：m是法向載荷的均值；a是法向載荷的振幅；B,m是軸向載荷的均值；B,a是軸向載荷的振幅；是法向載荷和軸向載荷之間的相位差。

接觸分析為高度非線性問題，接觸區(qū)有限元網(wǎng)格的質(zhì)量和單元尺寸直接影響到計算的收斂性和精度。因此對接觸區(qū)的網(wǎng)格進(jìn)行了加密，使用尺寸為10 μm× 10 μm的正方形網(wǎng)格可以滿足計算要求。將微動墊和試件的對應(yīng)區(qū)域定義為接觸，接觸行為遵循庫侖摩擦定律，定義滑動摩擦因數(shù)為=0.65。墊片和試件的材料均為2024-T351鋁合金，其彈性模量=74 100 MPa，泊松比=0.33。

2.2 有限元模擬的加載分析過程及參數(shù)設(shè)置

載荷分三步施加到模型上。第一步，施加一個小的法向載荷來建立微動墊和試件之間的接觸；第二步，將法向載荷的均值a應(yīng)用于模型；第三步，施加循環(huán)軸向載荷B和法向載荷來模擬多軸循環(huán)荷載工況。利用ABAQUS子程序DLOAD對模型施加動態(tài)載荷。由于微動狀態(tài)在第二個循環(huán)中是穩(wěn)定的，因此模擬了2個加載循環(huán)，并將每個循環(huán)分為100個子步驟進(jìn)行詳細(xì)分析。

該模擬的主要目的是研究可變法向載荷和相位差對接觸行為的影響，因此在所有模擬情況下都使用了一個水平的軸向載荷，應(yīng)力比=?1，B,a=100 MPa，以及循環(huán)法向載荷m=225.2 N/mm，a=77 N/mm（對應(yīng)從140 MPa到200 MPa的最大接觸應(yīng)力0）。采用3個相位差（=0°、45°、90°）和2個彈簧剛度（f=2 084、3 770 N/mm，分別對應(yīng)于彈簧的彈性模量模量s= 1 990、3 600 MPa）來評估多軸加載對微動運行狀態(tài)的影響。

3 有限元結(jié)果分析

根據(jù)2.2節(jié)中的加載設(shè)置，進(jìn)行了一系列有限元仿真模擬，輸出橫向彈簧左側(cè)邊緣的反作用力（如圖4所示）等于切向力，以便進(jìn)一步分析。

3.1 不同載荷條件下的Ft-D曲線特征

t曲線（t是切向力，是相對位移幅值）的軸為，軸為，是分析微動狀態(tài)最重要的動力學(xué)信息[27]。Vingsbo等[28]提出了“微動圖”的概念，闡明微動磨損參數(shù)與微動磨損之間的關(guān)系。Zhou等[29]系統(tǒng)地研究了t曲線與微動損傷之間的關(guān)系，提出了微動圖理論來解釋微動機(jī)理。他們提出t曲線具有3種典型的形狀（見圖5）：（1）線性t曲線，表示接觸中心處于黏著狀態(tài)，接觸邊緣出現(xiàn)微滑移，對應(yīng)部分滑移的微動運行狀態(tài)；（2）橢圓t曲線，表示接觸體間的相對運動由彈塑性變形來協(xié)調(diào)，但微動運行狀態(tài)仍為部分滑移；（3）平行四邊形的t曲線，對應(yīng)全局滑動的微動運行狀態(tài)[30]。這3種類型的t曲線基于恒定的法向加載條件。然而，隨著法向載荷的變化，t曲線會發(fā)生顯著變化。

圖5 典型的Ft-D曲線

如圖6所示，在剛度為2 084 N/mm、相位差為0°時，切向力加載曲線幾乎是線性的，運行狀態(tài)保持為部分滑移狀態(tài)，與圖5a相似。如圖7所示，f=3 770 N/mm、=0°時，由于剛度增加，在開局有不穩(wěn)定的滑移階段后，曲線穩(wěn)定為平行四邊形，表明此加載下有周期性的完全滑移階段，與圖5c相似。

圖6 微動墊夾具剛度較小時比例加載條件下的Ft-D曲線（kf =2 084 N/mm，φ= 0°）

f=3 770 N/mm、= 90°加載時，如圖8所示，在可變法向載荷和相位差的影響下，曲線兩端形狀發(fā)生了扭曲，近似處于梯形，產(chǎn)生這種狀態(tài)的原因在于，與恒定的法向載荷且無相位差的加載條件相比，完全滑移階段的切向力不再保持恒定，所以曲線兩端形成了不規(guī)則的形狀。

圖7 微動墊夾具剛度較大時比例加載條件下的Ft-D曲線（kf=3 770 N/mm，φ= 0°）

圖8 微動墊夾具剛度較大時非比例加載條件下的Ft-D曲線（kf=3 770 N/mm，φ= 90°）

3.2 Q-P曲線分析方法及討論

由于在加載參數(shù)可變的多軸載荷條件下，運行狀態(tài)非常復(fù)雜且難以分析，所以在本節(jié)中提出了一種-曲線分析方法來輔助分析這些復(fù)雜微動情況下的運行狀態(tài)特性。

3.2.1曲線分析模型

切向力歷程和t-曲線提供的信息對于多軸載荷條件下的微動運行狀態(tài)分析是有限的。基于阿蒙頓-庫侖定律，發(fā)現(xiàn)和的函數(shù)可以提供更多關(guān)于微動運行狀態(tài)的有用信息。

如果接觸狀態(tài)在加載循環(huán)期間保持部分滑移狀態(tài)，和的函數(shù)可設(shè)為：

顯然，和的關(guān)系服從橢圓方程。如果把線性曲線=±和-橢圓曲線放在同一個坐標(biāo)系中，它們之間有5個典型的條件（見圖9）：（1）不相交；（2）在一點相切；（3）在兩點相切；（4）在兩點相交；（5）在三點相交；（6）在四點相交。橢圓的中心坐標(biāo)、長軸和短軸的長度以及傾斜角影響線性曲線和橢圓曲線之間的關(guān)系。而所有這些因素都是由方程中的加載參數(shù)控制的。

在上面的討論中，沒有考慮滑動摩擦力的影響，而實際上，當(dāng)橢圓曲線越過=±線時，曲線的形狀會發(fā)生變化，即線性曲線=±是部分滑移與完全滑移的分界線。在實際情況下，橢圓和曲線的交叉部分不具有物理意義，這些部分的載荷路徑（圖9d—f中的虛線部分）會變成與線性曲線=±重合的線性曲線。因此，圖9a—c為部分滑移區(qū)，圖9d—f有穩(wěn)定的完全滑移區(qū)。

3.2.2 比例加載條件下的微動運行狀態(tài)

在比例加載條件下，-橢圓函數(shù)可以化為線性函數(shù)，定義載荷幅值比為amp=B,a/a，則函數(shù)進(jìn)一步可化為：

從公式（5）可以看出，線性曲線有一個中心坐標(biāo)(m,stiB,m)，斜率為stiamp，即m決定曲線的水平位置，sti、B,m和決定曲線的垂直位置，曲線的旋轉(zhuǎn)角度由sti、amp和決定。

如圖10所示，2種剛度曲線均為線性函數(shù)，并且隨著彈簧剛度增加，曲線的斜率增加，小剛度曲線完全包含在內(nèi)，則其完全為部分滑移，而大剛度的曲線與兩條邊都有相交，說明其運行狀態(tài)有完全滑移。

圖9 直線Q=±μP和Q-P橢圓曲線之間的可能位置

圖10 相位差φ=0°多軸加載條件下的Q-P曲線

3.2.3 非比例加載條件下的微動運行狀態(tài)

在非比例加載條件下，橢圓軸與坐標(biāo)系不平行，函數(shù)曲線會變得復(fù)雜。有相位差時，根據(jù)式（5），曲線可化為標(biāo)準(zhǔn)橢圓函數(shù)：

橢圓方程的中心坐標(biāo)為(m,stiB,m)，長軸由a決定，短軸由sti、B,a和決定。相位差為45°時，如圖11所示，大剛度曲線的短軸和長軸都比小剛度的大，且有一定的旋轉(zhuǎn)，小剛度橢圓始終保持在內(nèi)部，大剛度時橢圓發(fā)生逆時針旋轉(zhuǎn)以及短軸和長軸增加導(dǎo)致與=±相交于四點，形成穩(wěn)定的完全滑移階段。

圖11 相位差=45°多軸加載條件下的曲線

Fig.11 The-curve under multi-axial loading condition with phase difference=45°

如圖12所示，在加載相位為90°時，可以看出當(dāng)彈簧剛度增加時，橢圓的長軸和中心點幾乎保持不變，但是短軸明顯增加，所以能得到在此相位時剛度增加的影響主要在于使短軸增加，導(dǎo)致大剛度橢圓與=±上下兩邊都有相交，發(fā)生了完全滑移。

圖12 相位差φ=90°多軸加載條件下的Q-P曲線

4 試驗驗證

在PWS-25G疲勞試驗機(jī)上進(jìn)行了一系列試驗，以評估曲線分析方法的合理性。采用加載條件（B,a=10 kN，B,m=0 kN，a=1.54 kN，m=?4.5 kN，=0°、45°、90°）進(jìn)行微動疲勞試驗。相位差為0°時，圖13a、圖14a所示，接觸區(qū)只存在少量的表面塑性變形。相位差為45°時，圖13b接觸區(qū)域均覆蓋黑色的氧化層，并且從圖14b中可以看出接觸區(qū)存在大量的塑性劃痕及變形。相位差為90°時，如圖13c所示，接觸區(qū)域氧化層較圖13b顯著增加，圖14c中可以看出磨損的程度很深，范圍大。

圖13 不同相位差微動疲勞試驗?zāi)p狀態(tài)

圖14 磨損表面的三維形貌

結(jié)合圖15所示的-圖，可對上述微動狀態(tài)進(jìn)行進(jìn)一步的討論。=0°的加載條件下，-曲線近似直線，=45°時曲線近似為橢圓，處在=±曲線內(nèi)部。值得注意的是，其與= ?近似相切。=90°時，直線部分與=±曲線平行，其余部分近似為橢圓曲線。根據(jù)本文提出的分析方法可推斷，在相位差為0°和45°時，接觸狀態(tài)為部分滑移，相位差為90°時，接觸狀態(tài)為完全滑移。反映到圖13、圖14所示的微動接觸區(qū)形態(tài)上，可判斷圖13a所示的微動處在部分滑移狀態(tài)；圖13c所示的微動處在完全滑移狀態(tài)。圖13b較為特殊，從接觸區(qū)形貌分析，可見其既有塑性變形又有較多的氧化層（但較圖13c少）。結(jié)合-圖推測，因其-曲線與=?近似相交，在微動運行的有限局部（切點附近）接觸狀態(tài)會短時為完全滑移，進(jìn)而產(chǎn)生了較多的磨屑，但總體上該狀態(tài)仍為部分滑移。

圖15 不同相位差下微動疲勞試驗的Q-P曲線

5 結(jié)論

1）對于部分滑移狀態(tài)，微動墊和試件保持靜摩擦狀態(tài)，微動墊的位移與試樣同步，由于切向力是由微動墊夾具（“彈簧”）的變形產(chǎn)生的，因此它與相位差和法向載荷的變化無關(guān)，僅取決于軸向載荷。

2）對于完全滑移狀態(tài)，t曲線不是規(guī)則的平行四邊形，這是由于法向荷載以及相位差的變化，導(dǎo)致曲線變得極其復(fù)雜。反映到曲線上，表現(xiàn)為橢圓與=±相交，產(chǎn)生了完全滑移。曲線變?yōu)橛刹糠謾E圓曲線與直線=±組合成的類橢圓曲線。由此可詳細(xì)地了解加載循環(huán)中不同時期的微動運行狀態(tài)。

3）在比例加載條件下，曲線為直線函數(shù)。曲線的水平位置與m呈正相關(guān)，垂直位置與sti、B,m和呈正相關(guān)，曲線的旋轉(zhuǎn)角度與sti、amp和呈正相關(guān)。

4）在非比例加載條件下，曲線是橢圓函數(shù)。m、sti、B,m和影響橢圓的中心坐標(biāo)，而a、sti、B,a和影響長軸和短軸的長度，影響橢圓的斜率。

[1] 楊茂勝, 陳躍良, 郁大照, 等. 微動疲勞研究的現(xiàn)狀與展望[J]. 強(qiáng)度與環(huán)境, 2008, 35(6): 45-54.

YANG Mao-sheng, CHEN Yue-liang, YU Da-zhao, et al. Several Key Problems in Studying on Fretting Fatigue[J]. Structure & Environment Engineering, 2008, 35(6): 45-54.

[2] 劉道新, 何家文. 微動疲勞影響因素及鈦合金微動疲勞行為[J]. 航空學(xué)報, 2001, 22(5): 454-457.

LIU Dao-xin, HE Jia-wen. Review of Factors that Influence Fretting Fatigue (FF) and Investigation on Ff Behavior of Ti-Alloy[J]. Acta Aeronautica et Astronautica Sinica, 2001, 22(5): 454-457.

[3] ARAúJO J A, NOWELL D. The Effect of Rapidly Varying Contact Stress Fields on Fretting Fatigue[J]. International Journal of Fatigue, 2002, 24(7): 763-775.

[4] ARAúJO J A, CASTRO F C, MATOS I M, et al. Life Prediction in Multiaxial High Cycle Fretting Fatigue[J]. International Journal of Fatigue, 2020, 134: 105504.

[5] WATERHOUSE R B. Fretting fatigue[M]. London: Applied Science Publishers, 1981

[6] 何明鑒. 機(jī)械構(gòu)件的微動疲勞[M]. 北京: 國防工業(yè)出版社, 1994.

HE Ming-jian. Fretting Fatigue of Mechanical Comp-onents[M]. Beijing: National Defense Industry Press, 1994.

[7] NOWELL D, DINI D, HILLS D A. Recent Developments in the Understanding of Fretting Fatigue[J]. Engineering Fracture Mechanics, 2006, 73(2): 207-222.

[8] RUIZ C, BODDINGTON P B, CHEN K C. An Invest-igation of Fatigue and Fretting in a Dovetail Joint[J]. Experimental Mechanics, 1984, 24(3): 208-217.

[9] HOEPPNER D, ADIBNAZARI S, MOESSER M W. Literature Review and Preliminary Studies of Fretting and Fretting Fatigue Including Special Applications to Aircraft Joints[J]. Fracture Mechanics, 1994

[10] 古遠(yuǎn)興. 高低周復(fù)合載荷下燕尾榫結(jié)構(gòu)微動疲勞壽命研究[D]. 南京: 南京航空航天大學(xué), 2007.

GU Yuan-xing. Research on Fretting Fatigue Life of Dovetail Joint under HCF-LCF Load[D]. Nanjing: Nan-jing University of Aeronautics and Astronautics, 2007.

[11] 徐友良, 崔海濤, 陳偉. 基于臨界面法的燕尾榫連接結(jié)構(gòu)微動疲勞壽命預(yù)測[J]. 航空動力學(xué)報, 2013, 28(3): 489-493.

XU You-liang, CUI Hai-tao, CHEN Wei. Fretting Fatigue Life Prediction Based on Critical Plane Approach for Dovetail Joint[J]. Journal of Aerospace Power, 2013, 28(3): 489-493.

[12] 米雪, 唐攀, 沈平川, 等. 690合金管在不同法向載荷下的切向微動磨損性能研究[J]. 表面技術(shù), 2020, 49(11): 191-197.

MI Xue, TANG Pan, SHEN Ping-chuan, et al. Tangential Fretting Wear Characteristics of 690 Alloy Tubes under Different Normal Force[J]. Surface Technology, 2020, 49(11): 191-197.

[13] 景鵬飛, 俞樹榮, 宋偉, 等. 接觸載荷對TC4鈦合金微動磨損行為的影響[J]. 表面技術(shù), 2019, 48(11): 266-274.

JING Peng-fei, YU Shu-rong, SONG Wei, et al. Effect of Contact Load on Fretting Wear Behavior of TC4 Titanium Alloy[J]. Surface Technology, 2019, 48(11): 266-274.

[14] 石煒, 溫衛(wèi)東, 崔海濤. 軸向和法向應(yīng)力對DD3和DZ125鎳基合金微動疲勞性能的影響[J]. 機(jī)械工程材料, 2014, 38(11): 10-15.

SHI Wei, WEN Wei-dong, CUI Hai-tao. Effects of Axial and Normal Stresses on Fretting Fatigue Properties of DD3 and DZ125 Nickel-Based Alloys[J]. Materials for Mechanical Engineering, 2014, 38(11): 10-15.

[15] VAN PETEGHEM B, FOUVRY S, PETIT J. Effect of Variable Normal Force and Frequency on Fretting Wear Response of Ti-6Al-4V Contact[J]. Wear, 2011, 271(9-10): 1535-1542.

[16] FOUVRY S, ARNAUD P, MIGNOT A, et al. Contact Size, Frequency and Cyclic Normal Force Effects on Ti-6Al-4V Fretting Wear Processes: An Approach Combining Friction Power and Contact Oxygenation[J]. Tribology International, 2017, 113: 460-473.

[17] ABBASI F, MAJZOOBI G H. Effect of Contact Pressure on Fretting Fatigue Behavior under Cyclic Contact Loading[J]. Surface Review and Letters, 2017, 24(Supp02): 1850032.

[18] ABBASI F, MAJZOOBI G H. Effect of Out-of-Phase Loading on Fretting Fatigue Response of Al7075-T6 under Cyclic Normal Loading Using a New Testing Apparatus[J]. Engineering Fracture Mechanics, 2018, 188: 93-111.

[19] LIAO Jia-peng, TAN Ji-bo, WU Xin-qiang, et al. Effects of Normal Load on Fretting Corrosion Fatigue of Alloy 690 in 285?℃ Pure Water[J]. Corrosion Science, 2018, 141: 158-167.

[20] SUN Shou-yi, LI Lei, HE Kun, et al. Fretting Fatigue Damage Mechanism of Nickel-Based Single Crystal Superalloys at High Temperature[J]. International Journal of Mechanical Sciences, 2020, 186: 105894.

[21] CIAVARELLA M, DEMELIO G. A Review of Analytical Aspects of Fretting Fatigue, with Extension to Damage Parameters, and Application to Dovetail Joints[J]. International Journal of Solids and Structures, 2001, 38(10-13): 1791-1811.

[22] RAJASEKARAN R, NOWELL D. Fretting Fatigue in Dovetail Blade Roots: Experiment and Analysis[J]. Tribo-logy International, 2006, 39(10): 1277-1285.

[23] BARBER J R, DAVIES M, HILLS D A. Frictional Elastic Contact with Periodic Loading[J]. International Journal of Solids and Structures, 2011, 48(13): 2041- 2047.

[24] HOJJATI-TALEMI R, WAHAB M A, DE BAETS P. Finite Element Simulation of Phase Difference Effects on Fretting Fatigue Crack Nucleation Behaviour[J]. Procee-dings of the Institution of Mechanical Engineers, Part J: Journal of Engineering Tribology, 2014, 228(4): 470-479.

[25] 趙華, 金雪巖, 朱民昊, 等. 微動接觸應(yīng)力的數(shù)值分析[J]. 四川大學(xué)學(xué)報(工程科學(xué)版), 2003, 35(5): 32-36.

ZHAO Hua, JIN Xue-yan, ZHU Min-hao, et al. Nume-rical Analysis of Fretting Contact Stresses[J]. Journal of Sichuan University (Engineering Science Edition), 2003, 35(5): 32-36.

[26] 趙華, 周仲榮. 數(shù)值方法在微動疲勞研究中的應(yīng)用進(jìn)展[J]. 摩擦學(xué)學(xué)報, 2000, 20(4): 317-320.

ZHAO Hua, ZHOU Zhong-rong. Application and Progress of Numerical Methods in Analysis of Fretting Fatigue[J]. Tribology, 2000, 20(4): 317-320.

[27] 周仲榮, Leo Vincent. 微動磨損[M]. 北京: 科學(xué)出版社, 2002.

ZHOU Zhong-rong. Fretting wear[M]. Beijing: Science Press, 2002.

[28] VINGSBO O, S?DERBERG S. On Fretting Maps[J]. Wear, 1988, 126(2): 131-147.

[29] ZHOU Z R, VINCENT L. Cracking Induced by Fretting of Aluminium Alloys[J]. Journal of Tribology, 1997, 119(1): 36-42.

[30] ZHOU Z R, VINCENT L. Mixed Fretting Regime[J]. Wear, 1995, 181-183: 531-536.

Fretting Running State under Complex Alternating Loads

(School of Mechanical-electronic and Vehicle Engineering, Beijing University of Civil Engineering and Architecture, Beijing 100044, China)

The present common-used fretting fatigue models simplify the actual loading conditions of engineering cases and ignore the effect of variable normal load. In this paper, the contact behavior of the fretting fatigue model subjected to variable normal load and axial load was investigated. The running state of fretting under alternating load was studied. Acurve analysis method was proposed. The contact state under typical loading conditions was discussed, and the influence of different loading parameters on the fretting running state was revealed. It shows that the proposedcurves analysis method can provide guidance for the study of the fretting running state under actual engineering working conditions.

A two-dimensional (plane strain) model was established by using the finite element software ABAQUS to simulate the contact state of the fretting pad and the specimen under cyclic normal load and axial load. The model was divided into 3 main parts: the fretting pad, the specimen and the transverse spring. The normal load was applied to the upper edge of the fretting pad and the axial load was applied to the right edge of the specimen. The main purpose of this simulation was to study the effect of variable normal load and phase difference on the contact behavior, so a bulk axial load with stress ratio= ?1,B,a=100 MPa and cyclic normal loadm=225.2 N/mm anda=77 N/mm (which corresponds to maximum contact stress0from 140 MPa to 200 MPa). Three phase difference=0°, 45°, 90° and two spring stiffnessf=2 084 N/mm, 3 770 N/mm (corresponding to the Young’s modulus of the springs= 1 990 MPa, 3 600 MPa, respectively) were used to evaluate the effect of multi-axial loading on the micro-action fatigue operating condition, with the reaction force at the left edge of the output transverse spring equal to the tangential force.

A series of fretting fatigue test corresponding to the simulation was performed on a biaxial fretting fatigue testing system The biaxial fretting fatigue testing system has two hydraulic actuators and three load cells, through which the loading of variable normal loads, bulk loadBand the measurement of tangential forcescould be implemented. A series of tests were performed using loading conditions (B,a=10 kN;B,m=0 kN;a=1.54 kN;m= ?4.5 kN;=0°, 45°, 90°) and the loading loops with different phase differences were recorded to evaluate the rationality of thecurve analysis method.

For proportional loading conditions, thetcurve is a straight line for a small pad fixture stiffness, the curve changes to a regular parallelogram for a larger stiffness. In these conditions, thecurves are linear functions, and the slope of the curve increases as the increase of spring stiffness.

For non-proportional loading conditions, the shape of thetcurve is no longer a parallelogram and becomes very peculiar. Thecurves obey the ellipse function. The shape of the ellipse is determined by pad fixture stiffness, normal and bulk loading parameters, the phase difference etc.

The proposedcurves analysis method can provide a guidance for the study of the fretting fatigue running state, and also provides a method for further discussion of fretting fatigue or fretting wear behavior.

fretting; complex alternating load; dynamics; contact state; finite element simulation

O346.2

A

1001-3660(2023)01-0112-09

10.16490/j.cnki.issn.1001-3660.2023.01.012

2021–12–28；

2022–04–08

2021-12-28；

2022-04-08

國家自然科學(xué)基金青年基金項目（51905028）；北京市教育委員會科技計劃一般項目（KM202110016002）；北京建筑大學(xué)金字塔人才培養(yǎng)工程（JDYC20200323）

National Natural Science Foundation of China Youth Fund Project (51905028); Beijing Municipal Education Commission Science and Technology Plan General Project (KM202110016002); Beijing University of Civil Engineering and Architecture Pyramid Talent Training Project (JDYC20200323)

田帥（1996—），男，碩士生，主要研究方向微動疲勞損傷機(jī)理

TIAN Shuai (1996-), Male, Postgraduate, Research focus: fretting fatigue damage regime

李欣（1984—），男，博士，副教授，主要研究方向為機(jī)械結(jié)構(gòu)強(qiáng)度及疲勞耐久性分析

LI Xin (1984-), Male, Doctor, Associate professor, Research focus: mechanical structure strength and fatigue durability analysis

田帥, 李欣, 楊建偉, 等.復(fù)雜交變載荷條件下微動運行狀態(tài)[J]. 表面技術(shù), 2023, 52(1): 112-120.

TIAN Shuai, LI Xin, YANG Jian-wei, et al. Fretting Running State under Complex Alternating Loads[J]. Surface Technology, 2023, 52(1): 112-120.

責(zé)任編輯：萬長清