熊壘，林有希

（福州大學(xué)機(jī)械工程及自動(dòng)化學(xué)院，福建福州 350108）

0 引言

磨損作為摩擦學(xué)的一個(gè)重要方面，許多學(xué)者已經(jīng)進(jìn)行了非常多的理論和實(shí)驗(yàn)研究［1-2］。目前關(guān)于磨損的研究手段主要以實(shí)驗(yàn)為主，但是由于摩擦學(xué)系統(tǒng)本身龐大而復(fù)雜，其中涉及的不可控制變量太多，相關(guān)實(shí)驗(yàn)需要耗費(fèi)大量的人力物力，而且實(shí)驗(yàn)的重復(fù)性和可比性較弱，這給定量研究磨損行為造成了很大困難［3］，所以有部分學(xué)者開始嘗試運(yùn)用模擬仿真的手段來進(jìn)行磨損的研究［4］。模擬仿真的優(yōu)點(diǎn)在于其變量易于控制，針對(duì)一種摩擦磨損現(xiàn)象，可以刻意的選擇相關(guān)參數(shù)和約束，從而避免其他因素的干擾。同時(shí)模擬仿真的重復(fù)性強(qiáng)，可比性好，更有利于有關(guān)規(guī)律和機(jī)理的深入研究［5］。本文綜述了模擬仿真技術(shù)在磨損研究中的應(yīng)用現(xiàn)狀，并就當(dāng)前存在的問題提出了建議。

1 兩種主要磨損的產(chǎn)生機(jī)理

1.1 粘著磨損

當(dāng)聚合物在金屬表面滑動(dòng)時(shí)產(chǎn)生粘著，粘著磨損的重要特征是軟的聚合物向硬質(zhì)金屬表面的轉(zhuǎn)移:開始時(shí)聚合物轉(zhuǎn)移到金屬上，接著聚合物本體與附著在金屬表面的聚合物轉(zhuǎn)移膜之間產(chǎn)生粘著作用，然后在剪切力的作用下，粘著點(diǎn)處的聚合物從本體上撕扯下來并積累在金屬表面，逐漸形成轉(zhuǎn)移膜，隨著轉(zhuǎn)移膜不斷變厚，同時(shí)在循環(huán)載荷的反復(fù)作用下，轉(zhuǎn)移膜破裂，一部分碎片脫離摩擦表面，一部分加入其他的磨損形式中，這樣就形成材料逐漸損耗的過程［6］。轉(zhuǎn)移膜的生成改變了摩擦接觸的類型，從金屬與聚合物的摩擦接觸變成聚合物與聚合物的摩擦接觸;同時(shí)也改變了金屬表面及形貌。一般來說較軟的物質(zhì)易在較硬的表面形成粘著磨損，如聚合物/金屬摩擦副;物理性質(zhì)相近或者相同者也易形成粘著磨損，如金屬/金屬摩擦副［7］。

1.2 磨粒磨損

外界硬質(zhì)顆?；蛘邔?duì)磨表面上的硬凸起物或粗糙峰在摩擦過程中引起表面材料脫落的現(xiàn)象稱為磨粒磨損。磨粒磨損很少是一種單一的磨損機(jī)制引起的，而經(jīng)常是多種磨損機(jī)制綜合作用的結(jié)果。而且磨損機(jī)制通常隨著工況條件的變化而變化。磨粒磨損過程中材料的去除機(jī)制主要分為以下兩種［8］:

a）由塑性變形機(jī)制引起的去除過程

1）犁溝效應(yīng)。在載荷作用下磨粒被壓入材料表面，同時(shí)向前運(yùn)動(dòng)，受到擠壓的材料產(chǎn)生塑性流動(dòng)，向兩側(cè)隆起。這種過程不會(huì)直接導(dǎo)致材料的去除，但是在磨粒的反復(fù)作用下，材料產(chǎn)生多次變形，最后便會(huì)脫落，形成磨損。

2）微觀切削。在載荷作用下，磨粒在材料表面發(fā)生類似刨削一樣的過程，這種過程會(huì)直接導(dǎo)致材料的去除。

b）由斷裂機(jī)制引起的去除過程

這種磨損機(jī)制主要針對(duì)脆性材料而言，比如陶瓷、玻璃等。磨粒在壓力作用下會(huì)與接觸表面形成極高的接觸應(yīng)力，以致于超過脆性材料的極限，從而發(fā)生脆性斷裂。

2 磨粒磨損的模擬仿真

根據(jù)磨粒磨損的概念可知，表面粗糙度對(duì)磨粒磨損的影響較大。由于表面形貌隨機(jī)性強(qiáng)、形成過程復(fù)雜、影響因素多，所以不可能通過實(shí)時(shí)監(jiān)測(cè)的手段來觀察表面形貌的變化，現(xiàn)在對(duì)表面形貌的研究多采用模擬仿真的手段［9］。趙永春［10］采用離散元原理直觀的對(duì)表面粗糙度進(jìn)行了模擬，結(jié)果表明在剛性光滑平面與理想粗糙表面相互作用后，粗糙度變小，表面結(jié)構(gòu)也發(fā)生很大變化，隨著兩個(gè)接觸表面法向接近量的增大，粗糙表面脫落的顆粒數(shù)增加，磨損加劇。Zhang［11-12］基于隨機(jī)過程理論提出了描述二體磨損過程中表面形貌變化的磨損動(dòng)態(tài)過程理論，通過假設(shè)剛性表面為絕對(duì)剛體，建立了軟表面與剛性表面摩擦?xí)r，受微觀切削而引起的表面粗糙度變化和磨損量的預(yù)測(cè)模型。劉峰壁［13］考慮到實(shí)際過程中硬表面也會(huì)受到磨損，于是在此模型基礎(chǔ)上進(jìn)行修正，建立了兩表面均受到磨損的預(yù)測(cè)模型。馮偉等［14］利用表面輪廓曲線高度近似服從高斯分布這一統(tǒng)計(jì)學(xué)特點(diǎn)，分別用matlab和ansys建立了表面形貌的三維仿真模型，并用赫茲接觸模型計(jì)算了不同載荷下接觸應(yīng)力在表面的分布情況。

吳國清等［15］運(yùn)用Monte Carlo方法隨機(jī)創(chuàng)建磨粒形狀，采用打靶法隨機(jī)生成磨粒表面，采用網(wǎng)格剖分法記憶表面形貌，在三維空間內(nèi)模擬了多個(gè)磨粒同時(shí)參與磨損的隨機(jī)動(dòng)態(tài)磨損過程，并通過實(shí)驗(yàn)驗(yàn)證了模型的有效性。岑巖宏［16］利用三維有限元法動(dòng)態(tài)模擬剛性球形磨粒在彈-塑性平板上的壓入-滑動(dòng)過程，模擬了不同壓入深度下表面的磨損劃傷行為，結(jié)合實(shí)驗(yàn)數(shù)據(jù)建立了適合碳鋼系列材料塑性變形脊形貌與磨粒壓入深度和寬度的曲線關(guān)系。Jacobson［17］建立了兩維磨損面內(nèi)僅限于純切削過程的統(tǒng)計(jì)學(xué)模型，在模型中考慮了多個(gè)磨粒同時(shí)作用的影響，從而使它更接近于實(shí)際情況。文中利用該模型研究了磨粒尺寸、載荷、工件表面硬度對(duì)磨損率的影響，預(yù)測(cè)了接觸的磨粒數(shù)量和磨損表面形貌，并與其他研究試驗(yàn)結(jié)果進(jìn)行對(duì)比，驗(yàn)證了模型的有效性。Podra［18］建立錐形接觸時(shí)考慮表面形貌的磨粒磨損模型，用有限元和解析方法對(duì)磨損進(jìn)行了計(jì)算，但兩種方法計(jì)算結(jié)果相差較大，分析認(rèn)為產(chǎn)生差異的原因是在解析方法中沒有考慮彈性變形的影響。M.Barge［19］等著重研究了磨粒磨損中的犁溝現(xiàn)象，建立了較硬的圓頭形磨粒在較軟表面的犁溝切削模型，考察了載荷、材料和磨粒的硬度等因素對(duì)表面應(yīng)力分布、殘余應(yīng)力和磨損量的影響，并與Bowden和Tabor的犁溝切削模型進(jìn)行了對(duì)比。Anders S?derberg［20］通 過 對(duì) 磨 損 過 程 的 簡(jiǎn) 化，基 于Archard磨損模型和顯式歐拉積分求解方法，建立了制動(dòng)片在磨合期時(shí)的三維有限元模型，分析了磨合過程中接觸應(yīng)力的分布和磨損變化，發(fā)現(xiàn)在磨合初期，表面粗糙度較大，應(yīng)力分布非常不均勻，隨著磨合的進(jìn)行，應(yīng)力分布趨于均勻，磨損程度和應(yīng)力分布有關(guān)，即接觸應(yīng)力大的區(qū)域磨損較為嚴(yán)重。

3 粘著磨損的模擬仿真

粘著磨損中較有代表性的是Holm模型和Archard模型［21］，相比Holm模型，Archard模型更接近于粘著磨損的真實(shí)情況，因此大部分粘著磨損都是在此模型基礎(chǔ)上加以修正而成。

Sung-San Cho［22］等從分子作用力的角度出發(fā)，建立了半球形彈性體與剛性表面的粘著接觸有限元模型，計(jì)算了載荷、接觸半徑、脫粘力等相關(guān)參數(shù)的關(guān)系，結(jié)果表明較小的表面粗糙度能有效減小靜摩擦力和粘著摩擦力，但是會(huì)加劇粘著磨損。X.Yin和K.Komvopoulos［23］建立了粗糙表面粘著效應(yīng)的粘著磨損模型，以表面微凸體是否發(fā)生完全塑性變形為磨損發(fā)生準(zhǔn)則，將接觸面積、磨損率、磨損系數(shù)等目標(biāo)參數(shù)表示為法向載荷、表面粗糙度、材料的彈塑性變形參數(shù)、表面能、材料相容度等因素的函數(shù)，并采用典型的陶瓷/陶瓷、陶瓷/金屬、金屬/金屬摩擦副的材料參數(shù)對(duì)模型進(jìn)行計(jì)算，得到了相關(guān)影響因素與目標(biāo)參數(shù)之間的曲線關(guān)系，該模型比較完整的包含了粘著磨損過程中的各影響因素，對(duì)磨損的建模過程具有一定的借鑒意義。

綜上所述，數(shù)值仿真方法作為新的磨損研究手段和方法己引起摩擦學(xué)界的重視，并將發(fā)揮越來越重要的作用。目前關(guān)于磨損仿真研究的主要不足在于磨損仿真模型預(yù)測(cè)結(jié)果同實(shí)際情況相差較大;建立仿真模型時(shí)的簡(jiǎn)化處理導(dǎo)致模型精確度不夠高;利用Arhcard等經(jīng)典磨損模型進(jìn)行仿真計(jì)算時(shí)缺乏關(guān)于材料磨損率規(guī)律的準(zhǔn)確數(shù)據(jù)［24］。磨損仿真結(jié)果的可靠性和實(shí)用性在相當(dāng)大程度上取決于磨損模型的正確性和精確性。為了更好地發(fā)揮數(shù)值仿真方法在摩擦磨損研究中的作用，應(yīng)在深化磨損機(jī)理研究、摩擦學(xué)系統(tǒng)分析方法及磨損表面微觀分析的基礎(chǔ)上建立更加完善和切合實(shí)際的磨損模型，并通過大量實(shí)驗(yàn)研究確定關(guān)于磨損率規(guī)律的準(zhǔn)確數(shù)據(jù)。

4 溫度－應(yīng)力場(chǎng)仿真

在載荷作用下，摩擦接觸區(qū)域會(huì)產(chǎn)生大量的摩擦熱，摩擦熱嚴(yán)重影響材料的物理力學(xué)性能和化學(xué)性能，對(duì)磨損形式的改變有著直接影響。例如制動(dòng)器在制動(dòng)過程中溫度場(chǎng)分布不均，這將引起兩接觸體的熱彈性變形，從而改變了接觸面之間的壓力分布。反過來，接觸界面之間的壓力分布不均勻，使制動(dòng)器溫度場(chǎng)分布更加不均勻，形成了一種惡性的循環(huán)，這種現(xiàn)象被稱作熱彈性不穩(wěn)定（TEI）［25］。對(duì)樹脂基制動(dòng)片/金屬摩擦副而言，摩擦表面溫度分布不均勻性會(huì)導(dǎo)致制動(dòng)片一系列物化性能的變化，如熱彈性不穩(wěn)定性、熱分解、熱膨脹等。由于摩擦過程的特殊性，摩擦表面的溫度和應(yīng)力大小、變化規(guī)律、分布情況等都不可能實(shí)時(shí)測(cè)量，這為進(jìn)一步的研究工作造成困難，模擬仿真技術(shù)使其成為可能。

4.1 溫度場(chǎng)

摩擦副表面溫度的計(jì)算模型普遍采用Block溫度模型，即半無限體表面承受單一集中熱源，并提出了閃溫的概念。接著Jaeger發(fā)展了這一理論，闡述了矩形形狀的移動(dòng)熱源作用在半無限體表面上的數(shù)學(xué)模型［26］。Ling［27］運(yùn)用表面形貌的隨機(jī)模型來估算滑動(dòng)表面的瞬時(shí)溫升，Wang S［28］則對(duì)具有分形特征的粗糙表面的滑動(dòng)摩擦局部溫升進(jìn)行了研究，提出了界面溫度分布的分形理論。Tian和Kennedy［29］的研究指出，實(shí)際的滑動(dòng)接觸都是多點(diǎn)離散接觸且滑動(dòng)副物體尺寸都是有限的，在接觸微凸體處除了有局部閃現(xiàn)溫升作用外，還有名義表面溫升的作用。這要求考慮物體的有限形狀，滑動(dòng)接觸的重復(fù)特性以及物體間總的“熱量平衡”。

由上述可知，制動(dòng)摩擦溫度場(chǎng)問題是典型的非線性問題，滑動(dòng)摩擦表面的溫升由名義表面溫升和局部表面溫升（閃現(xiàn)溫度）組成。求解名義表面溫升一般都假設(shè)熱流分配系數(shù)為一常數(shù)，然后按照理論分析（或某個(gè)實(shí)驗(yàn)數(shù)據(jù)）得出一個(gè)確定的熱流分配系數(shù)。實(shí)際上，制動(dòng)過程中摩擦表面上產(chǎn)生的熱量在零件間的分配是變化的，與物體的尺寸大小，材料的熱物理性能及制動(dòng)時(shí)間有關(guān)，并隨材料的導(dǎo)熱性與接觸表面處的溫度梯度之比值而變化［30］。

在制動(dòng)摩擦名義表面溫升計(jì)算的熱模型中，名義接觸表面摩擦熱流邊界條件的確定是建立和求解的關(guān)鍵。K.Fridrich［31］建立了單個(gè)鋼微凸體在PEEK或碳纖維-PEEK復(fù)合材料上滑動(dòng)的模型，該模型假設(shè)單位面積上所產(chǎn)生的摩擦熱取決于接觸應(yīng)力p、摩擦系數(shù)μ和滑動(dòng)速度v，然后通過計(jì)算接觸點(diǎn)處的單位熱源強(qiáng)度來計(jì)算接觸溫度。Kennedy和Ling［32］最早運(yùn)用數(shù)值模擬的方法研究盤式制動(dòng)器的熱彈性不穩(wěn)定現(xiàn)象，但他的模型簡(jiǎn)化為軸對(duì)稱模型，忽略了移動(dòng)熱源的問題。

4.2 應(yīng)力場(chǎng)

摩擦過程中溫度場(chǎng)和應(yīng)力場(chǎng)是相互耦合、彼此影響的。M.Tirovic和A.J.Day［33］在詳細(xì)地研究了盤式制動(dòng)器的壓力分布情況后，認(rèn)為摩擦制動(dòng)器的摩擦襯片與對(duì)偶件表面之間產(chǎn)生的摩擦熱并不是均勻地分布在滑動(dòng)表面上的，而是取決于局部應(yīng)力的分布情況;而由于溫度效應(yīng)造成熱膨脹和表面熱裂紋等結(jié)果又會(huì)反過來影響應(yīng)力分布。唐旭晟在求解溫度-應(yīng)力耦合場(chǎng)的做法是:首先假設(shè)初始階段應(yīng)力是均勻分布的，求出溫度場(chǎng)的解，然后將溫度作為施加載荷，反過來求解溫度引起的熱應(yīng)力分布。Hasan Sofuoglu等［34］分別求解了考慮溫度和不考慮溫度影響下的表面應(yīng)力分布，發(fā)現(xiàn)在溫度影響下表面應(yīng)力更大，而較大的表面應(yīng)力主要分布在第一粗糙峰上，這是因?yàn)樵诘谝淮畏ㄏ嘟佑|時(shí)，主要是這些較“突出”的微凸體承受載荷，當(dāng)施加位移時(shí)，摩擦熱在這些微凸體上造成溫度累積效應(yīng)，從而產(chǎn)生熱彈塑性變形，所以表面應(yīng)力要更大一些。Gong和Komvopoulos［35］也分別計(jì)算了只考慮彈性接觸/彈塑性接觸和溫度影響下的熱彈塑性接觸的應(yīng)力大小，結(jié)果與 Hasan Sofuoglu 一致。Ji-Hoon Choi［36］建立了三維有限元盤式制動(dòng)器模型，模擬了反復(fù)制動(dòng)情況下盤/片的熱彈性接觸問題，深入分析了溫度、應(yīng)力及材料屬性之間的內(nèi)在關(guān)系。發(fā)現(xiàn)在初始階段應(yīng)力分布較均勻，隨著制動(dòng)的進(jìn)行應(yīng)力沿徑向方向發(fā)生明顯變化，局部出現(xiàn)應(yīng)力集中，應(yīng)力的變化主要是由于溫度升高導(dǎo)致表面變形而造成的，應(yīng)力分布和大小隨著溫度變化而變化，即所謂的熱應(yīng)力不穩(wěn)定性（TEI）。在研究材料屬性與應(yīng)力分布關(guān)系時(shí)，發(fā)現(xiàn)熱膨脹系數(shù)和彈性模量對(duì)應(yīng)力分布有很大影響，在仿真中將彈性模量減半后，應(yīng)力分布較之前更加均勻，這說明較軟的材料有利于增加實(shí)際接觸面積，使應(yīng)力分布更加均勻。相比各向同性材料，各向異性材料的熱彈性穩(wěn)定性更好，溫度和應(yīng)力的分布也更加均勻。

5 結(jié)語

如前所述，由于模擬仿真技術(shù)更易于變量的控制，與試驗(yàn)研究相比，具有可操作性、可重復(fù)性、可比性等優(yōu)點(diǎn)，因而越來越受到摩擦磨損領(lǐng)域的研究學(xué)者們的重視。但是也存在一些問題，如模擬結(jié)果與實(shí)驗(yàn)數(shù)據(jù)不吻合，建模時(shí)的簡(jiǎn)化處理影響模型精度，個(gè)別物理量難以通過實(shí)驗(yàn)得到驗(yàn)證等。為了更好的使用模擬仿真技術(shù)進(jìn)行研究，筆者覺得應(yīng)該從以下方面入手:首先，深化磨損機(jī)理的研究。只有真正弄清楚磨損機(jī)理，建立完善的磨損理論，才能從根本上指導(dǎo)模擬仿真的建模過程，才能保證模型的科學(xué)性和精確性。其次，盡量在模型中考慮更多的參量，過度的簡(jiǎn)化會(huì)使仿真精度下降;再次，加強(qiáng)實(shí)驗(yàn)驗(yàn)證的環(huán)節(jié)，針對(duì)所建立的模型設(shè)計(jì)專門的實(shí)驗(yàn)，只有通過實(shí)驗(yàn)驗(yàn)證才能最終保證模型的正確性。

［1］溫詩鑄.摩擦學(xué)原理［M］.北京:清華大學(xué)出版社，1990.

［2］黃平.摩擦學(xué)教程［M］.北京:高等教育出版社，2008.

［3］江親瑜.零件磨損過程及壽命預(yù)測(cè)的數(shù)值仿真［J］.潤(rùn)滑與密封，1997，29（6）:29-30.

［4］Yan L，XuJ J.Research on the simulation method of wear process［J］.Tribology，1999，19（1）:50-55.

［5］趙春華.基于摩擦學(xué)與動(dòng)力學(xué)的摩擦學(xué)系統(tǒng)狀態(tài)描述方法研究［D］.武漢，武漢理工大學(xué)，2005.

［6］B.布尚.摩擦學(xué)導(dǎo)論［M］.北京:機(jī)械工業(yè)出版社，2006.

［7］蘇立志，李志章.耐溫聚合物基復(fù)合材料的摩擦學(xué)性能與磨損機(jī)理研究進(jìn)展［J］.材料科學(xué)與工程，2000（70）:125-130.

［8］劉家浚.材料磨損機(jī)理及其耐磨性［M］.北京:清華大學(xué)出版社，1993.

［9］閆小青，謝志龍，樊保圣，等.基于表面粗糙度影響的摩擦聲發(fā)射特性研究［J］.Tribology，2011（31）:587-591.

［10］趙永春，劉小君，王偉，等.基于離散元法的表面形貌相互作用研究［J］.合肥工業(yè)大學(xué)學(xué)報(bào)，2010（33）:487-490.

［11］Zhang B，Xie Y B.Two body micro-cutting wear model，Part I:Two-dimensional roughness model［J］.Wea，1989，129:37-48.

［12］Zhang B，Xie Y B.Two body micro-cutting wear model，Part II:Three-dimensional roughness model［J］.Wear，1989，129:49-58.

［13］Liu F B，Li X E，Xie Y B.A new model for predicting roughness of rubbing surface in three bodyabrasive wear［J］.Chinese Journal of Mechanical Engineering，2000，19（1）:9-12.

［14］馮偉，嚴(yán)新平，周新聰.銷盤滑動(dòng)磨損試驗(yàn)的仿真建模研究［J］.中國機(jī)械工程，2005（16）:2141-2144.

［15］吳國清，張曉峰，方亮，等.兩體磨料磨損的三維動(dòng)態(tài)模擬［J］.Tribology，2000（20）:360-364.

［16］岑啟宏，孫錕，方亮，等.二體磨料磨損犁溝及脊的三維有限元?jiǎng)討B(tài)模擬［J］.Tribology，2004（24）:249-253.

［17］Jacobson S，Wallen P，Hogmark S.Fundamental aspects of abrasive wear studied by a new numerical simulation model［J］.Wear，2008，123:207-223.

［18］ Podra P，Andersson S.Finite element analysis of wear simulation of a conical spinning contact considering surface topography［J］.Wear，1999，224:13-21.

［19］M.Barge，G.Kermouche，P.Gilles，J.M.Bergheau.Experimental and numerical study of the ploughing part of abrasive wear［J］.Wear，2003（255）:30-37.

［20］Anders Sderberg，S ren Andersson.Simulation of wear and contact pressure distribution at the pad-to-rotor interface in a disc brake using general purpose nite element analysis software［J］.Wear，2009（267）:2243-2251.

［21］Valentin L.Popov.Contact Mechansics and Friction physical Principles and Applications［M］.北京:清華大學(xué)出版社，2011.

［22］Sung-San Cho，Seungho Park.Finite element modeling of adhesive contact using molecular potential［J］.Tribology，2004（37）:763-769.

［23］X.Yin，K.Komvopoulos.An adhesive wear model of fractal surfaces in normal contact［J］.Solid and Structure，2010（47）:912-921.

［24］汪選國，嚴(yán)新平，李濤生，等.磨損數(shù)值仿真技術(shù)的研究進(jìn)展［J］.摩擦學(xué)學(xué)報(bào)，2004（24）:188-192.

［25］高誠輝，黃健萌，林謝昭，等.盤式制動(dòng)器摩擦磨損熱動(dòng)力學(xué)研究進(jìn)展［J］.中國工程機(jī)械學(xué)報(bào)，2006（4）:83-88.

［26］ O.O.Evtushenko，E.H.Ivanyk，N.V.Horbachova.Analytic methods for thermal calculation of brakes（review）［J］.Materials Science，2000，36（6）:857-862.

［27］F.Ling.On temperature transients at sliding interface［J］.ASME Journal of Lubrication Technology，1969，7:397-405.

［28］ Wang S.A fractal theory of the interfacial temperature distribution in the slow sliding Regime［J］.ASME J of Trio.，1994，116-812.

［29］Xuefeng Tian，F(xiàn)rancis E Kennedy Jr.Contact surface temperature models for finite bodies in dry and boundary lubricated sliding［J］.ASME Journal of Tribology，1993，115:411-418.

［30］馬保吉，朱均.摩擦制動(dòng)器接觸表面溫度計(jì)算模型［J］.西安工業(yè)學(xué)院學(xué)報(bào)，1999，19（1）:35-39.

［31］K.Friedrich，J.Flock，K.Varadi and Z.Neder.Numerical and finite element contact and thermal analysis of real composite-steel surfaces in sliding contact［J］.Wear，1999，225-229:368-379.

［32］Kennedy，F(xiàn).E.，Ling，F(xiàn).F.A thermal，Thermoelastic，and Wear Simulation of a High-Emergy SlidingContact Problem［J］.ASME J.Lubr.Technol，1974，l97，497-508.

［33］J.Halling.Principles of Tribology［M］.北京:機(jī)械工業(yè)出版社，1975.

［34］Hasan Sofuoglu，Alaettin Ozer.Thermomechanical analysis of elastoplastic medium in sliding contact with fractal surface［M］.Tribology，2008（41）:783-796.

［35］ Gong Z-Q，Komvopoulos K.Thermo mechanical analysis of semi-infinite solid in sliding contact with a fractal surface.ASME J Tribol，2005;127:331 – 42.

［36］Ji-Hoon Choi，In Lee.Finite element analysis of transient thermoelastic behaviors in disk brakes［J］.Wear，2004（257）:47-58.