• <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)單元劃分及其意義
    亚洲第一av免费看| 久久精品久久久久久噜噜老黄| 国产亚洲午夜精品一区二区久久| 亚洲精品日本国产第一区| 久久99热6这里只有精品| 精品一品国产午夜福利视频| 久久99一区二区三区| 亚洲av男天堂| 成人无遮挡网站| 18禁动态无遮挡网站| 国产在线男女| 亚洲精品aⅴ在线观看| 亚洲av男天堂| 好男人视频免费观看在线| 热re99久久精品国产66热6| 中文字幕亚洲精品专区| 三上悠亚av全集在线观看 | 在线观看免费视频网站a站| 99国产精品免费福利视频| 内射极品少妇av片p| 国产探花极品一区二区| 国产精品一区二区在线观看99| 看十八女毛片水多多多| 久久久久久久久久久免费av| 久久国产精品大桥未久av | 国产成人免费无遮挡视频| 久久久久人妻精品一区果冻| av不卡在线播放| 丝袜在线中文字幕| 精品亚洲乱码少妇综合久久| 国产精品久久久久久久久免| a 毛片基地| 97精品久久久久久久久久精品| 水蜜桃什么品种好| 99热这里只有是精品在线观看| 国产精品国产三级国产专区5o| 最近的中文字幕免费完整| 久久久久网色| 最后的刺客免费高清国语| 色婷婷久久久亚洲欧美| 国产视频内射| 久久99热6这里只有精品| 内地一区二区视频在线| 久久精品国产亚洲网站| 晚上一个人看的免费电影| 国产69精品久久久久777片| 纯流量卡能插随身wifi吗| 人妻一区二区av| 午夜日本视频在线| 亚洲国产精品一区三区| 狠狠精品人妻久久久久久综合| 七月丁香在线播放| 国产男女超爽视频在线观看| videos熟女内射| 在线观看免费视频网站a站| 午夜免费男女啪啪视频观看| 岛国毛片在线播放| 亚洲精品aⅴ在线观看| 日本vs欧美在线观看视频 | 乱码一卡2卡4卡精品| 少妇被粗大的猛进出69影院 | 国产免费一区二区三区四区乱码| 国产片特级美女逼逼视频| 黄色毛片三级朝国网站 | 成年人免费黄色播放视频 | 国产 一区精品| 国产 一区精品| 少妇人妻 视频| 十分钟在线观看高清视频www | 亚洲av电影在线观看一区二区三区| 夫妻午夜视频| 国产美女午夜福利| a级片在线免费高清观看视频| 在线观看av片永久免费下载| 亚洲欧洲精品一区二区精品久久久 | 国产高清国产精品国产三级| 欧美成人精品欧美一级黄| 欧美人与善性xxx| 国产真实伦视频高清在线观看| 日韩三级伦理在线观看| 日韩一区二区视频免费看| 免费人成在线观看视频色| 欧美日本中文国产一区发布| 男人舔奶头视频| 色婷婷av一区二区三区视频| 日本色播在线视频| 国产黄频视频在线观看| 成人18禁高潮啪啪吃奶动态图 | 777米奇影视久久| 久久久久视频综合| 日韩强制内射视频| 久久午夜福利片| 亚洲国产色片| 午夜福利视频精品| 国产视频内射| 中文在线观看免费www的网站| 久久这里有精品视频免费| www.av在线官网国产| 插逼视频在线观看| 国产精品无大码| 一级毛片 在线播放| 丰满乱子伦码专区| 99九九在线精品视频 | 精品久久国产蜜桃| 97精品久久久久久久久久精品| 亚洲无线观看免费| 日韩不卡一区二区三区视频在线| 人妻少妇偷人精品九色| 欧美老熟妇乱子伦牲交| 99热网站在线观看| 欧美xxxx性猛交bbbb| 亚洲欧美成人精品一区二区| 亚洲美女视频黄频| 美女cb高潮喷水在线观看| 成年女人在线观看亚洲视频| av在线观看视频网站免费| av在线播放精品| 亚洲不卡免费看| 美女脱内裤让男人舔精品视频| 一本—道久久a久久精品蜜桃钙片| 在现免费观看毛片| 伊人久久精品亚洲午夜| xxx大片免费视频| 97精品久久久久久久久久精品| 麻豆成人av视频| 看非洲黑人一级黄片| 久久韩国三级中文字幕| 久久久久久久亚洲中文字幕| 国产精品福利在线免费观看| 国精品久久久久久国模美| av在线播放精品| 日产精品乱码卡一卡2卡三| 亚洲精品日韩av片在线观看| 少妇的逼好多水| 91精品国产国语对白视频| 亚洲国产欧美日韩在线播放 | 在线观看三级黄色| 精品一区二区三区视频在线| 国产伦精品一区二区三区四那| 午夜日本视频在线| 少妇裸体淫交视频免费看高清| 免费看日本二区| 国产欧美亚洲国产| 亚洲欧洲国产日韩| 精品卡一卡二卡四卡免费| 老女人水多毛片| 我的老师免费观看完整版| 激情五月婷婷亚洲| 久久99热6这里只有精品| 日韩,欧美,国产一区二区三区| 自拍偷自拍亚洲精品老妇| 国产男人的电影天堂91| av一本久久久久| 日日爽夜夜爽网站| 一区二区三区乱码不卡18| 久久99热6这里只有精品| 成人午夜精彩视频在线观看| 亚洲国产av新网站| 国产日韩欧美视频二区| 亚洲第一区二区三区不卡| av一本久久久久| 中文字幕制服av| 在线天堂最新版资源| 夫妻午夜视频| av在线app专区| 国产永久视频网站| 国精品久久久久久国模美| 国产欧美日韩一区二区三区在线 | 菩萨蛮人人尽说江南好唐韦庄| 国内精品宾馆在线| 亚洲中文av在线| 国产一区亚洲一区在线观看| 建设人人有责人人尽责人人享有的| 性高湖久久久久久久久免费观看| 成人特级av手机在线观看| 香蕉精品网在线| 精品熟女少妇av免费看| 啦啦啦中文免费视频观看日本| videos熟女内射| 国产极品天堂在线| 国产精品久久久久成人av| 亚洲精华国产精华液的使用体验| 日韩三级伦理在线观看| 欧美丝袜亚洲另类| 午夜视频国产福利| 亚洲av成人精品一区久久| 亚洲精品成人av观看孕妇| 美女中出高潮动态图| 国产永久视频网站| 18禁在线播放成人免费| 国产精品久久久久久久久免| a级一级毛片免费在线观看| 三级国产精品欧美在线观看| 大码成人一级视频| 久久久久国产精品人妻一区二区| 赤兔流量卡办理| av免费观看日本| 久久 成人 亚洲| 少妇精品久久久久久久| 在线看a的网站| av免费在线看不卡| 中文字幕av电影在线播放| 免费黄色在线免费观看| 91成人精品电影| 久久鲁丝午夜福利片| 只有这里有精品99| √禁漫天堂资源中文www| 99久久精品国产国产毛片| 亚洲欧美成人精品一区二区| 777米奇影视久久| 嘟嘟电影网在线观看| 亚洲欧美精品自产自拍| 久久 成人 亚洲| 午夜免费观看性视频| 亚洲欧洲国产日韩| 国产精品久久久久久精品古装| 九九在线视频观看精品| 欧美性感艳星| 久久精品国产鲁丝片午夜精品| 在线观看av片永久免费下载| 插阴视频在线观看视频| 欧美激情国产日韩精品一区| 啦啦啦在线观看免费高清www| 成人二区视频| 高清黄色对白视频在线免费看 | 少妇被粗大的猛进出69影院 | 五月开心婷婷网| 国产男女内射视频| 97在线人人人人妻| 国产熟女午夜一区二区三区 | 少妇被粗大的猛进出69影院 | 国产有黄有色有爽视频| 国产高清有码在线观看视频| 午夜免费男女啪啪视频观看| av免费在线看不卡| 亚洲精品乱码久久久v下载方式| 亚洲av国产av综合av卡| 免费观看的影片在线观看| 麻豆成人av视频| 2018国产大陆天天弄谢| 丝袜喷水一区| 99久久精品热视频| 亚洲精品中文字幕在线视频 | av卡一久久| 亚洲欧美成人精品一区二区| √禁漫天堂资源中文www| 十八禁高潮呻吟视频 | 丝瓜视频免费看黄片| 国产精品国产三级专区第一集| 最近手机中文字幕大全| 夫妻性生交免费视频一级片| 偷拍熟女少妇极品色| 欧美97在线视频| 成人国产麻豆网| 国产日韩欧美视频二区| 国产精品一区二区性色av| 一区二区三区四区激情视频| 久久久久久久久久成人| 91久久精品国产一区二区三区| 极品少妇高潮喷水抽搐| 亚洲美女视频黄频| 中文字幕制服av| 久久久久久久久久久久大奶| 久久久久网色| 制服丝袜香蕉在线| videossex国产| 一个人看视频在线观看www免费| 亚洲精品色激情综合| 中文字幕亚洲精品专区| 中文欧美无线码| 人妻少妇偷人精品九色| 国产精品一区二区性色av| 91精品伊人久久大香线蕉| 国产在线男女| 精品久久国产蜜桃| 精品视频人人做人人爽| 国产欧美另类精品又又久久亚洲欧美| 涩涩av久久男人的天堂| 精品久久久久久久久亚洲| 狂野欧美激情性bbbbbb| 伊人亚洲综合成人网| 午夜福利影视在线免费观看| av有码第一页| 国产综合精华液| 久久精品熟女亚洲av麻豆精品| 蜜臀久久99精品久久宅男| 一本—道久久a久久精品蜜桃钙片| 国产极品天堂在线| 18禁裸乳无遮挡动漫免费视频| 综合色丁香网| 国产亚洲一区二区精品| 免费大片黄手机在线观看| 日本av免费视频播放| 免费大片18禁| 香蕉精品网在线| 黄色一级大片看看| 我的女老师完整版在线观看| 亚洲精品第二区| 十八禁网站网址无遮挡 | av女优亚洲男人天堂| 日韩欧美一区视频在线观看 | 啦啦啦啦在线视频资源| 尾随美女入室| 国产熟女午夜一区二区三区 | 亚洲四区av| 国产色爽女视频免费观看| 在线观看av片永久免费下载| 综合色丁香网| 国产又色又爽无遮挡免| 国产爽快片一区二区三区| 91午夜精品亚洲一区二区三区| 亚洲精品国产av成人精品| 久久久国产精品麻豆| 视频中文字幕在线观看| 亚洲,欧美,日韩| av天堂久久9| a级一级毛片免费在线观看| av又黄又爽大尺度在线免费看| 内地一区二区视频在线| 免费播放大片免费观看视频在线观看| 中文字幕免费在线视频6| 91aial.com中文字幕在线观看| 欧美精品亚洲一区二区| av网站免费在线观看视频| 热99国产精品久久久久久7| 国产男女超爽视频在线观看| 成人无遮挡网站| 国产黄片美女视频| 一区二区三区免费毛片| 老司机影院毛片| 少妇猛男粗大的猛烈进出视频| 不卡视频在线观看欧美| 亚洲精品久久午夜乱码| 久久鲁丝午夜福利片| 观看美女的网站| 中文资源天堂在线| 偷拍熟女少妇极品色| 日韩在线高清观看一区二区三区| 成人漫画全彩无遮挡| 国产在线男女| 我要看日韩黄色一级片| 国产毛片在线视频| 久久久久久久精品精品| 久久久久国产网址| 久久 成人 亚洲| 欧美激情极品国产一区二区三区 | 狂野欧美激情性xxxx在线观看| 日韩一区二区视频免费看| 大片电影免费在线观看免费| av专区在线播放| 一边亲一边摸免费视频| a级一级毛片免费在线观看| 一边亲一边摸免费视频| av视频免费观看在线观看| 欧美精品国产亚洲| 秋霞伦理黄片| 久久久久人妻精品一区果冻| 国产精品人妻久久久久久| 伦理电影免费视频| 麻豆成人av视频| 久久综合国产亚洲精品| 久久鲁丝午夜福利片| 日本爱情动作片www.在线观看| 一个人看视频在线观看www免费| 综合色丁香网| 少妇高潮的动态图| 最后的刺客免费高清国语| 欧美区成人在线视频| 王馨瑶露胸无遮挡在线观看| 国产av码专区亚洲av| 日韩欧美精品免费久久| 大话2 男鬼变身卡| 最近2019中文字幕mv第一页| 免费人妻精品一区二区三区视频| 在线观看www视频免费| 久久99精品国语久久久| 极品教师在线视频| 亚洲国产av新网站| 精品亚洲成国产av| 中文精品一卡2卡3卡4更新| 在线观看一区二区三区激情| 大片电影免费在线观看免费| 国产精品国产av在线观看| 国产高清国产精品国产三级| 亚洲在久久综合| av国产精品久久久久影院| 麻豆成人av视频| 九色成人免费人妻av| 免费黄网站久久成人精品| 寂寞人妻少妇视频99o| 精品久久久久久久久av| 99热这里只有是精品在线观看| 少妇高潮的动态图| 九色成人免费人妻av| 如何舔出高潮| 亚洲电影在线观看av| 免费av中文字幕在线| 欧美丝袜亚洲另类| 日本免费在线观看一区| 亚洲精品国产成人久久av| 欧美成人精品欧美一级黄| 永久网站在线| 日本爱情动作片www.在线观看| 国产一区有黄有色的免费视频| 一级,二级,三级黄色视频| 美女福利国产在线| 人人妻人人看人人澡| 色婷婷av一区二区三区视频| 国产69精品久久久久777片| 日本-黄色视频高清免费观看| 亚洲欧美日韩卡通动漫| 在线播放无遮挡| 18禁在线播放成人免费| 少妇高潮的动态图| 一本大道久久a久久精品| 国产成人精品无人区| 国产 精品1| 亚洲av国产av综合av卡| 欧美成人精品欧美一级黄| 亚洲av成人精品一二三区| 新久久久久国产一级毛片| 日韩大片免费观看网站| 一级毛片aaaaaa免费看小| 丰满乱子伦码专区| 国产在线视频一区二区| 精品久久久精品久久久| 两个人的视频大全免费| 韩国高清视频一区二区三区| 视频区图区小说| 国产亚洲av片在线观看秒播厂| 大码成人一级视频| 亚洲欧美成人精品一区二区| 啦啦啦视频在线资源免费观看| 黑丝袜美女国产一区| 欧美 亚洲 国产 日韩一| 国内少妇人妻偷人精品xxx网站| 99热6这里只有精品| 成人黄色视频免费在线看| 亚洲人成网站在线播| av免费在线看不卡| 人妻 亚洲 视频| 亚洲国产精品一区三区| 三级经典国产精品| 女性被躁到高潮视频| 国产黄色免费在线视频| 一级av片app| 菩萨蛮人人尽说江南好唐韦庄| 午夜影院在线不卡| 成人18禁高潮啪啪吃奶动态图 | 18禁裸乳无遮挡动漫免费视频| 精品亚洲成a人片在线观看| 精品午夜福利在线看| 久久毛片免费看一区二区三区| 综合色丁香网| 美女脱内裤让男人舔精品视频| 国产片特级美女逼逼视频| 欧美亚洲 丝袜 人妻 在线| 中文字幕亚洲精品专区| 中国三级夫妇交换| 国产av一区二区精品久久| 一级毛片久久久久久久久女| 久久久久久久久久久免费av| 中文乱码字字幕精品一区二区三区| 婷婷色麻豆天堂久久| 亚洲国产欧美日韩在线播放 | 免费看不卡的av| 亚洲色图综合在线观看| 天堂中文最新版在线下载| 亚洲精品成人av观看孕妇| 久久国产乱子免费精品| 丝袜喷水一区| 我要看日韩黄色一级片| 欧美最新免费一区二区三区| 曰老女人黄片| 欧美日韩在线观看h| 亚洲精品乱码久久久v下载方式| 97超视频在线观看视频| 久久婷婷青草| 人人妻人人添人人爽欧美一区卜| 亚洲成人av在线免费| 久久久久精品性色| 91aial.com中文字幕在线观看| 久久久久久久久久久丰满| 日韩av免费高清视频| 在线精品无人区一区二区三| 国产在视频线精品| 亚洲av中文av极速乱| 国产91av在线免费观看| 国产亚洲欧美精品永久| 伊人久久精品亚洲午夜| 99久国产av精品国产电影| 亚洲精品久久久久久婷婷小说| 国产91av在线免费观看| 九色成人免费人妻av| 久久精品国产亚洲av涩爱| 午夜老司机福利剧场| 在线观看av片永久免费下载| 精品99又大又爽又粗少妇毛片| 国产精品不卡视频一区二区| 国产亚洲最大av| 在线精品无人区一区二区三| 高清欧美精品videossex| 蜜臀久久99精品久久宅男| 亚洲精品乱码久久久v下载方式| 九色成人免费人妻av| 日韩,欧美,国产一区二区三区| 一本大道久久a久久精品| 日产精品乱码卡一卡2卡三| 老熟女久久久| 伦理电影免费视频| 老司机影院毛片| 51国产日韩欧美| 日韩不卡一区二区三区视频在线| 国产伦精品一区二区三区四那| 久久97久久精品| 高清黄色对白视频在线免费看 | 亚洲一级一片aⅴ在线观看| 日韩在线高清观看一区二区三区| videossex国产| 婷婷色麻豆天堂久久| 久久精品国产a三级三级三级| 欧美国产精品一级二级三级 | 麻豆成人午夜福利视频| 国产午夜精品久久久久久一区二区三区| 国产日韩欧美亚洲二区| 欧美性感艳星| 免费观看在线日韩| 日本黄色片子视频| 国产一区亚洲一区在线观看| 一个人看视频在线观看www免费| 2021少妇久久久久久久久久久| 国产精品无大码| 午夜影院在线不卡| 免费观看性生交大片5| 午夜福利在线观看免费完整高清在| 另类精品久久| 自拍欧美九色日韩亚洲蝌蚪91 | 日韩免费高清中文字幕av| av在线观看视频网站免费| 精品亚洲成国产av| 欧美日韩视频高清一区二区三区二| av.在线天堂| 成人二区视频| 国产亚洲欧美精品永久| 伦理电影大哥的女人| 国产黄色视频一区二区在线观看| 免费播放大片免费观看视频在线观看| 成人午夜精彩视频在线观看| 久久午夜福利片| 免费高清在线观看视频在线观看| 极品少妇高潮喷水抽搐| 9色porny在线观看| 久久人妻熟女aⅴ| 国产一级毛片在线| 国产男女内射视频| 成人亚洲精品一区在线观看| 99热全是精品| 亚洲av.av天堂| 欧美+日韩+精品| av不卡在线播放| 最黄视频免费看| 国产在线男女| 一本久久精品| 国产熟女午夜一区二区三区 | 91午夜精品亚洲一区二区三区| 久久久久久久久久久丰满| 亚洲国产精品一区三区| 午夜福利在线观看免费完整高清在| 如何舔出高潮| √禁漫天堂资源中文www| 老司机影院毛片| 在线观看www视频免费| 国产亚洲91精品色在线| 看十八女毛片水多多多| 国产精品国产三级国产专区5o| 爱豆传媒免费全集在线观看| 国产欧美亚洲国产| 狂野欧美激情性bbbbbb| 色吧在线观看| 成年人免费黄色播放视频 | 久久这里有精品视频免费| av有码第一页| 久久人人爽人人爽人人片va| 国产av一区二区精品久久| 视频中文字幕在线观看| 国产黄片美女视频| 我要看黄色一级片免费的| 午夜福利网站1000一区二区三区| 观看美女的网站| 国产精品麻豆人妻色哟哟久久| 我要看日韩黄色一级片| 成人亚洲欧美一区二区av| 国产成人aa在线观看| 男女边摸边吃奶| 熟妇人妻不卡中文字幕| 五月玫瑰六月丁香| 少妇的逼水好多| 校园人妻丝袜中文字幕| 久久久国产精品麻豆| 大香蕉97超碰在线| 国产精品福利在线免费观看| 久久狼人影院| 午夜老司机福利剧场| 69精品国产乱码久久久| 国产成人午夜福利电影在线观看| 2021少妇久久久久久久久久久| 又大又黄又爽视频免费| 尾随美女入室| 丝袜脚勾引网站| 色吧在线观看| 国产精品一区二区三区四区免费观看| 久久人人爽人人片av| 日韩三级伦理在线观看| 成人亚洲精品一区在线观看| 妹子高潮喷水视频|