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

    高速直齒輪齒面有限元網(wǎng)格精準(zhǔn)離散方法

    2023-04-11 12:36:22黃一倫胡玉梅
    關(guān)鍵詞:切向速度輪齒齒面

    黃一倫,陳 旭,胡玉梅

    (1.重慶理工大學(xué) 車輛工程學(xué)院, 重慶 400054;2.重慶大學(xué) 機械傳動國家重點實驗室, 重慶 400044)

    有限元作為一種常用的齒輪嚙合仿真分析方法,被廣泛應(yīng)用于齒輪動力學(xué)性能的研究,其中齒輪齒面在進(jìn)行有限元網(wǎng)格離散時,網(wǎng)格的大小對于齒輪嚙合有限元仿真結(jié)果的正確性影響很大。在過去幾十年間,國內(nèi)外許多學(xué)者針對齒輪有限元網(wǎng)格離散的問題展開了研究。劉更等[1-2]通過程序編寫出不同的齒廓網(wǎng)格節(jié)點生成函數(shù),根據(jù)計算得到的節(jié)點坐標(biāo)自動生成輪齒有限元網(wǎng)格。穆慧勇等[3]將斜齒輪齒根中點拉應(yīng)力值取得平衡時的初始網(wǎng)格密度作為齒面齒向網(wǎng)格密度。張濤等[4]通過ANSYS的參數(shù)化程序語言,利用8節(jié)點六面體單元對齒輪進(jìn)行有限元離散。曹雪梅等[5]利用刨齒展成法精確計算出齒面節(jié)點坐標(biāo),在旋轉(zhuǎn)投影平面上劃分網(wǎng)格,最終映射到齒面上。李紅勛等[6]綜合計算時長、精度和收斂性等因素,將齒面網(wǎng)格細(xì)化,而輪轂基體網(wǎng)格設(shè)置較稀疏。Coy等[7]提出了一種考慮赫茲變形影響的齒面網(wǎng)格尺寸選擇方法。Vijayakar等[8]通過結(jié)合曲面積分法和傳統(tǒng)有限元分析方法研究了齒輪齒面的有限元網(wǎng)格離散問題。Gonzalez-Perez等[9-11]研究了考慮齒輪扭轉(zhuǎn)變形和軸偏轉(zhuǎn)影響時的有限元網(wǎng)格劃分方法。Gonzalez-Perez等[12]提出了一種基于多點約束的齒輪傳動有限元模型網(wǎng)格細(xì)化策略。

    綜上所述,國內(nèi)外學(xué)者對齒輪有限元網(wǎng)格離散研究已經(jīng)取得了豐碩成果,然而這些研究并未考慮齒面有限元網(wǎng)格尺寸在高轉(zhuǎn)速下的適應(yīng)性問題。隨著現(xiàn)代航空航天技術(shù)的不斷發(fā)展,高線速直齒圓柱齒輪作為航空發(fā)動機動力傳輸?shù)年P(guān)鍵零部件具有質(zhì)量輕、轉(zhuǎn)速高和載荷大等特性,其服役過程呈現(xiàn)強烈的非線性和瞬態(tài)性,所以對該類齒輪齒面有限元網(wǎng)格離散化的精度提出了更高要求。為此,本課題組用兩圓柱滾子曲面接觸有限元模型來模擬齒面接觸,通過有限元仿真分析不同轉(zhuǎn)速下兩圓柱滾子接觸產(chǎn)生的滑移能損失,提出一種齒面網(wǎng)格大小能夠適應(yīng)齒面滾滑速度以及曲率半徑變化的齒面有限元網(wǎng)格精準(zhǔn)離散方法。

    1 齒面嚙合簡化模型

    在齒輪嚙合過程中,輪齒之間的接觸屬于2個曲面之間接觸,而曲面有限元離散的基本思想是用有限個單元去近似模擬連續(xù)光滑的曲面體,要達(dá)到同等的離散化精度,曲率半徑越小的曲面,所需要的單元網(wǎng)格尺寸也越小。此外,高線速直齒輪轉(zhuǎn)速較高,齒面有限元網(wǎng)格存在較大的離心力,導(dǎo)致齒面有限元模型在參與嚙合時會產(chǎn)生單元網(wǎng)格之間的沖擊振動,這與實際的齒面接觸存在一定誤差。所以要想建立精確的齒面接觸有限元模型,關(guān)鍵問題是找到在不同轉(zhuǎn)速下齒面有限元網(wǎng)格尺寸與齒面曲率半徑之間的最佳倍數(shù)關(guān)系。

    滑移能損失是2個相對滑動的面碰撞所導(dǎo)致的能量交換,不可避免。曲面接觸的滑移能損失與摩擦系數(shù)及正壓力有關(guān),理論滑移能損失E可以表示為:

    E=f×F×Vt×t

    (1)

    式中:f為接觸摩擦系數(shù);F為垂直載荷,N;Vt為兩接觸面的相對滑動速度,m/s;t為時間,s。齒輪嚙合過程中的能量損失主要來源于輪齒之間的滑移能損失,所以輪齒之間的滑移能損失在一定程度上能夠反映出輪齒之間的接觸狀態(tài)。由于齒面上不同嚙合位置處的齒面壓力以及相對滑移速度均不相同,所以理論求解齒面滑移能損失存在一定難度。直齒輪的漸開線齒廓本質(zhì)上是由一系列曲率半徑不同的圓弧所組成的光滑曲線,并且齒面接觸區(qū)寬度遠(yuǎn)遠(yuǎn)小于接觸區(qū)域的曲率半徑,因而可以對嚙合齒面作適當(dāng)簡化。根據(jù)Weck等[13]試驗結(jié)果表明,在相同工作條件下,嚙合輪齒之間的接觸狀態(tài)可以簡化為2個圓柱滾子沿其母線方向上的接觸,如圖1所示。圖1中圓柱體滾子的半徑R1、R2分別代表齒輪嚙合位置處的齒面曲率半徑;n1和n2為圓柱滾子的轉(zhuǎn)速,上下圓柱滾子接觸點處的線速度分別代表主從輪齒齒面在嚙合位置處的切向速度;圓柱滾子垂直向下的壓力F代表齒面嚙合位置處的法向載荷;b代表齒寬。

    圖1 圓柱滾子曲面接觸模型簡圖

    2 直齒輪嚙合幾何分析

    圓柱滾子曲面接觸模型能精準(zhǔn)模擬高線速直齒圓柱齒輪齒面嚙合的前提是對其施加準(zhǔn)確的邊界條件以及設(shè)置合理的幾何參數(shù),其中包括圓柱滾子的曲率半徑、轉(zhuǎn)速以及垂直載荷,而這些邊界條件和幾何參數(shù)的值與高線速直齒圓柱齒輪的齒面曲率半徑、齒面切向速度、相對滑移速度以及齒面法向載荷的大小密切相關(guān)。本課題組根據(jù)直齒圓柱齒輪漸開線齒廓的嚙合過程,通過理論推導(dǎo)得出這些參數(shù)的計算方法。

    圖2為直齒輪齒廓嚙合過程示意圖。P為節(jié)點,K1點和K2點是嚙合線上任意位置嚙合點。N1N2為理論嚙合線,B1B2為實際嚙合線,主動輪從B1點進(jìn)入嚙合并推動從動輪的齒頂旋轉(zhuǎn),從動輪在B2點與主動輪齒頂脫離嚙合。

    圖2 齒輪齒廓嚙合過程簡圖

    根據(jù)齒廓嚙合過程中的幾何關(guān)系可知,齒面上任意嚙合點的曲率半徑等于齒面漸開線的發(fā)生線長度,即圖2中嚙合點距離N1點或N2點的長度,所以任意嚙合點在主、從動輪上的曲率半徑可以表示為:

    (2)

    (3)

    式中:ρ1為任意嚙合點在主動輪上的曲率半徑;ρ2為任意嚙合點在從動輪上的曲率半徑;α為分度圓壓力角;r1為主動輪分度圓半徑;r2為從動輪分度圓半徑;rk為任意嚙合點在主動輪上的半徑;rb1為主動輪基圓半徑。

    圖2中的V1k、V1t、V1n分別是嚙合點K2在主動輪上的絕對速度、切向速度以及法向速度,根據(jù)任意嚙合點K2在主動輪上的速度矢量三角形與直角三角形⊿K2O1N1的相似關(guān)系,主動輪上任意嚙合點沿齒面的切向速度V1t可以表示為:

    (4)

    式中:n1為主動輪轉(zhuǎn)速,r/min;ρ1為任意嚙合點K2在主動輪上的曲率半徑,mm;V1t為主動輪上任意嚙合點K2沿齒面的切向速度,m/s。

    圖2中的V2k、V2t、V2n分別是嚙合點K2在從動輪上的絕對速度、切向速度以及法向速度,根據(jù)任意嚙合點K2在從動輪上的速度矢量三角形與直角三角形⊿K2O2N2的相似關(guān)系,從動輪上任意嚙合點沿齒面的切向速度V2t可以表示為:

    (5)

    式中:n2為從動輪轉(zhuǎn)速,r/min;ρ2為任意嚙合點K2在從動輪上的曲率半徑,mm;V2t為從動輪上任意嚙合點K2沿齒面的切向速度,m/s。

    主、從動輪任意嚙合點處的相對滑移速度Vt可以表示為:

    Vt=V1t-V2t

    (6)

    利用式(2)(3)計算出主從動輪齒面的最大最小曲率半徑大小,并將其作為圓柱滾子曲率半徑的確定依據(jù)。利用式(4)—(6)計算出最大齒面切向速度和相對滑移速度作為圓柱滾子模型中滾子轉(zhuǎn)速的確定依據(jù)。

    此外,圓柱滾子的垂直載荷大小應(yīng)當(dāng)與航空高速重載直齒圓柱齒輪的齒面法向載荷相對應(yīng)。直齒圓柱齒輪齒輪從動輪齒受力分析如圖3所示,其中d2為從動輪分度圓直徑,T為扭矩。為了計算圖3中的齒面名義法向載荷Fn,將其在分度圓處分解為圓周力Ft和徑向力Fr,根據(jù)力平衡條件和各力之間的幾何關(guān)系,齒面名義法向載荷Fn可以表示為[14]:

    (7)

    式中:rb2為從動輪基圓半徑;P為從動輪輸出功率,kW。

    圖3 直齒圓柱齒輪輪齒受力分析

    根據(jù)式(7)可以求出齒面法向載荷的大小,并將其作為圓柱滾子模型中滾子垂直載荷的確定依據(jù)。

    3 建模與仿真分析

    以某在研航空發(fā)動機齒輪箱為研究對象,根據(jù)式(2)—(7)分別計算得出齒輪箱中所有直齒圓柱齒輪的齒面曲率半徑、齒面切向速度、齒面相對滑移速度以及齒面載荷,并最終總結(jié)得出的航空發(fā)動機高速重載齒輪齒面曲率半徑、齒面切向速度、齒面相對滑移速度以及齒面載荷的大致范圍如表1所示。基于此范圍確定圓柱滾子的曲率、轉(zhuǎn)速以及載荷等參數(shù)大小。

    由于相對滑移速度的大小隨著兩接觸面的切向速度變化而變化,所以為了定量分析齒面切向速度對主從輪齒接觸的影響,保持圓柱滾子接觸面之間的相對滑移速度Vt為20 m/s。直齒圓柱齒輪理想潤滑條件下的齒面摩擦系數(shù)為0.01~0.02[15],所以滾子之間的接觸摩擦系數(shù)固定為0.01。

    表1 航空發(fā)動機高速重載齒輪齒面曲率、速度和載荷等參數(shù)范圍

    分別建立4組接觸面切向速度不同的圓柱滾子曲面接觸模型,具體參數(shù)如表2所示。

    表2 圓柱滾子曲面接觸模型參數(shù)

    圖4為圓柱滾子曲面接觸有限元模型。分別對上下滾子施加大小不同的轉(zhuǎn)速,并在上滾子軸心線的節(jié)點上施加豎直向下的均布載荷,最后在上下滾子的接觸面之間建立接觸。

    圖4 圓柱滾子曲面接觸有限元模型

    將模型導(dǎo)入ANSYS軟件中使用LS-DYNA求解器仿真計算。對于線速度大于20 m/s的航空發(fā)動機直齒圓柱齒輪,單對輪齒的嚙合時間往往不超過0.001 s,所以將兩圓柱滾子在0.001 s時間內(nèi)接觸產(chǎn)生的滑移能仿真值與理論值之間的百分比誤差小于10%作為前提,對比分析得出在不同切向速度下曲面網(wǎng)格尺寸大小與曲率半徑之間的最佳倍數(shù)關(guān)系。

    圖5分別為表2中4組圓柱滾子曲面接觸模型在0.000 2、0.000 4、0.000 6、0.000 8、0.001 s時的滑移能損失誤差曲線。

    圖5 不同切向速度下的圓柱滾子滑移能誤差曲線

    根據(jù)以上仿真結(jié)果可知,隨著網(wǎng)格的細(xì)化,滑移能誤差百分比值呈現(xiàn)逐漸減小的趨勢,而隨著時間的推移,滑移能仿真值與理論值之間百分比誤差維持在小范圍內(nèi)波動,基本保持恒定。同時,在相同的網(wǎng)格細(xì)化程度下,滑移能誤差百分比與2個滾子接觸曲面之間最大切向速度的絕對值成正比關(guān)系,最大切向速度越大,則滑移能誤差百分比值越大。對于最大切向速度為25 m/s的滾子模型,接觸面網(wǎng)格尺寸為R/40~R/50時的滑移能百分比誤差相對較小,滿足建模精度要求。最大切向速度為40 m/s時,接觸面網(wǎng)格尺寸為R/70~R/80時的滑移能百分比誤差小于10%。最大切向速度為55 m/s時,接觸面網(wǎng)格的細(xì)化程度最少要達(dá)到R/110才能達(dá)到建模精度要求。而當(dāng)最大切向速度為70 m/s的滾子模型,接觸面網(wǎng)格的細(xì)化程度至少需要要達(dá)到R/140,才滿足誤差小于10%的工程精度要求。綜合誤差、齒輪齒面參數(shù)以及模型計算時間等多方面因素,總結(jié)出滿足工程建模精度的高速齒輪齒面網(wǎng)格細(xì)化準(zhǔn)則,如表3所示。

    表3 直齒圓柱齒輪齒面網(wǎng)格尺寸確定標(biāo)準(zhǔn)

    4 算例與驗證

    為了證明上述方法的有效性,以某在研航空發(fā)動機齒輪箱中的一對直齒圓柱齒輪副為例,根據(jù)表3中給出的齒面網(wǎng)格離散準(zhǔn)則劃分齒輪齒面有限元網(wǎng)格,并分別建立4組不同轉(zhuǎn)速工況下的齒輪嚙合有限元模型,對其進(jìn)行仿真分析。齒輪基本參數(shù)如表4所示。

    表4 航空發(fā)動機齒輪箱中某直齒圓柱齒輪副參數(shù)

    利用式(2)(3)可以計算出主從輪齒齒面的最小曲率半徑分別為24.7 mm和17.2 mm,利用式(4)—(6)計算出4種轉(zhuǎn)速下齒輪副輪齒嚙入到嚙出過程中的最大齒面切向速度,并根據(jù)表3確定出該齒輪副在4種轉(zhuǎn)速下的最佳齒面有限元網(wǎng)格尺寸,如表5所示。

    以表5中的第一組工況為例,說明齒輪副嚙合有限元模型的建模方法(其余工況除轉(zhuǎn)速、扭矩及齒面網(wǎng)格尺寸不同外,建模流程均相同)。圖6為主從輪齒齒面有限元模型。其中齒面體網(wǎng)格均使用六面體網(wǎng)格,齒廓線上的節(jié)點間距以及齒向方向的網(wǎng)格節(jié)點間距均等于表5中計算得出的齒面網(wǎng)格尺寸。

    表5 不同工況下的最佳齒面有限元網(wǎng)格尺寸

    圖6 主從輪齒齒面有限元模型

    如圖7所示,分別以對主從動輪施加轉(zhuǎn)速、轉(zhuǎn)矩及接觸為邊界條件,建立齒輪副嚙合有限元模型。最終導(dǎo)入ANSYS中利用LS-DYNA求解器進(jìn)行仿真求解。

    分別提取齒輪接觸力穩(wěn)定階段時4種不同轉(zhuǎn)速下節(jié)點位置處的齒面接觸應(yīng)力,如圖8所示。根據(jù)仿真結(jié)果可知,主從輪齒齒面接觸區(qū)的應(yīng)力呈較為理想的橢圓形分布,表5中工況1、工況2、工況3和工況4所對應(yīng)的最大齒面接觸應(yīng)力分別為608.6、545.2、597.2和666.5 MPa。

    圖7 齒輪副嚙合有限元模型

    圖8 4種工況下的齒面接觸壓應(yīng)力

    參考國家標(biāo)準(zhǔn)[16],齒輪使用系數(shù)、齒間載荷分配系數(shù)及齒向載荷分布系數(shù)均取1,而表5中的工況1、工況2、工況3和工況4所對應(yīng)的動載系數(shù)分別取1.65、1.44、1.21和1.06,根據(jù)赫茲公式計算得出的最大齒面接觸應(yīng)力理論值分別為557.8、597.3、662.9和706.5 MPa,與仿真值之間的誤差分別為9.1%、8.7%、9.9%和5.6%。說明上述高速直齒圓柱齒輪齒面有限元網(wǎng)格離散方法能保證不同轉(zhuǎn)速下的齒輪有限元仿真誤差小于10%。

    同理,相比于表5中的齒面網(wǎng)格尺寸,在工況1的條件下,分別采用一組更大(主從動輪齒面網(wǎng)格尺寸為0.35 mm)和一組更小(主從動輪齒面網(wǎng)格尺寸分別為0.095 mm)的網(wǎng)格尺寸對齒面進(jìn)行離散,并采用相同的建模方法建立齒輪嚙合有限元模型,并對其仿真分析。分別提取2種不同齒面網(wǎng)格尺寸下節(jié)點位置處的齒面接觸壓應(yīng)力,如圖9所示。根據(jù)仿真結(jié)果可知,大網(wǎng)格條件下所對應(yīng)的最大齒面接觸應(yīng)力為340.5 MPa,與理論值之間的誤差為39%,仿真精度明顯變差;而小網(wǎng)格條件下所對應(yīng)的最大齒面接觸應(yīng)力為541.3 MPa,與理論值之間的誤差為3%,仿真精度有所提高。如圖10所示,對比3種網(wǎng)格尺寸下的仿真誤差及仿真時間可知,齒面網(wǎng)格尺寸越小,有限元仿真精度越高,但同時仿真時間也會大幅延長。說明采用表3所推薦的齒面網(wǎng)格尺寸,既保證了齒輪有限元仿真的精度,也將仿真時長控制在了可接受范圍之內(nèi),證明上述高線速直齒圓柱齒輪齒面有限元網(wǎng)格離散方法是正確有效的。

    圖9 不同齒面網(wǎng)格尺寸下的齒面接觸壓應(yīng)力

    圖10 不同齒面網(wǎng)格尺寸下的齒面壓應(yīng)力誤差及仿真時長

    5 結(jié)論

    1) 當(dāng)最大齒面切向速度在0~25 m/s時,齒面網(wǎng)格尺寸應(yīng)當(dāng)取齒面最小曲率半徑的1/40。

    2) 當(dāng)最大齒面切向速度在25~40 m/s時,齒面網(wǎng)格尺寸應(yīng)當(dāng)取齒面最小曲率半徑的1/70。

    3) 當(dāng)最大齒面切向速度在40~55 m/s時,齒面網(wǎng)格尺寸應(yīng)當(dāng)取齒面最小曲率半徑的1/110。

    4) 當(dāng)最大齒面切向速度在55~70 m/s時,齒面網(wǎng)格尺寸應(yīng)當(dāng)取齒面最小曲率半徑的1/140。

    5) 采用上述齒面有限元網(wǎng)格離散方法所得到的有限元模型,其仿真結(jié)果與理論值誤差小于10%。

    猜你喜歡
    切向速度輪齒齒面
    曳引式電梯蝸輪輪齒失效原因分析及預(yù)防措施
    旋風(fēng)分離器內(nèi)氣相旋轉(zhuǎn)流不穩(wěn)定性的實驗研究
    基于NURBS理論的漸開線齒輪齒面修復(fù)
    重型機械(2020年3期)2020-08-24 08:31:46
    雙旋流耦合式旋流反應(yīng)器內(nèi)切向速度分布研究
    流體機械(2020年4期)2020-05-12 09:20:48
    齒輪發(fā)生隨機斷裂的原因和預(yù)防措施
    電梯蝸輪輪齒失效原因分析及警示
    基于BP神經(jīng)網(wǎng)絡(luò)的面齒輪齒面粗糙度研究
    高速動車組弧齒錐齒輪齒面疲勞點蝕失效分析
    風(fēng)力發(fā)電機齒輪箱輪齒斷裂原因分析
    重型機械(2016年1期)2016-03-01 03:42:12
    17CrNiMo6齒輪軸齒面剝落原因分析
    大型鑄鍛件(2015年1期)2016-01-12 06:33:01
    中文字幕高清在线视频| 亚洲精品美女久久av网站| 亚洲 国产 在线| 九色亚洲精品在线播放| av一本久久久久| 超色免费av| 精品国产一区二区三区久久久樱花| 男女午夜视频在线观看| 亚洲精品第二区| 美女国产高潮福利片在线看| 亚洲国产成人一精品久久久| 黄片小视频在线播放| 亚洲性夜色夜夜综合| 在线观看免费日韩欧美大片| √禁漫天堂资源中文www| 在线观看舔阴道视频| 动漫黄色视频在线观看| 蜜桃在线观看..| 亚洲性夜色夜夜综合| 成人黄色视频免费在线看| www.精华液| 国产在视频线精品| 黑人巨大精品欧美一区二区蜜桃| 一本一本久久a久久精品综合妖精| 满18在线观看网站| www.精华液| 一本大道久久a久久精品| 黄色 视频免费看| 欧美97在线视频| 90打野战视频偷拍视频| 亚洲一卡2卡3卡4卡5卡精品中文| 我的亚洲天堂| 狂野欧美激情性bbbbbb| 日本91视频免费播放| 美女高潮到喷水免费观看| 桃红色精品国产亚洲av| svipshipincom国产片| 天天躁日日躁夜夜躁夜夜| 日韩中文字幕欧美一区二区| 国产一卡二卡三卡精品| 国产黄频视频在线观看| 国产男女内射视频| 岛国在线观看网站| 国产欧美日韩一区二区三区在线| 日本五十路高清| 久久av网站| 国产成人av激情在线播放| 国产精品一区二区免费欧美 | 午夜免费成人在线视频| 成人免费观看视频高清| 99久久国产精品久久久| 99国产精品99久久久久| 好男人电影高清在线观看| 国产精品二区激情视频| 青草久久国产| 国产高清videossex| 超碰成人久久| 国产成人系列免费观看| 欧美xxⅹ黑人| 伊人亚洲综合成人网| 老司机影院成人| 亚洲avbb在线观看| 免费高清在线观看视频在线观看| 最近最新中文字幕大全免费视频| 女人久久www免费人成看片| 一二三四社区在线视频社区8| 欧美黑人欧美精品刺激| 亚洲精品美女久久av网站| 亚洲av电影在线进入| 99精品欧美一区二区三区四区| av在线播放精品| 天堂8中文在线网| 两个人看的免费小视频| 两人在一起打扑克的视频| 桃红色精品国产亚洲av| a在线观看视频网站| 如日韩欧美国产精品一区二区三区| 母亲3免费完整高清在线观看| 久久精品国产亚洲av高清一级| 免费观看人在逋| 男女免费视频国产| 精品国产乱码久久久久久男人| 人人妻人人澡人人爽人人夜夜| 香蕉丝袜av| 最近最新免费中文字幕在线| 亚洲精品日韩在线中文字幕| bbb黄色大片| 可以免费在线观看a视频的电影网站| 国产精品自产拍在线观看55亚洲 | 亚洲精品国产av蜜桃| 丰满饥渴人妻一区二区三| 国产精品一二三区在线看| 无遮挡黄片免费观看| 少妇被粗大的猛进出69影院| a级片在线免费高清观看视频| 大香蕉久久网| 人人妻人人添人人爽欧美一区卜| 菩萨蛮人人尽说江南好唐韦庄| 亚洲九九香蕉| 亚洲第一青青草原| av在线老鸭窝| 成人国产一区最新在线观看| av电影中文网址| 老司机午夜十八禁免费视频| 丝袜人妻中文字幕| 精品国产一区二区三区久久久樱花| 久久久欧美国产精品| 久久久久久久国产电影| www日本在线高清视频| 人人妻人人添人人爽欧美一区卜| 亚洲avbb在线观看| tocl精华| 男女下面插进去视频免费观看| av有码第一页| av在线老鸭窝| 热re99久久精品国产66热6| 亚洲综合色网址| 中文欧美无线码| 亚洲精品国产一区二区精华液| 久久热在线av| 狠狠狠狠99中文字幕| 免费观看人在逋| 交换朋友夫妻互换小说| 一级片免费观看大全| 在线观看免费日韩欧美大片| 亚洲人成电影免费在线| 亚洲av男天堂| 免费不卡黄色视频| 精品一区二区三区四区五区乱码| 精品国产一区二区久久| 精品久久久久久电影网| 国产区一区二久久| 一级片'在线观看视频| 人人澡人人妻人| 91精品国产国语对白视频| 欧美日韩国产mv在线观看视频| 亚洲一区中文字幕在线| 国产精品一二三区在线看| 久久久久久久大尺度免费视频| videos熟女内射| 亚洲av电影在线进入| 精品久久久久久电影网| 久久精品国产综合久久久| av线在线观看网站| 国产成人精品久久二区二区91| 亚洲国产av影院在线观看| av一本久久久久| 日日摸夜夜添夜夜添小说| 免费一级毛片在线播放高清视频 | 美国免费a级毛片| 国产精品免费视频内射| 午夜两性在线视频| 亚洲国产看品久久| 免费在线观看完整版高清| 国产精品自产拍在线观看55亚洲 | 制服人妻中文乱码| av一本久久久久| 亚洲成人手机| 亚洲午夜精品一区,二区,三区| 国产精品一二三区在线看| 97在线人人人人妻| 另类亚洲欧美激情| 国内毛片毛片毛片毛片毛片| 69精品国产乱码久久久| 欧美中文综合在线视频| 精品高清国产在线一区| 亚洲午夜精品一区,二区,三区| 亚洲欧洲精品一区二区精品久久久| 亚洲国产欧美在线一区| 桃红色精品国产亚洲av| 国产男女超爽视频在线观看| 久久久国产欧美日韩av| 老司机亚洲免费影院| 精品第一国产精品| 久久久精品国产亚洲av高清涩受| 最近中文字幕2019免费版| 亚洲国产成人一精品久久久| a 毛片基地| 啦啦啦在线免费观看视频4| 大码成人一级视频| 如日韩欧美国产精品一区二区三区| 欧美日本中文国产一区发布| 80岁老熟妇乱子伦牲交| 亚洲国产毛片av蜜桃av| 国产亚洲精品第一综合不卡| 无遮挡黄片免费观看| 亚洲成人免费av在线播放| 亚洲精品国产av蜜桃| 色婷婷久久久亚洲欧美| 曰老女人黄片| 亚洲熟女毛片儿| 亚洲精品中文字幕在线视频| 丝袜喷水一区| 久久毛片免费看一区二区三区| 岛国毛片在线播放| 中文字幕另类日韩欧美亚洲嫩草| 久久久久精品人妻al黑| 亚洲国产日韩一区二区| 在线精品无人区一区二区三| 女人久久www免费人成看片| 在线看a的网站| 日韩人妻精品一区2区三区| 亚洲av欧美aⅴ国产| 在线看a的网站| 日本av免费视频播放| 黄片大片在线免费观看| tocl精华| 大香蕉久久成人网| 亚洲一卡2卡3卡4卡5卡精品中文| 国产国语露脸激情在线看| 69精品国产乱码久久久| 别揉我奶头~嗯~啊~动态视频 | 午夜老司机福利片| 成人国产av品久久久| 王馨瑶露胸无遮挡在线观看| 免费少妇av软件| 免费在线观看黄色视频的| 久久国产精品人妻蜜桃| 每晚都被弄得嗷嗷叫到高潮| 1024香蕉在线观看| 高清视频免费观看一区二区| 美女脱内裤让男人舔精品视频| 精品国产乱码久久久久久男人| 日韩大码丰满熟妇| 男女边摸边吃奶| 黄频高清免费视频| 青草久久国产| 国产97色在线日韩免费| 淫妇啪啪啪对白视频 | 国产成人精品无人区| 男人操女人黄网站| 欧美少妇被猛烈插入视频| 色视频在线一区二区三区| 精品福利观看| 午夜福利影视在线免费观看| 久久久国产成人免费| 亚洲精品国产一区二区精华液| 一级毛片女人18水好多| 日韩 亚洲 欧美在线| 午夜福利乱码中文字幕| 国产精品av久久久久免费| 国产色视频综合| 国产成人啪精品午夜网站| 欧美乱码精品一区二区三区| 国产精品免费视频内射| 亚洲少妇的诱惑av| 久久av网站| 中国美女看黄片| 一本色道久久久久久精品综合| 国产精品一区二区在线不卡| 久久国产亚洲av麻豆专区| 另类精品久久| 天天影视国产精品| 欧美黑人欧美精品刺激| 满18在线观看网站| 亚洲国产av影院在线观看| 9色porny在线观看| 一二三四在线观看免费中文在| 欧美精品高潮呻吟av久久| 五月开心婷婷网| 最近最新中文字幕大全免费视频| 久久久久久亚洲精品国产蜜桃av| av在线播放精品| 亚洲三区欧美一区| 一区二区日韩欧美中文字幕| 欧美激情高清一区二区三区| 亚洲欧美成人综合另类久久久| 无限看片的www在线观看| 国产伦理片在线播放av一区| 人妻久久中文字幕网| 热re99久久精品国产66热6| 久久人妻福利社区极品人妻图片| av国产精品久久久久影院| 不卡一级毛片| 亚洲欧美精品自产自拍| 久久青草综合色| 在线亚洲精品国产二区图片欧美| 美女脱内裤让男人舔精品视频| 国产男女内射视频| 成人国产一区最新在线观看| 国产有黄有色有爽视频| 欧美日韩视频精品一区| 日韩一卡2卡3卡4卡2021年| 国产精品欧美亚洲77777| 法律面前人人平等表现在哪些方面 | 在线 av 中文字幕| 黄色 视频免费看| 日韩欧美一区视频在线观看| 天天添夜夜摸| 十八禁人妻一区二区| 国产男女内射视频| 91九色精品人成在线观看| 欧美亚洲日本最大视频资源| 十八禁网站网址无遮挡| 最新在线观看一区二区三区| 丝袜脚勾引网站| 美女高潮到喷水免费观看| 亚洲精品国产一区二区精华液| 亚洲精品国产区一区二| 这个男人来自地球电影免费观看| 不卡av一区二区三区| 欧美黑人欧美精品刺激| 老司机午夜福利在线观看视频 | 99九九在线精品视频| 午夜福利乱码中文字幕| 欧美日韩亚洲综合一区二区三区_| 午夜激情久久久久久久| 婷婷成人精品国产| 交换朋友夫妻互换小说| 91九色精品人成在线观看| 精品乱码久久久久久99久播| 18禁裸乳无遮挡动漫免费视频| 蜜桃国产av成人99| cao死你这个sao货| 久久青草综合色| 久热这里只有精品99| 日韩制服骚丝袜av| 12—13女人毛片做爰片一| 大陆偷拍与自拍| 亚洲精品国产色婷婷电影| 极品少妇高潮喷水抽搐| 精品国产乱码久久久久久小说| 九色亚洲精品在线播放| 人成视频在线观看免费观看| 我要看黄色一级片免费的| www.999成人在线观看| 午夜激情久久久久久久| 国产一区有黄有色的免费视频| 视频区图区小说| 亚洲专区国产一区二区| 最近最新中文字幕大全免费视频| 一本综合久久免费| 叶爱在线成人免费视频播放| 91精品国产国语对白视频| 精品乱码久久久久久99久播| av一本久久久久| 美女午夜性视频免费| 国产成人av教育| 久久久久久免费高清国产稀缺| 国产福利在线免费观看视频| 老司机深夜福利视频在线观看 | 不卡av一区二区三区| 久久久国产精品麻豆| 黄片播放在线免费| 99国产精品99久久久久| 女性生殖器流出的白浆| 黄色毛片三级朝国网站| 日本黄色日本黄色录像| 久久热在线av| 欧美黑人欧美精品刺激| 久久热在线av| 欧美激情 高清一区二区三区| 国产精品欧美亚洲77777| 9191精品国产免费久久| 亚洲精品一二三| av又黄又爽大尺度在线免费看| 乱人伦中国视频| 美女主播在线视频| 悠悠久久av| 亚洲av国产av综合av卡| 一二三四在线观看免费中文在| 国产福利在线免费观看视频| 一边摸一边抽搐一进一出视频| 成在线人永久免费视频| 亚洲av成人一区二区三| 午夜福利免费观看在线| 啦啦啦视频在线资源免费观看| 91精品三级在线观看| 久久ye,这里只有精品| 国产淫语在线视频| 午夜久久久在线观看| 国产亚洲av片在线观看秒播厂| 国产国语露脸激情在线看| 999久久久精品免费观看国产| 1024视频免费在线观看| 少妇 在线观看| 亚洲欧美精品自产自拍| 欧美精品啪啪一区二区三区 | 国产日韩一区二区三区精品不卡| 日韩欧美一区视频在线观看| 国产激情久久老熟女| 国产亚洲av片在线观看秒播厂| 日本撒尿小便嘘嘘汇集6| 女人被躁到高潮嗷嗷叫费观| 中文字幕人妻熟女乱码| 老司机靠b影院| 国产av又大| 热re99久久精品国产66热6| 丝袜人妻中文字幕| 久久久久国产精品人妻一区二区| 汤姆久久久久久久影院中文字幕| 一级黄色大片毛片| 日韩大码丰满熟妇| 亚洲三区欧美一区| 后天国语完整版免费观看| 久久久精品区二区三区| 亚洲美女黄色视频免费看| 动漫黄色视频在线观看| 后天国语完整版免费观看| 动漫黄色视频在线观看| 美女主播在线视频| 悠悠久久av| 波多野结衣一区麻豆| 老司机靠b影院| 制服诱惑二区| 蜜桃国产av成人99| bbb黄色大片| 天天躁夜夜躁狠狠躁躁| 日本wwww免费看| 嫩草影视91久久| 99久久国产精品久久久| 国产亚洲精品第一综合不卡| 中文字幕色久视频| 男女无遮挡免费网站观看| 婷婷丁香在线五月| 亚洲精品久久成人aⅴ小说| 自拍欧美九色日韩亚洲蝌蚪91| 午夜影院在线不卡| 国产亚洲欧美精品永久| 精品国产一区二区三区四区第35| 2018国产大陆天天弄谢| 在线观看免费视频网站a站| 国产激情久久老熟女| 精品人妻在线不人妻| 国产精品 欧美亚洲| 久久精品熟女亚洲av麻豆精品| √禁漫天堂资源中文www| 久热这里只有精品99| 国产精品欧美亚洲77777| 人妻 亚洲 视频| 成人影院久久| www.999成人在线观看| 精品少妇久久久久久888优播| 黄色视频不卡| 91精品伊人久久大香线蕉| 国产免费现黄频在线看| 天天躁日日躁夜夜躁夜夜| 日韩欧美免费精品| 亚洲va日本ⅴa欧美va伊人久久 | 91精品伊人久久大香线蕉| 亚洲av日韩精品久久久久久密| 日日爽夜夜爽网站| 天天添夜夜摸| 人妻人人澡人人爽人人| 亚洲国产精品一区二区三区在线| 97人妻天天添夜夜摸| 中文字幕另类日韩欧美亚洲嫩草| 欧美xxⅹ黑人| 高清av免费在线| 久久毛片免费看一区二区三区| 亚洲av电影在线观看一区二区三区| 精品乱码久久久久久99久播| 欧美精品av麻豆av| 久久狼人影院| 国产麻豆69| 免费黄频网站在线观看国产| 中国美女看黄片| 老熟妇乱子伦视频在线观看 | 精品熟女少妇八av免费久了| 国产亚洲精品第一综合不卡| 国产精品免费大片| 宅男免费午夜| 久久国产亚洲av麻豆专区| 精品国产一区二区三区久久久樱花| 欧美日韩国产mv在线观看视频| 亚洲七黄色美女视频| 色老头精品视频在线观看| 电影成人av| 热re99久久精品国产66热6| 麻豆国产av国片精品| 久久久久久免费高清国产稀缺| 亚洲一卡2卡3卡4卡5卡精品中文| 精品人妻1区二区| 久久久久视频综合| 日韩精品免费视频一区二区三区| 9色porny在线观看| 国产精品.久久久| 夜夜骑夜夜射夜夜干| 国产成人精品久久二区二区91| 亚洲精品第二区| 国产av又大| 自线自在国产av| av又黄又爽大尺度在线免费看| 日韩欧美一区视频在线观看| 看免费av毛片| 午夜福利在线免费观看网站| 手机成人av网站| 久久久久国内视频| 亚洲 国产 在线| 亚洲av男天堂| av欧美777| 成人18禁高潮啪啪吃奶动态图| 日韩,欧美,国产一区二区三区| 欧美日韩国产mv在线观看视频| 亚洲av日韩在线播放| 欧美成人午夜精品| 操出白浆在线播放| 精品亚洲成a人片在线观看| 婷婷丁香在线五月| kizo精华| 午夜福利视频精品| 亚洲欧美一区二区三区黑人| 亚洲精品一卡2卡三卡4卡5卡 | 老司机亚洲免费影院| 亚洲精华国产精华精| 热re99久久精品国产66热6| 国产国语露脸激情在线看| av在线老鸭窝| 亚洲国产精品一区三区| 视频区欧美日本亚洲| 日韩免费高清中文字幕av| 欧美精品人与动牲交sv欧美| 建设人人有责人人尽责人人享有的| 精品国产乱码久久久久久小说| 在线观看人妻少妇| 精品一区在线观看国产| 国产在线免费精品| 色老头精品视频在线观看| 免费人妻精品一区二区三区视频| 国产高清视频在线播放一区 | 国产又爽黄色视频| 亚洲av电影在线观看一区二区三区| 岛国毛片在线播放| 久久久久网色| 嫁个100分男人电影在线观看| 成人av一区二区三区在线看 | 在线永久观看黄色视频| av国产精品久久久久影院| 亚洲九九香蕉| 亚洲五月色婷婷综合| 日韩有码中文字幕| 久久精品熟女亚洲av麻豆精品| 亚洲情色 制服丝袜| 五月开心婷婷网| 成年人午夜在线观看视频| 国产色视频综合| 十八禁高潮呻吟视频| 欧美中文综合在线视频| 美女高潮到喷水免费观看| 欧美性长视频在线观看| 日韩大码丰满熟妇| 极品少妇高潮喷水抽搐| 一级毛片电影观看| 日韩,欧美,国产一区二区三区| 热99国产精品久久久久久7| 亚洲一区二区三区欧美精品| 成年av动漫网址| 日本猛色少妇xxxxx猛交久久| 亚洲精品美女久久久久99蜜臀| 王馨瑶露胸无遮挡在线观看| 亚洲精品日韩在线中文字幕| 伊人久久大香线蕉亚洲五| 中文精品一卡2卡3卡4更新| 亚洲精品美女久久久久99蜜臀| 国产在视频线精品| 99国产综合亚洲精品| 精品国产国语对白av| 亚洲精华国产精华精| 亚洲精品国产av蜜桃| 国产激情久久老熟女| 亚洲国产日韩一区二区| 男男h啪啪无遮挡| 色94色欧美一区二区| 99国产精品99久久久久| 亚洲精品美女久久av网站| 国产激情久久老熟女| 久久亚洲精品不卡| 2018国产大陆天天弄谢| 免费黄频网站在线观看国产| 亚洲成国产人片在线观看| 国产一区有黄有色的免费视频| 少妇的丰满在线观看| 亚洲视频免费观看视频| 自拍欧美九色日韩亚洲蝌蚪91| 日本欧美视频一区| 大片电影免费在线观看免费| 美女中出高潮动态图| 2018国产大陆天天弄谢| avwww免费| 国产精品一区二区在线不卡| 99国产精品一区二区蜜桃av | 日韩视频一区二区在线观看| 在线十欧美十亚洲十日本专区| 国产一区二区三区av在线| 亚洲国产毛片av蜜桃av| 黄色a级毛片大全视频| 国产成人欧美在线观看 | 国产一级毛片在线| 亚洲成国产人片在线观看| 深夜精品福利| 女性被躁到高潮视频| 久久久国产成人免费| 久久亚洲精品不卡| 99国产极品粉嫩在线观看| 两人在一起打扑克的视频| 久久亚洲国产成人精品v| 中文字幕精品免费在线观看视频| 一区二区三区精品91| 免费少妇av软件| 久久 成人 亚洲| 1024视频免费在线观看| av天堂在线播放| 精品第一国产精品| 91av网站免费观看| 精品人妻在线不人妻| 国产免费一区二区三区四区乱码| 一级,二级,三级黄色视频| 男女午夜视频在线观看| 久久热在线av| 久久99一区二区三区| 狠狠精品人妻久久久久久综合| 免费av中文字幕在线| 欧美激情高清一区二区三区| 亚洲精品在线美女| 熟女少妇亚洲综合色aaa.|