朱茂桃,瞿恒軍,周啟豪
(江蘇大學(xué)汽車(chē)與交通工程學(xué)院,江蘇鎮(zhèn)江 212013)
濕式多片離合器依靠多片摩擦副摩擦傳遞力矩。多片和密閉浸油的設(shè)計(jì)特點(diǎn)讓離合器具有接合力矩大、工作可靠、磨損穩(wěn)定等優(yōu)點(diǎn),同時(shí)也導(dǎo)致在換擋期間離合器多副摩擦產(chǎn)生大量的摩擦熱。濕式離合器主要的失效方式是熱失效,由于熱流分配的不均勻和熱應(yīng)力的作用,摩擦副出現(xiàn)表面燒蝕,熱應(yīng)力的作用最終導(dǎo)致翹曲變形,摩擦材料徹底失效。因此,摩擦副溫度場(chǎng)研究是離合器優(yōu)化設(shè)計(jì)必要的參考手段。
國(guó)內(nèi)外學(xué)者對(duì)摩擦副溫度場(chǎng)的問(wèn)題展開(kāi)了大量研究。三維多片摩擦副熱機(jī)耦合分析計(jì)算量大、耗時(shí)長(zhǎng),目前研究主要將問(wèn)題簡(jiǎn)化為三維單副模型、二維軸對(duì)稱(chēng)模型和半無(wú)限厚模型等,通過(guò)有限元法、差分法等數(shù)值方法進(jìn)行分析[1-5]。張金樂(lè)[6-7]利用Abaqus軟件將多片摩擦副簡(jiǎn)化為單副,進(jìn)行了三維模型的熱機(jī)耦合仿真,分析了油壓和轉(zhuǎn)速差等因素對(duì)摩擦副的溫升影響。此模型盤(pán)片上下面均關(guān)于中面對(duì)稱(chēng),未考慮盤(pán)片上下表面的差異和多組摩擦副片與片之間的差異。何澤銀[8]利用ANSYS軟件進(jìn)行了船用濕式離合器的瞬態(tài)熱結(jié)構(gòu)耦合分析,考慮了花鍵間的摩擦力作用,認(rèn)為不同摩擦副間的正壓力有所不同,最后得出了接排過(guò)程中各片摩擦副的溫度分布以及接觸壓力的瞬態(tài)變化過(guò)程。朱茂桃等[9-11]建立了汽車(chē)動(dòng)力學(xué)模型計(jì)算3種起步條件下干式雙離合器摩擦熱流,并通過(guò)計(jì)算流體仿真(CFD)分析了壓盤(pán)和摩擦片的溫度。上述研究認(rèn)為約束載荷和油缸活塞的正壓力是均勻分布的,忽略了實(shí)際約束的影響。在摩擦副實(shí)際工作過(guò)程中,外力輸入通常是活塞盤(pán)與卡簧,且活塞盤(pán)通常采用壓指的形式,這就導(dǎo)致了摩擦副表面壓力的不均勻,這些約束會(huì)加重各組元件間的溫度差異。
圖1所示為濕式離合器摩擦系統(tǒng)的示意圖,該模型包括3對(duì)摩擦副,即3個(gè)摩擦片和3個(gè)對(duì)偶鋼片,摩擦片由芯片和襯片粘接而成,摩擦副一側(cè)是活塞盤(pán),另一側(cè)依次是調(diào)整鋼片和卡簧。對(duì)偶鋼片作為被動(dòng)件通過(guò)花鍵與從動(dòng)軸連接,摩擦片作為主動(dòng)件通過(guò)花鍵與驅(qū)動(dòng)軸連接,摩擦副周向固定,可以軸向移動(dòng)。當(dāng)摩擦副需要接合傳遞動(dòng)力時(shí),液壓力推動(dòng)活塞盤(pán),逐步壓緊摩擦組件做軸向運(yùn)動(dòng),外殼體上的卡簧限制軸向位移,間隙消除,主從動(dòng)件之間的轉(zhuǎn)速差消除為0,完成接合。期間摩擦片與鋼片之間滑摩做功轉(zhuǎn)化為熱量在元件之間進(jìn)行傳遞。油液帶走熱量。表1給出了摩擦元件的材料和結(jié)構(gòu)參數(shù)。
圖1 濕式多片離合器摩擦系統(tǒng)結(jié)構(gòu)示意圖
表1 各元件結(jié)構(gòu)參數(shù)
考慮到對(duì)稱(chēng)結(jié)構(gòu)的特殊性,將模型處理為二維軸對(duì)稱(chēng)結(jié)構(gòu),如圖1所示,壓力大小為P0的壓力均勻分布在活塞表面,當(dāng)調(diào)整鋼片受到摩擦片5傳遞過(guò)來(lái)的壓力而與卡簧接觸時(shí),摩擦系統(tǒng)又會(huì)受到卡簧的集中載荷軸向約束。給各元件依次編號(hào)為鋼片0、摩擦片1、鋼片2、摩擦片3、鋼片4、摩擦片5和調(diào)整鋼片6,鋼片0到調(diào)整鋼片6的接觸表面依次編號(hào)為S0到S7。
接觸壓力是決定熱量分布的重要因素,因此參照?qǐng)D1所示結(jié)構(gòu)在有限元軟件Abaqus中建立多片離合器二維壓力模型,如圖2所示。分析卡簧、活塞盤(pán)作用下,各個(gè)摩擦面的接觸壓力分布。切分摩擦片,賦予襯片和基片不同的材料屬性。使用penalty函數(shù)定義接觸。模型采用CAX4I,4節(jié)點(diǎn)軸對(duì)稱(chēng)雙線性四邊形,非協(xié)調(diào)模式單元。考慮到主從面剛度差異,設(shè)置摩擦片網(wǎng)格0.8 mm,對(duì)偶鋼片網(wǎng)格1 mm,為了結(jié)果更接近真實(shí),獲得更準(zhǔn)確的壓力分布,鋼片和摩擦襯片在劃分網(wǎng)格時(shí)分別采用2∶2和2∶3的偏置比。
圖2 濕式多片離合器壓力模型示意圖
在仿真過(guò)程中,將均布?jí)毫Φ闹礟0設(shè)置為0.04 MPa,約束卡簧外側(cè)的6個(gè)自由度,考慮到壓力施加過(guò)程中,鋼片的花鍵齒處會(huì)發(fā)生彎曲變形,為了消除這種變形給仿真帶來(lái)的影響,對(duì)花鍵齒接合中心點(diǎn)進(jìn)行約束,約束其轉(zhuǎn)動(dòng)自由度,僅允許齒面接合處軸向位移。同樣地,約束摩擦片和活塞中面內(nèi)徑點(diǎn)除軸向以外的自由度。
將仿真結(jié)果處理為如圖3所示的3維可視化圖形,提取徑向節(jié)點(diǎn),導(dǎo)出節(jié)點(diǎn)數(shù)據(jù),利用3次樣條插值擬合法,得到接觸面徑向壓力曲線,如圖4所示。由圖3、4可知,多片摩擦副各面的壓力有差異,并不等于活塞壓力,且因?yàn)榭ɑ傻募休d荷的影響,各摩擦表面的壓力并不是均勻分布的。從軸向來(lái)看,接觸壓力的分布不同,在內(nèi)徑一側(cè),越靠近卡簧的接觸表面壓力越??;在外徑一側(cè),結(jié)果相反,越靠近卡簧的摩擦表面壓力越大。從徑向來(lái)看,最外端是花鍵齒接觸面,該面上接觸壓力為0,越靠近外徑處壓力越大。從圖4可以直觀看出,越靠近卡簧的摩擦表面壓力分布越不均勻,S6摩擦表面出現(xiàn)最大徑向壓差,達(dá)到0.173 MPa,壓力的不同間接導(dǎo)致摩擦副的溫度場(chǎng)不同,因此,活塞和卡簧約束對(duì)摩擦副溫度的影響不容忽視。
圖3 接觸壓力仿真
圖4 接觸面徑向壓力曲線
柱坐標(biāo)系下,二維熱傳導(dǎo)的控制方程為[12]
式中:ρ為材料密度;C為材料比熱容;λ為材料熱導(dǎo)率;T為溫度;r為徑向坐標(biāo);z為軸向坐標(biāo)。
摩擦接觸表面滿足第二類(lèi)邊界條件,熱流密度計(jì)算方程為:
式中:μ為摩擦因數(shù);r為徑向位置;對(duì)于熱量轉(zhuǎn)化率η,取值1;P為不同徑向位置和時(shí)間的壓力;ω表示摩擦副間轉(zhuǎn)速差,取值150 r/min。
摩擦因數(shù)μ是關(guān)于相對(duì)滑摩速度的函數(shù)[13],
式中:μk、μs是動(dòng)靜摩擦因數(shù),分別取 0.1、0.2;δ是常數(shù),取0.007?;^(guò)程符合Amonton定律,摩擦因數(shù)恒定不變。
表2為摩擦副的熱物性參數(shù),忽略摩擦片溝槽的影響,接觸面視為理想平面,則接觸面溫度相等。進(jìn)入各物體表面的熱通量按材料熱阻確定的熱流分配系數(shù)分配:
式中:下標(biāo)p、d分別表示鋼片和摩擦片;q表示熱通量;λ、ρ、C分別表示導(dǎo)熱系數(shù)、密度和比熱容。
鋼片和摩擦片各自分配的熱通量為:
表2 摩擦副熱物性參數(shù)
潤(rùn)滑油為摩擦副提供冷卻降溫的作用,潤(rùn)滑油性能參數(shù)如表3所示。
表3 潤(rùn)滑油參數(shù)(40℃)
內(nèi)外端面的對(duì)流換熱系數(shù)為
式中:下標(biāo) o表示外端面;i表示內(nèi)端面;Reo,i為圓柱面雷諾數(shù);c、n為常數(shù),c取值0.193,n取值0.618[12]。
鋼片無(wú)槽表面上的對(duì)流換熱系數(shù)[12]為
式中:re為當(dāng)量摩擦半徑,Rer為圓環(huán)面雷諾數(shù)。
油槽內(nèi)對(duì)流換熱系數(shù)計(jì)算式為[12]:
式中:Re為槽內(nèi)流體雷諾數(shù);Sg為槽橫截面積;Cg為槽截面周長(zhǎng);ηf和ηw為表征油液與溫度相關(guān)的物理性質(zhì),忽略溫度對(duì)油液的影響,取ηf/ηw=1。
由于接觸面被簡(jiǎn)化為理想平面,因此在摩擦面上使用等效對(duì)流換熱系數(shù)代替油槽的散熱作用:
式中:SG為溝槽總面積;ST為摩擦副摩擦面積。
整個(gè)模型對(duì)應(yīng)的邊界方程為
式中:Toil為油溫,即環(huán)境溫度;T0為初始時(shí)刻各盤(pán)片的內(nèi)部溫度。
對(duì)于非摩擦接觸面,如S0上只存在熱傳導(dǎo),熱傳導(dǎo)系數(shù)為λsp。SP和S7面只存在對(duì)流換熱。摩擦襯片與摩擦芯片粘結(jié)構(gòu)成了摩擦片整體,芯片與襯片間只存在熱傳導(dǎo),由于材料性質(zhì)不同,取摩擦片粘結(jié)面當(dāng)量熱物理性參數(shù):
式中:Zp和Zd分別表示鋼片軸向和摩擦片軸向;Ce、ρe和λe分別表示襯片和芯片粘接面的當(dāng)量比熱容、當(dāng)量密度和當(dāng)量熱導(dǎo)率;δ表示微元,這里取δzp=δzd。
有限差分法比有限元法更靈活快捷,概念直觀且通用性強(qiáng)。為保證多片溫度模擬的穩(wěn)定性,采用完全隱式的方法對(duì)熱傳導(dǎo)控制方程和邊界條件進(jìn)行離散,利用Taylor級(jí)數(shù)展開(kāi)得到式(12)的差分形式,將合并同類(lèi)項(xiàng),將方程變形,就得到了較為簡(jiǎn)潔的方程。熱傳導(dǎo)控制方程的差分形式為[14]:
寫(xiě)出邊界差分方程:
式中:a為熱擴(kuò)散系數(shù),a=λ/ρC;i為 r方向節(jié)點(diǎn)的位置編號(hào),從左到右依次增大;j為z方向節(jié)點(diǎn)的位置編號(hào),從下到上依增大。
利用Matlab編程,進(jìn)行數(shù)值求解和圖像處理。
前面進(jìn)行了接觸壓力的仿真,為了對(duì)比驗(yàn)證仿真的準(zhǔn)確性,進(jìn)行試驗(yàn)臺(tái)架的壓力分析(圖5)。試驗(yàn)中電機(jī)和制動(dòng)電機(jī)分別與摩擦片的主從動(dòng)件相連,通過(guò)調(diào)節(jié)電機(jī)轉(zhuǎn)速,獲得摩擦副的相對(duì)轉(zhuǎn)速。液壓站為摩擦副提供試驗(yàn)所需的壓力和潤(rùn)滑油,通過(guò)控制油壓調(diào)節(jié)摩擦副的接合。數(shù)據(jù)采集系統(tǒng)則采集各種信號(hào),通過(guò)傳感器收集轉(zhuǎn)速信號(hào)、壓力信號(hào)、油壓信號(hào)和轉(zhuǎn)矩信號(hào)。為方便鋼片溫度信號(hào)的采集,在滑摩期間使用制動(dòng)器令離合器從動(dòng)端鋼片保持靜止,摩擦元件進(jìn)行定速滑摩。在從動(dòng)軸上測(cè)量摩擦轉(zhuǎn)矩,主動(dòng)軸上測(cè)量相對(duì)轉(zhuǎn)速。
為了直接測(cè)量摩擦副的滑摩溫度,試驗(yàn)前選擇2、4號(hào)鋼片,布置熱電偶傳感器,考慮到鋼片的結(jié)構(gòu)特點(diǎn),沿著鋼片周向挖3個(gè)徑向長(zhǎng)孔,3個(gè)徑向長(zhǎng)孔深度分別為5、20、35 mm,對(duì)應(yīng)標(biāo)號(hào)分別為 A、B、C。
試驗(yàn)?zāi)M車(chē)輛緩速行駛過(guò)程,離合器還未處于完全接合狀態(tài),忽略半接合摩擦副熱彈變形的作用,將表面接觸壓力作為產(chǎn)熱主要因素。離合器相對(duì)轉(zhuǎn)速差100 r/min,接合壓力幅值0.04 MPa。
圖5 離合器試驗(yàn)臺(tái)架示意圖
圖6為定速滑摩試驗(yàn)中接合壓力隨時(shí)間變化曲線。由圖6可知,試驗(yàn)時(shí)長(zhǎng)共22 s。0.04 MPa的油壓并非一下實(shí)現(xiàn),還是從3 s左右時(shí)開(kāi)始上升,大概到7 s左右時(shí)達(dá)到設(shè)定值,隨后一直保持在0.04 MPa,直到22 s試驗(yàn)結(jié)束。
圖6 定速滑摩試驗(yàn)中接合壓力隨時(shí)間變化的曲線
根據(jù)熱電偶傳感器采集的溫度信息,整理成鋼片徑向側(cè)點(diǎn)處溫度隨時(shí)間的變化曲線,如圖7所示。
圖7 對(duì)偶鋼片徑向測(cè)點(diǎn)溫度隨時(shí)間變化的曲線
由圖7發(fā)現(xiàn):A4溫度明顯高于A2,分別為101、94℃;B4溫度略低于B2,分別為77、79℃;A4溫度略低于A2,分別為52、54℃。
試驗(yàn)證明,摩擦副溫度與軸向分布有關(guān),不同軸向位置處的摩擦件溫度不同,靠近活塞的2號(hào)鋼片溫度更低;在徑向上,溫度存在明顯的不同,內(nèi)徑溫度比外徑溫度低得明顯。
根據(jù)圖6,假定仿真的有效滑摩總時(shí)長(zhǎng)為20 s,截取3.0~23.0 s的試驗(yàn)區(qū)間,對(duì)試驗(yàn)壓力曲線做最小二乘擬合,得到接合壓力的幅值歷程量綱為一函數(shù):
式中:p為壓力;t為滑摩時(shí)間。
導(dǎo)入Abaqus徑向節(jié)點(diǎn)壓力數(shù)據(jù),擬合3次樣條函數(shù) p(r),最終獲得各接觸面壓力函數(shù) p(r,t),作為輸入賦予溫度模型。短時(shí)間接合內(nèi),忽略卡簧和花鍵齒有限的導(dǎo)熱和散熱作用,得到溫度結(jié)果:
1)由圖8可知,A4點(diǎn)仿真溫度明顯高于A2號(hào)鋼片,分別為98、97.1℃;B4點(diǎn)溫度略低于B2點(diǎn),分別為77.2、77.8℃;C4點(diǎn)溫度略低于C2點(diǎn),分別為51.8、53.8℃。B點(diǎn)和C點(diǎn)處的溫度與試驗(yàn)結(jié)果十分接近,A點(diǎn)處溫度與試驗(yàn)值誤差大約為3.3%。仿真結(jié)果和試驗(yàn)結(jié)果基本一致,曲線走向相互吻合。因此,仿真具有一定的正確性。
圖8 徑向測(cè)點(diǎn)溫度仿真曲線
2)由圖9、10可知,滑摩面的溫度要高于盤(pán)片內(nèi)部溫度和非產(chǎn)熱面。高溫區(qū)出現(xiàn)在S6面外徑上,因此靠近集中載荷一側(cè)的摩擦元件更易產(chǎn)生熱失效??拷钊粋?cè)的摩擦表面內(nèi)徑溫度比卡簧側(cè)的溫度更高,外徑溫度則更低。
圖9 滑摩面軸向溫度場(chǎng)
圖10 滑摩面徑向溫度場(chǎng)
載荷約束會(huì)使溫度場(chǎng)分布不均,為了避免元件翹曲變形,應(yīng)調(diào)整約束形式使鋼片表面徑向溫度場(chǎng)趨于均勻。通?;钊P(pán)的壓力施加方式有幾種,除了上述的施加方案,還有壓指式的活塞盤(pán),常見(jiàn)的壓指式活塞盤(pán)如圖11所示。
圖11 4種活塞施壓方案示意圖
為保證摩擦副傳扭能力不變,接合壓力按式(15)計(jì)算,仿真結(jié)果如圖12~15所示。
圖12 壓指中置滑摩溫度場(chǎng)
圖14 壓指內(nèi)置滑摩溫度場(chǎng)
圖15 內(nèi)外雙壓指滑摩溫度場(chǎng)
由圖12~15可知,當(dāng)單個(gè)壓指推動(dòng)活塞時(shí),壓指處的溫度更高,壓指內(nèi)置的軸向溫度差異明顯。當(dāng)存在2個(gè)壓指的集中載荷時(shí),在2個(gè)壓指處形成區(qū)域高溫,且靠近外徑的壓指處溫度最高。
表4是5種方案的溫度結(jié)果,其中(e)為活塞壓力均布的方案,仿真結(jié)果數(shù)據(jù)如圖9所示。由表4數(shù)據(jù)和圖12可知,當(dāng)卡簧布置于離合器外側(cè),壓指中置的(a)方案溫度場(chǎng)分布均勻,最高溫度最小,能有效避免摩擦表面燒蝕,而且該方案的徑向溫差最小,可以避免翹曲變形,所以相對(duì)來(lái)說(shuō),壓指中置的活塞盤(pán)結(jié)構(gòu)最優(yōu),能夠有效地延長(zhǎng)離合器壽命。
其次是雙壓指(d)和均布?jí)毫Γ╡)的方案,這2種方案效果相當(dāng),徑向溫差和最高溫度僅次于壓指中置(a)方案,相比(a)方案,這2種方案軸向溫差更小,軸向溫差小的優(yōu)點(diǎn)具體表現(xiàn)在當(dāng)摩擦副增多時(shí),對(duì)摩擦副的破壞更小。
壓指外置(c)和壓指內(nèi)置(b)的方案,最高溫度的值最大,徑向溫差也遠(yuǎn)大于其他3種方案。
表4 5種約束方案溫度差異 ℃
1)在活塞和卡簧的約束作用下,不同摩擦表面的壓力并不相同,而對(duì)于摩擦表面,壓力分布也是不均勻的,因此,要將離合器所有的摩擦副納入仿真,有必要考慮活塞和卡簧對(duì)溫度場(chǎng)的影響。
2)受約束影響,多片副高溫高壓區(qū)將靠近集中載荷,使靠近集中載荷元件產(chǎn)生明顯的徑向溫差。
3)在傳扭能力不變的前提下,分析了5種活塞約束下的多片副溫度場(chǎng),發(fā)現(xiàn)活塞盤(pán)壓力集中在中徑處導(dǎo)致的溫升和徑向溫差最小。
重慶理工大學(xué)學(xué)報(bào)(自然科學(xué))2021年4期