• <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
    国产高清视频在线播放一区| 亚洲熟女毛片儿| 久久久久国产一级毛片高清牌| 视频区欧美日本亚洲| 国产淫语在线视频| 精品国产亚洲在线| 中文字幕精品免费在线观看视频| www.熟女人妻精品国产| 亚洲精品乱久久久久久| 日韩有码中文字幕| 一级片免费观看大全| 亚洲一区中文字幕在线| 在线观看免费午夜福利视频| 俄罗斯特黄特色一大片| a级片在线免费高清观看视频| 在线十欧美十亚洲十日本专区| 亚洲精华国产精华精| 一二三四社区在线视频社区8| 一区二区三区乱码不卡18| 一级毛片精品| 91成人精品电影| 亚洲精品av麻豆狂野| 久久久久国内视频| 亚洲成av片中文字幕在线观看| 久久国产精品男人的天堂亚洲| 国产视频一区二区在线看| 国产日韩一区二区三区精品不卡| 亚洲avbb在线观看| 久久午夜亚洲精品久久| 巨乳人妻的诱惑在线观看| 免费在线观看视频国产中文字幕亚洲| 成年人黄色毛片网站| 亚洲少妇的诱惑av| 老汉色∧v一级毛片| 人人妻人人添人人爽欧美一区卜| 天堂中文最新版在线下载| 1024香蕉在线观看| 精品少妇一区二区三区视频日本电影| 男女之事视频高清在线观看| 国产日韩欧美视频二区| 免费一级毛片在线播放高清视频 | 超碰97精品在线观看| 久久热在线av| 国产黄频视频在线观看| 变态另类成人亚洲欧美熟女 | 18禁国产床啪视频网站| 久久人人97超碰香蕉20202| 一边摸一边抽搐一进一出视频| 美女扒开内裤让男人捅视频| 波多野结衣一区麻豆| 亚洲 欧美一区二区三区| 久久九九热精品免费| 在线天堂中文资源库| 国产成人系列免费观看| 一区二区三区激情视频| 精品久久久精品久久久| 丰满迷人的少妇在线观看| 两人在一起打扑克的视频| 免费看十八禁软件| 老司机影院毛片| kizo精华| 1024香蕉在线观看| 国产精品自产拍在线观看55亚洲 | 亚洲va日本ⅴa欧美va伊人久久| 丰满迷人的少妇在线观看| 欧美黄色淫秽网站| 99国产精品一区二区蜜桃av | 搡老乐熟女国产| 亚洲国产欧美在线一区| 国产精品秋霞免费鲁丝片| 久久久久久久久久久久大奶| 在线av久久热| 国产成人精品久久二区二区91| 久久久久久久精品吃奶| 日韩人妻精品一区2区三区| 久久av网站| 在线亚洲精品国产二区图片欧美| 久久国产亚洲av麻豆专区| 别揉我奶头~嗯~啊~动态视频| 黄片大片在线免费观看| 精品视频人人做人人爽| 久久天躁狠狠躁夜夜2o2o| 无遮挡黄片免费观看| 国产黄频视频在线观看| 国产亚洲av高清不卡| 色婷婷久久久亚洲欧美| 极品人妻少妇av视频| 久久精品成人免费网站| 国产一区二区激情短视频| 国产亚洲精品第一综合不卡| 国产精品亚洲一级av第二区| 九色亚洲精品在线播放| 中文字幕最新亚洲高清| 久久人妻福利社区极品人妻图片| 高清欧美精品videossex| √禁漫天堂资源中文www| 俄罗斯特黄特色一大片| 久久毛片免费看一区二区三区| 91成人精品电影| 亚洲 国产 在线| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品aⅴ一区二区三区四区| 日韩免费av在线播放| 亚洲中文字幕日韩| 色老头精品视频在线观看| 手机成人av网站| 色老头精品视频在线观看| 久久国产精品人妻蜜桃| 两人在一起打扑克的视频| 国产黄频视频在线观看| 精品福利永久在线观看| 久久人妻福利社区极品人妻图片| 欧美日本中文国产一区发布| 人妻 亚洲 视频| 国产在线观看jvid| 狠狠精品人妻久久久久久综合| 午夜两性在线视频| 国产精品免费视频内射| 欧美 亚洲 国产 日韩一| 日日夜夜操网爽| 制服人妻中文乱码| 国产成人精品无人区| 在线永久观看黄色视频| 亚洲成人免费电影在线观看| 国产精品熟女久久久久浪| 咕卡用的链子| 免费一级毛片在线播放高清视频 | 高清毛片免费观看视频网站 | 国产精品美女特级片免费视频播放器 | 国产不卡一卡二| 中文亚洲av片在线观看爽 | 成人18禁高潮啪啪吃奶动态图| 最新的欧美精品一区二区| 日本黄色日本黄色录像| 国产日韩欧美在线精品| 91麻豆精品激情在线观看国产 | 国产在线视频一区二区| 91精品国产国语对白视频| 日韩一区二区三区影片| 亚洲三区欧美一区| 亚洲欧美一区二区三区久久| 国产欧美日韩精品亚洲av| 国产主播在线观看一区二区| 久久久久精品人妻al黑| 岛国毛片在线播放| 日本欧美视频一区| 伊人久久大香线蕉亚洲五| 国产高清视频在线播放一区| 色在线成人网| 99re在线观看精品视频| 国产精品久久久久久精品电影小说| 热99国产精品久久久久久7| 香蕉久久夜色| 国产国语露脸激情在线看| 操美女的视频在线观看| 亚洲人成电影免费在线| 久久性视频一级片| 男人操女人黄网站| 免费在线观看黄色视频的| 黄片小视频在线播放| 婷婷成人精品国产| 色尼玛亚洲综合影院| 亚洲综合色网址| 一边摸一边抽搐一进一小说 | 另类亚洲欧美激情| 欧美在线黄色| 亚洲精品久久成人aⅴ小说| 久热这里只有精品99| 在线观看人妻少妇| 黄色成人免费大全| 国产亚洲av高清不卡| 日本黄色视频三级网站网址 | 啦啦啦在线免费观看视频4| 亚洲精品国产色婷婷电影| 欧美精品av麻豆av| av一本久久久久| av一本久久久久| 黑人欧美特级aaaaaa片| 黑人欧美特级aaaaaa片| 亚洲一区中文字幕在线| 丝瓜视频免费看黄片| 久久亚洲精品不卡| 女人被躁到高潮嗷嗷叫费观| 国产日韩欧美在线精品| 在线观看一区二区三区激情| 国产成人系列免费观看| 成年动漫av网址| 欧美变态另类bdsm刘玥| 亚洲,欧美精品.| 精品少妇黑人巨大在线播放| 99国产精品免费福利视频| 美女高潮到喷水免费观看| 视频区欧美日本亚洲| 国产有黄有色有爽视频| 男女下面插进去视频免费观看| 中亚洲国语对白在线视频| 国产成人影院久久av| 亚洲天堂av无毛| 一区二区av电影网| 少妇的丰满在线观看| www.精华液| 国产亚洲精品一区二区www | 国产欧美日韩精品亚洲av| 欧美av亚洲av综合av国产av| 欧美黑人欧美精品刺激| 国产成人精品久久二区二区91| 中文字幕人妻熟女乱码| 高清av免费在线| 超碰成人久久| 又紧又爽又黄一区二区| 亚洲国产欧美日韩在线播放| 婷婷成人精品国产| 在线观看人妻少妇| 亚洲久久久国产精品| cao死你这个sao货| 汤姆久久久久久久影院中文字幕| 性色av乱码一区二区三区2| 少妇 在线观看| 国产精品 欧美亚洲| 精品国产乱子伦一区二区三区| 国产精品久久久久成人av| 日本一区二区免费在线视频| 成人国语在线视频| 91老司机精品| 久热爱精品视频在线9| 国产91精品成人一区二区三区 | 欧美性长视频在线观看| 成人国产一区最新在线观看| 亚洲性夜色夜夜综合| 亚洲第一欧美日韩一区二区三区 | 日韩一区二区三区影片| 亚洲精品国产区一区二| 久久久久久久国产电影| 久久狼人影院| 91成年电影在线观看| 18在线观看网站| 国产精品久久电影中文字幕 | 亚洲一卡2卡3卡4卡5卡精品中文| 操美女的视频在线观看| 国产1区2区3区精品| 美女国产高潮福利片在线看| 成人特级黄色片久久久久久久 | 国产高清激情床上av| 婷婷成人精品国产| 人人妻,人人澡人人爽秒播| 人人妻人人爽人人添夜夜欢视频| 51午夜福利影视在线观看| 久久久精品国产亚洲av高清涩受| 美女扒开内裤让男人捅视频| 国产在线精品亚洲第一网站| 9色porny在线观看| 久久国产精品人妻蜜桃| 激情在线观看视频在线高清 | 2018国产大陆天天弄谢| aaaaa片日本免费| 久久ye,这里只有精品| 美女福利国产在线| 九色亚洲精品在线播放| 国产深夜福利视频在线观看| 午夜福利乱码中文字幕| 狠狠狠狠99中文字幕| 99久久精品国产亚洲精品| 一级毛片电影观看| 日韩免费高清中文字幕av| 成年人午夜在线观看视频| 亚洲欧美日韩另类电影网站| 国产在线一区二区三区精| 久热这里只有精品99| 五月开心婷婷网| 国产免费福利视频在线观看| 亚洲精品自拍成人| 免费在线观看黄色视频的| 久久人人爽av亚洲精品天堂| 一级a爱视频在线免费观看| 精品亚洲成国产av| 人人妻人人澡人人爽人人夜夜| 亚洲av欧美aⅴ国产| 亚洲成人国产一区在线观看| 欧美国产精品va在线观看不卡| 午夜视频精品福利| 日本精品一区二区三区蜜桃| 久久免费观看电影| 操出白浆在线播放| av福利片在线| 多毛熟女@视频| 午夜久久久在线观看| 久久国产亚洲av麻豆专区| 天堂俺去俺来也www色官网| 人人妻,人人澡人人爽秒播| 国产精品久久久久久精品电影小说| 又黄又粗又硬又大视频| 老司机福利观看| 亚洲成a人片在线一区二区| 国产一区二区在线观看av| 丁香六月欧美| 12—13女人毛片做爰片一| 韩国精品一区二区三区| 国产91精品成人一区二区三区 | 国产精品亚洲av一区麻豆| 国产97色在线日韩免费| 精品一区二区三区av网在线观看 | 国产男靠女视频免费网站| 午夜福利在线免费观看网站| 丝袜人妻中文字幕| 欧美在线一区亚洲| 桃红色精品国产亚洲av| 女人精品久久久久毛片| 色播在线永久视频| 日韩欧美国产一区二区入口| 午夜福利一区二区在线看| 视频区欧美日本亚洲| 黄频高清免费视频| 在线观看免费午夜福利视频| 国产精品一区二区在线不卡| 老司机亚洲免费影院| 日本五十路高清| 99在线人妻在线中文字幕 | 国产精品av久久久久免费| 欧美日韩一级在线毛片| 人妻久久中文字幕网| av不卡在线播放| 免费在线观看影片大全网站| 黄色毛片三级朝国网站| 亚洲成人免费电影在线观看| 国产极品粉嫩免费观看在线| 国产熟女午夜一区二区三区| 亚洲欧洲精品一区二区精品久久久| 欧美性长视频在线观看| 99久久人妻综合| 黄色视频不卡| 黄片大片在线免费观看| 欧美中文综合在线视频| 不卡av一区二区三区| 女人高潮潮喷娇喘18禁视频| 久久久水蜜桃国产精品网| 视频在线观看一区二区三区| 亚洲精品乱久久久久久| 下体分泌物呈黄色| 无遮挡黄片免费观看| 老司机影院毛片| 亚洲黑人精品在线| 欧美黄色淫秽网站| 欧美日韩国产mv在线观看视频| 精品亚洲成国产av| 亚洲第一av免费看| 十八禁网站网址无遮挡| 黄片大片在线免费观看| 中文欧美无线码| 午夜久久久在线观看| 老司机亚洲免费影院| 午夜福利视频在线观看免费| 午夜福利免费观看在线| 国产精品美女特级片免费视频播放器 | 两人在一起打扑克的视频| 免费一级毛片在线播放高清视频 | 精品高清国产在线一区| 成人永久免费在线观看视频 | 18禁黄网站禁片午夜丰满| 天天躁日日躁夜夜躁夜夜| 亚洲精品国产区一区二| 久久精品人人爽人人爽视色| 精品欧美一区二区三区在线| 精品久久久精品久久久| 757午夜福利合集在线观看| 色在线成人网| 欧美变态另类bdsm刘玥| 99riav亚洲国产免费| 精品国产乱子伦一区二区三区| www.自偷自拍.com| 国产一区二区激情短视频| 夜夜爽天天搞| 69精品国产乱码久久久| 免费久久久久久久精品成人欧美视频| 又黄又粗又硬又大视频| 美女午夜性视频免费| 亚洲av欧美aⅴ国产| 91大片在线观看| 午夜福利视频精品| 黄色成人免费大全| 两个人免费观看高清视频| 肉色欧美久久久久久久蜜桃| 色婷婷av一区二区三区视频| av片东京热男人的天堂| 国产1区2区3区精品| av天堂在线播放| 国产精品av久久久久免费| 国产激情久久老熟女| 下体分泌物呈黄色| tube8黄色片| 天天操日日干夜夜撸| 成人18禁高潮啪啪吃奶动态图| 建设人人有责人人尽责人人享有的| 999久久久国产精品视频| 1024香蕉在线观看| 欧美乱妇无乱码| 中文字幕高清在线视频| 欧美精品高潮呻吟av久久| 免费女性裸体啪啪无遮挡网站| 精品国产乱码久久久久久小说| 亚洲国产中文字幕在线视频| 咕卡用的链子| 天天躁狠狠躁夜夜躁狠狠躁| 人人妻人人澡人人爽人人夜夜| 日韩有码中文字幕| 久久热在线av| 自线自在国产av| 欧美日韩亚洲综合一区二区三区_| 两人在一起打扑克的视频| 王馨瑶露胸无遮挡在线观看| tube8黄色片| 欧美精品人与动牲交sv欧美| 91大片在线观看| 老熟妇仑乱视频hdxx| 99国产综合亚洲精品| av又黄又爽大尺度在线免费看| 在线观看66精品国产| 国产精品一区二区在线观看99| 99久久精品国产亚洲精品| 午夜福利影视在线免费观看| 999精品在线视频| av有码第一页| 桃花免费在线播放| 美女国产高潮福利片在线看| 好男人电影高清在线观看| 免费在线观看完整版高清| 性色av乱码一区二区三区2| 男人操女人黄网站| 亚洲成人手机| 亚洲美女黄片视频| 亚洲国产精品一区二区三区在线| 国产亚洲午夜精品一区二区久久| 国产亚洲一区二区精品| 国产一区二区激情短视频| 操美女的视频在线观看| 国产精品国产高清国产av | 日韩中文字幕欧美一区二区| 国产成人啪精品午夜网站| 成人亚洲精品一区在线观看| 国产亚洲欧美精品永久| 免费在线观看影片大全网站| 日本a在线网址| 亚洲欧美日韩高清在线视频 | 亚洲av国产av综合av卡| 国产免费视频播放在线视频| 一进一出好大好爽视频| 亚洲色图 男人天堂 中文字幕| xxxhd国产人妻xxx| 国产视频一区二区在线看| av有码第一页| 午夜成年电影在线免费观看| 天天躁夜夜躁狠狠躁躁| 黄频高清免费视频| 老司机午夜福利在线观看视频 | 欧美在线一区亚洲| 18禁国产床啪视频网站| 一区二区av电影网| 黄色成人免费大全| 日韩视频一区二区在线观看| 欧美国产精品一级二级三级| 欧美亚洲日本最大视频资源| 日本黄色日本黄色录像| 18禁裸乳无遮挡动漫免费视频| 麻豆乱淫一区二区| 国产免费现黄频在线看| www.熟女人妻精品国产| 亚洲精品在线观看二区| tocl精华| 国产精品免费视频内射| 国产亚洲欧美在线一区二区| 久久久国产一区二区| 少妇粗大呻吟视频| 国产成人精品无人区| 啦啦啦在线免费观看视频4| 国产xxxxx性猛交| 成人手机av| 欧美 亚洲 国产 日韩一| 人人妻人人爽人人添夜夜欢视频| 亚洲人成电影免费在线| 下体分泌物呈黄色| 日韩熟女老妇一区二区性免费视频| 婷婷丁香在线五月| tube8黄色片| 人妻久久中文字幕网| 国产高清激情床上av| 成年女人毛片免费观看观看9 | 一区在线观看完整版| 久久精品国产99精品国产亚洲性色 | 国产精品免费一区二区三区在线 | 亚洲久久久国产精品| 国产精品秋霞免费鲁丝片| 亚洲精品一卡2卡三卡4卡5卡| 视频区欧美日本亚洲| svipshipincom国产片| 韩国精品一区二区三区| 1024视频免费在线观看| 久久精品国产综合久久久| 欧美激情 高清一区二区三区| 国产区一区二久久| 另类精品久久| 在线看a的网站| 精品国内亚洲2022精品成人 | 亚洲专区字幕在线| www.999成人在线观看| 欧美日韩亚洲国产一区二区在线观看 | 亚洲av日韩在线播放| 国精品久久久久久国模美| 亚洲精品久久成人aⅴ小说| 国产精品秋霞免费鲁丝片| 亚洲欧美激情在线| 51午夜福利影视在线观看| 国精品久久久久久国模美| 国产av一区二区精品久久| aaaaa片日本免费| 黄网站色视频无遮挡免费观看| 亚洲一码二码三码区别大吗| 热re99久久国产66热| 日韩欧美一区视频在线观看| 日韩大片免费观看网站| 悠悠久久av| 免费人妻精品一区二区三区视频| 成年人午夜在线观看视频| 国产成人精品在线电影| 日日摸夜夜添夜夜添小说| 一本一本久久a久久精品综合妖精| 夫妻午夜视频| 巨乳人妻的诱惑在线观看| 婷婷丁香在线五月| 亚洲第一av免费看| 日本一区二区免费在线视频| 亚洲五月色婷婷综合| 亚洲精品中文字幕一二三四区 | 久久天堂一区二区三区四区| 人妻久久中文字幕网| 国产野战对白在线观看| 国产成人av教育| 999久久久精品免费观看国产| 一边摸一边抽搐一进一小说 | 中文字幕av电影在线播放| 亚洲男人天堂网一区| 免费av中文字幕在线| 国产精品久久久av美女十八| 国产一区二区在线观看av| 在线亚洲精品国产二区图片欧美| 亚洲五月婷婷丁香| 一级a爱视频在线免费观看| 黄频高清免费视频| 久久精品国产综合久久久| 天堂动漫精品| 国产精品一区二区在线不卡| 日韩欧美一区视频在线观看| 超碰97精品在线观看| a在线观看视频网站| 中文字幕精品免费在线观看视频| 黑丝袜美女国产一区| 午夜福利视频精品| 精品人妻熟女毛片av久久网站| 精品欧美一区二区三区在线| 999久久久国产精品视频| 久久九九热精品免费| 成年人午夜在线观看视频| 精品国产一区二区三区久久久樱花| 美女福利国产在线| 日本黄色日本黄色录像| 人人妻人人添人人爽欧美一区卜| av欧美777| 大码成人一级视频| 国产亚洲精品第一综合不卡| 国产欧美日韩综合在线一区二区| 不卡一级毛片| 丝袜人妻中文字幕| 曰老女人黄片| 免费观看人在逋| 夜夜骑夜夜射夜夜干| 黄色视频,在线免费观看| 国产99久久九九免费精品| 国产一区二区三区综合在线观看| 中文字幕最新亚洲高清| 日本精品一区二区三区蜜桃| 日日爽夜夜爽网站| 大陆偷拍与自拍| 99九九在线精品视频| 国产av精品麻豆| 久久热在线av| 免费在线观看完整版高清| 久久人妻av系列| 99精品欧美一区二区三区四区| 国产aⅴ精品一区二区三区波| 18禁国产床啪视频网站| 黄片播放在线免费| 91精品三级在线观看| 我的亚洲天堂| 国产一区二区三区综合在线观看| 国产99久久九九免费精品| 在线天堂中文资源库| 国产精品电影一区二区三区 | 免费少妇av软件| 成人三级做爰电影| 国产xxxxx性猛交| 免费久久久久久久精品成人欧美视频| 日韩精品免费视频一区二区三区| 国产精品免费大片| 久久久精品免费免费高清| av天堂久久9| 午夜福利,免费看| 黄色怎么调成土黄色| 日本黄色日本黄色录像| 91精品国产国语对白视频| 露出奶头的视频| 女人爽到高潮嗷嗷叫在线视频| 香蕉丝袜av| 五月天丁香电影| 欧美精品人与动牲交sv欧美| av免费在线观看网站| 18在线观看网站| 妹子高潮喷水视频|