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

    基于ANSYS的某型航空發(fā)動機(jī)軸承試驗(yàn)器支承剛度研究

    2014-07-09 18:34楊法立廖明夫王四季蔣云帆
    現(xiàn)代電子技術(shù) 2014年13期
    關(guān)鍵詞:結(jié)構(gòu)單元航空發(fā)動機(jī)

    楊法立+廖明夫+王四季+蔣云帆

    摘 要: 某型航空發(fā)動機(jī)軸承試驗(yàn)器采用雙“L”形安裝邊支承結(jié)構(gòu),采用ANSYS中solid45,solid95及solid187結(jié)構(gòu)單元來計(jì)算其支承剛度,并通過傳遞矩陣法分別計(jì)算軸承試驗(yàn)器的臨界轉(zhuǎn)速。通過與軸承試驗(yàn)器的實(shí)際臨界轉(zhuǎn)速對比,分析三種單元結(jié)構(gòu)對于該支承結(jié)構(gòu)剛度的計(jì)算準(zhǔn)確性,結(jié)果表明采用solid187結(jié)構(gòu)單元計(jì)算的臨界轉(zhuǎn)速誤差最小,可用于此種支承結(jié)構(gòu)的靜態(tài)剛度仿真計(jì)算。

    關(guān)鍵詞: 結(jié)構(gòu)單元; 支承剛度; 臨界轉(zhuǎn)速; 航空發(fā)動機(jī)

    中圖分類號: TN710?34;V231.96 文獻(xiàn)標(biāo)識碼: A 文章編號: 1004?373X(2014)13?0116?03

    Research on supporting stiffness of an aero?engine bearing tester based on ANSYS

    YANG Fa?li, LIAO Ming?fu, WANG Si?ji, JIANG Yun?fan

    (School of Power and Energy, Northwestern Polytechnical University, Xian 710072, China)

    Abstract: The supporting structure with the double ‘L shape assembling sides is adopted for an aero?engine bearing tester. The solid45, solid95 and solid187 structural units in ANSYS are used to calculate the supporting stiffness of the tester. The critical rotation speed of the bearing tester is calculated with the transfer matrix method. The calculation accuracy of three unit structures is analyzed for stiffness of the supporting structure by comparing it with the real critical rotation speed. The results show that the critical rotation speed got by solid187 has the smallest computation error, so the element type of solid187 can be applied to the simulation calculation of static stiffness of the supporting structure.

    Keywords: structural unit; supporting stiffness; critical rotation speed; aero?engine

    0 引 言

    在工程實(shí)踐中對于轉(zhuǎn)子臨界轉(zhuǎn)速的計(jì)算,已經(jīng)擁有了很多種計(jì)算方法。總結(jié)起來,幾乎所有的計(jì)算方法都要考慮到轉(zhuǎn)子的質(zhì)量、阻尼和剛度因素對臨界轉(zhuǎn)速的影響[1]。轉(zhuǎn)子支承剛度是影響轉(zhuǎn)子臨界轉(zhuǎn)速的重要因素之一[2?4],因此在進(jìn)行轉(zhuǎn)子設(shè)計(jì)時(shí),首先需要確定轉(zhuǎn)子的支承剛度,然后預(yù)估轉(zhuǎn)子的臨界轉(zhuǎn)速,保證設(shè)計(jì)轉(zhuǎn)子的動力學(xué)特性滿足設(shè)計(jì)要求。

    本文研究的為某型航空發(fā)動機(jī)軸承試驗(yàn)器,其軸承座結(jié)構(gòu)與真實(shí)發(fā)動機(jī)結(jié)構(gòu)相似,采用雙“L”形安裝邊與支承底座相連。為預(yù)估本轉(zhuǎn)子的臨界轉(zhuǎn)速,必須預(yù)估其支承結(jié)構(gòu)的支承剛度。本文研究了ANSYS軟件中solid45、solid95及solid187這三種結(jié)構(gòu)單元用于計(jì)算此類支承結(jié)構(gòu)的支承剛度的適用性,并通過與實(shí)驗(yàn)結(jié)果對比,分析了這三種結(jié)構(gòu)單元用于計(jì)算支承剛度的準(zhǔn)確性。

    1 支承剛度仿真計(jì)算

    軸承試驗(yàn)器采用的是雙支承單盤單轉(zhuǎn)子結(jié)構(gòu),轉(zhuǎn)子的兩個支承分別為:試驗(yàn)軸承端,用于研究高轉(zhuǎn)速下軸承的動力學(xué)特性;陪試軸承端,為轉(zhuǎn)子系統(tǒng)提供合理的支點(diǎn)位置和支承剛度。圖1所示為軸承試驗(yàn)器三維圖,圖中標(biāo)明了“L”形安裝邊的結(jié)構(gòu)位置。

    圖1 轉(zhuǎn)子試驗(yàn)器三維圖

    對與該轉(zhuǎn)子試驗(yàn)器需要計(jì)算其臨界轉(zhuǎn)速,以便在其動力學(xué)特性滿足試驗(yàn)要求。為了得到可靠合理的計(jì)算結(jié)果,必須對軸承試驗(yàn)器兩端的支承結(jié)構(gòu)進(jìn)行較為準(zhǔn)確剛度預(yù)估。

    1.1 支承結(jié)構(gòu)單元

    根據(jù)查閱相關(guān)文獻(xiàn),有學(xué)者針對solid45單元類型[5],用該單元對支承結(jié)構(gòu)進(jìn)行航空發(fā)動機(jī)機(jī)匣結(jié)構(gòu)的建模分析,該單元通過8個節(jié)點(diǎn)來定義,每個節(jié)點(diǎn)有3個沿著[xyz]方向平移的自由度;有的學(xué)者針對solid95單元模型[6],該單元是solid45的更高階單元,為三維20節(jié)點(diǎn)四面體,在保證精度的同時(shí),允許使用不規(guī)則的形狀,適用于模擬曲線的邊界,該單元具有空間的任意方向,同時(shí)該單元還具有可塑性,蠕動,應(yīng)力鋼化,大變形和大應(yīng)變能力。用該單元進(jìn)行軸承座結(jié)構(gòu)的分析;最新版本的ANSYS軟件,綜合原來實(shí)體模型新增了一種高階單元solid187,該單元是帶中間節(jié)點(diǎn)的四面體,具有二次位移模式可以更好的模擬不規(guī)則模型,理論上可以得到更高的計(jì)算精度。

    本文將三種結(jié)構(gòu)單元對兩端支承結(jié)構(gòu)進(jìn)行靜剛度仿真分析,首先在UG中對支承結(jié)構(gòu)進(jìn)行三維建模,然后將三維模型轉(zhuǎn)為PARASOLID格式導(dǎo)入ANSYS軟件中進(jìn)行靜態(tài)受力分析,采用智能網(wǎng)格劃分方法,選取網(wǎng)格精度等級1。

    由于各結(jié)構(gòu)力學(xué)參數(shù)基本一致,因此采用同一材料參數(shù),材料參數(shù)為:泊松比[μ=0.3;]彈性模量[E=][2.08×1011] Pa。

    1.2 試驗(yàn)軸承端剛度仿真計(jì)算

    試驗(yàn)軸承端用于研究軸承的結(jié)構(gòu)、工藝參數(shù)對軸承支承動力學(xué)特性的影響,其支承結(jié)構(gòu)形式與發(fā)動機(jī)高壓轉(zhuǎn)子前支承一致,因此結(jié)構(gòu)較為復(fù)雜,主要結(jié)構(gòu)有軸承座、安裝環(huán)、“L”形安裝邊、支承底座。對支承結(jié)構(gòu)進(jìn)行整體網(wǎng)格劃分后如圖2所示。

    圖2 試驗(yàn)軸承端網(wǎng)格劃分

    考慮到支承底座通過四個地腳螺栓與試驗(yàn)臺固定,因此在進(jìn)行計(jì)算時(shí),將支承底座與試驗(yàn)臺之間的配合面施加[UX=0,][UY=0,][Z=0]的約束。分別在軸承配合面上施加豎直向上和豎直向下1 500 N的載荷,通過變形量來計(jì)算支承結(jié)構(gòu)的剛度。

    圖3為sloid95結(jié)構(gòu)單元施加載荷后的變形圖,另外兩個結(jié)構(gòu)單元仿真結(jié)果圖與其相似,僅具體的變形量不同,三種結(jié)構(gòu)單元的仿真計(jì)算結(jié)果見表1。

    1.3 陪試軸承端剛度仿真計(jì)算

    陪試軸承端的作用主要為轉(zhuǎn)子系統(tǒng)提供合理的支點(diǎn)位置和支承剛度,為了保證試驗(yàn)端試驗(yàn)數(shù)據(jù)的可分析性,要求陪試軸承端的結(jié)構(gòu)簡單,以減少對試驗(yàn)結(jié)果的影響,主要結(jié)構(gòu)有軸承座、“L”形安裝邊、支承底座。對支承結(jié)構(gòu)進(jìn)行整體網(wǎng)格劃分后如圖4所示。

    圖3 試驗(yàn)軸承端仿真結(jié)果圖

    圖4 陪試軸承端網(wǎng)格劃分

    由于支承底座同樣通過四個地腳螺栓與試驗(yàn)臺固定,因此也將支承底座與試驗(yàn)臺之間的配合面施加[UX=]0,[UY=0,][Z=0]的約束,分別在軸承配合面上施加豎直向上和豎直向下[1 500]N的載荷,通過變形量來計(jì)算支承結(jié)構(gòu)的剛度。

    圖5為sloid187結(jié)構(gòu)單元施加載荷后的變形圖,另外兩個結(jié)構(gòu)單元仿真結(jié)果圖與其相似,僅具體的變形量不同,三種結(jié)構(gòu)單元的仿真計(jì)算結(jié)果見表2。

    圖5 陪試軸承端仿真結(jié)果圖

    2 動力學(xué)特性預(yù)估

    根據(jù)前面預(yù)估的軸承支承剛度值,采用傳遞矩陣法對轉(zhuǎn)子的動力學(xué)特性進(jìn)行預(yù)估,傳遞矩陣法的基本原理是,取不同的轉(zhuǎn)速值,循環(huán)進(jìn)行各軸段截面狀態(tài)參數(shù)的逐段推算,直至滿足轉(zhuǎn)軸另一端的邊界條件[7?8]。本文根據(jù)轉(zhuǎn)軸不同的截面直徑、支承位置及轉(zhuǎn)盤位置,將轉(zhuǎn)軸離散為14個軸段、2個彈性支承站和1個盤站,圖6所示為軸承試驗(yàn)器的計(jì)算模型。

    圖6 轉(zhuǎn)子計(jì)算模型

    將各個單元所計(jì)算剛度的平均值分別代入到計(jì)算模型中的彈性支承站,編程計(jì)算轉(zhuǎn)子臨界轉(zhuǎn)速,計(jì)算結(jié)果見表3。

    3 試驗(yàn)與仿真結(jié)果對比分析

    圖7所示為軸承試驗(yàn)器實(shí)物圖,軸承座上安裝速度傳感器測量軸承振動,盤上安裝位移傳感器測量轉(zhuǎn)子振動,通過測量各個通道的振動幅值來確定轉(zhuǎn)子的一階臨界轉(zhuǎn)速。

    圖8所示為軸承試驗(yàn)器增速過程,圖中顯示的為盤上的振動數(shù)據(jù),實(shí)驗(yàn)測得轉(zhuǎn)子的一階臨界轉(zhuǎn)速為5 222 r/min。

    圖7 軸承試驗(yàn)器實(shí)物圖

    將計(jì)算結(jié)果與實(shí)驗(yàn)測得的臨界轉(zhuǎn)速對比可知:采用solid45單元結(jié)構(gòu)進(jìn)行剛度仿真計(jì)算,得到的臨界轉(zhuǎn)速比實(shí)測值高398 r/min,誤差為7.6%;采用solid95結(jié)構(gòu)單元進(jìn)行剛度仿真計(jì)算,得到的臨界轉(zhuǎn)速比實(shí)測值高311,誤差為5.9%;采用solid187結(jié)構(gòu)單元進(jìn)行剛度仿真計(jì)算,得到的臨界轉(zhuǎn)速比實(shí)測值高168 r/min,誤差為3.2%。

    圖8 試驗(yàn)器增速Bode圖

    4 結(jié) 語

    通過仿真計(jì)算與實(shí)驗(yàn)對比分析可知,對于采用雙“L”安裝邊支承結(jié)構(gòu),利用solid87單元結(jié)構(gòu)進(jìn)行剛度預(yù)估,能夠得到較為準(zhǔn)確的臨界轉(zhuǎn)速計(jì)算結(jié)果,因此在今后的設(shè)計(jì)工作中可以采用solid187結(jié)構(gòu)單元對相似的支承結(jié)構(gòu)進(jìn)行剛度仿真計(jì)算。同時(shí)注意到采用這三種結(jié)構(gòu)單元計(jì)算出的臨界轉(zhuǎn)速全部偏大,若采用此種方法計(jì)算轉(zhuǎn)子的臨界轉(zhuǎn)速,需要在軸承試驗(yàn)器增速過程,注意臨界轉(zhuǎn)速的提前來到。

    參考文獻(xiàn)

    [1] 白中祥,吳偉亮.轉(zhuǎn)子支承系統(tǒng)的剛度對其臨界轉(zhuǎn)速的影響[J].機(jī)電設(shè)備,2012(5):46?49.

    [2] 洪杰,王華,肖大為,等.轉(zhuǎn)子支承動剛度對轉(zhuǎn)子動力特性的影響分析[J].航空發(fā)動機(jī),2008,34(1):23?27.

    [3] 程小勇,陳果,李成剛,等.航空發(fā)動機(jī)轉(zhuǎn)子系統(tǒng)模擬支承設(shè)計(jì)與剛度計(jì)算[J].航空計(jì)算技術(shù),2012,42(6):9?12.

    [4] 盛步云,張濤,丁毓峰,等.支承剛度對汽輪機(jī)轉(zhuǎn)子動力學(xué)特性的影響分析[J].機(jī)械設(shè)計(jì),2008,25(12):38?40.

    [5] 趙文濤,陳果,李瓊,等.航空發(fā)動機(jī)機(jī)匣支承動剛度有限元計(jì)算及驗(yàn)證[J].航空計(jì)算技術(shù),2011,41(5):34?38.

    [6] 何新榮,傅行軍.基于ANSYS的軸承?轉(zhuǎn)子系統(tǒng)動力特性研究[J].機(jī)械研究與應(yīng)用,2011(6):39?41.

    [7] 呂文林.航空發(fā)動機(jī)強(qiáng)度計(jì)算[M].北京:航空工業(yè)出版社,2000.

    [8] 徐龍祥.高速旋轉(zhuǎn)機(jī)械軸系動力學(xué)設(shè)計(jì)[M].北京:國防工業(yè)出版社,1994.

    [9] 馬輝,惠紅杰,唐玉生,等.不同單元計(jì)算轉(zhuǎn)子臨界轉(zhuǎn)速的對比分析[J].振動與沖擊,2012,31(z1):40?43.

    [10] 陳忠.滾動軸承及其支承的剛度計(jì)算[J].煤礦機(jī)械,2006,27(3):387?388.

    [11] 鐘芳明,龔建政,賀星.燃?xì)廨啓C(jī)支承剛度計(jì)算研究[J].船海工程,2012,41(1):41?44.

    [12] 王毅,吳立言,寒冰.ANSYS的兩種有限單元應(yīng)用研究[J].科學(xué)技術(shù)與工程,2007,7(6):955?958.

    1.2 試驗(yàn)軸承端剛度仿真計(jì)算

    試驗(yàn)軸承端用于研究軸承的結(jié)構(gòu)、工藝參數(shù)對軸承支承動力學(xué)特性的影響,其支承結(jié)構(gòu)形式與發(fā)動機(jī)高壓轉(zhuǎn)子前支承一致,因此結(jié)構(gòu)較為復(fù)雜,主要結(jié)構(gòu)有軸承座、安裝環(huán)、“L”形安裝邊、支承底座。對支承結(jié)構(gòu)進(jìn)行整體網(wǎng)格劃分后如圖2所示。

    圖2 試驗(yàn)軸承端網(wǎng)格劃分

    考慮到支承底座通過四個地腳螺栓與試驗(yàn)臺固定,因此在進(jìn)行計(jì)算時(shí),將支承底座與試驗(yàn)臺之間的配合面施加[UX=0,][UY=0,][Z=0]的約束。分別在軸承配合面上施加豎直向上和豎直向下1 500 N的載荷,通過變形量來計(jì)算支承結(jié)構(gòu)的剛度。

    圖3為sloid95結(jié)構(gòu)單元施加載荷后的變形圖,另外兩個結(jié)構(gòu)單元仿真結(jié)果圖與其相似,僅具體的變形量不同,三種結(jié)構(gòu)單元的仿真計(jì)算結(jié)果見表1。

    1.3 陪試軸承端剛度仿真計(jì)算

    陪試軸承端的作用主要為轉(zhuǎn)子系統(tǒng)提供合理的支點(diǎn)位置和支承剛度,為了保證試驗(yàn)端試驗(yàn)數(shù)據(jù)的可分析性,要求陪試軸承端的結(jié)構(gòu)簡單,以減少對試驗(yàn)結(jié)果的影響,主要結(jié)構(gòu)有軸承座、“L”形安裝邊、支承底座。對支承結(jié)構(gòu)進(jìn)行整體網(wǎng)格劃分后如圖4所示。

    圖3 試驗(yàn)軸承端仿真結(jié)果圖

    圖4 陪試軸承端網(wǎng)格劃分

    由于支承底座同樣通過四個地腳螺栓與試驗(yàn)臺固定,因此也將支承底座與試驗(yàn)臺之間的配合面施加[UX=]0,[UY=0,][Z=0]的約束,分別在軸承配合面上施加豎直向上和豎直向下[1 500]N的載荷,通過變形量來計(jì)算支承結(jié)構(gòu)的剛度。

    圖5為sloid187結(jié)構(gòu)單元施加載荷后的變形圖,另外兩個結(jié)構(gòu)單元仿真結(jié)果圖與其相似,僅具體的變形量不同,三種結(jié)構(gòu)單元的仿真計(jì)算結(jié)果見表2。

    圖5 陪試軸承端仿真結(jié)果圖

    2 動力學(xué)特性預(yù)估

    根據(jù)前面預(yù)估的軸承支承剛度值,采用傳遞矩陣法對轉(zhuǎn)子的動力學(xué)特性進(jìn)行預(yù)估,傳遞矩陣法的基本原理是,取不同的轉(zhuǎn)速值,循環(huán)進(jìn)行各軸段截面狀態(tài)參數(shù)的逐段推算,直至滿足轉(zhuǎn)軸另一端的邊界條件[7?8]。本文根據(jù)轉(zhuǎn)軸不同的截面直徑、支承位置及轉(zhuǎn)盤位置,將轉(zhuǎn)軸離散為14個軸段、2個彈性支承站和1個盤站,圖6所示為軸承試驗(yàn)器的計(jì)算模型。

    圖6 轉(zhuǎn)子計(jì)算模型

    將各個單元所計(jì)算剛度的平均值分別代入到計(jì)算模型中的彈性支承站,編程計(jì)算轉(zhuǎn)子臨界轉(zhuǎn)速,計(jì)算結(jié)果見表3。

    3 試驗(yàn)與仿真結(jié)果對比分析

    圖7所示為軸承試驗(yàn)器實(shí)物圖,軸承座上安裝速度傳感器測量軸承振動,盤上安裝位移傳感器測量轉(zhuǎn)子振動,通過測量各個通道的振動幅值來確定轉(zhuǎn)子的一階臨界轉(zhuǎn)速。

    圖8所示為軸承試驗(yàn)器增速過程,圖中顯示的為盤上的振動數(shù)據(jù),實(shí)驗(yàn)測得轉(zhuǎn)子的一階臨界轉(zhuǎn)速為5 222 r/min。

    圖7 軸承試驗(yàn)器實(shí)物圖

    將計(jì)算結(jié)果與實(shí)驗(yàn)測得的臨界轉(zhuǎn)速對比可知:采用solid45單元結(jié)構(gòu)進(jìn)行剛度仿真計(jì)算,得到的臨界轉(zhuǎn)速比實(shí)測值高398 r/min,誤差為7.6%;采用solid95結(jié)構(gòu)單元進(jìn)行剛度仿真計(jì)算,得到的臨界轉(zhuǎn)速比實(shí)測值高311,誤差為5.9%;采用solid187結(jié)構(gòu)單元進(jìn)行剛度仿真計(jì)算,得到的臨界轉(zhuǎn)速比實(shí)測值高168 r/min,誤差為3.2%。

    圖8 試驗(yàn)器增速Bode圖

    4 結(jié) 語

    通過仿真計(jì)算與實(shí)驗(yàn)對比分析可知,對于采用雙“L”安裝邊支承結(jié)構(gòu),利用solid87單元結(jié)構(gòu)進(jìn)行剛度預(yù)估,能夠得到較為準(zhǔn)確的臨界轉(zhuǎn)速計(jì)算結(jié)果,因此在今后的設(shè)計(jì)工作中可以采用solid187結(jié)構(gòu)單元對相似的支承結(jié)構(gòu)進(jìn)行剛度仿真計(jì)算。同時(shí)注意到采用這三種結(jié)構(gòu)單元計(jì)算出的臨界轉(zhuǎn)速全部偏大,若采用此種方法計(jì)算轉(zhuǎn)子的臨界轉(zhuǎn)速,需要在軸承試驗(yàn)器增速過程,注意臨界轉(zhuǎn)速的提前來到。

    參考文獻(xiàn)

    [1] 白中祥,吳偉亮.轉(zhuǎn)子支承系統(tǒng)的剛度對其臨界轉(zhuǎn)速的影響[J].機(jī)電設(shè)備,2012(5):46?49.

    [2] 洪杰,王華,肖大為,等.轉(zhuǎn)子支承動剛度對轉(zhuǎn)子動力特性的影響分析[J].航空發(fā)動機(jī),2008,34(1):23?27.

    [3] 程小勇,陳果,李成剛,等.航空發(fā)動機(jī)轉(zhuǎn)子系統(tǒng)模擬支承設(shè)計(jì)與剛度計(jì)算[J].航空計(jì)算技術(shù),2012,42(6):9?12.

    [4] 盛步云,張濤,丁毓峰,等.支承剛度對汽輪機(jī)轉(zhuǎn)子動力學(xué)特性的影響分析[J].機(jī)械設(shè)計(jì),2008,25(12):38?40.

    [5] 趙文濤,陳果,李瓊,等.航空發(fā)動機(jī)機(jī)匣支承動剛度有限元計(jì)算及驗(yàn)證[J].航空計(jì)算技術(shù),2011,41(5):34?38.

    [6] 何新榮,傅行軍.基于ANSYS的軸承?轉(zhuǎn)子系統(tǒng)動力特性研究[J].機(jī)械研究與應(yīng)用,2011(6):39?41.

    [7] 呂文林.航空發(fā)動機(jī)強(qiáng)度計(jì)算[M].北京:航空工業(yè)出版社,2000.

    [8] 徐龍祥.高速旋轉(zhuǎn)機(jī)械軸系動力學(xué)設(shè)計(jì)[M].北京:國防工業(yè)出版社,1994.

    [9] 馬輝,惠紅杰,唐玉生,等.不同單元計(jì)算轉(zhuǎn)子臨界轉(zhuǎn)速的對比分析[J].振動與沖擊,2012,31(z1):40?43.

    [10] 陳忠.滾動軸承及其支承的剛度計(jì)算[J].煤礦機(jī)械,2006,27(3):387?388.

    [11] 鐘芳明,龔建政,賀星.燃?xì)廨啓C(jī)支承剛度計(jì)算研究[J].船海工程,2012,41(1):41?44.

    [12] 王毅,吳立言,寒冰.ANSYS的兩種有限單元應(yīng)用研究[J].科學(xué)技術(shù)與工程,2007,7(6):955?958.

    1.2 試驗(yàn)軸承端剛度仿真計(jì)算

    試驗(yàn)軸承端用于研究軸承的結(jié)構(gòu)、工藝參數(shù)對軸承支承動力學(xué)特性的影響,其支承結(jié)構(gòu)形式與發(fā)動機(jī)高壓轉(zhuǎn)子前支承一致,因此結(jié)構(gòu)較為復(fù)雜,主要結(jié)構(gòu)有軸承座、安裝環(huán)、“L”形安裝邊、支承底座。對支承結(jié)構(gòu)進(jìn)行整體網(wǎng)格劃分后如圖2所示。

    圖2 試驗(yàn)軸承端網(wǎng)格劃分

    考慮到支承底座通過四個地腳螺栓與試驗(yàn)臺固定,因此在進(jìn)行計(jì)算時(shí),將支承底座與試驗(yàn)臺之間的配合面施加[UX=0,][UY=0,][Z=0]的約束。分別在軸承配合面上施加豎直向上和豎直向下1 500 N的載荷,通過變形量來計(jì)算支承結(jié)構(gòu)的剛度。

    圖3為sloid95結(jié)構(gòu)單元施加載荷后的變形圖,另外兩個結(jié)構(gòu)單元仿真結(jié)果圖與其相似,僅具體的變形量不同,三種結(jié)構(gòu)單元的仿真計(jì)算結(jié)果見表1。

    1.3 陪試軸承端剛度仿真計(jì)算

    陪試軸承端的作用主要為轉(zhuǎn)子系統(tǒng)提供合理的支點(diǎn)位置和支承剛度,為了保證試驗(yàn)端試驗(yàn)數(shù)據(jù)的可分析性,要求陪試軸承端的結(jié)構(gòu)簡單,以減少對試驗(yàn)結(jié)果的影響,主要結(jié)構(gòu)有軸承座、“L”形安裝邊、支承底座。對支承結(jié)構(gòu)進(jìn)行整體網(wǎng)格劃分后如圖4所示。

    圖3 試驗(yàn)軸承端仿真結(jié)果圖

    圖4 陪試軸承端網(wǎng)格劃分

    由于支承底座同樣通過四個地腳螺栓與試驗(yàn)臺固定,因此也將支承底座與試驗(yàn)臺之間的配合面施加[UX=]0,[UY=0,][Z=0]的約束,分別在軸承配合面上施加豎直向上和豎直向下[1 500]N的載荷,通過變形量來計(jì)算支承結(jié)構(gòu)的剛度。

    圖5為sloid187結(jié)構(gòu)單元施加載荷后的變形圖,另外兩個結(jié)構(gòu)單元仿真結(jié)果圖與其相似,僅具體的變形量不同,三種結(jié)構(gòu)單元的仿真計(jì)算結(jié)果見表2。

    圖5 陪試軸承端仿真結(jié)果圖

    2 動力學(xué)特性預(yù)估

    根據(jù)前面預(yù)估的軸承支承剛度值,采用傳遞矩陣法對轉(zhuǎn)子的動力學(xué)特性進(jìn)行預(yù)估,傳遞矩陣法的基本原理是,取不同的轉(zhuǎn)速值,循環(huán)進(jìn)行各軸段截面狀態(tài)參數(shù)的逐段推算,直至滿足轉(zhuǎn)軸另一端的邊界條件[7?8]。本文根據(jù)轉(zhuǎn)軸不同的截面直徑、支承位置及轉(zhuǎn)盤位置,將轉(zhuǎn)軸離散為14個軸段、2個彈性支承站和1個盤站,圖6所示為軸承試驗(yàn)器的計(jì)算模型。

    圖6 轉(zhuǎn)子計(jì)算模型

    將各個單元所計(jì)算剛度的平均值分別代入到計(jì)算模型中的彈性支承站,編程計(jì)算轉(zhuǎn)子臨界轉(zhuǎn)速,計(jì)算結(jié)果見表3。

    3 試驗(yàn)與仿真結(jié)果對比分析

    圖7所示為軸承試驗(yàn)器實(shí)物圖,軸承座上安裝速度傳感器測量軸承振動,盤上安裝位移傳感器測量轉(zhuǎn)子振動,通過測量各個通道的振動幅值來確定轉(zhuǎn)子的一階臨界轉(zhuǎn)速。

    圖8所示為軸承試驗(yàn)器增速過程,圖中顯示的為盤上的振動數(shù)據(jù),實(shí)驗(yàn)測得轉(zhuǎn)子的一階臨界轉(zhuǎn)速為5 222 r/min。

    圖7 軸承試驗(yàn)器實(shí)物圖

    將計(jì)算結(jié)果與實(shí)驗(yàn)測得的臨界轉(zhuǎn)速對比可知:采用solid45單元結(jié)構(gòu)進(jìn)行剛度仿真計(jì)算,得到的臨界轉(zhuǎn)速比實(shí)測值高398 r/min,誤差為7.6%;采用solid95結(jié)構(gòu)單元進(jìn)行剛度仿真計(jì)算,得到的臨界轉(zhuǎn)速比實(shí)測值高311,誤差為5.9%;采用solid187結(jié)構(gòu)單元進(jìn)行剛度仿真計(jì)算,得到的臨界轉(zhuǎn)速比實(shí)測值高168 r/min,誤差為3.2%。

    圖8 試驗(yàn)器增速Bode圖

    4 結(jié) 語

    通過仿真計(jì)算與實(shí)驗(yàn)對比分析可知,對于采用雙“L”安裝邊支承結(jié)構(gòu),利用solid87單元結(jié)構(gòu)進(jìn)行剛度預(yù)估,能夠得到較為準(zhǔn)確的臨界轉(zhuǎn)速計(jì)算結(jié)果,因此在今后的設(shè)計(jì)工作中可以采用solid187結(jié)構(gòu)單元對相似的支承結(jié)構(gòu)進(jìn)行剛度仿真計(jì)算。同時(shí)注意到采用這三種結(jié)構(gòu)單元計(jì)算出的臨界轉(zhuǎn)速全部偏大,若采用此種方法計(jì)算轉(zhuǎn)子的臨界轉(zhuǎn)速,需要在軸承試驗(yàn)器增速過程,注意臨界轉(zhuǎn)速的提前來到。

    參考文獻(xiàn)

    [1] 白中祥,吳偉亮.轉(zhuǎn)子支承系統(tǒng)的剛度對其臨界轉(zhuǎn)速的影響[J].機(jī)電設(shè)備,2012(5):46?49.

    [2] 洪杰,王華,肖大為,等.轉(zhuǎn)子支承動剛度對轉(zhuǎn)子動力特性的影響分析[J].航空發(fā)動機(jī),2008,34(1):23?27.

    [3] 程小勇,陳果,李成剛,等.航空發(fā)動機(jī)轉(zhuǎn)子系統(tǒng)模擬支承設(shè)計(jì)與剛度計(jì)算[J].航空計(jì)算技術(shù),2012,42(6):9?12.

    [4] 盛步云,張濤,丁毓峰,等.支承剛度對汽輪機(jī)轉(zhuǎn)子動力學(xué)特性的影響分析[J].機(jī)械設(shè)計(jì),2008,25(12):38?40.

    [5] 趙文濤,陳果,李瓊,等.航空發(fā)動機(jī)機(jī)匣支承動剛度有限元計(jì)算及驗(yàn)證[J].航空計(jì)算技術(shù),2011,41(5):34?38.

    [6] 何新榮,傅行軍.基于ANSYS的軸承?轉(zhuǎn)子系統(tǒng)動力特性研究[J].機(jī)械研究與應(yīng)用,2011(6):39?41.

    [7] 呂文林.航空發(fā)動機(jī)強(qiáng)度計(jì)算[M].北京:航空工業(yè)出版社,2000.

    [8] 徐龍祥.高速旋轉(zhuǎn)機(jī)械軸系動力學(xué)設(shè)計(jì)[M].北京:國防工業(yè)出版社,1994.

    [9] 馬輝,惠紅杰,唐玉生,等.不同單元計(jì)算轉(zhuǎn)子臨界轉(zhuǎn)速的對比分析[J].振動與沖擊,2012,31(z1):40?43.

    [10] 陳忠.滾動軸承及其支承的剛度計(jì)算[J].煤礦機(jī)械,2006,27(3):387?388.

    [11] 鐘芳明,龔建政,賀星.燃?xì)廨啓C(jī)支承剛度計(jì)算研究[J].船海工程,2012,41(1):41?44.

    [12] 王毅,吳立言,寒冰.ANSYS的兩種有限單元應(yīng)用研究[J].科學(xué)技術(shù)與工程,2007,7(6):955?958.

    猜你喜歡
    結(jié)構(gòu)單元航空發(fā)動機(jī)
    基于結(jié)構(gòu)單元的機(jī)床動態(tài)特性分析??
    鈣鈦礦結(jié)構(gòu)類型的功能材料的結(jié)構(gòu)單元和結(jié)構(gòu)演變
    某型發(fā)動機(jī)喘振故障的研究
    航空發(fā)動機(jī)小零件的標(biāo)準(zhǔn)化實(shí)踐與數(shù)字化管理
    航空發(fā)動機(jī)管理信息系統(tǒng)研究
    淺談航空發(fā)動機(jī)起動與點(diǎn)火系統(tǒng)
    太赫茲波段雙寬帶極化無關(guān)超材料吸波體的設(shè)計(jì)研究
    一種具有表面活性功能的聚合物及其制備方法和應(yīng)用
    兩個基于二噻吩乙烯結(jié)構(gòu)單元雙核釕乙烯配合物的合成,表征和性質(zhì)
    永安鎮(zhèn)油田永3斷塊沙二下河口壩儲層結(jié)構(gòu)單元劃分及其意義
    久99久视频精品免费| 91大片在线观看| 777米奇影视久久| 狂野欧美激情性xxxx| 欧美一级毛片孕妇| 国产成人啪精品午夜网站| av一本久久久久| 9色porny在线观看| 精品少妇久久久久久888优播| 91老司机精品| 国产高清激情床上av| 亚洲片人在线观看| 午夜福利一区二区在线看| avwww免费| 欧美黄色片欧美黄色片| 捣出白浆h1v1| 极品教师在线免费播放| 亚洲精品中文字幕一二三四区| 久久亚洲真实| 国产色视频综合| 亚洲人成伊人成综合网2020| 久久精品国产99精品国产亚洲性色 | 国产人伦9x9x在线观看| 成人影院久久| 99精品在免费线老司机午夜| 黑人操中国人逼视频| 午夜福利影视在线免费观看| 老司机在亚洲福利影院| 在线观看www视频免费| 午夜视频精品福利| 精品电影一区二区在线| 嫩草影视91久久| 国产成人系列免费观看| 久久中文看片网| 日本精品一区二区三区蜜桃| 午夜福利免费观看在线| 天堂中文最新版在线下载| 人成视频在线观看免费观看| 亚洲精品自拍成人| 久久香蕉激情| 亚洲成国产人片在线观看| 黄色丝袜av网址大全| 亚洲 欧美一区二区三区| 国产亚洲精品一区二区www | 亚洲av欧美aⅴ国产| 在线观看免费日韩欧美大片| 午夜精品久久久久久毛片777| 三级毛片av免费| 成人三级做爰电影| 国产精品久久电影中文字幕 | 老司机亚洲免费影院| 色老头精品视频在线观看| 热re99久久精品国产66热6| 18禁国产床啪视频网站| 中文字幕人妻熟女乱码| 亚洲成人手机| 天天影视国产精品| 国产精品 国内视频| 免费av中文字幕在线| 亚洲色图av天堂| 一区福利在线观看| 1024视频免费在线观看| 亚洲欧洲精品一区二区精品久久久| 777久久人妻少妇嫩草av网站| 无遮挡黄片免费观看| 黄色女人牲交| 久久久水蜜桃国产精品网| 国产成人精品在线电影| 午夜亚洲福利在线播放| 免费观看人在逋| 午夜亚洲福利在线播放| 国产熟女午夜一区二区三区| 亚洲av片天天在线观看| 国产一区有黄有色的免费视频| 亚洲九九香蕉| 欧美中文综合在线视频| 久久国产精品影院| 天堂√8在线中文| 麻豆国产av国片精品| 人妻丰满熟妇av一区二区三区 | 久久亚洲真实| 伦理电影免费视频| 久久婷婷成人综合色麻豆| 天堂√8在线中文| 国产欧美日韩综合在线一区二区| 好看av亚洲va欧美ⅴa在| 国产亚洲精品第一综合不卡| 18禁美女被吸乳视频| 欧美一级毛片孕妇| 水蜜桃什么品种好| www.999成人在线观看| 欧美色视频一区免费| 一级毛片女人18水好多| 精品亚洲成国产av| 日韩人妻精品一区2区三区| 色尼玛亚洲综合影院| 少妇 在线观看| 久久天堂一区二区三区四区| 美女 人体艺术 gogo| 91在线观看av| 国产精品香港三级国产av潘金莲| 欧美亚洲 丝袜 人妻 在线| 两个人看的免费小视频| av线在线观看网站| 国产成人精品久久二区二区91| 美女扒开内裤让男人捅视频| 热re99久久精品国产66热6| 国产欧美亚洲国产| 亚洲成国产人片在线观看| 美女国产高潮福利片在线看| 十八禁网站免费在线| 欧美日韩福利视频一区二区| 亚洲欧美一区二区三区久久| 精品久久久久久,| av网站免费在线观看视频| 久久人妻av系列| 欧美亚洲日本最大视频资源| 叶爱在线成人免费视频播放| 欧美一级毛片孕妇| 精品一区二区三区视频在线观看免费 | 亚洲国产毛片av蜜桃av| 久久精品成人免费网站| 久久人人爽av亚洲精品天堂| 12—13女人毛片做爰片一| 成人免费观看视频高清| 婷婷丁香在线五月| 妹子高潮喷水视频| 日韩欧美三级三区| 亚洲人成电影免费在线| 女人精品久久久久毛片| 欧美精品av麻豆av| 国产成人免费无遮挡视频| 欧美日韩精品网址| 亚洲成人免费电影在线观看| 国产亚洲精品第一综合不卡| 捣出白浆h1v1| 18在线观看网站| 可以免费在线观看a视频的电影网站| 人人妻人人添人人爽欧美一区卜| 国产精品国产高清国产av | 后天国语完整版免费观看| 美女扒开内裤让男人捅视频| 亚洲欧美一区二区三区久久| 国产高清视频在线播放一区| 少妇被粗大的猛进出69影院| 午夜福利免费观看在线| 亚洲熟女精品中文字幕| 多毛熟女@视频| 热99国产精品久久久久久7| ponron亚洲| 亚洲精品一二三| 中亚洲国语对白在线视频| x7x7x7水蜜桃| av国产精品久久久久影院| 777久久人妻少妇嫩草av网站| 亚洲美女黄片视频| 国产精品九九99| 国产成人影院久久av| 婷婷丁香在线五月| 好看av亚洲va欧美ⅴa在| 好男人电影高清在线观看| 精品乱码久久久久久99久播| 天天操日日干夜夜撸| 成人国语在线视频| 两个人看的免费小视频| av电影中文网址| 在线观看66精品国产| 国产欧美日韩一区二区三区在线| 亚洲精品美女久久久久99蜜臀| 成人黄色视频免费在线看| 黑人操中国人逼视频| 午夜福利乱码中文字幕| 男人操女人黄网站| 又黄又爽又免费观看的视频| 大香蕉久久网| 大码成人一级视频| 飞空精品影院首页| 亚洲av第一区精品v没综合| av片东京热男人的天堂| 高潮久久久久久久久久久不卡| 少妇 在线观看| 新久久久久国产一级毛片| 久久久久精品人妻al黑| 日本a在线网址| 国产成人欧美在线观看 | 麻豆国产av国片精品| 中文字幕精品免费在线观看视频| 亚洲熟妇中文字幕五十中出 | 香蕉久久夜色| 天堂动漫精品| 国产精品欧美亚洲77777| 成人av一区二区三区在线看| 母亲3免费完整高清在线观看| 黄色怎么调成土黄色| 国产日韩欧美亚洲二区| 欧美成狂野欧美在线观看| 日韩免费高清中文字幕av| 热99久久久久精品小说推荐| 三级毛片av免费| 九色亚洲精品在线播放| 99精国产麻豆久久婷婷| 精品一区二区三区av网在线观看| 国产真人三级小视频在线观看| 亚洲五月婷婷丁香| 欧美日韩瑟瑟在线播放| 亚洲av日韩精品久久久久久密| 国产精品免费大片| av中文乱码字幕在线| 日韩人妻精品一区2区三区| 真人做人爱边吃奶动态| 18禁国产床啪视频网站| 少妇的丰满在线观看| 欧美国产精品va在线观看不卡| 99国产极品粉嫩在线观看| 91精品三级在线观看| 成年人黄色毛片网站| 国产成人av激情在线播放| 亚洲中文字幕日韩| 欧美午夜高清在线| 十八禁人妻一区二区| 国产精品永久免费网站| 欧美日韩av久久| 欧美黑人精品巨大| 大型av网站在线播放| 熟女少妇亚洲综合色aaa.| www日本在线高清视频| 亚洲第一av免费看| 久久精品国产亚洲av高清一级| 大码成人一级视频| 曰老女人黄片| 80岁老熟妇乱子伦牲交| 亚洲黑人精品在线| www.精华液| 久久国产精品影院| 欧美激情 高清一区二区三区| 超色免费av| 9色porny在线观看| 老司机靠b影院| 日韩三级视频一区二区三区| 精品熟女少妇八av免费久了| 一级毛片高清免费大全| 人人妻,人人澡人人爽秒播| 国产一区二区三区在线臀色熟女 | av片东京热男人的天堂| 国产免费男女视频| 国产精品久久久av美女十八| 精品久久久久久久毛片微露脸| 成人国语在线视频| 成年人午夜在线观看视频| 亚洲片人在线观看| a在线观看视频网站| 久久热在线av| 久久久久久人人人人人| 国产男女超爽视频在线观看| 首页视频小说图片口味搜索| 久久精品国产综合久久久| 十八禁高潮呻吟视频| 精品国产乱码久久久久久男人| 在线av久久热| 丝袜美足系列| 伊人久久大香线蕉亚洲五| 久久香蕉国产精品| 日本撒尿小便嘘嘘汇集6| 脱女人内裤的视频| 成人特级黄色片久久久久久久| 宅男免费午夜| 在线观看免费日韩欧美大片| 美女国产高潮福利片在线看| 五月开心婷婷网| 欧美中文综合在线视频| 国产xxxxx性猛交| 久久这里只有精品19| 亚洲成人免费电影在线观看| 亚洲精品自拍成人| 欧美精品亚洲一区二区| 国产有黄有色有爽视频| 在线观看舔阴道视频| 男女床上黄色一级片免费看| 国产日韩一区二区三区精品不卡| 一进一出好大好爽视频| 女人被狂操c到高潮| 一边摸一边抽搐一进一出视频| 色在线成人网| 国产成人精品在线电影| 在线观看日韩欧美| 久久人人爽av亚洲精品天堂| 一边摸一边做爽爽视频免费| 日韩熟女老妇一区二区性免费视频| 久9热在线精品视频| 99国产精品99久久久久| 亚洲 国产 在线| 女性生殖器流出的白浆| 女人精品久久久久毛片| 亚洲欧美一区二区三区久久| 久久精品熟女亚洲av麻豆精品| 美女 人体艺术 gogo| 亚洲精品中文字幕在线视频| 成年人黄色毛片网站| 国产精品免费视频内射| 亚洲五月天丁香| 韩国av一区二区三区四区| 午夜福利欧美成人| 高清毛片免费观看视频网站 | 男人操女人黄网站| 午夜91福利影院| 亚洲国产看品久久| 久久精品亚洲av国产电影网| 美女福利国产在线| 欧美精品高潮呻吟av久久| 在线观看午夜福利视频| 啦啦啦视频在线资源免费观看| tube8黄色片| 在线观看午夜福利视频| 久久 成人 亚洲| 久久久国产欧美日韩av| 国产精品久久视频播放| 精品国产乱子伦一区二区三区| 丝瓜视频免费看黄片| 亚洲九九香蕉| av一本久久久久| 欧美+亚洲+日韩+国产| 国产伦人伦偷精品视频| av天堂久久9| 亚洲成人免费av在线播放| 国产成人影院久久av| 国精品久久久久久国模美| 男女午夜视频在线观看| 夜夜爽天天搞| 脱女人内裤的视频| 免费少妇av软件| 久久午夜亚洲精品久久| 亚洲午夜理论影院| 日韩精品免费视频一区二区三区| 男人的好看免费观看在线视频 | 久久久久精品国产欧美久久久| 在线观看www视频免费| 大陆偷拍与自拍| 欧洲精品卡2卡3卡4卡5卡区| 嫩草影视91久久| 男女免费视频国产| 亚洲欧美日韩高清在线视频| 91成人精品电影| 欧美亚洲 丝袜 人妻 在线| 亚洲av片天天在线观看| 欧美黑人精品巨大| 老司机午夜十八禁免费视频| 国产精品乱码一区二三区的特点 | 9热在线视频观看99| 天堂√8在线中文| 91麻豆av在线| 欧美 日韩 精品 国产| а√天堂www在线а√下载 | www.999成人在线观看| 女人被狂操c到高潮| 久久久精品区二区三区| 欧美激情高清一区二区三区| 正在播放国产对白刺激| 好男人电影高清在线观看| 电影成人av| 黄色毛片三级朝国网站| 日本欧美视频一区| 在线观看日韩欧美| 国产精品99久久99久久久不卡| 他把我摸到了高潮在线观看| 女警被强在线播放| 一边摸一边抽搐一进一出视频| 99在线人妻在线中文字幕 | 国产xxxxx性猛交| 一级作爱视频免费观看| 欧美日韩精品网址| 午夜福利免费观看在线| 正在播放国产对白刺激| 男女午夜视频在线观看| 亚洲五月婷婷丁香| 又黄又粗又硬又大视频| 丰满的人妻完整版| 日本wwww免费看| 亚洲欧美日韩另类电影网站| 人人妻人人澡人人爽人人夜夜| 午夜久久久在线观看| 黑人巨大精品欧美一区二区蜜桃| 亚洲精华国产精华精| 男男h啪啪无遮挡| 一进一出好大好爽视频| 国产单亲对白刺激| 国产亚洲av高清不卡| 露出奶头的视频| 国产成人啪精品午夜网站| 三上悠亚av全集在线观看| 久久久久久久久久久久大奶| 91精品国产国语对白视频| 一进一出抽搐gif免费好疼 | 亚洲中文日韩欧美视频| 国产精品自产拍在线观看55亚洲 | 国产精品久久视频播放| 热re99久久精品国产66热6| 波多野结衣一区麻豆| 国产精品综合久久久久久久免费 | 人人妻人人澡人人爽人人夜夜| 亚洲av成人一区二区三| 亚洲av日韩在线播放| 9热在线视频观看99| 久久天躁狠狠躁夜夜2o2o| 一边摸一边做爽爽视频免费| 久久香蕉激情| 女人高潮潮喷娇喘18禁视频| 国产成人一区二区三区免费视频网站| 亚洲欧美日韩另类电影网站| 91av网站免费观看| 精品少妇久久久久久888优播| 一区二区日韩欧美中文字幕| 久久ye,这里只有精品| 亚洲av美国av| 热re99久久国产66热| 免费高清在线观看日韩| 高清欧美精品videossex| 国产成人精品无人区| 欧美日韩精品网址| 香蕉国产在线看| 午夜精品国产一区二区电影| 日本撒尿小便嘘嘘汇集6| 精品久久久久久,| 国产男女内射视频| 热99久久久久精品小说推荐| 欧美最黄视频在线播放免费 | 久久99一区二区三区| 欧美日韩精品网址| 午夜精品在线福利| 嫩草影视91久久| 久久久国产成人精品二区 | 亚洲黑人精品在线| 建设人人有责人人尽责人人享有的| 一区二区三区精品91| 婷婷成人精品国产| 极品人妻少妇av视频| 亚洲国产毛片av蜜桃av| 正在播放国产对白刺激| 少妇被粗大的猛进出69影院| 精品午夜福利视频在线观看一区| 十八禁高潮呻吟视频| av福利片在线| 精品国产亚洲在线| 交换朋友夫妻互换小说| 午夜91福利影院| 女人久久www免费人成看片| 欧美另类亚洲清纯唯美| 国产男靠女视频免费网站| av福利片在线| 久久久精品区二区三区| 中亚洲国语对白在线视频| 欧美精品啪啪一区二区三区| 成人亚洲精品一区在线观看| 夜夜躁狠狠躁天天躁| 人妻 亚洲 视频| 国产精品亚洲av一区麻豆| 纯流量卡能插随身wifi吗| 两性夫妻黄色片| 成人特级黄色片久久久久久久| 亚洲精品自拍成人| 国产精品1区2区在线观看. | 亚洲一区二区三区欧美精品| 国产激情欧美一区二区| 91麻豆精品激情在线观看国产 | 中文字幕人妻熟女乱码| 国产精品国产高清国产av | 真人做人爱边吃奶动态| 91精品三级在线观看| 看片在线看免费视频| 国产伦人伦偷精品视频| 亚洲欧美激情在线| 男人操女人黄网站| 久久久国产欧美日韩av| 中文字幕最新亚洲高清| 91成年电影在线观看| 又黄又爽又免费观看的视频| 午夜精品在线福利| 色在线成人网| 国产aⅴ精品一区二区三区波| 国产一区二区激情短视频| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲五月婷婷丁香| 午夜免费成人在线视频| 午夜久久久在线观看| 亚洲熟妇中文字幕五十中出 | 亚洲成人免费av在线播放| 午夜成年电影在线免费观看| 91精品三级在线观看| 色播在线永久视频| 亚洲五月天丁香| 老熟妇仑乱视频hdxx| 女人精品久久久久毛片| 国产一卡二卡三卡精品| 在线观看免费高清a一片| 天堂中文最新版在线下载| 亚洲成人免费电影在线观看| 性少妇av在线| 亚洲视频免费观看视频| 少妇猛男粗大的猛烈进出视频| 国产伦人伦偷精品视频| 热re99久久精品国产66热6| av片东京热男人的天堂| 久久精品国产亚洲av香蕉五月 | 午夜亚洲福利在线播放| 精品国产亚洲在线| 首页视频小说图片口味搜索| 久久天堂一区二区三区四区| 啦啦啦视频在线资源免费观看| 高清黄色对白视频在线免费看| 亚洲中文字幕日韩| 久久久久久亚洲精品国产蜜桃av| 久久中文看片网| 黑人操中国人逼视频| 久久精品91无色码中文字幕| 亚洲av成人不卡在线观看播放网| 国产精品久久电影中文字幕 | 大片电影免费在线观看免费| 国产欧美日韩精品亚洲av| 国产亚洲精品一区二区www | 国产男靠女视频免费网站| 欧洲精品卡2卡3卡4卡5卡区| 成人免费观看视频高清| 亚洲精品一卡2卡三卡4卡5卡| 日本一区二区免费在线视频| 日韩欧美一区视频在线观看| 精品高清国产在线一区| 中国美女看黄片| 自线自在国产av| √禁漫天堂资源中文www| 国产精品二区激情视频| 午夜91福利影院| 丝袜人妻中文字幕| 久久精品国产a三级三级三级| 亚洲国产中文字幕在线视频| 精品人妻熟女毛片av久久网站| 亚洲va日本ⅴa欧美va伊人久久| 如日韩欧美国产精品一区二区三区| 在线视频色国产色| 久久久久久久久免费视频了| 香蕉丝袜av| 国产无遮挡羞羞视频在线观看| 中文字幕色久视频| 久久久水蜜桃国产精品网| 99精品在免费线老司机午夜| 亚洲精品国产精品久久久不卡| 高潮久久久久久久久久久不卡| 久久香蕉国产精品| 黄色怎么调成土黄色| 亚洲,欧美精品.| 国产精品一区二区免费欧美| 欧美黑人欧美精品刺激| 国产日韩一区二区三区精品不卡| 天堂动漫精品| 99国产精品99久久久久| 久久人妻福利社区极品人妻图片| 黄网站色视频无遮挡免费观看| 国产精品亚洲av一区麻豆| 精品乱码久久久久久99久播| 超碰成人久久| 波多野结衣一区麻豆| 久久狼人影院| 中文字幕av电影在线播放| cao死你这个sao货| 伊人久久大香线蕉亚洲五| 正在播放国产对白刺激| 欧美激情高清一区二区三区| 黄色 视频免费看| 国产精品免费大片| 免费高清在线观看日韩| 动漫黄色视频在线观看| 狂野欧美激情性xxxx| 国产精品免费一区二区三区在线 | 成人av一区二区三区在线看| 国产成人av教育| 女人高潮潮喷娇喘18禁视频| 国产亚洲欧美精品永久| 午夜免费观看网址| 免费观看精品视频网站| 国产成人精品在线电影| 女性被躁到高潮视频| 午夜成年电影在线免费观看| 国产欧美日韩一区二区三| 怎么达到女性高潮| 老司机影院毛片| 欧美日韩瑟瑟在线播放| 又大又爽又粗| 午夜视频精品福利| 国产亚洲欧美98| 一区福利在线观看| 两个人看的免费小视频| 欧美日韩中文字幕国产精品一区二区三区 | 成人av一区二区三区在线看| 亚洲欧美激情综合另类| 国产不卡一卡二| 欧美精品av麻豆av| 亚洲精品国产精品久久久不卡| 女人爽到高潮嗷嗷叫在线视频| 18禁美女被吸乳视频| 女同久久另类99精品国产91| 别揉我奶头~嗯~啊~动态视频| 免费黄频网站在线观看国产| 欧美国产精品一级二级三级| 欧美老熟妇乱子伦牲交| 久久国产亚洲av麻豆专区| 国产精品av久久久久免费| 97人妻天天添夜夜摸| 桃红色精品国产亚洲av| 老熟妇仑乱视频hdxx| 亚洲第一欧美日韩一区二区三区| 国产精品免费一区二区三区在线 | 国产又爽黄色视频| xxx96com| 极品教师在线免费播放| 国产精品1区2区在线观看. | 欧美乱码精品一区二区三区| 亚洲国产欧美一区二区综合|