顧 偉 張 博 丁 虎 陳立群
(上海大學(xué)力學(xué)與工程科學(xué)學(xué)院,上海 200444)
(長安大學(xué)理學(xué)院,西安 710064)
(上海大學(xué)上海市應(yīng)用數(shù)學(xué)和力學(xué)研究所,上海 200444)
(上海大學(xué)上海市能源工程力學(xué)重點實驗室,上海 200444)
在實際工程結(jié)構(gòu)中,葉片是最重要的傳動部件之一,如在渦輪機(jī)葉片、船用螺旋槳、直升機(jī)葉片以及風(fēng)力發(fā)電裝置中都有廣泛應(yīng)用.近年來,隨著人們對機(jī)械性能要求的提高,葉片的工作環(huán)境更加嚴(yán)峻,常需在復(fù)雜工況下(如高溫、摩擦等)切換轉(zhuǎn)速,極易使葉片產(chǎn)生變形,導(dǎo)致發(fā)生破壞,因此研究預(yù)變形變轉(zhuǎn)速葉片的動力學(xué)特性對于合理設(shè)計葉片結(jié)構(gòu)具有重要意義.
由于機(jī)械工程上的大量應(yīng)用,旋轉(zhuǎn)葉片的動力學(xué)特性長期以來都是國內(nèi)外學(xué)者關(guān)注的重點,許多學(xué)者采用理論分析,數(shù)值模擬以及實驗方法針對不同葉片模型,在線性或非線性框架下,研究了常轉(zhuǎn)速下旋轉(zhuǎn)葉片的動力學(xué)響應(yīng)規(guī)律.前期,一些學(xué)者[1-2]研究了在常轉(zhuǎn)速下,考慮旋轉(zhuǎn)葉片模型的建立問題,得出其動力學(xué)方程,通過不同方法求解其自由振動等問題.Yoo 等[3-5]通過混合坐標(biāo)系建模研究了常轉(zhuǎn)速下預(yù)扭懸臂梁受不同外激勵下的振動特性.蔡國平和洪嘉振[6]通過變結(jié)構(gòu)控制法對勻速轉(zhuǎn)動懸臂梁進(jìn)行主動控制,并且比較零次近似模型和一次近似模型的巨大差異.Yang 等[7]利用廣義哈密頓原理推導(dǎo)了常轉(zhuǎn)速下Euler-Bernoulli 梁橫向振動的非線性微分方程.姜靜和韓廣才[8]研究了常轉(zhuǎn)速葉片在彎扭耦合下由于質(zhì)量偏心的動力學(xué)分析.Inoue 等[9]研究了實際問題中定轉(zhuǎn)速旋轉(zhuǎn)的風(fēng)力發(fā)電機(jī)受重力和風(fēng)力影響的振動特性.Ramakrishnan 和Feeny[10]進(jìn)一步在模型建立過程中加入非線性項來研究其動力學(xué)特性,隨后Ramakrishnan 和Feeny[11]建立了強(qiáng)迫馬蒂厄方程來考慮風(fēng)力葉片的共振,李國強(qiáng)等[12]通過風(fēng)洞試驗研究翼型橫擺震蕩.趙國威和吳志剛[13]在文獻(xiàn)[6]的基礎(chǔ)上進(jìn)一步用一次近似模型考慮將軸向和橫向變形產(chǎn)生的耦合項代入應(yīng)變能中,考慮其對動力學(xué)特性的影響.楊鄂川等[14]研究了常轉(zhuǎn)速下,葉片受開口裂紋作用下的振動特性,觀察梁的一階、二階固有頻率隨不同裂紋深度和位置的變化.Oh 和Yoo[15]考慮了拉伸、彎曲和扭轉(zhuǎn)的耦合作用下,常轉(zhuǎn)速預(yù)扭轉(zhuǎn)葉片的振動分析.Tian 等[16]提出一種改進(jìn)的變分方法,推導(dǎo)了旋轉(zhuǎn)梁的自由和瞬態(tài)振動特性.李炳強(qiáng)等[17]研究了在常轉(zhuǎn)速下,轉(zhuǎn)子--葉片耦合系統(tǒng)在主共振條件下的振動響應(yīng)和動態(tài)穩(wěn)定性.吳吉等[18]基于絕對節(jié)點坐標(biāo)法研究懸臂曲梁的純彎曲問題和帶有剛?cè)狁詈系男D(zhuǎn)柔性非線性動力學(xué)響應(yīng)問題.
由于轉(zhuǎn)子圓盤葉片系統(tǒng)經(jīng)常會受外部擾動,考慮轉(zhuǎn)速隨時間變化更加普遍也更貼近實際.例如當(dāng)發(fā)動機(jī)啟動、變速、停機(jī)等過程中,轉(zhuǎn)子的輸入與輸出功率(扭矩)失衡.轉(zhuǎn)子在轉(zhuǎn)動過程中將伴隨扭振,從而產(chǎn)生了速度脈動,使得葉片轉(zhuǎn)速不再是一個常數(shù).近年來,有越來越多的學(xué)者將目光轉(zhuǎn)移到變轉(zhuǎn)速葉片的動力學(xué)行為的研究中.Kammer 和Schlack[19]首先將角速度表示為穩(wěn)態(tài)值和微小擾動的疊加,探究變轉(zhuǎn)速對梁振動的影響.Young[20]研究了變轉(zhuǎn)速下,預(yù)扭轉(zhuǎn)變截面梁的動態(tài)響應(yīng),利用多尺度法處理方程,并研究阻尼系數(shù)對不穩(wěn)定區(qū)域的影響.Yang 和Tsao[21]研究了非定常轉(zhuǎn)速下,預(yù)扭轉(zhuǎn)葉片振動與穩(wěn)定性分析.張偉等[22]針對1/2 亞諧共振和1:3 內(nèi)共振,通過數(shù)值模擬得到二維、三維相圖,分析了風(fēng)速變化對變轉(zhuǎn)速葉片振動的影響.Yao 等[23]考慮了在高溫超聲波氣流變轉(zhuǎn)速葉片的動力學(xué)響應(yīng).隨后,Yao 等[24]研究了葉片在不同轉(zhuǎn)速下的非線性振動和穩(wěn)態(tài)響應(yīng).Georgiades[25]研究了在非定常轉(zhuǎn)速彈性連續(xù)體中傳動軸動力學(xué)行為.張偉等[26]建立了變截面葉片模型,分析在變轉(zhuǎn)速下葉片受氣動力和離心力影響的非線性動力學(xué)行為.Hu 等[27]提出一種新型的三維預(yù)扭梁單元公式來捕捉變轉(zhuǎn)速梁高度耦合的振動特性.
近期,有學(xué)者報道了葉片結(jié)構(gòu)的預(yù)變形將嚴(yán)重影響葉片的動力學(xué)響應(yīng).Zhang 和Li[28]研究了在常轉(zhuǎn)速下,預(yù)扭轉(zhuǎn)葉片在2:1 內(nèi)共振下的受迫振動.Kang 等[29]基于載荷增量法,提出一種新的預(yù)變形方法來補(bǔ)償葉片變形,通過修正剛度矩陣考慮非線性葉片剛度的影響.Zhang 等[30]研究在2:1 內(nèi)共振下旋轉(zhuǎn)預(yù)扭轉(zhuǎn)梁的主共振情況,詳細(xì)研究了不同參數(shù)變化對幅頻響應(yīng)曲線的影響.
通過文獻(xiàn)調(diào)研發(fā)現(xiàn),較多文獻(xiàn)考慮變轉(zhuǎn)速情況時往往不考慮預(yù)變形作用,而文獻(xiàn)[28]主要考慮常轉(zhuǎn)速下預(yù)變形葉片受外激勵作用.同時考慮變轉(zhuǎn)速和預(yù)變形兩種因素共同作用下的葉片動力學(xué)行為相關(guān)的研究還尚未見報道.本文主要研究變轉(zhuǎn)速預(yù)變形葉片在2:1 內(nèi)共振下受參數(shù)激勵的響應(yīng)并考慮立方非線性項對系統(tǒng)的影響.
如圖1 所示為一長度為L的變轉(zhuǎn)速預(yù)變形細(xì)長葉片.其一端固結(jié)在轉(zhuǎn)軸半徑為r的圓盤上,另一端自由,考慮葉片轉(zhuǎn)速由一定常轉(zhuǎn)速?1和一簡諧變化的微小擾動?2sin(!t)疊加而成,即?(t)=?1+?2sin(!t).葉片在安裝時與圓盤有一定的角度,稱為安裝角Ψ.葉片的總預(yù)扭角為Θ.由于受溫度梯度影響,葉片產(chǎn)生預(yù)變形u20和u30,分別表示弦向和翼向,忽略軸向預(yù)變形(軸向預(yù)變形遠(yuǎn)小于其他兩個方向),利用旋轉(zhuǎn)坐標(biāo)可得到在熱梯度作用下預(yù)扭葉片沿截面兩主軸的曲率函數(shù),見文獻(xiàn)[28].為了描述葉片振動過程中的構(gòu)形,本文建立了兩套坐標(biāo)系,一套慣性坐標(biāo)系XYZ固定在圓盤中心,一套旋轉(zhuǎn)坐標(biāo)系xyz連接在旋轉(zhuǎn)葉片根部上,x沿著未變形葉片的軸線方向(左右),y為弦向(前后),z為翼向(上下).
圖1 變轉(zhuǎn)速預(yù)變形葉片示意圖Fig.1 Diagram of pre-deformed blade with varying speed
為了推導(dǎo)出葉片的解析模型,本文作出如文獻(xiàn)[5,28]中假設(shè),根據(jù)Von-Karman 應(yīng)變--位移公式
上式中,u1,u2,u3表示為梁軸線上任一點處位移,通過上述假設(shè),變形能可寫成如下形式
由于離心力作用產(chǎn)生的軸向收縮勢能可表達(dá)為
可得總勢能為
旋轉(zhuǎn)葉片的動能可表達(dá)如下
其中,頂標(biāo) 表示對時間求微分,在此文中對于外力做功,只考慮黏性阻尼力的作用,阻尼系數(shù)為cd.
利用假設(shè)模態(tài)法,將位移函數(shù)分解成關(guān)于時間和空間的函數(shù),空間函數(shù)表達(dá)如下
上式中參數(shù)βj,可見文獻(xiàn)[31].
將式(4)和式(5)代入Lagrange 動力學(xué)方程中,為了使方程更具普遍性,引入無量綱參數(shù),如下
對方程進(jìn)行處理,為了書寫的簡潔性,將式(7)中無量綱參數(shù)中上劃線符號忽略不寫,可得下列無量綱化動力學(xué)方程
上式方程中下劃線代表與文獻(xiàn)[28]一文中受迫振動不同項,方程(8)控制著軸向運(yùn)動,方程(9)、方程(10)控制著弦向和翼向,從上式中可以發(fā)現(xiàn)加上預(yù)變形,方程中多了平方項和立方項,令預(yù)變形項為零、轉(zhuǎn)速變?yōu)槌^D(zhuǎn)速并且不考慮外力項,方程就退化成文獻(xiàn)[5]的模型.
正如文獻(xiàn)[5]中所說,對于細(xì)長歐拉梁,其軸向第一階固有頻率遠(yuǎn)高于彎曲振動固有頻率,葉片系統(tǒng)前幾階模態(tài)主要表現(xiàn)為彎曲變形,因此忽略軸向收縮變形與彎曲變形耦合作用,上式只留下方程(9)和(10),同時將方程中雙下劃線部分(跟q1j相關(guān)項)略去,并通過模態(tài)矩陣[30]對線性部分進(jìn)行解耦,可得下式
上式方程中各系數(shù)可見附錄,ηijk,ξijk參數(shù)可見文獻(xiàn)[32]對方程(11)進(jìn)行重刻度
并將非線性方程穩(wěn)態(tài)解展開成下式
其中,T0=t為快時間尺度,T1=為慢時間尺度,將方程(13)代入方程(11)中,分別取階可得
從方程(14)中可求得
其中,Ai是關(guān)于T1的復(fù)函數(shù),cc表示前面各項的復(fù)共軛項.Ai和Bi表達(dá)式如下
其中,ai和ζi分別為幅值和相角,均為慢時間尺度T1的函數(shù),將式(16)代入式(15)中,可得
文獻(xiàn)[28]的圖3 和圖4 詳細(xì)研究了旋轉(zhuǎn)預(yù)變形葉片隨轉(zhuǎn)速、預(yù)變形程度變化過程中,系統(tǒng)前兩階固有頻率的變化規(guī)律.發(fā)現(xiàn)在特定的組合參數(shù)下,系統(tǒng)具有發(fā)生2:1 內(nèi)共振的可能性.為研究2:1 內(nèi)共振下,非線性振動系統(tǒng)發(fā)生一階次諧波共振的振動特性,引入兩個解諧參數(shù)σ1和σ2,用σ1表示接近的程度,用σ2表示接近的程度
從方程(18)中消除久期項,分離實部和虛部可得如下的演化方程
其中,上標(biāo)I 表示虛部,R 表示實部
當(dāng)a1,a2,與時間T1無關(guān)時,可以求解得到系統(tǒng)的穩(wěn)態(tài)響應(yīng).因此方程(20)可令等式左邊為零,并求解對應(yīng)的非線性代數(shù)方程組即得到系統(tǒng)的穩(wěn)態(tài)響應(yīng),解分為兩種情況:(i)a1=0,a2≠0 稱為單模態(tài)解;(ii)a1≠0,a2≠0 稱為雙模態(tài)解,系統(tǒng)的穩(wěn)定性可以通過Lyapunov 準(zhǔn)則判斷,具體可見文獻(xiàn)[32-33].
圖2~圖4 分別考慮了溫度梯度、阻尼系數(shù)以及轉(zhuǎn)速擾動幅值對模態(tài)頻率響應(yīng)的影響,圖中實線表示穩(wěn)定區(qū)域,虛線表示不穩(wěn)定區(qū)域,由上述圖可以得知在2:1 內(nèi)共振中,模態(tài)頻率響應(yīng)曲線中會同時出現(xiàn)軟特性和硬特性.由圖2 可知,隨著溫度梯度的增大,一階、二階模態(tài)頻率響應(yīng)曲線有向外張開的趨勢且峰值增大,其不穩(wěn)定區(qū)域變大;在圖3 中,隨著阻尼系數(shù)的增大,曲線仍有向外張開的趨勢,但峰值降低,表明當(dāng)有阻尼約束時,對系統(tǒng)振動有較大的抑制作用,由圖3 可見,阻尼系數(shù)由0.1 變?yōu)?.3,響應(yīng)峰值降低近3 倍,同時不穩(wěn)定區(qū)域變??;在圖4 中,隨著擾動轉(zhuǎn)速的增大,曲線呈現(xiàn)出向內(nèi)閉合的趨勢且峰值增大,并且不同于溫度梯度影響,在不同擾動轉(zhuǎn)速在同一頻率下,其幅值不會相交.
圖2 溫度梯度變化對前兩階模態(tài)頻率響應(yīng)曲線的影響Fig.2 The effect of temperature gradient on the responses of the first two modal frequency
圖3 阻尼變化對前兩階模態(tài)頻率響應(yīng)曲線的影響Fig.3 The effect of damping on the responses of the first two modal frequency
圖4 轉(zhuǎn)速擾動幅值變化對前兩階模態(tài)頻率響應(yīng)曲線的影響Fig.4 The effect of amplitude of revolution disturbance on the responses of the first two modal frequency
為了驗證本文建模的正確性,忽略本文建立的動力學(xué)模型中預(yù)變形和轉(zhuǎn)速的影響,將本文模型得到的前四階固有頻率與其他文獻(xiàn)進(jìn)行了對比,如表1 表示.
表1 中與文獻(xiàn)[5,34]進(jìn)行了比較,發(fā)現(xiàn)前四階固有頻率最大誤差不超過0.1%.
為了驗證多尺度法得到的系統(tǒng)穩(wěn)態(tài)解的正確性,本節(jié)采用Runge-Kutta 法對系統(tǒng)動力學(xué)方程(11)進(jìn)行數(shù)值積分.在數(shù)值積分中,仿真時間設(shè)定足夠長,以確保系統(tǒng)處于穩(wěn)態(tài)響應(yīng),同時,數(shù)值積分初值在單模態(tài)穩(wěn)定解附近選取.
表1 前四階固有頻率的比較(γ=0,κ=0:25,δ=2)Table 1 Comparison of the lowest four natural frequencies(γ=0,κ=0:25,δ=2)
通過正向、反向掃頻研究立方項在2:1 內(nèi)共振中對方程的影響程度,結(jié)果如下.
圖5 中圓圈代表對原系統(tǒng)數(shù)值積分中不考慮立方非線性項影響,加號代表考慮立方非線性項,實線代表解析解,紅色虛線代表不穩(wěn)定區(qū)域.可以發(fā)現(xiàn)在2:1 內(nèi)共振下,立方非線性對系統(tǒng)響應(yīng)影響很小,僅在完全內(nèi)共振處有細(xì)微偏差.同時,在正向和反向掃頻中,發(fā)現(xiàn)有兩處跳躍點,在頻響曲線的多解頻帶中,系統(tǒng)有兩個漸進(jìn)穩(wěn)定解和一個不穩(wěn)定解,由于在實驗中只能實現(xiàn)漸進(jìn)穩(wěn)定運(yùn)動,因此在正反兩個方向掃頻激勵中會出現(xiàn)圖示路徑的跳躍現(xiàn)象,并且跳躍點附近雙模態(tài)解吸引域較小,穩(wěn)態(tài)解吸引到單模態(tài)解中,振幅發(fā)生突變.解析解和數(shù)值解吻合較好,從側(cè)面證明多尺度處理的正確性.
此外,利用控制變量的思想,探究了擾動速度γ2和解諧參數(shù)σ2的改變對系統(tǒng)動力響應(yīng)的影響.
圖6 中,當(dāng)固定σ2=0:8 時(圖6(a),圖6(b),圖6(c)),隨著γ2的增大,系統(tǒng)經(jīng)歷了單周期擬周期單周期的變化歷程.類似的,當(dāng)固定γ2=0:000 9 時(圖6(d),圖6(b),圖6(e),圖6(f)),系統(tǒng)經(jīng)歷了單周期擬周期單周期多周期的變化歷程.可以看出系統(tǒng)動力學(xué)行為與σ2,γ2緊密相關(guān),隨參數(shù)的細(xì)微變化,系統(tǒng)動力學(xué)響應(yīng)會出現(xiàn)分岔等現(xiàn)象.
圖5 正反向掃頻激勵數(shù)值仿真與解析解的對比Fig.5 Comparison of numerical simulation and analytical solution with forward sweep and backward sweep
圖6 不同γ2 和σ2 時系統(tǒng)的相圖Fig.6 The phase diagrams of system for different γ2 and σ2
圖6 不同γ2 和σ2 時系統(tǒng)的相圖(續(xù))Fig.6 The phase diagrams of system for different γ2 and σ2 (continued)
圖6 不同γ2 和σ2 時系統(tǒng)的相圖(續(xù))Fig.6 The phase diagrams of system for different γ2 and σ2 (continued)
本文主要研究了參數(shù)激勵下,預(yù)變形變轉(zhuǎn)速葉片的非線性動力學(xué)特性,研究不同參數(shù)變化對系統(tǒng)模態(tài)頻率響應(yīng)的影響,并考慮立方項在2:1 內(nèi)共振時對方程的影響,最后利用數(shù)值積分與解析解進(jìn)行驗證.結(jié)果發(fā)現(xiàn):隨著溫度梯度增大時,模態(tài)頻率響應(yīng)曲線有向外張開的趨勢,同時幅值增大,不穩(wěn)定區(qū)域增大;當(dāng)轉(zhuǎn)速擾動幅值增大時,幅值同樣增大,不穩(wěn)定區(qū)域增大,不同于溫度梯度,其模態(tài)頻率響應(yīng)曲線向內(nèi)閉合,并且不同擾動轉(zhuǎn)速在同一固有頻率下,其幅值不會相交;阻尼系數(shù)改變對系統(tǒng)的幅值影響最大,當(dāng)阻尼系數(shù)由0.1 變?yōu)?.3,其幅值降低近3 倍,同時,不穩(wěn)定區(qū)域減?。涣⒎巾椩?:1 內(nèi)共振中對方程的影響可近似忽略,利用原系統(tǒng)數(shù)值積分對解析解進(jìn)行驗證,發(fā)現(xiàn)兩者吻合較好.
附錄A
式(11)中各矩陣及其表達(dá)式如下