研究馬守東,郭立新,孫德華
(東北大學(xué) 機(jī)械工程與自動化學(xué)院,遼寧 沈陽 110819)
?
考慮時變速度的盤式制動系統(tǒng)的顫振特性
研究馬守東,郭立新,孫德華
(東北大學(xué) 機(jī)械工程與自動化學(xué)院,遼寧 沈陽 110819)
制動過程中的黏滑與顫振是一個與摩擦密切相關(guān)的復(fù)雜動力學(xué)問題.為了研究制動盤與摩擦片的黏滑與顫振現(xiàn)象,建立了制動系統(tǒng)考慮制動盤旋轉(zhuǎn)自由度的五自由度非線性動力學(xué)模型.利用數(shù)值仿真方法,通過改變制動過程中的制動壓力和初速度,研究速度時變的動力學(xué)特性,從而探索其黏滑與顫振現(xiàn)象產(chǎn)生的相關(guān)原因.研究結(jié)果表明:制動盤和摩擦片的切向振幅隨著制動初速度的增加,呈先增大后減小的趨勢;低速時制動盤和摩擦片切向?qū)a(chǎn)生明顯的黏滑顫振,隨著制動初速度的增加,黏滑顫振逐漸減小;隨著制動壓力的增加,切向和法向的振動有增強(qiáng)趨勢,系統(tǒng)由純滑動逐漸進(jìn)入黏滑運(yùn)動并產(chǎn)生顫振.
盤式制動系統(tǒng); 黏滑; 制動顫振; 非線性動力學(xué); 數(shù)值仿真
制動系統(tǒng)是車輛安全性的重要組成部分.如果車輛制動系統(tǒng)設(shè)計不合理,會使車輛在低速制動工況下出現(xiàn)黏滑和顫振現(xiàn)象[1-2],并伴有大量噪聲,嚴(yán)重影響舒適度.
國內(nèi)外學(xué)者對制動系統(tǒng)的振動摩擦問題進(jìn)行了一系列研究[3-5],楊鳳紅等[6]對汽車盤式制動系統(tǒng)建立了二自由度的非光滑動力學(xué)模型,并通過數(shù)值模擬方法研究了制動盤的轉(zhuǎn)速對系統(tǒng)動力學(xué)行為的影響,結(jié)果表明:盤式制動系統(tǒng)中存在著混沌運(yùn)動和周期加倍等復(fù)雜非線性動力學(xué)現(xiàn)象.SHIN等[7]建立了二自由度的動力學(xué)模型,通過改變參數(shù)研究系統(tǒng)的穩(wěn)定性,研究發(fā)現(xiàn)單一增加制動盤或摩擦片的阻尼會使系統(tǒng)變得不穩(wěn)定,伴隨大量噪聲.KINKAID等[8]闡述了盤式制動器噪聲的一種新理論,即多自由度摩擦振動的瞬態(tài)動力學(xué)問題.他們通過研究一組制動塊和轉(zhuǎn)子組成的四自由度模型的動力學(xué)特性,發(fā)現(xiàn)當(dāng)滑動速度在制動方向上很低時,系統(tǒng)滑動速度方向沿徑向有一個突變,這種突變將引起復(fù)雜的振動.但他們大多以二自由度模型的制動盤和摩擦片為研究對象,只研究了摩擦片和制動盤水平方向的振動,并且研究的是制動速度不變、穩(wěn)定狀態(tài)下的振動特性,由于制動時法向方向的振動會引起摩擦力的變化,從而導(dǎo)致一些非線性現(xiàn)象,所以制動系統(tǒng)切向和法向的振動都值得研究.
因此,本文在考慮制動盤和摩擦片切向和法向剛度和阻尼基礎(chǔ)上,又考慮了制動盤的旋轉(zhuǎn)因素,建立了五自由度汽車制動系統(tǒng)的非線性動力學(xué)模型,主要研究系統(tǒng)隨時間變化、制動盤與摩擦片在其相對速度變化過程中的黏滑和顫振現(xiàn)象.
將盤式制動系統(tǒng)簡化成如圖1所示的五自由度動力學(xué)模型.
圖1 制動系統(tǒng)動力學(xué)模型Fig.1 Dynamic model of braking systems
上述模型中考慮制動盤和摩擦片切向、法向的振動和制動盤的旋轉(zhuǎn),由牛頓第二定律可得該制動系統(tǒng)的動力學(xué)方程:
(1)
式中:Fμ為動摩擦力;h為摩擦片與制動盤的中心距;J為制動盤的轉(zhuǎn)動慣量.
系統(tǒng)的運(yùn)動形式由黏滑運(yùn)動時的靜摩擦力和純滑動時的動摩擦力決定.對于黏滑運(yùn)動,靜摩擦力小于等于最大靜摩擦力,Fs≤μsN,并且靜摩擦力等于制動盤和摩擦片之間切向的作用力與反作用力之和[9].
因此,摩擦力Fμ的表達(dá)式為
(2)
式中:Fs為兩接觸表面黏著時的靜摩擦力;μs為靜摩擦因數(shù);μ(vr)為動摩擦因數(shù)[10];vr和FN分別為摩擦片與制動盤間的相對滑動速度和接觸壓力.
(3)
(4)
(5)
(6)
(7)
限于篇幅原因,僅研究制動初速度和制動壓力的變化對盤式制動系統(tǒng)特性振動的影響,阻尼、剛度等因素對制動系統(tǒng)振動特性的影響將另行討論.
2.1 制動初速度對制動顫振的影響
為了與現(xiàn)有文獻(xiàn)[11-12]結(jié)果具有一定的對比性,取γ=0.03,=0.000 013 463 7,δ=0.05,μs=0.6,m1=9 kg,m2=9 kg,c1=c2=c3=c4=1 N·s·m-1,k1=k2=k3=k4=63 N·m-1,N=40 N,分別取初始角速度ω0=5,10,100 rad·s-1進(jìn)行仿真,觀察其對應(yīng)的位移、相圖和相對速度變化特征,研究其對制動顫振的影響.
圖2,3,4分別為制動初速度ω0=5,10,100 rad·s-1的時間歷程及其相軌跡.從摩擦片和制動盤的水平位移曲線可以看出,制動盤和摩擦片切線方向的振幅是隨著角速度的增加呈現(xiàn)出先增大后減小的趨勢.從制動盤和摩擦片法向相圖可以看出,隨著角速度的增加,法向方向振幅不變并一直處于無黏滯周期運(yùn)動狀態(tài),最終達(dá)到穩(wěn)定,可見制動初速度對制動盤和摩擦片豎直方向的振動沒有影響.從制動盤角速度曲線可以看出,隨著制動初速度的增加達(dá)到穩(wěn)定的時間增加,并且在低速時制動角速度會有較大波動,這也會引起摩擦片和制動盤切線方向的較大振動.從制動盤和摩擦片的相圖與相對速度曲線可以看出,低速時產(chǎn)生黏滑現(xiàn)象,在位移圖像中反映了存在顫振現(xiàn)象,隨著制動速度的增加,黏滑現(xiàn)象逐漸減小直至消失.
圖2 ω0=5 rad·s-1的時間歷程及其相軌跡Fig.2 Time history and phase trajectory when ω0=5 rad·s-1
圖3 ω0=10 rad·s-1的時間歷程及其相軌跡Fig.3 Time history and phase trajectory when ω0=10 rad·s-1
2.2 制動壓力對制動顫振的影響
研究在其他參數(shù)保持不變的情況下,通過改變制動壓力來研究其對制動顫振的影響(m1=9 kg,m2=9 kg,c1=c2=c3=c4=1 N·s·m-1,k1=k2=k3=k4=63 N·m-1,ω0=10 rad·s-1,N取10,100 N進(jìn)行仿真分析).
圖5,6分別為制動壓力為N=10,100 N時的時間歷程及其相軌跡.從制動盤和摩擦片的位移曲線可以看出,隨著制動壓力的增加,制動盤和摩擦片切向和法向的振幅都增大.從角速度曲線可以看出,隨著制動壓力的增加角速度趨近零時產(chǎn)生的波動也增加,從而導(dǎo)致制動盤和摩擦片切向的顫振.從相圖和相對速度圖可以看出,隨著制動壓力的增加,法向振動一直做無黏滯周期運(yùn)動并逐漸趨于穩(wěn)定,制動盤和摩擦片切線方向由純滑動逐漸進(jìn)入黏滑運(yùn)動狀態(tài),并且隨著制動壓力的增大,黏滑階段持續(xù)的時間也逐漸增加.
圖4 ω0=100 rad·s-1的時間歷程及其相軌跡Fig.4 Time history and phase trajectory when ω0=100 rad·s-1
圖5 N=10 N的時間歷程及其相軌跡Fig.5 Time history and phase trajectory when N=10 N
(1) 用數(shù)值仿真方法研究了制動初速度和制動壓力對制動盤和摩擦片制動穩(wěn)定性的影響,仿真結(jié)果表明:在低速時,系統(tǒng)存在明顯的黏滯現(xiàn)象,并且角速度有很大波動,從而導(dǎo)致切線方向產(chǎn)生較大的顫振;隨著制動初速度的增加,切向顫振減弱,黏滑現(xiàn)象也逐漸減弱并進(jìn)入純滑動.
圖6 N=100 N的時間歷程及其相軌跡Fig.6 Time history and phase trajectory when N=100 N
(2) 隨著制動壓力的增加,制動時間減小,角速度趨近零時產(chǎn)生的波動增大,制動盤和摩擦片切向的黏滑顫振越明顯,法向振幅也增加.這為合理地操縱車輛、減少制動過程中產(chǎn)生的黏滑顫振提供了依據(jù).
[1] KANG J Y,KROUSGRILL C M.Comprehensive stability analysis of disc brake vibrations including gyroscopic,negative friction slope and mode-coupling mechanisms[J].Journal of Sound and Vibration,2009,324(1):387-407.
[2] IBRAHIMR A,MADHAVAN S,QIAO S L,et al. Experimental investigation of friction-induced noise in disc brake systems[J].International Journal of Vehicle Design,2000,23(3):218-240.
[3] NAKAE T,RYU T,SUEOKA A,et al.Squeal and chatter phenomena generated in a mountain bike disc brake[J].Journal of Sound and Vibration,2011,330(10):2138-2149.
[4] OBERST S,LAI J C S.Nonlinear transient and chaotic interactions in disc brake squeal[J].Journal of Sound and Vibration,2015,342(4):272-289.
[5] 劉獻(xiàn)棟,任增杰,王海霞,等.盤式制動器摩擦特性及制動尖叫測試與分析[J].振動、測試與診斷,2013,33(5):746-750.
LIU X D,REN Z G,WANG H X,et al.Disc brake friction characteristics,braking scream test and analysis[J].Journal of Vibration,Measurement & Diagnosis,2013,33(5):746-750.
[6] 楊鳳紅,張偉.盤式制動系統(tǒng)的非線性動力學(xué)分析[J].振動與沖擊,2009,28(6):93-99.
YANG F H,ZHANG W.Nonlinear dynamical analysis for a disc braking system[J].Journal of Vibration and Shock,2009,28(6):93-99.
[7] SHIN K,BRENNAN M J,HARRIS C J,et al.Analysis of disc brake noise using a two-degree-of-freedom model[J].Journal of Sound and Vibration,2002,254(5):837-848.
[8] KINKAID N M,OREILLY O M,PAPADOPOULOS P.On the transient dynamics of multi-degree-of-freedom friction oscillator:a new mechanism for disc brake noise[J].Journal of Sound and Vibration,2005,287(4):901-917.
[9] SHIN K,BRENNAN M J.Analysis of disc brake noise using a two-degree-of-freedom model[J].Journal of Sound and Vibration,2002,254(5):837-848.
[10] POPP K,STELTER P.Stick-slip vibrations and chaos[J].Philosophical Transactions of the Royal Society of London,1990,332(1624):89-105.
[11] 宋春元.車輛制動系統(tǒng)顫振研究[D].成都:西南交通大學(xué),2006.
SONG C Y.Investigation into chatter of vehicle brake system[D].Chengdu:Southwest Jiaotong University,2006.
[12] 李博.剎車系統(tǒng)干摩擦非光滑動力學(xué)分析[D].天津:天津大學(xué),2007.
LI B.Analysis of non-smooth dynamics of brake system with dry-friction[D].Tianjin:Tianjin University,2007.
Flutter property study on disc brake system with time-varied velocities
MA Shoudong,GUO Lixin,SUN Dehua
(School of Mechanical Engineering & Automation, Northeastern University, Shenyang 110819,Liaoning,China)
In general, the stick-slip and flutter during braking process are associated with a complex dynamic problem, which is closely related to friction. In order to study the stick-slip and flutter between the brake and friction disks, the nonlinear dynamic model of five degree freedoms is first established including the rotational freedom of brake disc. Then, the dynamic properties of braking system are investigated through changing brake pressure and initial velocity by using numerical simulation software. As such, the reasons for stick-slip and flutter can be speculated. Finally, it is shown from results that the tangent vibration amplitude increases initially, and thereafter decreases with increase of brake disc initial velocity. When the initial velocity of brake disc is low, the serious stick-slip and flutter occur. With the increase of initial velocity, it gradually weakens. When the tangent and normal vibration amplitude increase with increase of braking pressure, the stick-slip and flutter of braking system occur.
braking system; stick-slip; braking flutter; nonlinear dynamics; numerical simulation
國家自然科學(xué)基金資助項目(51275082)
馬守東(1990-),男,碩士.E-mail:msd12142838@163.com
TH 113.1
A
1672-5581(2017)01-0047-05