• <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中文字幕| 色5月婷婷丁香| 亚洲激情五月婷婷啪啪| 国产亚洲午夜精品一区二区久久 | 亚洲成色77777| 成人黄色视频免费在线看| 一级a做视频免费观看| 久久国内精品自在自线图片| 乱码一卡2卡4卡精品| 免费看不卡的av| 国产高潮美女av| 亚洲精品乱久久久久久| 婷婷色综合www| 在线亚洲精品国产二区图片欧美 | 午夜福利视频精品| 直男gayav资源| 伊人久久精品亚洲午夜| 99热国产这里只有精品6| 亚洲精品久久午夜乱码| 精品99又大又爽又粗少妇毛片| 久久久久久久久久成人| 哪个播放器可以免费观看大片| 偷拍熟女少妇极品色| 免费黄色在线免费观看| av在线亚洲专区| 国产极品天堂在线| 亚洲av一区综合| 免费在线观看成人毛片| 特级一级黄色大片| 色综合色国产| 日本色播在线视频| 免费观看无遮挡的男女| 日本三级黄在线观看| 免费大片黄手机在线观看| 亚洲成人精品中文字幕电影| 大片电影免费在线观看免费| 亚洲精品第二区| av专区在线播放| 两个人的视频大全免费| 亚洲精品国产色婷婷电影| 亚洲最大成人手机在线| 国产一区有黄有色的免费视频| 欧美区成人在线视频| www.色视频.com| 国产免费一区二区三区四区乱码| 国产一区有黄有色的免费视频| 精品熟女少妇av免费看| 国产精品三级大全| 色视频www国产| 日韩成人伦理影院| av免费观看日本| 日本爱情动作片www.在线观看| 婷婷色av中文字幕| 免费av不卡在线播放| 亚洲欧美成人综合另类久久久| 免费观看无遮挡的男女| 男女那种视频在线观看| 日韩一本色道免费dvd| 麻豆国产97在线/欧美| 简卡轻食公司| a级一级毛片免费在线观看| 中文天堂在线官网| 偷拍熟女少妇极品色| 超碰97精品在线观看| 搞女人的毛片| 精品一区二区免费观看| 好男人在线观看高清免费视频| 丝袜喷水一区| 在线观看av片永久免费下载| av在线天堂中文字幕| 日本黄色片子视频| 九草在线视频观看| 国产精品国产三级国产av玫瑰| 日本一二三区视频观看| 成人美女网站在线观看视频| 汤姆久久久久久久影院中文字幕| 久久精品久久精品一区二区三区| 伊人久久国产一区二区| 全区人妻精品视频| 亚洲最大成人av| 韩国av在线不卡| 日韩成人伦理影院| 精品一区二区免费观看| 亚洲人与动物交配视频| 亚洲欧美日韩东京热| 五月伊人婷婷丁香| 插阴视频在线观看视频| 亚洲丝袜综合中文字幕| 日韩欧美精品v在线| 亚洲欧美日韩卡通动漫| 高清毛片免费看| 九九在线视频观看精品| 国产成人福利小说| 免费看av在线观看网站| 欧美日韩综合久久久久久| 亚洲av免费在线观看| 国产淫片久久久久久久久| 成人午夜精彩视频在线观看| 亚洲精品日韩av片在线观看| 婷婷色综合www| 免费人成在线观看视频色| av专区在线播放| 永久免费av网站大全| 国产亚洲午夜精品一区二区久久 | 视频区图区小说| 亚洲高清免费不卡视频| 精品视频人人做人人爽| 亚洲天堂av无毛| 国产免费视频播放在线视频| av在线老鸭窝| 亚洲真实伦在线观看| 视频中文字幕在线观看| 别揉我奶头 嗯啊视频| 免费观看av网站的网址| 嫩草影院精品99| 欧美3d第一页| 人体艺术视频欧美日本| 激情五月婷婷亚洲| 久久久久久久大尺度免费视频| 色综合色国产| 色哟哟·www| 久久99热6这里只有精品| 亚洲在久久综合| 男女下面进入的视频免费午夜| 97热精品久久久久久| 真实男女啪啪啪动态图| 久久热精品热| 嫩草影院新地址| 亚洲精品乱码久久久v下载方式| 国产熟女欧美一区二区| 久久久精品欧美日韩精品| 一级av片app| 看十八女毛片水多多多| 麻豆国产97在线/欧美| 国产乱来视频区| 在线精品无人区一区二区三 | 国产成人免费无遮挡视频| 日日啪夜夜爽| 亚洲精品亚洲一区二区| 97热精品久久久久久| 黄色欧美视频在线观看| 街头女战士在线观看网站| 国产高清有码在线观看视频| 夜夜爽夜夜爽视频| 欧美日韩视频精品一区| 男插女下体视频免费在线播放| 美女被艹到高潮喷水动态| 国产亚洲5aaaaa淫片| 亚洲色图综合在线观看| 成人一区二区视频在线观看| 婷婷色综合大香蕉| 亚洲av免费高清在线观看| 亚洲av中文字字幕乱码综合| 国产一区二区在线观看日韩| 老司机影院成人| 午夜免费男女啪啪视频观看| 99热网站在线观看| 亚洲国产色片| 黄色配什么色好看| 一级片'在线观看视频| a级毛色黄片| 看十八女毛片水多多多| 精品一区在线观看国产| 久久精品国产亚洲av天美| 色网站视频免费| 一级片'在线观看视频| 国产精品一二三区在线看| 久久女婷五月综合色啪小说 | 免费黄网站久久成人精品| 全区人妻精品视频| 深爱激情五月婷婷| 亚洲精品一区蜜桃| 日韩一本色道免费dvd| 午夜亚洲福利在线播放| 免费黄频网站在线观看国产| 黄片无遮挡物在线观看| 久久国内精品自在自线图片| 国产精品不卡视频一区二区| 我的女老师完整版在线观看| 久久久久性生活片| 亚洲国产高清在线一区二区三| 久久久久久国产a免费观看| 性色av一级| 国产黄片美女视频| 久久影院123| 日本三级黄在线观看| 纵有疾风起免费观看全集完整版| 视频中文字幕在线观看| 精品久久久噜噜| 婷婷色麻豆天堂久久| 国产高清有码在线观看视频| 日韩一区二区视频免费看| 欧美日韩在线观看h| 插阴视频在线观看视频| 九草在线视频观看| 国产精品爽爽va在线观看网站| 国产精品久久久久久久久免| 午夜福利网站1000一区二区三区| av卡一久久| 国产色婷婷99| 精品久久久久久久久亚洲| 亚洲国产最新在线播放| 内射极品少妇av片p| 少妇人妻久久综合中文| 国产成人aa在线观看| 国产欧美亚洲国产| 久久久精品94久久精品| 又爽又黄a免费视频| 男女啪啪激烈高潮av片| 97热精品久久久久久| 午夜视频国产福利| 成人国产av品久久久| 国产成人精品久久久久久| 舔av片在线| 99久久人妻综合| 亚洲av中文av极速乱| 国产成人a∨麻豆精品| 国产探花在线观看一区二区| 搡女人真爽免费视频火全软件| 中文乱码字字幕精品一区二区三区| 欧美人与善性xxx| 97超视频在线观看视频| 国产成人aa在线观看| 国产毛片a区久久久久| 亚洲av福利一区| 性色av一级| 国产熟女欧美一区二区| 一区二区三区四区激情视频| 欧美高清成人免费视频www| 国产精品av视频在线免费观看| 亚洲精品日本国产第一区| 一个人看视频在线观看www免费| 亚洲av男天堂| 晚上一个人看的免费电影| 日本-黄色视频高清免费观看| 成人漫画全彩无遮挡| 尤物成人国产欧美一区二区三区| 黄片无遮挡物在线观看| 最近2019中文字幕mv第一页| 亚洲av日韩在线播放| 嫩草影院新地址| 夜夜爽夜夜爽视频| 亚洲av成人精品一二三区| 简卡轻食公司| 在线天堂最新版资源| 亚洲最大成人av| 精品人妻一区二区三区麻豆| 日产精品乱码卡一卡2卡三| 一级毛片我不卡| 久久久久久久久久成人| 丰满乱子伦码专区| 伊人久久国产一区二区| av免费观看日本| 麻豆成人午夜福利视频| 国产高清不卡午夜福利| 欧美zozozo另类| 国产探花在线观看一区二区| 99久久中文字幕三级久久日本| 欧美丝袜亚洲另类| 日韩欧美 国产精品| 91在线精品国自产拍蜜月| 久久精品熟女亚洲av麻豆精品| 欧美日韩精品成人综合77777| 国产精品久久久久久精品电影小说 | 男女无遮挡免费网站观看| 五月伊人婷婷丁香| 久久久久性生活片| 在线观看av片永久免费下载| 热99国产精品久久久久久7| 国产成人freesex在线| 亚洲内射少妇av| 中文字幕免费在线视频6| 国内少妇人妻偷人精品xxx网站| 午夜福利视频精品| 国产一区二区三区综合在线观看 | 国产老妇伦熟女老妇高清| 国产白丝娇喘喷水9色精品| 久久精品综合一区二区三区| 在线观看人妻少妇| 日韩不卡一区二区三区视频在线| 亚洲成人一二三区av| 国产乱人视频| 久久久色成人| 99精国产麻豆久久婷婷| 一二三四中文在线观看免费高清| 尾随美女入室| 成人特级av手机在线观看| 一级爰片在线观看| 又大又黄又爽视频免费| 国产精品蜜桃在线观看| 免费观看a级毛片全部| 国产精品99久久99久久久不卡 | 久久人人爽人人片av| 成年女人看的毛片在线观看| 六月丁香七月| 久久久久国产网址| 欧美性猛交╳xxx乱大交人| 国产免费又黄又爽又色| 99热网站在线观看| 久久精品国产a三级三级三级| 国产 一区精品| 精品一区二区三区视频在线| av.在线天堂| 亚洲aⅴ乱码一区二区在线播放| 欧美老熟妇乱子伦牲交| a级毛片免费高清观看在线播放| 涩涩av久久男人的天堂| 视频区图区小说| 如何舔出高潮| 日韩欧美精品v在线| 51国产日韩欧美| 日本免费在线观看一区| 亚洲精品影视一区二区三区av| 欧美日本视频| 欧美另类一区| 80岁老熟妇乱子伦牲交| 久久久久久久午夜电影| 成人一区二区视频在线观看| 国产av不卡久久| 别揉我奶头 嗯啊视频| 欧美日韩视频精品一区| 成人鲁丝片一二三区免费| 交换朋友夫妻互换小说| 18禁裸乳无遮挡免费网站照片| 成人亚洲精品一区在线观看 | 国产高清有码在线观看视频| 搞女人的毛片| 久久久久久久精品精品| 18+在线观看网站| 亚洲一级一片aⅴ在线观看| 亚洲国产成人一精品久久久| 欧美日韩视频精品一区| 美女视频免费永久观看网站| 成人黄色视频免费在线看| 亚洲欧美日韩另类电影网站 | av国产免费在线观看| 青春草国产在线视频| 男女啪啪激烈高潮av片| 一本一本综合久久| 国产免费福利视频在线观看| 日本一二三区视频观看| 欧美日韩一区二区视频在线观看视频在线 | 91精品国产九色| 亚洲精品自拍成人| 精品人妻熟女av久视频| 在线天堂最新版资源| 丰满少妇做爰视频| 99精国产麻豆久久婷婷| 成人无遮挡网站| 成年版毛片免费区| 在线观看免费高清a一片| 久久精品国产亚洲av涩爱| 欧美性猛交╳xxx乱大交人| 人妻 亚洲 视频| 亚洲av免费在线观看| 尤物成人国产欧美一区二区三区| 精品一区二区免费观看| 伊人久久国产一区二区| 99久久中文字幕三级久久日本| 尤物成人国产欧美一区二区三区| 99热网站在线观看| 乱码一卡2卡4卡精品| 九九久久精品国产亚洲av麻豆| 亚洲av不卡在线观看| 国产精品偷伦视频观看了| 国产美女午夜福利| 国产黄色免费在线视频| 女人被狂操c到高潮| 亚洲精品久久久久久婷婷小说| 97精品久久久久久久久久精品| 身体一侧抽搐| 亚洲国产精品成人久久小说| 久久热精品热| 成年免费大片在线观看| 深夜a级毛片| 99热这里只有是精品在线观看| 中文字幕av成人在线电影| 在现免费观看毛片| 亚洲精品aⅴ在线观看| 国内精品宾馆在线| 国产永久视频网站| 亚洲成人久久爱视频| 边亲边吃奶的免费视频| 国产成人aa在线观看| 成人毛片a级毛片在线播放| 国产精品熟女久久久久浪| 制服丝袜香蕉在线| 水蜜桃什么品种好| 国产亚洲5aaaaa淫片| 韩国av在线不卡| 亚洲av国产av综合av卡| 热re99久久精品国产66热6| 亚洲欧美日韩无卡精品| 国产一级毛片在线| 国产一区二区在线观看日韩| 丝瓜视频免费看黄片| 99热国产这里只有精品6| 欧美极品一区二区三区四区| 亚洲欧美成人综合另类久久久| 亚洲不卡免费看| 九九久久精品国产亚洲av麻豆| 国产在线一区二区三区精| 99久久精品热视频| 老师上课跳d突然被开到最大视频| 99久久精品热视频| 日本一本二区三区精品| 精品久久久久久久久av| 亚洲精华国产精华液的使用体验| 偷拍熟女少妇极品色| 国产伦在线观看视频一区| 久久综合国产亚洲精品| www.av在线官网国产| 18禁在线无遮挡免费观看视频| av天堂中文字幕网| 婷婷色综合www| 亚洲第一区二区三区不卡| 免费观看的影片在线观看| 日本黄大片高清| 日韩 亚洲 欧美在线| 日韩人妻高清精品专区| 99热这里只有是精品在线观看| 国产精品福利在线免费观看| 特大巨黑吊av在线直播| 亚洲美女视频黄频| av在线老鸭窝| 夫妻午夜视频| 亚洲国产精品999| 午夜免费鲁丝| 成人亚洲精品一区在线观看 | 人妻系列 视频| 2022亚洲国产成人精品| 美女视频免费永久观看网站| 搡女人真爽免费视频火全软件| 日日啪夜夜爽| 亚洲欧洲日产国产| 国产伦理片在线播放av一区| 91在线精品国自产拍蜜月| 99久久九九国产精品国产免费| 久久韩国三级中文字幕| 肉色欧美久久久久久久蜜桃 | 久久国产乱子免费精品| 欧美极品一区二区三区四区| 日韩制服骚丝袜av| 99热这里只有是精品50| 亚洲国产欧美人成| 人人妻人人爽人人添夜夜欢视频 | 免费播放大片免费观看视频在线观看| 尤物成人国产欧美一区二区三区| 午夜激情久久久久久久| 人人妻人人看人人澡| 超碰av人人做人人爽久久| 亚洲最大成人中文| 天天躁日日操中文字幕| 免费观看av网站的网址| 麻豆乱淫一区二区| 国产一区二区在线观看日韩| videos熟女内射| 免费不卡的大黄色大毛片视频在线观看| 国产探花在线观看一区二区| 国产精品精品国产色婷婷| 黄色配什么色好看| 高清欧美精品videossex| 亚洲精品aⅴ在线观看| 免费在线观看成人毛片| 少妇丰满av| 99热这里只有是精品50| 亚洲精品国产色婷婷电影| 国产精品99久久久久久久久| 亚洲国产精品成人综合色| 别揉我奶头 嗯啊视频| 亚洲av.av天堂| videos熟女内射| 国产成人福利小说| 中文资源天堂在线| 大片免费播放器 马上看| 免费观看在线日韩| 韩国av在线不卡| av国产免费在线观看| 久久久久久久久久久丰满| 在线观看三级黄色| 亚洲国产精品999| 久久韩国三级中文字幕| 亚洲欧美精品自产自拍| 国产中年淑女户外野战色| 久久99蜜桃精品久久| 亚洲美女搞黄在线观看| 九色成人免费人妻av| 欧美老熟妇乱子伦牲交| 91精品伊人久久大香线蕉| 日本免费在线观看一区| 99热全是精品| www.色视频.com| 成人欧美大片| 久久久久久久精品精品| 插逼视频在线观看| 久久久久国产精品人妻一区二区| eeuss影院久久| 在线观看人妻少妇| av一本久久久久| 久久久久国产网址| 日本欧美国产在线视频| 国产一区二区三区av在线| 亚洲av在线观看美女高潮| 在线观看三级黄色| 国产精品99久久久久久久久| 午夜福利视频1000在线观看| 神马国产精品三级电影在线观看| 国产爽快片一区二区三区| 黑人高潮一二区| 亚洲欧美成人精品一区二区| 亚洲精品影视一区二区三区av| 永久免费av网站大全| 欧美3d第一页| 七月丁香在线播放| 亚洲人成网站在线观看播放| 欧美激情在线99| 国产免费一级a男人的天堂| 欧美3d第一页| 成人国产av品久久久| 亚洲精品日韩在线中文字幕| 狠狠精品人妻久久久久久综合| 国产欧美另类精品又又久久亚洲欧美| 成年女人在线观看亚洲视频 | 夜夜看夜夜爽夜夜摸| 少妇 在线观看| 亚洲欧洲国产日韩| 秋霞伦理黄片| 一区二区三区四区激情视频| av在线蜜桃| 亚洲成人一二三区av| 国产精品一及| 国产精品国产av在线观看| 女人十人毛片免费观看3o分钟| 成人毛片60女人毛片免费| 国产精品伦人一区二区| 免费不卡的大黄色大毛片视频在线观看| 免费看不卡的av| 91精品伊人久久大香线蕉| 日韩av不卡免费在线播放| 狠狠精品人妻久久久久久综合| 在线观看三级黄色| 欧美人与善性xxx| 国产在视频线精品| 日本与韩国留学比较| 国产高清三级在线| 99九九线精品视频在线观看视频| 美女脱内裤让男人舔精品视频| 免费在线观看成人毛片| 中国三级夫妇交换| 伦精品一区二区三区| 欧美成人一区二区免费高清观看| 亚洲激情五月婷婷啪啪| 国产在线一区二区三区精| 久久久精品欧美日韩精品| 国产老妇女一区| 一级黄片播放器| www.色视频.com| 18+在线观看网站| 亚洲va在线va天堂va国产| 亚洲精品久久久久久婷婷小说| 国产在线一区二区三区精| 成人特级av手机在线观看| 欧美最新免费一区二区三区| 国产精品一二三区在线看| 国产精品一及| 人妻一区二区av| 久久国内精品自在自线图片| 日韩中字成人| 国产欧美日韩精品一区二区| 国产免费视频播放在线视频| 22中文网久久字幕| 99久久人妻综合| 欧美性感艳星| 永久网站在线| 国产一级毛片在线| 少妇的逼好多水| av福利片在线观看| 国产视频内射| av免费在线看不卡| 麻豆成人av视频| 久久国内精品自在自线图片| 少妇被粗大猛烈的视频| 成人欧美大片| 亚洲精品影视一区二区三区av| 神马国产精品三级电影在线观看| 又黄又爽又刺激的免费视频.| 99九九线精品视频在线观看视频| 国产av码专区亚洲av| 熟女人妻精品中文字幕| 草草在线视频免费看| 18+在线观看网站| 最近手机中文字幕大全| 日本与韩国留学比较| 国产成人a区在线观看| 久久99热这里只频精品6学生| 久久久久久九九精品二区国产| av线在线观看网站| 中国国产av一级| 99视频精品全部免费 在线| 男人狂女人下面高潮的视频| 国产亚洲午夜精品一区二区久久 | 国产美女午夜福利| 亚洲怡红院男人天堂| 日日摸夜夜添夜夜添av毛片| 好男人在线观看高清免费视频| 日韩精品有码人妻一区| 日韩欧美一区视频在线观看 | 精品少妇久久久久久888优播| 两个人的视频大全免费| 蜜桃久久精品国产亚洲av| 男女国产视频网站| 肉色欧美久久久久久久蜜桃 |