【正文】 遼寧工程技 術(shù)大學(xué)畢業(yè)設(shè)計(jì)（論文） 1 1 緒論 問(wèn)題的提出 輪盤是發(fā)動(dòng)機(jī)的重要零件之一 ,亦是主要受力零件 .輪盤結(jié)構(gòu)復(fù)雜 ,轉(zhuǎn)速高 ,其 工作條件 十分惡劣 ,設(shè)計(jì)要求高、難度大 .輪盤不僅在高轉(zhuǎn)速下承受較大的外負(fù)荷 ,由于結(jié)構(gòu)要求 ,大部分輪盤都開有中心 孔或其它通氣孔 ,還開有榫槽以便于與葉片聯(lián)接 ,造成輪盤局部應(yīng)力 集中 .渦輪盤的徑向、軸向還有較 大的溫差 ,承受很大的熱應(yīng)力 ,承受工作轉(zhuǎn)速和溫度循環(huán)變化的疲勞應(yīng)力 .渦輪發(fā)動(dòng)機(jī)常見(jiàn)的輪盤故 障之一就是輪盤由于共振而產(chǎn)生的開裂破壞 .如某型發(fā)動(dòng)機(jī)壓氣機(jī)九級(jí)盤曾發(fā)生輪緣爆破破壞 ,造成重 大故障 ,其 原因 是 輪盤產(chǎn)生了 一階傘 型振動(dòng)的共振 ,導(dǎo)致輪緣產(chǎn)生疲勞損傷而爆破 .由于日益追求發(fā) 動(dòng)機(jī)重量輕和壽命長(zhǎng)， 輪盤結(jié)構(gòu)日趨輕型化 ,輪盤做得很薄 。 因此在高負(fù)荷、長(zhǎng)壽命的工作情況 下 ,輪盤因振動(dòng)疲勞而斷裂的故障有所增多 。 旋轉(zhuǎn)輪盤在高速旋轉(zhuǎn)狀態(tài)下會(huì)表現(xiàn)出與低速或非旋轉(zhuǎn)狀態(tài)下迥異的力學(xué)性能。許多轉(zhuǎn)子機(jī)械受剛度特性和非保守效應(yīng)的阻尼特性的影響，其運(yùn)轉(zhuǎn)工況下的動(dòng)態(tài)性能與靜止時(shí)相比有很大的差異。在進(jìn)行高速旋轉(zhuǎn)機(jī)械的轉(zhuǎn)子系統(tǒng)動(dòng)力設(shè)計(jì)時(shí)，需要對(duì)轉(zhuǎn)動(dòng)部件進(jìn)行模態(tài)分析，求解出其固有頻率和相應(yīng)的模態(tài)振型。通過(guò)合理的設(shè)計(jì)使其工作轉(zhuǎn)速盡量遠(yuǎn)離轉(zhuǎn)子系統(tǒng)的固有頻率以避免由 于共振而產(chǎn)生開裂破壞。而對(duì)于高速部件，工作時(shí)由于受到離心力的影響，其固有頻率跟靜止時(shí)相比會(huì)有一定的變化。為此，在進(jìn)行模態(tài)分析時(shí)需要考慮離心力引起的預(yù)應(yīng)力的影響。如果忽略了預(yù)應(yīng)力的影響就會(huì)使原有的系統(tǒng)遭到破壞，降低系統(tǒng)的穩(wěn)定性和可靠性。因此掌握高速旋轉(zhuǎn)輪盤運(yùn)轉(zhuǎn)工況各種轉(zhuǎn)速下的動(dòng)力特性，對(duì)于了解不平衡振動(dòng)響應(yīng)以及可能出現(xiàn)的不穩(wěn)定性等均有實(shí)際意義。 國(guó)內(nèi)外研究現(xiàn)狀 輪盤是旋轉(zhuǎn)機(jī)械中的常見(jiàn)部件，如水泵葉輪、風(fēng)機(jī)葉輪輪盤、汽輪機(jī)轉(zhuǎn)子上的輪盤等，它的強(qiáng)度、振動(dòng)是設(shè)計(jì)、加工中重點(diǎn)考慮的問(wèn)題之一。由于 輪盤結(jié)構(gòu)復(fù)雜 , 因此，在利用有限元分析計(jì)算時(shí)， 建模是很重要的一個(gè)環(huán)節(jié)。 文獻(xiàn) [1]采用在 UG 上進(jìn)行幾何建模，運(yùn)用NASTRAN 建立有限元模型，對(duì)輪盤進(jìn)行靜強(qiáng)度分析和模態(tài)分析。 for ( MSC公司推出的最新的運(yùn)行于 windows環(huán)境的有限元分析工具。它可 以在 windows95 或 windowsNT 下運(yùn)行。 可以進(jìn)行靜力、模態(tài)、屈曲、熱傳導(dǎo)、 動(dòng)力學(xué)、非線性和設(shè)計(jì)優(yōu)化等各種類型分析 。 擁有綜合全面和先進(jìn)的網(wǎng) 格劃分技術(shù)，根據(jù)不同的幾何模型提 供了多種不同的生成和定義有限元模型工具，包括 : 丁文文：預(yù)應(yīng)力對(duì)高速旋轉(zhuǎn)輪盤模態(tài)的影響分析 2 網(wǎng)格劃分、邊界和載荷定義。 預(yù)應(yīng)力技術(shù)能在不增加零件尺寸及材料性能的前提下，提高其承載能力，使用性能及工作壽命。它是使零件在承受工作載荷之前被施加一定的反向預(yù)應(yīng)力，當(dāng)其承受工作載荷時(shí)，工作應(yīng)力與預(yù)應(yīng)力疊減，從而使零件的實(shí)際應(yīng)力得以減小，承載能力提高，變形量減小，壽命延長(zhǎng)。對(duì)于以離心力為其主要工作荷載的高速旋轉(zhuǎn)零件。如離心壓縮葉機(jī)，高速泵葉輪，高速旋轉(zhuǎn)刀具等，也可以通過(guò)預(yù)應(yīng)力處理來(lái)提高其承載能力，延長(zhǎng)壽命。 文獻(xiàn) [2]從基本理論出發(fā)，對(duì)該類零件在一定轉(zhuǎn)速下 的動(dòng)態(tài)固有頻率作三維 有限元分析。經(jīng)過(guò)理論分析和實(shí)例證明，由于離心力的影響，旋轉(zhuǎn)輪盤 低階模態(tài)的動(dòng)態(tài)固有頻率增加較大。另外，輪盤的厚度越薄，轉(zhuǎn)速越高，其影響 也越 大。 文獻(xiàn) [3]采用激光全息干涉法對(duì)某型發(fā)動(dòng)機(jī)一級(jí)壓氣機(jī)盤進(jìn)行振動(dòng)試驗(yàn)研究 ,在 0～ 15 000Hz 頻率范圍內(nèi)獲得了 7 個(gè)振動(dòng)模 態(tài)的靜頻及振型 ,并根據(jù)試驗(yàn)結(jié)果研究該輪盤的行波振動(dòng) ,分析輪盤的共振特性及臨界轉(zhuǎn)速。 文獻(xiàn) [4]用有限元等數(shù)值分析的方法，探討了葉片動(dòng)力剛化對(duì)其振動(dòng)特性的影響。利用國(guó)際流行的有限元分析軟件 ANSYS，對(duì)某型壓 氣機(jī)的盤鼓在工作過(guò)程中的應(yīng) 力狀態(tài)以及振動(dòng)模態(tài)進(jìn)行了仿真計(jì)算。 文獻(xiàn) [5]針對(duì) 600MW 汽輪發(fā)電機(jī)組低壓轉(zhuǎn)子 軸承系統(tǒng) ,建立了整軸模型和低壓轉(zhuǎn)子模型 ,采用 有限元分析軟件 ANSYS 進(jìn)行模態(tài)分析 時(shí) 提出使用 Lagrange 乘子法來(lái)處理約束方程， 計(jì)算汽輪機(jī)轉(zhuǎn)子的固有頻率和臨界轉(zhuǎn)速 ,從而避免工作轉(zhuǎn)速達(dá)到 臨界轉(zhuǎn)速產(chǎn)生共振現(xiàn)象，最后通過(guò)臨界轉(zhuǎn)速和振型圖分析了低壓轉(zhuǎn)子的特性。 文獻(xiàn) [6]簡(jiǎn)述了通過(guò)預(yù)應(yīng)力提高旋轉(zhuǎn)零件承載能力與壽命的基本原理，給出了最低及最高兩個(gè)極限試驗(yàn)轉(zhuǎn)速，并對(duì)試驗(yàn)設(shè)備的技術(shù)性能提出要求。 隨著制造水平的提高以及對(duì)新型軸承結(jié)構(gòu)的理論探 索 ,外存儲(chǔ)器正朝著超大容量、高數(shù)據(jù)傳輸率和小尺寸方向發(fā)展 .這意味著存儲(chǔ)器主軸轉(zhuǎn)速將越來(lái)越快 ,盤片的厚度越來(lái)越薄 ,因此 ,盤 /軸系統(tǒng)高速旋轉(zhuǎn)時(shí) ,盤片的彈性變形主要源于旋轉(zhuǎn)時(shí)的軸對(duì)稱膜應(yīng)力 離心力 ,且盤片的變形將對(duì)整個(gè)系統(tǒng)的動(dòng)態(tài)特性及性能影響最大 .國(guó)內(nèi)曾用理論計(jì)算和實(shí)驗(yàn)方法研究過(guò)磁盤在不旋轉(zhuǎn)時(shí)的自由振動(dòng)和動(dòng)態(tài)特性 ,并給出了部分振動(dòng)模態(tài)圖 . 文獻(xiàn) [7]研究給出了盤 /軸系統(tǒng)高速旋轉(zhuǎn)時(shí)盤片振動(dòng)模態(tài)的計(jì)算方法 ,并開發(fā)了振動(dòng)模態(tài)的仿真程序 ,給出了部分仿真結(jié)果 . 文獻(xiàn) [8]中以磁懸浮軸承轉(zhuǎn)子系統(tǒng)的組成及工作原理為基礎(chǔ) ,在 ANSYS9. 0 中建立了磁懸浮轉(zhuǎn)子的三維有限元模型。采用 Subspace 法計(jì)算了前 4 階固有頻率和振型 ,并與試驗(yàn)?zāi)＿|寧工程技 術(shù)大學(xué)畢業(yè)設(shè)計(jì)（論文） 3 態(tài)分析結(jié)果進(jìn)行了對(duì)比。結(jié)果顯示 ,基于有限元法得到的數(shù)據(jù)與試驗(yàn)?zāi)B(tài)分析基本一致 ,所建立的有限元模型為柔性轉(zhuǎn)子系統(tǒng)設(shè)計(jì)提供了一定的理論依據(jù)。 文獻(xiàn) [9]建立了葉片 輪盤耦合系統(tǒng)振動(dòng)固有特性計(jì)算分析模型 ,利用有限元通用程序NASTARN 的循環(huán) 對(duì)稱分析功能和幾何非線性分析功能 ,對(duì)某艦用燃?xì)廨啓C(jī)葉片 輪盤結(jié)構(gòu)進(jìn)行計(jì)算分析。計(jì)算中考慮了旋轉(zhuǎn) 離心力及材料參數(shù)隨溫度變化對(duì)固有特性的影響 ,得到的諧波共振頻率與實(shí)測(cè)結(jié)果 基本吻合 。 本文的主要工作 旋轉(zhuǎn)輪盤是許多機(jī)械裝置最基本的也是最重要的零件之一。輪盤在高速旋轉(zhuǎn)狀態(tài)下由于離心力的影響，其固有頻率與靜止時(shí)相比也會(huì)發(fā)生一定的變化，因此在對(duì)輪盤進(jìn)行模態(tài)分析不可忽略預(yù)應(yīng)力的影響。本文主要做了以下工作： 1） 本文首先介紹了有限元法以及 ANSYS 軟件，建立了高速旋轉(zhuǎn)輪盤的有限元模 型。 2） 介紹靜力學(xué)分析基本理論，分析了輪盤在相對(duì)穩(wěn)定的離心力載荷的作用下的靜 應(yīng)力及位移的問(wèn)題。 3） 在此靜應(yīng)力分析的基礎(chǔ)上，基于旋轉(zhuǎn)體動(dòng)力學(xué)方程，略去初位移剛度矩陣，建立了輪 盤的有限元形式動(dòng)力學(xué)方程。并基于此方程作了模態(tài)分析，詳細(xì)分析輪盤在旋轉(zhuǎn)時(shí)的振動(dòng)特性。利用 ANSYS 軟件，直接求出輪盤在靜 止?fàn)?態(tài)和 不同 旋轉(zhuǎn) 速度 狀 態(tài)下整體結(jié)構(gòu)的固有頻率、振型，并給出振型的彩色云圖和動(dòng)畫顯示；通過(guò)對(duì)固有頻率、振型的分析得出預(yù)應(yīng)力對(duì)高速旋轉(zhuǎn)輪盤模態(tài)的影響規(guī)律。 丁文文：預(yù)應(yīng)力對(duì)高速旋轉(zhuǎn)輪盤模態(tài)的影響分析 4 2 高速旋轉(zhuǎn)輪盤模態(tài)分析基本理論與方法 在實(shí)際工作中， 輪盤經(jīng)常要受到隨時(shí)間變化的動(dòng)載荷的作用。當(dāng)所受的動(dòng)載荷較大 時(shí) ,或者雖然不大 ,當(dāng) 作用力的頻率與輪盤的某一階固有頻率接近時(shí) , 輪盤 將產(chǎn)生強(qiáng)烈的共振 ,從而引起很高的動(dòng)應(yīng)力 , 造成輪盤 強(qiáng)度破壞或變形。因此 , 有必要對(duì)輪盤 作動(dòng)力學(xué)模態(tài) 分析 , 也就是要分析輪盤 的固有頻率、主振型 , 以及在動(dòng)載荷作用下的動(dòng)撓度、動(dòng)應(yīng)力等等。 模態(tài)分析概述 模態(tài)分析的定義 一般地，以振動(dòng)理論為基礎(chǔ)、以模態(tài)參數(shù)為目的的分析方法，稱為模態(tài)分析。更確切地說(shuō)，模態(tài)分析是研究系統(tǒng)物理參數(shù)模型、模態(tài)參數(shù)模型和非參數(shù)模型的關(guān)系，并通過(guò)一定手段確定這些系統(tǒng)模型理論及其應(yīng)用的一門學(xué)科 [10]。通過(guò)模態(tài)分析得到的結(jié)果，可以幫助結(jié)構(gòu)設(shè)計(jì)者設(shè)計(jì)合理結(jié)構(gòu)以避開共振頻率。 模態(tài)分析的實(shí)質(zhì)是一種坐 標(biāo)變換。其目的在于把原在物理坐標(biāo)系統(tǒng)中描述的響應(yīng)向量，放到所謂“模態(tài)坐標(biāo)系統(tǒng)”中描述。這一坐標(biāo)系統(tǒng)的每一個(gè)基向量恰是振動(dòng)系統(tǒng)的一個(gè)特征向量。運(yùn)用這一坐標(biāo)系統(tǒng)的好處是：利用各特征向量之間的正交特性，可使描述響應(yīng)向量的各個(gè)坐標(biāo)均可單獨(dú)求解。換句話說(shuō)，在這一坐標(biāo)系統(tǒng)下，振動(dòng)方程是一組互無(wú)耦合的方程，每一個(gè)坐標(biāo)均可單獨(dú)求解。 根據(jù)研究模態(tài)分析的手段和方法的不同，模態(tài)分析分為理論模態(tài)分析和實(shí)驗(yàn)?zāi)B(tài)分析兩種。理論模態(tài)分析是指以線性振動(dòng)理論為基礎(chǔ)，研究激勵(lì)、系統(tǒng)和響應(yīng)三者的關(guān)系；實(shí)驗(yàn)?zāi)B(tài)分析又稱模態(tài)分析的實(shí)驗(yàn)過(guò)程，是解析 模態(tài)分析的逆過(guò)程。 計(jì)算模態(tài)分析實(shí)際上是一種理論建模過(guò)程，主要是運(yùn)用有限元法對(duì)振動(dòng)結(jié)構(gòu)進(jìn)行離散，建立系統(tǒng)特征值問(wèn)題的數(shù)學(xué)模型，用各種近似方法求解系統(tǒng)特征值和特征向量。 模態(tài)分析在工程中的應(yīng)用 1） 模態(tài)分析在結(jié)構(gòu)性能評(píng)價(jià)中的直接應(yīng)用 根據(jù)模態(tài)分析的結(jié)果，即模態(tài)頻率、模態(tài)振型、模態(tài)阻尼等模態(tài)參數(shù)，對(duì)被測(cè)結(jié)構(gòu)進(jìn)行直接的動(dòng)態(tài)性能評(píng)估。對(duì)一般結(jié)構(gòu)，要求各階模態(tài)頻率遠(yuǎn)離工作頻率，或工作頻率不落在某階模態(tài)的半功率帶寬內(nèi)；對(duì)結(jié)構(gòu)貢獻(xiàn)較大的振型，應(yīng)使其不影響結(jié)構(gòu)正常工作為最佳。這是模態(tài)分析的直接 應(yīng)用，已成為工程界的基本方法。 2） 模態(tài)分析在結(jié)構(gòu)動(dòng)態(tài)設(shè)計(jì)中的應(yīng)用 遼寧工程技 術(shù)大學(xué)畢業(yè)設(shè)計(jì)（論文） 5 以模態(tài)分析為基礎(chǔ)的結(jié)構(gòu)動(dòng)態(tài)設(shè)計(jì)，是近年來(lái)振動(dòng)工程界開展的最廣泛的研究領(lǐng)域之一。傳統(tǒng)的結(jié)構(gòu)設(shè)計(jì)是以經(jīng)驗(yàn)和反復(fù)實(shí)測(cè)為主要手段。因?yàn)楸M管依據(jù)模態(tài)分析結(jié)果和響應(yīng)實(shí)驗(yàn)容易判斷出初步結(jié)構(gòu)的性能缺陷，但在結(jié)構(gòu)修改問(wèn)題上卻往往茫然無(wú)所知，設(shè)計(jì)工程師只能依據(jù)經(jīng)驗(yàn)和現(xiàn)有條件進(jìn)行反復(fù)修改和實(shí)測(cè)，有時(shí)甚至將原設(shè)計(jì)完全推翻重新設(shè)計(jì)。這樣就大大減緩了設(shè)計(jì)速度，設(shè)計(jì)質(zhì)量也難以達(dá)到最優(yōu)。 3） 模態(tài)分析在故障診斷和狀態(tài)檢測(cè)中的應(yīng)用 利用模態(tài)分析得到的模態(tài)參數(shù)等結(jié)果進(jìn)行故障判 別，日益成為一種有效而實(shí)用的故障診斷和安全檢驗(yàn)方法。如根據(jù)模態(tài)頻率的變化判斷裂紋的出現(xiàn)，根據(jù)振型的分析判斷裂紋的位置，根據(jù)轉(zhuǎn)子支撐系統(tǒng)阻尼的改變判斷和預(yù)報(bào)轉(zhuǎn)子的失穩(wěn)。 4） 模態(tài)分析在聲控中的應(yīng)用 聲音控制包括利用振動(dòng)和抑制振動(dòng)兩個(gè)方面。抑制結(jié)構(gòu)的輻射噪聲，在很多問(wèn)題中都很突出。模態(tài)分析為分析噪聲產(chǎn)生的原因及治理提供了有效的方法。 高速旋轉(zhuǎn)輪盤模態(tài)分析的基本理論 固有頻率 根據(jù)達(dá)朗貝爾原理 , 只要在研究對(duì)象所受的外力中加入慣性力 , 就可以像建立靜力學(xué)平衡方程那樣去建立動(dòng)力學(xué)方程 [11]。一個(gè)多自由度系統(tǒng)有阻尼振動(dòng)方程如下 : ??M ??x +??C ??x? +??K ??x =? ?)(tP （ 21） 式 中 ??M ： 系統(tǒng)質(zhì)量矩 ??C ： 系統(tǒng)阻尼矩陣 ??K ：系統(tǒng)剛度矩陣 ??x ：廣義坐標(biāo)的列陣 ??x? ：代表廣義速度 ??x ：代表廣義加速度 ? ?)(tP ：所受外部載荷 因?yàn)榻Y(jié)構(gòu)的固有頻率和振型與所受外力 ? ?)(tP 無(wú)關(guān) , 而小阻 尼對(duì)固有頻率和振型影響并不大。因此常用無(wú)阻尼自由度振動(dòng)方程求結(jié)構(gòu)的頻率和振型 , 即上式簡(jiǎn)化為下式 : 丁文文：預(yù)應(yīng)力對(duì)高速旋轉(zhuǎn)輪盤模態(tài)的影響分析 6 ??M ??x +??K ??x =0 （ 22） 由于彈性體的自由振動(dòng)總可以分解為一系列的簡(jiǎn)諧振動(dòng)的疊加。因此可設(shè)其特解 為 )sin ( ?? ?? tAx jj （ j=1， 2， … ， n） （ 23） 此特解表示系統(tǒng)內(nèi)各個(gè)坐標(biāo)偏離平衡值時(shí)均以同一頻率 ? 和同一初相角 ? 作不同振幅的簡(jiǎn)諧運(yùn)動(dòng)。上式也可以寫作矩陣形式： ? ? ? ? )sin ( ?? ?? tAx （ 24） 其 中 jAA? 為各坐標(biāo)振幅組成的 n階列陣。 將上式代入方程 （ 22） ，化作矩陣 K 和 M 的廣義本征值問(wèn)題 (??K 2? ??M )??A =0 （ 25） ??A 有非零解的的充分與必要條件為系數(shù)行列式等于零 ? ? ? ? 02 ?? MK ? （ 26） 展開后得到 2? 的 n 次代數(shù)方程，即系統(tǒng)的本征方程 2 2 ( 1 ) 211 0nn nna a a? ? ?? ?? ? ? ? ? （ 27） 對(duì)于平衡狀態(tài)穩(wěn)定的正定系統(tǒng)，各坐標(biāo)只能在平衡位置附近作微幅簡(jiǎn)諧振動(dòng)。方程（ 27） 存在 2? 的 n 個(gè)正實(shí)根，即系統(tǒng)的本征值。每個(gè)本征值所對(duì)應(yīng)的 i? ( 1,2, , )in? 為系統(tǒng)的 n 個(gè)固有頻率，將其由小到大按順序排列 為 321 ??? ?? ? n?? （ 28） 因此本征方程 （ 27） 也可稱作頻率方程。與單自由