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

    迷宮密封-滾動(dòng)軸承-懸臂轉(zhuǎn)子系統(tǒng)非線性動(dòng)力學(xué)特性分析

    2020-05-13 14:31羅躍綱王鵬飛王晨勇徐昊
    振動(dòng)工程學(xué)報(bào) 2020年2期
    關(guān)鍵詞:有限元

    羅躍綱 王鵬飛 王晨勇 徐昊

    摘要:對(duì)于帶有迷宮密封的航空發(fā)動(dòng)機(jī)轉(zhuǎn)子系統(tǒng)氣流激振問(wèn)題,基于有限元理論,應(yīng)用非線性滾動(dòng)軸承支承力模型以及Muzynska密封力模型建立了兩個(gè)滾動(dòng)軸承支承的迷宮密封一懸臂轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)模型,并運(yùn)用Newmark-β數(shù)值積分法求解得到系統(tǒng)在不同轉(zhuǎn)速、偏心量和密封結(jié)構(gòu)參數(shù)下的動(dòng)力學(xué)響應(yīng)特征。研究結(jié)果表明,系統(tǒng)在一定轉(zhuǎn)速范圍內(nèi)作周期一運(yùn)動(dòng),隨著轉(zhuǎn)速的升高系統(tǒng)發(fā)生失穩(wěn)并作擬周期運(yùn)動(dòng);適當(dāng)增大偏心量會(huì)導(dǎo)致轉(zhuǎn)子在共振區(qū)出現(xiàn)偏心力所引起的短暫的混沌運(yùn)動(dòng);增大密封間隙會(huì)使系統(tǒng)在高轉(zhuǎn)速區(qū)重新回歸周期一運(yùn)動(dòng),而且失穩(wěn)區(qū)域也隨之減小;適當(dāng)提高密封長(zhǎng)度,系統(tǒng)僅表現(xiàn)為周期一運(yùn)動(dòng),但繼續(xù)增大密封長(zhǎng)度,懸臂端承受密封圓盤的重量也將提高,失穩(wěn)轉(zhuǎn)速提前;另外還分析了失穩(wěn)轉(zhuǎn)速和密封力的影響因素及其影響規(guī)律,為轉(zhuǎn)子系統(tǒng)的密封激振故障診斷及密封結(jié)構(gòu)優(yōu)化設(shè)計(jì)提供一定的理論依據(jù)。

    關(guān)鍵詞:非線性振動(dòng);懸臂轉(zhuǎn)子系統(tǒng);迷宮密封;密封力;有限元

    中圖分類號(hào):0322;0347.6文獻(xiàn)標(biāo)志碼:A 文章編號(hào):1004-4523(2020)02-0256-09

    DOI:10.16385/j.cnki.issn.1004-4523.2020.02.005

    引言

    迷宮密封是普遍安裝在現(xiàn)代航空發(fā)動(dòng)機(jī)、汽輪機(jī)等旋轉(zhuǎn)機(jī)械結(jié)構(gòu)中的有效封嚴(yán)結(jié)構(gòu),它作為一種非接觸式密封,具有結(jié)構(gòu)簡(jiǎn)單、耗能小、使用壽命長(zhǎng)、無(wú)需潤(rùn)滑等特點(diǎn),其作用是減少軸端與各級(jí)問(wèn)的流體泄漏損失。對(duì)于帶有迷宮密封的轉(zhuǎn)子系統(tǒng),由于工作轉(zhuǎn)速的提高、轉(zhuǎn)子柔性增大和高參數(shù)密封致使密封激振作用極易發(fā)生,并導(dǎo)致轉(zhuǎn)子失穩(wěn)。因此,為加強(qiáng)該類系統(tǒng)的運(yùn)行穩(wěn)定性與工作安全性,研究含有密封激振力作用下的轉(zhuǎn)子系統(tǒng)動(dòng)力學(xué)特征并分析一些典型參數(shù)影響規(guī)律有著重要的意義。

    多年以來(lái),國(guó)內(nèi)外許多專家學(xué)者在含有密封的轉(zhuǎn)子動(dòng)力學(xué)領(lǐng)域作了大量研究,比如在求解密封動(dòng)力特性系數(shù)并分析其影響因素方面,wang等通過(guò)應(yīng)用單控制體模型及攝動(dòng)法對(duì)含有迷宮密封的轉(zhuǎn)子系統(tǒng)進(jìn)行動(dòng)力學(xué)建模并對(duì)其進(jìn)行計(jì)算;文獻(xiàn)[2-3]利用cFX-TAscflow流體動(dòng)力學(xué)軟件計(jì)算了密封轉(zhuǎn)子動(dòng)力系數(shù),并研究了它的影響因素等。通過(guò)這些研究進(jìn)而得到能夠反映密封性能的泄漏量以及穩(wěn)定性的影響因素,并提出了減少密封泄漏、抑制密封激振的改進(jìn)措施,例如sun等、冀大偉等、陳堯興等還分別考慮了安裝渦流制動(dòng)器、防旋板及非均勻進(jìn)氣等情況的影響。此外Childs通過(guò)實(shí)驗(yàn)對(duì)比發(fā)現(xiàn)在定子上安裝密封齒要比安裝于轉(zhuǎn)子的穩(wěn)定性更好;wang等還發(fā)現(xiàn)位于密封腔的流體泄漏對(duì)于系統(tǒng)失穩(wěn)轉(zhuǎn)速也存在很大影響。

    目前的文獻(xiàn)中多以汽輪機(jī)轉(zhuǎn)子為研究對(duì)象,或?qū)⒛P秃?jiǎn)化為Jeffcott轉(zhuǎn)子,并以滑動(dòng)軸承作為支承。而對(duì)于以航空發(fā)動(dòng)機(jī)轉(zhuǎn)子為代表的帶有滾動(dòng)軸承支承的懸臂轉(zhuǎn)子密封結(jié)構(gòu)的動(dòng)力學(xué)特性還研究較少。本文以某型渦扇航空發(fā)動(dòng)機(jī)轉(zhuǎn)子為研究對(duì)象,為突出主要矛盾將葉片簡(jiǎn)化為圓盤,考慮了陀螺效應(yīng)的影響,采用Muszynska密封力模型和滾動(dòng)軸承支承力模型建立密封一滾動(dòng)軸承一懸臂轉(zhuǎn)子系統(tǒng)有限元模型,采用Newmark-β數(shù)值算法對(duì)系統(tǒng)進(jìn)行模擬仿真分析,研究不同轉(zhuǎn)速、偏心量和密封結(jié)構(gòu)參數(shù)等對(duì)系統(tǒng)的影響,并得到了失穩(wěn)轉(zhuǎn)速和密封力的影響因素及其影響規(guī)律,為密封一轉(zhuǎn)子系統(tǒng)的氣流激振故障診斷及結(jié)構(gòu)優(yōu)化設(shè)計(jì)提供一定的理論依據(jù)。

    1迷宮密封一懸臂轉(zhuǎn)子一滾動(dòng)軸承系統(tǒng)模型

    1.1轉(zhuǎn)子系統(tǒng)的動(dòng)力學(xué)模型

    1.2密封力模型

    本文采用的是Muszynska密封力模型,如圖4所示為某一密封腔的徑向截面剖視圖,其中靠近轉(zhuǎn)子處的流體周向角速度為ω,靠近定子處的流體周向角速度則降為0,引入τ來(lái)表示流體周向平均流速比,則密封腔中平均流速可用τω來(lái)表示。其中密封力對(duì)轉(zhuǎn)子的擾動(dòng)反力隨流體一起以平均流速τω繞轉(zhuǎn)子旋轉(zhuǎn),因此密封力的旋轉(zhuǎn)效應(yīng)是致使系統(tǒng)失穩(wěn)的主因。

    2.1轉(zhuǎn)速的影響

    選定的密封結(jié)構(gòu)尺寸為:密封問(wèn)隙c=0.06mm,壓差△P=0.2MPa,軸向流速Va=30m/s;密封圓盤半徑R=100mm,圓盤長(zhǎng)度l=25mm,偏心量me=500g·mm;圖6和7分別為無(wú)/有密封力作用下系統(tǒng)轉(zhuǎn)速n從2000r/min升高到20000r/min過(guò)程中的分岔圖,可見(jiàn)無(wú)密封力時(shí)整個(gè)升速過(guò)程系統(tǒng)做穩(wěn)定的周期一運(yùn)動(dòng),有密封力后轉(zhuǎn)速在9800r/min時(shí)發(fā)生失穩(wěn);圖8為有密封力作用時(shí),轉(zhuǎn)速在5000r/min及16000r/min情況下的軸心運(yùn)動(dòng)軌跡圖、時(shí)域波形圖、頻率響應(yīng)圖和Poincar6圖。由圖8(a)可以看出,在轉(zhuǎn)速較低時(shí),密封圓盤處的軸心軌跡呈現(xiàn)為橢圓形,時(shí)域波形運(yùn)行平穩(wěn),頻譜中存在工頻以及一定的2倍頻分量,Poincar6圖中僅存在孤立的單點(diǎn),這些特征表明當(dāng)前轉(zhuǎn)子運(yùn)行平穩(wěn),密封力的作用并不明顯,系統(tǒng)作周期一運(yùn)動(dòng);當(dāng)轉(zhuǎn)速上升至16000r/min時(shí)的響應(yīng)如圖8(b)所示,軸心運(yùn)動(dòng)軌跡表現(xiàn)出含有多個(gè)圓環(huán)相互嵌套成橢圓的特征,頻率譜中呈現(xiàn)出工頻和幅值較高的1.1倍頻率,Poincar6圖表現(xiàn)為若干個(gè)離散的點(diǎn)所組成的封閉圓環(huán),這些現(xiàn)象說(shuō)明在系統(tǒng)失穩(wěn)以后由于氣流激振作用轉(zhuǎn)子作擬周期運(yùn)動(dòng)。

    2.2偏心量的影響

    圖9是偏心量1500g·mm時(shí)轉(zhuǎn)子升速過(guò)程的分岔圖,與圖6對(duì)比可以發(fā)現(xiàn)整個(gè)升速過(guò)程中有兩個(gè)失穩(wěn)區(qū),其中一個(gè)發(fā)生于轉(zhuǎn)速較低的3000r/min,是由于臨近共振區(qū)偏心力的作用增強(qiáng)而引起的,另一個(gè)發(fā)生在11200r/min,這是由于氣流激振力所引起的密封失穩(wěn)。圖10表示轉(zhuǎn)速穩(wěn)定在3200r/min時(shí)的系統(tǒng)響應(yīng),此時(shí)的軸心運(yùn)動(dòng)軌跡曲線比較紊亂,頻譜圖中在4倍頻內(nèi)出現(xiàn)了連續(xù)的分頻譜成分,Poincare映射圖上呈現(xiàn)出了很多散亂的點(diǎn),結(jié)合圖8可以發(fā)現(xiàn)轉(zhuǎn)子偏心量適當(dāng)增大后在3000-4600r/min的共振區(qū)附近出現(xiàn)了偏心力所引起的短暫的混沌運(yùn)動(dòng)。雖然在整體上均先后經(jīng)歷了周期一運(yùn)動(dòng)、擬周期運(yùn)動(dòng)、周期一運(yùn)動(dòng)的轉(zhuǎn)變,但發(fā)生分岔失穩(wěn)的轉(zhuǎn)速

    2.3密封間隙的影響

    圖11(a)為密封問(wèn)隙c=0.14mm時(shí)豎直向轉(zhuǎn)速分岔圖,與圖6對(duì)比發(fā)現(xiàn),密封問(wèn)隙量由0.06mm提高至0.14mm,系統(tǒng)發(fā)生了豐富的動(dòng)力學(xué)行為,在2000-11400r/rain時(shí)作周期一運(yùn)動(dòng),在11400-17600r/rain時(shí)作擬周期運(yùn)動(dòng),而當(dāng)轉(zhuǎn)子轉(zhuǎn)速超過(guò)17600r/min時(shí)轉(zhuǎn)子系統(tǒng)發(fā)生跳躍并重新回到周期一運(yùn)動(dòng),而且此時(shí)y方向的振動(dòng)位移量明顯減小;通過(guò)圖11不難看出密封問(wèn)隙逐漸增加后,系統(tǒng)向后推延,而重新回歸穩(wěn)態(tài)運(yùn)行的轉(zhuǎn)速亦有所提前,系統(tǒng)的失穩(wěn)區(qū)域逐漸減小。另外在轉(zhuǎn)速較低時(shí)密封問(wèn)隙對(duì)系統(tǒng)動(dòng)力學(xué)行為的影響并不明顯,隨著轉(zhuǎn)速的不斷升高密封問(wèn)隙的影響作用不斷增強(qiáng),這是由于密封問(wèn)隙增大導(dǎo)致密封腔內(nèi)部流體泄漏量增大,而轉(zhuǎn)速的升高進(jìn)一步帶動(dòng)了氣體的流動(dòng),進(jìn)而加劇密封處流體泄漏量,致使密封激振力的影響迅速減弱。

    2.4密封長(zhǎng)度的影響

    以迷宮密封的齒腔寬度不變?yōu)榍疤?,通過(guò)增加密封腔的個(gè)數(shù)來(lái)增大有效密封長(zhǎng)度。圖12(a)和(b)分別表示在密封長(zhǎng)度50mm和100mm兩種典型情形下的系統(tǒng)分岔圖。對(duì)比圖6與12(a),密封長(zhǎng)度由25mm伸長(zhǎng)到50mm,系統(tǒng)于工作轉(zhuǎn)速范圍內(nèi)將不再呈現(xiàn)失穩(wěn)分岔以及擬周期運(yùn)動(dòng),僅表現(xiàn)為穩(wěn)定的同步運(yùn)動(dòng),而且系統(tǒng)的振動(dòng)位移有所減小,這與張恩杰等的結(jié)果是一致的;但當(dāng)密封長(zhǎng)度增大到100mm時(shí),轉(zhuǎn)子的失穩(wěn)轉(zhuǎn)速又將提前,這表明適當(dāng)增大密封長(zhǎng)度利于系統(tǒng)平穩(wěn)運(yùn)行,但過(guò)大的密封長(zhǎng)度會(huì)導(dǎo)致密封圓盤過(guò)于龐大,懸臂轉(zhuǎn)子的伸出端將承受較大的密封圓盤重量,反而不利于轉(zhuǎn)子運(yùn)行的穩(wěn)定,甚至?xí)?duì)支承處的滾動(dòng)軸承造成一定的損害。

    2.5系統(tǒng)失穩(wěn)轉(zhuǎn)速的影響因素分析

    失穩(wěn)轉(zhuǎn)速是反映密封一轉(zhuǎn)子系統(tǒng)穩(wěn)定性能的重要指標(biāo),圖13分別分析了軸向流速、密封壓力差、轉(zhuǎn)子偏心量、密封問(wèn)隙、密封半徑和長(zhǎng)度這6種因素對(duì)于失穩(wěn)轉(zhuǎn)速的關(guān)系曲線,其中圖13(a)為軸向流速在10-30m/s情況下與失穩(wěn)轉(zhuǎn)速的關(guān)系,不難發(fā)現(xiàn)軸向流速對(duì)失穩(wěn)轉(zhuǎn)速的變化表現(xiàn)為非線性關(guān)系,隨著軸向流速的不斷增加,失穩(wěn)轉(zhuǎn)速不斷減小,而且其變化率逐漸降低,即系統(tǒng)發(fā)生失穩(wěn)的情況會(huì)提前到來(lái),不利于系統(tǒng)穩(wěn)定;如圖13(b)所示,進(jìn)出口密封壓力差從0.1MPa增大到0.9MPa時(shí)對(duì)失穩(wěn)轉(zhuǎn)速的影響近似表現(xiàn)為線性關(guān)系,伴隨壓差逐漸升高,失穩(wěn)臨界轉(zhuǎn)速亦隨之向后推遲,這有益于密封轉(zhuǎn)子的穩(wěn)定運(yùn)行;圖13(c)和(d)分別表示轉(zhuǎn)子偏心量由500g·mm增至2500g·mm以及密封問(wèn)隙在0.06-0.22mm范圍內(nèi)的失穩(wěn)轉(zhuǎn)速影響關(guān)系圖,可以發(fā)現(xiàn)伴隨著二者的提高,失穩(wěn)轉(zhuǎn)速有升高的趨向,這表明適當(dāng)增大偏心和問(wèn)隙可提升系統(tǒng)穩(wěn)定性,但過(guò)大偏心極易致使轉(zhuǎn)定子問(wèn)產(chǎn)生碰摩故障,而且擴(kuò)大問(wèn)隙也會(huì)導(dǎo)致流體泄漏增多、密封性能下降,所以工程實(shí)際結(jié)構(gòu)設(shè)計(jì)中需要充分考慮其他因素的共同影響情況;圖13(e)和(f)為密封半徑以及長(zhǎng)度對(duì)于失穩(wěn)轉(zhuǎn)速的關(guān)系曲線,可見(jiàn)密封圓盤半徑對(duì)失穩(wěn)轉(zhuǎn)速的影響近似為線性的,增加密封半徑會(huì)導(dǎo)致系統(tǒng)失穩(wěn)轉(zhuǎn)速降低,系統(tǒng)穩(wěn)定性變差,而且改變密封半徑還會(huì)影響轉(zhuǎn)子系統(tǒng)其他結(jié)構(gòu)尺寸,提高轉(zhuǎn)子設(shè)計(jì)工作量,因此宜將密封件安裝于小半徑轉(zhuǎn)子處;當(dāng)有效密封長(zhǎng)度由10mm伸長(zhǎng)至100mm時(shí),其對(duì)失穩(wěn)轉(zhuǎn)速的影響呈現(xiàn)出先適度減小后急劇升高到最高點(diǎn)再下降的趨勢(shì),密封長(zhǎng)度取25mm時(shí)對(duì)應(yīng)的轉(zhuǎn)速達(dá)到最低點(diǎn),而在50-70mm時(shí)又迅速升高,這說(shuō)明增加有效密封長(zhǎng)度可以使密封腔內(nèi)平均流速降低,推遲氣流激振所引起的轉(zhuǎn)速失穩(wěn)閾值,有利于該系統(tǒng)在正常轉(zhuǎn)速區(qū)域內(nèi)平穩(wěn)運(yùn)轉(zhuǎn);但繼續(xù)延長(zhǎng)密封長(zhǎng)度,懸臂轉(zhuǎn)子所承受的密封圓盤的重量也會(huì)加大,反而不利于系統(tǒng)穩(wěn)定運(yùn)行。

    2.6密封力的影響因素分析

    圖14是工作轉(zhuǎn)速分別選取2500,4000和8000r/min三種情況下軸向流速、密封壓力差、密封問(wèn)隙、轉(zhuǎn)子偏心量以及密封半徑等因素對(duì)水平方向和豎直方向密封力的影響規(guī)律總結(jié),從整體上看轉(zhuǎn)速對(duì)密封力的影響是不斷增強(qiáng)的。由圖14(a)可以看出,隨著軸向流速由5m/s增大至30m/s,在轉(zhuǎn)速較低時(shí)水平方向的密封力逐漸減小,轉(zhuǎn)速較高時(shí)水平方向密封力負(fù)向減小,而豎直方向密封力整體呈負(fù)向減小趨勢(shì);由圖14(b)可看出,提高密封壓力差,水平方向的密封分力負(fù)向增加,且轉(zhuǎn)速越高其增幅越明顯,而在高轉(zhuǎn)速下密封壓差的增大使水平向密封力增速逐漸減緩,而其對(duì)豎直向密封分力的影響則與當(dāng)前轉(zhuǎn)子運(yùn)行的轉(zhuǎn)速有關(guān):在2500r/min時(shí)豎直向密封分力隨壓差升高而增大,而在4000r/rain時(shí)壓差對(duì)它并無(wú)太大的影響,轉(zhuǎn)速升高至8000r/rain其隨壓差的升高表現(xiàn)為先負(fù)向增加后負(fù)向降低;增大密封問(wèn)隙尺寸對(duì)密封激振力的影響如圖14(c)所示,當(dāng)轉(zhuǎn)速較低時(shí)水平方向分力負(fù)向增長(zhǎng),豎直方向的分力則為正向增長(zhǎng),但是當(dāng)轉(zhuǎn)速較高時(shí)二者均改變?yōu)樨?fù)向減小,而且隨著問(wèn)隙的不斷增大,密封力的變化率逐漸降低;圖14(d)為轉(zhuǎn)子不平衡量對(duì)密封力的影響關(guān)系曲線,隨著偏心量增大,水平方向和豎直方向的密封力均呈現(xiàn)出負(fù)向增長(zhǎng)的趨勢(shì);在圖14(e)中,當(dāng)密封半徑由50mm提高至100mm時(shí),2500r/rain時(shí)其對(duì)水平向密封分力的作用不太明顯,然而在轉(zhuǎn)速升高到4000r/rain時(shí)呈現(xiàn)出先負(fù)向增長(zhǎng)后又負(fù)向降低的變化規(guī)律,當(dāng)轉(zhuǎn)速高達(dá)8000r/rain時(shí)水平與豎直向的密封分力均負(fù)向增長(zhǎng),而低轉(zhuǎn)速區(qū)的豎直向密封分力則有正向升高的趨勢(shì)。

    3結(jié)論

    (1)通過(guò)對(duì)有/無(wú)密封力作用下系統(tǒng)動(dòng)力學(xué)特性對(duì)比發(fā)現(xiàn):密封力在轉(zhuǎn)速較低時(shí)的作用并不明顯,此時(shí)轉(zhuǎn)子作穩(wěn)定的周期一運(yùn)動(dòng);隨著轉(zhuǎn)速升高,當(dāng)超過(guò)一定閾值后系統(tǒng)失穩(wěn)并作擬周期運(yùn)動(dòng)。

    (2)適當(dāng)增大偏心量會(huì)導(dǎo)致轉(zhuǎn)子在共振區(qū)出現(xiàn)短暫的偏心力引起的混沌運(yùn)動(dòng);增加密封問(wèn)隙會(huì)使位于高轉(zhuǎn)速區(qū)的系統(tǒng)重新回歸周期一運(yùn)動(dòng),且系統(tǒng)失穩(wěn)區(qū)域隨著問(wèn)隙的增加而不斷收縮,這是高轉(zhuǎn)速導(dǎo)致密封氣流泄漏量增加,密封激振力的影響程度減弱而引起的。

    (3)適當(dāng)提高密封有效長(zhǎng)度,系統(tǒng)在工作轉(zhuǎn)速區(qū)將不再出現(xiàn)擬周期運(yùn)動(dòng),僅表現(xiàn)為單一的周期一運(yùn)動(dòng),但繼續(xù)增大密封長(zhǎng)度,懸臂端承受密封圓盤的重量也將升高,系統(tǒng)失穩(wěn)轉(zhuǎn)速提前,不利于系統(tǒng)的穩(wěn)定運(yùn)行。

    (4)密封長(zhǎng)度對(duì)失穩(wěn)轉(zhuǎn)速的影響最大,而密封問(wèn)隙對(duì)其的影響相對(duì)較小。增大密封壓差、偏心量、密封問(wèn)隙會(huì)使失穩(wěn)轉(zhuǎn)速升高,有利于系統(tǒng)的平穩(wěn)運(yùn)行;提高密封半徑和軸向流速將不利于轉(zhuǎn)子工作穩(wěn)定;提升轉(zhuǎn)速和偏心量,密封力亦隨之增加;提高軸向流速會(huì)使密封力減小,提高密封壓差易導(dǎo)致水平方向密封力負(fù)向增大;密封問(wèn)隙和密封半徑對(duì)密封力的影響比較復(fù)雜,均與當(dāng)前轉(zhuǎn)速有關(guān)。

    猜你喜歡
    有限元
    有限元基礎(chǔ)與應(yīng)用課程專業(yè)賦能改革與實(shí)踐
    基于有限元的Q345E鋼補(bǔ)焊焊接殘余應(yīng)力的數(shù)值模擬
    定方位射孔器燃爆過(guò)程有限元計(jì)算模型及模擬結(jié)果分析
    將有限元分析引入材料力學(xué)組合變形的教學(xué)探索
    大型壓濾機(jī)主梁的彈性力學(xué)計(jì)算及有限元分析設(shè)計(jì)
    大型壓濾機(jī)主梁的彈性力學(xué)計(jì)算及有限元分析設(shè)計(jì)
    顯式Runge—Kutta局部間斷Galerkin方法的穩(wěn)定性分析
    起重機(jī)吊臂振動(dòng)模態(tài)的有限元分析
    基于短纖維增強(qiáng)的復(fù)合氣壓砂輪基體性能研究
    有限元編程方法 第5版
    伊人久久国产一区二区| 日日啪夜夜爽| 午夜老司机福利片| 1024香蕉在线观看| 亚洲精品第二区| 国产亚洲最大av| 亚洲精品自拍成人| 丝瓜视频免费看黄片| 久久人人爽av亚洲精品天堂| 日韩,欧美,国产一区二区三区| 一二三四在线观看免费中文在| 操出白浆在线播放| 亚洲国产最新在线播放| 免费av中文字幕在线| 可以免费在线观看a视频的电影网站 | 老汉色∧v一级毛片| 精品久久蜜臀av无| 18禁观看日本| 制服人妻中文乱码| 精品国产露脸久久av麻豆| 精品国产一区二区三区久久久樱花| 又黄又粗又硬又大视频| 亚洲av福利一区| 亚洲一级一片aⅴ在线观看| 又粗又硬又长又爽又黄的视频| 国产精品秋霞免费鲁丝片| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩视频精品一区| 你懂的网址亚洲精品在线观看| 91精品三级在线观看| 国产色婷婷99| 少妇精品久久久久久久| 亚洲欧美清纯卡通| 日韩大片免费观看网站| 国产欧美亚洲国产| 丝袜在线中文字幕| 亚洲欧美中文字幕日韩二区| 国产97色在线日韩免费| 欧美日本中文国产一区发布| 午夜精品国产一区二区电影| 久久精品久久久久久久性| 亚洲精品国产色婷婷电影| 一级毛片电影观看| 这个男人来自地球电影免费观看 | 这个男人来自地球电影免费观看 | 一本大道久久a久久精品| 国产成人精品久久二区二区91 | 制服人妻中文乱码| 欧美精品一区二区大全| 亚洲专区中文字幕在线 | 巨乳人妻的诱惑在线观看| 免费看av在线观看网站| 一本久久精品| 亚洲国产精品一区三区| 精品久久蜜臀av无| 国产黄色免费在线视频| 久久99一区二区三区| 80岁老熟妇乱子伦牲交| 极品少妇高潮喷水抽搐| 女人爽到高潮嗷嗷叫在线视频| 中文字幕高清在线视频| 国产成人av激情在线播放| 两个人看的免费小视频| 日本av手机在线免费观看| 免费观看性生交大片5| 中文字幕av电影在线播放| 亚洲精品乱久久久久久| 久久性视频一级片| 成人亚洲欧美一区二区av| 免费少妇av软件| tube8黄色片| 亚洲国产精品国产精品| 一边摸一边做爽爽视频免费| 亚洲欧美中文字幕日韩二区| 777久久人妻少妇嫩草av网站| 老司机在亚洲福利影院| 久久人人爽av亚洲精品天堂| 久久青草综合色| 一二三四中文在线观看免费高清| 精品国产国语对白av| 亚洲av电影在线观看一区二区三区| 国产精品香港三级国产av潘金莲 | 丰满迷人的少妇在线观看| 丝袜人妻中文字幕| 亚洲精品美女久久久久99蜜臀 | 亚洲精品日本国产第一区| 两个人看的免费小视频| 最近最新中文字幕大全免费视频 | 亚洲人成网站在线观看播放| 美女中出高潮动态图| 国产精品久久久av美女十八| 日韩视频在线欧美| 久热这里只有精品99| 国产av一区二区精品久久| 国产日韩欧美亚洲二区| 久久久久精品性色| 日韩精品免费视频一区二区三区| 91老司机精品| 亚洲伊人久久精品综合| 国产精品成人在线| 久久久久久久久久久久大奶| 亚洲成人av在线免费| 欧美激情极品国产一区二区三区| 另类亚洲欧美激情| 乱人伦中国视频| 如何舔出高潮| 天天躁夜夜躁狠狠躁躁| 少妇 在线观看| 欧美国产精品va在线观看不卡| 国产免费视频播放在线视频| 亚洲综合色网址| 国产精品成人在线| 精品国产乱码久久久久久男人| 丝袜喷水一区| 免费日韩欧美在线观看| 老司机影院成人| 成年人午夜在线观看视频| 免费看av在线观看网站| 精品亚洲成国产av| 国产精品久久久av美女十八| av视频免费观看在线观看| 成年人免费黄色播放视频| 制服丝袜香蕉在线| 国产在视频线精品| 成人影院久久| av天堂久久9| 国产男女超爽视频在线观看| 99热全是精品| 亚洲欧美成人综合另类久久久| 亚洲精品乱久久久久久| 丰满乱子伦码专区| 黄片播放在线免费| 亚洲精品国产av成人精品| 亚洲成人一二三区av| 嫩草影视91久久| 亚洲国产毛片av蜜桃av| 日本91视频免费播放| 日韩精品免费视频一区二区三区| 国产精品二区激情视频| 久久精品人人爽人人爽视色| 97精品久久久久久久久久精品| 伊人久久大香线蕉亚洲五| 午夜福利影视在线免费观看| 亚洲av电影在线观看一区二区三区| 色视频在线一区二区三区| 免费不卡黄色视频| 99热网站在线观看| 久久99精品国语久久久| 欧美激情高清一区二区三区 | 亚洲精品第二区| 男人舔女人的私密视频| 亚洲人成网站在线观看播放| 亚洲中文av在线| 亚洲欧美一区二区三区黑人| 免费人妻精品一区二区三区视频| 看非洲黑人一级黄片| 大码成人一级视频| 国产福利在线免费观看视频| 欧美xxⅹ黑人| 一区二区三区精品91| 中文乱码字字幕精品一区二区三区| 国产精品无大码| 老司机影院毛片| 日本黄色日本黄色录像| 老汉色av国产亚洲站长工具| 亚洲精品美女久久av网站| 国产成人精品福利久久| 黄色一级大片看看| 亚洲中文av在线| 欧美日韩亚洲综合一区二区三区_| 午夜免费男女啪啪视频观看| 国产精品香港三级国产av潘金莲 | xxxhd国产人妻xxx| 国产日韩一区二区三区精品不卡| 日韩制服骚丝袜av| 看非洲黑人一级黄片| 日韩av免费高清视频| 丁香六月欧美| 一二三四在线观看免费中文在| 日日啪夜夜爽| 欧美人与善性xxx| 国产成人精品久久久久久| 欧美黑人欧美精品刺激| 久久国产精品男人的天堂亚洲| 午夜免费男女啪啪视频观看| 国产毛片在线视频| 欧美日韩综合久久久久久| 色精品久久人妻99蜜桃| 亚洲精品自拍成人| av福利片在线| 中文字幕色久视频| 91老司机精品| 国产老妇伦熟女老妇高清| 国产亚洲av高清不卡| 国产深夜福利视频在线观看| 日韩精品免费视频一区二区三区| 黄色视频不卡| 日韩大码丰满熟妇| av片东京热男人的天堂| 欧美激情极品国产一区二区三区| 中文字幕人妻丝袜一区二区 | 超色免费av| 少妇人妻久久综合中文| 久久久久人妻精品一区果冻| 黑人欧美特级aaaaaa片| 在线观看免费高清a一片| www.熟女人妻精品国产| 亚洲一级一片aⅴ在线观看| 美女扒开内裤让男人捅视频| 国产亚洲最大av| 卡戴珊不雅视频在线播放| 欧美国产精品va在线观看不卡| 欧美成人精品欧美一级黄| 久久狼人影院| 亚洲在久久综合| 久久久久久久久久久久大奶| 一区二区av电影网| 男女国产视频网站| 午夜福利视频精品| 9热在线视频观看99| 18禁动态无遮挡网站| 久久婷婷青草| 十八禁高潮呻吟视频| 在线天堂中文资源库| 国产一卡二卡三卡精品 | 日本欧美视频一区| 80岁老熟妇乱子伦牲交| 国产在线视频一区二区| 精品国产一区二区三区久久久樱花| 国产女主播在线喷水免费视频网站| 只有这里有精品99| 爱豆传媒免费全集在线观看| 人人澡人人妻人| 看免费av毛片| www日本在线高清视频| 国产精品久久久久久精品古装| 性色av一级| 午夜精品国产一区二区电影| avwww免费| 国产精品二区激情视频| 亚洲精品自拍成人| 亚洲av福利一区| 久久 成人 亚洲| 黄色视频在线播放观看不卡| 免费黄网站久久成人精品| 人人妻人人爽人人添夜夜欢视频| 亚洲精品一区蜜桃| 欧美在线黄色| 亚洲国产精品999| 国产av码专区亚洲av| 久久久久久久久久久久大奶| e午夜精品久久久久久久| 久久这里只有精品19| 久久久久久久久免费视频了| 久久久久精品性色| 久久久国产欧美日韩av| 亚洲中文av在线| 国产 一区精品| 亚洲精品国产色婷婷电影| 女性被躁到高潮视频| 新久久久久国产一级毛片| 国产在线免费精品| 久久性视频一级片| 欧美精品av麻豆av| 色婷婷久久久亚洲欧美| 曰老女人黄片| 中文字幕制服av| 亚洲av在线观看美女高潮| 男女床上黄色一级片免费看| 国产又色又爽无遮挡免| 满18在线观看网站| 日韩一本色道免费dvd| 男人操女人黄网站| 成人手机av| 日本午夜av视频| 国产亚洲午夜精品一区二区久久| 成人18禁高潮啪啪吃奶动态图| a级毛片黄视频| 一级黄片播放器| 日韩一区二区视频免费看| 久久99一区二区三区| 国产97色在线日韩免费| netflix在线观看网站| 纯流量卡能插随身wifi吗| 人妻一区二区av| 成人亚洲欧美一区二区av| 色精品久久人妻99蜜桃| 一区在线观看完整版| 免费在线观看完整版高清| 最近中文字幕2019免费版| 午夜福利网站1000一区二区三区| 亚洲一码二码三码区别大吗| 在线观看免费日韩欧美大片| 精品一区二区免费观看| 国产高清国产精品国产三级| 国产成人免费无遮挡视频| 亚洲五月色婷婷综合| 纯流量卡能插随身wifi吗| 日韩电影二区| 国产欧美亚洲国产| 国产精品国产三级专区第一集| 国产精品 欧美亚洲| 国产一区二区 视频在线| www日本在线高清视频| 日本91视频免费播放| 男女无遮挡免费网站观看| 欧美黑人精品巨大| 亚洲第一区二区三区不卡| 菩萨蛮人人尽说江南好唐韦庄| 一二三四中文在线观看免费高清| 午夜精品国产一区二区电影| 成人亚洲精品一区在线观看| 少妇人妻久久综合中文| 久久久久久久精品精品| 老司机亚洲免费影院| 免费黄频网站在线观看国产| 人体艺术视频欧美日本| 久久久久视频综合| 国精品久久久久久国模美| 亚洲精品第二区| 亚洲精品久久成人aⅴ小说| 制服丝袜香蕉在线| 高清av免费在线| 97人妻天天添夜夜摸| 一区福利在线观看| 男的添女的下面高潮视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲成人手机| 十八禁高潮呻吟视频| 国产av精品麻豆| 亚洲精品乱久久久久久| 不卡av一区二区三区| 免费少妇av软件| 欧美精品人与动牲交sv欧美| 精品一区在线观看国产| 黄网站色视频无遮挡免费观看| 深夜精品福利| 一区二区三区四区激情视频| 中文字幕人妻丝袜制服| 最黄视频免费看| 日日撸夜夜添| 亚洲av男天堂| 搡老乐熟女国产| 嫩草影院入口| 成人18禁高潮啪啪吃奶动态图| 亚洲国产最新在线播放| 天美传媒精品一区二区| 国产黄色视频一区二区在线观看| 99久久人妻综合| 免费观看人在逋| 国产一区二区在线观看av| 9191精品国产免费久久| 免费黄频网站在线观看国产| 亚洲国产精品一区二区三区在线| 国产男女内射视频| 狂野欧美激情性bbbbbb| 考比视频在线观看| 香蕉国产在线看| 大片免费播放器 马上看| tube8黄色片| 婷婷色麻豆天堂久久| 久久97久久精品| 在线观看人妻少妇| 国产女主播在线喷水免费视频网站| 免费av中文字幕在线| 五月天丁香电影| 日韩精品免费视频一区二区三区| 大话2 男鬼变身卡| 一级片'在线观看视频| 男男h啪啪无遮挡| 老司机靠b影院| 国产不卡av网站在线观看| 色精品久久人妻99蜜桃| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲自偷自拍图片 自拍| tube8黄色片| 国产av码专区亚洲av| 成人国产av品久久久| 日韩大码丰满熟妇| 国产人伦9x9x在线观看| 天堂俺去俺来也www色官网| 免费看不卡的av| 国产精品成人在线| 久久人人爽av亚洲精品天堂| 婷婷成人精品国产| 熟女av电影| 国产熟女午夜一区二区三区| av国产精品久久久久影院| 在线观看人妻少妇| 大话2 男鬼变身卡| 精品第一国产精品| 99久久人妻综合| 亚洲一码二码三码区别大吗| 成人影院久久| 成年av动漫网址| 欧美激情极品国产一区二区三区| av片东京热男人的天堂| 哪个播放器可以免费观看大片| 黄片播放在线免费| 亚洲国产欧美日韩在线播放| 久久久欧美国产精品| 只有这里有精品99| 天天躁日日躁夜夜躁夜夜| 免费黄频网站在线观看国产| 人体艺术视频欧美日本| 国产极品天堂在线| 啦啦啦中文免费视频观看日本| 伊人久久国产一区二区| 男女下面插进去视频免费观看| 黑人欧美特级aaaaaa片| 久久久久久人妻| av.在线天堂| 国产亚洲最大av| 波多野结衣av一区二区av| 日本色播在线视频| 午夜福利影视在线免费观看| h视频一区二区三区| 日本午夜av视频| 亚洲精品国产av成人精品| 日韩一区二区视频免费看| 久久国产精品大桥未久av| 久久久精品国产亚洲av高清涩受| 久久97久久精品| 欧美日韩亚洲高清精品| 日本黄色日本黄色录像| 一区在线观看完整版| 午夜免费鲁丝| 免费看av在线观看网站| 亚洲成av片中文字幕在线观看| 欧美亚洲 丝袜 人妻 在线| 黄片播放在线免费| 自线自在国产av| 交换朋友夫妻互换小说| 亚洲精品久久午夜乱码| av在线老鸭窝| 街头女战士在线观看网站| 婷婷色综合www| 精品少妇黑人巨大在线播放| 亚洲一级一片aⅴ在线观看| 建设人人有责人人尽责人人享有的| 中文字幕另类日韩欧美亚洲嫩草| 男女国产视频网站| 在线观看一区二区三区激情| 亚洲精品,欧美精品| 精品卡一卡二卡四卡免费| av国产精品久久久久影院| av在线老鸭窝| 欧美日韩视频精品一区| av不卡在线播放| 啦啦啦 在线观看视频| 啦啦啦在线观看免费高清www| 蜜桃国产av成人99| 亚洲国产日韩一区二区| 亚洲成国产人片在线观看| 97精品久久久久久久久久精品| 一区二区日韩欧美中文字幕| 国产又爽黄色视频| 如日韩欧美国产精品一区二区三区| 色视频在线一区二区三区| 精品人妻一区二区三区麻豆| 满18在线观看网站| 91老司机精品| 精品亚洲成a人片在线观看| 日韩av不卡免费在线播放| 日日爽夜夜爽网站| 波多野结衣av一区二区av| 日韩大码丰满熟妇| 99久久99久久久精品蜜桃| 国产精品香港三级国产av潘金莲 | 18禁裸乳无遮挡动漫免费视频| 精品一区二区三区av网在线观看 | 国产成人免费无遮挡视频| 国产人伦9x9x在线观看| 精品久久蜜臀av无| 久久久久久免费高清国产稀缺| 欧美日韩综合久久久久久| 亚洲欧洲精品一区二区精品久久久 | 亚洲国产最新在线播放| 亚洲第一青青草原| 日韩中文字幕视频在线看片| 韩国高清视频一区二区三区| 你懂的网址亚洲精品在线观看| 我的亚洲天堂| 精品一区二区三卡| 最近最新中文字幕免费大全7| a级毛片黄视频| 亚洲男人天堂网一区| 啦啦啦在线观看免费高清www| 少妇人妻精品综合一区二区| 国产在线免费精品| 自拍欧美九色日韩亚洲蝌蚪91| 侵犯人妻中文字幕一二三四区| 伦理电影大哥的女人| 亚洲成国产人片在线观看| 建设人人有责人人尽责人人享有的| 欧美激情极品国产一区二区三区| 女人精品久久久久毛片| 国产精品国产三级国产专区5o| 欧美日本中文国产一区发布| 国产精品三级大全| 亚洲精品一区蜜桃| 日本午夜av视频| 亚洲精品成人av观看孕妇| 美女高潮到喷水免费观看| 国产精品 欧美亚洲| videos熟女内射| 成人毛片60女人毛片免费| 最近中文字幕高清免费大全6| 80岁老熟妇乱子伦牲交| 十分钟在线观看高清视频www| 中文乱码字字幕精品一区二区三区| 久久久国产精品麻豆| 国产一区二区 视频在线| 欧美精品高潮呻吟av久久| av福利片在线| 日韩一区二区三区影片| 男女下面插进去视频免费观看| 欧美少妇被猛烈插入视频| 亚洲熟女毛片儿| 男的添女的下面高潮视频| 国语对白做爰xxxⅹ性视频网站| 国产一区二区激情短视频 | 久久久久久久大尺度免费视频| 电影成人av| 高清欧美精品videossex| 欧美日本中文国产一区发布| svipshipincom国产片| 精品人妻一区二区三区麻豆| 亚洲天堂av无毛| 97精品久久久久久久久久精品| 欧美日韩一区二区视频在线观看视频在线| 日本av手机在线免费观看| av视频免费观看在线观看| 亚洲一级一片aⅴ在线观看| 久久久久久人妻| 两性夫妻黄色片| a级毛片黄视频| 亚洲七黄色美女视频| 美国免费a级毛片| 亚洲,欧美精品.| 亚洲综合色网址| 999久久久国产精品视频| 麻豆av在线久日| 亚洲四区av| 丝袜美腿诱惑在线| 王馨瑶露胸无遮挡在线观看| 91精品三级在线观看| 亚洲情色 制服丝袜| 一区福利在线观看| 久久久国产一区二区| 搡老岳熟女国产| 爱豆传媒免费全集在线观看| 国产乱人偷精品视频| 久久精品久久久久久噜噜老黄| 国产精品蜜桃在线观看| 在线观看免费日韩欧美大片| 国产精品久久久久久精品电影小说| 国产精品av久久久久免费| 亚洲精品久久午夜乱码| 欧美老熟妇乱子伦牲交| 大片电影免费在线观看免费| 国产精品久久久av美女十八| 99久久精品国产亚洲精品| 狠狠精品人妻久久久久久综合| 欧美人与性动交α欧美精品济南到| 麻豆精品久久久久久蜜桃| 只有这里有精品99| 97精品久久久久久久久久精品| 交换朋友夫妻互换小说| av不卡在线播放| 黑丝袜美女国产一区| 99九九在线精品视频| 精品酒店卫生间| 黑人欧美特级aaaaaa片| 久久av网站| 中国国产av一级| 深夜精品福利| 丰满迷人的少妇在线观看| 十八禁高潮呻吟视频| 欧美日韩一级在线毛片| 悠悠久久av| 精品国产乱码久久久久久小说| 婷婷成人精品国产| 欧美亚洲 丝袜 人妻 在线| 韩国av在线不卡| 丝袜人妻中文字幕| 国产乱来视频区| 十八禁人妻一区二区| 久久久精品区二区三区| 制服诱惑二区| 另类精品久久| 久久久亚洲精品成人影院| 久久久久网色| 精品国产一区二区三区久久久樱花| xxxhd国产人妻xxx| 深夜精品福利| 一本色道久久久久久精品综合| 婷婷成人精品国产| 国产精品偷伦视频观看了| 美女中出高潮动态图| 各种免费的搞黄视频| 国产淫语在线视频| xxx大片免费视频| 亚洲国产毛片av蜜桃av| 午夜激情久久久久久久| 欧美av亚洲av综合av国产av | 下体分泌物呈黄色| 精品福利永久在线观看| 在线 av 中文字幕| 久久久久久久久久久久大奶| 一级a爱视频在线免费观看| 最新的欧美精品一区二区| 国产亚洲av高清不卡| 国产精品无大码| 晚上一个人看的免费电影|