• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    航空發(fā)動(dòng)機(jī)失諧葉盤動(dòng)態(tài)特性研究進(jìn)展*

    2016-05-30 06:59:27姚建堯王建軍
    航空制造技術(shù) 2016年21期
    關(guān)鍵詞:王建軍葉盤模態(tài)

    姚建堯,高 陽,王建軍

    (1.重慶大學(xué)航空航天學(xué)院,重慶 400044;2.中國航空工業(yè)集團(tuán)公司貴州航空發(fā)動(dòng)機(jī)研究所,貴陽 550081;3.北京航空航天大學(xué)能源與動(dòng)力工程學(xué)院,北京 100191)

    姚建堯

    重慶大學(xué)航空航天學(xué)院研究員,主要從事航空發(fā)動(dòng)機(jī)結(jié)構(gòu)強(qiáng)度、振動(dòng)和可靠性,流固耦合理論及應(yīng)用等方面的研究。

    葉盤結(jié)構(gòu)是航空發(fā)動(dòng)機(jī)中的關(guān)鍵構(gòu)件,其工作狀態(tài)直接影響到整機(jī)的結(jié)構(gòu)完整性和工作可靠性[1]。隨著結(jié)構(gòu)設(shè)計(jì)、材料和加工工藝等方面的不斷進(jìn)步,現(xiàn)代航空發(fā)動(dòng)機(jī)的壓氣機(jī)/風(fēng)扇中廣泛采用整體葉盤(環(huán)),與高效氣動(dòng)設(shè)計(jì)相結(jié)合,使壓氣機(jī)/風(fēng)扇的效率提高、流動(dòng)損失降低,發(fā)動(dòng)機(jī)整機(jī)的性能參數(shù)顯著改善。先進(jìn)的氣動(dòng)和結(jié)構(gòu)設(shè)計(jì)使得葉盤所承受的氣動(dòng)負(fù)荷增加,而結(jié)構(gòu)本身卻變得更為輕巧,因此其振動(dòng)問題也更加突出,嚴(yán)重地威脅著航空發(fā)動(dòng)機(jī)的結(jié)構(gòu)強(qiáng)度和壽命[2]。

    葉盤結(jié)構(gòu)以往的設(shè)計(jì)和分析中,通常假設(shè)各扇區(qū)是一致的。由于材料分散性、加工誤差、使用中的不均勻磨損、抑制顫振而引入錯(cuò)頻等因素,葉盤各扇區(qū)間總是不可避免地存在微小的差別,稱為失諧[3]。失諧的存在破壞了結(jié)構(gòu)的周期對稱性,某些情況下會(huì)造成明顯的振動(dòng)局部化和振動(dòng)幅值放大現(xiàn)象,使得葉盤結(jié)構(gòu)的結(jié)構(gòu)完整性和可靠性問題更為嚴(yán)峻。

    葉盤結(jié)構(gòu)的失諧問題一直受到工程界和學(xué)術(shù)界的關(guān)注,國內(nèi)外學(xué)者從理論、試驗(yàn)和工程應(yīng)用等方面都進(jìn)行了深入的研究[3-8]。近幾十年來,隨著理論和計(jì)算水平的不斷進(jìn)步,失諧葉盤的研究也經(jīng)歷了從簡化模型到高保真模型、從單級(jí)葉盤到多級(jí)葉盤、從假設(shè)失諧到真實(shí)失諧的過程,建模和計(jì)算分析精度不斷提高,為葉盤結(jié)構(gòu)設(shè)計(jì)及其安全性和可靠性的保障提供了有力的支持。本文將從建模、減縮和分析評價(jià)3個(gè)方面對失諧葉盤結(jié)構(gòu)動(dòng)態(tài)特性的國內(nèi)外相關(guān)研究進(jìn)展進(jìn)行總結(jié)。

    失諧葉盤建模技術(shù)

    合理、準(zhǔn)確地建模是葉盤動(dòng)態(tài)特性分析的基礎(chǔ),基于葉盤結(jié)構(gòu)的特點(diǎn),失諧葉盤結(jié)構(gòu)分析模型主要有集中參數(shù)模型、連續(xù)參數(shù)模型和高保真有限元模型3類[9]。在最近的研究中也報(bào)道了精確幾何失諧和多級(jí)葉盤結(jié)構(gòu)的建模技術(shù)。下面分別對其進(jìn)行介紹。

    1 集中參數(shù)模型

    集中參數(shù)模型采用集中質(zhì)量模擬葉片和輪盤的質(zhì)量,采用彈簧來模擬葉片剛度、輪盤剛度和耦合剛度,如圖1所示[10-12]。該類模型規(guī)模小,能很好地反映葉盤的一族或少數(shù)幾族諧調(diào)和失諧的動(dòng)態(tài)特性,廣泛應(yīng)用于失諧葉盤的機(jī)理研究,并取得良好的效果。

    葉盤集中參數(shù)模型的建模過程是一個(gè)典型的反問題(參數(shù)識(shí)別),需要通過已知的模態(tài)或響應(yīng)特性來獲取集中參數(shù),并使得集中參數(shù)模型的響應(yīng)特性與高保真有限元模型接近。Mignolet和Lin[13]利用最小二乘估計(jì)(Least square estimation)、最大似然估計(jì)(Maximum likelihood estimation)以及兩者的組合等3種方法建立起已知的失諧響應(yīng)和未知的集中參數(shù)的關(guān)系,并進(jìn)行求解來獲取系統(tǒng)中的參數(shù)。同樣地,也可以利用隨機(jī)模態(tài)剛度法(Radom modal stiffness approach)或最大似然估計(jì)方法,通過模態(tài)數(shù)據(jù)來獲取集中參數(shù)模型中的參數(shù)[14-15]。集中參數(shù)模型也應(yīng)用至葉盤結(jié)構(gòu)的魯棒性分析[16-17]以及多級(jí)葉盤失諧動(dòng)態(tài)特性分析中[18-19]。

    2 連續(xù)參數(shù)模型

    葉盤結(jié)構(gòu)的集中參數(shù)模型精度有限,難以模擬葉片和葉盤結(jié)構(gòu)的復(fù)雜的振動(dòng)及耦合振動(dòng)特性。為了提高動(dòng)態(tài)特性的分析精度,也有學(xué)者采用連續(xù)參數(shù)模型模擬葉盤。在這類模型中,通常采用梁模擬葉片,采用板模擬輪盤,而葉片間的耦合仍用集中參數(shù)的彈簧進(jìn)行模擬,典型結(jié)構(gòu)如圖2[20-21]所示。

    連續(xù)參數(shù)模型可以方便地考慮剪切變形、旋轉(zhuǎn)慣性、離心效應(yīng)、科氏力等[22-24]。采用該模型也可以方便地引入葉盤裂紋,考察其在旋轉(zhuǎn)狀態(tài)下的耦合振動(dòng)特性[25-27]。葉先磊等[28]采用該類模型研究了大小葉片整體葉盤的耦合振動(dòng)特性,其結(jié)果與有限元模型的計(jì)算結(jié)果具有良好的一致性。

    與集中參數(shù)模型相比,基于梁和板構(gòu)件的連續(xù)參數(shù)模型能夠提高模擬精度,但由于該模型仍難以模擬實(shí)際的葉片和輪盤的復(fù)雜振動(dòng)形式,計(jì)算分析精度仍然較低,因此也多用于失諧機(jī)理的研究。

    3 高保真有限元模型

    圖1 葉盤結(jié)構(gòu)集中參數(shù)模型Fig.1 Lumped-parameter model for bladed disk

    圖2 葉盤結(jié)構(gòu)連續(xù)參數(shù)模型Fig.2 Beam blade model for bladed disk

    前面介紹的兩種建模方法都是針對葉盤結(jié)構(gòu)動(dòng)態(tài)特性而提出的簡化模型,雖然能反映出葉盤結(jié)構(gòu)的基本振動(dòng)特征,但并不能精確描述某一特定葉盤的振動(dòng)特性,如模態(tài)特性和響應(yīng)水平等。隨著硬件和軟件水平的不斷提高,葉盤結(jié)構(gòu)的高保真有限元模型得到越來越廣泛的應(yīng)用,目前已成為通用做法。失諧葉盤有限元建模的關(guān)鍵在于如何將失諧引入至有限元模型中,目前主要通過改變?nèi)~片材料屬性和改變?nèi)~片幾何形狀來模擬和實(shí)現(xiàn)失諧。

    通過改變?nèi)~片材料屬性(通常為彈性模量)來模擬失諧是最簡便的方法,這種失諧類型通常被稱為“比例失諧(proportional mistuning)[29]”,所謂比例失諧是指失諧量與諧調(diào)情況的質(zhì)量或剛度矩陣成比例,失諧前后固有頻率變化很小,模態(tài)振型也與諧調(diào)時(shí)基本一致。這種方式能很好地反映失諧造成的各葉片間的頻率差,模擬加工誤差或使用磨損所造成的失諧,在失諧葉盤動(dòng)態(tài)特性分析中得到廣泛的應(yīng)用[30-32]。

    另一種實(shí)現(xiàn)失諧的方法是改變?nèi)~片的幾何形狀,這種方式同時(shí)改變了系統(tǒng)的質(zhì)量和剛度矩陣。在這種情況下,不僅模態(tài)頻率變化較大,而且模態(tài)振型也會(huì)發(fā)生較大的變化,因此被稱為“非比例失諧(nonproportional mistuning)[29]”。非比例失諧與實(shí)際情況更為接近,可用來模擬加工誤差、裂紋、掉角、外物打傷等所造成的失諧。Lim等[33]采用改變?nèi)~片單元位置和形狀的方法模擬了外物打傷造成的單個(gè)葉片大幾何失諧的情況(圖3[33]);Capiez-Lernout等[34-35]則通過改變?nèi)~片截面前緣和尾緣的扭轉(zhuǎn)角的方法模擬了加工誤差造成的幾何失諧;張輝有和王紅建[36]則通過改變?nèi)~片厚度來模擬葉盤的幾何失諧。

    4 幾何失諧葉盤結(jié)構(gòu)的建模

    在前面介紹的各種模型中,失諧多是通過假設(shè)的形式引入到葉盤中的。研究表明,葉盤中實(shí)際存在的幾何失諧同時(shí)改變了名義模型的剛度和質(zhì)量的分布,某些情況下其振動(dòng)特性與比例失諧葉盤有顯著差別[37],因此需在分析模型中準(zhǔn)確體現(xiàn)幾何失諧的影響。隨著測量和建模技術(shù)的進(jìn)步,近年來提出了通過葉片表面坐標(biāo)測量獲取精確幾何失諧的方法,主要有接觸式測量和非接觸式測量兩種方法。

    在接觸式測量方面,Sinha等[38]與普惠公司合作,采用坐標(biāo)測量機(jī)獲取葉片表面坐標(biāo)。通過對測量數(shù)據(jù)的恰當(dāng)正交分解(Proper Orthogonal Decomposition,POD)獲取幾何失諧的數(shù)據(jù)特征,結(jié)合CAD和CAE軟件建立有限元模型。

    非接觸式的光學(xué)測量是獲取葉片表面幾何失諧的另一種主要手段。與接觸式測量方法相比,該方法更為快捷方便,適應(yīng)性更廣,但精度要比接觸式稍低。Kaszynski等[39]實(shí)現(xiàn)了利用光學(xué)掃描高精度測量整體葉盤幾何坐標(biāo)的自動(dòng)化過程,利用測量數(shù)據(jù)重建了幾何失諧葉盤的三維模型,進(jìn)而通過有限元分析獲取其動(dòng)態(tài)特性(圖4[39])。在其后續(xù)研究中,還提出一種直接關(guān)聯(lián)名義有限元模型和幾何測量數(shù)據(jù)的網(wǎng)格變形法[40-41]以快速精確地建立幾何失諧葉盤的有限元模型。

    5 多級(jí)葉盤結(jié)構(gòu)的建模

    在多級(jí)葉盤的動(dòng)態(tài)分析中,有限元模型占有主導(dǎo)地位。Bladh等[42]建立了簡化的兩級(jí)失諧葉盤的有限元模型,研究了盤的柔性以及盤間耦合對失諧動(dòng)態(tài)特性的影響。葛長闖等[43-44]以某型發(fā)動(dòng)機(jī)兩級(jí)壓氣機(jī)葉盤為研究對象,建立有限元模型并對其自由振動(dòng)特性進(jìn)行分析,并基于應(yīng)變能定義了新的模態(tài)局部化因子。D’Souza等[45-46]研究了兩級(jí)葉盤的概率動(dòng)態(tài)特性,并首次考慮了氣動(dòng)作用對多級(jí)葉盤失諧動(dòng)態(tài)特性的影響。需要指出的是,多級(jí)葉盤的有限元模型規(guī)模一般都很大,需要對其進(jìn)行減縮處理以提高分析效率。

    葉盤結(jié)構(gòu)模型減縮技術(shù)

    圖3 失諧葉盤的有限元模型Fig.3 Finite element model for bladed disk

    圖4 非接觸式幾何失諧測量系統(tǒng)Fig.4 Optical 3D canning system for geometrical mistuning

    由于失諧破壞了葉盤結(jié)構(gòu)的周期對稱性,通常情況下須采用完整模型進(jìn)行分析,加之實(shí)際工程中的葉盤幾何形狀復(fù)雜、模型規(guī)模大,給計(jì)算帶來很大困難。特別是對于多級(jí)葉盤,即使是在諧調(diào)的情況下,由于兩級(jí)葉片數(shù)一般不存在公約數(shù),因此在分析中往往需要采用整體模型,這導(dǎo)致多級(jí)葉盤有限元模型的規(guī)模會(huì)更大。并且,失諧具有很強(qiáng)的隨機(jī)性,很多情況下需要通過蒙特卡洛模擬對多個(gè)樣本的分析以獲取其統(tǒng)計(jì)特征,進(jìn)一步增加了計(jì)算負(fù)荷。

    為此,國內(nèi)外學(xué)者研究了各種模型減縮技術(shù)和高效的求解方法,以實(shí)現(xiàn)高效、高保真的模擬和分析葉盤結(jié)構(gòu)動(dòng)態(tài)特性的目的。目前,廣泛應(yīng)用的模型減縮方法主要有基于部件的減縮、基于模態(tài)族的減縮和基于周期對稱性的減縮等3大類。

    1 基于部件的減縮技術(shù)

    基于系統(tǒng)部件的模態(tài)綜合法(Component Mode Synthesis,CMS)[47-48]是結(jié)構(gòu)動(dòng)力學(xué)中常用的子結(jié)構(gòu)模態(tài)綜合法,是將復(fù)雜的整體結(jié)構(gòu)分成較小、較易于處理的動(dòng)態(tài)子結(jié)構(gòu)的一種方法。對于葉盤結(jié)構(gòu),常見的是將葉片和輪盤分別處理為子結(jié)構(gòu),兩個(gè)子結(jié)構(gòu)之間用固定界面、自由界面或混合界面連接。這類方法在單級(jí)葉盤的失諧分析中已得到廣泛的應(yīng)用,具體可參考相關(guān)文獻(xiàn)[6]和專著[49],這里重點(diǎn)介紹該方法在幾何失諧葉盤中的擴(kuò)展應(yīng)用。

    Lim等[33,50]針對外物打傷造成的單個(gè)葉片大幾何失諧的情況進(jìn)行研究,提出了基于系統(tǒng)部件的雜交模態(tài)綜合法,將失諧葉盤分為諧調(diào)葉盤(諧調(diào)子結(jié)構(gòu))和相應(yīng)的葉片失諧部分(失諧子結(jié)構(gòu)),進(jìn)而進(jìn)行子結(jié)構(gòu)的模態(tài)減縮。

    Beck等[51-53]提出了適用于失諧葉盤動(dòng)態(tài)特性分析的Craig-Bampton部件模態(tài)綜合法(C-B CMS),該方法可將葉片與輪盤分別作為子結(jié)構(gòu)(共n+1個(gè))或者將每個(gè)扇區(qū)作為子結(jié)構(gòu)(共n個(gè))。對于比例失諧的情況,各葉片采用與諧調(diào)葉片相同的減縮基;而對于存在幾何偏差的非比例失諧的情況,各葉片則需采用不同的減縮基。在其后續(xù)的研究中,該模型減縮方法還被應(yīng)用至雙流道[54-56]葉盤的動(dòng)態(tài)特性分析中(圖5[55])。

    2 基于模態(tài)族的減縮技術(shù)

    葉盤結(jié)構(gòu)是一種典型的周期結(jié)構(gòu),即使在失諧狀態(tài)下,其模態(tài)仍具有成族出現(xiàn)的特征。利用這一特征,Yang和Griffin[57]提出了使用部分諧調(diào)模態(tài)(subset of nominal system modes,SNM)作為基向量進(jìn)行模型減縮;Feiner和Griffin[58]進(jìn)一步對基向量進(jìn)行減縮,僅使用一族模態(tài)對葉盤進(jìn)行減縮,稱為基礎(chǔ)失諧模型(Fundamental Model of Mistuning,FMM),除穩(wěn)態(tài)頻響分析外,該模型可用于失諧葉盤瞬態(tài)響應(yīng)分析[59]和失諧識(shí)別[60-62]等。在Martel等[63-64]提出的漸近失諧模型(Asymptotic Mistuning Model,AMM)中,使用的則是能夠反映系統(tǒng)動(dòng)態(tài)特性的最小模態(tài)子集,提高了這種基于模態(tài)族的減縮方法的適用性(圖6[63])。

    圖5 雙通道整體葉盤及其子結(jié)構(gòu)劃分Fig.5 Dual flow-path integrally bladed rotor and its cyclic sector C-B CMS partitionment

    前面已經(jīng)說明,幾何失諧屬于非比例失諧,使得葉片的振型發(fā)生較為顯著的變化,這就使得前面介紹的模態(tài)減縮方法的精度受到一定的影響。Ganine等[65-66]針對幾何失諧葉盤,通過靜模態(tài)補(bǔ)償(static mode compensation,SMC)對模態(tài)減縮方法進(jìn)行改進(jìn),將大的幾何失諧和小的隨機(jī)失諧采用不同的減縮基分別處理,而后再將兩者疊加。該方法可準(zhǔn)確高效地預(yù)測幾何失諧葉盤的模態(tài)和響應(yīng)等動(dòng)態(tài)特性。

    圖6 FMM和AMM失諧降階模型的模態(tài)選取Fig.6 Choices of basic vectors of FMM and AMM

    為了解決模態(tài)減縮方法在幾何失諧葉盤動(dòng)態(tài)特性分析中的精度問題,Sinha[67-68]提出一種改進(jìn)的模態(tài)域分析(modified modal domain analysis,MMDA)方法。該方法首先對葉片幾何失諧數(shù)據(jù)進(jìn)行[67]POD分析,以獲取失諧特征(圖7);進(jìn)而將具有這些幾何失諧特征的葉盤進(jìn)行諧調(diào)情況下的模態(tài)分析,與名義設(shè)計(jì)諧調(diào)葉盤的模態(tài)共同組成模態(tài)減縮基。在后續(xù)的研究中,該方法也應(yīng)用到多級(jí)幾何失諧葉盤[69]、工業(yè)葉盤[70]、失諧識(shí)別[71]、受迫響應(yīng)統(tǒng)計(jì)特性分析[72]等。

    3 基于周期對稱性的減縮技術(shù)

    葉盤結(jié)構(gòu)是一種典型的循環(huán)周期對稱結(jié)構(gòu),在單級(jí)葉盤的諧調(diào)分析中,僅可采用一個(gè)扇區(qū)進(jìn)行分析。失諧和多級(jí)葉盤雖然破壞了嚴(yán)格意義上的周期對稱性,但仍可認(rèn)為是廣義的周期對稱結(jié)構(gòu),可利用該性質(zhì)對分析模型進(jìn)行減縮。

    姚建堯等[73-74]提出了一種利用周期對稱性求解失諧葉盤瞬態(tài)響應(yīng)的計(jì)算方法。該方法將失諧部分作為非線性激勵(lì)移到運(yùn)動(dòng)方程右端,左端仍保持諧調(diào)的形式,仍可僅用單扇區(qū)進(jìn)行求解分析。失諧非線性激勵(lì)可與結(jié)構(gòu)本身的非線性(摩擦、間隙、大變形等)和非線性氣動(dòng)載荷等一并處理,提高了分析效率。王培屹和李琳[75]基于周期對稱的諧波平衡法,提出了利用單扇區(qū)模型計(jì)算失諧葉盤穩(wěn)態(tài)頻率響應(yīng)的方法,該方法的關(guān)鍵仍在于對失諧激勵(lì)的處理方式。

    周期對稱性在多級(jí)葉盤的建模中有著更為廣泛的應(yīng)用。多級(jí)葉盤周期對稱建模中最大的挑戰(zhàn)在于每級(jí)葉盤的葉片數(shù)量不同,需要對級(jí)間連接部分進(jìn)行特殊處理。Song等[76]提出了多級(jí)葉盤的逐級(jí)模態(tài)綜合法,將每一級(jí)作為一個(gè)部件,在每一級(jí)中都可僅用一個(gè)扇區(qū)進(jìn)行建模分析;在級(jí)間邊界處理上,利用級(jí)間鼓筒的軸對稱性,通過傅里葉基向量來實(shí)現(xiàn)級(jí)間的坐標(biāo)轉(zhuǎn)換,以滿足不同級(jí)間的邊界條件(圖8[76])。在后續(xù)研究中,該方法被應(yīng)用至多級(jí)葉盤的系統(tǒng)參數(shù)識(shí)別中[77]。Laxalde等[78]采用類似的方法處理級(jí)間連接鼓筒,但級(jí)間使用的是連續(xù)的網(wǎng)格。需要說明的是,上述兩個(gè)研究都是針對多級(jí)諧調(diào)葉盤進(jìn)行的。

    Bhartiya和Sinha[69]同樣將諧調(diào)多級(jí)葉盤的每一級(jí)作為子結(jié)構(gòu)分別進(jìn)行分析,在邊界處理時(shí),分別選取連接鼓筒處的約束模態(tài)和自由模態(tài)作為減縮基,并利用之前研究的MMDA方法對幾何失諧多級(jí)葉盤進(jìn)行動(dòng)態(tài)特性分析。

    失諧葉盤結(jié)構(gòu)動(dòng)態(tài)特性分析評價(jià)

    葉盤結(jié)構(gòu)中的失諧對其動(dòng)態(tài)特性產(chǎn)生非常顯著的影響,為定量衡量失諧的影響,國內(nèi)外學(xué)者定義了一系列無量綱參數(shù),如模態(tài)局部化因子、幅值放大因子、模態(tài)參與因子等。此外,人為失諧、科氏力、非線性等因素也會(huì)對失諧葉盤的動(dòng)態(tài)特性產(chǎn)生更為復(fù)雜的影響。

    1 失諧葉盤動(dòng)態(tài)特性的評價(jià)

    失諧對葉盤結(jié)構(gòu)動(dòng)態(tài)特性最為明顯的影響表現(xiàn)為模態(tài)局部化和振動(dòng)傳遞局部化,為定量描述這些現(xiàn)象,分別定義了模態(tài)局部化因子和幅值放大因子。王建軍等[79]基于實(shí)際葉盤結(jié)構(gòu)失諧振動(dòng)局部化的基本特征和物理含義,分別定義了基于模態(tài)位移、模態(tài)應(yīng)力和模態(tài)應(yīng)變能的模態(tài)局部化因子,得到評價(jià)失諧葉盤結(jié)構(gòu)模態(tài)局部化程度的定量方法。在這些定義中,都需要將失諧模態(tài)與相應(yīng)的諧調(diào)模態(tài)進(jìn)行比較,為了使兩者的對應(yīng)更具有物理意義,姚建堯等[80-81]提出振型的節(jié)徑譜的概念,通過失諧振型的主節(jié)徑成分確定與之對應(yīng)的諧調(diào)振型。

    在失諧對葉盤響應(yīng)的定量評價(jià)方面,廣泛采用的是幅值放大因子,其定義為失諧最大響應(yīng)(位移或應(yīng)力)與諧調(diào)最大響應(yīng)之比。這種定義物理意義明確,計(jì)算簡便,在工程和學(xué)術(shù)界得到廣泛應(yīng)用[18,32,72,82-84]。

    圖7 葉片幾何失諧的POD特征Fig.7 Blade thickness for each POD feature

    圖8 多級(jí)葉盤的周期對稱模型Fig.8 Cyclic model for multistage bladed disk

    與單級(jí)葉盤相比,多級(jí)葉盤模態(tài)局部化特性的評價(jià)要更為復(fù)雜,主要需考慮級(jí)間耦合的問題。Sternchuss和Balmes[85]通過應(yīng)變能定義了各級(jí)的模態(tài)參與因子(relative participation factor),以衡量各階模態(tài)是以單級(jí)振動(dòng)為主導(dǎo)還是以耦合振動(dòng)為主導(dǎo)。D’Souza和Epureanu[86]則分別利用應(yīng)變能比(strain energy ratio)和模態(tài)置信因子對多級(jí)葉盤的模態(tài)進(jìn)行分類評價(jià)。多級(jí)葉盤響應(yīng)評價(jià)相對簡潔,一般仍采用與單級(jí)葉盤的相同的幅值放大因子[69]。

    此外,Rotea和D’Amato[11]和姚建堯等[16-17]也試圖從系統(tǒng)魯棒性的角度評價(jià)失諧葉盤的動(dòng)態(tài)特性,但目前只能應(yīng)用于簡單的集中參數(shù)模型。Nikolic等[87]提出了通過增加失諧范圍來提高葉盤結(jié)構(gòu)魯棒性的設(shè)計(jì)理念;Mbaye等[88]提出利用錯(cuò)頻手段提高葉盤結(jié)構(gòu)的魯棒性。

    2 人為失諧的影響

    Kaza和Kielb[89-90]的研究表明,適當(dāng)程度的人為失諧可以有效提高壓氣機(jī)的顫振裕度。目前,引入人為失諧或錯(cuò)頻以抑制顫振已成為發(fā)動(dòng)機(jī)設(shè)計(jì)的常用技術(shù)手段[91-93]。同時(shí),某些人為失諧也能有效地抑制隨機(jī)失諧對葉盤響應(yīng)的不利影響。Choi等[94]研究了AB兩種葉片類型人為失諧排列方式對隨機(jī)失諧響應(yīng)放大的抑制效果,并應(yīng)用至離心葉盤的設(shè)計(jì)中。Kenyon和Griffin[82]研究了更為復(fù)雜的諧波形式的人為失諧對隨機(jī)失諧葉盤響應(yīng)的影響。Murthy和Mignolet[95]研究了干摩擦阻尼器人為失諧對葉盤的影響。姚建堯等[80-81]則利用節(jié)徑譜和模態(tài)激勵(lì)因子的概念,很好地解釋了隨機(jī)失諧的閾值效應(yīng)和人為失諧的作用機(jī)理。袁惠群[96]、李巖[97]和趙天宇等[98]分別采用蟻群算法、離散遺傳粒子群算法和模擬退火算法對人為失諧排列方式進(jìn)行優(yōu)化。

    3 科氏力的影響

    葉盤在高速旋轉(zhuǎn)的工作狀態(tài)下,科氏力會(huì)對其動(dòng)態(tài)特性產(chǎn)生較為顯著的影響。李永強(qiáng)等[99]基于薄殼理論考慮了科氏力對旋轉(zhuǎn)葉片動(dòng)頻的影響,但并未考慮對輪盤的影響。Huang和Kuang[24]采用連續(xù)參數(shù)模型研究了考慮科氏力效應(yīng)下的失諧葉盤的模態(tài)局部化,結(jié)果表明,科氏力可能會(huì)加劇失諧葉盤的模態(tài)局部化現(xiàn)象。Nikolic等[100]研究了科氏力和失諧共同作用下葉盤的響應(yīng)特性,并進(jìn)行了相應(yīng)的試驗(yàn)驗(yàn)證,結(jié)果表明,某些情況下科氏力使得失諧響應(yīng)顯著增大。辛健強(qiáng)等[101]采用有限元法對考慮科氏力的工業(yè)葉盤進(jìn)行分析,結(jié)果同樣表明,科氏力對失諧葉盤的模態(tài)和響應(yīng)特性會(huì)產(chǎn)生顯著的影響。

    4 非線性的影響

    在失諧葉盤動(dòng)態(tài)特性分析中,也考慮了摩擦、裂紋、大變形等非線性因素的影響。

    干摩擦阻尼器是航空發(fā)動(dòng)機(jī)中最常見的阻尼器之一,其對失諧葉盤動(dòng)態(tài)特性,特別是響應(yīng)水平有著顯著的影響。Poudou和Pierre[102]提出了帶干摩擦阻尼器的葉盤結(jié)構(gòu)響應(yīng)計(jì)算的頻域-時(shí)域混合方法。Cha和Sinha[103]研究了白噪聲和窄帶激勵(lì)下帶摩擦阻尼器的失諧葉盤響應(yīng)的統(tǒng)計(jì)特性。王紅建[104]和賀爾銘等[105]研究了干摩擦阻尼失諧對葉盤受迫響應(yīng)的影響,結(jié)果表明,干摩擦散亂失諧會(huì)導(dǎo)致阻尼件減振性能的降低。Petrov等[106-111]系統(tǒng)研究了考慮摩擦和間隙等非線性因素情況下的失諧葉盤的建模和分析技術(shù)。

    裂紋是葉片中常見的缺陷形式,同樣是引起失諧的主要因素之一。Kuang和Huang[25]以及Hou[27]通過連續(xù)參數(shù)模型考慮了葉片裂紋位置以及大小對失諧葉盤模態(tài)局部化的影響,但并未考慮裂紋的非線性性質(zhì)。Saito等[112]建立了裂紋葉片的有限元模型,并考慮了裂紋開合造成的非線性,基于模態(tài)綜合法建立了裂紋失諧葉盤的降階模型;在其后續(xù)研究中,還將其擴(kuò)展至多級(jí)葉盤[113]和裂紋檢測中[114]。Wang等[115]利用類似的方法對失諧離心葉盤進(jìn)行了分析。

    大變形和大位移所引起的幾何非線性對失諧葉盤有一定的影響。Capiez-lernout等[116-119]研究了考慮幾何非線性的高負(fù)荷離心壓氣機(jī)葉盤的建模、降階和響應(yīng)分析技術(shù)。

    結(jié)論和展望

    本文重點(diǎn)回顧了失諧葉盤的建模、降階以及動(dòng)態(tài)特性分析評價(jià)等方面的相關(guān)研究進(jìn)展,著重對多級(jí)葉盤和幾何失諧葉盤的最新研究進(jìn)行了介紹。需要指出的是,由于論文篇幅有限,這里并未涉及葉盤結(jié)構(gòu)流固耦合、失諧識(shí)別、試驗(yàn)技術(shù)等重要領(lǐng)域。

    從前述的研究進(jìn)展可以看出,隨著計(jì)算機(jī)軟硬件和建模技術(shù)的不斷提高,失諧葉盤高保真模型已得到廣泛的應(yīng)用,特別是基于坐標(biāo)測量的精細(xì)建模技術(shù)也已于近年來開始應(yīng)用。但由于失諧葉盤的隨機(jī)性,更為關(guān)注其動(dòng)態(tài)特性的統(tǒng)計(jì)特征,因此高效精確的降階技術(shù)仍是研究的重要方向之一。為了更加準(zhǔn)確地預(yù)測失諧響應(yīng)及其概率分布特征,進(jìn)行失諧葉盤流固耦合的概率特性分析也是未來重要的研究方向之一。

    除理論研究外,如何將失諧分析和設(shè)計(jì)的理念在工程中加以應(yīng)用,指導(dǎo)航空發(fā)動(dòng)機(jī)的設(shè)計(jì)、評估和使用是未來需要重點(diǎn)關(guān)注的問題。

    [1]李其漢,王延榮.航空發(fā)動(dòng)機(jī)結(jié)構(gòu)強(qiáng)度設(shè)計(jì)問題[M].上海: 上海交通大學(xué)出版社,2014.

    LI Qihan,WANG Yanrong.The design problem of aero-engine structure strength[M].Shanghai: Shanghai Jiao Tong University Press,2014.

    [2]李其漢,王延榮,王建軍.航空發(fā)動(dòng)機(jī)葉片高循環(huán)疲勞失效研究[J].航空發(fā)動(dòng)機(jī),2003,29(4): 16-18.

    LI Qihan,WANG Yanrong,WANG Jianjun.Investigation of high cycle fatigue failures for the aero engine blades[J].Aeroengine,2003,29(4):16-18.

    [3]王建軍,李其漢,朱梓根.失諧葉片-輪盤結(jié)構(gòu)系統(tǒng)振動(dòng)局部化問題的研究進(jìn)展[J].力學(xué)進(jìn)展,2000,30(4): 517-528.

    WANG JianJun,LI Qihan,ZHU ZiGen.Vibratory localization of mistuned bladed disk assemblies—a review[J].Advances in Mechanics,2000,30(4): 517-528.

    [4]SRINIVASAN A V.Flutter and resonant vibration characteristics of engine blades[J].Journal of Engineering for Gas Turbines and Power,1997,119(10): 742-775.

    [5]SLATER J C,MINKIEWICZ G R,BLAIR A J.Forced response of bladed disk assemblies-a survey[J].Shock and Vibration Digest,1999,31(1): 17-24.

    [6]CASTANIER M P,PIERRE C.Modeling and analysis of mistuned bladed disk vibration: status and emerging directions[J].Journal of Propulsion and Power,2006,22(2):384-396.

    [7]臧朝平,蘭海強(qiáng).失諧葉盤結(jié)構(gòu)振動(dòng)問題研究新進(jìn)展[J].航空工程進(jìn)展, 2011,2(2): 133-142.

    ZANG Chaoping,LAN Haiqiang.Advances in research vibration problem of mistuned blisk assemblies[J].Advances in Aeronautical Science and Engineering,2011,2(2): 133-142.

    [8]白斌,白廣忱,童曉晨,等.整體葉盤結(jié)構(gòu)失諧振動(dòng)的國內(nèi)外研究狀況[J].航空動(dòng)力學(xué)報(bào),2014,29(1): 91-103.

    BAI Bin,BAI Guangchen,TONG Xiaochen,et al.Research on vibration problem of integral mistuned bladed disk assemblies at home and abroad[J].Journal of Aerospace Power,2014,29(1): 91-103.

    [9]王建軍,許建東,李其漢.失諧葉片-輪盤結(jié)構(gòu)振動(dòng)局部化的分析模型[J].汽輪機(jī)技術(shù), 2004,46(4): 256-259.

    WANG Jianjun,XU Jiandong,LI Qihan.Analytical models of mistuned bladed disk assembles—a review[J].Turbine Technology,2004,46(4): 256-259.

    [10]CASTANIER M P,PIERRE C.Using Intentional mistuning in the design of turbomachinery rotors[J].AIAA Journal,2002,40(10): 2077-2086.

    [11]ROTEA M,D’AMATO F.New tools for analysis and optimization of mistuned bladed disks[C]//38th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit,2002.

    [12]王建軍,姚建堯,李其漢.剛度隨機(jī)失諧葉盤結(jié)構(gòu)概率模態(tài)特性分析[J].航空動(dòng)力學(xué)報(bào),2008,23(2): 256-262.

    WANG Jianjun,YAO Jianyao,LI Qihan.Probability characteristics of vibratory mode of bladed disk assemblies with random stiffness mistuning[J].Journal of Aerospace Power,2008,23(2): 256-262.

    [13]MIGNOLET M P,LIN C.Identification of structural parameters in mistuned bladed disks[J].Journal of Vibration and Acoustics,1997,119(3): 428-438.

    [14]MIGNOLET M P,RIVAS-GUERRA A J,DELOR J P.Identification of mistuning characteristics of bladed disks from free response data—part I[J].Journal of Engineering for Gas Turbines and Power,2001,123(2): 395-403.

    [15]RIVAS-GUERRA A J,MIGNOLET M P,DELOR J P,et al.Identification of mistuning characteristics of bladed disks from free response data—part II[J].Journal of Engineering for Gas Turbines and Power,2001,123(2): 404-411.

    [16]姚建堯,王建軍,李其漢.失諧周期結(jié)構(gòu)的魯棒分析模型[J].應(yīng)用力學(xué)學(xué)報(bào),2009,26(3): 426-431.

    YAO Jianyao,WANG Jianjun,LI Qihan.Robustness analysis model for mistuned periodic structures[J].Chinese Journal of Applied Mechanics,2009,26(3): 426-431.

    [17]姚建堯,王建軍,李其漢.失諧葉盤結(jié)構(gòu)魯棒性能分析[J].航空動(dòng)力學(xué)報(bào),2010,25(7): 1634-1639.

    YAO Jianyao,WANG Jianjun,LI Qihan.Robust performance analysis of mistuned bladed disks[J].Journal of Aerospace Power,2010,25(7):1634-1639.

    [18]廖海濤,王建軍,李其漢.多級(jí)葉盤結(jié)構(gòu)隨機(jī)失諧響應(yīng)特性分析[J].振動(dòng)與沖擊,2011,30(3): 22-29.

    LIAO Haitao,WANG Jianjun,LI Qihan.Mistuned forced response characteristics analysis of mistuned multi-stages bladed disks[J].Journal of Vibration And Shock,2011,30(3): 22-29.

    [19]SINHA A.Reduced-order model of a mistuned multi-stage bladed rotor[C]//Proceedings of GT2007,ASME Turbo Expo 2007:Power for Land,Sea and Air,2007.

    [20]TURCOTTE J S,HOLLKAMP J J,GORDON R W.Vibration of a mistuned bladed-disk assembly using structurally damped beams[J].AIAA Journal,1998,36(12): 2225-2228.

    [21]TOMIOKA T,KOBAYASHI Y,YAMADA G.Analysis of free vibration of rotating disk–blade coupled systems by using artificial springs and orthogonal polynomials[J].Journal of Sound and Vibration,1996,191(1): 53-73.

    [22]RZADKOWSKI R.The general model of free vibrations of mistuned bladed discs,part i:Theory[J].Journal of Sound and Vibration,1994,173(3): 377-393.

    [23]RZADKOWSKI R.The general model of free vibrations of mistuned bladed discs,part ii: Numerical results[J].Journal of Sound and Vibration,1994,173(3): 395-413.

    [24]HUANG B W,KUANG J H.Mode locallization in a rotating mistuned turbo disk with Coriolis effects[J].International Journal of Mechanical Sciences,2001,43: 1643-1660.

    [25]KUANG J H,HUANG B W.The effect of blade crack on mode localization in rotating bladed disks[J].Journal of Sound and Vibration,1999,227(1): 85-103.

    [26]HUANG B W,KUNG H K,KUANG J H.Stability in a twisted periodic blade system with cracks[J].AIAA Journal,2006,44(7): 1436-1444.

    [27]HOU J.Cracking-induced mistuning in bladed disks[J].AIAA Journal,2006,44(11):2542-2546.

    [28]葉先磊,王建軍,朱梓根,等.大小葉片結(jié)構(gòu)連續(xù)參數(shù)模型和振動(dòng)模態(tài)[J].航空動(dòng)力學(xué)報(bào),2005,20(1): 66-72.

    YE Xianlei,WANG Jianjun,ZHU Zigen,et al.Continuous parameter model and vibration modal of splitter vane rotor[J].Journal of Aerospace Power,2005,20(1): 66-72.

    [29]LIM S H,BLADH R,CASTANIER M P,et al.A compact,generalized component mode mistuning representation for mode bladed disk vibration[C]//44th AIAA/ASME/ASCE/AHS Structures,Structural Dynamics,and Materials Conference.2003.

    [30]BLADH R,PIERRE C,CASTANIER M P.Dynamic response predictions for a mistuned industrial turbomachinery rotor using reducedorder modeling[J].Journal of Engineering for Gas Turbines and Power,2002,124(4): 311-324.

    [31]于長波,王建軍,李其漢.錯(cuò)頻葉盤結(jié)構(gòu)的概率模態(tài)局部化特性分析[J].航空動(dòng)力學(xué)報(bào),2009,24(9): 2040-2045.

    YU Changbo,WANG Jianjun,LI Qihan.Probability characteristics for vibratory mode of detuned bladed disk assemblies[J].Journal of Aerospace Power,2009,24(9): 2040-2045.

    [32]于長波,王建軍,李其漢.失諧葉盤結(jié)構(gòu)的概率響應(yīng)局部化特性[J].航空動(dòng)力學(xué)報(bào),2010,25(9): 2006-2012.

    YU Changbo,WANG Jianjun,LI Qihan.Probability characteristics for response localization of mistuned bladed disk assemblies[J].Journal of Aerospace Power,2009,25(9): 2006-2012.

    [33]LIM S H,BLADH R,CASTANIER M P,et al.Compact,generalized component mode mistuning representation for modeling bladed disk vibration[J].AIAA Journal,2007,45(9): 2285-2298.

    [34]CAPIEZ-LERNOUT E,SOIZE C.Nonparametric modeling of random uncertainties for dynamic response of mistuned balded disks[J].Journal of Engineering for Gas Turbines and Power,2004,126(7): 610-618.

    [35]CAPIEZ-LERNOUT E,SOIZE C,LOMBARD J P,et al.Blade manufacturing tolerances definition for a mistuned industrial bladed disk[J].Journal of Engineering for Gas Turbines and Power,2005,127(7): 621-628.

    [36]張輝有,王紅建.一種基于葉盤結(jié)構(gòu)幾何失諧的降階分析方法[J].航空工程進(jìn)展,2014,5(4): 481-486.

    ZHANG Huiyou,WANG Hongjian.A reduced order model of bladed disks with geometric mistuning[J].Advances in Aeronautical Science and Engineering,2014,5(4): 481-486.

    [37]BECK J A,BROWN J M,SLATER J C,et al.Probabilistic mistuning assessment using nominal and geometry based mistuning methods[J].Journal of Turbomachinery,2013,135(5): 51004.

    [38]SINHA A,HALL B,CASSENTI B,et al.Vibratory parameters of blades from coordinate measurement machine data[J].Journal of Turbomachinery,2008,130(1): 523-531.

    [39]KASZYNSKI A A,BECK J A,BROWN J M.Uncertainties of an automated optical 3d geometry measurement,modeling,and analysis process for mistuned integrally bladed rotor reverse engineering[J].Journal of Engineering for Gas Turbines and Power,2013,135(10): 1025-1034.

    [40]KASZYNSKI A A,BECK J A,BROWN J M.Automated finite element model mesh updating scheme applicable to mistuning analysis[C]//ASME Turbo Expo 2014: Turbine Technical Conference and Exposition.American Society of Mechanical Engineers,2014.

    [41]KASZYNSKI A A,BECK J A,BROWN J M.Experimental validation of a mesh quality optimized morphed geometric mistuning model[C]//ASME Turbo Expo 2015: Turbine Technical Conference and Exposition.2015.

    [42]BLADH R,CASTANIER M P,PIERRE C.Effects of multi-stage coupling and disk flexibility on mistuned bladed disk dynamics[J].Journal of Engineering for Gas Turbines & Power,2001,125(1): V004T033A040.

    [43]葛長闖,王建軍,劉永泉.2級(jí)葉片-輪盤系統(tǒng)模態(tài)特性研究[J].航空發(fā)動(dòng)機(jī),2009,35(5): 19-23.

    GE Changchuang,WANG Jianjun,LIU Yongquan.Investigation of characteristics of twostage blade-disk system modal[J].Aeroengine,2009,35(5): 19-23.

    [44]葛長闖,王建軍,劉永泉.失諧多級(jí)整體葉盤振動(dòng)模態(tài)特性定量評價(jià)方法研究[J].航空發(fā)動(dòng)機(jī),2012,38(1): 25-28.

    GE Changchuang,WANG Jianjun,LIU Yongquan.Quantitative assessment method of vibration mode characteristics for mistuned multistage blisk[J].Aeroengine,2012,38(1): 25-28.

    [45]D’SOUZA K X,EPUREANU B I.A statistical characterization of the effects of mistuning in multistage bladed disks[J].Journal of Engineering for Gas Turbines and Power,2012,134(1): n37-n47.

    [46]D’SOUZA K X,JUNG C,EPUREANU B I.Analyzing mistuned multistage turbomachinery rotors with aerodynamic effects[J].Journal of Fluids and Structures,2013,42: 388-400.

    [47]CRAIG R,CHANG C.Free-interface methods of substructure coupling for dynamic analysis[J].AIAA Journal,1976,14(11): 1633-1635.

    [48]CRAIG R R.Coupling of substructures for dynamic analyses: an overview[C]// Proceedings of AIAA/ASME/ASCE/AHS/ASC Structures,Structural Dynamics,and Materials Conference and Exhibit.2000.

    [49]王建軍,李其漢.航空發(fā)動(dòng)機(jī)失諧葉盤振動(dòng)減縮模型與應(yīng)用[M].北京: 國防工業(yè)出版社,2009.

    WANG Jianjun,LI Qihan.Methods and applications of reduction modeling for mistuned bladed disk vibration in aero-engine[M].Beijing:National Defend Industry Press ,2009.

    [50]LIM S,CASTANIER M,PIERRE C.Vibration modeling of bladed disks subject to geometric mistuning and design changes[C]//45th AIAA/ASME/ASCE/AHS/ASC Structures,Structural Dynamics and Material Comference.2004.

    [51]BECK J A,BROWN J M,SLATER J C,et al.Probabilistic mistuning assessment using nominal and geometry based mistuning methods[J].Journal of Turbomachinery,2012,135(5): 1085-1097.

    [52]BECK J A,BROWN J M,CROSS C J,et al.Component-mode reduced-order models for geometric mistuning of integrally bladed rotors[J].AIAA Journal,2014,52(7): 1345-1356.

    [53]BECK J A,BROWN J M,KASZYNSKI A A,et al.Geometric mistuning reduced-order models for integrally bladed rotors with mistuned disk–blade boundaries[J].Journal of Turbomachinery,2015,137(7): 71001.

    [54]BECK J A,BROWN J M,KASZYNSKI A A,et al.Mistuned response prediction of dual flow-path integrally bladed rotors with geometric mistuning[J].Journal of Engineering for Gas Turbines and Power,2015,137(6): 62501.

    [55]BECK J A,BROWN J M,SCOTTEMUAKPOR O E,et al.Dynamic response characteristics of dual flow-path integrally bladed rotors[J].Journal of Sound and Vibration,2015,336: 150-163.

    [56]BECK J A,SCOTT-EMUAKPOR O E,BROWN J M,et al.Validation of geometric mistuning reduced-order models for single and dual flow-path integrally bladed rotors[C]//17th AIAA Non-Deterministic Approaches Conference.2015.

    [57]YANG M T,GRIFFIN J H.A reduced-order model of mistuning using a subset of nominal system modes[J].Journal of Engineering for Gas Turbines and Power,2001,123(10): 893-900.

    [58]FEINER D M,GRIFFIN J H.A fundamental model of mistuning for a single family of modes[J].Journal of Turbomachinery,2002,124(10): 597-605.

    [59]AYERS J P,FEINER D M,GRIFFIN J H.A reduced-order model for transient analysis of bladed disk forced response[J].Journal of Turbomachinery,2006,128(7): 466-473.

    [60]FEINER D M,GRIFFIN J H.Mistunign identification of bladed disks using a fundamental mistuning model—part ii:application[J].Journal of Turbomachinery,2004,126(1): 159-165.

    [61]FEINER D M,GRIFFIN J H.Mistuning identification of bladed disks using a fundamental mistuning model—part i: theory[J].Journal of Turbomachinery,2004,126(1): 150-158.

    [62]王帥,王建軍,李其漢.一種基于模態(tài)減縮技術(shù)的整體葉盤結(jié)構(gòu)失諧識(shí)別方法[J].航空動(dòng)力學(xué)報(bào),2009,24(3): 662-669.

    WANG Shuai,WANG Jianjun,LI Qihan.Mistuning identification of integrally bladed disk based on the modal reduced technique[J].Journal of Aerospace Power,2009,24(3): 662-669.

    [63]MARTEL C,CORRAL R.Asymptotic description of maximum mistuning amplification of bladed disk forced response[J].Journal of Engineering for Gas Turbines and Power,2009,131(2): 435-446.

    [64]KHEMIRI O,MARTEL C,CORRAL R.Quantitative validation of the asymptotic mistuning model using a high fidelity bladed disk model[C]// ASME Turbo Expo 2010: Power for Land,Sea,and Air.2010.

    [65]GANINE V,LEGRAND M,PIERRE C,et al.A reduction technique for mistuned bladed disks with superposition of large geometric mistuning and small model uncertainties[C]//Proceedings of the 12th International Symposium on Transport Phenomena and Dynamics of Rotating Machinery.2008.

    [66]GANINE V,LEGRAND M,MICHALSKA H,et al.A sparse preconditioned iterative method for vibration analysis of geometrically mistuned bladed disks[J].Computers & Structures,2009,87(5): 342-354.

    [67]SINHA A.Reduced-order model of a bladed rotor with geometric mistuning[C]//Proceedings of ASME Turbo Expo 2007: Power for Land,Sea and Air.2007.

    [68]SINHA A.Reduced-order model of a bladed rotor with geometric mistuning[J].Journal of Turbomachinery,2009,131(3): 310-314.

    [69]BHARTIYA Y,SINHA A.Reduced order model of a multistage bladed rotor with geometric mistuning via modal analyses of finite element sectors[J].Journal of Turbomachinery,2012,134(4): 410-411.

    [70]VISHWAKARMA V,SINHA A,BHARTIYA Y,et al.Modified modal domain analysis of a bladed rotor using coordinate measurement machine data on geometric mistuning[J].Journal of Engineering for Gas Turbines and Power,2015,137(4): 110-121.

    [71]BHARTIYA Y,SINHA A.Geometric mistuning identification of integrally bladed rotors using modified modal domain analysis[J].Journal of Engineering for Gas Turbines and Power,2014,136(12): 205-215.

    [72]VISHWAKARMA V,SINHA A.Forced response statistics of a bladed rotor with geometric mistuning[J].AIAA Journal,2015,53(9): 2776-2781.

    [73]姚建堯,王建軍.周期對稱性在葉盤結(jié)構(gòu)瞬態(tài)響應(yīng)求解中的應(yīng)用[J].航空動(dòng)力學(xué)報(bào),2011,26(2): 385-391.

    YAO Jianyao,WANG Jianjun.Application of cyclic symmetry property in transient forced vibration of bladed disks[J].Journal of Aerospace Power,2011,26(2): 385-391.

    [74]YAO J Y,ZHU W X,HU N,et al.Application of symmetry property for transient response analysis of mistuned bladed disks[C]//52nd AIAA/SAE/ASEE Joint Propulsion Conference.2016.

    [75]王培屹,李琳.用于失諧葉盤動(dòng)力學(xué)特性分析的減縮計(jì)算方法[J].航空動(dòng)力學(xué)報(bào),2014,29(6): 1395-1402.

    WANG Peiyi,LI Lin.Reduced order computational method for analysis of mistuning bladed disk dynamics characteristic[J].Journal of Aerospace Power,2014,29(6): 1395-1402.

    [76]SONG S H,CASTANIER M P,PIERRE C.Multi-stage modeling of turbine engine rotor vibration[C]// ASME 2005 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference.2005.

    [77]SONG S H,CASTANIER M P,PIERRE C.System identification of multistage turbine engine rotors[C]// Proceedings of GT2007,ASME Turbo Expo 2007: Power for Land,Sea and Air.2007.

    [78]LAXALDE D,LOMBARD J P,THOUVEREZ F.Dynamics of multistage bladed disks systems[J].Journal of Engineering for Gas Turbines and Power,2007,129(5): 1058-1064.

    [79]王建軍,于長波,姚建堯,等.失諧葉盤振動(dòng)模態(tài)局部化定量描述方法[J].推進(jìn)技術(shù),2009,30(4): 457-461.

    WANG Jianjun,YU Changbo,YAO Jianyao,et al.Vibratory mode localization factors of mistuned bladed disk assemblies[J].Journal of Propulsion Technology,2009,30(4): 457-461.

    [80]姚建堯,王建軍,李其漢.基于振型節(jié)徑譜的失諧葉盤結(jié)構(gòu)動(dòng)態(tài)特性評價(jià)[J].推進(jìn)技術(shù),2011,32(5): 645-653.

    YAO Jianyao,WANG Jianjun,LI Qihan.Dynamic characteristics assessment of mistuned bladed disk using nodal diameter spectrum of mode shapes[J].Journal of Propulsion Technology,2011,32(5): 645-653.

    [81]YAO J Y,WANG J J,Li Q H.Improved modal localization and excitation factors for understanding mistuned bladed disk response[J].Journal of Propulsion and Power,2011,27(1): 50-60.

    [82]KENYON J A,GRIFFIN J H.Forced response of turbine engine bladed disks and sensitivity to harmonic mistuning[J].Journal of Engineering for Gas Turbines and Power,2003,125(1): 113-120.

    [83]LIAO H T,WANG J J,YAO J Y,et al.Mistuning forced response characteristics analysis of mistuned bladed disks[J].Journal of Engineering for Gas Turbines and Power,2010,132(12): 1113-1122.

    [84]KHEMIRI O,MARTEL C,CORRAL R.Forced response of mistuned bladed disks:quantitative validation of the asymptotic description[J].Journal of Propulsion and Power,2014,30(2): 397-406.

    [85]STERNCHUSS A,BALMES E.Reduction of multistage rotor models using cyclic modeshapes[C]// Proceedings of GT2007,ASME Turbo Expo 2007: Power for Land,Sea and Air.2007.

    [86]D’SOUZA K X,EPUREANU B I.A statistical characterization of the effects of mistuning in multistage bladed disks[J].Journal of Engineering for Gas Turbines and Power,2012,134(1): 1137-1147.

    [87]NIKOLIC M,PETROV E P,EWINS D J.Robust strategies for forced response reduction of bladed discs on large mistuning concept[C]// Proceedings of GT2007,ASME Turbo Expo 2007: Power for Land,Sea and Air.2007.

    [88]MBAYE M,SOIZE C,OUSTY J,et al.Robust analysis of design in vibration of turbomachines[J].Journal of Turbomachinery,2013,135(2): 210-218.

    [89]KAZA K R V,KIELB R E.Flutter and response of a mistuned cascade in incompressible flow[J].AIAA Journal,1982,20(8): 1120-1127.

    [90]KAZA K R V,KIELB R E.Flutter of turbofan rotors with mistuned blades[J].AIAA Journal,1984,22(11): 1618-1625.

    [91]徐建民,宋兆泓.錯(cuò)頻葉片轉(zhuǎn)子的穩(wěn)定性分析[J].北京航空航天大學(xué)學(xué)報(bào),1986(4): 21-32.

    XU Jianmin,SONG Zhaohong.The stability analysis of the rotor of the fault frequency[J].Journal of Beijing University of Aeronautics and Astronautics,1986(4): 21-32.

    [92]鄭赟,王靜.錯(cuò)頻對葉片的氣動(dòng)彈性穩(wěn)定性影響[J].航空動(dòng)力學(xué)報(bào),2013,28(5):1029-1036.

    ZHENG Yun,WANG Jing.Influence of frequency mistuning on aeroelastic stability of blade[J].Journal of Aerospace Power,2013,28(5):1029-1036.

    [93]徐可寧,王延榮,劉金龍.壓氣機(jī)轉(zhuǎn)子錯(cuò)頻葉盤結(jié)構(gòu)振動(dòng)響應(yīng)分析[J].燃?xì)鉁u輪試驗(yàn)與研究,2013,26(3): 6-12.

    XU Kening,WANG Yanrong,LIU Jinlong.Vibration response analysis on mistuned bladed disk assembly of compressor rotors[J].Gas Turbine Experiment and Research,2013,26(3): 6-12.

    [94]CHOI B K,LENTZ J,MIGNOLET M P.Optimization of intentional mistuning patterns for the reduction of the forced response effects of unintentional mistuning: formulation and assessment[J].Antipode,2002,16(3): 47-59.

    [95]MURTHY R,MIGNOLET M P.On the benefits of intentional mistuning of friction dampers to reduce the response of tuned and mistuned bladed disks[C]// ASME Turbo Expo 2013: Turbine Technical Conference and Exposition.2013.

    [96]袁惠群,張亮,韓清凱,等.基于蟻群算法的航空發(fā)動(dòng)機(jī)失諧葉片減振排布優(yōu)化分析[J].振動(dòng)與沖擊,2012,31(11): 169-172.

    YUAN Huiqun,ZHANG Liang,HAN Qingkai,et al.Optimization of mistuning blades arrangement for vibration absorption in an aeroengine based on artificial ant colony algorithm[J].Journal of Vibration and Shock,2012,31(11):169-172.

    [97]李巖,袁惠群,梁明軒.基于改進(jìn)DPSO算法的航空發(fā)動(dòng)機(jī)失諧葉片排序[J].東北大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,33(4): 569-572.

    LI Yan,YUAN Huiqun,LIANG Mingxuan.Mistuned blade sorting based on improved dpso algorithm for aero-engine[J].Journal of Northeastern University(Natural Science),2013,33(4): 569-572.

    [98]趙天宇,袁惠群,楊文軍,等.非線性摩擦失諧葉片排序并行退火算法[J].航空動(dòng)力學(xué)報(bào),2016,31(5): 1053-1064.

    ZHAO Tianyu,YUAN Huiqun,YANG Wenjun.Mistuned blade arrangement with non-linear friction based on parallel annealing algorithm[J].Journal of Aerospace Power,2016,31(5): 1053-1064.

    [99]李永強(qiáng),李健,郭星輝.科氏力對旋轉(zhuǎn)葉片動(dòng)頻的影響[J].振動(dòng)與沖擊,2006,25(1): 79-83.

    LI Yongqiang,LI Jian,GUO Xinghui.Effect of coriolis force on dynamic frequency of high speed spinning blade[J].Journal of Vibration and Shock,2006,25(1): 79-83.

    [100]NIKOLIC M,PETROV E P,EWINS D J.Coriolis forces in forced response analysis of mistuned bladed disks[J].Journal of Turbomachinery,2007,129(10): 730-739.

    [101]辛健強(qiáng),王建軍,李其漢.科氏力對失諧葉盤振動(dòng)特性的影響分析[J].推進(jìn)技術(shù),2011,32(5): 637-644.

    XIN Jianqiang,WANG Jianjun,LI Qihan.Analysis on vibration characteristics of mistuned bladed disk assemblies with Coriolis force[J].Journal of Propulsion Technology,2011,32(5):637-644.

    [102]POUDOU O,PIERRE C.Hybrid frequency-time domain methods forthe analysis of complex structural systems with dry friction damping[C]//44th AIAA/ASME/ASCE/AHS Structures,Structural Dynamics,and Materials Conference.2003.

    [103]CHA D,SINHA A.Statistics of responses of a mistuned and frictionally damped bladed disk assembly subjected to white noise narrow band excitations[J].Probabilistic Engineering Mechanics,2006,21: 384-396.

    [104]王紅建,賀爾銘,余仕俠.具有干摩擦散亂失諧的葉盤受迫響應(yīng)特性[J].航空動(dòng)力學(xué)報(bào),2006,21(4): 711-715.

    WANG Hongjian,HE Erming,YU Shixia.Forced response characteristics of bladed disks with disordered dry friction[J].Journal of Aerospace Power,2006,21(4): 711-715.

    [105]賀爾銘,王紅建,余仕俠.含有非線性摩擦阻尼的失諧葉盤系統(tǒng)受迫響應(yīng)研[J].機(jī)械強(qiáng)度,2007,29(4): 532-535.

    HE Erming,WANG Hongjian,YU Shixia.Investigation of forced response of mistuned bladed disks with non-linear friction damping[J].Journal of Mechanical Strength,2007,29(4):532-535.

    [106]PETROV E P.A method for use of cyclic symmetry properties in analysis of nonlinear multiharmonic vibrations of bladed disks[J].Journal of Turbomach,2004,126(1):175-183.

    [107]PETROV E P,EWINS D J.Generic friction models for time-domain vibration analysis of bladed disks[J].Journal of Turbomachinery,2004,126(1): 184-192.

    [108]PETROV E P,EWINS D J.Methods for analysis of nonlinear multiharmonic vibrations of mistuned bladed disks with scatter of contact interface characteristics[J].Journal of Turbomachinery,2005,127(1): 128-136.

    [109]PETROV E P,ZACHARIADIS Z I,BERETTA A,et al.A study of nonlinear vibrations in a frictionally damped turbine bladed disk with comprehensive modeling of aerodynamic effects[J].Journal of Engineering for Gas Turbines& Power,2013,135(3): 1239-1251.

    [110]PETROV E P,DEPARTMENT M E,FRICTION,et al.Analysis of flutter-induced limit cycle oscillations in gas-turbine structures with[J].Journal of Turbomachinery,2014,134:61011-61018.

    [111]SEVER I A,PETROV E E,EWINS D J.Experimental and numerical investigation of rotating bladed disk forced response using under-platform friction dampers[C]// Proceedings of GT2007,ASME Turbo Expo 2007: Power for Land,Sea and Air.2007.

    [112]SAITO A,CASTANIER M P,PIERRE C.Effects of a cracked blade on mistuned turbine engine rotor vibration[J].Journal of Vibration and Acoustics,2009,131(6): 49-62.

    [113]D’SOUZA K,SAITO A,EPUREANU B I.Reduced-order modeling for nonlinear analysis of cracked mistuned multistage bladed-disk systems[J].AIAA Journal,2012,50(2): 304-312.

    [114]JUNG C,SAITO A,EPUREANU B I.Detection of cracks in mistuned bladed disks using reduced-order models and vibration data[J].Journal of Vibration and Acoustics.2012,134(6): 1-10.

    [115]WANG S,ZI Y,LI B,et al.Reduced-order modeling for mistuned centrifugal impellers with crack damages[J].Journal of Sound and Vibration,2014,333(25): 6979-6995.

    [116]CAPIEZ-LERNOUT E,SOIZE C,MBAYE M.Geometric nonlinear dynamic analysis of uncertain structures with cyclic symmetry-Application to a mistuned industrial bladed disk[C]// International Conference on Uncertainty in Structural Dynamics.2014.

    [117]CAPIEZ-LERNOUT E E L,SOIZE C,MBAYE M.Mistuning analysis and uncertainty quantification of an industrial bladed disk with geometrical nonlinearity[J].Journal of Sound and Vibration,2015,356: 124-143.

    [118]CAPIEZ-LERNOUT E E L,SOIZE C,MBAYE M.Uncertainty quantification for an industrial mistuned bladed disk with geometrical nonlinearities[C]//ASME Turbo Expo 2015: Turbine Technical Conference and Exposition.2015.

    [119]CAPIEZ-LERNOUT E E L,SOIZE C,MBAYE M.Uncertainty quantification in nonlinear structural dynamics for mistuned bladed disks[C]//UNCECOMP 2015,1st ECCOMAS Thematic International Conference on Uncertainty Quantification in Computational Sciences and Engineering.2015.

    猜你喜歡
    王建軍葉盤模態(tài)
    Fe-doped ZnS film fabricated by electron beam evaporation and its application as saturable absorber for Er:ZBLAN fiber laser*
    憫蠶
    位卑未敢忘憂國——記上海馳愛機(jī)械有限公司總經(jīng)理王建軍
    王建軍接替宋麗萍任深交所總經(jīng)理
    國內(nèi)多模態(tài)教學(xué)研究回顧與展望
    某型航空發(fā)動(dòng)機(jī)整體葉盤強(qiáng)度分析*
    基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識(shí)別
    由單個(gè)模態(tài)構(gòu)造對稱簡支梁的抗彎剛度
    基于FE-ERSM航空發(fā)動(dòng)機(jī)葉盤結(jié)構(gòu)可靠性研究
    級(jí)間接觸耦合的失諧葉盤模態(tài)局部化問題研究
    飞空精品影院首页| 1024视频免费在线观看| 夜夜夜夜夜久久久久| 9色porny在线观看| 亚洲欧美日韩另类电影网站| 男女无遮挡免费网站观看| 岛国毛片在线播放| 午夜精品久久久久久毛片777| 欧美精品av麻豆av| 欧美精品亚洲一区二区| 精品亚洲成国产av| 免费观看av网站的网址| 精品欧美一区二区三区在线| 两性午夜刺激爽爽歪歪视频在线观看 | 国产深夜福利视频在线观看| 国产亚洲欧美在线一区二区| 麻豆国产av国片精品| av天堂久久9| 亚洲国产精品成人久久小说| 国产在线观看jvid| 一本一本久久a久久精品综合妖精| 电影成人av| 777米奇影视久久| 精品久久久精品久久久| 无遮挡黄片免费观看| 久久久水蜜桃国产精品网| 伊人亚洲综合成人网| 久久久久国产一级毛片高清牌| 精品高清国产在线一区| 97在线人人人人妻| 狠狠婷婷综合久久久久久88av| 欧美人与性动交α欧美精品济南到| 久久影院123| av免费在线观看网站| 美女午夜性视频免费| 成人亚洲精品一区在线观看| 又大又爽又粗| 欧美另类一区| 国产精品久久久久久精品电影小说| 少妇粗大呻吟视频| 成人亚洲精品一区在线观看| 女警被强在线播放| 啦啦啦在线免费观看视频4| 999久久久精品免费观看国产| 亚洲成人国产一区在线观看| 亚洲色图 男人天堂 中文字幕| 搡老乐熟女国产| 一二三四在线观看免费中文在| 久久精品熟女亚洲av麻豆精品| 欧美激情高清一区二区三区| 一本色道久久久久久精品综合| 欧美日韩国产mv在线观看视频| 国产精品影院久久| 欧美xxⅹ黑人| 免费不卡黄色视频| 五月天丁香电影| 国产野战对白在线观看| 狠狠狠狠99中文字幕| 中文字幕色久视频| 99国产精品一区二区三区| 无遮挡黄片免费观看| 国产亚洲精品久久久久5区| 成人国语在线视频| 欧美日韩av久久| av欧美777| 天天操日日干夜夜撸| 欧美日韩亚洲国产一区二区在线观看 | 99国产精品一区二区蜜桃av | av片东京热男人的天堂| 国产在视频线精品| 久久亚洲精品不卡| 亚洲五月婷婷丁香| 久久国产精品人妻蜜桃| 建设人人有责人人尽责人人享有的| 亚洲精品一二三| 啪啪无遮挡十八禁网站| 亚洲情色 制服丝袜| 久久久精品区二区三区| 国产欧美亚洲国产| 亚洲免费av在线视频| 久久精品国产综合久久久| 精品亚洲乱码少妇综合久久| 免费在线观看完整版高清| 亚洲五月婷婷丁香| 免费在线观看日本一区| 男女午夜视频在线观看| 精品福利永久在线观看| 久久国产精品人妻蜜桃| av线在线观看网站| 老司机影院毛片| 99热网站在线观看| 欧美乱码精品一区二区三区| 人妻 亚洲 视频| 久久免费观看电影| 人人澡人人妻人| 婷婷成人精品国产| 日韩欧美免费精品| 高清在线国产一区| 午夜视频精品福利| 美女脱内裤让男人舔精品视频| 极品人妻少妇av视频| 搡老熟女国产l中国老女人| 涩涩av久久男人的天堂| 精品福利永久在线观看| 久久热在线av| 老司机影院成人| 国产成人系列免费观看| 人人妻人人澡人人爽人人夜夜| 亚洲五月婷婷丁香| 中文字幕高清在线视频| 男人添女人高潮全过程视频| 青草久久国产| 中文字幕av电影在线播放| 欧美+亚洲+日韩+国产| 精品久久久久久久毛片微露脸 | 成人国产一区最新在线观看| 久久精品国产亚洲av高清一级| 午夜久久久在线观看| 久久国产精品影院| 国产精品.久久久| 日韩制服骚丝袜av| 不卡av一区二区三区| 两个人免费观看高清视频| 亚洲精品久久久久久婷婷小说| 精品一区二区三区av网在线观看 | 免费少妇av软件| 国产精品久久久久久人妻精品电影 | 制服人妻中文乱码| 国产无遮挡羞羞视频在线观看| 日韩,欧美,国产一区二区三区| 18禁黄网站禁片午夜丰满| 丰满少妇做爰视频| 久久毛片免费看一区二区三区| 亚洲五月婷婷丁香| netflix在线观看网站| 国产又爽黄色视频| 制服诱惑二区| 亚洲专区字幕在线| 午夜免费观看性视频| 黄片大片在线免费观看| 热re99久久国产66热| 亚洲成人免费av在线播放| 男女之事视频高清在线观看| 国产精品一区二区免费欧美 | 欧美一级毛片孕妇| 午夜两性在线视频| 19禁男女啪啪无遮挡网站| 夜夜夜夜夜久久久久| 黄色视频在线播放观看不卡| 美女高潮喷水抽搐中文字幕| 成人国语在线视频| 老汉色∧v一级毛片| 亚洲av日韩精品久久久久久密| 国产97色在线日韩免费| 黄网站色视频无遮挡免费观看| 国产精品.久久久| 可以免费在线观看a视频的电影网站| 狂野欧美激情性bbbbbb| 国产精品一区二区在线观看99| 可以免费在线观看a视频的电影网站| 久久久久国产精品人妻一区二区| 久久久国产精品麻豆| 中文字幕色久视频| 精品国产超薄肉色丝袜足j| 女人久久www免费人成看片| 国产97色在线日韩免费| 丝袜在线中文字幕| 永久免费av网站大全| 成人18禁高潮啪啪吃奶动态图| 亚洲国产成人一精品久久久| 纵有疾风起免费观看全集完整版| 黄色毛片三级朝国网站| 精品人妻在线不人妻| 欧美成狂野欧美在线观看| 真人做人爱边吃奶动态| 精品国产一区二区三区四区第35| 欧美日韩福利视频一区二区| 精品福利观看| 叶爱在线成人免费视频播放| 久久国产精品人妻蜜桃| 精品人妻1区二区| 母亲3免费完整高清在线观看| 91老司机精品| 在线观看免费视频网站a站| 一区在线观看完整版| 天天躁日日躁夜夜躁夜夜| 国产免费av片在线观看野外av| 久9热在线精品视频| 咕卡用的链子| av视频免费观看在线观看| 久久九九热精品免费| 一本久久精品| 亚洲精品国产精品久久久不卡| 少妇的丰满在线观看| 久久人人爽人人片av| videosex国产| 精品第一国产精品| 久久国产精品男人的天堂亚洲| 汤姆久久久久久久影院中文字幕| 国产1区2区3区精品| 亚洲少妇的诱惑av| 日韩人妻精品一区2区三区| 日本av手机在线免费观看| 国产亚洲av片在线观看秒播厂| 1024视频免费在线观看| 不卡av一区二区三区| 国产欧美日韩一区二区三 | 欧美日韩av久久| 99精品欧美一区二区三区四区| 久久久精品国产亚洲av高清涩受| 亚洲五月色婷婷综合| www.999成人在线观看| 99热网站在线观看| 丰满饥渴人妻一区二区三| 久久人妻福利社区极品人妻图片| 看免费av毛片| 国产不卡av网站在线观看| 国产成人av教育| 欧美日韩av久久| 欧美日韩中文字幕国产精品一区二区三区 | 别揉我奶头~嗯~啊~动态视频 | 一区二区三区精品91| 美女福利国产在线| 一区二区日韩欧美中文字幕| 国产日韩欧美亚洲二区| 日韩视频一区二区在线观看| 9色porny在线观看| 男女无遮挡免费网站观看| 女警被强在线播放| 久久久久久免费高清国产稀缺| 久久人妻熟女aⅴ| 欧美久久黑人一区二区| 青春草视频在线免费观看| 欧美人与性动交α欧美精品济南到| 国产精品熟女久久久久浪| 亚洲第一青青草原| 中文精品一卡2卡3卡4更新| 亚洲精品中文字幕在线视频| 亚洲一码二码三码区别大吗| 亚洲国产欧美网| 黄片小视频在线播放| 极品人妻少妇av视频| 亚洲五月色婷婷综合| 亚洲性夜色夜夜综合| 法律面前人人平等表现在哪些方面 | 69av精品久久久久久 | 男人舔女人的私密视频| 国产一区二区激情短视频 | 人人妻人人爽人人添夜夜欢视频| 欧美另类亚洲清纯唯美| 99国产极品粉嫩在线观看| 亚洲五月婷婷丁香| 97精品久久久久久久久久精品| 国产一区二区三区在线臀色熟女 | 久久精品人人爽人人爽视色| 久久国产精品人妻蜜桃| 午夜影院在线不卡| 免费不卡黄色视频| 韩国精品一区二区三区| 后天国语完整版免费观看| 在线十欧美十亚洲十日本专区| 五月开心婷婷网| 高清在线国产一区| 久久亚洲精品不卡| 人妻久久中文字幕网| 18禁观看日本| 日本av免费视频播放| 热99国产精品久久久久久7| 十八禁网站网址无遮挡| 久久热在线av| 午夜福利一区二区在线看| 国产精品亚洲av一区麻豆| 欧美亚洲日本最大视频资源| 在线观看免费日韩欧美大片| 久久久久网色| 免费观看人在逋| 精品国内亚洲2022精品成人 | 我要看黄色一级片免费的| 人妻一区二区av| av又黄又爽大尺度在线免费看| 国产成人欧美| 又大又爽又粗| 国产精品99久久99久久久不卡| 国产男女超爽视频在线观看| 亚洲情色 制服丝袜| 午夜福利乱码中文字幕| 中文字幕色久视频| 无遮挡黄片免费观看| 女警被强在线播放| 国产无遮挡羞羞视频在线观看| 满18在线观看网站| netflix在线观看网站| 亚洲国产日韩一区二区| 2018国产大陆天天弄谢| 在线观看一区二区三区激情| 99久久国产精品久久久| 999久久久国产精品视频| 亚洲激情五月婷婷啪啪| 天天躁日日躁夜夜躁夜夜| 欧美激情极品国产一区二区三区| 啦啦啦中文免费视频观看日本| 国产91精品成人一区二区三区 | 狂野欧美激情性bbbbbb| 老鸭窝网址在线观看| 女人爽到高潮嗷嗷叫在线视频| 午夜福利,免费看| 成人国产av品久久久| 午夜福利乱码中文字幕| 777米奇影视久久| 日本五十路高清| 国产熟女午夜一区二区三区| 黄网站色视频无遮挡免费观看| 国产精品 国内视频| 国产精品熟女久久久久浪| 黄色视频,在线免费观看| 久久 成人 亚洲| 91精品三级在线观看| 欧美少妇被猛烈插入视频| 国产亚洲精品第一综合不卡| 性少妇av在线| 久久精品国产a三级三级三级| 如日韩欧美国产精品一区二区三区| 午夜激情av网站| 最近最新中文字幕大全免费视频| 老司机影院成人| 中文字幕另类日韩欧美亚洲嫩草| 高清欧美精品videossex| 久久久国产精品麻豆| 国产日韩一区二区三区精品不卡| 精品国产一区二区三区四区第35| 亚洲精品国产一区二区精华液| 精品亚洲成国产av| 少妇的丰满在线观看| 极品少妇高潮喷水抽搐| 成人国产av品久久久| kizo精华| 成人影院久久| 亚洲av美国av| 夫妻午夜视频| 亚洲精品久久成人aⅴ小说| 亚洲欧美成人综合另类久久久| 国产高清视频在线播放一区 | 99久久人妻综合| 18禁黄网站禁片午夜丰满| 久久久久久人人人人人| 黄色片一级片一级黄色片| 欧美亚洲日本最大视频资源| 99精品久久久久人妻精品| 亚洲精品在线美女| 后天国语完整版免费观看| 妹子高潮喷水视频| 午夜精品久久久久久毛片777| 亚洲精品中文字幕在线视频| 别揉我奶头~嗯~啊~动态视频 | 国产精品亚洲av一区麻豆| videos熟女内射| 菩萨蛮人人尽说江南好唐韦庄| 在线亚洲精品国产二区图片欧美| 久久九九热精品免费| 丰满饥渴人妻一区二区三| 欧美亚洲 丝袜 人妻 在线| 99热全是精品| tocl精华| 欧美 日韩 精品 国产| 69精品国产乱码久久久| 成人国产一区最新在线观看| 免费不卡黄色视频| 搡老岳熟女国产| 各种免费的搞黄视频| 男女午夜视频在线观看| 亚洲国产欧美在线一区| 亚洲国产看品久久| av一本久久久久| 久久亚洲精品不卡| 久久午夜综合久久蜜桃| 精品少妇黑人巨大在线播放| 精品国内亚洲2022精品成人 | 老司机靠b影院| 欧美成狂野欧美在线观看| 亚洲精品国产区一区二| 99香蕉大伊视频| 悠悠久久av| 在线观看免费午夜福利视频| 亚洲精品一区蜜桃| 亚洲第一欧美日韩一区二区三区 | 久久久欧美国产精品| 国产精品欧美亚洲77777| 999久久久精品免费观看国产| 亚洲美女黄色视频免费看| 亚洲精品在线美女| 国产欧美日韩综合在线一区二区| 国产成人欧美| 18禁黄网站禁片午夜丰满| 国产深夜福利视频在线观看| 首页视频小说图片口味搜索| 久久精品熟女亚洲av麻豆精品| 免费久久久久久久精品成人欧美视频| 国产欧美日韩一区二区三 | 91成人精品电影| 亚洲国产av影院在线观看| 日韩欧美国产一区二区入口| 免费不卡黄色视频| 国产真人三级小视频在线观看| 欧美另类亚洲清纯唯美| 18禁国产床啪视频网站| 一级毛片电影观看| 91精品三级在线观看| 亚洲精品av麻豆狂野| 国产三级黄色录像| 精品国产国语对白av| 91av网站免费观看| 97人妻天天添夜夜摸| 最近最新免费中文字幕在线| 亚洲中文字幕日韩| 69精品国产乱码久久久| 日本五十路高清| 少妇的丰满在线观看| 日韩 欧美 亚洲 中文字幕| 日韩中文字幕视频在线看片| 欧美在线黄色| 国产精品一区二区在线不卡| 日韩三级视频一区二区三区| 丝袜美足系列| 国产伦人伦偷精品视频| 一二三四在线观看免费中文在| 一区二区三区精品91| 亚洲国产欧美日韩在线播放| 亚洲av男天堂| 欧美黑人精品巨大| 国产精品一区二区在线不卡| 大香蕉久久网| 国产精品久久久久成人av| 交换朋友夫妻互换小说| 国产视频一区二区在线看| 国产一区二区三区av在线| 女人被躁到高潮嗷嗷叫费观| 精品久久久久久电影网| 青春草视频在线免费观看| 亚洲精品国产色婷婷电影| 99re6热这里在线精品视频| 日本av手机在线免费观看| 色94色欧美一区二区| 天天躁日日躁夜夜躁夜夜| 国产精品av久久久久免费| 亚洲自偷自拍图片 自拍| 国产福利在线免费观看视频| 亚洲午夜精品一区,二区,三区| 一级片免费观看大全| 欧美成人午夜精品| 亚洲国产精品成人久久小说| 亚洲精品日韩在线中文字幕| 老熟妇仑乱视频hdxx| 亚洲免费av在线视频| 精品一品国产午夜福利视频| 青春草亚洲视频在线观看| 国产欧美亚洲国产| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜日韩欧美国产| 免费女性裸体啪啪无遮挡网站| 亚洲va日本ⅴa欧美va伊人久久 | 在线十欧美十亚洲十日本专区| 1024视频免费在线观看| 三上悠亚av全集在线观看| 亚洲五月婷婷丁香| 国产在线一区二区三区精| 亚洲免费av在线视频| 精品人妻在线不人妻| 亚洲欧美日韩高清在线视频 | 首页视频小说图片口味搜索| 久久毛片免费看一区二区三区| 高潮久久久久久久久久久不卡| 19禁男女啪啪无遮挡网站| 午夜福利影视在线免费观看| 免费看十八禁软件| 伊人久久大香线蕉亚洲五| 99香蕉大伊视频| 免费黄频网站在线观看国产| 2018国产大陆天天弄谢| 国产精品二区激情视频| 一本久久精品| 十八禁人妻一区二区| 亚洲国产成人一精品久久久| 亚洲国产欧美在线一区| 深夜精品福利| 9191精品国产免费久久| 日韩人妻精品一区2区三区| 悠悠久久av| 国产精品亚洲av一区麻豆| 他把我摸到了高潮在线观看 | 亚洲精品国产av蜜桃| 狂野欧美激情性xxxx| 少妇被粗大的猛进出69影院| 人人妻人人添人人爽欧美一区卜| 午夜福利免费观看在线| 国产精品一区二区精品视频观看| 麻豆国产av国片精品| tube8黄色片| 成人18禁高潮啪啪吃奶动态图| 欧美日韩中文字幕国产精品一区二区三区 | 91精品国产国语对白视频| 亚洲av日韩在线播放| 热re99久久精品国产66热6| 国产精品秋霞免费鲁丝片| 国产又爽黄色视频| 亚洲欧美日韩高清在线视频 | 99热国产这里只有精品6| 一个人免费看片子| 9色porny在线观看| 国产有黄有色有爽视频| 亚洲熟女毛片儿| 欧美日韩福利视频一区二区| 巨乳人妻的诱惑在线观看| 人人澡人人妻人| 国产有黄有色有爽视频| 国产精品一二三区在线看| 性少妇av在线| 亚洲精品一二三| netflix在线观看网站| 色精品久久人妻99蜜桃| 纯流量卡能插随身wifi吗| 精品久久蜜臀av无| 十八禁网站免费在线| 老司机深夜福利视频在线观看 | av在线app专区| 多毛熟女@视频| 亚洲精品美女久久av网站| 色婷婷久久久亚洲欧美| 欧美成狂野欧美在线观看| 丝袜喷水一区| 丰满迷人的少妇在线观看| 国产精品av久久久久免费| 一级黄色大片毛片| 99久久国产精品久久久| 日韩电影二区| av国产精品久久久久影院| 国产成人精品在线电影| 热99国产精品久久久久久7| tocl精华| 美女中出高潮动态图| 成人影院久久| 国产在线视频一区二区| 色播在线永久视频| 亚洲av国产av综合av卡| 国产精品影院久久| 制服人妻中文乱码| 少妇裸体淫交视频免费看高清 | 亚洲激情五月婷婷啪啪| 国产在线一区二区三区精| 成人影院久久| 亚洲九九香蕉| 亚洲欧美激情在线| 桃花免费在线播放| 99精品欧美一区二区三区四区| av不卡在线播放| 69精品国产乱码久久久| av欧美777| 国产成人精品久久二区二区免费| 亚洲精品自拍成人| 精品亚洲成a人片在线观看| 国产精品一区二区免费欧美 | 久久久久久免费高清国产稀缺| 各种免费的搞黄视频| 日本黄色日本黄色录像| www日本在线高清视频| 一区在线观看完整版| 成年动漫av网址| 中文字幕av电影在线播放| 老鸭窝网址在线观看| 乱人伦中国视频| 精品第一国产精品| 亚洲七黄色美女视频| 三上悠亚av全集在线观看| 欧美少妇被猛烈插入视频| 成在线人永久免费视频| 另类精品久久| 91麻豆av在线| 久久精品人人爽人人爽视色| 日韩一区二区三区影片| 十分钟在线观看高清视频www| www.999成人在线观看| 曰老女人黄片| 精品国产一区二区久久| 国产又色又爽无遮挡免| 国产精品久久久久成人av| 最黄视频免费看| 一区福利在线观看| 午夜福利视频精品| 久久精品成人免费网站| 国产成人av激情在线播放| 国产福利在线免费观看视频| 嫩草影视91久久| 精品第一国产精品| 永久免费av网站大全| 乱人伦中国视频| 青青草视频在线视频观看| 交换朋友夫妻互换小说| 伊人久久大香线蕉亚洲五| 精品一区在线观看国产| 国产亚洲精品久久久久5区| 老熟妇仑乱视频hdxx| 欧美精品av麻豆av| 日本a在线网址| 伊人久久大香线蕉亚洲五| 欧美久久黑人一区二区| 69av精品久久久久久 | 国产成人精品在线电影| 性色av一级| 午夜久久久在线观看| 国内毛片毛片毛片毛片毛片| 国产视频一区二区在线看| 爱豆传媒免费全集在线观看| 国产三级黄色录像| 日韩中文字幕欧美一区二区| 大码成人一级视频| 亚洲av电影在线进入|