兩物體間的摩擦生熱分析

1、兩物體相對轉(zhuǎn)動過程中的摩擦生熱分析摘要：有限元分析方法是隨計算機(jī)的發(fā)展而迅速發(fā)展起來的在計算數(shù)學(xué)、計算力學(xué)和計算工程科學(xué)領(lǐng)域的先進(jìn)計算方法。其中的熱分析可以識別出系統(tǒng)或部件的溫度分布及其他熱物理參數(shù)，為系統(tǒng)的結(jié)構(gòu)分析以及結(jié)構(gòu)特性的優(yōu)化設(shè)計提供依據(jù)。文章利用有限元軟件對兩物體相對轉(zhuǎn)動過程中的摩擦生熱進(jìn)行了有限元熱-結(jié)構(gòu)耦合分析，計算了模型的溫度場以及熱應(yīng)力場，并通過對結(jié)果的分析對零件結(jié)構(gòu)進(jìn)行了說明，為該結(jié)構(gòu)的實際應(yīng)用提供了設(shè)計依據(jù)。關(guān)鍵字：有限元 ANSYS 熱結(jié)構(gòu)耦合分析1. 有限元思想及ANSYS簡介有限元法是隨著電子計算機(jī)的發(fā)展而迅速發(fā)展起來的一種現(xiàn)代計算方法，是20世紀(jì)50年代首先在連續(xù)

2、力學(xué)領(lǐng)域飛機(jī)結(jié)構(gòu)靜、動態(tài)特性分析中應(yīng)用的一種有效的數(shù)值分析方法，隨后很快就廣泛地用于求解熱傳導(dǎo)、電磁場、流體力學(xué)等連續(xù)性問題。1.1. 有限元法簡介1.1.1. 有限元法的基本思想有限元分析計算的思路和作法可歸納如下：1. 物理離散化將某個工程結(jié)構(gòu)離散為由各種連接單元組成的計算模型，這一步稱作單元剖分。離散后單元與單元之間利用單元的節(jié)點相互連接起來。單元節(jié)點的設(shè)置、性質(zhì)、數(shù)目等應(yīng)視問題的性質(zhì)，描述變形形態(tài)要根據(jù)需要和計算精度而定。所以有限元法中分析的結(jié)構(gòu)已不是原有的物體或結(jié)構(gòu)物，而是同樣的材料由眾多單元以一定方式連接成的離散物體。一般情況，單元劃分越細(xì)則描述變形情況越精確，即越接近實際變形，但

3、計算量越大。2. 單元特性分析（1）選擇未知量模式在有限元法中，選擇節(jié)點位移作為基本未知量時稱為位移法；選擇節(jié)點力作為基本未知量時稱為力法；取一部分節(jié)點力和一部分節(jié)點位移作為基本未知量時稱為混合法。位移法易于實現(xiàn)計算自動化，所以在有限單元法中位移法應(yīng)用范圍最廣。（2）分析單元的力學(xué)性質(zhì)根據(jù)單元的材料、形狀、尺寸、節(jié)點數(shù)目、位置及其含義等，找出單元節(jié)點力和節(jié)點位移的關(guān)系式，這是單元分析中的關(guān)鍵一步，也是有限元法的基本步驟之一。（3）計算等效節(jié)點力對于實際的連續(xù)體，力是從單元的公共邊界傳遞到另一個單元中去的。因而，這種作用在單元邊界上的表面力、體積力或集中力都需要等效的移到節(jié)點上去，也就是用等效的

4、節(jié)點力來替代所有作用在單元上的力。3. 單元組集利用結(jié)構(gòu)力的平衡條件和邊界條件把各個單元按照原來的結(jié)構(gòu)重新連接起來，形成整體的有限元方程。4. 求解未知節(jié)點位移求解有限元方程式（上式）得出位移?？梢钥闯觯邢拊ǖ幕舅枷胧恰耙环忠缓稀保质菫榱诉M(jìn)行單元分析，合則是為了對整體結(jié)構(gòu)進(jìn)行綜合分析。1.1.2. 有限元法的基本要素構(gòu)成有限元系統(tǒng)的3個基本要素是節(jié)點、單元和自由度：節(jié)點（Node）：節(jié)點是構(gòu)成有限元系統(tǒng)的基本對象，也就是整個工程系統(tǒng)中的最基本點。它包含了坐標(biāo)位置以及具有物理意義的自由度信息。單元（Element）：單元是由節(jié)點與節(jié)點相連而成，是構(gòu)成有限元系統(tǒng)的基礎(chǔ)。一個有限元系統(tǒng)中必須

5、有至少一個以上的單元。單元和單元之間由各個節(jié)點相互連接。自由度（DOF，Degree Of Freedom）：包括系統(tǒng)自由度和節(jié)點自由度。整個系統(tǒng)的自由度，在分析中需要進(jìn)行適當(dāng)?shù)募s束，系統(tǒng)的每個節(jié)點都有各自的節(jié)點坐標(biāo)系和對應(yīng)的節(jié)點自由度，對于不同的單元上的節(jié)點，具有不同的自由度。1.1.3. 有限元法的分析步驟有限元分析是物理現(xiàn)象（幾何及載荷工況）的模擬，是對真實情況的數(shù)值近似。通過分析對象劃分網(wǎng)絡(luò)，求解有限個數(shù)值來近似模擬真實環(huán)境的無限個未知量。ANSYS分析過程中包含三個主要步驟。1. 創(chuàng)建有限元模型（1）創(chuàng)建或讀入幾何模型（2）定義材料屬性（3）劃分網(wǎng)格（節(jié)點及單元）2. 施加載荷并求解

6、（1）施加載荷及載荷選項、設(shè)定約束條件（2）求解3. 查看結(jié)果（1）查看分析結(jié)果（2）檢驗結(jié)果（分析是否正確）1.2. ANSYS軟件簡介1.2.1. ANSYS軟件的發(fā)展歷史ANSYS是一種融結(jié)構(gòu)、熱、流體、電磁和聲學(xué)于一體的大型CAE通用有限元分析軟件，可廣泛用于核工業(yè)、鐵道、石油化工、航空航天、機(jī)械制造、能源、汽車交通、國防軍工、電子、土木工程、造船、生物醫(yī)學(xué)、輕工、地礦、水利，以及日用家電等一般工業(yè)及科學(xué)研究。該軟件可在大多數(shù)計算機(jī)及操作系統(tǒng)中運行。從PC到工作站，直至巨型計算機(jī)，ANSYS文件在其所有的產(chǎn)品系列和工作平臺上均兼容。ANSYS是第一個集成計算機(jī)流體動力學(xué)功能的軟件，也是

7、唯一一個包括多物理場分析功能軟件。ANSYS是Analysis SYStem的縮寫，是一種廣泛性的商業(yè)套裝工程分析軟件。它由世界上著名的有限元分析軟件公司ANSYS開發(fā)，它能與大多數(shù)CAD軟件結(jié)合使用，實現(xiàn)數(shù)據(jù)共享和交換，如AutoCAD、I-DEAS、Pro/Engineer、NASTRAN、Alogor等，是實現(xiàn)現(xiàn)代產(chǎn)品設(shè)計中的高級CAD工具之一。該軟件從1971年的2.0版本至現(xiàn)在的12.0版本，已有近40年的歷史。目前已有許多國際化大公司以ANSYS作為其標(biāo)準(zhǔn)。1.2.2. ANSYS軟件的基本功能ANSYS的基本功能有：結(jié)構(gòu)靜力分析、結(jié)構(gòu)動力學(xué)分析、結(jié)構(gòu)非線性分析、動力學(xué)分析、熱分析

8、、電磁場分析、計算流體動力學(xué)分析、聲場分析、壓電分析等，高級功能有多物理場耦合分析、優(yōu)化設(shè)計、拓?fù)鋬?yōu)化等。熱分析用于計算一個系統(tǒng)或部件的溫度分布及其它熱物理參數(shù)，如熱量的獲取或損失、熱梯度、熱流密度（熱通量）等。熱分析在許多工程應(yīng)用中扮演著重要角色，如內(nèi)燃機(jī)、渦輪機(jī)、換熱器、管路系統(tǒng)、電子元件、鍛造、鑄造等。ANSYS熱分析基于能量守恒原理的熱平衡方程，用有限元法計算各節(jié)點的溫度，并導(dǎo)出其它熱物理參數(shù)。ANSYS熱分析包括熱傳導(dǎo)、熱對流及熱輻射三種熱傳遞方式。此外，還可以分析相變、有內(nèi)熱源、接觸熱阻等問題。ANSYS熱分析主要分為兩大類，一是穩(wěn)態(tài)傳熱，是指系統(tǒng)的溫度場不隨時間變化；二是瞬態(tài)傳熱

9、，是指系統(tǒng)的溫度場隨時間明顯變化。ANSYS中與熱相關(guān)的耦合場分析主要有熱結(jié)構(gòu)耦合、熱流體耦合、熱電耦合、熱磁耦合以及熱電磁結(jié)構(gòu)耦合等。2. 熱結(jié)構(gòu)耦合分析的有限元法2.1. 熱分析基本知識2.1.1. 熱傳遞的方式如上文所述，熱傳遞的方式主要有熱傳導(dǎo)、熱對流和熱輻射三種方式。在絕大多數(shù)情況下，我們分析的熱傳導(dǎo)問題都帶有對流和/或輻射邊界條件。1、熱傳導(dǎo)熱傳導(dǎo)可以定義為完全接觸的兩個物體之間或一個物體的不同部分之間由于溫度梯度而引起的內(nèi)能的交換。熱傳導(dǎo)遵循傅里葉定律：，其中，為熱流密度（W/m2），Km為導(dǎo)熱系數(shù)（W/m·），為沿向德溫度帝都，負(fù)號表示熱量流向溫度降低的方向。2、熱對

10、流熱對流是指固體的表面和與它周圍接觸的流體之間，由于溫差的存在引起的熱量的交換。熱對流可以分為兩類：自然對流和強(qiáng)制對流。對流一般作為面邊界條件施加。熱對流用牛頓冷卻方程來描述：，其中，為對流換熱系數(shù)（或稱膜傳熱系數(shù)、給熱系數(shù)、膜系數(shù)等），為固定表面的溫度，為周圍流體的溫度。3、熱輻射熱輻射指物體發(fā)射電磁能，并被其它物體吸收轉(zhuǎn)變?yōu)闊岬臒崃拷粨Q過程。物體溫度越高，單位時間輻射的熱量越多。熱傳導(dǎo)和熱對流都需要有傳熱介質(zhì)，而熱輻射無須任何介質(zhì)。實質(zhì)上，在真空中的熱輻射效率最高。在工程中通?？紤]兩個或兩個以上物體之間的輻射，系統(tǒng)中每個物體同時輻射并吸收熱量。它們之間的凈熱量傳遞可以用斯蒂芬波爾茲曼方程來

11、計算：，其中，Q為熱流率，為吸射率（黑度），為斯蒂芬-波爾茲曼常數(shù)，約為5.67 BTU/h·in2·K4，A1為輻射面1的面積，F(xiàn)12為由輻射面2的形狀系數(shù)，T1為輻射面1的絕對溫度，T2為輻射面2的絕對溫度。由上式可以看出，包含熱輻射的熱分析是高度非線性的。在ANSYS中將輻射按平面現(xiàn)象處理（體都假設(shè)為不透明的）。2.1.2. 熱力學(xué)第一定律熱力學(xué)第一定律是熱分析的理論依據(jù)，也成能量守恒定律，即對于一個封閉的系統(tǒng)（沒有質(zhì)量地流入和流出）：，其中，Q為熱量；W為作功；為系統(tǒng)內(nèi)能；為系統(tǒng)動能；為系統(tǒng)勢能。對于大多數(shù)工程傳熱問題：。通?？紤]沒有做功：W=0，則：；對于穩(wěn)態(tài)熱分

12、析：=0，即流入系統(tǒng)的熱量等于流出的熱量；對于瞬態(tài)熱分析：，即流入或流出的熱傳遞速率q等于系統(tǒng)的內(nèi)能的變化。將其應(yīng)用到一個微元體上，就可以得到熱傳導(dǎo)的控制微分方程。2.1.3. 熱分析的控制方程熱傳導(dǎo)的控制微分方程為：，其中，其中，為媒介傳導(dǎo)速率。2.2. 熱分析的有限元法熱分析一般可以分為穩(wěn)態(tài)熱分析、瞬態(tài)熱分析與非線性熱分析、熱輻射分析、相變分析、CFD分析以及與溫度場有關(guān)的耦合場分析。如果系統(tǒng)的凈熱流率為0，即流入系統(tǒng)的熱量加上系統(tǒng)自身產(chǎn)生的熱量等于流出系統(tǒng)的熱量：q流入+q生成-q流出=0，則系統(tǒng)處于熱穩(wěn)態(tài)。在穩(wěn)態(tài)熱分析中任一節(jié)點的溫度不隨時間變化。瞬態(tài)傳熱過程是指一個系統(tǒng)的加熱或冷卻過

13、程。在這個過程中系統(tǒng)的溫度、熱流率、熱邊界條件以及系統(tǒng)內(nèi)能隨時間都有明顯變化。ANSYS熱分析的邊界條件或初始條件可分為七種：溫度、熱流率、熱流密度、對流、輻射、絕熱、生熱。熱分析涉及到的單元有大約40種，其中純粹用于熱分析的有14種，它們?nèi)绫?所示。表2 熱分析單元列表單元類型名稱 說明 線性 LINK32 LINK33 LINK34 LINK31 兩維二節(jié)點熱傳導(dǎo)單元三維二節(jié)點熱傳導(dǎo)單元二節(jié)點熱對流單元二節(jié)點熱輻射單元二維實體PLANE55 PLANE77 PLANE35 PLANE75 PLANE78 四節(jié)點四邊形單元八節(jié)點四邊形單元三節(jié)點三角形單元四節(jié)點軸對稱單元八節(jié)點軸對稱單元三維實

14、體SOLID87 SOLID70 SOLID90 六節(jié)點四面體單元八節(jié)點六面體單元二十節(jié)點六面體單元殼 SHELL57 四節(jié)點四邊形殼單元點 MASS71 節(jié)點質(zhì)量單元2.3. 熱結(jié)構(gòu)耦合分析的有限元法熱-結(jié)構(gòu)耦合場分析是指在有限元分析的過程中考慮了溫度和應(yīng)力兩種物理場的交叉作用和相互影響。耦合場分析主要有兩種方法：序貫耦合方法和直接耦合方法。熱結(jié)構(gòu)耦合問題是結(jié)構(gòu)分析中通常遇到的一類耦合分析問題。由于結(jié)構(gòu)溫度場的分布不均會引起結(jié)構(gòu)的熱應(yīng)力，或者結(jié)構(gòu)部件在高溫環(huán)境中工作，材料受到溫度的影響會發(fā)生性能的改變，這些都是進(jìn)行結(jié)構(gòu)分析時需要考慮的因素。為此需要先進(jìn)行相應(yīng)的熱分析，然后再進(jìn)行結(jié)構(gòu)分析。因此

15、，熱結(jié)構(gòu)耦合分析是指求解溫度場對結(jié)構(gòu)中應(yīng)力、應(yīng)變和位移等物理量影響的分析類型。對于熱結(jié)構(gòu)耦合分析，在ANSYS中通常采用順序耦合分析方法，即先進(jìn)行熱分析求得結(jié)構(gòu)的溫度場，然后再進(jìn)行結(jié)構(gòu)分析，且將前面得到的溫度場作為體載荷加到結(jié)構(gòu)中，求解結(jié)構(gòu)的應(yīng)力分布。1. 兩物體相對轉(zhuǎn)動過程中的摩擦生熱分析實例在ANSYS中，兩物體由于摩擦產(chǎn)生的總熱流率由下式計算：，其中，F(xiàn)HTG為摩擦生熱的能力轉(zhuǎn)化因子（默認(rèn)為1）；為等效摩擦應(yīng)力；為兩物體的相對滑動速率。接觸面的熱流率為：，其中，為接觸面所得到的熱流率；FWGT為目標(biāo)面和接觸面熱量分配權(quán)因子（默認(rèn)值為0.5）。目標(biāo)面的熱流率：，其中為目標(biāo)面所得到的熱流率。

16、1.1. 問題描述及分析一銅塊在鋼環(huán)上滑動，鋼環(huán)固定，其材料參數(shù)見表1，鋼塊和銅塊間的摩擦因素為0.2，滑塊的角速度分為三組，分別為0.000666rad/s、0.00333rad/s、0.00666rad/s，計算時間為10s，計算鋼塊和銅塊由于摩擦產(chǎn)生的溫度場，以及鋼塊和銅塊的應(yīng)力分布，初始溫度為20，分析時，溫度采用，其他單位采用國際單位制。表1 銅塊和鋼塊的材料參數(shù)表材料溫度/彈性模量/GPa密度/(kg/m3傳熱系數(shù)/(W/m·比熱/(J/kg·熱膨脹系數(shù)/(1/泊松比銅201038900383390鋼202067800460該實例屬于熱結(jié)構(gòu)耦合場分析，屬于旋轉(zhuǎn)摩

17、擦生熱問題，選用耦合場三維六面體SOLID5八節(jié)點單元進(jìn)行分析，將角速度載荷轉(zhuǎn)化為切向位移載荷施加在銅滑塊上。1.2. 建立幾何模型ANSYS軟件的幾何建模主要有兩種形式，一是利用軟件本身的建模功能建立需要的模型；二是導(dǎo)入通過其他三維CAD軟件建好的幾何模型。由于ANSYS本身的建模功能相對較弱，其又與主流三維CAD軟件（如CATIA、UG、PRO/E等）擁有已經(jīng)非常成熟的接口，因此一般均采用導(dǎo)入模型的方法。由于本例中的兩物體結(jié)構(gòu)相對簡單，因此采用了軟件本身的建模方法，幾何模型如圖1所示。圖1 幾何模型1.3. 創(chuàng)建有限元模型及網(wǎng)格劃分在對問題進(jìn)行有限元分析時，首先要做的就是針對問題建立適當(dāng)?shù)?/p>

18、有限元模型，模型要與結(jié)構(gòu)系統(tǒng)的集合外形基本一致。建立有限元模型有兩種方法，即直接法和間接法。直接法是直接根據(jù)結(jié)構(gòu)的幾何外形建立節(jié)點和元素，不需要再進(jìn)行網(wǎng)格劃分，適用于比較簡單的結(jié)構(gòu)系統(tǒng)。而間接法是通過點、線、面、體積，先建立有限元幾何模型，再進(jìn)行實體網(wǎng)格劃分，以完成有限元模型的建立，適用于較復(fù)雜的結(jié)構(gòu)。有限元模型包括節(jié)點、單元、材料屬性、實常數(shù)、邊界條件，以及其他用來表達(dá)這個物理系統(tǒng)的特征。本實例采用的是間接法，選用了耦合場三維六面體SOLID5八節(jié)點單元進(jìn)行分析，建立適當(dāng)?shù)木W(wǎng)格密度之后，進(jìn)行了網(wǎng)格劃分，模型如圖2所示。并選取鋼環(huán)內(nèi)圈表面和銅滑塊外側(cè)表面為兩相對轉(zhuǎn)動物體的接觸面，建立如圖3所示

19、的接觸。圖2 有限元模型劃分網(wǎng)格 圖3 有限元模型建立接觸1.4. 加載求解根據(jù)實際情況給模型施加適當(dāng)?shù)膲毫d荷、溫度載荷、約束條件以及位移載荷。對于位移載荷，將轉(zhuǎn)動在微小位移上看成是一個移動，因此根據(jù)不同的轉(zhuǎn)速將此位移轉(zhuǎn)化為0.002、0.01、002三種位移分別進(jìn)行加載計算。其中熱-結(jié)構(gòu)耦合分析時溫度邊界條件的加載過程是將前面穩(wěn)態(tài)溫度場得到的溫度結(jié)果文件（后綴名為.Rth）讀入熱-結(jié)構(gòu)耦合分析的網(wǎng)格模型中，模型中的每個節(jié)點都加載上了各自的溫度值，溫差的存在使節(jié)點之間產(chǎn)生熱應(yīng)力，熱應(yīng)力再和機(jī)械應(yīng)力進(jìn)行耦合，最后得到模型中節(jié)點的綜合應(yīng)力。三種轉(zhuǎn)速情況下的溫度場分布云圖及其對應(yīng)的應(yīng)力分布云圖如

20、圖4、圖5、圖6、圖7、圖8、圖9所示。1.5. 結(jié)果分析通過ANSYS熱結(jié)構(gòu)耦合分析的方法，得出了兩物體相對轉(zhuǎn)動時摩擦生熱后的溫度場分布情況以及相應(yīng)的應(yīng)力分布情況，并對不同轉(zhuǎn)速時的溫度場及應(yīng)力分布情況進(jìn)行了對比說明，為了避免產(chǎn)生因摩擦生熱影響機(jī)械結(jié)構(gòu)的可靠性，可采取如下措施：（1）選擇適當(dāng)?shù)霓D(zhuǎn)速。由上文圖表所示，轉(zhuǎn)速越高，其溫度場的最高溫度就越高，且高溫集中現(xiàn)象越明顯，相應(yīng)的應(yīng)力分布基值就越大，那么該結(jié)構(gòu)的可靠性就越差。（2）注意磨損。本模型的摩擦實質(zhì)上屬于滑動摩擦，可以認(rèn)為是一個滑動摩擦副。工程中，對于長時間處于高壓、高速運動狀態(tài)的相互接觸的滑動摩擦副來說，接觸區(qū)的溫度升高會導(dǎo)致粘著磨損現(xiàn)

21、象的發(fā)生，嚴(yán)重影響著摩擦副的正常運動。在摩擦過程中，因表層材料的變形或破裂而耗掉的能量大部分轉(zhuǎn)變成熱，從而引起表面溫度的升高。在載荷作用下，摩擦副表面間因有相對滑動而產(chǎn)生摩擦，摩擦生熱致使局部產(chǎn)生很高的溫升，有可能形成瞬時過熱，一個瞬時高溫的熱點會導(dǎo)致表面上相應(yīng)節(jié)點材料狀態(tài)的改變，造成摩擦表面的焊聯(lián)作用。隨后在分離的瞬時，連結(jié)點被撕裂。這樣，摩擦副表面將產(chǎn)生局部的初期粘著損傷，隨著多次損傷的積累，將引發(fā)明顯的粘著磨損。因此必須注意摩擦磨損現(xiàn)象。（3）加強(qiáng)關(guān)鍵部位結(jié)構(gòu)強(qiáng)度。由上文圖表所示，該模型的最高應(yīng)力集中在鋼環(huán)兩端支腳外端以及銅塊外側(cè)面，因此在進(jìn)行結(jié)構(gòu)設(shè)計時應(yīng)該加強(qiáng)這兩部位的結(jié)構(gòu)強(qiáng)度，以免發(fā)生失效。2. 總結(jié)本文對有限元思想進(jìn)行了簡要介紹，說明了ANSYS軟件在做熱分析，特別是熱結(jié)構(gòu)耦合分析方面的方法。通過實例的操作，熟悉了ANSYS軟件常用功能的操作使用，并對熱結(jié)構(gòu)耦合方面的知識有了較為深入的了解