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

    基于微分進化算法的非對稱弧齒錐齒輪約束多目標優(yōu)化

    2014-03-21 05:03:14錢學毅遲建華
    圖學學報 2014年1期
    關(guān)鍵詞:誘導法齒面微分

    錢學毅,遲建華,吳 雙

    (1.武夷學院電子工程系,福建 武夷山 354300;2.大連開發(fā)區(qū)大開產(chǎn)業(yè)發(fā)展研究中心,遼寧 大連 116600;3.浙江經(jīng)貿(mào)職業(yè)技術(shù)學院應(yīng)用工程系,浙江 杭州 310018)

    對弧齒錐齒輪傳動進行約束多目標優(yōu)化設(shè)計是提高產(chǎn)品設(shè)計質(zhì)量、縮短產(chǎn)品設(shè)計周期、提升產(chǎn)品市場競爭力的有效手段。傳統(tǒng)的格利森制弧齒錐齒輪是按照齒面兩側(cè)具有相同嚙合性能的原則進行對稱設(shè)計。但在工程實踐中,大部分傳動裝置如汽車等正轉(zhuǎn)一側(cè)齒面的使用率大大高于另一側(cè)。因此,為充分利用材料,研究弧齒錐齒輪強度的非對稱設(shè)計(在反轉(zhuǎn)齒面用一個小的壓力角、在正轉(zhuǎn)齒面使用一個較大的壓力角)即可以提高齒輪的齒根彎曲強度、齒面接觸強度和膠合強度又可以有效地避免齒頂過分變尖的現(xiàn)象,具有工程實用意義。文獻[1-2]提出了設(shè)計雙壓力角非對稱齒廓齒輪的思想,用以在發(fā)揮大壓力角優(yōu)點的同時又可避開輪齒兩側(cè)壓力角同時增大會導致齒頂強度和柔度下降,抗沖擊能力降低,加載時容易崩齒等缺點。文獻[3]提出了在準雙曲面齒輪輪坯設(shè)計中,摒棄強度的對稱設(shè)計,修正正轉(zhuǎn)用齒面的壓力角,以增加正轉(zhuǎn)用齒面的強度,減少齒根拉應(yīng)力。

    本文在上述設(shè)計思想的基礎(chǔ)上,根據(jù)齒輪嚙合原理和現(xiàn)代摩擦學理論,建立了以弧齒錐齒輪傳動齒面上瞬時接觸線方向與相對滑動速度之間夾角的余切值最小、傳動總體積最小和齒面誘導法曲率主值最小為目標函數(shù)的約束多目標優(yōu)化設(shè)計數(shù)學模型。借助于改進的微分進化多目標優(yōu)化算法,利用MATLAB編制優(yōu)化程序,通過范例對非對稱弧齒錐齒輪傳動進行兩目標和三目標優(yōu)化設(shè)計。并分析了齒面上瞬時接觸線方向與相對滑動速度之間夾角及齒面誘導法曲率主值對齒面接觸強度和膠合強度的影響,為其綜合性能的進一步優(yōu)化提供基礎(chǔ)。

    1 齒面上瞬時接觸線方向與相對滑動速度之間的夾角

    接觸線方向與相對運動速度之間的夾角θ在齒輪嚙合中是一個很重要的角度。它對于兩齒面在相對運動時形成潤滑油膜的條件有很大的影響。根據(jù)彈性流體動力潤滑理論中的Dowson-Higginson線接觸最小油膜厚度公式[4]

    式中:hmin—最小油膜厚度(mm)

    a—潤滑油的黏壓系數(shù)(MPa-1)

    η0—潤滑油的動力黏度(MPa·s)

    U—卷吸速度(mm/s)

    R—當量曲率半徑(mm)

    L—兩齒面接觸線的長度(mm)

    —齒輪材料綜合彈性模量(MPa)

    W—齒面間的法向載荷(N)

    從上式中可見,最小彈流油膜厚度與卷吸速度(即形成油膜的有效速度)有關(guān),卷吸速度越大,最小膜厚越大,潤滑條件越好,抗膠合能力和抗磨損能力越強。卷吸速度等于相對滑動速度在與接觸線垂直方向上的分量。因此,為了有利于形成油膜,改善潤滑條件,接觸線方向與相對運動速度之間的夾角θ應(yīng)盡可能的大。依據(jù)齒輪嚙合原理

    式中:α—法向壓力角(o)

    β—螺旋角(o)

    顯然,由式(2)計算的值cotθ越小越好。

    2 弧齒錐齒輪齒面誘導法曲率主值

    根據(jù)微分幾何理論和齒輪嚙合原理,兩相嚙合齒面在切點沿同一個方向的法曲率之差稱為該方向的誘導法曲率,是相對曲率,它刻劃了兩嚙合齒面在該方向上相對彎曲的程度,也就是切點沿同一方向的鄰近點相對離開的快慢程度。共軛齒面的誘導法曲率是衡量齒輪傳動質(zhì)量的一個重要指標,該值越小,說明兩接觸齒面在這個方向上相對彎曲的程度越弱,當量曲率半徑越大,Hertz彈性接觸應(yīng)力越小,接觸強度越高,同時Hertz應(yīng)力的減小還使齒間最小油膜厚度增加,抗膠合承載能力得以提高。

    過嚙合齒面上的某一切點,沿不同的方向,具有不同的誘導法曲率,必有一個最大值和一個最小值,其最大值和最小值稱為誘導法曲率主值。取得主值的方向稱為主方向,兩個主方向具有正交性。

    弧齒錐齒輪在傳動時是線接觸,兩齒面在接觸線方向τ的誘導法曲率等于零,這是它的一個主方向。誘導法曲率的另一主方向與接觸線方向τ垂直,設(shè)這個方向的誘導法曲率主值為kσ,則kσ是最大值,用式(3)計算[5]:式中:σk—誘導法曲率主值(mm-1)

    δ1,δ2—小輪、大輪的分度錐角(o)L1—小輪分度圓錐的錐距(mm)其他符號意義同前。

    3 微分進化約束多目標優(yōu)化算法

    微分進化算法(Differential Dvolution, DE)是近來提出的簡單而高效的多目標進化算法,是一種采用浮點矢量編碼的全局空間中進行隨機搜索的優(yōu)化算法。目前國內(nèi)外對于多目標微分進化算法的研究尚處于起步階段,相關(guān)研究主要集中在算法性能的改進與完善上。國外具有代表性的有Abbass[6]的自適應(yīng)Pareto微分進化算法,Madavan[7]及Xue等[8]學者的基于Pareto最優(yōu)概念的多目標微分進化算法。國內(nèi)的研究主要有張利彪等[9-11]學者提出的基于極大極小距離密度的多目標微分進化算法。國內(nèi)外學者的研究在一定程度上提高了多目標微分進化算法的性能,但這些研究還僅限于對數(shù)學測試函數(shù)的研究,還不能直接用來進行實際工程的優(yōu)化設(shè)計,尚需要進行改進。

    本文在對現(xiàn)有的微分進化多目標優(yōu)化技術(shù)進一步研究的基礎(chǔ)上,針對其在機械工程設(shè)計中存在的問題,提出了相應(yīng)的策略,改進了現(xiàn)有的微分進化多目標優(yōu)化技術(shù),給出了基于微分進化技術(shù)、適用于機械工程領(lǐng)域的約束多目標優(yōu)化算法。

    3.1 現(xiàn)有微分進化算法在機械設(shè)計中存在的問題

    機械設(shè)計中的多目標優(yōu)化設(shè)計往往是在復雜而苛刻的約束條件下進行。當利用現(xiàn)存的算法進行迭代時,會發(fā)生下面的問題:

    (1)在進行變異和交叉操作產(chǎn)生的新個體有可能在可行域內(nèi),也有可能在不可行域內(nèi),而支配關(guān)系的比較和極大極小距離密度的計算都是在目標函數(shù)空間中進行的,不可行域內(nèi)的個體產(chǎn)生的目標函數(shù)完全有可能支配可行域內(nèi)的個體產(chǎn)生的目標函數(shù),其密度值的計算也同樣具有不確定性,這樣刪除的很可能將是有用的個體,保留的將是無用的個體。最后迭代出的結(jié)果很可能毫無意義。

    (2)若變異和交叉操作只在可行域內(nèi)進行,則每產(chǎn)生一個新個體都要進行可行性檢驗,滿足約束條件則保留,不滿足約束條件則罰掉,再重新搜索,直到滿足約束條件為止。設(shè)計實踐表明,復雜而苛刻的約束條件不但會使收斂速度慢得讓人難以接受,而且約束邊界上的目標個體還很容易使程序陷入死循環(huán)。例如,當進行齒輪設(shè)計時,若將齒數(shù)限定在20~50,當某一目標個體的齒數(shù)恰好是20或50時,用DE的變異和交叉產(chǎn)生新個體的操作就會使程序陷入死循環(huán),導致優(yōu)化設(shè)計無法進行。

    3.2 基于極大極小距離密度實用多目標微分進化算法

    為解決前述多目標微分進化算法所產(chǎn)生的問題,使迭代高速有效,避免程序進入死循環(huán),提出基于極大極小距離密度實用多目標微分進化算法如下:

    (1)設(shè)置初始參數(shù):向量維度D、初始群體規(guī)模NP、縮放因子F、交叉常量CR、最終可行解規(guī)模N;

    (2)在可行域內(nèi)隨機生成初始群體P,規(guī)模為NP;

    (3)將群體P中的非劣解加入Pareto候選解集中;

    (4)在無約束條件下,用DE的變異和交叉操作,對群體P和Pareto候選解集中的每一目標個體x都生成一個新個體x′;

    (5)對每一個新個體x′都進行可行性檢驗,滿足約束條件的保留,否則直接罰掉;

    (6)將所有目標個體和滿足約束條件的新個體一起構(gòu)成群體T;

    (7)如果Pareto候選解集的規(guī)模沒有達到規(guī)定的大小,將所獲得的非劣解直接加入到Pareto候選解集中,直到達到候選解集的設(shè)定規(guī)模;否則,比較個體間的支配關(guān)系,保留支配個體,丟棄被支配個體,直到達到候選解集的規(guī)模;如果所有被支配個體都丟棄以后,候選解集的規(guī)模仍超過設(shè)定值N,則計算Pareto候選解集中每一個體的極大極小距離密度,并刪除密度值最大的個體,直到達到候選解集的設(shè)定規(guī)模;

    (8)如果滿足中止條件,停止迭代,否則返回(4)。

    4 范 例

    非對稱弧齒錐齒輪傳動的承載能力和傳動性能與下面的兩個因素密切相關(guān):

    (1)接觸線方向與相對運動速度之間的夾角θ;

    (2)齒面誘導法曲率主值。接觸線方向與相對運動速度之間的夾角θ越接近于90o,齒面誘導法曲率主值越小,承載能力越高,傳動性能越好。因此非對稱弧齒錐齒輪傳動優(yōu)化設(shè)計是在滿足約束條件下尋求cotθ最小及齒面誘導法曲率主值最小,同時還要尋求齒輪傳動總體積最?。ㄗ罱?jīng)濟),是一個約束多目標優(yōu)化問題,多目標沒有公度性,不宜用傳統(tǒng)的多目標優(yōu)化方法來處理。這是因為傳統(tǒng)的多目標優(yōu)化方法是將多個目標通過加權(quán)求和轉(zhuǎn)化為單目標問題進而進行優(yōu)化。這樣做有以下缺點:

    ① cotθ無量綱,齒面誘導法曲率主值(mm-1),總體積(mm3),各目標單位不一致,無可比性;

    ② 多個目標加權(quán)值的分配帶有較大的主觀色彩,難以反映客觀真貌;

    ③ 優(yōu)化目標僅為多個目標的加權(quán)和,優(yōu)化過程中各目標的優(yōu)度進展不可操作;

    ④ 多個目標之間通過決策變量相互制約,存在相互矛盾,致使加權(quán)目標函數(shù)的拓樸結(jié)構(gòu)十分復雜。改進的微分進化技術(shù)在解決此類問題中得到成功應(yīng)用。

    4.1 已知數(shù)據(jù)

    主動齒輪的轉(zhuǎn)速n1=960r/min,被動齒輪的轉(zhuǎn)速n2=320r/min,傳遞的功率P=100kW;選用齒輪油HL-30潤滑,其運動黏度ν100=30mm2/s;油池溫度θoil=50℃;兩輪齒面粗糙度分別為Ra1=Ra2=0.8μm。

    4.2 選取向量維度D

    選取7個設(shè)計變量(D=7)

    式中,z1為主動輪齒數(shù);mt為端面模數(shù);α、β意義同前;xt為主動輪切向變位系數(shù);x為主動輪高變位系數(shù);b為齒輪寬度(mm)。

    4.3 建立目標函數(shù)

    (1)將齒輪傳動總體積最小作為一個追求的目標,即

    式中,dm1、dm2分別為兩齒輪中點的分度圓直徑(mm)。

    (2)將齒面上瞬時接觸線方向與相對滑動速度之間夾角的余切值最小作為第二個追求的目標,即

    (3)將齒面誘導法曲率主值最小作為同時追求的第3個目標,即

    4.4 確立約束條件

    參照文獻[12-13]及工程實際經(jīng)驗確立膠合強度約束、接觸強度和彎曲強度可靠性約束、齒頂厚約束、齒數(shù)約束、模數(shù)約束、壓力角約束、螺旋角約束、變位系數(shù)約束、齒寬約束等。

    4.5 優(yōu)化計算及結(jié)果遴選

    用MATLAB進行程序設(shè)計計算,設(shè)置初始參數(shù):D=7、NP=50、F=0.3、CR=0.5、N=100,循環(huán)200次,結(jié)果得到100組可行解及相應(yīng)的目標函數(shù)值。以f1(X)和f2(X)為目標函數(shù)的兩目標優(yōu)化非劣最優(yōu)目標域如圖1所示。根據(jù)程序運行結(jié)果,從中選取2個目標函數(shù)都接近最小值,綜合性能最佳的一組結(jié)果,本例中選取圖1中左下邊界點(決策者也可根據(jù)自己的偏好選擇其他的結(jié)果)作為最終方案:X*=[z1,mt,α,β,xt,x,b]T=[12,4, 24.6202, 39.4660, 0.7541, 0.2268, 50]T,總體積V=1.5237×105mm3,齒面上瞬時接觸線方向與相對滑動速度之間夾角θ的余切值f2(X)=cotθ=2.9154,θ=18.93°。以f1(X)、f2(X)和f3(X)為目標函數(shù)的3目標優(yōu)化非劣最優(yōu)目標域如圖2所示。最終方案為:X*=[z1,mt,α,β,xt,x,b]T=[12, 6, 25.8226, 39.8978, 0.8567, 0.5570,17]T,總體積V=2.0980×105mm3,齒面上瞬時接觸線方向與相對滑動速度之間夾角θ的余切值f2(X)=cotθ=2.7459,θ=20.01°,齒面誘導法曲率主值f3(X)=kσ=0.0409mm-1。

    圖1 兩目標優(yōu)化非劣最優(yōu)目標域

    圖2 三目標優(yōu)化非劣最優(yōu)目標域

    4.6 用粒子群算法進行優(yōu)化并與微分進化算法做比較

    用粒子群優(yōu)化算法(particle swarm optimization, PSO)[14-15],參照文獻[16-18]中提出的尋優(yōu)原則和搜索方法對本算例進行優(yōu)化。粒子群規(guī)模為100,進行200次迭代,所得非劣最優(yōu)目標域如圖3所示。比較圖1和圖3可以看出,微分進化算法個體的多樣性和收斂性都好于粒子群算法。優(yōu)化實踐表明,微分進化的收斂速度更是遠遠大于粒子群算法,群體規(guī)模都是100,同樣都是MATLAB編程,同一臺計算機上都是運行200次循環(huán),微分進化僅用了35分鐘,而粒子群算法則用了20.8個小時,這是由于:

    (1)粒子群算法受約束條件的限制,隨著循環(huán)次數(shù)的增加,在可行域內(nèi)搜索新粒子的速度會越來越慢,而DE算法在尋求新個體的過程中不考慮約束條件,不必反復迭代;

    (2)PSO算法在尋求新粒子時需要考慮粒子的歷史痕跡,而DE算法不必顧及個體進化的歷史痕跡,只需考慮當前狀態(tài)即可,相對來講,DE優(yōu)化程序簡單,計算量小,PSO程序復雜,計算量大。

    圖3 PSO非劣最優(yōu)目標域

    5 結(jié) 論

    (1)本文給出了基于極大極小距離密度多目標微分進化算法的改進策略,并將改進的微分進化算法用于弧齒錐齒輪的多目標優(yōu)化設(shè)計,成功地驗證了改進算法的有效性,該方法比傳統(tǒng)的多目標優(yōu)化方法及目前常用的其他現(xiàn)代優(yōu)化方法都具有明顯的科學合理性和實用性。

    (2)本文給出的弧齒錐齒輪傳動系統(tǒng)多目標優(yōu)化設(shè)計方法和得到的設(shè)計參數(shù),能有效地提高產(chǎn)品的設(shè)計質(zhì)量,提升產(chǎn)品的綜合經(jīng)濟技術(shù)指標。該方法進一步拓展了微分進化算法的工程應(yīng)用領(lǐng)域,對于其他復雜機械傳動系統(tǒng)的約束多目標優(yōu)化設(shè)計也具有一定的參考價值。

    [1]Litvin F L, Fuentes A, Howkins M.Designing,generation and TCA of new type of asymmetric face-gear drive with modified geometry [J].Computer Methods in Applied Mechanics and Engineering,2001, 190(43-44): 5837-5865.

    [2]Difrancesco G, Marini S.Structural analysis of asymmetrical teeth: reduction of size and weight [J].Gear Technology, 1997, 14(5): 121-125.

    [3]張金良, 方宗德, 鄧效忠, 吳 斌.準雙曲面齒輪強度的非對稱設(shè)計[J].機械科學與技術(shù), 2005, 24(5):568-602.

    [4]溫詩鑄, 黃 平.摩擦學原理[M].北京: 清華大學出版社, 2002: 73-75.

    [5]吳序堂.齒輪嚙合原理[M].西安: 西安交通大學出版社, 2009: 293-296.

    [6]Abbass H A.The self-adaptive pareto differential evolution algorithm [C]//Proc of the Congress on Evolutionary Compution (CEC’2002).Piscataway:IEEE Service Center, 2002: 831-836.

    [7]Madavan N K.Multiobjective optimization using a pareto differential evolution approach [C]//Proc of the Congress on Evolution Computation (CEC’2002).Piscataway: IEEE Service Center, 2002: 1145-1150.

    [8]Xue F, Sanderson A C, Graves R J.Pareto-based multi-objective differentil evolution [C]//Proc of the 2003 Congress on Evolutionary Computation(CEC’2003).Piscataway, NJ: IEEE press, 2003:862-869.

    [9]張利彪, 許相莉, 馬 銘, 孫彩堂, 周春光.基于微分進化求解多目標優(yōu)化問題中的退化現(xiàn)象[J].吉林大學學報(工學版), 2009, 39(4): 1041-1046.

    [10]張利彪.基于粒子群和微分進化的優(yōu)化算法研究[D].長春: 吉林大學計算機科學與技術(shù)學院, 2007:92-110.

    [11]張利彪, 周春光, 馬 銘, 孫彩堂.基于極大極小距離密度的多目標微分進化算法[J].計算機研究與發(fā)展, 2007, 44(1): 177-184.

    [12]秦大同, 邢子坤, 王建宏.基于動力學和可靠性的風力發(fā)電齒輪傳動系統(tǒng)參數(shù)優(yōu)化設(shè)計[J].機械工程學報, 2008, 44(7): 24-31.

    [13]齒輪手冊編委會.齒輪手冊[M].北京: 機械工業(yè)出版社, 2005: 5-9-5-134.

    [14]Kennedy J, Eberhart R.Particle swarm optimization [C]//IEEE International Conference on Neural Networks:Perth, Australia.Piscataway.NJ: IEEE Service Center,1995: 1942-1948.

    [15]Eberhart R, Kennedy J.A new optimizer using particle swarm theory [C]//Roc.Sixth International Symposium on Micro Machine and Human Science:Nagoya, Japan.Piscataway, NJ: IEEE Service Center,1995: 39-43.

    [16]魏靜萱.解決單目標和多目標優(yōu)化問題的進化算法[D].西安: 西安電子科技大學, 2009: 73-77.

    [17]張 敏.約束優(yōu)化和多目標優(yōu)化的進化算法研究[D].合肥: 中國科學技術(shù)大學, 2008: 9-15.

    [18]張利彪, 周春光, 馬 銘, 劉小華.基于粒子群算法求解多目標優(yōu)化問題[J].計算機研究與發(fā)展,2004, 41(7): 1287-1290.

    猜你喜歡
    誘導法齒面微分
    擬微分算子在Hp(ω)上的有界性
    基于NURBS理論的漸開線齒輪齒面修復
    重型機械(2020年3期)2020-08-24 08:31:46
    純種米曲霉發(fā)酵制備陽江豆豉及其抗氧化活性研究
    上下解反向的脈沖微分包含解的存在性
    誘導法在中學道德與法治課堂教學中的運用
    談思品教育在小學美術(shù)教學中的有效滲透
    基于BP神經(jīng)網(wǎng)絡(luò)的面齒輪齒面粗糙度研究
    借助微分探求連續(xù)函數(shù)的極值點
    高速動車組弧齒錐齒輪齒面疲勞點蝕失效分析
    17CrNiMo6齒輪軸齒面剝落原因分析
    大型鑄鍛件(2015年1期)2016-01-12 06:33:01
    国产免费福利视频在线观看| 色精品久久人妻99蜜桃| 国产亚洲av高清不卡| 亚洲情色 制服丝袜| 操出白浆在线播放| 精品一品国产午夜福利视频| 国产av国产精品国产| 亚洲人成电影观看| 亚洲精品久久成人aⅴ小说| 色老头精品视频在线观看| 欧美一级毛片孕妇| 国产男女超爽视频在线观看| 欧美日韩黄片免| 中文字幕高清在线视频| 别揉我奶头~嗯~啊~动态视频| 久久毛片免费看一区二区三区| 欧美精品啪啪一区二区三区| www.精华液| 欧美亚洲 丝袜 人妻 在线| 国产一卡二卡三卡精品| www日本在线高清视频| 日韩大片免费观看网站| 久久精品人人爽人人爽视色| 亚洲成人免费av在线播放| 一区二区三区国产精品乱码| 99热网站在线观看| 久久国产精品大桥未久av| 嫩草影视91久久| 香蕉久久夜色| av又黄又爽大尺度在线免费看| 免费女性裸体啪啪无遮挡网站| 每晚都被弄得嗷嗷叫到高潮| 亚洲色图综合在线观看| 婷婷丁香在线五月| 母亲3免费完整高清在线观看| 精品人妻熟女毛片av久久网站| avwww免费| 国产又爽黄色视频| 999精品在线视频| 久久精品国产亚洲av高清一级| 亚洲美女黄片视频| 99九九在线精品视频| 亚洲综合色网址| 久久久久国产一级毛片高清牌| 别揉我奶头~嗯~啊~动态视频| 香蕉丝袜av| 少妇精品久久久久久久| 又紧又爽又黄一区二区| 天堂8中文在线网| 久久这里只有精品19| 国内毛片毛片毛片毛片毛片| 国产精品香港三级国产av潘金莲| 日韩中文字幕欧美一区二区| av国产精品久久久久影院| 久久久久久亚洲精品国产蜜桃av| 啪啪无遮挡十八禁网站| 18禁国产床啪视频网站| 午夜视频精品福利| 老司机深夜福利视频在线观看| 国产97色在线日韩免费| 国产成人精品久久二区二区91| tube8黄色片| 九色亚洲精品在线播放| 50天的宝宝边吃奶边哭怎么回事| 免费人妻精品一区二区三区视频| 啦啦啦 在线观看视频| 国产成人欧美| 咕卡用的链子| 首页视频小说图片口味搜索| 黄频高清免费视频| 精品视频人人做人人爽| 老司机在亚洲福利影院| 精品国产乱子伦一区二区三区| 亚洲人成伊人成综合网2020| 1024视频免费在线观看| 午夜福利在线免费观看网站| 国产精品麻豆人妻色哟哟久久| 色老头精品视频在线观看| 丝袜美足系列| 欧美久久黑人一区二区| 欧美一级毛片孕妇| 久久久精品94久久精品| 色尼玛亚洲综合影院| 亚洲色图 男人天堂 中文字幕| 国产男女超爽视频在线观看| 欧美精品高潮呻吟av久久| 超碰成人久久| 正在播放国产对白刺激| 欧美精品一区二区免费开放| 国产三级黄色录像| 在线亚洲精品国产二区图片欧美| 亚洲欧美一区二区三区久久| 欧美精品啪啪一区二区三区| 看免费av毛片| 如日韩欧美国产精品一区二区三区| 午夜福利视频在线观看免费| 五月天丁香电影| 在线av久久热| 一级片免费观看大全| 欧美精品人与动牲交sv欧美| 免费女性裸体啪啪无遮挡网站| 精品一区二区三卡| 亚洲熟女毛片儿| 一级毛片女人18水好多| 欧美黄色淫秽网站| 成年人黄色毛片网站| 十八禁网站免费在线| 精品亚洲成a人片在线观看| 欧美性长视频在线观看| 国产精品欧美亚洲77777| 亚洲伊人久久精品综合| 丝袜人妻中文字幕| 国产有黄有色有爽视频| 欧美激情 高清一区二区三区| 亚洲精品在线观看二区| 亚洲精品一卡2卡三卡4卡5卡| 大片电影免费在线观看免费| 国产黄色免费在线视频| 免费在线观看完整版高清| 精品一品国产午夜福利视频| 国产深夜福利视频在线观看| 侵犯人妻中文字幕一二三四区| 两性午夜刺激爽爽歪歪视频在线观看 | 性少妇av在线| 天堂中文最新版在线下载| 十八禁网站免费在线| 91麻豆精品激情在线观看国产 | 99国产精品一区二区蜜桃av | 涩涩av久久男人的天堂| 午夜精品国产一区二区电影| 黑人巨大精品欧美一区二区蜜桃| 三上悠亚av全集在线观看| 中文亚洲av片在线观看爽 | 女警被强在线播放| 亚洲成人国产一区在线观看| 久久人人97超碰香蕉20202| 老熟妇仑乱视频hdxx| 老司机亚洲免费影院| 99精品在免费线老司机午夜| 色尼玛亚洲综合影院| 国产精品欧美亚洲77777| 9色porny在线观看| 一级毛片电影观看| 精品久久蜜臀av无| 久久久久久久久久久久大奶| kizo精华| 少妇精品久久久久久久| 999精品在线视频| 妹子高潮喷水视频| 中文亚洲av片在线观看爽 | 三级毛片av免费| 夜夜爽天天搞| 国产免费现黄频在线看| 1024香蕉在线观看| 日韩成人在线观看一区二区三区| 亚洲精品中文字幕在线视频| 99国产精品一区二区三区| 人人澡人人妻人| 国产精品一区二区精品视频观看| 美女视频免费永久观看网站| 在线观看免费视频日本深夜| 国产精品 欧美亚洲| 久久久国产一区二区| 色老头精品视频在线观看| 首页视频小说图片口味搜索| 午夜老司机福利片| 免费观看a级毛片全部| 久久天躁狠狠躁夜夜2o2o| 午夜福利在线观看吧| 国产福利在线免费观看视频| 狠狠狠狠99中文字幕| 男女免费视频国产| 精品国产一区二区三区久久久樱花| 亚洲精品粉嫩美女一区| 国产不卡av网站在线观看| 国产真人三级小视频在线观看| 精品亚洲成a人片在线观看| 久久天堂一区二区三区四区| 久久亚洲真实| 久久这里只有精品19| 在线观看免费午夜福利视频| 国产有黄有色有爽视频| av福利片在线| 一区二区三区精品91| 色婷婷久久久亚洲欧美| 十分钟在线观看高清视频www| 精品国产乱码久久久久久男人| 国产亚洲精品第一综合不卡| 国产亚洲av高清不卡| 亚洲成人手机| 人妻久久中文字幕网| 后天国语完整版免费观看| 90打野战视频偷拍视频| 色吧在线观看| av中文乱码字幕在线| 在线观看免费视频日本深夜| 9191精品国产免费久久| 亚洲欧洲精品一区二区精品久久久| 91麻豆av在线| 99久久精品国产亚洲精品| 俄罗斯特黄特色一大片| 18禁国产床啪视频网站| 欧美在线黄色| 日韩人妻高清精品专区| 中文在线观看免费www的网站| 国产成人精品久久二区二区91| 午夜免费观看网址| 美女免费视频网站| 国产三级中文精品| 成人特级黄色片久久久久久久| 天堂av国产一区二区熟女人妻| 亚洲第一电影网av| 精品日产1卡2卡| 最新在线观看一区二区三区| a级毛片在线看网站| av在线天堂中文字幕| 久久中文看片网| 91av网一区二区| av在线蜜桃| 亚洲国产精品久久男人天堂| 久久精品国产99精品国产亚洲性色| 欧美色欧美亚洲另类二区| 在线观看午夜福利视频| 51午夜福利影视在线观看| 亚洲性夜色夜夜综合| 18禁国产床啪视频网站| 久久国产乱子伦精品免费另类| 又黄又爽又免费观看的视频| 欧美色欧美亚洲另类二区| 又大又爽又粗| 欧美日韩乱码在线| 我的老师免费观看完整版| 99热这里只有是精品50| 国产精品亚洲av一区麻豆| 90打野战视频偷拍视频| 国产高潮美女av| 一级毛片女人18水好多| 色尼玛亚洲综合影院| aaaaa片日本免费| 欧美黄色片欧美黄色片| 亚洲国产精品999在线| 2021天堂中文幕一二区在线观| 午夜两性在线视频| 在线观看66精品国产| 757午夜福利合集在线观看| 在线国产一区二区在线| 国产精品电影一区二区三区| 丁香六月欧美| 欧美一级毛片孕妇| 一进一出抽搐动态| 国产亚洲精品久久久com| 国产精品 国内视频| 免费大片18禁| 日韩欧美三级三区| av中文乱码字幕在线| 亚洲成人免费电影在线观看| 亚洲 欧美一区二区三区| 日韩欧美国产在线观看| 亚洲一区高清亚洲精品| 亚洲18禁久久av| 亚洲自拍偷在线| 色在线成人网| 白带黄色成豆腐渣| 久久热在线av| 久久国产精品人妻蜜桃| 天堂影院成人在线观看| 国产激情久久老熟女| 国产成人精品久久二区二区91| 免费看美女性在线毛片视频| 五月伊人婷婷丁香| av天堂在线播放| 丁香欧美五月| 国产三级在线视频| 免费在线观看影片大全网站| 国产精品98久久久久久宅男小说| 日韩国内少妇激情av| 成人一区二区视频在线观看| 国产av麻豆久久久久久久| 天堂网av新在线| 亚洲在线观看片| 欧美极品一区二区三区四区| 热99re8久久精品国产| 久久久国产成人免费| 亚洲狠狠婷婷综合久久图片| 亚洲18禁久久av| 久久天堂一区二区三区四区| 一进一出好大好爽视频| 久久久久亚洲av毛片大全| 三级国产精品欧美在线观看 | 精品久久久久久久人妻蜜臀av| 久久国产精品人妻蜜桃| 色视频www国产| 国产成人aa在线观看| 白带黄色成豆腐渣| 中文字幕av在线有码专区| 高清在线国产一区| 亚洲第一欧美日韩一区二区三区| 日本撒尿小便嘘嘘汇集6| 久久精品国产99精品国产亚洲性色| 少妇的丰满在线观看| 很黄的视频免费| 国产美女午夜福利| aaaaa片日本免费| 麻豆成人午夜福利视频| 亚洲人成伊人成综合网2020| 亚洲欧美激情综合另类| 女人被狂操c到高潮| 国产单亲对白刺激| 蜜桃久久精品国产亚洲av| 99国产精品99久久久久| 窝窝影院91人妻| 天堂网av新在线| 亚洲人成网站高清观看| 99热精品在线国产| 在线看三级毛片| 女人被狂操c到高潮| 成人18禁在线播放| 国产亚洲精品综合一区在线观看| 国产精品亚洲一级av第二区| 日韩中文字幕欧美一区二区| 99视频精品全部免费 在线 | 国产成人精品久久二区二区免费| 国产激情偷乱视频一区二区| 亚洲成a人片在线一区二区| 色尼玛亚洲综合影院| 亚洲欧美日韩高清在线视频| 国产成人av激情在线播放| 色综合婷婷激情| 两性午夜刺激爽爽歪歪视频在线观看| 久久国产精品人妻蜜桃| av片东京热男人的天堂| 免费看十八禁软件| 夜夜爽天天搞| 后天国语完整版免费观看| 日本撒尿小便嘘嘘汇集6| 97超级碰碰碰精品色视频在线观看| 久久久久国产一级毛片高清牌| 高清在线国产一区| 日本三级黄在线观看| 美女扒开内裤让男人捅视频| e午夜精品久久久久久久| 99久久综合精品五月天人人| 国产综合懂色| 欧美一级a爱片免费观看看| 国产成人av激情在线播放| 看免费av毛片| 久久久久免费精品人妻一区二区| 亚洲午夜理论影院| 黑人巨大精品欧美一区二区mp4| 亚洲电影在线观看av| 全区人妻精品视频| 999久久久国产精品视频| 日韩成人在线观看一区二区三区| 婷婷六月久久综合丁香| 亚洲 国产 在线| 精品久久久久久,| 久久精品国产清高在天天线| 久久久久久久午夜电影| 哪里可以看免费的av片| 中文字幕熟女人妻在线| 亚洲精品中文字幕一二三四区| 国产成人精品久久二区二区91| 国产精品 国内视频| 久久99热这里只有精品18| 久久中文看片网| 女生性感内裤真人,穿戴方法视频| 一a级毛片在线观看| 男女午夜视频在线观看| 999久久久国产精品视频| 女人高潮潮喷娇喘18禁视频| 国产99白浆流出| 神马国产精品三级电影在线观看| 一a级毛片在线观看| 国产欧美日韩精品亚洲av| 叶爱在线成人免费视频播放| 哪里可以看免费的av片| 午夜亚洲福利在线播放| 国产精华一区二区三区| 麻豆av在线久日| 欧美日韩福利视频一区二区| x7x7x7水蜜桃| 午夜免费激情av| 毛片女人毛片| 久久久久九九精品影院| 国产成人福利小说| 午夜激情福利司机影院| 国产高潮美女av| 狂野欧美白嫩少妇大欣赏| 成人无遮挡网站| 欧美在线一区亚洲| 男插女下体视频免费在线播放| 国产精品女同一区二区软件 | 97人妻精品一区二区三区麻豆| 国产91精品成人一区二区三区| 香蕉丝袜av| 亚洲五月天丁香| 国产精品美女特级片免费视频播放器 | 黑人欧美特级aaaaaa片| 欧美不卡视频在线免费观看| 一a级毛片在线观看| 99久久无色码亚洲精品果冻| 久9热在线精品视频| 成年女人毛片免费观看观看9| 久久久精品大字幕| 色尼玛亚洲综合影院| 精品久久久久久久人妻蜜臀av| 日韩国内少妇激情av| 1000部很黄的大片| 亚洲成人精品中文字幕电影| 女生性感内裤真人,穿戴方法视频| 最近最新中文字幕大全免费视频| 操出白浆在线播放| 人妻久久中文字幕网| 波多野结衣巨乳人妻| 亚洲人成网站在线播放欧美日韩| 午夜两性在线视频| 美女大奶头视频| 19禁男女啪啪无遮挡网站| 嫁个100分男人电影在线观看| av片东京热男人的天堂| 亚洲av成人精品一区久久| 搞女人的毛片| 不卡一级毛片| 日本 av在线| 久久久精品欧美日韩精品| 欧美性猛交╳xxx乱大交人| 久久午夜亚洲精品久久| 欧美日韩亚洲国产一区二区在线观看| 国内久久婷婷六月综合欲色啪| 国产激情久久老熟女| 国产成+人综合+亚洲专区| 国产乱人视频| 日韩免费av在线播放| 麻豆成人av在线观看| 亚洲欧美日韩无卡精品| 精品乱码久久久久久99久播| 长腿黑丝高跟| 中国美女看黄片| 熟女少妇亚洲综合色aaa.| 看黄色毛片网站| 成人av一区二区三区在线看| 国产av不卡久久| 午夜激情欧美在线| 午夜成年电影在线免费观看| 高潮久久久久久久久久久不卡| 欧美不卡视频在线免费观看| 欧美日韩黄片免| 久久久国产欧美日韩av| 中文字幕人妻丝袜一区二区| 99热6这里只有精品| 性色avwww在线观看| 丁香六月欧美| 人人妻,人人澡人人爽秒播| 88av欧美| 在线观看美女被高潮喷水网站 | 久久香蕉国产精品| 免费人成视频x8x8入口观看| 最新在线观看一区二区三区| 成人欧美大片| 亚洲欧美精品综合一区二区三区| 午夜a级毛片| 俄罗斯特黄特色一大片| 麻豆av在线久日| 欧美日韩亚洲国产一区二区在线观看| 成熟少妇高潮喷水视频| 两性午夜刺激爽爽歪歪视频在线观看| ponron亚洲| 欧美又色又爽又黄视频| 免费看美女性在线毛片视频| 非洲黑人性xxxx精品又粗又长| 国产高清三级在线| 三级男女做爰猛烈吃奶摸视频| 国产精品久久久av美女十八| 99国产精品一区二区蜜桃av| 国产真实乱freesex| 黄色丝袜av网址大全| 成年版毛片免费区| 欧美又色又爽又黄视频| 色综合站精品国产| 狂野欧美白嫩少妇大欣赏| 日本一本二区三区精品| 亚洲色图av天堂| www日本黄色视频网| 51午夜福利影视在线观看| 久久久久久久精品吃奶| 午夜福利在线观看吧| 老司机福利观看| 最新美女视频免费是黄的| 欧美乱色亚洲激情| 成人国产综合亚洲| 99国产精品一区二区三区| 男人舔奶头视频| 九九热线精品视视频播放| 国产成人影院久久av| 法律面前人人平等表现在哪些方面| 91在线观看av| 久久久久久人人人人人| 一二三四社区在线视频社区8| 两性夫妻黄色片| 国产精品久久电影中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 又黄又粗又硬又大视频| 久久人人精品亚洲av| 欧美黑人欧美精品刺激| 两性夫妻黄色片| 久久热在线av| 夜夜看夜夜爽夜夜摸| 无遮挡黄片免费观看| 亚洲性夜色夜夜综合| 国产久久久一区二区三区| 岛国在线观看网站| 最好的美女福利视频网| 最近视频中文字幕2019在线8| 最近在线观看免费完整版| 一区二区三区国产精品乱码| 神马国产精品三级电影在线观看| 久久久久久国产a免费观看| 狂野欧美白嫩少妇大欣赏| 人人妻人人看人人澡| 国产一区二区三区视频了| 日韩三级视频一区二区三区| 国产成人aa在线观看| 午夜日韩欧美国产| 国产高潮美女av| 亚洲av免费在线观看| 麻豆av在线久日| 久久国产精品人妻蜜桃| 999久久久精品免费观看国产| 久久人人精品亚洲av| 亚洲性夜色夜夜综合| 最近在线观看免费完整版| 久久久精品大字幕| 国产精品99久久久久久久久| 又粗又爽又猛毛片免费看| 久久精品国产亚洲av香蕉五月| 一个人免费在线观看电影 | 在线观看免费午夜福利视频| 日日夜夜操网爽| 男人舔女人的私密视频| 偷拍熟女少妇极品色| 天堂影院成人在线观看| 搡老岳熟女国产| 国产麻豆成人av免费视频| 亚洲精品在线美女| 国产一区二区在线av高清观看| 中文在线观看免费www的网站| 99re在线观看精品视频| 女警被强在线播放| 国产精品久久久久久久电影 | 韩国av一区二区三区四区| 欧美日韩综合久久久久久 | 欧美3d第一页| 国产精品亚洲美女久久久| 国产不卡一卡二| 亚洲aⅴ乱码一区二区在线播放| 老司机午夜福利在线观看视频| 不卡av一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 欧美黑人巨大hd| 99久久成人亚洲精品观看| 久久香蕉精品热| 亚洲中文av在线| 国产日本99.免费观看| 天天一区二区日本电影三级| 久久精品国产亚洲av香蕉五月| 午夜免费激情av| 亚洲国产中文字幕在线视频| 女同久久另类99精品国产91| 国产一区二区在线av高清观看| 伦理电影免费视频| 老汉色∧v一级毛片| 中文字幕人妻丝袜一区二区| 99在线视频只有这里精品首页| 日日摸夜夜添夜夜添小说| 麻豆久久精品国产亚洲av| 国产一区二区激情短视频| 黄色片一级片一级黄色片| 免费在线观看视频国产中文字幕亚洲| 国产男靠女视频免费网站| 日韩av在线大香蕉| 国内精品久久久久精免费| 亚洲精品在线观看二区| 九色成人免费人妻av| 久久这里只有精品19| 亚洲自偷自拍图片 自拍| 午夜影院日韩av| 亚洲五月天丁香| 国产精品电影一区二区三区| 麻豆久久精品国产亚洲av| 88av欧美| 久久这里只有精品19| 成年女人毛片免费观看观看9| 免费无遮挡裸体视频| www.999成人在线观看| 国产真实乱freesex| netflix在线观看网站| 国产真人三级小视频在线观看| 欧美在线一区亚洲| 免费在线观看日本一区| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久久大精品| 蜜桃久久精品国产亚洲av| 12—13女人毛片做爰片一| 在线观看美女被高潮喷水网站 | 亚洲精品一卡2卡三卡4卡5卡| 久久婷婷人人爽人人干人人爱| 后天国语完整版免费观看| 亚洲一区二区三区不卡视频| 精品乱码久久久久久99久播| bbb黄色大片| 淫妇啪啪啪对白视频| 国产aⅴ精品一区二区三区波| 全区人妻精品视频| 19禁男女啪啪无遮挡网站| 桃色一区二区三区在线观看|