• <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版
    别揉我奶头 嗯啊视频| 搡老岳熟女国产| 国产精品美女特级片免费视频播放器| h日本视频在线播放| 久久6这里有精品| 中国美白少妇内射xxxbb| 99九九线精品视频在线观看视频| 日韩一区二区视频免费看| 美女高潮喷水抽搐中文字幕| 观看免费一级毛片| 国产精品免费一区二区三区在线| 亚洲天堂国产精品一区在线| 国产私拍福利视频在线观看| 精品久久久久久成人av| 午夜福利18| xxxwww97欧美| 亚洲精品日韩av片在线观看| 91精品国产九色| 一a级毛片在线观看| 九九爱精品视频在线观看| 久久这里只有精品中国| 中文字幕免费在线视频6| 国产 一区精品| 欧美成人免费av一区二区三区| www日本黄色视频网| 亚洲久久久久久中文字幕| 九九久久精品国产亚洲av麻豆| 在线播放国产精品三级| 精品国产三级普通话版| 天堂影院成人在线观看| 国产男人的电影天堂91| 国产av麻豆久久久久久久| 亚洲av日韩精品久久久久久密| 国产精品人妻久久久久久| 日本免费a在线| 欧美不卡视频在线免费观看| 国模一区二区三区四区视频| 97热精品久久久久久| 免费观看人在逋| 日韩欧美三级三区| 午夜福利在线在线| 亚洲经典国产精华液单| 久久国产精品人妻蜜桃| 欧美一区二区国产精品久久精品| 成人高潮视频无遮挡免费网站| 97热精品久久久久久| 中文字幕av成人在线电影| 亚洲精品456在线播放app | 我的女老师完整版在线观看| 99精品在免费线老司机午夜| 听说在线观看完整版免费高清| 精品国产三级普通话版| 黄片wwwwww| 日本免费一区二区三区高清不卡| 欧美激情久久久久久爽电影| 日日夜夜操网爽| 人妻久久中文字幕网| 性欧美人与动物交配| 久久久久久久久大av| 国产精品人妻久久久久久| 国模一区二区三区四区视频| 久久精品人妻少妇| 亚洲男人的天堂狠狠| 91麻豆av在线| 久久精品人妻少妇| 免费看日本二区| 国产精品国产三级国产av玫瑰| 看黄色毛片网站| 一个人观看的视频www高清免费观看| 精品一区二区三区视频在线观看免费| x7x7x7水蜜桃| 久久精品国产亚洲av香蕉五月| 精品国内亚洲2022精品成人| 综合色av麻豆| 又紧又爽又黄一区二区| 无人区码免费观看不卡| 国内揄拍国产精品人妻在线| 啦啦啦观看免费观看视频高清| 免费观看人在逋| 不卡一级毛片| 特大巨黑吊av在线直播| 国产高清视频在线观看网站| 国产探花极品一区二区| 真人一进一出gif抽搐免费| 国产精品无大码| 日本一本二区三区精品| 日本免费a在线| 一进一出抽搐gif免费好疼| 亚洲av电影不卡..在线观看| 国产男人的电影天堂91| 天天一区二区日本电影三级| 精品一区二区三区视频在线| 日本五十路高清| 十八禁国产超污无遮挡网站| 九九热线精品视视频播放| 精品不卡国产一区二区三区| 欧美丝袜亚洲另类 | aaaaa片日本免费| 精品乱码久久久久久99久播| 欧美激情国产日韩精品一区| 一本一本综合久久| 久久精品国产自在天天线| 日韩欧美在线乱码| 最后的刺客免费高清国语| 亚洲精品一区av在线观看| 淫妇啪啪啪对白视频| 欧美zozozo另类| 神马国产精品三级电影在线观看| 日本在线视频免费播放| 黄色视频,在线免费观看| 成年女人看的毛片在线观看| 欧美成人免费av一区二区三区| 亚洲无线在线观看| 全区人妻精品视频| 精品99又大又爽又粗少妇毛片 | 亚洲精品乱码久久久v下载方式| a级毛片免费高清观看在线播放| 小说图片视频综合网站| 亚洲av不卡在线观看| 亚洲无线在线观看| 国产精品久久久久久av不卡| 亚洲四区av| 久久久国产成人精品二区| 内地一区二区视频在线| 日韩欧美在线二视频| 国产精品精品国产色婷婷| 一进一出抽搐gif免费好疼| 久久人妻av系列| 亚洲精华国产精华精| 成人三级黄色视频| 国产亚洲精品综合一区在线观看| 男人和女人高潮做爰伦理| 久久精品国产亚洲av香蕉五月| 少妇熟女aⅴ在线视频| 亚洲成a人片在线一区二区| 女生性感内裤真人,穿戴方法视频| 国产熟女欧美一区二区| 免费黄网站久久成人精品| 91麻豆精品激情在线观看国产| 久久久久久久久久成人| 亚洲人成伊人成综合网2020| 99热只有精品国产| 国产一区二区亚洲精品在线观看| 亚州av有码| 亚洲综合色惰| 久久久成人免费电影| 亚洲美女视频黄频| 久久精品久久久久久噜噜老黄 | 欧美丝袜亚洲另类 | 国产aⅴ精品一区二区三区波| 中亚洲国语对白在线视频| 热99在线观看视频| 亚洲经典国产精华液单| 别揉我奶头 嗯啊视频| 亚洲精华国产精华精| 他把我摸到了高潮在线观看| 亚洲成人中文字幕在线播放| 无遮挡黄片免费观看| 亚洲自拍偷在线| 美女免费视频网站| 18禁黄网站禁片免费观看直播| 亚洲精品日韩av片在线观看| 一a级毛片在线观看| 国产白丝娇喘喷水9色精品| 亚洲天堂国产精品一区在线| 亚洲av免费高清在线观看| 天堂√8在线中文| ponron亚洲| 成人精品一区二区免费| 久久久久性生活片| 一夜夜www| 草草在线视频免费看| 在线天堂最新版资源| 亚洲在线自拍视频| 别揉我奶头~嗯~啊~动态视频| 长腿黑丝高跟| 久久久午夜欧美精品| 国产免费男女视频| 色av中文字幕| 中亚洲国语对白在线视频| 国产精品免费一区二区三区在线| 999久久久精品免费观看国产| 中文字幕久久专区| 在线免费观看不下载黄p国产 | av黄色大香蕉| 亚洲欧美日韩高清专用| 琪琪午夜伦伦电影理论片6080| 亚洲精品一区av在线观看| 成人高潮视频无遮挡免费网站| 亚洲美女搞黄在线观看 | 欧美一区二区国产精品久久精品| 精品人妻1区二区| 亚洲最大成人中文| 又爽又黄a免费视频| 亚洲欧美日韩高清专用| 国产精品精品国产色婷婷| 国产一区二区在线观看日韩| 小说图片视频综合网站| 免费观看精品视频网站| 狂野欧美白嫩少妇大欣赏| 精品免费久久久久久久清纯| 国产伦精品一区二区三区视频9| 永久网站在线| 亚洲成av人片在线播放无| 蜜桃亚洲精品一区二区三区| 欧美xxxx黑人xx丫x性爽| 国产精品人妻久久久影院| 亚洲精品国产成人久久av| 精品久久久久久久久亚洲 | 日韩强制内射视频| 欧美bdsm另类| 欧美三级亚洲精品| 色综合婷婷激情| 国产久久久一区二区三区| 亚洲国产色片| 观看美女的网站| 国产私拍福利视频在线观看| 91午夜精品亚洲一区二区三区 | 成年女人看的毛片在线观看| 午夜日韩欧美国产| 中文字幕免费在线视频6| 1024手机看黄色片| 亚洲性夜色夜夜综合| 日韩欧美精品v在线| 一个人看的www免费观看视频| 国产又黄又爽又无遮挡在线| 一本一本综合久久| or卡值多少钱| 亚洲成人中文字幕在线播放| 亚洲专区国产一区二区| 亚洲av成人精品一区久久| 免费av观看视频| 国产免费一级a男人的天堂| 国产精品一区二区三区四区免费观看 | 国产又黄又爽又无遮挡在线| 看十八女毛片水多多多| 日韩一本色道免费dvd| 久久中文看片网| 中亚洲国语对白在线视频| 国产av一区在线观看免费| 国内精品一区二区在线观看| 国产极品精品免费视频能看的| 日本五十路高清| 国产私拍福利视频在线观看| 国产精品免费一区二区三区在线| 久久欧美精品欧美久久欧美| 日本成人三级电影网站| 日日摸夜夜添夜夜添av毛片 | 一区二区三区激情视频| 欧美日韩国产亚洲二区| 九色国产91popny在线| 国产精品久久久久久久电影| 久久人人爽人人爽人人片va| 色精品久久人妻99蜜桃| 亚洲性夜色夜夜综合| 美女高潮喷水抽搐中文字幕| 91av网一区二区| 国产不卡一卡二| 国产精品久久久久久av不卡| 亚洲av不卡在线观看| 老熟妇乱子伦视频在线观看| 精品人妻视频免费看| 欧美黑人巨大hd| 日本三级黄在线观看| 国产视频一区二区在线看| 国产亚洲精品综合一区在线观看| 尤物成人国产欧美一区二区三区| xxxwww97欧美| 国产精品一区www在线观看 | 老女人水多毛片| 老司机福利观看| 亚洲成a人片在线一区二区| 亚洲国产精品久久男人天堂| 联通29元200g的流量卡| 日本与韩国留学比较| www.色视频.com| 国内毛片毛片毛片毛片毛片| 久久九九热精品免费| 看片在线看免费视频| 欧美人与善性xxx| 精品一区二区三区视频在线观看免费| 精品人妻视频免费看| 少妇的逼好多水| 十八禁网站免费在线| 97超视频在线观看视频| 91麻豆精品激情在线观看国产| 亚洲四区av| 免费av观看视频| 欧美日韩综合久久久久久 | 国产白丝娇喘喷水9色精品| 亚洲va日本ⅴa欧美va伊人久久| 国产亚洲精品av在线| 波多野结衣巨乳人妻| 日日夜夜操网爽| 黄色欧美视频在线观看| 成人高潮视频无遮挡免费网站| 亚洲精品日韩av片在线观看| 天堂动漫精品| 一区二区三区四区激情视频 | 国产精品亚洲一级av第二区| 国产91精品成人一区二区三区| 国产亚洲精品久久久com| 亚洲中文日韩欧美视频| 亚洲真实伦在线观看| 国产精品久久久久久精品电影| 国产高清激情床上av| 国产精品av视频在线免费观看| 88av欧美| 国产中年淑女户外野战色| 国产精品一区二区免费欧美| 又黄又爽又刺激的免费视频.| 又黄又爽又刺激的免费视频.| 国产高清三级在线| 一区二区三区激情视频| 精品不卡国产一区二区三区| 麻豆成人av在线观看| 女同久久另类99精品国产91| 精品不卡国产一区二区三区| 三级男女做爰猛烈吃奶摸视频| 免费av不卡在线播放| 麻豆成人av在线观看| 国产熟女欧美一区二区| 久久精品91蜜桃| 1000部很黄的大片| 国产高清视频在线播放一区| 免费在线观看影片大全网站| 国产免费一级a男人的天堂| 国产高潮美女av| 亚洲av中文av极速乱 | 在线国产一区二区在线| 日本色播在线视频| 亚洲精华国产精华精| 99riav亚洲国产免费| 国内精品久久久久精免费| 美女免费视频网站| 久久久久久久久久黄片| 十八禁网站免费在线| 久久久久久国产a免费观看| a级毛片a级免费在线| 日本免费a在线| 一区二区三区四区激情视频 | 女的被弄到高潮叫床怎么办 | 又粗又爽又猛毛片免费看| 精品人妻1区二区| 久久国产乱子免费精品| 国产主播在线观看一区二区| 免费电影在线观看免费观看| 久久精品影院6| 国产精品自产拍在线观看55亚洲| 三级毛片av免费| 亚洲国产精品合色在线| 国产精品人妻久久久影院| 欧美激情在线99| 男人的好看免费观看在线视频| 久久久国产成人免费| 两人在一起打扑克的视频| 女的被弄到高潮叫床怎么办 | 男人舔奶头视频| 18禁黄网站禁片免费观看直播| 久久99热这里只有精品18| 美女高潮喷水抽搐中文字幕| 亚洲国产欧美人成| 一个人免费在线观看电影| 久久午夜福利片| 简卡轻食公司| 成年女人毛片免费观看观看9| 日本在线视频免费播放| 免费在线观看日本一区| 亚洲成av人片在线播放无| 国产毛片a区久久久久| 日韩中文字幕欧美一区二区| 久久热精品热| 五月伊人婷婷丁香| 在线免费观看的www视频| 亚洲美女黄片视频| www日本黄色视频网| 亚洲一区二区三区色噜噜| 97碰自拍视频| 婷婷六月久久综合丁香| 亚洲精品粉嫩美女一区| 嫩草影院新地址| 国产精华一区二区三区| 韩国av在线不卡| 欧美国产日韩亚洲一区| 日韩一区二区视频免费看| 国内精品一区二区在线观看| 国产一区二区三区视频了| 村上凉子中文字幕在线| 在线免费观看的www视频| 国产精品精品国产色婷婷| www.www免费av| 免费在线观看成人毛片| 亚洲国产精品久久男人天堂| 国产激情偷乱视频一区二区| 白带黄色成豆腐渣| 97碰自拍视频| 精品一区二区免费观看| 舔av片在线| 最后的刺客免费高清国语| 很黄的视频免费| 国产 一区 欧美 日韩| 成人永久免费在线观看视频| 国产视频一区二区在线看| 国产男靠女视频免费网站| 亚洲av.av天堂| 制服丝袜大香蕉在线| 日韩亚洲欧美综合| 亚洲av电影不卡..在线观看| 国产精品野战在线观看| 波野结衣二区三区在线| 国产av在哪里看| 国产成年人精品一区二区| 男女之事视频高清在线观看| 久9热在线精品视频| 亚洲国产精品久久男人天堂| netflix在线观看网站| 久久精品91蜜桃| 男人舔女人下体高潮全视频| 国产精品日韩av在线免费观看| 欧美xxxx黑人xx丫x性爽| 成年免费大片在线观看| 美女被艹到高潮喷水动态| 婷婷亚洲欧美| 午夜免费成人在线视频| 亚州av有码| 一个人观看的视频www高清免费观看| 国产成人av教育| 十八禁国产超污无遮挡网站| 免费看光身美女| 亚洲精品在线观看二区| 亚洲五月天丁香| 亚洲国产日韩欧美精品在线观看| 国产男人的电影天堂91| 国产人妻一区二区三区在| 亚洲,欧美,日韩| 欧美bdsm另类| 春色校园在线视频观看| 国产在线男女| 少妇丰满av| 欧美成人一区二区免费高清观看| 精品日产1卡2卡| 搞女人的毛片| 国产av在哪里看| a级毛片免费高清观看在线播放| 一本一本综合久久| 最近最新中文字幕大全电影3| 亚洲av五月六月丁香网| 91av网一区二区| 国内毛片毛片毛片毛片毛片| 三级毛片av免费| av专区在线播放| 亚洲成人久久性| 免费高清视频大片| 免费观看的影片在线观看| 少妇丰满av| 十八禁网站免费在线| 一本一本综合久久| 欧美日韩亚洲国产一区二区在线观看| 色5月婷婷丁香| 国产三级在线视频| 中文字幕熟女人妻在线| 亚洲成人久久爱视频| 人妻丰满熟妇av一区二区三区| 日本-黄色视频高清免费观看| 我的女老师完整版在线观看| 日韩欧美精品免费久久| 久久久久久久久久久丰满 | 身体一侧抽搐| 国产亚洲91精品色在线| 国产淫片久久久久久久久| 人妻丰满熟妇av一区二区三区| av.在线天堂| 国产精品国产三级国产av玫瑰| 性色avwww在线观看| 老女人水多毛片| 日本-黄色视频高清免费观看| 97超视频在线观看视频| 亚洲四区av| 免费人成在线观看视频色| 国产爱豆传媒在线观看| 亚洲熟妇熟女久久| 在线免费十八禁| 国产伦精品一区二区三区四那| 国产精品综合久久久久久久免费| 国产黄a三级三级三级人| 熟女人妻精品中文字幕| 人人妻人人看人人澡| 欧美又色又爽又黄视频| 一级毛片久久久久久久久女| 免费看av在线观看网站| 成人高潮视频无遮挡免费网站| 少妇丰满av| 最后的刺客免费高清国语| 国产高清三级在线| 99久久久亚洲精品蜜臀av| 国产精品三级大全| aaaaa片日本免费| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲成人中文字幕在线播放| 毛片一级片免费看久久久久 | 色综合婷婷激情| 亚洲欧美日韩无卡精品| 亚洲黑人精品在线| 乱码一卡2卡4卡精品| 麻豆久久精品国产亚洲av| 亚洲av日韩精品久久久久久密| 99热这里只有是精品在线观看| 不卡视频在线观看欧美| 久久99热这里只有精品18| 一本一本综合久久| 看黄色毛片网站| 亚洲精品成人久久久久久| 午夜老司机福利剧场| 在线天堂最新版资源| 午夜免费激情av| 国产一区二区三区视频了| 麻豆成人午夜福利视频| 22中文网久久字幕| 欧美一区二区亚洲| 亚洲图色成人| 日韩欧美 国产精品| 免费看美女性在线毛片视频| 亚洲国产精品久久男人天堂| 国内精品美女久久久久久| 91久久精品国产一区二区三区| 久久草成人影院| 在线天堂最新版资源| 综合色av麻豆| 亚洲美女搞黄在线观看 | 中国美白少妇内射xxxbb| 亚洲成人中文字幕在线播放| 俺也久久电影网| 人人妻人人看人人澡| 男女下面进入的视频免费午夜| 给我免费播放毛片高清在线观看| 看黄色毛片网站| 搡女人真爽免费视频火全软件 | 99国产极品粉嫩在线观看| 日韩av在线大香蕉| 99热精品在线国产| 国产 一区精品| 嫁个100分男人电影在线观看| 欧美高清成人免费视频www| 久久亚洲精品不卡| 亚洲成av人片在线播放无| 久久亚洲精品不卡| 在线看三级毛片| 男女那种视频在线观看| 免费看a级黄色片| 日韩av在线大香蕉| 窝窝影院91人妻| 色噜噜av男人的天堂激情| 久9热在线精品视频| 亚洲最大成人中文| АⅤ资源中文在线天堂| 亚洲av不卡在线观看| .国产精品久久| or卡值多少钱| 乱码一卡2卡4卡精品| 村上凉子中文字幕在线| 美女cb高潮喷水在线观看| 熟妇人妻久久中文字幕3abv| 五月玫瑰六月丁香| 特大巨黑吊av在线直播| 久久久久久九九精品二区国产| 岛国在线免费视频观看| 啦啦啦韩国在线观看视频| 国产一区二区在线观看日韩| 久99久视频精品免费| 亚洲专区国产一区二区| 极品教师在线免费播放| 精品久久久久久久久久免费视频| 成人三级黄色视频| 午夜爱爱视频在线播放| 亚洲专区国产一区二区| 偷拍熟女少妇极品色| 国产色婷婷99| 亚洲国产欧洲综合997久久,| 久久国产精品人妻蜜桃| 久久久久久国产a免费观看| 中文字幕av在线有码专区| 九九热线精品视视频播放| 成人精品一区二区免费| 久久精品综合一区二区三区| 国产av麻豆久久久久久久| 精品久久久久久久久av| 中文字幕久久专区| 九色成人免费人妻av| 免费人成在线观看视频色| 国产精品人妻久久久影院| 亚洲图色成人| 欧美日韩精品成人综合77777| 非洲黑人性xxxx精品又粗又长| 999久久久精品免费观看国产| 99久久无色码亚洲精品果冻| 欧美+日韩+精品| 亚洲精品国产成人久久av| 成人综合一区亚洲| 国产私拍福利视频在线观看| 校园春色视频在线观看| 午夜福利18| 日本撒尿小便嘘嘘汇集6| 午夜福利高清视频| 国内精品宾馆在线| 欧美日韩综合久久久久久 | 亚洲七黄色美女视频| 在线观看舔阴道视频| 琪琪午夜伦伦电影理论片6080| 欧美一级a爱片免费观看看| 国产精品久久久久久精品电影| 一卡2卡三卡四卡精品乱码亚洲| 亚洲专区中文字幕在线| 美女高潮的动态| 搡女人真爽免费视频火全软件 | 日本欧美国产在线视频|