• <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免费看| 国产又色又爽无遮挡免费看| 婷婷六月久久综合丁香| 亚洲精品久久成人aⅴ小说| 女人被躁到高潮嗷嗷叫费观| 国产99白浆流出| 国产亚洲精品一区二区www| 成人永久免费在线观看视频| 韩国av一区二区三区四区| 女人爽到高潮嗷嗷叫在线视频| 香蕉久久夜色| 免费在线观看视频国产中文字幕亚洲| 岛国在线观看网站| 无限看片的www在线观看| 丝袜美足系列| 91在线观看av| 欧美日本视频| 美女国产高潮福利片在线看| 亚洲三区欧美一区| 不卡一级毛片| 日本五十路高清| 美女午夜性视频免费| 99热只有精品国产| 色播亚洲综合网| 亚洲成国产人片在线观看| 91九色精品人成在线观看| 午夜精品国产一区二区电影| 国产欧美日韩精品亚洲av| 一本久久中文字幕| 夜夜躁狠狠躁天天躁| 欧美日本视频| 又紧又爽又黄一区二区| 国产视频一区二区在线看| 欧美成人午夜精品| 精品久久久精品久久久| 色婷婷久久久亚洲欧美| 精品福利观看| 国产精品乱码一区二三区的特点 | 日本黄色视频三级网站网址| 亚洲电影在线观看av| 亚洲成人免费电影在线观看| 精品不卡国产一区二区三区| 亚洲人成网站在线播放欧美日韩| 9热在线视频观看99| 色播在线永久视频| 夜夜躁狠狠躁天天躁| 一级毛片精品| 国产高清视频在线播放一区| 久久久久久大精品| 别揉我奶头~嗯~啊~动态视频| 国产成人av激情在线播放| 亚洲三区欧美一区| 亚洲中文av在线| 在线观看免费午夜福利视频| 国产成+人综合+亚洲专区| 色在线成人网| 香蕉丝袜av| 极品人妻少妇av视频| 波多野结衣巨乳人妻| 亚洲专区中文字幕在线| 免费高清在线观看日韩| 桃色一区二区三区在线观看| 欧美 亚洲 国产 日韩一| 男女午夜视频在线观看| 日韩欧美免费精品| 99久久综合精品五月天人人| 九色亚洲精品在线播放| 美女高潮喷水抽搐中文字幕| 久久久久国内视频| 怎么达到女性高潮| 欧美日本亚洲视频在线播放| 一个人观看的视频www高清免费观看 | 老汉色av国产亚洲站长工具| 精品无人区乱码1区二区| 亚洲国产欧美网| 97人妻天天添夜夜摸| 伊人久久大香线蕉亚洲五| 国产精品一区二区三区四区久久 | 91av网站免费观看| 欧洲精品卡2卡3卡4卡5卡区| 不卡一级毛片| 变态另类丝袜制服| 久久久国产成人精品二区| 亚洲一区二区三区不卡视频| 在线观看免费视频网站a站| 啪啪无遮挡十八禁网站| 国产精品一区二区在线不卡| 国产日韩一区二区三区精品不卡| 精品少妇一区二区三区视频日本电影| 午夜激情av网站| 亚洲精品在线美女| 成人免费观看视频高清| 99久久精品国产亚洲精品| 国产伦一二天堂av在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 久久天躁狠狠躁夜夜2o2o| 国产av一区在线观看免费| 亚洲专区国产一区二区| 宅男免费午夜| 久9热在线精品视频| 国产精品电影一区二区三区| 亚洲自拍偷在线| 亚洲欧美一区二区三区黑人| 级片在线观看| 国产99白浆流出| 老汉色av国产亚洲站长工具| 国产成人欧美在线观看| 欧美日韩瑟瑟在线播放| 精品人妻在线不人妻| 日本 欧美在线| 级片在线观看| 午夜视频精品福利| 亚洲一码二码三码区别大吗| 国产精品一区二区免费欧美| 亚洲伊人色综图| 极品教师在线免费播放| 成人三级做爰电影| 12—13女人毛片做爰片一| 久久久国产精品麻豆| 女人精品久久久久毛片| 欧美成人午夜精品| 国产精品精品国产色婷婷| 很黄的视频免费| 两个人看的免费小视频| 亚洲精品av麻豆狂野| 一夜夜www| 久久久久久亚洲精品国产蜜桃av| 看片在线看免费视频| 国产精品综合久久久久久久免费 | 久久久久久久精品吃奶| 日本一区二区免费在线视频| 国产片内射在线| 又黄又爽又免费观看的视频| 日韩欧美国产一区二区入口| 日本五十路高清| 亚洲人成77777在线视频| 亚洲色图 男人天堂 中文字幕| 久久久国产成人精品二区| 91在线观看av| av中文乱码字幕在线| 99国产综合亚洲精品| 国产国语露脸激情在线看| 最近最新中文字幕大全电影3 | 日韩欧美在线二视频| 美国免费a级毛片| 在线观看免费视频日本深夜| 人人妻人人澡欧美一区二区 | 黄色丝袜av网址大全| 一本综合久久免费| 欧美一区二区精品小视频在线| 欧美午夜高清在线| 日韩视频一区二区在线观看| 黄片大片在线免费观看| 99久久精品国产亚洲精品| 看免费av毛片| 少妇 在线观看| 国产成人系列免费观看| 久久久精品欧美日韩精品| 十八禁人妻一区二区| 免费在线观看日本一区| 久久亚洲精品不卡| 操出白浆在线播放| 免费在线观看视频国产中文字幕亚洲| 亚洲精品美女久久久久99蜜臀| 很黄的视频免费| 欧美不卡视频在线免费观看 | 亚洲精品久久成人aⅴ小说| 青草久久国产| 久久精品亚洲精品国产色婷小说| 午夜福利,免费看| 99久久国产精品久久久| 一二三四在线观看免费中文在| 亚洲久久久国产精品| 非洲黑人性xxxx精品又粗又长| 国产伦一二天堂av在线观看| 一进一出好大好爽视频| 日本vs欧美在线观看视频| 日韩精品中文字幕看吧| 女性生殖器流出的白浆| 色尼玛亚洲综合影院| 丝袜在线中文字幕| 免费观看人在逋| 日韩欧美国产在线观看| 国产野战对白在线观看| 青草久久国产| 亚洲自偷自拍图片 自拍| 精品人妻在线不人妻| av福利片在线| 亚洲最大成人中文| 88av欧美| 成人国语在线视频| 亚洲精品在线美女| 黄频高清免费视频| 免费在线观看日本一区| 久热这里只有精品99| av福利片在线| 9热在线视频观看99| 99在线人妻在线中文字幕| 女人被狂操c到高潮| 久久国产精品影院| or卡值多少钱| 国产片内射在线| 成人国产一区最新在线观看| 国产成人影院久久av| 久久精品国产综合久久久| 两个人免费观看高清视频| 国产熟女午夜一区二区三区| av免费在线观看网站| 欧美黑人欧美精品刺激| 欧洲精品卡2卡3卡4卡5卡区| 午夜影院日韩av| 色精品久久人妻99蜜桃| 熟妇人妻久久中文字幕3abv| 国产一区二区激情短视频| 校园春色视频在线观看| 国产精品久久久久久亚洲av鲁大| 在线观看免费日韩欧美大片| 欧美日韩瑟瑟在线播放| 免费高清在线观看日韩| 久久青草综合色| 美女大奶头视频| 国产亚洲av嫩草精品影院| 9191精品国产免费久久| 欧美日本中文国产一区发布| 成人特级黄色片久久久久久久| av视频免费观看在线观看| 精品国产一区二区久久| 久久久久久免费高清国产稀缺| 极品教师在线免费播放| 日日爽夜夜爽网站| 午夜福利,免费看| 中文字幕色久视频| 国产精品影院久久| 两性夫妻黄色片| av电影中文网址| 精品乱码久久久久久99久播| 法律面前人人平等表现在哪些方面| 日本免费一区二区三区高清不卡 | 国产午夜福利久久久久久| 性色av乱码一区二区三区2| 高潮久久久久久久久久久不卡| 国产成人系列免费观看| 老司机午夜十八禁免费视频| 欧美日韩黄片免| 男女下面插进去视频免费观看| 久久 成人 亚洲| 午夜福利一区二区在线看| x7x7x7水蜜桃| 国产一区二区三区视频了| 亚洲av熟女| 波多野结衣巨乳人妻| 久久天堂一区二区三区四区| 色综合欧美亚洲国产小说| 纯流量卡能插随身wifi吗| 亚洲精华国产精华精| 午夜日韩欧美国产| 长腿黑丝高跟| 高清在线国产一区| av在线播放免费不卡| 女警被强在线播放| 国产一级毛片七仙女欲春2 | 天堂影院成人在线观看| 精品一品国产午夜福利视频| 中文字幕色久视频| 涩涩av久久男人的天堂| 精品国产一区二区久久| 午夜福利高清视频| 波多野结衣一区麻豆| 免费无遮挡裸体视频| 黄色视频,在线免费观看| 我的亚洲天堂| 婷婷六月久久综合丁香| 欧美另类亚洲清纯唯美| 亚洲成人精品中文字幕电影| 日韩一卡2卡3卡4卡2021年| 黄网站色视频无遮挡免费观看| 国产伦人伦偷精品视频| 欧美日本视频| av有码第一页| 久久人人爽av亚洲精品天堂| 国产精品久久久人人做人人爽| 性欧美人与动物交配| 午夜精品国产一区二区电影| 麻豆国产av国片精品| 超碰成人久久| 美女午夜性视频免费| 亚洲成av人片免费观看| 在线国产一区二区在线| 亚洲国产毛片av蜜桃av| 中国美女看黄片| 成人18禁高潮啪啪吃奶动态图| 曰老女人黄片| 神马国产精品三级电影在线观看 | 亚洲熟女毛片儿| 亚洲久久久国产精品| 美女 人体艺术 gogo| 亚洲午夜理论影院| 丁香欧美五月| 日本 欧美在线| 黄片播放在线免费| 黑人巨大精品欧美一区二区蜜桃| 国产av精品麻豆| 日本在线视频免费播放| 午夜福利高清视频| 欧美日韩乱码在线| 50天的宝宝边吃奶边哭怎么回事| 国产精品,欧美在线| 校园春色视频在线观看| 日韩 欧美 亚洲 中文字幕| 精品一区二区三区四区五区乱码| 亚洲全国av大片| 两性夫妻黄色片| 一级a爱视频在线免费观看| 亚洲av第一区精品v没综合| 一级毛片高清免费大全| 色综合站精品国产| 丁香六月欧美| 免费在线观看完整版高清| 亚洲国产精品久久男人天堂| 久久亚洲精品不卡| √禁漫天堂资源中文www| 美女大奶头视频| 日本vs欧美在线观看视频| 国产精品综合久久久久久久免费 | 亚洲av五月六月丁香网| 国产精品av久久久久免费| 中文字幕高清在线视频| 久久久久久国产a免费观看| 夜夜躁狠狠躁天天躁| 国产一卡二卡三卡精品| 久久影院123| 欧美在线黄色| 亚洲欧美精品综合久久99| 欧美日本亚洲视频在线播放| 女人精品久久久久毛片| 久久国产乱子伦精品免费另类| 久久久精品国产亚洲av高清涩受| 乱人伦中国视频| 69精品国产乱码久久久| 欧美一区二区精品小视频在线| 亚洲成av人片免费观看| 老汉色∧v一级毛片| 又大又爽又粗| 一级黄色大片毛片| 欧美成人一区二区免费高清观看 | 国产欧美日韩一区二区精品| 亚洲电影在线观看av| 日韩大码丰满熟妇| 一边摸一边做爽爽视频免费| 午夜福利视频1000在线观看 | 热99re8久久精品国产| 老熟妇仑乱视频hdxx| 欧美人与性动交α欧美精品济南到| 亚洲三区欧美一区| 欧美成狂野欧美在线观看| 中文亚洲av片在线观看爽| 亚洲精品久久国产高清桃花| 亚洲五月色婷婷综合| 啪啪无遮挡十八禁网站| 少妇 在线观看| 色精品久久人妻99蜜桃| 丝袜美足系列| 国产主播在线观看一区二区| 可以在线观看的亚洲视频| 欧美日韩精品网址| 亚洲av第一区精品v没综合| 久久精品亚洲精品国产色婷小说| 51午夜福利影视在线观看| 欧美乱码精品一区二区三区| 日韩欧美免费精品| 国产激情欧美一区二区| 成人av一区二区三区在线看| 亚洲精品一区av在线观看| 久久人人爽av亚洲精品天堂| 黄色a级毛片大全视频| 99久久99久久久精品蜜桃| 亚洲av电影不卡..在线观看| 韩国av一区二区三区四区| 国产一区二区三区视频了| 男人舔女人下体高潮全视频| 极品教师在线免费播放| 国产成人av教育| 国产精品野战在线观看| 麻豆久久精品国产亚洲av| 日韩欧美一区二区三区在线观看| 欧美日韩亚洲综合一区二区三区_| 黄色片一级片一级黄色片| 中文字幕av电影在线播放| 久久婷婷成人综合色麻豆| 国产私拍福利视频在线观看| 免费看十八禁软件| 日本三级黄在线观看| 国产乱人伦免费视频| 日本 av在线| 非洲黑人性xxxx精品又粗又长| 变态另类成人亚洲欧美熟女 | 日韩欧美在线二视频| 最好的美女福利视频网| 国产熟女午夜一区二区三区| 欧美成人午夜精品| 亚洲黑人精品在线| 看黄色毛片网站| 久久狼人影院| 亚洲一卡2卡3卡4卡5卡精品中文| 老熟妇乱子伦视频在线观看| 淫秽高清视频在线观看| 亚洲av电影不卡..在线观看| 美女大奶头视频| 人妻久久中文字幕网| 精品一品国产午夜福利视频| 久久午夜综合久久蜜桃| 99在线人妻在线中文字幕| 91精品国产国语对白视频| 露出奶头的视频| 精品卡一卡二卡四卡免费| 俄罗斯特黄特色一大片| 欧美激情久久久久久爽电影 | 亚洲精品中文字幕在线视频| 中文字幕人妻熟女乱码| 天天躁夜夜躁狠狠躁躁| 亚洲一区高清亚洲精品| 午夜a级毛片| 午夜成年电影在线免费观看| 亚洲第一青青草原| 中文字幕久久专区| 久久久国产精品麻豆| 91九色精品人成在线观看| 欧美日韩亚洲综合一区二区三区_| 欧美 亚洲 国产 日韩一| 欧美在线一区亚洲| 夜夜看夜夜爽夜夜摸| 成年人黄色毛片网站| 一进一出抽搐gif免费好疼| av福利片在线| 村上凉子中文字幕在线| www.精华液| 91麻豆精品激情在线观看国产| 淫妇啪啪啪对白视频| 国产欧美日韩一区二区三| 少妇的丰满在线观看| 一区二区三区高清视频在线| 18美女黄网站色大片免费观看| 老司机深夜福利视频在线观看| 久久香蕉国产精品| 欧美久久黑人一区二区| avwww免费| 最新在线观看一区二区三区| 亚洲国产欧美日韩在线播放| 一边摸一边抽搐一进一小说| 香蕉久久夜色| АⅤ资源中文在线天堂| 亚洲第一电影网av| 可以免费在线观看a视频的电影网站| 亚洲人成77777在线视频| 成人永久免费在线观看视频| 日韩精品免费视频一区二区三区| 欧美久久黑人一区二区| 国产午夜精品久久久久久| 色播在线永久视频| 日日干狠狠操夜夜爽| 777久久人妻少妇嫩草av网站| 99riav亚洲国产免费| 国产99白浆流出| 日韩有码中文字幕| 97人妻精品一区二区三区麻豆 | 亚洲成国产人片在线观看| 精品一品国产午夜福利视频| 夜夜躁狠狠躁天天躁| 久久精品人人爽人人爽视色| 又黄又爽又免费观看的视频| 91国产中文字幕| 自线自在国产av| 成人国语在线视频| 一进一出好大好爽视频| 免费一级毛片在线播放高清视频 | 亚洲五月色婷婷综合| 国产欧美日韩一区二区精品| 国产亚洲精品第一综合不卡| 国产激情久久老熟女| 88av欧美| 91九色精品人成在线观看| 香蕉丝袜av| 亚洲一区二区三区不卡视频| or卡值多少钱| 视频在线观看一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 最好的美女福利视频网| 别揉我奶头~嗯~啊~动态视频| 国产欧美日韩一区二区三区在线| 精品熟女少妇八av免费久了| 一区二区三区高清视频在线| 中文字幕人妻熟女乱码| 女人被躁到高潮嗷嗷叫费观| 国产成人影院久久av| 国产精品1区2区在线观看.| 日本a在线网址| 熟妇人妻久久中文字幕3abv| 极品教师在线免费播放| 精品高清国产在线一区| 久久香蕉国产精品| 国产精品电影一区二区三区| 一二三四社区在线视频社区8| 91大片在线观看| 久热这里只有精品99| 午夜久久久在线观看| 亚洲国产看品久久| 老司机午夜福利在线观看视频| 国产野战对白在线观看| 两人在一起打扑克的视频| av天堂在线播放| 亚洲欧美激情综合另类| 国产精品久久久久久亚洲av鲁大| 99久久国产精品久久久| 成人国产一区最新在线观看| 国产麻豆成人av免费视频| 天天躁夜夜躁狠狠躁躁| 国产欧美日韩综合在线一区二区| 色哟哟哟哟哟哟| 久久九九热精品免费| 亚洲熟女毛片儿| 午夜福利视频1000在线观看 | 日本免费a在线| 岛国在线观看网站| 最近最新免费中文字幕在线| 黑人巨大精品欧美一区二区mp4| 久久久久国内视频| 激情视频va一区二区三区| 校园春色视频在线观看| 国产真人三级小视频在线观看| 国产成人欧美在线观看| 亚洲成人精品中文字幕电影| 久久久国产成人精品二区| 免费少妇av软件| 国产亚洲欧美精品永久| 精品人妻在线不人妻| 亚洲欧洲精品一区二区精品久久久| 婷婷六月久久综合丁香| 级片在线观看| 免费在线观看完整版高清| 亚洲欧美日韩高清在线视频| 亚洲精品美女久久久久99蜜臀| 人人澡人人妻人| 欧美日本中文国产一区发布| 日日爽夜夜爽网站| 成人免费观看视频高清| 在线观看舔阴道视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲专区字幕在线| 欧美激情高清一区二区三区| 午夜福利影视在线免费观看| 免费观看人在逋| 亚洲午夜精品一区,二区,三区| 亚洲国产欧美一区二区综合| 国产精品一区二区精品视频观看| 母亲3免费完整高清在线观看| 亚洲五月色婷婷综合| 国产精华一区二区三区| 老司机靠b影院| 女警被强在线播放| 最好的美女福利视频网| 久久久久国产一级毛片高清牌| 99re在线观看精品视频| 日本在线视频免费播放| 亚洲成人国产一区在线观看| 高潮久久久久久久久久久不卡| 美女高潮到喷水免费观看| 久久精品国产综合久久久| 一夜夜www| 国语自产精品视频在线第100页| 18禁裸乳无遮挡免费网站照片 | 9191精品国产免费久久| 中文字幕另类日韩欧美亚洲嫩草| 日韩精品中文字幕看吧| 色精品久久人妻99蜜桃| 91字幕亚洲| 中文字幕最新亚洲高清| 午夜亚洲福利在线播放| 999精品在线视频| 精品午夜福利视频在线观看一区| 亚洲精品美女久久久久99蜜臀| 美女国产高潮福利片在线看| 午夜老司机福利片| 亚洲成人国产一区在线观看| 免费在线观看视频国产中文字幕亚洲| 久久香蕉激情| 久久精品国产亚洲av香蕉五月| 男男h啪啪无遮挡| 十八禁人妻一区二区| 国产在线观看jvid| 精品国内亚洲2022精品成人| 日本 欧美在线| 欧美日本亚洲视频在线播放| 国内精品久久久久精免费| 色播亚洲综合网| 午夜免费成人在线视频| 国产成人av激情在线播放| 欧美在线一区亚洲| av天堂在线播放| 50天的宝宝边吃奶边哭怎么回事| 免费人成视频x8x8入口观看| 国产精品影院久久| 91老司机精品| 一个人观看的视频www高清免费观看 | 国产不卡一卡二| 国产成人精品久久二区二区免费| 视频在线观看一区二区三区| 亚洲国产毛片av蜜桃av| 国产精品野战在线观看| 日日干狠狠操夜夜爽| 黄网站色视频无遮挡免费观看| 黄色片一级片一级黄色片| 亚洲自拍偷在线|