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

    齒廓修形人字齒輪副時變嚙合剛度計算方法

    2021-05-17 05:32:08林騰蛟陳夢寒
    振動與沖擊 2021年9期
    關鍵詞:人字修形齒廓

    林騰蛟,陳夢寒,楊 金

    (重慶大學 機械傳動國家重點實驗室,重慶 400044)

    人字齒輪具有高重合度、低軸向力、高承載能力和平穩(wěn)傳動等優(yōu)點,廣泛應用于航空航天、大型船舶等重要傳動裝置中。人字齒輪運轉過程中,輪齒交替嚙合,導致嚙合沖擊及接觸剛度變化,影響傳動系統(tǒng)動態(tài)性能。時變嚙合剛度是齒輪系統(tǒng)振動噪聲的重要激勵因素,齒廓修形在有效改善輪齒嚙合性能的同時,對其時變嚙合剛度也將產生較大影響,因此開展修形人字齒輪副時變嚙合剛度解析算法研究對人字齒輪系統(tǒng)動力學精準分析具有重要的理論及工程意義。

    近年來,國內外學者在輪齒修形及時變嚙合剛度研究方面取得了諸多成果。Litvin等[1-2]改變標準齒條刀具切削刃,采用修形曲線代替齒條直線齒廓,獲得修形齒輪齒面方程;唐進元等[3-5]采用有限元法研究了修形直齒輪的嚙合剛度;Sun等[6]在基于切片法的嚙合剛度算法中引入有限元法與解析算法間誤差的修正系數(shù),建立了修形直齒輪副嚙合剛度準確計算模型;Wan等[7-8]采用積分勢能法計算了斜齒輪的嚙合剛度;魏靜等[9-11]將齒頂修形等效為輪齒誤差,提出了齒廓修形齒輪的嚙合剛度計算方法;王成等[12-13]將齒頂修形量表示為嚙合線長度的函數(shù),建立了考慮齒頂修形的直齒輪嚙合剛度計算模型。在人字齒輪嚙合剛度計算方面,王峰等[14]基于齒面LTCA方法,以嚙合剛度波動幅值為目標,對人字齒輪進行了齒面修形優(yōu)化設計;郭家舜等[15]將齒廓修形作為沿嚙合線方向的時變間隙,研究了修形參數(shù)對人字齒輪副動載系數(shù)幅值的影響規(guī)律。以上文獻取得了重要研究成果,但主要針對修形直、斜齒輪嚙合剛度,解析算法中多將修形量等效為誤差或間隙等因素;目前涉及人字齒輪嚙合剛度的研究報道較少,且均未同時考慮修形輪齒的齒廓幾何特征、重合度變化以及退刀槽對人字齒輪嚙合剛度的影響。

    本文以人字齒輪副為研究對象,提出一種計及齒廓修形參數(shù)和退刀槽寬度的人字齒輪時變嚙合剛度計算方法。建立含齒頂和齒根修形的人字齒輪齒面方程,分析輪齒修形和受載變形對齒輪重合度和接觸狀態(tài)的影響;基于勢能法和數(shù)值積分公式,推導人字齒輪端面嚙合剛度和軸向嚙合剛度,并研究退刀槽寬度、修形參數(shù)、輸入轉矩等對人字齒輪時變嚙合剛度的影響規(guī)律。

    1 齒廓修形人字齒輪齒面方程推導

    1.1 刀具端面齒廓方程

    輪齒齒廓修形包括齒頂和齒根修形,為建立修形齒廓的齒面方程,采用假想刀具加工輪齒,假想刀具齒廓由4段曲線組成,如圖1所示,其中直線1(CD段)、曲線2(DE段)、曲線3(BC段)和橢圓弧4(AB段)分別包絡出被加工齒輪的漸開線區(qū)、齒頂修形區(qū)、齒根修形區(qū)和過渡圓弧區(qū),曲線2、3的方程為y=axn,a為修形曲線系數(shù),由修形參數(shù)確定,n為修形曲線的階次。

    圖1 假想刀具端面齒廓Fig.1 Transverse profile of hypothetical rack cutter

    圖1中S0(x0Oy0)為假想刀具端面齒廓坐標系,Sa(xaOaya)和Sb(xbObyb)分別為曲線2、3的局部坐標系。通過坐標變換,可得端面齒廓上各曲線對應的方程RS0,i(i=1~4),分別表示為

    (1)

    式中:l1為CD上任意點到D點的距離;αt為端面壓力角;“+”號用于刀具右齒面,“-”號用于刀具左齒面。

    (2)

    RS0,3(Δ2,h2,n)=

    (3)

    (4)

    式中:xOe和yOe為橢圓弧圓心坐標;c1和d1分別為橢圓弧長半軸和短半軸;θ為橢圓弧段角度變量。

    圖2為假想刀具與被加工齒輪坐標系位置關系,其中假想刀具1加工出人字齒輪的右旋齒面,刀具2加工出左旋齒面,由坐標變換得到S1(x1O1y1)坐標系中刀具齒面方程RS1,i

    圖2 假想刀具與被加工齒輪坐標系Fig.2 Coordinate systems of hypothetical rack cutter and machined gear

    RS1,i=M10RS0,i

    (5)

    式中:M10為S0到S1的坐標變換矩陣。

    1.2 齒輪齒面方程

    根據(jù)假想刀具與被加工齒輪在齒面上任意點的正確嚙合條件n·v=0,推導出含齒廓修形的人字齒輪在S2(x2O2y2)坐標系中的齒面方程RS2,i

    (6)

    式中:RS1,ix、RS1,iy、RS1,iz分別為S1坐標系中刀具齒面方程x、y、z分量;φi為形成被加工齒輪齒面時S2坐標系的轉動角度;r為被加工齒輪的節(jié)圓半徑。

    2 齒廓修形人字齒輪時變嚙合剛度計算原理

    人字齒輪由左右兩側斜齒輪和退刀槽構成,在工程實踐中人字齒輪的支撐方式通常為主動輪采用軸向可浮動的軸承支撐,從動輪采用能承受少量軸向載荷的軸承支撐,以滿足左右兩側斜齒輪均載要求,而航空航天、大型船舶等重要傳動裝置中人字齒輪制造加工精度較高,嚙合過程中人字齒輪軸向竄動量較小,齒面嚙合狀態(tài)良好,本文在研究人字齒輪時變嚙合剛度時假定齒面嚙合為理想接觸,在嚙合過程中左右兩側斜齒輪同時參與嚙合。輪齒嚙合剛度是隨參與嚙合齒數(shù)和輪齒接觸位置周期變化的時變函數(shù),齒輪接觸線長度具有時變性,由于人字齒輪的接觸線與軸線不平行,單側斜齒輪在嚙合過程中受軸向力作用,因此齒輪嚙合剛度需按端面嚙合剛度和軸向嚙合剛度分別計算。

    當人字齒輪存在退刀槽時,退刀槽的變形會引起輪齒嚙合剛度變化,因此計算人字齒輪副嚙合剛度時還需要計及退刀槽的影響。

    將輪齒變形等效為沿嚙合線方向的彈性變形,在嚙合力F作用下,綜合考慮輪齒發(fā)生的端面變形、軸向變形以及包含退刀槽在內的基體變形而儲存的彈性勢能可表示為[16]

    (7)

    因此嚙合周期內一對輪齒的嚙合剛度,即單齒嚙合剛度表示為

    (8)

    式中:j=1、2,分別表示主、從動輪;Utb、ktb為端面彎曲勢能和剛度,Uts、kts為端面剪切勢能和剛度,Uta、kta為徑向壓縮勢能和剛度,Utf、ktf為基體徑向勢能和剛度,Uab、kab為軸向彎曲勢能和剛度,Uat、kat為軸向扭轉勢能和剛度,Uaf、kaf為基體軸向勢能和剛度,Uh、kh為赫茲接觸勢能和剛度。

    2.1 端面嚙合剛度

    齒輪端面嚙合剛度包括嚙合力F作用下的輪齒彎曲剛度、剪切剛度、徑向壓縮剛度、赫茲接觸剛度和基體徑向剛度五部分。計算人字齒輪端面嚙合剛度時,基于勢能法和切片思想,將輪齒沿齒寬方向切片,見圖3,通過計算各薄片嚙合剛度,求和可得齒輪端面嚙合剛度。修形薄片輪齒的懸臂梁模型如圖4所示。

    圖3 輪齒切片模型Fig.3 Slice model of tooth

    圖4 齒廓修形薄片輪齒懸臂梁模型Fig.4 Cantilever beam model of tooth slice with profile modification

    采用切片法時,厚度為Δz的每一薄片輪齒中端面彎曲勢能、端面剪切勢能和端面壓縮勢能分別表示為

    (9)

    (10)

    (11)

    式中:Mt為嚙合力F作用下輪齒端面所受彎矩,Mt=Fcosβ[cosαz(x(αz)-x)-sinαzy(αz)];x(αz)、y(αz)分別為薄片輪齒嚙合點坐標值,由齒面方程(6)計算可得;nS2,x(αz)、nS2,y(αz)為薄片輪齒嚙合點處齒面法矢的x、y分量;Δz為薄片輪齒齒寬,Δz=(2B+Bt)/N,B為單側斜齒輪齒寬;Bt為退刀槽寬度;N為人字齒輪沿齒寬方向切分的薄片份數(shù);αz為每片薄片輪齒嚙合點處的端面壓力角,αz=arctan(nS2,x(αz)/nS2,y(αz));β為螺旋角;E、G表示齒輪材料的彈性模量和剪切模量。

    對每一薄片齒輪剛度沿人字齒輪的齒寬方向求和,得到人字齒輪的端面彎曲剛度、端面剪切剛度、徑向壓縮剛度分別為

    (12)

    (13)

    (14)

    對于每一薄片輪齒,由于所處軸向位置不同,輪齒嚙合位置也不同,各薄片輪齒上嚙合點的坐標值x(αz)、y(αz)和壓力角αz是關于軸向位置z和嚙合時間t的函數(shù),如圖5所示,其應滿足

    圖5 任意時刻薄片輪齒嚙合狀態(tài)示意圖Fig.5 Sketch map of meshing state of tooth slice at any time

    (15)

    (16)

    (17)

    0≤l(t)≤εαpbt

    (18)

    式中:“+”號用于主動輪,“-”號用于從動輪;v=2πnjrbj,nj為齒輪的轉速,rbj為齒輪基圓半徑;εα為齒輪端面重合度;rsj表示端面嚙合起始點到齒輪圓心的距離,rtj表示任意時刻端面嚙合點到齒輪圓心的距離;βb為基圓螺旋角;pbt為端面基圓齒距。

    當薄片輪齒上的αz、x(αz)、y(αz)不滿足式(15)~(18)時,此薄片輪齒未參與嚙合。

    人字齒輪赫茲接觸勢能和接觸剛度分別為

    (19)

    (20)

    式中:μ為齒輪材料泊松比,L為接觸線長度。

    對每一薄片齒輪基體徑向勢能和基體徑向剛度沿齒寬方向進行求和,即可得到齒輪基體的徑向勢能和徑向剛度,其表達式分別為

    (21)

    (22)

    式中:L*,M*,P*,Q*為尺寸系數(shù)[17],uf,Sf如圖4中所示,可表示為

    (23)

    式中,Z為齒輪的齒數(shù)。

    2.2 軸向嚙合剛度

    人字齒輪軸向嚙合剛度包括輪齒的軸向彎曲剛度、軸向扭轉剛度以及基體的軸向剛度。假設接觸線上所受的嚙合力相等,齒輪在軸向力作用下繞Y軸的彎矩Ma和繞X軸的扭矩MT,如圖6所示,可分別表示為

    圖6 人字齒輪軸向嚙合剛度計算模型Fig.6 Axial mesh stiffness calculation model of double helical gear

    (24)

    式中,Ns表示參與嚙合的輪齒薄片數(shù)量。

    由于單側斜齒輪的軸向力會使輪齒產生軸向彎曲變形和扭轉變形,因此人字齒輪輪齒的軸向嚙合剛度同樣包括軸向彎曲剛度和軸向扭轉剛度。嚙合力F作用下,輪齒軸向彎曲勢能和軸向扭轉勢能表示為

    (25)

    (26)

    (27)

    基于勢能法,輪齒的軸向彎曲剛度和軸向扭轉剛度可表示為

    (28)

    (29)

    在兩側斜齒輪的嚙合力作用下,含退刀槽的人字齒輪基體部分會產生彎曲變形,基體的軸向彎曲勢能可表示為

    (30)

    式中:Maf為作用于齒輪基體的彎矩,Iaf為基體距離齒輪中心x1處的截面慣性矩,其表達式分別為

    (31)

    (32)

    式中:rf為齒根圓半徑;r0為齒輪軸孔半徑;rt為退刀槽外圓半徑。

    因此人字齒輪基體的軸向剛度可表示為

    (33)

    3 考慮受載變形的修形人字齒輪嚙合特性

    齒輪修形參數(shù)與傳遞載荷有關,在實際運行工況下,合理的輪齒修形可有效降低嚙合沖擊,提升齒輪系統(tǒng)運轉平穩(wěn)性,但輪齒修形和受載變形會導致齒輪副重合度的變化,進而影響齒輪時變嚙合剛度,即齒輪副在嚙合過程中同時嚙合的輪齒在單位齒寬上的單齒嚙合剛度之和,因此有必要對輪齒受載變形下的齒廓修形人字齒輪嚙合狀態(tài)和時變嚙合剛度進行分析。

    3.1 修形人字齒輪副的嚙合起止點

    為研究受載變形下齒廓修形人字齒輪副的嚙合狀態(tài),需確定實際載荷作用下輪齒的端面嚙合起止點。人字齒輪同時有多對輪齒參與嚙合,基于修形齒輪的嚙合力F和單齒嚙合剛度k,參與嚙合的各對輪齒沿嚙合線方向的變形δ(αz)可表示為

    (34)

    式中:km(αz)為第m對參與嚙合的單齒嚙合剛度;Nc為參與嚙合的輪齒對數(shù)。

    表1 人字齒輪副基本參數(shù)Tab.1 Basic parameters of double helical gear pair

    為驗證嚙合起止點計算方法的準確性,根據(jù)文中推導的齒廓修形人字齒輪齒面方程,利用ANSYS軟件建立修形人字齒輪副靜力接觸有限元分析模型,對主動輪內圈節(jié)點施加轉矩,并約束其軸向和徑向位移;對從動輪內圈節(jié)點施加全約束,求解得到輪齒的接觸狀態(tài)如圖8所示,提取主動輪嚙合起止點處結點坐標,得出輪齒嚙合起止點到主動輪圓心的距離分別為38.004 mm和41.226 mm,與圖7解析法計算結果對比,兩者誤差分別為0.32%和0.47%,驗證了實際嚙合起止點計算方法的準確性。

    圖7 修形齒輪嚙合起止點計算流程圖Fig.7 Flow chart of calculating meshing start and end points of gear with modification

    3.2 修形人字齒輪副的重合度

    圖9給出了修形人字齒輪副的實際端面重合度計算示意圖,圖中B2、B1為未修形時齒輪端面嚙合起止點,C2、C1為修形后齒輪端面嚙合起止點。

    圖9 齒廓修形齒輪副端面重合度計算示意圖Fig.9 Sketch map of contact ratio calculation of gears with profile modification

    修形人字齒輪的端面重合度可表示為

    (35)

    式中,u為傳動比。

    由于本文只對人字齒輪進行了齒廓修形,對人字齒輪的縱向重合度不產生影響,則縱向重合度為[18]

    (36)

    因此齒廓修形人字齒輪副的總重合度為

    (a) 等效應力云圖

    (37)

    3.3 人字齒輪副時變嚙合剛度

    由于輪齒受載變形后修形人字齒輪接觸狀態(tài)和重合度發(fā)生變化,輪齒的接觸線長度也會隨之變化,考慮輪齒受載變形后修形人字齒輪的時變接觸線長度L可表示為

    L(F,Δ,h,n)=

    (38)

    L(F,Δ,h,n)=

    (39)

    考慮輪齒受載變形時,根據(jù)輪齒的實際嚙合起止點,齒輪副嚙合剛度計算的積分域應滿足

    (40)

    為了驗證含修形人字齒輪嚙合剛度解析算法的正確性,利用文中所給算法和有限元數(shù)值方法分別計算修形人字齒輪副的嚙合剛度,兩者的計算結果如圖10所示,圖中橫坐標T為無量綱時間,表示輪齒從嚙入到嚙出的整個過程。

    圖10 人字齒輪副的嚙合剛度對比曲線Fig.10 Comparison curves of double helical gears mesh stiffness

    對比人字齒輪副時變嚙合剛度解析算法和有限元數(shù)值法的計算結果,兩者吻合良好,嚙合剛度變化規(guī)律一致,兩者計算所得人字齒輪副嚙合剛度的均值分別為17.885 N/(mm·μm)和18.570 N/(mm·μm),誤差為3.83%,在工程允許范圍內。

    4 人字齒輪嚙合剛度影響因素分析

    針對表1所示的人字齒輪副,研究退刀槽寬度、修形量、修形長度、修形曲線階次、輸入轉矩對人字齒輪副重合度和時變嚙合剛度的影響規(guī)律。

    4.1 退刀槽寬度

    在修形參數(shù)和輸入轉矩保持不變情況下,保持單側輪齒寬度不變,選取退刀槽寬度Bt的變化范圍為0~40 mm,計算可得不同退刀槽寬度時人字齒輪副的重合度和嚙合剛度變化規(guī)律如表2和圖11所示。隨著退刀槽寬度的增大,人字齒輪的基體剛度增大,齒輪副嚙合剛度隨之增大;由于改變退刀槽寬度并不影響人字齒輪副的重合度,嚙合剛度波動值變化不大。

    表2 不同退刀槽寬度下的人字齒輪重合度及嚙合剛度Tab.2 Contact ratio and mesh stiffness of double helical gears with different width of run-out groove

    4.2 修形量

    其他參數(shù)保持不變,選取主、從動輪齒頂和齒根修形量變化范圍為0~15 μm,計算得到不同修形量時人字齒輪副的重合度和嚙合剛度變化規(guī)律如表3和圖12所示。齒廓修形后齒頂和齒根修形區(qū)齒厚減薄,輪齒剛度減弱,在輸入轉矩作用下齒輪部分修形區(qū)參與嚙合,齒輪重合度有所減小,致使修形齒輪嚙合剛度小于未修形齒輪。隨著修形量的增大,輪齒齒厚削弱程度增加,輪齒修形區(qū)參與嚙合的區(qū)域減小,齒輪副重合度減小,嚙合剛度均值減小,嚙合剛度波動值隨之變化,當齒輪副總重合度接近整數(shù)時,嚙合剛度波動值相對較大,而當端面重合度接近整數(shù)時,嚙合剛度波動值相對較小。當修形量過大時,齒輪副端面重合度小于1,可能引起嚙合沖擊,降低傳動的平穩(wěn)性。

    圖12 不同修形量時人字齒輪副的嚙合剛度Fig.12 Mesh stiffness of double helical gears with different profile modification amount

    表3 不同修形量時人字齒輪副的重合度及嚙合剛度Tab.3 Contact ratio and mesh stiffness of double helical gears with different profile modification mount

    4.3 修形長度

    其他參數(shù)保持不變,選取主、從動輪的齒頂和齒根修形長度變化范圍為0.6~1.2 mm,計算得到不同修形長度時人字齒輪副的重合度和嚙合剛度變化規(guī)律如表4和圖13所示。隨著修形長度的增大,齒廓修形區(qū)域增大,相同載荷作用下,輪齒嚙合區(qū)域減小,齒輪副重合度降低,導致接觸線長度減少,嚙合剛度均值減小,嚙合剛度波動值隨之變化,當修形長度為1 mm時,齒輪副總重合度最接近于整數(shù)3,嚙合剛度的波動值較大。

    表4 不同修形長度時人字齒輪副的重合度及嚙合剛度Tab.4 Contact ratio and mesh stiffness of double helical gears with different profile modification length

    圖13 不同修形長度時人字齒輪副的嚙合剛度Fig.13 Mesh stiffness of double helical gears with different profile modification length

    4.4 修形曲線階次

    其他參數(shù)保持不變,選取主、從動輪的齒頂和齒根修形曲線階次分別為2、3、4和6,計算可得不同修形曲線階次時人字齒輪副的重合度和嚙合剛度變化規(guī)律如表5和圖14所示。隨著修形曲線的階次增大,修形區(qū)域的齒廓曲率增大,輪齒齒厚削弱程度減小,相同載荷作用下,齒輪副的重合度增大,嚙合剛度均值隨之增大。當修形曲線階次大于4時,由于修形量一般為微米量級,而修形曲線曲率變化很小,齒輪副重合度變化也很小,此時嚙合剛度均值和波動值隨修形曲線階次的變化不大。

    表5 不同修形曲線階次時人字齒輪副的重合度及嚙合剛度Tab.5 Contact ratio and mesh stiffness of double helical gears with different profile modification curve

    圖14 不同修形曲線階次時人字齒輪副的嚙合剛度Fig.14 Mesh stiffness of double helical gears with different profile modification curve

    4.5 輸入轉矩

    其他參數(shù)保持不變,選取輸入轉矩變化范圍為250~875 N·m,計算可得不同輸入轉矩下人字齒輪副的重合度和嚙合剛度變化規(guī)律如表6和圖15所示。隨著輸入轉矩的增加,輪齒受載變形增大,齒廓修形區(qū)參與嚙合的區(qū)域增大,齒輪副的重合度隨之增加,輪齒接觸線長度變長,時變嚙合剛度均值增大,嚙合剛度波動值隨之變化。當輸入轉矩為500 N·m時,齒輪副重合度接近整數(shù)3,嚙合剛度波動值較大;當輸入轉矩較小時,輪齒修形量偏大,齒輪副端面重合度小于1,可能引起嚙合沖擊,降低傳動的平穩(wěn)性。當輸入轉矩大于750 N·m時,嚙合輪齒接觸變形較大,齒廓修形區(qū)幾乎完全進入嚙合,齒輪副重合度不再變化,人字齒輪副嚙合剛度不再受輸入轉矩影響。

    表6 不同輸入轉矩下人字齒輪副的重合度及嚙合剛度Tab.6 Contact ratio and mesh stiffness of double helical gears under different input torque

    圖15 不同輸入轉矩下人字齒輪副的嚙合剛度Fig.15 Mesh stiffness of double helical gears under different input torque

    5 結 論

    (1) 考慮輪齒修形和受載變形對齒輪副重合度和嚙合狀態(tài)的影響,提出了一種計入齒廓修形參數(shù)和退刀槽寬度的人字齒輪嚙合剛度計算方法,與有限元法計算結果對比分析,驗證了解析算法的正確性。

    (2) 齒廓修形后,齒頂和齒根修形區(qū)齒厚減薄,輪齒剛度減弱,受載時人字齒輪部分修形區(qū)參與嚙合,端面重合度減小,修形齒輪嚙合剛度較未修形齒輪有所減小。

    (3) 人字齒輪副退刀槽寬度增大,嚙合剛度均值增大,嚙合剛度波動值變化不大;修形人字齒輪副的重合度及嚙合剛度均值隨修形量、修形長度的增加而減小,隨修形曲線階次的增加而增大,輸入轉矩增大時,嚙合剛度均值先增大而后不再變化;修形量、修形長度和輸入轉矩對齒輪嚙合剛度波動值有較大影響,當總重合度接近整數(shù)時,嚙合剛度波動值較大。

    猜你喜歡
    人字修形齒廓
    一種新型的昆蟲仿生齒廓設計
    風電齒輪螺旋角修形方式的探討
    機械工程師(2021年1期)2021-01-22 02:57:54
    荷葉杯·鴻雁淡描人字
    基于ISO 1328-1:2013的齒廓偏差評定方法研究*
    機電工程(2020年1期)2020-03-04 06:06:48
    齒輪修形在直升機傳動系統(tǒng)中的應用研究
    機械工程師(2018年7期)2018-07-30 09:05:34
    斜齒輪對角修形設計研究
    漸開線齒廓鏈輪梳齒刀的設計
    六圓弧齒廓螺旋齒輪及其嚙合特性
    Summer Clothes
    贊美吧,寫下人字的人
    回族文學(2016年6期)2017-01-11 20:07:52
    免费看a级黄色片| 一个人免费在线观看电影 | 成人三级黄色视频| 一个人看的www免费观看视频| 国产1区2区3区精品| 免费在线观看日本一区| 国产精品 欧美亚洲| 老司机午夜十八禁免费视频| 久久久久国内视频| 免费观看人在逋| 亚洲精品国产精品久久久不卡| 国产高清视频在线观看网站| 亚洲真实伦在线观看| 国内少妇人妻偷人精品xxx网站 | 亚洲狠狠婷婷综合久久图片| 一个人免费在线观看电影 | 国产av麻豆久久久久久久| 高清在线国产一区| 99热这里只有是精品50| 成年免费大片在线观看| 99久久精品国产亚洲精品| tocl精华| 男插女下体视频免费在线播放| 国产精品久久视频播放| 久久久精品欧美日韩精品| 长腿黑丝高跟| 亚洲,欧美精品.| 国产欧美日韩一区二区三| 中文字幕人成人乱码亚洲影| 久久精品人妻少妇| 一区二区三区高清视频在线| 久久国产精品影院| 亚洲va日本ⅴa欧美va伊人久久| 久久午夜亚洲精品久久| 午夜福利在线观看吧| 欧美黑人巨大hd| 国产高清有码在线观看视频| 长腿黑丝高跟| 两性夫妻黄色片| 成人18禁在线播放| 国产亚洲av高清不卡| 精华霜和精华液先用哪个| 一个人免费在线观看的高清视频| 老司机午夜十八禁免费视频| 欧美绝顶高潮抽搐喷水| 搡老岳熟女国产| 露出奶头的视频| 麻豆国产av国片精品| 国产久久久一区二区三区| 老司机深夜福利视频在线观看| 国产亚洲av高清不卡| 久久99热这里只有精品18| 精品免费久久久久久久清纯| 成人亚洲精品av一区二区| 欧美黑人欧美精品刺激| 一本久久中文字幕| 91麻豆av在线| 欧美不卡视频在线免费观看| 男人舔奶头视频| 国产精品九九99| 亚洲成av人片免费观看| 亚洲在线自拍视频| 一进一出抽搐gif免费好疼| 国产精品,欧美在线| 国产精品 国内视频| 亚洲精品456在线播放app | 久久婷婷人人爽人人干人人爱| 午夜亚洲福利在线播放| 精品午夜福利视频在线观看一区| 国产午夜精品论理片| 此物有八面人人有两片| 哪里可以看免费的av片| 午夜激情福利司机影院| 好看av亚洲va欧美ⅴa在| 日本免费一区二区三区高清不卡| 亚洲无线在线观看| 一进一出抽搐动态| 欧美绝顶高潮抽搐喷水| 啦啦啦免费观看视频1| 日韩欧美免费精品| 校园春色视频在线观看| 国产精品国产高清国产av| 精品久久久久久,| 久久久精品欧美日韩精品| 国产麻豆成人av免费视频| 国产精品电影一区二区三区| 97碰自拍视频| 噜噜噜噜噜久久久久久91| bbb黄色大片| 一级毛片女人18水好多| 又粗又爽又猛毛片免费看| 久久久国产精品麻豆| 国产成人aa在线观看| 脱女人内裤的视频| 婷婷精品国产亚洲av在线| 久久久久久久午夜电影| 久99久视频精品免费| 日韩欧美国产在线观看| 国产一区二区三区视频了| 在线观看一区二区三区| 欧美成狂野欧美在线观看| svipshipincom国产片| 亚洲欧美日韩东京热| 日本五十路高清| 成年女人永久免费观看视频| 无人区码免费观看不卡| 亚洲av中文字字幕乱码综合| 99国产极品粉嫩在线观看| 久久久久久久午夜电影| 色综合欧美亚洲国产小说| 国内精品一区二区在线观看| 国产一区二区在线观看日韩 | 最近最新免费中文字幕在线| 男女那种视频在线观看| 国产一区二区三区视频了| 亚洲精品美女久久av网站| 国产美女午夜福利| 特级一级黄色大片| 国产精品一及| 国产成人aa在线观看| netflix在线观看网站| 日本五十路高清| 性色av乱码一区二区三区2| 国产精品亚洲av一区麻豆| 久久精品亚洲精品国产色婷小说| 欧美av亚洲av综合av国产av| 夜夜爽天天搞| 夜夜躁狠狠躁天天躁| 悠悠久久av| 国产午夜精品久久久久久| 女警被强在线播放| 亚洲最大成人中文| 丰满人妻一区二区三区视频av | 免费在线观看成人毛片| 色老头精品视频在线观看| 久久久成人免费电影| 久久久久国内视频| 免费一级毛片在线播放高清视频| 色在线成人网| 欧美成人一区二区免费高清观看 | 国产亚洲欧美在线一区二区| 91麻豆av在线| av在线蜜桃| 国产一区二区激情短视频| 很黄的视频免费| 亚洲欧美日韩东京热| 久久99热这里只有精品18| 美女cb高潮喷水在线观看 | 极品教师在线免费播放| 色综合欧美亚洲国产小说| 在线a可以看的网站| 精品乱码久久久久久99久播| 中文字幕精品亚洲无线码一区| 欧美乱码精品一区二区三区| 久久99热这里只有精品18| 给我免费播放毛片高清在线观看| 日韩欧美在线乱码| 国产欧美日韩一区二区三| 国产淫片久久久久久久久 | 非洲黑人性xxxx精品又粗又长| 久久精品国产综合久久久| 国产99白浆流出| 久久中文看片网| 岛国在线免费视频观看| 久久久国产成人精品二区| 欧美绝顶高潮抽搐喷水| 97人妻精品一区二区三区麻豆| 亚洲性夜色夜夜综合| 国产美女午夜福利| 真人做人爱边吃奶动态| 亚洲人成网站在线播放欧美日韩| 女人高潮潮喷娇喘18禁视频| 久久中文看片网| 久久久久久久午夜电影| 久久亚洲真实| 青草久久国产| 国产成人av激情在线播放| 国产高潮美女av| 琪琪午夜伦伦电影理论片6080| 精品欧美国产一区二区三| 色尼玛亚洲综合影院| 成人精品一区二区免费| 女同久久另类99精品国产91| 成人国产一区最新在线观看| 精品国产三级普通话版| 一二三四社区在线视频社区8| 一夜夜www| 中国美女看黄片| 很黄的视频免费| 日韩大尺度精品在线看网址| 无限看片的www在线观看| 操出白浆在线播放| 欧美日韩福利视频一区二区| 国产毛片a区久久久久| 大型黄色视频在线免费观看| 国产精品一及| 两性夫妻黄色片| 国产激情欧美一区二区| 麻豆久久精品国产亚洲av| 亚洲,欧美精品.| 欧美成狂野欧美在线观看| 免费在线观看日本一区| 人妻夜夜爽99麻豆av| 国产三级黄色录像| 后天国语完整版免费观看| 亚洲精品乱码久久久v下载方式 | 丁香六月欧美| 免费av不卡在线播放| av天堂在线播放| 中文字幕人成人乱码亚洲影| 热99在线观看视频| 色综合站精品国产| 男插女下体视频免费在线播放| 99久国产av精品| 欧美日韩国产亚洲二区| 男女午夜视频在线观看| 久久国产乱子伦精品免费另类| 99久久精品一区二区三区| 婷婷精品国产亚洲av| 最好的美女福利视频网| 国产精品爽爽va在线观看网站| 久9热在线精品视频| 国产精品一区二区免费欧美| 三级男女做爰猛烈吃奶摸视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲av免费在线观看| 天天添夜夜摸| 国内毛片毛片毛片毛片毛片| 久久久国产成人免费| 99精品在免费线老司机午夜| 老熟妇乱子伦视频在线观看| 神马国产精品三级电影在线观看| 中文亚洲av片在线观看爽| 伊人久久大香线蕉亚洲五| 香蕉av资源在线| 狠狠狠狠99中文字幕| 久久久久性生活片| 伦理电影免费视频| 在线a可以看的网站| 99精品欧美一区二区三区四区| 色哟哟哟哟哟哟| 99精品在免费线老司机午夜| 男人舔女人的私密视频| 久久国产乱子伦精品免费另类| 国产精品爽爽va在线观看网站| 国内精品久久久久久久电影| 天天一区二区日本电影三级| 国产野战对白在线观看| 久久天堂一区二区三区四区| 99riav亚洲国产免费| 国产精品一区二区免费欧美| 人妻夜夜爽99麻豆av| 国产精品一区二区三区四区久久| 真人一进一出gif抽搐免费| 国产综合懂色| 小说图片视频综合网站| 一区福利在线观看| 中文字幕熟女人妻在线| 97超级碰碰碰精品色视频在线观看| 国产精品乱码一区二三区的特点| av中文乱码字幕在线| 国产av不卡久久| 午夜免费激情av| 久久精品国产99精品国产亚洲性色| 国产亚洲精品综合一区在线观看| 三级国产精品欧美在线观看 | 这个男人来自地球电影免费观看| 亚洲精品乱码久久久v下载方式 | 国产黄片美女视频| 美女大奶头视频| 热99在线观看视频| 中国美女看黄片| 久久午夜综合久久蜜桃| 网址你懂的国产日韩在线| 婷婷精品国产亚洲av在线| 免费在线观看视频国产中文字幕亚洲| 不卡av一区二区三区| 熟女电影av网| 男人舔奶头视频| 波多野结衣高清作品| 国产亚洲精品久久久com| 伦理电影免费视频| 亚洲精品在线美女| 美女高潮喷水抽搐中文字幕| 麻豆av在线久日| 精品久久久久久成人av| 少妇熟女aⅴ在线视频| 身体一侧抽搐| 日韩三级视频一区二区三区| 亚洲人成电影免费在线| 啦啦啦韩国在线观看视频| 国产乱人视频| 88av欧美| 亚洲色图 男人天堂 中文字幕| xxx96com| 久久精品国产亚洲av香蕉五月| 最好的美女福利视频网| 一本一本综合久久| 亚洲国产高清在线一区二区三| 亚洲无线在线观看| 免费一级毛片在线播放高清视频| 老汉色av国产亚洲站长工具| 亚洲av成人精品一区久久| 巨乳人妻的诱惑在线观看| 国产精品美女特级片免费视频播放器 | 怎么达到女性高潮| 欧美午夜高清在线| 在线播放国产精品三级| 国产精品av视频在线免费观看| 久久久久性生活片| 国产高潮美女av| 高清毛片免费观看视频网站| 女同久久另类99精品国产91| 99国产精品一区二区蜜桃av| www日本在线高清视频| 视频区欧美日本亚洲| 久久人妻av系列| 美女午夜性视频免费| 中文字幕精品亚洲无线码一区| 日本黄大片高清| 亚洲黑人精品在线| 亚洲av五月六月丁香网| 亚洲黑人精品在线| 白带黄色成豆腐渣| 欧美成人性av电影在线观看| 国产精品乱码一区二三区的特点| 欧美性猛交╳xxx乱大交人| 啦啦啦观看免费观看视频高清| 免费av毛片视频| bbb黄色大片| 国产1区2区3区精品| 99视频精品全部免费 在线 | 一进一出抽搐动态| 搞女人的毛片| 欧美+亚洲+日韩+国产| 久久久久久国产a免费观看| 黄色片一级片一级黄色片| 成人欧美大片| av在线蜜桃| 淫秽高清视频在线观看| 免费看日本二区| 欧美日韩精品网址| 深夜精品福利| 欧美日韩福利视频一区二区| 国产亚洲精品一区二区www| 国产三级在线视频| av视频在线观看入口| 一本精品99久久精品77| 精品电影一区二区在线| 国产精品日韩av在线免费观看| 夜夜夜夜夜久久久久| 亚洲午夜精品一区,二区,三区| 黄色视频,在线免费观看| 一级毛片精品| 亚洲精品国产精品久久久不卡| 欧美一级a爱片免费观看看| 国产精品,欧美在线| www.熟女人妻精品国产| 久久这里只有精品中国| 久久久久精品国产欧美久久久| 国产真人三级小视频在线观看| 欧美3d第一页| 女人高潮潮喷娇喘18禁视频| 美女高潮的动态| 亚洲专区国产一区二区| 97碰自拍视频| 757午夜福利合集在线观看| 国内少妇人妻偷人精品xxx网站 | 热99re8久久精品国产| 最近最新免费中文字幕在线| 麻豆av在线久日| 91老司机精品| 亚洲电影在线观看av| 男女做爰动态图高潮gif福利片| 91久久精品国产一区二区成人 | 国产一区二区在线av高清观看| 国产麻豆成人av免费视频| 不卡av一区二区三区| 热99在线观看视频| 可以在线观看的亚洲视频| 啦啦啦韩国在线观看视频| 亚洲男人的天堂狠狠| 国产毛片a区久久久久| 美女免费视频网站| 一本一本综合久久| 波多野结衣巨乳人妻| 国产主播在线观看一区二区| 久久性视频一级片| 免费看光身美女| 成年女人永久免费观看视频| 精品99又大又爽又粗少妇毛片 | 搞女人的毛片| 在线免费观看的www视频| 久久久国产成人免费| 91在线观看av| 国产成人系列免费观看| 亚洲精品中文字幕一二三四区| 日本三级黄在线观看| ponron亚洲| 日本免费a在线| 91老司机精品| 成人三级黄色视频| 精品一区二区三区视频在线 | 国产一区在线观看成人免费| 俄罗斯特黄特色一大片| 欧美日韩亚洲国产一区二区在线观看| 国产高潮美女av| 搡老妇女老女人老熟妇| 美女高潮的动态| 午夜两性在线视频| 久久精品国产亚洲av香蕉五月| 欧美日韩一级在线毛片| 久久久久久久久中文| 欧美在线黄色| 国产精品久久视频播放| 亚洲中文字幕日韩| 成人av在线播放网站| 91字幕亚洲| 99久久精品国产亚洲精品| 亚洲欧美日韩无卡精品| 欧美极品一区二区三区四区| 亚洲激情在线av| 中文在线观看免费www的网站| 欧美午夜高清在线| 看免费av毛片| 亚洲性夜色夜夜综合| 成人永久免费在线观看视频| 色播亚洲综合网| 变态另类成人亚洲欧美熟女| 久久久色成人| 欧美又色又爽又黄视频| 色视频www国产| 日韩成人在线观看一区二区三区| 蜜桃久久精品国产亚洲av| 国产精品av视频在线免费观看| 级片在线观看| 亚洲精品粉嫩美女一区| 一本综合久久免费| 老鸭窝网址在线观看| av福利片在线观看| 青草久久国产| 此物有八面人人有两片| 久久久久国产一级毛片高清牌| 精品国产三级普通话版| 亚洲av成人精品一区久久| 麻豆国产97在线/欧美| 国模一区二区三区四区视频 | 国产单亲对白刺激| 欧美激情在线99| 日本黄大片高清| 国产成人福利小说| 久久久国产成人精品二区| 999久久久精品免费观看国产| 国产精品免费一区二区三区在线| 999久久久精品免费观看国产| 亚洲性夜色夜夜综合| 欧美3d第一页| 国产久久久一区二区三区| www国产在线视频色| 熟妇人妻久久中文字幕3abv| 18禁观看日本| 岛国在线免费视频观看| 成人一区二区视频在线观看| 亚洲国产精品999在线| 一本久久中文字幕| 久久久久久久久中文| 亚洲国产精品999在线| 日本一二三区视频观看| 中文字幕av在线有码专区| 热99在线观看视频| 99久久国产精品久久久| 国产亚洲av高清不卡| 长腿黑丝高跟| 免费看日本二区| 身体一侧抽搐| 免费大片18禁| 欧美一级a爱片免费观看看| 亚洲国产精品sss在线观看| 国产真实乱freesex| 久久精品国产清高在天天线| 久久久国产成人免费| 成人高潮视频无遮挡免费网站| 18禁国产床啪视频网站| 国产成人av教育| 久久久久国产精品人妻aⅴ院| 亚洲av成人av| 女人被狂操c到高潮| 我的老师免费观看完整版| 免费人成视频x8x8入口观看| 国内精品久久久久精免费| 丁香欧美五月| 午夜精品久久久久久毛片777| 岛国在线观看网站| 哪里可以看免费的av片| 色综合欧美亚洲国产小说| 亚洲一区二区三区色噜噜| 69av精品久久久久久| 国产欧美日韩精品一区二区| 国产一区二区在线观看日韩 | 亚洲人成网站在线播放欧美日韩| 岛国视频午夜一区免费看| 精品午夜福利视频在线观看一区| 成人av一区二区三区在线看| 女生性感内裤真人,穿戴方法视频| 亚洲avbb在线观看| 日韩欧美国产一区二区入口| 亚洲精华国产精华精| 国产欧美日韩精品一区二区| 国产一区在线观看成人免费| 男人的好看免费观看在线视频| 国产欧美日韩精品一区二区| 久久婷婷人人爽人人干人人爱| 熟女人妻精品中文字幕| 激情在线观看视频在线高清| 日韩欧美三级三区| 亚洲欧美精品综合久久99| 亚洲成人免费电影在线观看| 日韩欧美免费精品| 看免费av毛片| av欧美777| 久久这里只有精品19| 亚洲色图av天堂| 嫩草影院精品99| 最新美女视频免费是黄的| 中文字幕人成人乱码亚洲影| 国产精品久久视频播放| 亚洲av成人精品一区久久| 美女高潮的动态| 中文亚洲av片在线观看爽| 久久久国产精品麻豆| 中文在线观看免费www的网站| av中文乱码字幕在线| 久久天躁狠狠躁夜夜2o2o| 波多野结衣巨乳人妻| 五月玫瑰六月丁香| 一卡2卡三卡四卡精品乱码亚洲| 国产三级在线视频| 好男人在线观看高清免费视频| 免费在线观看视频国产中文字幕亚洲| 亚洲av美国av| 亚洲精品乱码久久久v下载方式 | 免费在线观看视频国产中文字幕亚洲| 久久久国产精品麻豆| 国产人伦9x9x在线观看| 国产精品影院久久| 成人av在线播放网站| 又紧又爽又黄一区二区| 最近最新免费中文字幕在线| 日韩精品中文字幕看吧| 美女黄网站色视频| 国产成人av激情在线播放| 亚洲激情在线av| 日韩欧美在线二视频| 一级a爱片免费观看的视频| 国产欧美日韩一区二区精品| 日本黄色视频三级网站网址| 欧美激情在线99| 久久久久国产精品人妻aⅴ院| 久久中文字幕人妻熟女| 国产又黄又爽又无遮挡在线| 欧美在线黄色| 国产97色在线日韩免费| 久久久久亚洲av毛片大全| 久久午夜亚洲精品久久| 一级黄色大片毛片| 欧美zozozo另类| 国产伦在线观看视频一区| 九九热线精品视视频播放| 啦啦啦观看免费观看视频高清| 欧美午夜高清在线| 成年女人毛片免费观看观看9| 国产精品久久久av美女十八| 国产成人aa在线观看| 日本黄大片高清| 亚洲成av人片免费观看| 国产黄片美女视频| 国产精品精品国产色婷婷| 成人欧美大片| 国产精品日韩av在线免费观看| 两个人看的免费小视频| 婷婷六月久久综合丁香| 国产视频内射| 国产精品久久久久久亚洲av鲁大| 成人av在线播放网站| 欧美中文日本在线观看视频| 日韩人妻高清精品专区| 99国产精品一区二区蜜桃av| 在线十欧美十亚洲十日本专区| 亚洲精品美女久久久久99蜜臀| 窝窝影院91人妻| 亚洲熟妇熟女久久| 老司机午夜福利在线观看视频| 亚洲真实伦在线观看| 中文字幕精品亚洲无线码一区| 午夜免费成人在线视频| 成年免费大片在线观看| 身体一侧抽搐| 国产麻豆成人av免费视频| 久久久久亚洲av毛片大全| 欧美色视频一区免费| 久久久久久九九精品二区国产| 高清毛片免费观看视频网站| 动漫黄色视频在线观看| 欧美日韩乱码在线| 国产精品99久久99久久久不卡| 亚洲国产精品成人综合色| 欧美日韩精品网址| 91在线精品国自产拍蜜月 | 欧洲精品卡2卡3卡4卡5卡区| 亚洲美女视频黄频| 国产又黄又爽又无遮挡在线| 天堂动漫精品| 国产探花在线观看一区二区| 19禁男女啪啪无遮挡网站| 中文字幕久久专区|