一種正壓力隨時間變化的微滑移干摩擦模型

上官博 徐自力 肖俊峰 高松 于飛龍

摘要： 在Csaba摩擦模型的基礎(chǔ)上發(fā)展了一種能夠考慮接觸面摩擦運(yùn)動過程中正壓力隨時間變化的微滑移干摩擦模型。模型將局部微滑移運(yùn)動與整體宏滑移運(yùn)動進(jìn)行了結(jié)合，采用在每個時間迭代上逐步分析判斷摩擦微元運(yùn)動狀態(tài)的方法來描述摩擦運(yùn)動過程中的黏滯、局部微動滑移、整體宏滑移以及分離狀態(tài)，能夠克服Csaba摩擦模型只能描述局部微滑移摩擦力、并且正壓力分布不變的不足。對單諧波和多諧波外力作用下無分離狀態(tài)及存在分離狀態(tài)的干摩擦力進(jìn)行了模擬。計算了凸肩阻尼結(jié)構(gòu)模擬葉片在復(fù)雜激勵下的非線性振動響應(yīng)，研究了系統(tǒng)參數(shù)對干摩擦阻尼減振效果的影響，為工程中研究具有干摩擦阻尼結(jié)構(gòu)葉片的接觸面摩擦力提供了理論依據(jù)。關(guān)鍵詞： 轉(zhuǎn)子動力學(xué)；葉片；干摩擦力； 微滑移模型； 變正壓力

中圖分類號： V231.96 ； TK263.3文獻(xiàn)標(biāo)志碼： A文章編號： 1004-4523（2016）03-0444-08

DOI：10.16385/j.cnki.issn.10044523.2016.03.010

引言

葉片是汽輪機(jī)、燃?xì)廨啓C(jī)等透平機(jī)械實現(xiàn)能量轉(zhuǎn)換的關(guān)鍵部件，常因振動產(chǎn)生疲勞破壞而導(dǎo)致機(jī)組故障。由于附加干摩擦阻尼結(jié)構(gòu)可以有效降低葉片的振動水平，減少葉片振動疲勞破壞，干摩擦阻尼結(jié)構(gòu)在透平機(jī)械中已得到了廣泛的應(yīng)用，因此其摩擦減振特性的研究具有十分重要的理論意義和應(yīng)用價值。常見的葉片干摩擦阻尼結(jié)構(gòu)形式包括：圍帶、拉筋、凸肩、燕尾型葉根等等。由于干摩擦接觸往往是較為復(fù)雜的非線性接觸，葉片摩擦阻尼結(jié)構(gòu)設(shè)計的困難主要在于如何描述接觸界面的摩擦運(yùn)動狀態(tài)，以及更精確的計算接觸面的非線性干摩擦力。

宏滑移模型（macroslip）和微滑移模型（microslip）是干摩擦阻尼結(jié)構(gòu)研究中最常用到的兩種摩擦力模型。宏滑移模型又稱為整體滑移模型，假設(shè)接觸面上所有點(diǎn)的變形是均勻的，接觸面上所有接觸點(diǎn)同時滑動或黏滯，接觸面上的運(yùn)動狀況可通過一個點(diǎn)進(jìn)行描述。在法向接觸正壓力較小或者接觸面積較小的情況下，應(yīng)用這種原理簡單的模型對接觸面進(jìn)行簡化也可以得到較理想的結(jié)果[14]。為解決正壓力較大時整體滑移模型的局限性，準(zhǔn)確描述接觸面上應(yīng)力分布不均，采用多個接觸點(diǎn)來描述接觸面摩擦接觸特性的微動滑移模型（部分滑移模型）[510]，逐漸被學(xué)者們采用。微動滑移模型將接觸界面離散成多個接觸點(diǎn)對，分別求解各個接觸點(diǎn)對的運(yùn)動狀態(tài)以及由其產(chǎn)生的摩擦約束力，較為精細(xì)地描述摩擦接觸面的摩擦狀態(tài)。何尚文等人[11]采用具有彈性剪切層的微滑移模型對BG型葉片緣板阻尼器的減振特性進(jìn)行了研究。徐自力[1213]等人采用Csaba微滑移摩擦模型研究了結(jié)構(gòu)參數(shù)對葉片減振效應(yīng)影響，研究結(jié)果為阻尼器設(shè)計提供了理論依據(jù)。

本文在Csaba微滑移摩擦模型的基礎(chǔ)上發(fā)展了一種能夠考慮接觸面摩擦運(yùn)動過程中正壓力隨時間變化的微滑移干摩擦模型。該模型通過在每個時間迭代上逐步分析摩擦接觸面離散微元運(yùn)動狀態(tài)的方法得到接觸面的摩擦接觸力，可以描述摩擦表面的黏滯滑動共存狀態(tài)和法向接觸正壓力的非均勻分布。模型還將局部微滑移運(yùn)動與整體宏滑移運(yùn)動有效結(jié)合，克服了原模型只能描述局部微滑移摩擦力的不足。模擬了單諧波和多諧波外力作用下無分離狀態(tài)及存在分離狀態(tài)的干摩擦力遲滯回線，并計算了凸肩阻尼結(jié)構(gòu)模擬葉片在復(fù)雜激勵下的非線性振動響應(yīng)。

1變正壓力微滑移干摩擦模型[*2]1.1模型描述兩接觸面間的干摩擦阻尼通過圖1所示的一個壓在剛性平面上的矩形板來模擬，矩形板的彈性模量為E，橫截面積為A，長度為l，接觸面的摩擦系數(shù)為μ。矩形板右端作用著外力F，右端位移為u。矩形板上端作用有隨時間變化的法向分布載荷，將其表示為二元函數(shù)q（x，t），t為時間變量。為了更精細(xì)地求解摩擦運(yùn)動過程中矩形板的干摩擦力，采用數(shù)值方法將矩形劃分為若干離散單元，單元個數(shù)為n，節(jié)點(diǎn)數(shù)為n+1。

4結(jié)論

（1） 本文發(fā)展的變正壓力微滑移干摩擦模型，是在Csaba摩擦模型的基礎(chǔ)上考慮了接觸面法向正壓力分布隨時間變化對接觸面干摩擦力的影響，不僅克服了原模型中正壓力分布始終保持不變的不足，更突破了原模型只能描述小運(yùn)動位移的限制，將其擴(kuò)展到宏滑移階段。該模型可以描述摩擦接觸表面的黏滯狀態(tài)、微動滑移狀態(tài)、整體滑移狀態(tài)以及分離狀態(tài)，適用于變正壓力的復(fù)雜摩擦運(yùn)動部件。

（2） 通過對單諧波和多諧波外力作用下接觸面干摩擦力遲滯回線的模擬，發(fā)現(xiàn)本文發(fā)展的變正壓力微滑移干摩擦模型能夠彌補(bǔ)YANG模型無法考慮微滑移摩擦所消耗的那部分接觸面阻尼，這一結(jié)果對接觸面摩擦阻尼耗能的評估至關(guān)重要。

（3） 研究發(fā)現(xiàn)存在一個最優(yōu)初始正壓力區(qū)間，在該區(qū)間內(nèi)葉片的減振效果較好，該結(jié)論可以指導(dǎo)制造和安裝過程中的裕度設(shè)定。外激勵增加導(dǎo)致阻尼器剛度效果減弱從而使系統(tǒng)共振頻率減小。

參考文獻(xiàn)：

[1]Petrov E P， Ewins D J. Advanced modeling of underplatform friction dampers for analysis of bladed disk vibration[J]. Journal of Turbomachinery， 2007， 129： 143—150.

[2]Yang B D， Menq C H. Characterization of contact kinematics and application to the design of wedge dampers in turbomachinery blading： part 1-stickslip contact kinematics[J]. Journal of Engineering for Gas Turbines and Power， 1998， 120： 410—417.

[3]徐自力， 上官博， 吳其林，等. 松裝葉片燕尾型葉根干摩擦力模型及振動響應(yīng)分析[J]. 西安交通大學(xué)學(xué)報， 2009， 43 （7）： 1—6.

Xu Zili， Shangguan Bo， Wu Qilin， et al. Dry Friction force model and dynamic response analysis of blades with loosely assembled dovetail attachment[J]. Journal of Xian Jiaotong University， 2009， 43 （7）： 1—6.

[4]Chen J J， Yang B D， Menq C H. Periodic forced response of structures having threedimensional frictional constraints[J]. Journal of Sound and Vibration， 2000， 229 （4）： 775—792.

[5]Iwan W D. On a class of models for the yielding behavior of continuous and composite systerms[J]. ASME Journal of Applied Mechanics， 1967， 89： 612—617.

[6]Menq C H， Bielak J， Griffin J H. The influence of microslip on vibratory response， part I： a new microslip model[J]. Journal of Sound and Vibration， 1986， 107（2）： 279—293.

[7]Menq C H， Bielak J， Griffin J H. The influence of microslip on vibratory response， part II： a comparision with experimental results[J]. Journal of Sound and Vibration， 1986， 107（2）： 295—307.

[8]Csaba G. Modeling microslip friction damping and its influence on turbine blade vibrations[D]. Linkoping， Sweden： Linkoping University， 1998.

[9]Csaba G. Forced response analysis in time and frequency domains of a tuned bladed disk with friction dampers[J]. Journal of Sound and Vibration， 1998， 214 （3）： 395—412.

[10]Sextro W. The caculation of the forced response of shrouded blades with friction contacts and its experimental verification[C]. ASME Paper 2000GT540， Munich， Germany， 2000.

[11]何尚文， 任興民， 秦衛(wèi)陽， 等. 基于微滑移模型的BG型葉片緣板阻尼器減振特性研究[J]. 振動與沖擊， 2012， 31（1）：44—48.

He Shangwen， Ren Xingmin， Qin Weiyang， et al. Reduced vibration characteristics of a BG platform damper based on a microslip model[J]. Journal of Vibration and Shock， 2012， 31（1）： 44—48.

[12]徐自力， 常東鋒， 上官博. 微滑移離散模型及在干摩擦阻尼葉片振動分析中的應(yīng)用[J]. 機(jī)械科學(xué)與技術(shù)， 2007， 26 （10）： 1304—1307.

Xu Zili， Chang Dongfeng， Shangguan Bo. Onebar microslip discrete model and its application to vibration analysis of blade with dry friction damper[J]. Mechanical Science and Technology for Aerospace Engineering， 2007， 26 （10）： 1304—1307.

[13]徐自力， 常東鋒， 劉雅琳. 基于微滑移解析模型的干摩擦阻尼葉片穩(wěn)態(tài)響應(yīng)分析[J]. 振動工程學(xué)報， 2008， 21 （5）： 505—510.

Xu Zili， Chang Dongfeng， Liu Yalin. Forced response analysis of blade system with dry friction damper using onebar microslip analytic model. Journal of Vibration Engineering， 2008， 21 （5）： 505—510.

[14]Yang B D， Chu M L， Menq C H. Stickslipseparation analysis and nonlinear stiffness and damping characterization of friction contacts having variable normal load[J]. Journal of Sound and Vibration， 1998， 210（2）：461—481.