電流變液靜壓導軌系統(tǒng)靜動態(tài)特性研究

胡均平，劉成沛，郭 勇，梁向京

(1. 中南大學 機電工程學院,長沙 410083; 2. 湖南科技大學 機電工程學院,湖南 湘潭 411201)

?

電流變液靜壓導軌系統(tǒng)靜動態(tài)特性研究

胡均平1，劉成沛1，郭 勇2，梁向京1

(1. 中南大學 機電工程學院,長沙 410083; 2. 湖南科技大學 機電工程學院,湖南 湘潭 411201)

研究了一種以電流變液為潤滑介質(zhì)的靜壓導軌系統(tǒng)靜動態(tài)特性。利用液阻網(wǎng)絡(luò)理論，在分析電場強度與電流變液黏性關(guān)系，綜合考慮關(guān)聯(lián)參數(shù)的基礎(chǔ)上，建立了電流變液靜壓導軌系統(tǒng)的CFD計算數(shù)值模型。根據(jù)攝動理論，利用Fluent動網(wǎng)格技術(shù)，研究了在不同工作變量下電流變液靜壓導軌的靜動態(tài)特性。研究結(jié)果表明：電場強度、負載因素與初始壓力比等工作變量均會對電流變液靜壓導軌系統(tǒng)性能帶來影響；電場強度增強會使系統(tǒng)流量減小，阻尼系數(shù)增大，但不會對系統(tǒng)靜剛度帶來影響；改變切削力載荷頻率，導軌系統(tǒng)會出現(xiàn)共振現(xiàn)象，增強電場強度可顯著提高系統(tǒng)動剛度，抑制共振時的最大振幅。實驗結(jié)果與仿真數(shù)值一致，為電流變液靜壓導軌的設(shè)計應(yīng)用提供了理論參考依據(jù)。

靜壓導軌；電流變液；電場強度；靜動態(tài)特性；動網(wǎng)格；關(guān)聯(lián)參數(shù)

液體靜壓導軌以其高精度、高剛度、高阻尼吸振性、低摩擦和長壽命等技術(shù)優(yōu)勢，在超精密加工領(lǐng)域獲得了越來越廣泛的應(yīng)用。使用電流變液作為潤滑介質(zhì)取代液壓油，是智能材料在潤滑領(lǐng)域的重要應(yīng)用，可以顯著提高液體靜壓支承系統(tǒng)的性能。目前針對電流變液特性及其在潤滑領(lǐng)域的應(yīng)用方面：KO等[1]介紹了一種負電流變液的材料機理、性能及應(yīng)用范圍；BOUZIDANCE等[2]研究了基于電流變液的靜壓徑向軸承因轉(zhuǎn)子轉(zhuǎn)動不平衡引起的慣性振動的抑制方法；AOYAMA等[3]分析了電流變液在機床工作臺與軸承等關(guān)鍵區(qū)域的應(yīng)用方案。靜壓支承系統(tǒng)主要性能參數(shù)的求解方法主要有：從雷諾方程出發(fā)，采用解析方法進行推導[4-5]；使用有限差分法進行數(shù)值求解；根據(jù)攝動理論，利用CFD軟件Fluent動網(wǎng)格技術(shù)求解剛度與阻尼系數(shù)的方法[6]。在靜壓支承系統(tǒng)靜動態(tài)特性研究方面：CHEN等[7]研究了靜壓主軸的靜動態(tài)特性；THOMAS等[8]研究了不同油腔結(jié)構(gòu)參數(shù)對靜壓徑向軸承的等效剛度與阻尼的影響情況；DU等[9]在考慮油液可壓縮性與慣性效應(yīng)的基礎(chǔ)上建立了靜壓導軌系統(tǒng)的質(zhì)量、剛度、阻尼等效模型；楊小高等[10]對基于固定節(jié)流的靜壓徑向軸承進行了動態(tài)特性分析；于賀春等[11]靜壓氣體軸承轉(zhuǎn)子系統(tǒng)提出了流場與轉(zhuǎn)子的準動態(tài)動力學特性耦合分析方法。另外導軌系統(tǒng)結(jié)構(gòu)參數(shù)對油液溫度場的影響[12]以及油液溫度因素對靜壓支承系統(tǒng)性能的影響[13]也經(jīng)常出現(xiàn)在文獻之中。但是基于電流變液為潤滑介質(zhì)的靜壓導軌系統(tǒng)的靜動態(tài)特性研究領(lǐng)域卻鮮有報道。本文作者在分析電流變液特性與電場強度關(guān)系，并綜合考慮關(guān)聯(lián)參數(shù)的基礎(chǔ)上，建立了基于關(guān)聯(lián)參數(shù)的CFD計算數(shù)值模型；根據(jù)攝動理論，利用Fluent動網(wǎng)格技術(shù)研究了不同工作變量對靜壓導軌系統(tǒng)性能參數(shù)的影響規(guī)律；以及為在切削力載荷下提高系統(tǒng)動剛度并減小導軌共振振幅而選取合適的電場強度提供理論依據(jù)。

1 電流變液靜壓導軌系統(tǒng)分析建模

1.1 電流變液靜壓導軌系統(tǒng)原理

圖1 電流變液靜壓導軌系統(tǒng)原理圖Fig.1 Scheme of hydrostatic guideway system based on electrorheological fluid

圖1是電流變液靜壓導軌系統(tǒng)原理圖。電流變液由液壓泵泵出，分別經(jīng)毛細管節(jié)流器進入上下油腔，在導軌溜板與固定座之間形成油膜，以抵抗溜板與工件重量和外載荷的作用。它屬于毛細管節(jié)流器補償?shù)暮銐汗┯烷]式電流變液靜壓導軌。

該電流變液靜壓導軌上下油腔之間尺寸參數(shù)不一致，但其結(jié)構(gòu)型式均為為矩形油腔，如圖2所示。

圖2 矩形油腔結(jié)構(gòu)圖Fig.2 Structure of rectangle oil pocket

靜壓導軌在工作過程中受到的外載荷分為兩大類：靜態(tài)載荷與動態(tài)載荷。靜態(tài)載荷主要分為固定重量(導軌溜板)與可變重量(工件與夾具)，亦可統(tǒng)稱為負載因素；動態(tài)載荷主要指加工過程中的切削力載荷。

F=Mg+mg+f(t)

(1)

式中：F為靜壓導軌總負載，N；M為導軌溜板質(zhì)量，kg；m為工件與夾具質(zhì)量，kg；g為重力加速度，m/s2；f(t)為切削力載荷，N。

定義負載比α為靜態(tài)載荷與固定質(zhì)量(導軌溜板)之比，用來定量地表征負載因素。

(2)

在初始設(shè)計狀態(tài)，導軌在固定質(zhì)量以及上下油腔作用力下保持靜態(tài)平衡，由力平衡條件可知：

n(Pr0Ae2-Pr0Ae1)=Mg

(3)

Aei=(Li+li)(Bi+bi)/4

(4)

式中：Pr0為初始設(shè)計狀態(tài)各油腔壓力，Pa；Aei為油腔有效承載面積，m2；n為油腔對數(shù)；Li為油墊總長度，m；Bi為油墊總寬度，m；li為油腔長度，m；bi為油腔寬度，m；其中下標i為1時指下油腔及相關(guān)參數(shù)，i為2時指上油腔及相關(guān)參數(shù)，下同。

初始壓力比β0是靜壓導軌系統(tǒng)重要的初始設(shè)計參數(shù)，其定義為為初始設(shè)計狀態(tài)下油腔壓力Pr0與供油壓力Ps之比。調(diào)節(jié)Ps大小可改變β0。

(5)

在負載因素影響的任意狀態(tài)，根據(jù)力平衡條件：

n(Pr2Ae2-Pr1Ae1)=αMg

(6)

值得注意的是，當α=1時，Pr1=Pr2=Pr0。此時初始設(shè)計狀態(tài)下的式(3)是式(6)的特殊情況。

根據(jù)油路液阻網(wǎng)絡(luò)理論[14]，由連續(xù)性條件，液體靜壓導軌系統(tǒng)有如下關(guān)系：

(7)

對于毛細管節(jié)流器，進油液阻Rgi表達式為:

(8)

對于矩形油腔，出油液阻Rhi表達式為:

(9)

(10)

同樣在初始設(shè)計狀態(tài)下Δh=0，則h1=h2=h0。由式(5)及式(7)～式(10)可求得：

(11)

由式(7)可知，如果使初始狀態(tài)下上下油膜厚度相等，須使RgiAbi=C，C為常數(shù)。所以根據(jù)初始壓力比β0與C可確定初始膜厚h0。

由式(6)～式(10)可知，在負載因素下油膜厚度的變化量Δh，又因為Δh值較小，故略去它的高次方項。所以油膜厚度變化量Δh可表達為:

(12)

由式(12)可知，在負載因素作用下，Δh是負載比α的函數(shù)；且當α=1時，Δh=0。同時在可求得任意狀態(tài)下，導軌上下油膜厚度：

hi=h0-(-1)iΔh

(13)

進一步，任意狀態(tài)下導軌上下油腔壓力：

(14)

圖3 電流變液表觀黏度與電場強度關(guān)系Fig.3 Apparent viscosity versus electric field of ER fluid

由式(12)～式(14)及圖3可知，hi、Pri與μ均是E、α及β0的函數(shù)，前三者構(gòu)成了一組關(guān)聯(lián)參數(shù)Θ。將工作變量記作Γ={E,α,β0}，則關(guān)聯(lián)參數(shù)Θ為：

Θ={μ,hi,Pri}=f(Γ)

(15)

關(guān)聯(lián)參數(shù)Θ中的各變量不可隨意組合，它們是聯(lián)動變化的?？赏ㄟ^對Γ賦值以求得Θ，為電流變液靜壓導軌系統(tǒng)建模提供參數(shù)依據(jù)。

1.2 電流變液靜壓導軌系統(tǒng)模型

基于動網(wǎng)格技術(shù)的導軌性能參數(shù)計算方法是在計算流體動力學(Computational Fluid Dynamics, CFD)軟件Fluent的基礎(chǔ)上，采用自定義程序?qū)崿F(xiàn)導軌間隙擾動功能，結(jié)合差分計算模型，計算導軌油膜性能參數(shù)的一種新方法。但該文獻既未從整個系統(tǒng)的角度來考慮因工作變量不同所導致關(guān)聯(lián)參數(shù)變化的聯(lián)動效應(yīng)，也未從電流變液作為特殊潤滑介質(zhì)的角度來研究靜壓支承系統(tǒng)的模型。

靜壓導軌在某一關(guān)聯(lián)參數(shù)下的平衡位置受到小擾動作用(位移擾動與速度擾動)時，油膜力會發(fā)生變化。根據(jù)攝動理論，對油膜力進行一階泰勒級數(shù)線性化處理，油膜力變化與位移擾動和速度擾動之間的關(guān)系[17]可以表達為:

(16)

這樣在任一關(guān)聯(lián)參數(shù)Θ下導軌系統(tǒng)靜剛度與阻尼的計算公式為:

(17)

所以求解式(17)關(guān)鍵是要解決四個技術(shù)問題：其一是為電流變液選取適當?shù)挠嬎隳Ｐ停黄涠墙⒒陉P(guān)聯(lián)參數(shù)Θ的數(shù)學模型；其三是采用動網(wǎng)格技術(shù)施加小擾動并界定擾動范圍；其四是求解計算模型得到油膜壓力分布與油膜力。

1)選取電流變液計算模型。當未加電場時，電流變液為牛頓流體，可參照一般油液計算模型；施加電場后，電流變液呈現(xiàn)為賓漢體，其模型為:

(18)

FLUENT提供了一種擴展的Herschel-Bulkley模型，它綜合反映了流體中的賓漢及指數(shù)行為效應(yīng)[18]：

(19)

當k=μ、λ=1且(τ0/μ0)→0時，F(xiàn)LUENT提供的Herschel-Bulkley模型便與文中需要的賓漢模型一致，可以用作電流變液的計算模型。該模型可在軟件中材料特性的黏度欄選取。

2)關(guān)聯(lián)參數(shù)Θ須運用到計算模型中去。其中油膜厚度hi決定了上下油膜流場三維模型的結(jié)構(gòu)；油液黏度μ決定了電流變液的參數(shù)；油腔壓力Pri決定了模型的壓力邊界條件。由于關(guān)聯(lián)參數(shù)Θ是工作變量集合Γ={E,α,β0}的函數(shù)，因此不同的工作變量Γ就會導致完全不同的計算模型。關(guān)聯(lián)參數(shù)Θ與計算模型的關(guān)系如圖4所示。

圖4 關(guān)聯(lián)參數(shù)與計算模型關(guān)系圖Fig.4 Relationship between correlation parameters and computational model

3)在位移擾動與速度擾動情況下，油膜幾何形狀隨油膜發(fā)生變化，必須引入動網(wǎng)格技術(shù)來進行處理。此時油膜與溜板之間的邊界層網(wǎng)格與油膜內(nèi)部網(wǎng)格均產(chǎn)生變形。油膜網(wǎng)格運動的計算模型采用的是基于線性假設(shè)的彈性光順Smoothing模型，其核心思想是基于胡克定律，網(wǎng)格線類似于彈簧，兩端節(jié)點作彈性移動的網(wǎng)格運動方式。油膜厚度擾動程序建立在Fluent軟件提供的用戶自定義函數(shù)(User Define Function, UDF)的基礎(chǔ)上，其功能通過宏DEFINE_PROFILE實現(xiàn)[19]：

((BOTTOMWALL_profile 2 point)

(time 0t)

(v_zVhVh))

在該UDF程序中，t指的是網(wǎng)格運動終結(jié)時間，Vh是指油膜厚度網(wǎng)格運動速度。參照熊萬里等提出的原則，即通過對比研究動網(wǎng)格瞬態(tài)計算和穩(wěn)態(tài)計算在平衡點的變化，分析位移擾動與速度擾度對油膜力的影響，進而界定其取值范圍：

位移擾動(求靜剛度)：Vh≤1 μm/s；Vht≤0.2 μm

速度擾動(求阻尼)：Vh≥0.1 mm/s；Vht≤0.2 μm

4)在用SolidWorks建立的三維油膜模型文件x_t的基礎(chǔ)上，采用前處理軟件ICEM CFD劃分網(wǎng)格，油膜厚度方向網(wǎng)格為10層。采用Fluent軟件求解模型，不同電場強度及負載因素下的參數(shù)取值遵從關(guān)聯(lián)參數(shù)Θ。經(jīng)過迭代數(shù)值計算，得到油膜流場分布。導軌系統(tǒng)油膜力W可進一步求得：

(20)

由上述的動網(wǎng)格瞬態(tài)計算最終可求得關(guān)聯(lián)參數(shù)下的導軌系統(tǒng)靜剛度KΘ與阻尼CΘ系數(shù)。另外通過穩(wěn)態(tài)計算還可求得電流變液靜壓導軌系統(tǒng)總流量Q：

(21)

1.3 導軌系統(tǒng)動態(tài)特性分析

電流變液靜壓導軌的油膜相當于帶有阻尼的液體彈簧，故文中的導軌系統(tǒng)也可以簡化為一個單自由度的質(zhì)量—剛度—阻尼系統(tǒng)。假設(shè)作用在電流變液靜壓導軌上的載荷為波動載荷，此時考慮不同關(guān)聯(lián)參數(shù)Θ下的系統(tǒng)振動模型為：

(22)

式中:fm為切削力幅值，N；ω為切削力頻率，Hz。

式(22)是一個典型的具有黏性阻尼系統(tǒng)的單自由度受迫振動。如果忽略瞬態(tài)振動而僅僅研究其穩(wěn)態(tài)行為，則導軌系統(tǒng)的動態(tài)響應(yīng)方程為:

Δh(t)=Bsin(ωt-φ)

(23)

式中:B為導軌系統(tǒng)振動振幅，m；φ為系統(tǒng)位移落后于切削力的相位角，rad。

由式(23)可求得靜壓導軌系統(tǒng)的振幅放大因子：

(24)

導軌系統(tǒng)動剛度KdΘ可由其靜剛度求得：

(25)

2 系統(tǒng)靜動態(tài)特性仿真結(jié)果分析

為了定量地分析工作變量Γ (含電場強度、初始壓力比及負載比)對系統(tǒng)性能參數(shù)的影響規(guī)律，以某精密機械研究所超精密加工系統(tǒng)靜壓導軌油墊為例進行分析。其基本參數(shù)見表1。

除了表1所列舉的基本參數(shù)外，導軌系統(tǒng)還有一些基礎(chǔ)參數(shù)：如供油壓力Ps=2.5 MPa、導軌溜板質(zhì)量M=500 kg等。另外還有前文所述的關(guān)于電流變液的特性參數(shù)等。

表1 靜壓導軌系統(tǒng)基本參數(shù)Tab.1 Parameters of hydrostatic guideway system

2.1 電場強度與初始壓力比的影響

流量Q、油膜靜剛度KΘ與阻尼CΘ是電流變液靜壓導軌系統(tǒng)的重要性能參數(shù)。本小節(jié)均是在初始設(shè)計狀態(tài)(負載比α=1)下，研究不同電場強度與初始壓力比對性能參數(shù)影響的。

圖5為系統(tǒng)流量Q隨電場強度E和初始壓力比β0變化曲線。圖中顯示系統(tǒng)流量Q隨電場強度E和初始壓力比β0的增大而減小。這是由于當電場強度增大時，電流變液黏度會增加，而系統(tǒng)流量與黏度成反比例，變會減小。當β0=1時，系統(tǒng)流量Q減小為0，這是由于導軌達到了供油壓力的承載極限。

圖5 系統(tǒng)流量隨電場強度和初始壓力比變化曲線Fig.5 Flow rate versus pressure ratio for different electric fields

圖6為系統(tǒng)靜剛度K隨電場強度E和初始壓力比β0變化曲線。圖中顯示電場強度E對電流變液靜壓導軌系統(tǒng)靜剛度K沒有影響，因為使用毛細管補償?shù)膶к夓o剛度均與油液黏度無關(guān)，自然也與電場強度無關(guān)；靜剛度K隨初始壓力比β0增大先變大，后減小，并且在β0=0.67處達到峰值，此時K=2 989 N/μm，這與BOUZIDANCE等分析靜壓軸承時得到的結(jié)論一致。

圖7為系統(tǒng)阻尼C隨電場強度E和初始壓力比β0變化曲線。圖中顯示系統(tǒng)阻尼C隨電場強度E和初始壓力比β0的增大而增大，由于阻尼C的變化幅度較大，縱坐標采用對數(shù)坐標。初始壓力比與電場強度分別通過對初始膜厚h0及電流變液黏度μ的影響來實現(xiàn)對系統(tǒng)阻尼影響的。

圖6 系統(tǒng)靜剛度隨電場強度和初始壓力比變化曲線Fig.6 Stiffness versus pressure ratio for different electric fields

圖7 系統(tǒng)阻尼隨電場強度和初始壓力比變化曲線Fig.7 Damping versus pressure ratio for different electric fields

2.2 電場強度與負載比的影響

負載比α也會對電流變液靜壓導軌系統(tǒng)性能參數(shù)造成影響。本小節(jié)均是在初始設(shè)計狀態(tài)(初始壓力比β0=0.6)下，研究不同電場強度與負載比對性能參數(shù)影響的。

圖8為系統(tǒng)靜剛度K隨電場強度E和負載比α變化曲線。圖中顯示系統(tǒng)靜剛度K隨負載比α的增大而減小，這是由于膜厚變化量Δh是負載比α的函數(shù)，而系統(tǒng)靜剛度會隨導軌偏心程度的增大而減?。凰栽趯к壍脑O(shè)計和使用階段，應(yīng)盡量降低負載比的大小，以使導軌工作在較大的系統(tǒng)靜剛度狀態(tài)之中。電場強度E依舊對系統(tǒng)靜剛度K沒有影響。

圖8 系統(tǒng)靜剛度隨電場強度和負載比變化曲線Fig.8 Stiffness versus load factor for different electric fields

圖9為系統(tǒng)阻尼C隨電場強度E和負載比α變化曲線。圖中顯示系統(tǒng)阻尼C隨負載比α的增大而變大；同樣由于電場強度對電流變液黏度的影響，使得系統(tǒng)阻尼C隨電場強度E的增大而變大。

圖9 系統(tǒng)阻尼隨電場強度和負載比變化曲線Fig.9 Damping versus load factor for different electric fields

另外由圖6～圖9可知，對于毛細管節(jié)流器補償?shù)碾娏髯円红o壓導軌，電場強度E對系統(tǒng)靜剛度K變化沒有作用，而對系統(tǒng)阻尼C有著顯著地影響。因此調(diào)節(jié)電場強度可以在不改變系統(tǒng)靜剛度的條件下調(diào)整系統(tǒng)的動態(tài)特性，具有明顯的優(yōu)勢。

2.3 電場強度對動態(tài)特性影響

從前文可知，電場強度E會對電流變液靜壓導軌系統(tǒng)動態(tài)特性帶來影響。本小節(jié)均是在初始設(shè)計狀態(tài)(負載比α=1、初始壓力比β0=0.6)下，研究不同電場強度對系統(tǒng)動態(tài)特性參數(shù)影響的。

圖10為系統(tǒng)振幅放大因子B/Be隨電場強度E和切削力頻率ω變化曲線。圖中顯示在切削力頻率ω增大過程中，振幅放大因子B/Be先變大后減小，存在共振頻率2 830 Hz；增加電場強度可以顯著降低共振過程中的振幅幅值。

圖10 振幅放大因子隨電場強度和切削力頻率變化曲線Fig.10 Vibration amplitude magnification due to machining force frequency for different electric fields

圖11為系統(tǒng)動剛度KdΘ隨電場強度E和切削力頻率ω變化曲線。圖中顯示系統(tǒng)動剛度KdΘ隨切削頻率ω的增大先變大后減小，且在共振頻率附近存在著最小的動剛度數(shù)值。增大電場強度E可以顯著提高共振頻率附近的系統(tǒng)動剛度，如20倍的局部放大圖Ⅰ所示。

圖11 系統(tǒng)動剛度隨電場強度和切削力頻率變化曲線Fig.11 Dynamic stiffness due to machining force frequency for different electric fields

3 實驗研究

實驗平臺是在中國航空工業(yè)集團北京航空精密機械研究所國防科技重點預研項目——“Nanosys-300 非球面超精密加工系統(tǒng)”的基礎(chǔ)上改造而成。

給導軌施加波動載荷，利用德國米依公司的非接觸式位移傳感器capaNCDT6300來測量導軌運動，并通過數(shù)據(jù)采集卡將數(shù)據(jù)傳遞到計算機中。利用HVG5000提供直流電源(電壓0～5 kV/mm，電流0～10 mA)施加在電流變液靜壓導軌上：導軌固定座電極接負極，溜板電極接正極。另外，可以通過改變工件或夾具質(zhì)量改變負載比α，調(diào)整溢流閥壓力調(diào)節(jié)初始壓力比β0的方式來得到導軌不同工作狀態(tài)下的振動響應(yīng)。

表2為電流變液靜壓導軌系統(tǒng)在不同電場強度E下，最大振幅B的實驗值與數(shù)值解比較。實驗中導軌系統(tǒng)處于初始設(shè)計狀態(tài)(負載比α=1、初始壓力比β0=0.6)；切削力f(t)=500sin(2 830 t)，切削力最大振幅為500 N，切削力頻率等于共振頻率大小，即2 830 Hz。從表中可以看出電流變液靜壓導軌系統(tǒng)最大振幅隨電場強度的增大而減小，且數(shù)值解均能較好地與實驗值相符合，差異維持在0.63%～1.42%之間。其差異是由于在導軌油膜流場數(shù)值模擬過程中忽略了油液慣性力所致。實驗結(jié)果證實了文章研究方法及結(jié)果的正確性。

表2 最大振幅B的實驗值與數(shù)值解比較Tab.2 Comparison between test and computational value of vibration amplitude B

4 結(jié) 論

(1) 電場是通過改變電流變液黏度來影響靜壓導軌系統(tǒng)靜動態(tài)特性的。電場強度E越大，系統(tǒng)流量Q就越小，而系統(tǒng)阻尼系數(shù)C則越大；但電場作用對靜壓導軌系統(tǒng)靜剛度K沒有影響。

(2)在電流變液靜壓導軌系統(tǒng)中，負載因素與初始壓力比依然會對系統(tǒng)性能參數(shù)帶來影響。負載比α越大，靜剛度K就會越小，而阻尼C則會變大；初始壓力比β0越大，流量Q就越小，而阻尼C會變大，系統(tǒng)靜剛度K則會隨初始壓力比先大后小，出現(xiàn)峰值為2 989 N/μm。

(3)在切削力載荷頻率ω逐漸增大的過程中，導軌系統(tǒng)動剛度KdΘ先減小后增大，系統(tǒng)振幅放大因子B/Be則先增大后減小，在共振頻率2 830 Hz處存在峰值；增加電場強度E可以有效地增大系統(tǒng)動剛度KdΘ，減小共振振幅。

(4) 在綜合考慮電場強度E對電流變液特性的影響的基礎(chǔ)上，建立了基于關(guān)聯(lián)參數(shù)的導軌系統(tǒng)CFD求解計算模型。研究了在不同工作變量下電流變液靜壓導軌的靜動態(tài)特性。實驗結(jié)果與仿真數(shù)值一致，為電流變液靜壓導軌的設(shè)計應(yīng)用提供了理論參考依據(jù)，具有較好的工程應(yīng)用價值。

[1] KO Y G, CHOI U S. Negative electrorheological fluids[J]. Journal of Rheology, 2013, 57(6): 1655-1667.

[2] BOUZIDANCE A, THOMAS M. An electrorheological hydrostatic journal bearing for controlling rotor vibration[J]. Computers and Structures, 2008, 86(3/4/5): 463-472.

[3] AOYAMA T, INASAKI I. Application of electro-rheological fluid dampers to machine tool elements[J]. Annals of CIRP, 1997, 46(1): 309-312.

[4] 葉紅玲, 鄭小龍, 沈靜嫻, 等. 液體靜壓導軌轉(zhuǎn)臺軸向振動的動力學建模與分析[J]. 工程力學, 2012, 29(3): 218-225. YE Hongling, ZHENG Xiaolong, SHEN Jingxian, et al. Dynamic modeling and analysis of axial vibration of the hydrostatic slide turntable[J]. Engineering Mechanics,2012, 29(3): 218-225.

[5] 趙建華. 龍門車銑中心液體靜壓導軌性能的理論分析與實驗研究[D]. 秦皇島: 燕山大學, 2013: 15-57.

[6] 熊萬里, 侯志泉. 基于動網(wǎng)格模型的液體動靜壓軸承剛度阻尼計算方法[J]. 機械工程學報, 2012, 48(23): 118-126. XIONG Wanli, HOU Zhiquan. Method for calculating stiffness and damping coefficients for hybrid bearings based on dynamic mesh model[J]. Journal of Mechanical Engineering, 2012, 48(23): 118-126.

[7] CHEN Dongju, FAN Jinwei, ZHANG Feihu, et al. Dynamic and static characteristics of a hydrostatic spindle for machine tools[J]. Journal of Manufacturing Systems, 2012, 31(1): 26-33.

[8] BOUZIDANE A, THOMAS M. Equivalent stiffness and damping investigation of a hydrostatic journal bearing[J]. Tribology Transactions, 2007, 50(2): 257-267.

[9] DU Yikang, MAO Kuanmin, ZHU Yaming, et al. Dynamic modeling of hydrostatic guideway considering compressibility and inertial effect[J]. Frontiers of Mechanical Engineering, 2015, 10(1): 78-88.

[10] 楊小高, 王勇勤, 江桂云, 等. 基于固定節(jié)流的徑向滑動軸承動態(tài)特性[J]. 中南大學學報(自然科學版), 2014, 45(9): 2993-3000. YANG Xiaogao, WANG Yongqin, JIANG Guiyun, et al. Dynamic characteristics of hydrostatic fixed throttled journal bearing[J]. Journal of Central South University (Science and Technology), 2014, 45(9): 2993-3000.

[11] 于賀春, 馬文琦, 趙廣, 等. 船舶增壓器靜壓氣體軸承-轉(zhuǎn)子系統(tǒng)動力學特性研究[J]. 振動與沖擊, 2011, 30(12):1-6. YU Hechun, MA Wenqi, ZHAO Guang, et al. Dynamic characteristics of an aerostatic bearing-rotor system with a ship turbocharger[J]. Journal of Vibration and Shock,2011,30(12): 1-6.

[12] SHAO Junpeng. The effect of oil cavity depth on temperature filed in heavy hydrostatic thrust bearing[J]. Journal of Hydrodynamics, 2011, 23(5): 676-680.

[13] MARUIS P. Research concerning the influence of the lubricant temperature on the functioning of hydrostatic guideway systems[J]. Journal of Engineering Studies and Research, 2012, 18(4): 71-77.

[14] 丁振乾. 流體靜壓支承設(shè)計[M]. 上海: 上海科學技術(shù)出版社, 1989: 35-66.

[15] CHO M S, CHOI H J, JHON M S. Shear stress analysis of a semiconducting polymer based electrorheological fluid system[J]. Polymer, 2005,46(25):11484-11488.

[16] CHO M S, CHOI H J, JHON M S. Shear stress analysis of a semiconducting polymer based electrorheological fluid system[J]. Polymer, 2005, 46(25): 11484-11488.

[17] ROWE W B, CHONG F S. Computation of dynamic force coefficients for hybrid (hydrostatic/hydrodynamic) journal bearings by the finite disturbance and perturbation techniques[J]. Tribology International, 1986, 19(5): 260-271.

[18] GERTZOS K P, NIKOLAKOPOULOS P G, PAPADOPOULOS C A. CFD analysis of journal bearing hydrodynamic lubrication by Bingham lubricant[J]. Tribology International, 2008, 41(12): 1190-1204.

[19] 朱紅鈞. ANSYS 14.5 熱流固耦合實戰(zhàn)指南[M]. 北京：人民郵電出版社, 2014: 15-86.

Static and dynamic characteristics of a hydrostatic guideway system based on electrorheological fluid

HU Junping1, LIU Chengpei1, GUO Yong2, LIANG Xiangjing1

(1. College of Mechanical and Electrical Engineering, Central South University, Changsha 410083, China; 2. College of Mechanical and Electrical Engineering, Hunan University of Science and Technology, Xiangtan 411201, China)

The static and dynamic characteristics of a hydrostatic guideway system based on the lubricating medium of electrorheological fluid were studied. A numerical model for fluid dynamics computation was established according to the hydraulic resistance network theory and in comprehensive consideration of correlated parameters and the influence of electric field intensity on the characteristics of electrorheological fluid. The static and dynamic characteristics of the hydrostatic guideway based on the electrorheological fluid with different working variables were studied by using the fluent dynamic mesh technology and the perturbation theory. The results show that the static and dynamic characteristics of the hydrostatic guideway system can be influenced by the electric filed intensity, load factor and initial pressure ratio. The system flow rate will decrease and the damping coefficient will increase with the system stiffness remaining unchanged, when the electric filed intensity is enhanced under the condition of unchanged initially designed parameters and load factors. The resonance phenomenon of the guideway system may occur when the cutting force frequency is changed. The resonance vibration can be suppressed with the dynamic stiffness being increased significantly by enhancing the electric filed intensity. The simulation results are in agreement with the experimental results, which provides a theoretical reference to the design and application of hydrostatic guideways.

hydrostatic guideway; electrorheological fluid; electric field intensity; static and dynamic characteristics; dynamic mesh; correlated parameters

國家自然科學基金資助項目(51175518)；湖南省科技型中小企業(yè)技術(shù)創(chuàng)新基金(12C26214305029)

2015-07-28 修改稿收到日期：2015-11-02

胡均平 男，博士，教授， 1965年6月生

劉成沛 男，博士生， 1987年6月生 E-mail：liuchengpei2012@163.com

TH137

A

10.13465/j.cnki.jvs.2016.22.004