董　赟　馮瑞成　剡昌鋒　段早琦

1.東南大學(xué),南京,211189　　2.蘭州理工大學(xué),蘭州,730050

?

一種計(jì)算輪轂疲勞損傷的數(shù)值模擬新算法

董赟1,2馮瑞成2剡昌鋒2段早琦1

1.東南大學(xué),南京,2111892.蘭州理工大學(xué),蘭州,730050

摘要：針對(duì)輪轂結(jié)構(gòu)復(fù)雜和周期載荷的特點(diǎn),采用剛性元法簡(jiǎn)化輪轂有限元模型, 基于材料疲勞損傷規(guī)律、應(yīng)變等效假設(shè)理論和有限元方法,提出了一種計(jì)算旋轉(zhuǎn)輪轂結(jié)構(gòu)疲勞損傷的新算法。該算法采用擬合得到的損傷演化材質(zhì)參數(shù)和ANSYS的二次開發(fā)工具對(duì)輪轂的疲勞損傷過程進(jìn)行數(shù)值計(jì)算,在得到裂紋全壽命的同時(shí)，給出裂紋群擴(kuò)展路徑，真實(shí)地模擬了輪轂動(dòng)態(tài)旋轉(zhuǎn)損傷進(jìn)程。試驗(yàn)研究證明了該方法的可行性和有效性。

關(guān)鍵詞:輪轂建模;旋轉(zhuǎn)加載;累計(jì)損傷;疲勞壽命

0引言

在實(shí)際工作狀態(tài)下,輪轂承受旋轉(zhuǎn)交變載荷，疲勞破壞是輪轂的主要破壞形式之一，因此，預(yù)測(cè)輪轂的疲勞損傷具有十分重要的科學(xué)與工程意義。目前有關(guān)輪轂疲勞壽命預(yù)測(cè)的方法有名義應(yīng)力法、局部應(yīng)力應(yīng)變法、能量法和場(chǎng)強(qiáng)法等基于應(yīng)力應(yīng)變疲勞理論的預(yù)測(cè)方法[1-7]。學(xué)者們多年的研究取得了一定的成果，對(duì)輪轂的疲勞壽命預(yù)測(cè)起了重要的作用，但這些方法均假設(shè)疲勞損傷累積到一定程度便發(fā)生斷裂失效,未能全面描述疲勞裂紋從微觀缺陷發(fā)展至宏觀裂紋直至斷裂的規(guī)律,所得的結(jié)果與實(shí)際壽命之間存在誤差,難以真正應(yīng)用于工程實(shí)際。而損傷力學(xué)理論可以解釋材料或構(gòu)件在交變載荷作用下裂紋萌生和擴(kuò)展的規(guī)律[8]。因此，基于材料的疲勞損傷進(jìn)程，研究輪轂疲勞壽命的預(yù)測(cè)是一種新的途徑。

本文針對(duì)輪轂結(jié)構(gòu)和工況的特點(diǎn),利用ANSYS完成了輪轂的三維幾何建模、網(wǎng)格劃分及載荷約束,通過Cerig命令在加載軸端建立一主節(jié)點(diǎn)，將該節(jié)點(diǎn)與其他螺栓孔受力節(jié)點(diǎn)通過剛性梁耦合,然后在主節(jié)點(diǎn)上加旋轉(zhuǎn)載荷,可以簡(jiǎn)化和模擬輪轂的實(shí)際工作過程。在疲勞壽命分析時(shí),從考慮材料疲勞損傷進(jìn)程的角度出發(fā)，結(jié)合疲勞損傷理論和有限元分析方法，建立了適合計(jì)算輪轂疲勞壽命的損傷力學(xué)-有效應(yīng)力-有限元法,采用擬合得到的損傷演化材質(zhì)參數(shù),利用ANSYS的二次開發(fā)工具APDL語言，對(duì)輪轂疲勞損傷過程進(jìn)行數(shù)值模擬計(jì)算，在得到輪轂裂紋萌生壽命和擴(kuò)展壽命的同時(shí)，給出了不同循環(huán)載荷下裂紋的擴(kuò)展路徑。最后，通過對(duì)輪轂彎曲疲勞試驗(yàn)結(jié)果與有限元數(shù)值仿真結(jié)果進(jìn)行對(duì)比驗(yàn)證，證明了此方法的可行性和有效性。

1輪轂有限元模型的簡(jiǎn)化

輪轂是由鋼板沖壓成形后的輪輞和輪輻經(jīng)焊接而成的薄壁結(jié)構(gòu)；試驗(yàn)對(duì)象為帶加載軸、法蘭盤及螺栓的整體輪轂。由于輪轂結(jié)構(gòu)比較復(fù)雜，所含的單元數(shù)較多，在強(qiáng)度計(jì)算過程中,計(jì)算速度較低。本文考慮了輪轂工作過程中受到的彎矩作用，通過剛性元法對(duì)輪轂?zāi)Ｐ瓦M(jìn)行了簡(jiǎn)化。

由于輪轂在實(shí)際工況下產(chǎn)生的疲勞為高周疲勞,因此輪轂發(fā)生的變形基本都是小變形。結(jié)合輪轂幾何模型和實(shí)際加載情況,在距輪轂輪輻面中心540 mm處建立一主節(jié)點(diǎn),定義為mass21剛性單元,然后通過Cerig命令把該節(jié)點(diǎn)與輪輻螺栓接觸面的所有節(jié)點(diǎn)耦合,再對(duì)輪轂輪輞進(jìn)行邊界約束,最后給主節(jié)點(diǎn)(加載軸端)施加額定載荷(F=3400 N),通過剛性單元進(jìn)行載荷傳遞。這種方法在端面額外增加了一定的剛度，一般適用于小變形分析。簡(jiǎn)化后的模型加載圖見圖1。圖2為受載后的應(yīng)力分布圖。單元在加載過程中應(yīng)具有承受大變形、大應(yīng)變的能力，并能在不喪失精確性的前提下忍受不規(guī)則形狀，結(jié)合單元“殺死”功能,本文選取SHELL143單元。

圖1　受Y負(fù)方向載荷的輪轂加載圖

圖2　受Y負(fù)方向載荷的輪轂應(yīng)力分布

2輪轂疲勞損傷

輪轂的疲勞損傷實(shí)質(zhì)是微裂紋的萌生及發(fā)展過程。在旋轉(zhuǎn)循環(huán)載荷作用下，材料的局部細(xì)觀組織發(fā)生了不可逆的塑性變形，非線性累計(jì)損傷增加，疲勞強(qiáng)度顯著減小。根據(jù)以上損傷特征，Chaboche[9]提出一個(gè)疲勞損傷累積模型來描述材料的劣化過程：

(1)

(2)

(3)

(4)

3損傷力學(xué)-有效應(yīng)力-有限元法

3.1損傷力學(xué)-有效應(yīng)力-有限元法計(jì)算格式的建立

由式(1)可知,每加載一次,第i單元產(chǎn)生的損傷增量為

(5)

算法具體步驟如下：

(6)

此時(shí),單元總損傷為

(7)

(8)

(9)

(4)重復(fù)步驟(3)，第r次加載后的有效應(yīng)力為

(10)

單元總損傷為

(11)

當(dāng)損傷最大單元總損傷

(12)

時(shí)單元發(fā)生破壞,此時(shí)“殺死”該單元。其中，Nf1為首個(gè)單元破壞時(shí)所對(duì)應(yīng)的加載次數(shù)。

此時(shí)第i單元的總損傷為

(13)

當(dāng)m個(gè)單元(本文m取36)被“殺死”時(shí),所對(duì)應(yīng)的加載次數(shù)為

(14)

由步驟(3)可知，對(duì)主節(jié)點(diǎn)沿α=15N°的方向上施加旋轉(zhuǎn)載荷時(shí)，在0°~360°范圍內(nèi)每加載24次,相當(dāng)于輪轂旋轉(zhuǎn)一周,所以輪轂的疲勞壽命可由下式表示:

f=[k]+1

式中,f為輪轂的疲勞壽命;M為周加載次數(shù),本文M取24;符號(hào)[]表示向下取整。

利用步驟(1)~步驟(4)的疲勞損傷有限元方法對(duì)采用剛性元簡(jiǎn)化后的輪轂進(jìn)行疲勞損傷預(yù)測(cè)過程中，由于施加旋轉(zhuǎn)載荷后計(jì)算的有效真應(yīng)力考慮了累計(jì)損傷值，并結(jié)合了“殺死”單元功能，從而有效描述了疲勞裂紋從微觀缺陷發(fā)展至宏觀裂紋直至斷裂的全進(jìn)程。

3.2算例驗(yàn)證

在用ANSYS對(duì)車輪軸端進(jìn)行加載時(shí),只能沿坐標(biāo)軸方向加載,所以應(yīng)對(duì)加載值分別沿Y軸和Z軸進(jìn)行載荷分解,即

(15)

式中,α為在OYZ平面內(nèi)加載方向從Y軸正方向開始逆時(shí)針旋轉(zhuǎn)的角位移;|F|為載荷矢量的模;FY、FZ分別為對(duì)Y軸和Z軸進(jìn)行分解的載荷值，值為負(fù)時(shí)表示加載方向和坐標(biāo)軸正方向相反。

表1　SAPH440鋼的疲勞損傷參數(shù)

圖3　受載方向與Y軸成0°時(shí)的加載圖和應(yīng)力分布

圖4　受載方向與Y軸成45°時(shí)的加載圖和應(yīng)力分布

圖5　受載方向與Y軸成165°時(shí)的加載圖和應(yīng)力分布

圖6　受載方向與Y軸成315°時(shí)的加載圖和應(yīng)力分布

通過對(duì)輪轂求解分析可知，當(dāng)旋轉(zhuǎn)加載時(shí)，大應(yīng)力單元主要分布在螺栓孔和通風(fēng)孔周圍，因此，通過應(yīng)力大小迭代求出單元損傷后，總損傷最先達(dá)到破壞值的單元應(yīng)在螺栓孔或通風(fēng)孔附近。

令θ為OYZ平面內(nèi)加載方向與Y軸正方向沿逆時(shí)針旋轉(zhuǎn)在0°~360°內(nèi)的夾角,則θ與角位移α之間的關(guān)系可表示為

α=(f-1)×360°+θ

(16)

圖7~圖12所示為不同加載角度θ、不同加載次數(shù)N下所對(duì)應(yīng)的損傷應(yīng)力場(chǎng)。

圖7　N=1 080 711(f=45 030周次，θ=255°)時(shí)的等效應(yīng)力圖

圖8　N=1 104 580(f=46 025周次,θ=60°)時(shí)的等效應(yīng)力圖

圖9　N=1 124 967 (f=46 874周次,θ=225°)時(shí)的等效應(yīng)力圖

圖10　N=1 175 434(f=48 977周次,θ=150°)時(shí)的等效應(yīng)力圖

圖11　N=1 252 026 (f=52 168周次,θ=270°)時(shí)的等效應(yīng)力圖

圖12　N=1 279 226 (f=53 302周次,θ=30°)時(shí)的等效應(yīng)力圖

通過對(duì)輪轂進(jìn)行旋轉(zhuǎn)循環(huán)加載,利用損傷有限元法進(jìn)行疲勞損傷計(jì)算可以看出,當(dāng)加載次數(shù)為1 080 711(疲勞壽命為45 030周次,θ=255°)時(shí)(圖7),輪轂螺栓孔處產(chǎn)生初始裂紋,即第1個(gè)單元發(fā)生破壞,此時(shí)的壽命為輪轂的萌生壽命；當(dāng)加載次數(shù)為1 104 580(疲勞壽命為46 025周次,θ=60°)時(shí)(圖8),初始裂紋進(jìn)行擴(kuò)展，有3個(gè)單元發(fā)生破壞；當(dāng)加載次數(shù)為1 124 967周次(疲勞壽命為46 874周次，θ=225°)時(shí)(圖9),通風(fēng)孔處有一單元發(fā)生破壞;當(dāng)加載次數(shù)為1 175 434(疲勞壽命f=48 977周次，θ=150°)時(shí)(圖10),裂紋進(jìn)一步擴(kuò)展，螺栓孔處有10個(gè)單元發(fā)生破壞，通風(fēng)孔處有7個(gè)單元發(fā)生破壞;當(dāng)加載次數(shù)為1 252 026次(疲勞壽命f=52 168周次,θ=270°)時(shí)(圖11),螺栓孔處又有一條新的裂紋產(chǎn)生，同時(shí)第一條裂紋擴(kuò)展至14個(gè)單元，第二條裂紋擴(kuò)展至15個(gè)單元；當(dāng)加載次數(shù)為1 279 226(疲勞壽命f=53 302周次,θ=30°)時(shí)(圖12),輪轂完全破壞,此時(shí)的疲勞壽命為輪轂的全壽命。

由計(jì)算結(jié)果可知，輪轂輪輻部位在產(chǎn)生裂紋的過程中，當(dāng)主裂紋產(chǎn)生并發(fā)展時(shí)，又有次生微裂紋產(chǎn)生，最后產(chǎn)生數(shù)條裂紋致使輪轂發(fā)生破壞。

4輪轂疲勞試驗(yàn)

根據(jù)GB/T5334-2005輪轂性能的要求和動(dòng)態(tài)彎曲疲勞試驗(yàn)的條件要求,將10只以熱軋帶鋼SAPH440生產(chǎn)的N1型鋼圈以及M12緊固螺栓、加載軸和法蘭在標(biāo)準(zhǔn)試驗(yàn)裝置[12]下做彎曲疲勞破壞試驗(yàn)。

4.1輪轂失效的判定準(zhǔn)則

隨著輪轂旋轉(zhuǎn)次數(shù)的增多,輪轂輪輻出現(xiàn)了破壞的跡象。根據(jù)國標(biāo)[13]中給出的判斷輪轂破壞的依據(jù)，輪轂在不能繼續(xù)承受載荷、原始裂紋產(chǎn)生擴(kuò)展或出現(xiàn)應(yīng)力導(dǎo)致輪轂斷面出現(xiàn)裂紋時(shí),可以斷定輪轂已經(jīng)失效。本試驗(yàn)在輪轂?zāi)巢课涣鸭y長(zhǎng)度為4 cm時(shí),停止加載,視為輪轂破壞。

4.2試驗(yàn)結(jié)果

對(duì)多個(gè)輪轂進(jìn)行破壞試驗(yàn),結(jié)果表明輪轂的破壞多發(fā)生在輪輻螺栓孔及通風(fēng)孔附近，部分輪轂的破壞圖見圖13和圖14。疲勞壽命試驗(yàn)結(jié)果和計(jì)算結(jié)果的比較如表2所示。

圖13　試樣1的裂紋分布

圖14　試樣3的裂紋分布

試樣號(hào)試驗(yàn)全壽命f1計(jì)算萌生壽命f2計(jì)算全壽命f3誤差|f3-f1|f1(%)1658232658003630004340005860006900007320008430009450001058000平均值58262450305330219191557384167241888.5

4.3試驗(yàn)結(jié)果與計(jì)算結(jié)果的比較分析

比較疲勞試驗(yàn)圖(圖13、圖14)和計(jì)算全壽命裂紋圖(圖12)可知，本文建立的損傷力學(xué)-有效應(yīng)力-有限元法對(duì)用剛性元簡(jiǎn)化后的輪轂?zāi)M計(jì)算的疲勞損傷裂紋與原型輪轂破壞試驗(yàn)結(jié)果非常相近，計(jì)算全壽命與原型輪轂試驗(yàn)全壽命平均值的誤差水平為8.5%，說明此方法是可行的。

通過對(duì)主節(jié)點(diǎn)的旋轉(zhuǎn)加載模擬了輪轂的動(dòng)態(tài)彎曲疲勞試驗(yàn)過程。輪轂在受到旋轉(zhuǎn)載荷加載時(shí)，危險(xiǎn)部位在螺栓孔和通風(fēng)孔附近，試驗(yàn)結(jié)果和計(jì)算結(jié)果吻合較好。

圖7、圖12及表2的結(jié)果表明,裂紋的萌生壽命超過全壽命的80%，這是因?yàn)檩嗇灣跏紦p傷很小，應(yīng)力分布比較均勻，初始裂紋產(chǎn)生后，同一截面上的單元幾乎同時(shí)被破壞,所以擴(kuò)展壽命非常短。

從理論上講,只要邊界條件正確，作用的載荷取得準(zhǔn)確，當(dāng)網(wǎng)格無限細(xì)化時(shí)，有限元解收斂于精確的解析解[11],但是由于輪轂結(jié)構(gòu)復(fù)雜，本文在劃分網(wǎng)格時(shí)考慮了計(jì)算速度問題，所以網(wǎng)格尺寸較粗，而當(dāng)單元?jiǎng)澐州^粗時(shí)，隨著單元性質(zhì)的劣化以及模型邊界條件的變化，會(huì)使部分單元的應(yīng)力場(chǎng)和應(yīng)變場(chǎng)發(fā)生突變，造成新的應(yīng)力集中，同時(shí)本文采用“殺死”后的單元模擬缺陷，由于單元“殺死”后沒有任何承載能力，與實(shí)際情況相比高估了損傷值，且本文在疲勞損傷有限元分析中沒有考慮材料的循環(huán)強(qiáng)化作用，所以以上這些都會(huì)使計(jì)算結(jié)果較試驗(yàn)結(jié)果偏小，導(dǎo)致表2中計(jì)算結(jié)果基本都是試驗(yàn)結(jié)果的下限。

仿真結(jié)果和試驗(yàn)結(jié)果存在的誤差可能還與以下因素有關(guān):①在計(jì)算輪轂應(yīng)力分析和疲勞壽命分析時(shí)采用彈性有限元模型，這與真實(shí)輪轂材料的本構(gòu)關(guān)系存在一定的誤差；②通過剛性元法對(duì)輪轂?zāi)Ｐ瓦M(jìn)行簡(jiǎn)化以及每隔15°對(duì)主節(jié)點(diǎn)進(jìn)行旋轉(zhuǎn)加載都與整體輪轂(帶法蘭盤和加載軸)的彎曲疲勞試驗(yàn)存在偏差；③本文沒有考慮旋轉(zhuǎn)的輪轂還會(huì)受到殘余應(yīng)力、螺栓預(yù)緊力以及離心力等作用的影響。

5結(jié)論

(1)由于輪轂在實(shí)際工況下產(chǎn)生的疲勞為高周疲勞，所以輪轂發(fā)生的變形基本都是小變形，結(jié)合輪轂幾何模型和實(shí)際加載情況，利用適合小變形情況的Cerig命令在輪轂輪輻面中心線上建立一主節(jié)點(diǎn)，然后把該節(jié)點(diǎn)與輪輻螺栓接觸面的所有節(jié)點(diǎn)進(jìn)行耦合，再對(duì)輪轂輪輞進(jìn)行邊界約束，最后給主節(jié)點(diǎn)(加載臂端)施加旋轉(zhuǎn)載荷，通過剛性單元進(jìn)行載荷傳遞，從而簡(jiǎn)化了輪轂整體的幾何模型。

(2)從考慮材料疲勞損傷進(jìn)程的角度出發(fā)，結(jié)合疲勞損傷理論和有限元分析方法，建立了適合計(jì)算輪轂等旋轉(zhuǎn)機(jī)構(gòu)疲勞壽命的損傷有限元方法。通過對(duì)加載臂端主節(jié)點(diǎn)進(jìn)行旋轉(zhuǎn)加載，將輪轂在工作過程中產(chǎn)生的損傷間接地用有效真應(yīng)力值表征，得到了輪轂的使用壽命和裂紋擴(kuò)展路徑。

(3)利用輪轂動(dòng)態(tài)彎曲疲勞試驗(yàn)?zāi)M輪轂在承受彎曲載荷作用下高速運(yùn)轉(zhuǎn)過程，將試驗(yàn)結(jié)果和數(shù)值仿真結(jié)果進(jìn)行對(duì)比，結(jié)果表明，本文建立的損傷力學(xué)-有效應(yīng)力-有限元法可以較準(zhǔn)確地預(yù)測(cè)輪轂的疲勞壽命和裂紋擴(kuò)展路徑。因此，結(jié)合損傷力學(xué)對(duì)輪轂疲勞壽命及破壞過程進(jìn)行預(yù)測(cè)是可行的。

參考文獻(xiàn):

[1]田軍.貨車車輪輻板孔裂紋及其運(yùn)用安全性研究[D].北京:北京交通大學(xué),2014.

[2]劉雄輝,李德源,呂文閣,等.水平軸風(fēng)力機(jī)輪轂強(qiáng)度和疲勞壽命分析方法研究[J].太陽能學(xué)報(bào),2012,33(5):763-768.

Liu Xionghui,Li Deyuan,Lü Wen’ge, et al.Research on Analysis Approach of Strength and Fatigue Life of Horizontal Axis Wind Turbine Hub[J].Acta Energiae Solaris Sinica,2012,33(5):763-768.

[3]Karandikar H M,Fuchs W .Fatigue Life Prediction for Wheels by Simulation of the Rotating Bending Test[J].SAE Paper,900140．

[4]Wang Ping,Cao Jiayong,Zhang Zhinan.Finite Element Analysis for Strength of Wheel Hubs of Large Wind Turbines[J].East China Electric Power,2009,37(7):1206-1208．

[5]胡衛(wèi)杰.鎂合金車輪彎曲疲勞壽命預(yù)測(cè)及動(dòng)態(tài)特性分析[D].秦皇島：燕山大學(xué)，2014.

[6]Kocabicak U,Firat M.A Simple Approach for Multiaxial Fatigue Damage FEM Postprocessing[J].Materials and Design,2004,25(6):73-82．

[7]Ridha R A.Finite Element Analysis of Automotive Wheels[J].SAE Paper,760085．

[8]董赟,蔡敢為,鄭戰(zhàn)光,等.基于損傷力學(xué)的疲勞壽命和裂紋擴(kuò)展的數(shù)值分析[J].中國機(jī)械工程,2010,21(20):2412-2415.

Dong Yun,Cai Ganwei,Zheng Zhanguang,et al.The Numerical Analysis of Fatigue Life and Crack Growth Based on Damage Mechanic[J].China Mechanical Engineering,2010,21(20):2412-2415.

[9]Chaboche J L.Continuum Damage Mechanics Part Ⅱ-Damage Growth,Crack Initiation and Crack Growth[J].ASME J. Appl. Mech.,1988,55(1):65-72.

[10]金磊,夏才初.理論流變力學(xué)模型中蠕變損傷的研究方法與問題[J].巖石力學(xué)與工程學(xué)報(bào),2012,31(1):3006-3014.

Jin Lei,Xia Caichu.Study Methods for Creep Damage in Theoretical Rheological Models and Some Problems[J].Chinese Journal of Rock Mechanics and Engineering,2012,31(1):3006-3014.

[11]李曉偉,黃佩珍.內(nèi)壓、外載下晶內(nèi)微裂紋演化的有限元分析[J].應(yīng)用力學(xué)學(xué)報(bào),2014,31(4):530-534.

Li Xiaowei,Huang Peizhen.The Evolution of Intragranular Microcracks under Internal Pressure and External Load Simulated by Finite Element Method[J].Chinese Journal of Applied Mechanics,2014,31(4):530-534.

[12]馬曉康，王宏祥，白 洋．汽車車輪動(dòng)態(tài)彎曲疲勞試驗(yàn)的有限元模擬[J].遼寧工業(yè)大學(xué)學(xué)報(bào),2014,34(6):386-389．

Ma Xiaokang,Wang Hongxiang,Bai Yang.Simulating Analysis of Fatigue Bending Test for Automobile Wheel[J].Simulating Analysis of Fatigue Bending Test for Automobile Wheel,2014,34(6):386-389.

[13]中華人民共和國汽車標(biāo)準(zhǔn)化技術(shù)委員會(huì).QC-T221-1997,汽車輕合金車輪的性能要求和試驗(yàn)方法[S].北京:中華人民共和國機(jī)械工業(yè)部,1997.

(編輯王艷麗)

A New Numerical Simulation Method for Calculating Fatigue Damage of Wheels

Dong Yun1,2Feng Ruicheng2Yan Changfeng2Duan Zaoqi1

1.Southeast University,Nanjing,211189 2.Lanzhou University of Technology,Lanzhou,730050

Abstract:According to the characteristics of the complicated structure and cyclic loading of wheels,a new numerical simulation method for calculating fatigue damage of wheels was proposed based on two ways,one was the fatigue damage coupling theory of strain equivalence and the finite element analysis method,the other was the law of material fatigue damage and the rigid method. Combining constants and properties of the fatigue damage model with fitting,a subroutine to calculate damage and to predict fatigue life on ANSYS platform was developed, while the initiation and growth of cracks of rim were simulated by damage developing, the crack propagation paths were given under the different cycles as the fatigue life was obtained.At last,the comparison between fatigue tests and calculation results indicates that this method is proper,feasible and of high accuracy.

Key words:wheel rim modeling;rotating load; accumulated damage; fatigue life

基金項(xiàng)目:國家自然科學(xué)基金資助項(xiàng)目(51165018);甘肅省自然科學(xué)基金資助項(xiàng)目(148RJZA008);甘肅省高等學(xué)?？蒲许?xiàng)目(2014A-033)

收稿日期：2015-03-04

中圖分類號(hào)：U463.343DOI：10.3969/j.issn.1004-132X.2015.23.011

作者簡(jiǎn)介：董赟,男,1984年生。東南大學(xué)機(jī)械工程學(xué)院博士研究生,蘭州理工大學(xué)機(jī)電工程學(xué)院講師。主要研究方向?yàn)榻Y(jié)構(gòu)強(qiáng)度、微納摩擦學(xué)。馮瑞成,男,1976年生。蘭州理工大學(xué)機(jī)電工程學(xué)院副教授。剡昌鋒,男,1974年生。蘭州理工大學(xué)機(jī)電工程學(xué)院教授、博士研究生導(dǎo)師。段早琦,男,1988年生。東南大學(xué)機(jī)械工程學(xué)院博士研究生。