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

    錐度和不對(duì)中誤差對(duì)止推箔片軸承承載特性的影響

    2023-10-17 14:43:22連軍偉張茂森金志磊賀雷
    軸承 2023年10期
    關(guān)鍵詞:承載力

    連軍偉,張茂森,金志磊,賀雷

    (1.中國(guó)航天科技集團(tuán)有限公司 低溫液體推進(jìn)技術(shù)實(shí)驗(yàn)室,北京 100076;2.北京航天動(dòng)力研究所,北京 100076)

    1 概述

    人類探測(cè)太空的能力主要受空間推進(jìn)技術(shù)的限制,針對(duì)高速度增量的航天任務(wù),傳統(tǒng)化學(xué)推進(jìn)技術(shù)的比沖較低,太陽(yáng)能電池陣—蓄電池技術(shù)嚴(yán)重依賴太陽(yáng)光照,深空生存能力較差[1]??臻g核電源系統(tǒng)具有能量密度大,容易實(shí)現(xiàn)大功率(數(shù)千瓦至數(shù)兆瓦)供電,機(jī)動(dòng)性和隱蔽性好,不依賴太陽(yáng)光照等優(yōu)勢(shì),是深空探測(cè)不可替代的空間電源[1]。

    動(dòng)態(tài)轉(zhuǎn)換的放射性同位素發(fā)電系統(tǒng)在空間核電源系統(tǒng)中轉(zhuǎn)換效率最高,其中研究較多且已應(yīng)用的是閉式布雷頓循環(huán)的動(dòng)力轉(zhuǎn)換系統(tǒng)(圖1),旋轉(zhuǎn)單元是該系統(tǒng)的核心組件,其結(jié)構(gòu)和性能直接影響系統(tǒng)的優(yōu)劣。止推箔片軸承作為一種柔性自適應(yīng)動(dòng)壓軸承,在旋轉(zhuǎn)單元中承受軸向載荷,與滾動(dòng)軸承相比,省略了復(fù)雜的冷卻潤(rùn)滑系統(tǒng),使系統(tǒng)更加清潔和高效。

    (a) 閉式布雷頓循環(huán)系統(tǒng)原理

    1854年,法國(guó)科學(xué)家提出將氣體作為潤(rùn)滑劑的可能,1886年法國(guó)科學(xué)家雷諾推導(dǎo)出潤(rùn)滑膜內(nèi)壓力分布的雷諾方程[2],對(duì)流體潤(rùn)滑的認(rèn)識(shí)上升到理論高度。止推箔片軸承作為階梯型動(dòng)壓軸承的衍生結(jié)構(gòu)軸承,具有波箔剛度非線性的特點(diǎn),其理論研究晚于徑向箔片軸承且相關(guān)研究較少。1983年,文獻(xiàn)[3]對(duì)止推箔片軸承的靜態(tài)特性進(jìn)行了詳細(xì)的計(jì)算分析,計(jì)算模型考慮了波箔剛度的非線性分布,第1次使用柔度系數(shù)α分析箔片軸承靜特性,該方法一直沿用至今。2000年,文獻(xiàn)[4]采用有限差分法與有限元相結(jié)合的方法,將可壓縮流體的雷諾方程求解與彈性變形耦合起來(lái)進(jìn)行止推箔片軸承分析。2017年,文獻(xiàn)[5]采用非等溫雷諾方程評(píng)估了轉(zhuǎn)子靜態(tài)和動(dòng)態(tài)角度不對(duì)中對(duì)止推箔片軸承性能的影響,并深入分析了軸承的承載力、摩擦功、氣膜剛度和阻尼系數(shù)的變化情況。

    20世紀(jì)80年代,文獻(xiàn)[6]應(yīng)用圓柱薄殼線性理論推導(dǎo)出波箔變形方程,這是國(guó)內(nèi)最早關(guān)于箔片軸承中波箔剛度的計(jì)算模型。文獻(xiàn)[7]采用有限單元法對(duì)平箔式止推箔片軸承進(jìn)行求解,并分析了箔片軸承靜態(tài)特性的影響因素。文獻(xiàn)[8]對(duì)止推箔片軸承進(jìn)行了啟停試驗(yàn),分析了箔片材料和表面涂層對(duì)軸承啟停磨損的影響。文獻(xiàn)[9]采用有限差分法得出,推力盤傾斜狀態(tài)下軸承的軸向承載力和黏性摩擦力矩均隨傾角的增大而增大,且傾角越大靜態(tài)特性越大。文獻(xiàn)[10]基于牛頓迭代法建立偏角不對(duì)中時(shí)止推箔片軸承的潤(rùn)滑模型,與線性替換法的比較說(shuō)明牛頓迭代法具有同樣的適用性。文獻(xiàn)[11]基于CFX的有限元計(jì)算與MATLAB的有限差分求解動(dòng)壓氣體軸承氣膜壓力,二者計(jì)算結(jié)果一致,隨著偏心率的增加CFX有限元的計(jì)算結(jié)果比有限差分法計(jì)算的略小。文獻(xiàn)[12]采用ANSYS和CFX軟件計(jì)算了超臨界二氧化碳止推箔片軸承動(dòng)態(tài)特性,得出動(dòng)態(tài)剛度和阻尼大小主要取決于氣膜厚度和箔片結(jié)構(gòu),且隨氣膜厚度的增大而減小。隨著更多類型箔片軸承的研究,箔片軸承的性能逐漸被熟知,但仍缺乏精準(zhǔn)的理論分析模型和實(shí)用的設(shè)計(jì)規(guī)范[13]。

    綜上可知,目前在止推箔片軸承的誤差研究方面,主要著重于軸承不對(duì)中誤差對(duì)其靜態(tài)性能的影響,而錐度誤差對(duì)軸承承載特性的影響以及不對(duì)中誤差下軸承動(dòng)態(tài)特性變化的研究較少。在生產(chǎn)和裝配過(guò)程中,錐度和不對(duì)中誤差均難以避免,不對(duì)中誤差主要有平行不對(duì)中、角不對(duì)中和綜合不對(duì)中。由于旋轉(zhuǎn)單元轉(zhuǎn)子系統(tǒng)的徑向設(shè)計(jì)間隙約束及推力盤徑向尺寸大于止推軸承高度,分析過(guò)程不考慮平行不對(duì)中的影響,本文所提及的不對(duì)中特指角不對(duì)中或主軸傾斜。本文在柱坐標(biāo)系下建立止推箔片軸承的雷諾方程,并采用牛頓迭代法和有限差分法對(duì)雷諾方程和氣膜厚度方程進(jìn)行耦合求解,分析生產(chǎn)和裝配過(guò)程引起的徑向錐度誤差和不對(duì)中誤差對(duì)止推箔片軸承承載特性的影響。

    2 理論分析

    止推箔片軸承主要由頂箔、波箔和底板組成。波箔和頂箔組成多個(gè)扇形瓦,每個(gè)瓦的度數(shù)和傾斜面占比相同。止推箔片軸承結(jié)構(gòu)和箔片參數(shù)示意圖如圖2所示,β為扇形瓦的瓦角,βp為扇形瓦節(jié)線位置角,σh為扇形瓦傾斜面高度,ω為推力盤角速度,h2為設(shè)計(jì)間隙。

    (a) 止推箔片軸承 (b) 箔片參數(shù)

    2.1 壓力控制雷諾方程推導(dǎo)

    推導(dǎo)等溫條件下的雷諾方程時(shí)進(jìn)行如下假設(shè):1)忽略氣體慣性力;2)潤(rùn)滑氣體介質(zhì)視為牛頓流體,符合牛頓黏性定律;3)氣膜間隙內(nèi)為層流;4)忽略氣膜法向流動(dòng);5)潤(rùn)滑介質(zhì)視為理想氣體。

    經(jīng)推導(dǎo)得到柱坐標(biāo)下潤(rùn)滑介質(zhì)為理想氣體的雷諾方程,即

    (1)

    式中:r為徑向坐標(biāo);h為氣膜厚度;P為氣膜壓力;θ為周向坐標(biāo);μ為氣體動(dòng)力黏度。

    (2)

    式中:Pa為環(huán)境壓力;R2為軸承外半徑;Λ為軸承數(shù)。

    2.2 氣膜厚度方程

    每個(gè)扇形瓦的氣膜由傾斜面氣膜和水平面氣膜組成,根據(jù)圖2的結(jié)構(gòu)參數(shù)得到單個(gè)扇形瓦的氣膜厚度方程,即

    h=h2+h(θ)+w(r,θ) ,

    (3)

    式中:w(r,θ)為箔片變形后構(gòu)建的變形區(qū)域氣膜厚度;h(θ)為扇形瓦傾斜面氣膜厚度。

    定義b為斜面占比,b=βP/β,氣膜厚度的表達(dá)式為

    (4)

    在工程應(yīng)用過(guò)程中存在加工誤差,推力盤不可避免地產(chǎn)生一定的徑向錐度βz,將直接影響氣膜厚度分布。錐度為負(fù)數(shù)時(shí)表示推力盤出現(xiàn)內(nèi)錐狀態(tài),為正數(shù)時(shí)出現(xiàn)外錐狀態(tài), 按內(nèi)錐度和外錐度這2種形式構(gòu)建徑向錐度模型(圖3),Ri為錐度半徑,R1為軸承內(nèi)半徑,Δh為錐度引起的氣膜厚度變化。錐度βz的定義為

    (a) 內(nèi)錐度

    (5)

    建立考慮推力盤徑向錐度的氣膜厚度方程,即

    h0=h2+h(θ)+w(r,θ)+Δh。

    (6)

    盡管設(shè)計(jì)及生產(chǎn)過(guò)程要嚴(yán)格控制形位公差,但在裝配過(guò)程仍會(huì)存在不對(duì)中現(xiàn)象[14-15],而不對(duì)中誤差會(huì)導(dǎo)致止推箔片軸承氣膜間隙發(fā)生變化,從而引起軸承承載特性的改變。如圖4所示,建立不對(duì)中模型和空間坐標(biāo)系。

    (a) 坐標(biāo)系及扇形瓦標(biāo)記

    假設(shè)轉(zhuǎn)子與止推箔片軸承的夾角(不對(duì)中角度)為φ,在空間坐標(biāo)下可分解為繞x軸的θx和繞y軸的θy??紤]轉(zhuǎn)子不對(duì)中的氣膜厚度方程為

    h1=h2+h(θ)+w(r,θ)-rcosθtanφ。

    (7)

    2.3 波箔剛度等效模型

    波箔單位面積剛度Kb的表達(dá)式為

    (8)

    式中:Eb為波箔的彈性模量;t為波箔厚度;l為半波箔長(zhǎng)度;s為波箔單位長(zhǎng)度;ν為波箔的泊松比。

    2.4 網(wǎng)格劃分及邊界條件

    圖5 單個(gè)扇形瓦的網(wǎng)格劃分

    2.5 流固耦合計(jì)算流程

    采用有限差分法求解(2)式得到氣膜壓力,通過(guò)迭代由(9)式判斷是否收斂?;贛ATLAB平臺(tái)按圖6程序邏輯進(jìn)行編譯計(jì)算。

    圖6 止推箔片軸承計(jì)算程序邏輯圖

    (9)

    3 計(jì)算與分析

    3.1 程序驗(yàn)證

    (10)

    式中:W為承載力。

    本文程序計(jì)算與文獻(xiàn)[3]的對(duì)比結(jié)果如圖7所示,最大誤差為2.53%,兩者吻合度較好,說(shuō)明本文止推箔片軸承數(shù)學(xué)模型的準(zhǔn)確性較高。

    圖7 本文程序與文獻(xiàn)[3]計(jì)算結(jié)果對(duì)比

    3.2 錐度誤差的影響

    空間核電動(dòng)力轉(zhuǎn)換系統(tǒng)的止推箔片軸承結(jié)構(gòu)參數(shù)及其他參數(shù)見表1。

    表1 止推箔片軸承結(jié)構(gòu)參數(shù)及其他參數(shù)

    錐度一定時(shí),各扇形瓦氣膜厚度和氣膜壓力在同一半徑和周向角度上相同,為提高計(jì)算精度并減少計(jì)算時(shí)間, 選取單個(gè)扇形瓦進(jìn)行分析。 在額定設(shè)計(jì)間隙下,內(nèi)錐(Ri=R1)和外錐(Ri=R2)狀態(tài)下量綱一的氣膜壓力和氣膜厚度分布如圖8和圖9所示:當(dāng)內(nèi)錐度為-4×10-4時(shí),最大氣膜壓力增大18.4%,對(duì)應(yīng)軸向承載力提高47.8%;當(dāng)外錐度為4×10-4時(shí),最大氣膜壓力減小12.7%,對(duì)應(yīng)軸向承載力降低33.8%。以圖9b中氣膜厚度A和B位置為基準(zhǔn),內(nèi)錐度為-4×10-4時(shí)B區(qū)域的氣膜厚度減小至10.6 μm,外錐度為4×10-4時(shí)A區(qū)域的氣膜厚度增大至29.4 μm。綜上可知:內(nèi)錐度減小了平均氣膜厚度,提高了承載力;外錐度增大了平均氣膜厚度,降低了承載力。

    (a) -4×10-4

    (a) -4×10-4

    設(shè)計(jì)間隙不變時(shí)軸向承載力和摩擦功隨錐度的變化以及軸向載荷恒定(200 N)時(shí)最小氣膜厚度和量綱一的最大氣膜壓力隨錐度的變化如圖10所示:設(shè)計(jì)間隙不變時(shí),軸向承載力和摩擦功隨著錐度的增加逐漸減小,由于錐度的增加增大了平均氣膜厚度,削弱了動(dòng)壓效應(yīng)進(jìn)而減小了軸向承載力和摩擦功。在軸向載荷恒定時(shí),隨著錐度絕對(duì)值的增加,最小氣膜厚度逐漸減小, 最小氣膜厚度達(dá)到0.67 μm,若繼續(xù)減小至零時(shí)程序計(jì)算很難收斂,由于模型未考慮表面粗糙度, 工程應(yīng)用時(shí)可認(rèn)為此時(shí)軸承已處于臨界碰磨狀態(tài)。說(shuō)明恒定載荷工況時(shí)只有允許一定范圍的錐度誤差才能保證軸承正常運(yùn)行。最大氣膜壓力與最小氣膜厚度正好相反,說(shuō)明軸向載荷恒定時(shí)只要錐度絕對(duì)值增加就會(huì)增大最大氣膜壓力。此外,內(nèi)錐度和外錐度引起的最大氣膜壓力和最小氣膜厚度均關(guān)于錐度βz=0對(duì)稱。

    圖10 止推箔片軸承靜態(tài)參數(shù)隨錐度的變化

    3.3 不對(duì)中的影響

    為簡(jiǎn)化分析模型,掌握不對(duì)中對(duì)軸承性能的影響,只選取繞某一軸傾斜產(chǎn)生的不對(duì)中進(jìn)行分析。對(duì)中和不對(duì)中時(shí)止推箔片軸承的氣膜壓力、氣膜厚度和箔片變形如圖11所示:由于不對(duì)中的產(chǎn)生,軸承各瓦的氣膜壓力和氣膜厚度分布不均勻,箔片產(chǎn)生的變形量也不同;2#瓦最大氣膜壓力(2.458)最大,高于對(duì)中狀態(tài)下各瓦的最大氣膜壓力(1.805),5#瓦的最大氣膜壓力最小(1.494);由于氣膜壓力作用,最大變形量為(1.205 μm)同樣出現(xiàn)在2#瓦上,是5#瓦最大變形量(0.412 μm)的2.925倍;由于箔片變形構(gòu)建了新的氣膜承載區(qū)域,氣膜壓力越高時(shí)箔片變形量越大,體現(xiàn)了止推箔片軸承較好的自適應(yīng)性。

    (a) 不對(duì)中狀態(tài)(θx=0.01°)

    止推箔片軸承和各瓦的承載力和摩擦功隨不對(duì)中角度的變化如圖12所示:軸承的承載力和摩擦功隨著不對(duì)中角度的增大而增加,但各瓦不同;當(dāng)推力盤按圖4坐標(biāo)系傾斜時(shí),x軸上方的3個(gè)扇形瓦平均氣膜厚度減小,動(dòng)壓效應(yīng)增強(qiáng),氣膜壓力升高,承載力和摩擦功隨不對(duì)中角度的增大而增加,2#瓦平均氣膜厚度最小,其承載力最大且增速最快,1#瓦次之,3#瓦增速最慢,x軸下方3個(gè)扇形瓦則相反;不對(duì)中造成各瓦的承載力和摩擦發(fā)熱嚴(yán)重不均衡,不對(duì)中角度越大時(shí)各瓦的差異越大。

    (a) 軸承的承載力和摩擦功

    由于不對(duì)中的存在,止推箔片軸承不僅產(chǎn)生沿軸向的位移剛度,還存在沿徑向的彎曲剛度。圖4中軸承產(chǎn)生繞x軸的彎矩Mx; 由于1#, 6#瓦與3#, 4#瓦承載力的差異及2#, 5#瓦關(guān)于y軸不對(duì)稱,產(chǎn)生繞y軸的彎矩My,則軸承在3個(gè)自由度上有9個(gè)剛度系數(shù)和9個(gè)阻尼系數(shù),組成的剛度系數(shù)矩陣K和阻尼系數(shù)矩陣C分別為

    (11)

    (12)

    剛度系數(shù)和阻尼系數(shù)隨不對(duì)中角度的變化如圖13所示:當(dāng)繞x軸傾斜時(shí),軸向主剛度系數(shù)和主阻尼系數(shù)比軸向交叉項(xiàng)的大,隨著不對(duì)中角度的增大,軸向主剛度系數(shù)和主阻尼系數(shù)變化較明顯;在繞x軸方向,交叉項(xiàng)Kzθx和Czθx隨不對(duì)中角度的變化較明顯,且幅值遠(yuǎn)高于其他項(xiàng),主要原因?yàn)椴粚?duì)中角度增加,承載力快速增大,造成承載力在該方向的交叉彎曲剛度系數(shù)和阻尼系數(shù)增大;在繞y軸方向,交叉項(xiàng)Czθy隨不對(duì)中角度變化較明顯且為負(fù)數(shù),剛度系數(shù)主項(xiàng)Kθyθy的變化最大,當(dāng)θx為0.016°時(shí),主項(xiàng)Kθyθy增幅為28.5%,交叉項(xiàng)Kzθy增幅為22.1%,說(shuō)明在繞y軸方向,隨著不對(duì)中角度的增加,由彎矩My和承載力引起的剛度變化占主導(dǎo)。上述現(xiàn)象說(shuō)明不對(duì)中角度對(duì)止推箔片軸承的剛度系數(shù)和阻尼系數(shù)的影響較大,而剛度系數(shù)和阻尼系數(shù)直接影響轉(zhuǎn)子系統(tǒng)的穩(wěn)定性。

    4 結(jié)論

    基于牛頓迭代法和有限差分法,耦合求解考慮徑向錐度和不對(duì)中誤差影響的氣膜厚度方程及壓力控制的雷諾方程,對(duì)閉式布雷頓循環(huán)系統(tǒng)的止推箔片軸承進(jìn)行承載特性影響因素分析,得到如下結(jié)論:

    1)在設(shè)計(jì)間隙一定的情況下,內(nèi)錐度會(huì)減小平均氣膜厚度,進(jìn)而增大軸承承載力,當(dāng)內(nèi)錐度為-4×10-4時(shí),承載力增大47.8%;外錐度會(huì)增大平均氣膜厚度,進(jìn)而降低承載力,當(dāng)外錐度為4×10-4時(shí),承載力降幅達(dá)到33.8%。在軸向載荷為200 N時(shí),最小氣膜厚度隨著錐度絕對(duì)值的增加而減小,當(dāng)錐度大于2.25×10-3時(shí),軸承進(jìn)入臨界碰磨狀態(tài)。生產(chǎn)過(guò)程應(yīng)嚴(yán)格控制推力盤的錐度,避免因加工誤差過(guò)大而造成軸承失效。

    2)不對(duì)中會(huì)改變軸承各瓦的氣膜壓力和氣膜厚度分布,當(dāng)不對(duì)中角度θx=0.01°時(shí),最大氣膜壓力增加36.2%。不對(duì)中角度造成各瓦的承載力嚴(yán)重不均衡,對(duì)軸承剛度系數(shù)和阻尼系數(shù)的影響較明顯,直接影響轉(zhuǎn)子系統(tǒng)的穩(wěn)定性。隨著不對(duì)中角度的增大,軸承的摩擦功呈指數(shù)形式增加,導(dǎo)致軸承的發(fā)熱量急速上升,裝配過(guò)程應(yīng)盡可能將不對(duì)中誤差降到最低。

    猜你喜歡
    承載力
    再生混凝土抗剪鍵接縫受剪性能及承載力計(jì)算
    基于單樁豎向承載力計(jì)算分析研究
    CFRP-PCP板加固混凝土梁的抗彎承載力研究
    南通沿海開發(fā)水資源承載力的實(shí)證研究
    PVA-ECC抗剪加固帶懸臂RC梁承載力計(jì)算研究
    耐火鋼圓鋼管混凝土柱耐火極限和承載力
    基于SAP2000的光伏固定支架結(jié)構(gòu)承載力分析
    FRP筋混凝土板抗彎承載力理論研究
    K形偏心支撐鋼框架的彈性抗側(cè)剛度與極限承載力
    潛艇極限承載力計(jì)算與分析
    1000部很黄的大片| 少妇被粗大猛烈的视频| 久久精品夜夜夜夜夜久久蜜豆| 欧美性猛交黑人性爽| av视频在线观看入口| 桃色一区二区三区在线观看| 日日摸夜夜添夜夜添av毛片 | 亚洲不卡免费看| 99热这里只有是精品50| 亚洲va日本ⅴa欧美va伊人久久| 夜夜看夜夜爽夜夜摸| 婷婷六月久久综合丁香| 亚洲无线在线观看| 一级黄片播放器| 可以在线观看毛片的网站| 亚洲人成网站在线播放欧美日韩| 97人妻精品一区二区三区麻豆| 少妇被粗大猛烈的视频| 欧美成人a在线观看| 国产日本99.免费观看| 日本一本二区三区精品| 欧美日韩黄片免| 香蕉av资源在线| 女人十人毛片免费观看3o分钟| 亚洲成a人片在线一区二区| 色综合婷婷激情| 激情在线观看视频在线高清| 精品久久久久久久久亚洲 | 成人性生交大片免费视频hd| 亚洲人成网站高清观看| 欧美日本视频| 欧美不卡视频在线免费观看| 国产一区二区三区视频了| 99国产极品粉嫩在线观看| 国产一级毛片七仙女欲春2| 亚洲,欧美,日韩| 亚洲中文字幕日韩| 日韩免费av在线播放| 日本免费一区二区三区高清不卡| 午夜福利18| 欧美bdsm另类| 天堂动漫精品| 亚洲自偷自拍三级| 国产一区二区三区在线臀色熟女| 欧美黑人巨大hd| 精品欧美国产一区二区三| 精品人妻熟女av久视频| 如何舔出高潮| 欧美潮喷喷水| 欧美不卡视频在线免费观看| 午夜福利高清视频| 此物有八面人人有两片| 最新中文字幕久久久久| 在线免费观看的www视频| 黄色视频,在线免费观看| 无人区码免费观看不卡| 人人妻人人看人人澡| 1000部很黄的大片| 99在线人妻在线中文字幕| 99在线视频只有这里精品首页| 男女那种视频在线观看| 亚洲精品亚洲一区二区| 亚洲最大成人av| 中文字幕精品亚洲无线码一区| 免费在线观看日本一区| 久久热精品热| 偷拍熟女少妇极品色| 国产视频内射| 亚洲国产精品合色在线| 他把我摸到了高潮在线观看| 中亚洲国语对白在线视频| 蜜桃久久精品国产亚洲av| 国语自产精品视频在线第100页| 一a级毛片在线观看| 嫩草影院新地址| 女人被狂操c到高潮| 老司机深夜福利视频在线观看| 午夜影院日韩av| 免费大片18禁| x7x7x7水蜜桃| 久久精品国产亚洲av涩爱 | 久久精品国产亚洲av天美| 欧美绝顶高潮抽搐喷水| 亚洲精品日韩av片在线观看| 国产精品久久视频播放| 看免费av毛片| 最好的美女福利视频网| 亚洲aⅴ乱码一区二区在线播放| 午夜福利成人在线免费观看| 夜夜看夜夜爽夜夜摸| 波多野结衣高清作品| 国产精品嫩草影院av在线观看 | 久久久精品大字幕| 亚洲精品在线观看二区| 久久亚洲精品不卡| 搡老岳熟女国产| 我的女老师完整版在线观看| 欧美精品国产亚洲| 搡老妇女老女人老熟妇| 毛片女人毛片| 亚洲国产精品久久男人天堂| 午夜视频国产福利| 黄色配什么色好看| 99久久久亚洲精品蜜臀av| 久久久成人免费电影| 综合色av麻豆| 日韩高清综合在线| 精品午夜福利视频在线观看一区| 好看av亚洲va欧美ⅴa在| 在线观看av片永久免费下载| 国产私拍福利视频在线观看| 一个人免费在线观看的高清视频| 国产精品久久视频播放| 99热这里只有精品一区| 在线观看舔阴道视频| 少妇高潮的动态图| 国产精品一区二区性色av| 精品国产三级普通话版| 极品教师在线免费播放| 亚洲av一区综合| 欧美性感艳星| 一边摸一边抽搐一进一小说| 51国产日韩欧美| 中文在线观看免费www的网站| a级一级毛片免费在线观看| 看十八女毛片水多多多| 国产乱人视频| 国产免费av片在线观看野外av| 俄罗斯特黄特色一大片| 国产成人aa在线观看| 中文字幕人妻熟人妻熟丝袜美| 一本综合久久免费| 日韩精品青青久久久久久| 国产成年人精品一区二区| 亚洲av熟女| 国产又黄又爽又无遮挡在线| 人人妻人人澡欧美一区二区| 国模一区二区三区四区视频| 看片在线看免费视频| 免费看a级黄色片| 成人午夜高清在线视频| 亚洲18禁久久av| 一进一出抽搐动态| 伦理电影大哥的女人| 日本黄色视频三级网站网址| 欧美性感艳星| 日本三级黄在线观看| 精品一区二区三区人妻视频| 狠狠狠狠99中文字幕| 级片在线观看| 日韩欧美精品v在线| 啪啪无遮挡十八禁网站| 日日夜夜操网爽| 亚洲精品在线观看二区| 午夜免费激情av| 成人欧美大片| 婷婷六月久久综合丁香| 一级作爱视频免费观看| 成人国产综合亚洲| 精品久久国产蜜桃| 亚洲第一欧美日韩一区二区三区| 一区二区三区激情视频| 日本免费a在线| 久久伊人香网站| 日韩欧美 国产精品| 偷拍熟女少妇极品色| 色尼玛亚洲综合影院| 三级国产精品欧美在线观看| 女人十人毛片免费观看3o分钟| 九色国产91popny在线| 亚洲精品亚洲一区二区| 老鸭窝网址在线观看| 精品国内亚洲2022精品成人| 婷婷六月久久综合丁香| 国产视频内射| 国产亚洲精品久久久久久毛片| 久久久久久国产a免费观看| 极品教师在线视频| 午夜福利在线观看免费完整高清在 | 久久久久久大精品| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 岛国在线免费视频观看| 精品久久久久久久人妻蜜臀av| 色噜噜av男人的天堂激情| 久久久久国内视频| 国产精品野战在线观看| 他把我摸到了高潮在线观看| 午夜激情福利司机影院| 亚洲人成网站在线播放欧美日韩| 一个人免费在线观看电影| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 一本一本综合久久| 国内精品一区二区在线观看| 国产精品三级大全| 久久草成人影院| 国内精品美女久久久久久| 18禁在线播放成人免费| 夜夜躁狠狠躁天天躁| av中文乱码字幕在线| 中亚洲国语对白在线视频| 赤兔流量卡办理| 国产主播在线观看一区二区| 99热这里只有是精品在线观看 | 国产成人av教育| 激情在线观看视频在线高清| 极品教师在线免费播放| 99精品在免费线老司机午夜| 他把我摸到了高潮在线观看| xxxwww97欧美| av在线蜜桃| 精品免费久久久久久久清纯| 成人精品一区二区免费| 国产高清视频在线观看网站| 日本三级黄在线观看| 91字幕亚洲| 亚洲aⅴ乱码一区二区在线播放| 99精品在免费线老司机午夜| 在线观看一区二区三区| 99热这里只有是精品在线观看 | 一级黄片播放器| 欧美最新免费一区二区三区 | 十八禁国产超污无遮挡网站| eeuss影院久久| 亚洲欧美日韩高清专用| 黄色一级大片看看| 黄色视频,在线免费观看| 久久人妻av系列| 久久精品91蜜桃| 99国产综合亚洲精品| 亚洲天堂国产精品一区在线| 日本五十路高清| 午夜老司机福利剧场| 国产熟女xx| 色吧在线观看| 亚洲av电影不卡..在线观看| 亚洲真实伦在线观看| 美女cb高潮喷水在线观看| 欧美乱色亚洲激情| 欧美性猛交黑人性爽| av黄色大香蕉| 午夜福利18| 1024手机看黄色片| 亚洲中文字幕日韩| 色精品久久人妻99蜜桃| 又黄又爽又免费观看的视频| 美女xxoo啪啪120秒动态图 | 国产精品爽爽va在线观看网站| 在线播放国产精品三级| 亚洲乱码一区二区免费版| 国产精华一区二区三区| 黄色日韩在线| 欧美性感艳星| 日韩欧美在线乱码| 久久精品影院6| 久久精品国产亚洲av天美| 自拍偷自拍亚洲精品老妇| 国产精品嫩草影院av在线观看 | 国内精品美女久久久久久| 亚洲aⅴ乱码一区二区在线播放| 中文字幕高清在线视频| 欧美精品国产亚洲| 男女那种视频在线观看| 久久精品国产清高在天天线| 久久草成人影院| 成人国产一区最新在线观看| 麻豆国产av国片精品| 欧美最新免费一区二区三区 | 成人欧美大片| 高清日韩中文字幕在线| 蜜桃久久精品国产亚洲av| 最好的美女福利视频网| 一级作爱视频免费观看| 身体一侧抽搐| 亚洲五月婷婷丁香| 国产视频一区二区在线看| 精品免费久久久久久久清纯| 亚洲一区二区三区不卡视频| 久久精品综合一区二区三区| 久久久精品大字幕| 亚洲av美国av| 午夜久久久久精精品| 91av网一区二区| 久99久视频精品免费| 又紧又爽又黄一区二区| 一夜夜www| 九色国产91popny在线| 少妇丰满av| 精品无人区乱码1区二区| 成年女人看的毛片在线观看| 97热精品久久久久久| 成人永久免费在线观看视频| 嫩草影院精品99| 亚洲三级黄色毛片| 亚洲av成人av| 观看免费一级毛片| 他把我摸到了高潮在线观看| 丰满人妻熟妇乱又伦精品不卡| 日韩欧美国产一区二区入口| 18禁黄网站禁片午夜丰满| 亚洲中文字幕日韩| 俺也久久电影网| 亚洲av中文字字幕乱码综合| 亚洲成人中文字幕在线播放| 欧美丝袜亚洲另类 | 老熟妇仑乱视频hdxx| 99riav亚洲国产免费| 亚洲国产精品合色在线| 岛国在线免费视频观看| 国产精品乱码一区二三区的特点| 在线观看美女被高潮喷水网站 | 国产精品一区二区性色av| 琪琪午夜伦伦电影理论片6080| 成人av在线播放网站| 美女被艹到高潮喷水动态| 国产精品影院久久| 制服丝袜大香蕉在线| 国产成人啪精品午夜网站| 不卡一级毛片| 嫩草影院精品99| 精品久久久久久成人av| 999久久久精品免费观看国产| 久久久久久久午夜电影| 亚洲成a人片在线一区二区| 伊人久久精品亚洲午夜| 高清日韩中文字幕在线| 国产中年淑女户外野战色| 欧美+日韩+精品| 丰满的人妻完整版| 亚洲成人久久爱视频| 欧美乱色亚洲激情| 精品乱码久久久久久99久播| aaaaa片日本免费| 夜夜看夜夜爽夜夜摸| 亚洲黑人精品在线| av黄色大香蕉| 草草在线视频免费看| 中出人妻视频一区二区| 日韩欧美精品v在线| 99久久精品国产亚洲精品| 国内少妇人妻偷人精品xxx网站| 高清毛片免费观看视频网站| 亚洲国产高清在线一区二区三| 午夜福利在线观看吧| 欧美激情久久久久久爽电影| 久久6这里有精品| 特级一级黄色大片| 日本a在线网址| 欧美+日韩+精品| 乱人视频在线观看| 欧美激情国产日韩精品一区| 日韩欧美免费精品| 特级一级黄色大片| 成人午夜高清在线视频| 国产乱人伦免费视频| 丝袜美腿在线中文| 亚洲精品一卡2卡三卡4卡5卡| 久久6这里有精品| 亚洲最大成人中文| 成人午夜高清在线视频| 国产乱人伦免费视频| 亚洲av熟女| 757午夜福利合集在线观看| 一级毛片久久久久久久久女| 在线十欧美十亚洲十日本专区| 精品免费久久久久久久清纯| 激情在线观看视频在线高清| 一区二区三区免费毛片| 日本黄色视频三级网站网址| 欧美色视频一区免费| av在线蜜桃| 久久久国产成人免费| 精品一区二区三区人妻视频| 能在线免费观看的黄片| av天堂在线播放| 亚洲人成伊人成综合网2020| 日本与韩国留学比较| 99热这里只有是精品在线观看 | 欧美国产日韩亚洲一区| 久久中文看片网| 又黄又爽又刺激的免费视频.| 亚洲激情在线av| 级片在线观看| 一二三四社区在线视频社区8| 日韩大尺度精品在线看网址| 日本精品一区二区三区蜜桃| 丁香六月欧美| 嫩草影院精品99| 中文资源天堂在线| 亚洲七黄色美女视频| 欧美激情久久久久久爽电影| 久久午夜福利片| 一级黄色大片毛片| 男人的好看免费观看在线视频| 欧美三级亚洲精品| 99在线视频只有这里精品首页| 一进一出抽搐gif免费好疼| 熟女人妻精品中文字幕| 免费看a级黄色片| 久久精品国产99精品国产亚洲性色| 午夜影院日韩av| 美女高潮的动态| 亚洲欧美日韩无卡精品| 日韩欧美精品免费久久 | 国产精品美女特级片免费视频播放器| 波多野结衣高清无吗| 一级黄片播放器| 国产午夜精品久久久久久一区二区三区 | 麻豆久久精品国产亚洲av| 精品一区二区三区视频在线| 亚洲av成人不卡在线观看播放网| 51国产日韩欧美| 婷婷精品国产亚洲av在线| 亚洲成a人片在线一区二区| 亚洲欧美日韩无卡精品| 国产久久久一区二区三区| 男女下面进入的视频免费午夜| 丰满人妻熟妇乱又伦精品不卡| 国产精品,欧美在线| 久久久精品大字幕| 国产精品爽爽va在线观看网站| 国产午夜精品论理片| 日日干狠狠操夜夜爽| 网址你懂的国产日韩在线| 亚洲精华国产精华精| 国产老妇女一区| 男人和女人高潮做爰伦理| 99久久久亚洲精品蜜臀av| 午夜激情欧美在线| www日本黄色视频网| 日本在线视频免费播放| 欧美精品国产亚洲| 91在线观看av| 日本免费a在线| 亚洲色图av天堂| 十八禁人妻一区二区| 国内精品美女久久久久久| 国产麻豆成人av免费视频| 亚洲精品一区av在线观看| 老司机午夜福利在线观看视频| 日本 av在线| 九色国产91popny在线| 亚洲内射少妇av| 精品久久久久久久人妻蜜臀av| 免费大片18禁| 变态另类成人亚洲欧美熟女| 亚洲av免费高清在线观看| 日本免费一区二区三区高清不卡| 国产成人欧美在线观看| 男女做爰动态图高潮gif福利片| 日本在线视频免费播放| 哪里可以看免费的av片| 韩国av一区二区三区四区| 嫩草影院新地址| 在线观看66精品国产| 午夜精品在线福利| 人人妻,人人澡人人爽秒播| 在线看三级毛片| 国产一区二区在线观看日韩| 亚洲精品乱码久久久v下载方式| 草草在线视频免费看| 久久人妻av系列| 亚洲精品粉嫩美女一区| 伊人久久精品亚洲午夜| 男人的好看免费观看在线视频| 精品人妻偷拍中文字幕| 久久精品国产99精品国产亚洲性色| 久久婷婷人人爽人人干人人爱| 淫妇啪啪啪对白视频| 淫秽高清视频在线观看| 日韩精品中文字幕看吧| 一本一本综合久久| 午夜福利18| 欧美国产日韩亚洲一区| 91麻豆av在线| www日本黄色视频网| 中文字幕人妻熟人妻熟丝袜美| 他把我摸到了高潮在线观看| 免费一级毛片在线播放高清视频| 国内毛片毛片毛片毛片毛片| 成年版毛片免费区| av女优亚洲男人天堂| 日韩欧美 国产精品| 国产淫片久久久久久久久 | 一本久久中文字幕| 午夜福利在线在线| 久久精品国产亚洲av天美| 免费高清视频大片| 亚洲午夜理论影院| 久久午夜福利片| 午夜福利高清视频| 亚洲av电影不卡..在线观看| 日韩欧美三级三区| 真人做人爱边吃奶动态| 一个人免费在线观看电影| 欧美日韩综合久久久久久 | 真人一进一出gif抽搐免费| 亚洲 国产 在线| 老司机午夜福利在线观看视频| 99久久精品热视频| 国产激情偷乱视频一区二区| 久久人人爽人人爽人人片va | 亚洲自偷自拍三级| 成人永久免费在线观看视频| 18禁裸乳无遮挡免费网站照片| ponron亚洲| 老司机午夜福利在线观看视频| 国产精品免费一区二区三区在线| 中文字幕高清在线视频| 欧美黑人巨大hd| 国产精品综合久久久久久久免费| 国产视频内射| 日韩欧美一区二区三区在线观看| 小蜜桃在线观看免费完整版高清| 日韩欧美免费精品| 欧美最新免费一区二区三区 | 亚洲专区国产一区二区| 三级毛片av免费| 国产不卡一卡二| 国产又黄又爽又无遮挡在线| 亚洲狠狠婷婷综合久久图片| 国产乱人伦免费视频| 禁无遮挡网站| 日韩成人在线观看一区二区三区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 精品久久久久久久久久久久久| 亚洲av日韩精品久久久久久密| 国内精品一区二区在线观看| 日韩欧美精品免费久久 | 日本三级黄在线观看| 久久久久久国产a免费观看| 婷婷亚洲欧美| 欧美色视频一区免费| 一级a爱片免费观看的视频| 亚洲在线自拍视频| 69人妻影院| 我的老师免费观看完整版| 男人狂女人下面高潮的视频| 亚洲精品色激情综合| 亚洲一区二区三区不卡视频| 青草久久国产| 一级a爱片免费观看的视频| 午夜精品一区二区三区免费看| 性色av乱码一区二区三区2| 精品不卡国产一区二区三区| 亚洲成a人片在线一区二区| 国产精品免费一区二区三区在线| 国产伦一二天堂av在线观看| 给我免费播放毛片高清在线观看| 免费看a级黄色片| 久久人妻av系列| 少妇裸体淫交视频免费看高清| 999久久久精品免费观看国产| 亚洲精品456在线播放app | av女优亚洲男人天堂| 少妇丰满av| 亚洲av中文字字幕乱码综合| 久久午夜福利片| 亚洲真实伦在线观看| 九色国产91popny在线| 在线观看美女被高潮喷水网站 | 日本 欧美在线| 国产精品久久久久久亚洲av鲁大| 天堂影院成人在线观看| 我要看日韩黄色一级片| 精品国产亚洲在线| 精品午夜福利视频在线观看一区| 有码 亚洲区| 波野结衣二区三区在线| 国内精品久久久久久久电影| 亚洲,欧美精品.| 日本熟妇午夜| 国产真实乱freesex| 国产主播在线观看一区二区| 99久久99久久久精品蜜桃| 日韩高清综合在线| 亚洲国产精品sss在线观看| 脱女人内裤的视频| 丝袜美腿在线中文| 亚洲av成人不卡在线观看播放网| 乱码一卡2卡4卡精品| 免费观看人在逋| 一区二区三区高清视频在线| 亚洲avbb在线观看| 色综合站精品国产| 久久香蕉精品热| 国产精品伦人一区二区| 国产国拍精品亚洲av在线观看| 黄色一级大片看看| 精品无人区乱码1区二区| 丰满人妻一区二区三区视频av| 99国产精品一区二区三区| 久9热在线精品视频| 国产伦精品一区二区三区视频9| 欧美日本视频| 真人一进一出gif抽搐免费| 亚洲片人在线观看| 亚洲成av人片在线播放无| 俺也久久电影网| 欧美在线黄色| 人人妻人人看人人澡| 免费观看的影片在线观看| 国产一区二区在线观看日韩| 午夜日韩欧美国产| 一本久久中文字幕| eeuss影院久久| 免费av观看视频| 毛片女人毛片| 日本 av在线| 国产午夜精品论理片| 午夜精品久久久久久毛片777| 久久久成人免费电影| 国产成+人综合+亚洲专区| 成人午夜高清在线视频| 成年女人看的毛片在线观看| 两人在一起打扑克的视频|