郭鐵能，李 玲，蔡力鋼，劉志峰，趙永勝，楊 坤

(北京工業(yè)大學(xué) 機械工程與應(yīng)用電子技術(shù)學(xué)院，北京 100124)

基于頻響函數(shù)辨識機械結(jié)合部動態(tài)參數(shù)的研究

郭鐵能，李 玲，蔡力鋼，劉志峰，趙永勝，楊 坤

(北京工業(yè)大學(xué) 機械工程與應(yīng)用電子技術(shù)學(xué)院，北京 100124)

進行機械結(jié)構(gòu)的動力分析與動態(tài)優(yōu)化設(shè)計時，為獲得比較準確的結(jié)構(gòu)邊界條件，結(jié)合部參數(shù)辨識成為其中一項關(guān)鍵技術(shù)。結(jié)合工程實際，應(yīng)用頻響函數(shù)辨識結(jié)合部剛度與阻尼參數(shù)的方法進行參數(shù)識別，對能否獲取完備頻響函數(shù)進行了研究，并作出相應(yīng)的求解策略。該方法避開了對頻響函數(shù)直接求逆，運用最小二乘原理將矛盾方程轉(zhuǎn)化為定解方程，保證了數(shù)值計算的穩(wěn)定。一般情況下，結(jié)合部處絕對阻尼值要遠遠小于絕對剛度值，為了更準確地辨識出阻尼值，需采用對阻尼進行二次辨識。該方法具有一定的工程實用價值，而且辨識過程簡單。算例證實了該方法具有很高的辨識精度。

結(jié)合部;頻響函數(shù);參數(shù)辨識

對結(jié)構(gòu)進行動態(tài)分析時，結(jié)合部的影響不可忽略。許多研究表明:結(jié)合部處的總?cè)岫日紮C床總?cè)岫鹊?0%［1］;機床上出現(xiàn)振動問題有60%以上是源自結(jié)合部［2］;結(jié)合部的阻尼占機床總阻尼的 90%以上［3］。這些結(jié)論說明，建立準確的結(jié)構(gòu)動力學(xué)模型關(guān)鍵在于能否獲得合理的結(jié)合部動力學(xué)參數(shù)。然而，影響結(jié)合部動態(tài)特性的因素很多，且比較復(fù)雜，很多學(xué)者已開展了大量的實驗和理論研究，但仍有很多問題沒有得到根本性的解決，成為制約整機結(jié)構(gòu)動靜態(tài)特性建模、工程分析計算與優(yōu)化的一個瓶頸。所以，尋找一種比較合理的結(jié)合部動態(tài)參數(shù)辨識方法，以獲得比較準確的結(jié)構(gòu)邊界條件，是進一步發(fā)展結(jié)構(gòu)動力分析與動態(tài)優(yōu)化設(shè)計的先決條件，是機械動力學(xué)由單件分析走向整機分析所要解決的關(guān)鍵問題。

近年來，國內(nèi)外學(xué)者已提出了許多辨識方法［4－12］，總體可分為兩種，一種是基于模型的辨識方法［4－7］，通過研究實驗數(shù)據(jù)和通過有限元法建立的理論模型相結(jié)合去辨識結(jié)合部參數(shù)。但是，結(jié)構(gòu)的阻尼特性一般很難通過有限元方法準確獲得，這為結(jié)合部的辨識帶來了難度。另一種是基于純實驗的辨識方法［8－12］，通過實驗數(shù)據(jù)去辨識結(jié)合部參數(shù)。這種方法最大的缺陷是，測量時隨機噪聲不可避免。隨機噪聲的存在，使得結(jié)合部的辨識變得更加困難。

以國內(nèi)某龍頭企業(yè)承擔(dān)的國家高檔數(shù)控機床與基礎(chǔ)制造裝備重大專項產(chǎn)品為研究與應(yīng)用對象，與國家高檔數(shù)控機床重大專項產(chǎn)品對接，在綜合前人研究成果的基礎(chǔ)上，假設(shè)結(jié)合部一端為固定剛體(圖1)，采用頻響函數(shù)辨識機械結(jié)合部剛度與阻尼參數(shù)的方法，對實驗中能否獲取完備頻響函數(shù)進行了研究。在辨識過程中避開了對頻響函數(shù)直接求逆，保證了數(shù)值計算的穩(wěn)定，其結(jié)果表明該方法具有很高的辨識精度。

1 識別原理

1.1 結(jié)合部動態(tài)參數(shù)與頻響函數(shù)的關(guān)系

機床中常見結(jié)合部一般可以用如圖1所示的動力學(xué)模型表示，即假定結(jié)構(gòu)一端被完全固定，另一端A自由，結(jié)合面處采用線性彈簧和阻尼器來模擬其特性。辨識的主要任務(wù)是要確定圖中剛度與阻尼參數(shù)。

子結(jié)構(gòu)?A在未聯(lián)接的情況下，動力學(xué)方程可表示為:

圖1 結(jié)合部模型Fig.1 Structure connected with joints

式中:M為結(jié)構(gòu)A的質(zhì)量矩陣;C為結(jié)構(gòu)A的阻尼矩陣;K為結(jié)構(gòu)A的剛度矩陣;F為結(jié)構(gòu)A的力向量陣;Xa為結(jié)構(gòu)A與支承部上假定的彈簧、阻尼器相聯(lián)的運動坐標(a=1，2，…，n)，即聯(lián)接坐標;

XA為結(jié)構(gòu)A除去Xa以外的其余運動坐標，即非聯(lián)接坐標。

這時，結(jié)構(gòu)A在未聯(lián)接狀態(tài)下的動剛度Z可表示:

式中H為結(jié)構(gòu)A在自由狀態(tài)下的頻響函數(shù)矩陣。

當結(jié)構(gòu)A通過線性彈簧和阻尼聯(lián)接于某一剛體后，其系統(tǒng)的動剛度～Z為:

由式(2)與式(3)得:

式(4)表達了待識別參數(shù)P與頻響函數(shù)的數(shù)學(xué)關(guān)系。顯然，按上式識別不僅需要全部的H和，而且需要對H和矩陣求逆。

為避免直接對矩陣求逆，對式(4)兩邊分別左乘～H、右乘H有:

下面根據(jù)式(5)建立辨識P的基本公式和辨識結(jié)合部剛度與阻尼參數(shù)的基本算法。

1.2 基本辨識表達式

假定圖1中被聯(lián)結(jié)構(gòu)在未聯(lián)與已聯(lián)狀態(tài)下，其頻響函數(shù)H和均由獲得。因測試需要，它要求被聯(lián)結(jié)構(gòu)能從原系統(tǒng)中分離出來，對于不可分離的被聯(lián)結(jié)構(gòu)必須在組裝前對其做未聯(lián)狀態(tài)下的激振實驗，以便得到矩陣H。

將式(5)按分塊矩陣展開并完成矩陣相乘

取上式任一分塊矩陣對應(yīng)相等，均可得到求解P的關(guān)系式，即:

以上四式即為基本辨識公式。

1.3 基本辨識方法

在實際現(xiàn)場測試過程中，由于機械系統(tǒng)的配置，整體工藝的要求等原因，有時不是所有的頻響函數(shù)都可以通過測試獲得。為了使問題變得更加實際可解，本文將對能否獲取完備頻響函數(shù)分別進行了討論，并作出相應(yīng)的求解策略。

1.3.1 基于完備頻響函數(shù)辨識方法

假定被聯(lián)結(jié)構(gòu)處于未聯(lián)接狀態(tài)，聯(lián)接前后的頻響函數(shù)矩陣均通過測試獲得。由式(7)得:

式中右邊第l行第m列的元素具有如下形式:

注意，由于頻響函數(shù)中，原點頻響函數(shù)比交叉點頻響函數(shù)可信度高，故取式(10)對角線元素對應(yīng)相等，可獲一組求解pi(i=1，2，…，n)的線性方程組

注意到上式是一組復(fù)數(shù)域線性方程組，設(shè)ER、EI分別為E的實、虛部，ER、EI分別為F的實、虛部，又已知P的實虛部分別為K'J與jωC'J，按實、虛部對應(yīng)相等的原則，上式又可表示為:

式(12)即為一組求解K'J、C'J的實線性方程組。

同樣，因頻響函數(shù)采樣頻率ω可以很多，故求解K'J、C'J的方程也很多，設(shè)ω有m個，則由式(12)有:

上式為一組矛盾方程，由最小二乘原理將其轉(zhuǎn)化為定解方程:

采用高斯列主元素法解上式，即可辨識出K'J、C'J。

一般來說按式(13)辨識的K'J是有相當精度的，但C'J的精度可能要差些，特別是K'J比C'J在量級上大得較多時更是如此，這主要是因為頻響函數(shù)中阻尼的信息往往較剛度信息弱。為了更準確地辨識出阻尼，這里建議對阻尼進行二次辨識。重新回到式(11)，并將已辨識的K'J代入式(11)，得:

然后仿照前面的方法即可再單獨辨識C'J。

1.3.2 基于不完備頻響函數(shù)的辨識方法

首先假定被聯(lián)結(jié)構(gòu)在已聯(lián)與未聯(lián)兩種狀態(tài)下，所有非聯(lián)接坐標點均為可測點，即AA與HAA各元素全可測;另外，注意到(1)的定義，只要我們假定的彈簧、阻尼器不安排在轉(zhuǎn)角等難測自由度上，則未聯(lián)狀態(tài)下與各聯(lián)接坐標相關(guān)的頻響函數(shù)Haa與HAa總是可測的，但已聯(lián)狀態(tài)下頻響函數(shù)aa與Aa(或HaA)中可能存在部分難測元素。這時，對難測的部分矩陣元素采用如下公式計算。

式(14)與式(15)即為補充公式，具體補充方法如下:

2 算例

結(jié)合工程實際，以圖2所示常見模型為分析模型來說明本文方法的實施過程。假定被連接結(jié)構(gòu)在未聯(lián)與已聯(lián)狀態(tài)下均沒有建立有限元模型，其頻響函數(shù)已知，即所有參數(shù)為已知:

圖2 實例系統(tǒng)Fig.2 example of the system

待識別的參數(shù)為 k1、k2和 c1、c2。

則 k1、k2、c1、c2未聯(lián)接時的動力學(xué)方程為:

而系統(tǒng)在已聯(lián)狀態(tài)下的動力學(xué)方程為:

由兩種狀態(tài)下的運動方程，在給定的測試采樣頻率下可計算出系統(tǒng)準確的Z、與H。本例采樣頻率ω取100個，頻率范圍為5 Hz～500 Hz，對未知系統(tǒng)，可根據(jù)聯(lián)接狀態(tài)下系統(tǒng)的幾階主要低頻固有頻率來確定測試頻帶范圍，且建議采樣點固有頻率附近多取些，這有利于反映結(jié)構(gòu)的真實動態(tài)特性。為仿真計算，先對本例系統(tǒng)H、的理論值添加一定水平的隨機誤差，將其作為頻響函數(shù)，本例添加的誤差水平對原點hii、ii與跨點hijij(i≠j)分別控制在5%與6%以內(nèi)，這比較符合誤差分布。同時，對原準確的系統(tǒng)物理參數(shù)添加5%左右的隨機誤差，再形成動剛度矩陣Z。在進行阻尼辨識中，運用了二次辨識的方法，辨識結(jié)果如表1誤差在8.7%以內(nèi)，表明采用該方法是行之有效的。

表1 仿真實例辨識結(jié)果Tab.1 Identification results of simulation example

3 結(jié)論

(1)本文應(yīng)用了基于頻響函數(shù)辨識結(jié)合部剛度與阻尼的方法，對能否獲取完備頻響函數(shù)分別進行了研究，并作出相應(yīng)的求解策略。辨識過程避免了對頻響函數(shù)矩陣求逆運算，且利用的是較低維矩陣關(guān)系，從而提高了辨識精度。

(2)文章是以結(jié)合部一端為固定剛體為前提，而進行的辨識分析，故在采用該方法時，應(yīng)將一端固定。但如果兩端均不能固定，則應(yīng)分別建立兩部件在自由狀態(tài)的有限元模型，或分別對兩部件在自由狀態(tài)下測試其頻響函數(shù)，再在整體狀態(tài)下測試頻響函數(shù)，仿照本文原理，可導(dǎo)出類似的辨識公式。

(3)采用本文的辨識方法可進行系列化實驗研究，以獲得可供實用的經(jīng)驗公式與圖表。并將辨識得到的結(jié)合部等效動力學(xué)參數(shù)，一方面可以判定結(jié)合部等效剛度和阻尼的特性及其隨外加載荷作用的變化規(guī)律;另一方面可獲得典型工況下結(jié)合部完整的動力學(xué)模型，為動力分析、動力修改及動態(tài)優(yōu)化設(shè)計奠定基礎(chǔ)，為產(chǎn)品的快速數(shù)字化設(shè)計提供理論與技術(shù)支持。

［1］Connolly R，Thomley R H.The Significance of Joint on the Overall Deflection of Machine Tool Structures 6thMTDR Conf，1965.

［2］Burdekin M，Baek N A.Cowley.Analysis of the local deformation in machine joints［J］.Mech.Eng.Sci，1979，34(3):15－17.

［3］Vafaei S，Rahnejat H，Aini R.Vibration monitoring of high speed spindlesusingspectralanalysistechniques［J］Machine Tools＆ Manufacture，2002，42:1223－1234.

［4］Nalitolela N G.A.Mass or Stiffness Addition Technique for Structure Parameter Updating［J］.The International Journal of Analytical and Experimental Modal Analysis.1992，7(3):157－168.

［5］HjelmstadK， MutualD. ResidualEnergyMethodfor Parameter Estimation in Structures［J］.Journal of Structural Engineering，1992，118(1):223 －242.

［6］李潤方，韓 西，等.齒輪傳動系統(tǒng)結(jié)合部動力學(xué)參數(shù)識別［J］.中國機械工程.2001，12(12):1333－1335.

［7］廖伯瑜，周新民，尹志宏.現(xiàn)代機械動力學(xué)及其工程應(yīng)用［M］.機械工業(yè)出版社.2004.

［8］Liu W，Ewins D J.Substructure synthesis via elastic media［J］.Journal of Sound and Vibration，2002，257(2):361－379.

［9］Damjan C，Miha B.Identification of the dynamic properties of joints using frequency-response functions［J］.Journal of Sound and Vibration，2008，317(1－2):158－174.

［10］ Mascardelli B A.Park S S，F(xiàn)reiheit T.Substructure coupling of micro end mills to aid in the suppression of chatter［J］.Journal of Manufacturing Science and Engineering，Transactionsofthe ASME， 2008， 130(l)0110101－01101012.

［11］ Ren Y，Beards C F.Identification of Effective'Linear Joints Using Coupling and Joint Identification Techniques［J］.Journal of Vibration and Acoustics.1998，120，331 －338.

［12］ Wang J H.Identification of Parameters of Structural Joints by Use of Noise Contaminated FRFS［J］.Journal of Sound and Vibration，1990，142:261 －277.

Identifying mechanical joint dynamic parameters based on measured frequency response functions

GUO Tie-neng，LI Ling，CAI Li-gang，LIU Zhi-feng，ZHAO Yong-sheng，Yang Kun

(College of Mechanical Engineering and Applied Electronics Technology，Beijing University of Technology，Beijing，100124，China)

In order to obtain a more accurate boundary condition of structure in dynamic analysis and optimization design，the joint parameter identification is one of the key technologs.In combination with engineering practice，the method based on frequency response function(FRF)to identify the characteristic parameters of the joint was proposed.Whether the complete FRF can be obtained or not was studied and the corresponding solving strategy was presented.The method avoids the matrix inversion directly，and transforms contradictory equation to equation with definite solution by using least-square principle，which ensures the stability of numerical calculation.Considering the damping information is inferior to stiffness information in identification process，in order to identify the damping parameter more accurately，a successive identification method was further proposed.Example analysis confirms that the method is of high identification accuracy.

joint;frequency response function;parameters identification

TH113;TH123

A

國家重大專項(2009ZX04014－033);國家自然科學(xué)基金(50905004)

2010－03－15 修改稿收到日期:2010－04－12

郭鐵能 男，博士，講師，1975年5月生

李 玲 男，博士生，1981年7月生