焦衛(wèi)東， 蔣永華， 施繼忠， 王曉燕

(浙江師范大學(xué)工學(xué)院， 浙江 金華 321004)

引 言

現(xiàn)有的裂紋模型可以分為開裂紋與呼吸裂紋兩大類，開裂紋模型適用于裂紋面受恒定方向拉應(yīng)力作用的情形，其效應(yīng)是導(dǎo)致轉(zhuǎn)軸剛度的局部定值削減[1]；呼吸裂紋模型的裂紋面承受拉、壓應(yīng)力的交替作用，裂紋周期性開/合(或呼吸)，激起轉(zhuǎn)軸剛度的周期性變化[2]。Darpe 等推導(dǎo)了包含橫向裂紋(即裂紋面同時垂直于轉(zhuǎn)軸軸線及其基面)轉(zhuǎn)子的柔度參數(shù)與剛度矩陣，并研究了在裂紋非線性呼吸行為作用下轉(zhuǎn)子的縱向、彎曲與扭轉(zhuǎn)耦合振動特性[3]。在已有的研究基礎(chǔ)上，Darpe 應(yīng)用應(yīng)變能理論將橫向裂紋推廣至更一般的斜裂紋情形，并與橫向裂紋進(jìn)行了對比。相比于橫向裂紋，斜裂紋會導(dǎo)致剛度矩陣中不同參數(shù)的更多耦合，進(jìn)而造成轉(zhuǎn)子彎曲、扭轉(zhuǎn)甚至是縱向振動的耦合[4]。Bachschmid 等指出，通常情況下轉(zhuǎn)子裂紋面垂直于轉(zhuǎn)軸軸線方向(即橫向)，但在大扭矩和強彎矩載荷作用下裂紋會沿著螺旋方向擴展，從而形成螺旋裂紋或斜裂紋。并以精細(xì)的三維非線性模型為基礎(chǔ)，建立起更為簡單的裂紋轉(zhuǎn)子一維有限元模型，據(jù)此對橫裂紋與斜裂紋的呼吸機理、彎扭聯(lián)合作用下裂紋的擴展行為進(jìn)行數(shù)值仿真與實驗驗證[5]。Han 等認(rèn)為齒輪嚙合力導(dǎo)致的大扭矩往往引發(fā)裂紋的斜向擴展，因此對含有斜裂紋的齒輪-轉(zhuǎn)子-軸承系統(tǒng)的動力學(xué)特性進(jìn)行了數(shù)值仿真研究[6]。Ebrahimi 等建立了橫向開裂紋連續(xù)彎曲振動模型，運動方程基于漢密爾頓原理構(gòu)建，采用改進(jìn)的伽遼金方法求解，所得的響應(yīng)結(jié)果與有限元法呈現(xiàn)出良好的一致性[7]。國內(nèi)對裂紋轉(zhuǎn)子動力學(xué)也開展了一定的研究。閆明等采用有限元方法和最大周向應(yīng)力準(zhǔn)則，分析計算了熱疲勞作用下斜裂紋的應(yīng)力強度因子和開裂角的周期性變化規(guī)律[8]。劉長利等對橫裂紋與斜裂紋轉(zhuǎn)子的變剛度特性進(jìn)行了對比[9]，并進(jìn)一步研究了包含雙盤、雙呼吸型裂紋的轉(zhuǎn)子系統(tǒng)的一些典型的非線性動力學(xué)特性[10]。夏恒恒等分析了斜裂紋轉(zhuǎn)子在裂紋全部張開狀態(tài)下剛度的影響因素[11]。

現(xiàn)有的絕大部分研究主要是基于斷裂力學(xué)中的應(yīng)變能理論，并采用 Euler 梁或Timoshenko 梁單元進(jìn)行建模[3-4, 6, 9-11]?；诹簡卧Ｅc應(yīng)變能理論的有限元法物理意義明確、理論基礎(chǔ)扎實，因此獲得廣泛應(yīng)用。但是，這些研究絕大部分針對的是簡單的橫向裂紋或裂紋面傾斜于轉(zhuǎn)軸軸線而垂直于轉(zhuǎn)軸基面的斜裂紋，這類裂紋稱為橫-斜或直-斜裂紋更為合適，很少見到對含有空間任意斜裂紋(即裂紋面與轉(zhuǎn)軸軸線和轉(zhuǎn)軸基面均不垂直)轉(zhuǎn)子的動力學(xué)特性進(jìn)行研究。本文著眼于探究空間任意斜裂紋轉(zhuǎn)子的多自由度耦合振動機理，采用基于梁單元建模與應(yīng)變能理論的有限元法，對轉(zhuǎn)子旋轉(zhuǎn)過程中由于裂紋的非線性呼吸行為導(dǎo)致的剛度變化以及不同方向剛度參數(shù)的交叉耦合特性進(jìn)行研究，并系統(tǒng)分析裂紋方位角以及裂紋深度等參數(shù)對轉(zhuǎn)子剛度特性的影響。裂紋轉(zhuǎn)子建模采用考慮剪切變形效應(yīng)的 Timoshenko 梁單元，并考慮縱向、彎曲與扭轉(zhuǎn)全部六個方向的自由度。

1 空間任意斜裂紋轉(zhuǎn)子剛度計算

1.1 裂紋轉(zhuǎn)子建模

轉(zhuǎn)軸長度為L= 0.7 m，直徑為D= 0.015 m，轉(zhuǎn)子盤質(zhì)量為m= 1 kg。裂紋位于轉(zhuǎn)子盤右端靠近盤位置。轉(zhuǎn)子-軸承系統(tǒng)采用Timoshenko梁單元進(jìn)行建模，共劃分為14個單元，包含15個節(jié)點，每個節(jié)點考慮了縱向、彎曲和扭轉(zhuǎn)所有6個自由度，如圖1(a)所示。從便于描述裂紋動力學(xué)效應(yīng)的角度，對于裂紋轉(zhuǎn)子軸段進(jìn)行單獨建模。

圖1 裂紋轉(zhuǎn)子的有限元模型及其坐標(biāo)系統(tǒng)(橫裂紋：xyz；橫-斜裂紋：x′y′z′；空間任意斜裂紋：x″y″z″)Fig.1 Finite element models of cracked rotors and their coordinate systems (transverse: xyz, transverse-slant: x′y′z′ and arbitrary spatial slant: x″y″z″)

為了具體說明空間任意斜裂紋、橫-斜裂紋以及橫裂紋相互之間的空間變換關(guān)系，圖1(b)～(d)對比給出了各自的坐標(biāo)系統(tǒng)，其中圖1(e)中xoy面為轉(zhuǎn)子的基面。在圖1(e)中，進(jìn)一步描述了各坐標(biāo)系統(tǒng)的空間變換關(guān)系?？梢钥吹?，由最基本的橫裂紋(xyz)出發(fā)，其裂紋面繞z坐標(biāo)軸旋轉(zhuǎn)一定角度(θ1)，即得橫-斜裂紋(x′y′z′)；橫-斜裂紋面繞其坐標(biāo)軸y′ 再旋轉(zhuǎn)一定角度(θ2)，即得空間任意斜裂紋，其坐標(biāo)系統(tǒng)為x″y″z″。厘清三類裂紋坐標(biāo)系統(tǒng)的空間變換關(guān)系，是裂紋有限元建模及動力學(xué)參數(shù)計算的必要條件。

1.2 裂紋單元柔度參數(shù)及剛度矩陣計算

設(shè)轉(zhuǎn)軸半徑為R，梁元長度為l。裂紋中心距離單元左端為x，裂紋深度為a。 單元承受剪力P2，P3和P8，P9，彎矩P5，P6和P11，P12，軸向力P1和P7以及扭矩P4和P10的作用，如圖1所示。令ui和Pi分別為沿著第i個坐標(biāo)方向的節(jié)點位移與節(jié)點力。根據(jù) Castinglianos 定理，裂紋單元的柔度參數(shù)表達(dá)為

(1)

式中U0為無裂紋轉(zhuǎn)軸單元的應(yīng)變能，Uc為裂紋導(dǎo)致的外加應(yīng)變能。

裂紋單元節(jié)點位移ui為

(2)

(3)

式中A= πR2為轉(zhuǎn)子軸橫截面積，E為楊氏彈性模量，G為剛性模量，I為轉(zhuǎn)子軸截面面積矩，I0為截面極慣矩，αs為Timoshenko 梁剪切系數(shù)。

由裂紋導(dǎo)致的外加應(yīng)變能計算式為

(4)

式(4)中的各個 SIF 計算式如下：

(1)張開模式(Opening Mode)

(5)

(6)

(7)

(8)

(9)

(10)

(2)滑移模式(Sliding Mode)

(11)

(12)

(13)

(14)

(15)

(16)

(3)剪開模式(Tearing Mode)

(17)

(18)

(19)

(20)

(21)

(22)

圖2 包含裂紋閉合線的 AS 型裂紋截面Fig.2 Cross-section of AS type crack including CCL

聯(lián)立以上各式，可得各柔度參數(shù)如下：

(23)

(24)

(25)

8Rbβsin2θ1sin2θ1sin2θ2sin2θ2+

(26)

(27)

(28)

(29)

(30)

4βsin4θ1sin2θ2sin2θ2)F1F2+(Rbsin2θ1sin2θ1sin22θ2+

(31)

(32)

(33)

(34)

4βsin3θ1sin22θ2)F1F2+(Rbsinθ1sin2θ1sin4θ2+

(35)

(36)

(37)

4βsin2θ1sin2θ1sin2θ2sin2θ2)F1F+(Rbsin22θ1sin22θ2+

4βsin5θ1sin2θ2sin2θ2)F1F2+(Rbsin3θ1sin2θ1sin22θ2+

(38)

(39)

(40)

4βsin5θ1sin2θ2sin2θ2)F1F2+(2Rbsin3θ1sin2θ1sin22θ2+

(41)

(42)

(43)

如圖1(b)，考慮裂紋單元各節(jié)點位移qi，i= 1～12 的靜平衡條件，有 {q1-12}T=T{q1-6}T，變換矩陣T參見文獻(xiàn) [3-4]。從而，裂紋單元的剛度矩陣為Kc=TG-1TT，其中G= {gij}，i,j=1～6，各柔度參數(shù)gij可以根據(jù)式(23)～(43)求得。

2 數(shù)值仿真結(jié)果

隨著轉(zhuǎn)子軸旋轉(zhuǎn)，裂紋發(fā)生交替性開/合，導(dǎo)致裂紋轉(zhuǎn)子剛度的周期性變化。為了仿真裂紋轉(zhuǎn)子的變剛度特性，在一個旋轉(zhuǎn)周期內(nèi)(360°)將裂紋邊均分為100個點，每個點確定一個 CCL 位置，由此確定不同的積分限。從而，可將式(29)計算所得的剛度參數(shù)變化視為是 CCL 位置的函數(shù)。

裂紋轉(zhuǎn)子的變剛度特性不僅取決于裂紋的類型，而且與裂紋的方位角(θ1與/或θ2)和裂紋深度a有關(guān)。接下來，將對此作詳細(xì)研究。

2.1 不同類型裂紋變剛度特性的對比

圖3 三類不同裂紋的剛度參數(shù)對比(橫裂紋：T；橫-斜裂紋：TS；空間任意斜裂紋：AS)Fig.3 Comparison on stiffness coefficients of three different types of crack (transverse: T； transverse-slant: TS and arbitrary spatial slant: AS)

由圖3(a)可見，在轉(zhuǎn)子整個旋轉(zhuǎn)周期中，空間任意斜裂紋(AS)轉(zhuǎn)子相比于橫裂紋和橫-斜裂紋(T 與 TS)轉(zhuǎn)子，具有最大的軸向與垂直剪切剛度(k11與k33)。在旋轉(zhuǎn)初期，AS 型裂紋轉(zhuǎn)子的扭轉(zhuǎn)剛度(k44)略小于 TS 型裂紋轉(zhuǎn)子；在旋轉(zhuǎn)末期，AS 型裂紋轉(zhuǎn)子的彎曲剛度(k55與k66)略小于 T 型裂紋轉(zhuǎn)子。不過，從整體上來看其剛度差不多也是最大的。但是，在轉(zhuǎn)子旋轉(zhuǎn)的大部分過程中 AS 型裂紋轉(zhuǎn)子的水平剪切剛度(k22)明顯小于其他兩類裂紋轉(zhuǎn)子，雖然在旋轉(zhuǎn)末期其剛度有明顯提高。此外，AS 型與 TS 型裂紋轉(zhuǎn)子的扭轉(zhuǎn)剛度 (k44) 很接近，彎曲剛度(k55與k66)的變化趨勢相近，表明在扭轉(zhuǎn)和彎曲方向上裂紋面方位角θ1是影響轉(zhuǎn)子剛度特性的主要因素。而在其他方向上，方位角θ2則發(fā)揮了重要的作用，導(dǎo)致AS型裂紋轉(zhuǎn)子具有與其他兩類裂紋轉(zhuǎn)子明顯不同的、更加復(fù)雜的變剛度特性，尤其在水平剪切方向上(k22)表現(xiàn)得更加明顯，剛度曲線的對稱性被嚴(yán)重破壞，剛度值明顯下降，整體小于其他兩類裂紋轉(zhuǎn)子。

交叉耦合剛度參數(shù)kij，i=j，是造成裂紋轉(zhuǎn)子不同方向振動響應(yīng)相互耦合的根本原因。對于 AS 型裂紋轉(zhuǎn)子，剛度參數(shù)彼此交叉耦合的現(xiàn)象更加顯著。例如，在水平剪切-垂直剪切(k23)、水平剪切-扭轉(zhuǎn)(k24)、水平剪切-垂直彎曲(k25)以及垂直剪切-水平彎曲(k36)等方向，均出現(xiàn)了強烈的耦合現(xiàn)象，如圖3(b)所示。

2.2 方位角對AS型裂紋變剛度特性的影響

保持裂紋面方位角θ1= 45° 不變，θ2在 30° 到 90°之間變化。在不同方位角的共同作用下，AS 型裂紋轉(zhuǎn)子的剛度特性曲線如圖4所示。需要注意的是，當(dāng)θ2= 90°時AS型裂紋退化為TS型裂紋，具體可以參考圖1(e)加以理解。

圖4 方位角(θ2)對裂紋轉(zhuǎn)子變剛度特性的影響Fig.4 Influence of oriented angle (θ2) on stiffness variation characteristics of cracked rotor

圖5 裂紋深度比(ā = )對裂紋轉(zhuǎn)子變剛度特性的影響Fig.5 Influence of crack depth ratio (ā = on stiffness variation characteristics of cracked rotor

由圖4可見，隨著θ2的增大，縱向(k11)、垂直剪切方向(k33)、扭轉(zhuǎn)方向(k44)以及垂直彎曲方向(k55)的剛度值單調(diào)下降。但是，在水平彎曲方向(k66)以及水平剪切方向(k22)，這種趨勢不復(fù)存在，剛度變化曲線彼此出現(xiàn)了明顯交叉，表明 AS 型裂紋轉(zhuǎn)子的裂紋面方位角θ1與θ2之間的交互作用效應(yīng)。顯然，在轉(zhuǎn)子整個旋轉(zhuǎn)周期的不同轉(zhuǎn)角位置，這種效應(yīng)對轉(zhuǎn)子剛度的影響是不同的，因此是動態(tài)的，也是非線性的。

2.3 裂紋深度對 AS 型裂紋變剛度特性的影響

進(jìn)一步地，保持 AS 型裂紋轉(zhuǎn)子的裂紋面方位角θ1= 45°，θ2= 60°不變，裂紋深度比 ā 從0.1變化到0.5，裂紋轉(zhuǎn)子的剛度特性如圖5所示。對應(yīng)于不同裂紋深度的剛度特性曲線，采用與圖4相同的線型按照ā遞增的順序進(jìn)行描述。

由圖5可見，隨著ā的增大，AS 型裂紋轉(zhuǎn)子各方向的剛度參數(shù)kii，i= 1～6 均單調(diào)下降。在轉(zhuǎn)子整個旋轉(zhuǎn)周期的不同轉(zhuǎn)角位置，剛度值是不斷變化的，但是對應(yīng)于不同裂紋深度的轉(zhuǎn)子變剛度特性曲線幾乎具有相同的走向和變化趨勢。這表明，裂紋深度參數(shù)對轉(zhuǎn)子剛度的影響雖然也是動態(tài)的，但是具有明顯的線性特征。

3 結(jié)論與討論

基于材料力學(xué)的應(yīng)力理論以及斷裂力學(xué)的應(yīng)變能理論，采用 Timoshenko 梁單元模型對裂紋單元進(jìn)行單獨建模，并考慮了縱向、彎曲以及扭轉(zhuǎn)所有六個方向的自由度，推導(dǎo)了含有空間任意斜(AS)裂紋轉(zhuǎn)子單元的柔度參數(shù)，進(jìn)而根據(jù)節(jié)點位移的靜平衡條件導(dǎo)出了裂紋單元剛度矩陣。在此基礎(chǔ)上，對比研究了不同類型裂紋(T：橫向、TS：橫-斜與AS：空間任意斜)轉(zhuǎn)子的剛度特性，并進(jìn)一步研究了裂紋面方位角以及裂紋深度參數(shù)對AS型裂紋轉(zhuǎn)子變剛度特性的影響。研究結(jié)果表明，AS型裂紋轉(zhuǎn)子由于受到裂紋面方位角θ1與θ2的交互作用，剛度特性明顯不同于其他兩類裂紋轉(zhuǎn)子。這種交互作用效應(yīng)，不僅體現(xiàn)在剛度參數(shù)的量值上，還體現(xiàn)在剛度曲線的變化趨勢上，并且具有明顯的非線性特征。由此導(dǎo)致的結(jié)果是，不同方向上 AS 型裂紋轉(zhuǎn)子的剛度具有更廣泛、更強烈且更加復(fù)雜的交叉耦合特性。

裂紋轉(zhuǎn)子剛度矩陣中的交叉耦合參數(shù)，可以解釋由裂紋引發(fā)的轉(zhuǎn)子在不同方向上的運動耦合機理。根據(jù)Darpe 的研究，對于 TS 型裂紋轉(zhuǎn)子，其幾乎所有的交叉耦合剛度參數(shù)均大于 T 型裂紋轉(zhuǎn)子，特別是當(dāng)涉及到扭轉(zhuǎn)方向時(例如，k14，k24，k34，k45和k46)。因此，Darpe將扭轉(zhuǎn)與彎曲、縱向振動的耦合即縱-彎-扭耦合振動作為研究重點[4]。對于 AS 型裂紋轉(zhuǎn)子，鑒于其更廣泛、更復(fù)雜的變剛度特性，有必要對轉(zhuǎn)子的縱-彎-扭耦合振動以及縱-彎、彎-彎(包括水平和垂直方向)耦合振動進(jìn)行全面研究，以獲取描述 AS 型裂紋轉(zhuǎn)子振動特性的典型特征，用于轉(zhuǎn)子裂紋故障的辨識。此外，無論是 T，TS 還是 AS 型裂紋轉(zhuǎn)子，其剛度參數(shù)與裂紋深度參數(shù)之間均具有良好的線性關(guān)系，這或許可以為裂紋故障的定量診斷提供一定的理論依據(jù)。

[1] Papadopoulos C A, Dimarogonas A D. Coupled longitudinal and bending vibrations of a rotating shaft with an open crack[J]. Journal of Sound and Vibration, 1987, 117(1): 81—93.

[2] Chen Y S, Zhang H B. Review and prospect on the research of dynamics of complete aero-engine systems[J]. Acta Aeronautica et Astronautica Sinica, 2011, 32(8): 1371—1391.

[3] Darpe A K, Gupta K, Chawla A. Coupled bending, longitudinal and torsional vibrations of a cracked rotor[J]. Journal of Sound and Vibration, 2004, 269(1): 33—60.

[4] Darpe A K. Coupled vibrations of a rotor with slant crack[J]. Journal of Sound and Vibration, 2007, 305(1-2): 172—193.

[5] Bachschmid N, Pennacchi P, Tanzi E. Some remarks on breathing mechanism, on non-linear effects and on slant and helicoidal cracks[J]. Mechanical Systems and Signal Processing, 2008, 22(4): 879—904.

[6] Han Q K, Zhao J S, Chu F L. Dynamic analysis of a geared rotor system considering a slant crack on the shaft[J]. Journal of Sound and Vibration, 2012, 331(26): 5803—5823.

[7] Ebrahimi A, Heydarin M, Behzad M. A continuous vibration theory for rotors with an open edge crack[J]. Journal of Sound and Vibration, 2014, 333(15): 3522—3535.

[8] 閆 明, 張義民, 何雪浤, 等. 熱疲勞斜裂紋應(yīng)力強度因子有限元分析[J]. 東北大學(xué)學(xué)報(自然科學(xué)版), 2011，32(5): 720—723.

Run Ming, Zhang Yi-min, He Xue-hong, et al. Time-varying stiffness characteristics of rotating shaft with transverse and slant crack[J]. Journal of Northeastern University (Natural Science), 2011，32(5): 720—723.

[9] 劉長利, 李 誠, 周邵萍, 等. 直斜裂紋轉(zhuǎn)軸的時變剛度特性研究[J]. 振動與沖擊, 2011，30(3)：165—170.

Liu Chang-li, Li Cheng, Zhou Shao-ping, et al. Time-varying stiffness characteristics of rotating shaft with transverse and slant crack[J]. Journal of Vibration and Shock, 2011，30(3): 165—170.

[10] 劉長利, 周邵萍, 江 君, 等. 雙盤雙呼吸型裂紋轉(zhuǎn)子的非線性動力學(xué)特性[J]. 振動、測試與診斷，2012，32(S1): 136—140.

Liu Chang-li, Zhou Shao-ping, Jiang Jun, et al. Nonlinear dynamics analysis of double-disc rotor with two breathing cracks[J]. Journal of Vibration, Measurement & Diagnosis, 2012，32(S1): 136—140.

[11] 夏恒恒, 李志農(nóng), 肖堯先. 斜裂紋轉(zhuǎn)子全開狀態(tài)的剛度影響因素分析[J]. 南昌航空大學(xué)學(xué)報(自然科學(xué)版), 2014，28(3): 22—28.

Xia Heng-heng, Li Zhi-nong, Xiao Yao-xian. Stiffness influencing factors analysis of the fully opened slant crack rotor[J]. Journal of Nanchang Hangkong University (Natural Science) , 2014，28(3): 22—28.

[12] 孫訓(xùn)方, 方孝淑, 陸耀洪. 材料力學(xué)[M].第3版.北京: 高等教育出版社, 2012: 152—161.

Sun Xun-fang, Fang Xiao-shu, Lu Yao-hong. Mechanics of Materials[M]. 3rd ed. Beijing: Higher Education Press, 2012: 152—161.