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

    基于響應(yīng)面法的轎車(chē)主減速齒輪可靠性靈敏度研究

    2016-06-23 03:29:01梁明軒王曉林袁惠群于印鑫
    中國(guó)機(jī)械工程 2016年4期

    梁明軒  王曉林  袁惠群  于印鑫

    1. 中國(guó)計(jì)量學(xué)院, 杭州, 310018 2. 東北大學(xué), 沈陽(yáng), 110819 3.昆明理工大學(xué), 昆明,650500

    基于響應(yīng)面法的轎車(chē)主減速齒輪可靠性靈敏度研究

    梁明軒1王曉林1袁惠群2于印鑫3

    1. 中國(guó)計(jì)量學(xué)院, 杭州, 3100182. 東北大學(xué), 沈陽(yáng), 110819 3.昆明理工大學(xué), 昆明,650500

    摘要:為獲得某轎車(chē)主減速齒輪尺寸參數(shù)對(duì)其傳動(dòng)可靠性的影響,提出采用改進(jìn)拉丁超立方抽樣(LHS)獲取隨機(jī)尺寸樣本的方法;通過(guò)顯式有限元仿真構(gòu)建齒面最大接觸應(yīng)力完全二次響應(yīng)面(RSM),研究了可靠度對(duì)參數(shù)的靈敏度。結(jié)果表明主減速齒輪可靠度對(duì)壓力角和模數(shù)的均值與方差較為敏感,應(yīng)合理控制其取值。

    關(guān)鍵詞:主減速齒輪;改進(jìn)拉丁超立方抽樣;完全二次響應(yīng)面;可靠性靈敏度

    0引言

    汽車(chē)變速器是影響汽車(chē)性能的關(guān)鍵部件之一,其內(nèi)部主減速齒輪副處于常嚙合狀態(tài),承受著整車(chē)的全部驅(qū)動(dòng)負(fù)載,且工況復(fù)雜多變,一旦出現(xiàn)故障不容易替換,該齒輪副的動(dòng)力穩(wěn)健性將會(huì)直接影響到整車(chē)的耐久性和可靠性等性能指標(biāo)[1-3]。影響齒輪副接觸應(yīng)力的因素有很多,采用蒙特卡羅(MonteCarlo)模擬齒輪接觸應(yīng)力分布時(shí),需將不確定因素作為隨機(jī)變量建立概率模型,此時(shí)樣本抽樣方法對(duì)于近似響應(yīng)面模型的建立至關(guān)重要[4]。近20年來(lái),拉丁超立方抽樣方法已廣泛應(yīng)用于預(yù)測(cè)模式的參數(shù)不確定性分析中[5-7]。為了用較少的樣本有效地估計(jì)響應(yīng)變量的均值,McKay等[5]提出了中值拉丁超立方抽樣(medianLatinhypercubesampling,MLHS)方法,它可以在抽取較少樣本的情況下,獲得較高的計(jì)算精度。Jirutitijaroen等[6]提出了含隨機(jī)排序法的拉丁超立方抽樣(LHSwithrandompermutation,LHSRP)方法,該方法需要準(zhǔn)備足夠數(shù)量的素?cái)?shù)并考慮素?cái)?shù)的根逆轉(zhuǎn)換,以上方法均沒(méi)有考慮抽樣參數(shù)的概率分布尾部特征。

    目前,機(jī)械零件可靠性靈敏度分析大多以單個(gè)一體化結(jié)構(gòu)的振動(dòng)響應(yīng)量作為極限狀態(tài)函數(shù)[8-11],對(duì)于含有齒輪耦合接觸非線(xiàn)性的可靠度靈敏性分析在文獻(xiàn)中鮮見(jiàn)。文獻(xiàn)[8]采用不含交叉項(xiàng)的純二次多項(xiàng)式響應(yīng)面法對(duì)翼梁結(jié)構(gòu)的強(qiáng)度可靠度進(jìn)行計(jì)算,獲得了機(jī)翼腹板和上下緣線(xiàn)等基本變量參數(shù)對(duì)結(jié)構(gòu)的可靠性靈敏度。文獻(xiàn)[9]將6σ可靠性方法與協(xié)同優(yōu)化相結(jié)合研究了汽車(chē)車(chē)橋輕量化問(wèn)題,表明協(xié)同優(yōu)化方法得到的最優(yōu)解提高了系統(tǒng)的可靠性。王新剛等[10]重點(diǎn)考慮零件可靠度和失效率受變載荷和強(qiáng)度退化的影響,分析了管形截面扭桿在時(shí)變載荷作用下的動(dòng)態(tài)可靠性靈敏度。

    事實(shí)上,由于設(shè)計(jì)公差、制造誤差和嚙合磨損等因素影響,汽車(chē)主減速齒輪副尺寸分布具有一定的隨機(jī)性,此時(shí)按確定性設(shè)計(jì)的齒輪副仍有可能發(fā)生疲勞失效。為體現(xiàn)參數(shù)分布的重要性,本文首先改進(jìn)了LHS抽樣,利用顯式動(dòng)態(tài)有限元程序仿真獲得了不同參數(shù)樣本下主減速齒輪副嚙合三維接觸應(yīng)力,從而構(gòu)建了齒面接觸極限狀態(tài)完全二次多項(xiàng)式顯式響應(yīng)面,計(jì)算了主減速齒輪副在隨機(jī)尺寸參數(shù)影響下的可靠度,并對(duì)影響嚙合可靠性的基本隨機(jī)變量的靈敏度進(jìn)行了分析,該方法和研究結(jié)果能較好地用于轎車(chē)主減速齒輪副結(jié)構(gòu)優(yōu)化。

    1改進(jìn)的拉丁超立方重要性抽樣

    x1,in=Fi-1((n-0.5)/N)

    (1)

    采用子區(qū)間隨機(jī)抽樣的LHSRP法樣本點(diǎn)x2,in可表示為

    x2,in=Fi-1(ri,n/N+(n-1)/N)

    (2)

    式中,ri,n為[0,1]內(nèi)的獨(dú)立偽隨機(jī)數(shù),n=1,2,…,N。

    當(dāng)采用LHIS法抽樣時(shí),抽取概率密度曲線(xiàn)期望左邊子區(qū)間的右邊界點(diǎn),抽取期望值右邊子區(qū)間的左邊界點(diǎn),這樣就保證抽取更多的期望值附近樣本點(diǎn),更接近實(shí)際樣本分布特征,其樣本點(diǎn)x3,in可以表示為

    (3)

    采用MLHS法、LHSRP法和LHIS法的抽樣過(guò)程如圖1a和圖1b所示,對(duì)于樣本Xi服從期望為4、標(biāo)準(zhǔn)差為2的正態(tài)分布的6次抽樣結(jié)果如圖1c和圖1d所示??梢钥闯鯩LHS法和LHSRP法抽樣更側(cè)重于在子區(qū)間內(nèi)取點(diǎn),體現(xiàn)了“分層”的思想,而LHIS法的關(guān)鍵在于選取子區(qū)間內(nèi)趨近期望值的邊界頂點(diǎn),兼顧了“分層”和“重要性”的思想,從而使抽樣結(jié)果能更好地體現(xiàn)樣本實(shí)際分布,具有較高的收斂性。

    (a)MLHS、LHSRP法概率分布曲線(xiàn)

    (b)LHIS法概率密度曲線(xiàn)

    (c)MLHS、LHSRP法的抽樣結(jié)果(d)本文LHIS法的抽樣結(jié)果圖1 MLHS、LHSRP和LHIS方法的樣點(diǎn)比較圖

    2主減速齒輪三維嚙合應(yīng)力響應(yīng)面

    2.1響應(yīng)面模型理論

    響應(yīng)面法(responsesurfacemethodology,RSM)是通過(guò)一系列確定性的試驗(yàn)擬合一個(gè)響應(yīng)面來(lái)模擬真實(shí)極限狀態(tài)曲面的方法。響應(yīng)面法作為經(jīng)典的代理模型方法包括線(xiàn)性響應(yīng)面、加權(quán)非線(xiàn)性響應(yīng)面、高精度響應(yīng)面以及Kriging響應(yīng)面等[11-12],本文采用顯式完全二次多項(xiàng)式響應(yīng)面擬合隱式極限狀態(tài)函數(shù),便于通過(guò)求可靠度顯式表達(dá)式的偏導(dǎo)數(shù)獲得系統(tǒng)的可靠性靈敏度。假設(shè)n維隨機(jī)輸入變量x和響應(yīng)向量y之間具有如下近似關(guān)系:

    (4)

    式中,a0、ai、aij(i=1,2,…,n; j=i,i+1,…,n)為待定系數(shù)。

    對(duì)輸入變量x的k個(gè)樣本值(x1,x2,…,xk)進(jìn)行仿真試驗(yàn),得到k個(gè)響應(yīng)量(y1,y2,…,yk),則樣本點(diǎn)試驗(yàn)值與近似響應(yīng)面值之間的誤差可以表示為

    (5)

    利用最小二乘法對(duì)式(5)進(jìn)行回歸分析,使誤差項(xiàng)最小,則可以得到

    (6)

    通過(guò)對(duì)式(6)進(jìn)行求解可以得到系統(tǒng)基本隨機(jī)變量和結(jié)構(gòu)響應(yīng)值的近似關(guān)系式中各系數(shù)的估計(jì)值,從而得到響應(yīng)面函數(shù)的顯式表達(dá)式。

    2.2基于顯式有限元法的齒輪嚙合仿真

    某型國(guó)產(chǎn)變速器主減速器輸入齒輪齒數(shù)z1=15,輸出齒輪齒數(shù)z2=65,一擋工況輸入轉(zhuǎn)速w=95.993 rad/s,額定扭矩T=1545.818N·m。在主減速齒輪副的尺寸范圍內(nèi),假設(shè)齒寬、螺旋角、壓力角和模數(shù)服從正態(tài)分布,各變量的取值全部集中在[μ-3σ, μ+3σ],樣本分布均滿(mǎn)足3σ準(zhǔn)則,其均值和標(biāo)準(zhǔn)差如表1所示。

    表1 主減速齒輪尺寸變量均值和標(biāo)準(zhǔn)差

    齒面出現(xiàn)最大接觸應(yīng)力是齒輪出現(xiàn)點(diǎn)蝕失效的直接原因,為了構(gòu)建接觸應(yīng)力響應(yīng)面模型,首先利用LHIS法對(duì)主減速齒輪副尺寸隨機(jī)變量進(jìn)行抽樣,得到仿真樣本點(diǎn),利用ABAQUS/Explicit顯式動(dòng)態(tài)有限元程序仿真得到每組樣本點(diǎn)齒面最大接觸應(yīng)力值。

    圖2為某型主減速齒輪副有限元模型,分別在兩齒輪軸心建立參考點(diǎn)RP-1與RP-2,建立如圖2所示坐標(biāo)系,將齒輪內(nèi)圈與參考點(diǎn)完全耦合,通過(guò)參考點(diǎn)施加邊界條件和載荷,約束參考點(diǎn)x、y、z方向位移以及繞x、y軸的旋轉(zhuǎn),保留繞z軸旋轉(zhuǎn)自由度;在RP-1點(diǎn)上施加轉(zhuǎn)速ω,在RP-2上施加負(fù)載扭矩T,仿真時(shí)間為20 ms,取齒面最大應(yīng)力作為響應(yīng)真實(shí)值。

    圖3為樣本初始取值時(shí)齒面最大接觸應(yīng)力云圖,可以看出在動(dòng)態(tài)嚙合過(guò)程中齒面最大接觸應(yīng)力為947.62 MPa。

    圖 2 ABAQUS/Explicit中主減速齒輪有限元模型

    圖 3 初始樣本值下齒輪接觸應(yīng)力云圖

    利用式(4)~式(6)進(jìn)行數(shù)據(jù)回歸分析可以得到主減速齒輪副接觸應(yīng)力響應(yīng)面模型,其系數(shù)如表2所示。

    表2 響應(yīng)面函數(shù)系數(shù)項(xiàng)

    3主減速器齒輪嚙合可靠性及靈敏度分析

    3.1可靠性分析及可靠度

    假設(shè)Pmax為主減速器齒輪副齒面許用接觸應(yīng)力強(qiáng)度因子,G(X)是接觸狀態(tài)響應(yīng)面函數(shù),齒輪接觸存在兩種狀態(tài):G(X)≥Pmax失效狀態(tài);G(X)

    μG=E(G(X))=μG(μ1,μ2,…,μk,D1,D2,…,Dk)

    (7)

    DG=D(G(X))=DG(μ1,μ2,…,μk,D1,D2,…,Dk)

    (8)

    由此可以得到狀態(tài)函數(shù)的可靠性指標(biāo)表達(dá)式為

    (9)

    對(duì)于任意分布的狀態(tài)函數(shù)可以采用蒙特卡羅法模擬得到其分布特征參數(shù),若G(X)服從正態(tài)分布,則系統(tǒng)可靠度可表示為R=ψ(β),其中ψ(·)為標(biāo)準(zhǔn)正態(tài)分布[10]。取設(shè)計(jì)要求的主減速齒輪嚙合極限接觸強(qiáng)度因子Pmax=985 MPa。圖4為利用蒙特卡羅法通過(guò)2000次模擬得到的主減速齒面最大接觸應(yīng)力響應(yīng)狀態(tài)函數(shù)G(X)概率直方圖和正態(tài)概率圖。從圖4可以看出G(X)服從正態(tài)分布,得到狀態(tài)函數(shù)均值和標(biāo)準(zhǔn)差為:μG=949.4799 MPa,σG=351.6672 MPa。通過(guò)式(9)計(jì)算得到β=2.70,從而得到其可靠度R=0.9965。

    (a)狀態(tài)函數(shù)G(X)概率直方圖

    (b)狀態(tài)函數(shù)G(X)正態(tài)概率圖圖 4 主減速齒輪副最大接觸應(yīng)力函數(shù)概率分布圖

    3.2主減速器接觸可靠性靈敏度分析

    參數(shù)的靈敏度分析是確定各參數(shù)對(duì)模式預(yù)測(cè)結(jié)果的貢獻(xiàn)大小,均值靈敏度反映了樣本均值大小對(duì)可靠性的影響程度,方差靈敏度反映了樣本參數(shù)波動(dòng)性對(duì)可靠性的影響程度,它有利于在進(jìn)一步的研究過(guò)程中有效減小模式預(yù)測(cè)參數(shù)的不確定性程度。主減速齒輪接觸強(qiáng)度可靠度對(duì)基本隨機(jī)變量的均值和方差的靈敏度為

    (10)

    (11)

    利用式(10)和式(11)將μG和DG對(duì)矩陣μ和D中各元素求偏導(dǎo)數(shù),略去求解過(guò)程,得到

    4結(jié)論

    (1)利用改進(jìn)的拉丁超立方抽樣對(duì)某型主減速齒輪副基本結(jié)構(gòu)尺寸進(jìn)行了抽樣,采用完全二次多項(xiàng)式顯式響應(yīng)面法對(duì)其接觸可靠性進(jìn)行了分析。

    (2)在主減速齒輪副基本尺寸樣本服從正態(tài)分布的情況下,由蒙特卡羅法模擬得到的最大齒面接觸應(yīng)力服從正態(tài)分布。主減速齒輪副模數(shù)分布是影響其可靠性的最主要因素;壓力角和模數(shù)均值的增大會(huì)降低其接觸可靠度,齒寬和螺旋角的增大會(huì)增大其可靠度;基本尺寸隨機(jī)波動(dòng)對(duì)其接觸可靠性影響相對(duì)較小。該方法可為轎車(chē)主減速齒輪嚙合傳動(dòng)優(yōu)化設(shè)計(jì)提供參考。

    參考文獻(xiàn):

    [1]Rao S S, Tjandra M .Reliability-based Design of Automotive Transmission Systems[J]. Reliability Engineering & System Safety, 1994, 46(2): 159-169.

    [2]李昌, 韓興. 基于響應(yīng)面法齒輪嚙合傳動(dòng)可靠性靈敏度分析[J]. 航空動(dòng)力學(xué)報(bào), 2011, 26(3): 711-715.

    Li Chang, Han Xing. Analysis of Reliability Sensitivity for Gear Engagement Based on Response Surface Methods[J]. Journal of Aerospace Power, 2011, 26(3): 711-715.

    [3]尚蛟, 馮慧華, 魯守衛(wèi), 等. 轎車(chē)變速器主減速器齒輪副動(dòng)態(tài)性能仿真[J]. 振動(dòng)工程學(xué)報(bào), 2010, 23(6): 665-669.

    Shang Jiao, Feng Huihua, Lu Shouwei. et al. Simulation Research on Dynamic Performance for Gear Pair of Main Gear Box in Transmission of a Car[J]. Journal of Vibration Engineering, 2010, 23(6): 665-669.

    [4]Hosseini S A A, Khadem S E. Vibration and Reliability of a Rotating Beam with Random Properties under Random Excitation[J]. International Journal of Mechanical Sciences, 2007, 49(12): 1377-1388.

    [5]Mckay M D, Beckman R J, Conover W J. A Comparison of Three Methods for Selecting Values of Input Variables in the Analysis of Output from a Computer Code[J]. Technometrics, 1979, 21(2): 239-245.

    [6]Jirutitijaroen P, Singh C. Comparison of Simulation Methods for Power System Reliability Indexes and Their Distributions[J]. IEEE Transactions on Power Systems, 2008, 23(2): 486-493.

    [7]Echard B, Gayton N, Lemaire M. AK-MCS: an Active Learning Reliability Method Combing Kriging and Monte Carlo Simiulation[J]. Structural Safety, 2011, 33(2): 145-154.[8]高宗戰(zhàn), 劉志群, 姜志峰, 等. 飛機(jī)翼梁結(jié)構(gòu)強(qiáng)度可靠性靈敏度分析[J]. 機(jī)械工程學(xué)報(bào), 2010, 46(14): 194-198.

    Gao Zongzhan, Liu Zhiqun, Jiang Zhifeng, et al. Strength Reliability and Parameter Sensitivity Analysis for Airfoil Spar Structure[J]. Journal of Mechanical Engineering, 2010, 46(14): 194-198.

    [9]姜武華, 王其東, 何煒, 等. 基于可靠性的車(chē)橋結(jié)構(gòu)多學(xué)科優(yōu)化[J]. 汽車(chē)工程, 2012, 34(12): 1136-1140.

    Jiang Wuhua, Wang Qidong, He Wei, et al. Reliability-based Multidisciplinary Optimization of Vehicle Axle Structure[J]. Automotive Engineering, 2012, 34(12): 1136-1140.

    [10]王新剛, 張義民, 王寶艷. 機(jī)械零部件的動(dòng)態(tài)可靠性靈敏度分析[J]. 機(jī)械工程學(xué)報(bào), 2010,46(10): 188-193.

    Wang Xingang, Zhang Yimin, Wang Baoyan. Dynamic Reliability Sensitivity Analysis of Mechanical Components[J]. Journal of Mechanical Engineering, 2010,46(10): 188-193.

    [11]呂震宙, 宋述芳, 李紅雙, 等. 結(jié)構(gòu)機(jī)構(gòu)可靠性及可靠性靈敏度分析[M]. 北京: 科學(xué)出版社, 2009.

    [12]呂震宙, 趙潔, 岳珠峰. 機(jī)械可靠性分析的高精度響應(yīng)面法[J]. 應(yīng)用數(shù)學(xué)和力學(xué), 2007, 28(1):16-24.

    Lü Zhenzhou, Zhao Jie, Yue Zhufeng. Advanced Response Surface Method for Mechanical Reliability Analysis[J]. Applied Mathematics and Mechanics, 2007, 28(1):16-24.

    (編輯王艷麗)

    ReliabilitySensitivityforMainReductionGearsofaCarBasedonResponseSurfaceMethods

    LiangMingxuan1WangXiaolin1YuanHuiqun2YuYinxin2

    1.ChinaJiliangUniversity,Hangzhou, 3100182.NortheasternUniversity,Shenyang,1108193.KunmingUniversityofScienceandTechnology,Kunming,650500

    Abstract:In order to obtain the effects of size parameters on reliability of the main reduction gears in a car, an improved LHS was presented to acquire the random size samples. The full quadratic RSM of the maximum contact stress on tooth surface was built through explicit finite element simulations to study reliability sensitivity to the key parameters. The results show that the mean and variance of pressure angle and modulus which are more sensitive to reliability of the gears must be selected reasonably in design stage.

    Key words:main reduction gear; improved Latin hypercube sampling; full quadratic response surface model; reliability sensitivity

    收稿日期:2015-05-04

    基金項(xiàng)目:國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863計(jì)劃)資助項(xiàng)目(2012AA040104);浙江省“儀器科學(xué)與技術(shù)”重中之重學(xué)科開(kāi)放基金資助項(xiàng)目

    中圖分類(lèi)號(hào):TK122

    DOI:10.3969/j.issn.1004-132X.2016.04.023

    作者簡(jiǎn)介:梁明軒,男,1986年生。中國(guó)計(jì)量學(xué)院機(jī)電工程學(xué)院講師。研究方向?yàn)辇X輪箱傳動(dòng)系統(tǒng)動(dòng)力學(xué)及其可靠性、轉(zhuǎn)子系統(tǒng)非線(xiàn)性動(dòng)力學(xué)及振動(dòng)控制。發(fā)表論文6篇。王曉林,男,1957年生。中國(guó)計(jì)量學(xué)院機(jī)電工程學(xué)院教授。袁惠群,男,1954年生。東北大學(xué)理學(xué)院教授、博士研究生導(dǎo)師。于印鑫,男,1983年生。昆明理工大學(xué)機(jī)電工程學(xué)院講師。

    爱豆传媒免费全集在线观看| 久久中文字幕一级| 女人被躁到高潮嗷嗷叫费观| 王馨瑶露胸无遮挡在线观看| 两人在一起打扑克的视频| 色精品久久人妻99蜜桃| 啦啦啦 在线观看视频| 三上悠亚av全集在线观看| 成人18禁高潮啪啪吃奶动态图| 国产熟女午夜一区二区三区| 丝袜美腿诱惑在线| 大型av网站在线播放| 成人手机av| 亚洲精品日韩在线中文字幕| 国产精品成人在线| 国产精品秋霞免费鲁丝片| 一级片'在线观看视频| 黑人猛操日本美女一级片| 中国国产av一级| √禁漫天堂资源中文www| 国产高清国产精品国产三级| 各种免费的搞黄视频| av天堂久久9| 精品国产国语对白av| av国产久精品久网站免费入址| 精品国产国语对白av| 国产一区有黄有色的免费视频| 一区二区三区四区激情视频| 婷婷丁香在线五月| 后天国语完整版免费观看| 波野结衣二区三区在线| 日韩精品免费视频一区二区三区| 亚洲国产欧美日韩在线播放| 手机成人av网站| 高清不卡的av网站| av网站在线播放免费| 香蕉国产在线看| 狠狠婷婷综合久久久久久88av| a级片在线免费高清观看视频| 国产高清视频在线播放一区 | 99久久综合免费| 91老司机精品| 激情五月婷婷亚洲| 精品一区二区三区四区五区乱码 | 日本vs欧美在线观看视频| 精品熟女少妇八av免费久了| 制服诱惑二区| 欧美日韩视频精品一区| 国产伦人伦偷精品视频| 啦啦啦 在线观看视频| 国产伦理片在线播放av一区| 日韩视频在线欧美| 日本91视频免费播放| 亚洲 国产 在线| 性色av一级| 亚洲五月色婷婷综合| 一区二区三区激情视频| 另类精品久久| 男女午夜视频在线观看| 超碰97精品在线观看| 少妇 在线观看| 丰满少妇做爰视频| 日日爽夜夜爽网站| 91精品伊人久久大香线蕉| 国产视频首页在线观看| 国产一区亚洲一区在线观看| 人人妻人人澡人人看| 久久99精品国语久久久| 亚洲七黄色美女视频| 国产伦人伦偷精品视频| 中文字幕另类日韩欧美亚洲嫩草| 欧美黄色片欧美黄色片| 丁香六月天网| 国产精品麻豆人妻色哟哟久久| a 毛片基地| 日本欧美国产在线视频| 日韩 亚洲 欧美在线| 18禁观看日本| 考比视频在线观看| 美女午夜性视频免费| 精品亚洲成a人片在线观看| 嫩草影视91久久| 悠悠久久av| 免费看av在线观看网站| 五月天丁香电影| 美国免费a级毛片| 一本综合久久免费| 91九色精品人成在线观看| 精品国产国语对白av| 精品卡一卡二卡四卡免费| av国产精品久久久久影院| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产中文字幕在线视频| 各种免费的搞黄视频| 黄色片一级片一级黄色片| 两性夫妻黄色片| 免费人妻精品一区二区三区视频| 超色免费av| 日韩一本色道免费dvd| 美女大奶头黄色视频| 亚洲九九香蕉| 波野结衣二区三区在线| 国产亚洲欧美在线一区二区| 女性生殖器流出的白浆| 99国产精品免费福利视频| 尾随美女入室| 深夜精品福利| 深夜精品福利| 一本一本久久a久久精品综合妖精| 欧美黑人精品巨大| 日韩精品免费视频一区二区三区| 国产淫语在线视频| 国产xxxxx性猛交| 美女午夜性视频免费| 欧美在线黄色| 成年动漫av网址| 两个人免费观看高清视频| 高清视频免费观看一区二区| 啦啦啦 在线观看视频| 国产精品九九99| 亚洲国产精品一区二区三区在线| 亚洲激情五月婷婷啪啪| 欧美精品av麻豆av| 老汉色∧v一级毛片| 亚洲av日韩精品久久久久久密 | 热99国产精品久久久久久7| 国产真人三级小视频在线观看| 1024视频免费在线观看| 精品国产国语对白av| 日本欧美视频一区| 如日韩欧美国产精品一区二区三区| 久久ye,这里只有精品| 蜜桃国产av成人99| 国产男女超爽视频在线观看| 男女边吃奶边做爰视频| 国产免费视频播放在线视频| 精品国产一区二区三区四区第35| 制服诱惑二区| 夜夜骑夜夜射夜夜干| 成年动漫av网址| 一区二区三区精品91| 亚洲精品av麻豆狂野| 国产亚洲av片在线观看秒播厂| 人成视频在线观看免费观看| 色视频在线一区二区三区| 亚洲av在线观看美女高潮| 黄频高清免费视频| 精品一区在线观看国产| 亚洲av在线观看美女高潮| 高清黄色对白视频在线免费看| 久久影院123| 欧美 日韩 精品 国产| 叶爱在线成人免费视频播放| 欧美大码av| 国产爽快片一区二区三区| 叶爱在线成人免费视频播放| 日韩av免费高清视频| 国产精品一区二区精品视频观看| 大香蕉久久网| 男女国产视频网站| 亚洲图色成人| 亚洲五月婷婷丁香| 女性被躁到高潮视频| 热99久久久久精品小说推荐| 亚洲免费av在线视频| 一级a爱视频在线免费观看| 国产不卡av网站在线观看| 精品国产乱码久久久久久男人| 两人在一起打扑克的视频| 欧美日韩av久久| 精品人妻熟女毛片av久久网站| 亚洲伊人久久精品综合| 欧美精品亚洲一区二区| 午夜福利在线免费观看网站| 丰满迷人的少妇在线观看| 狠狠精品人妻久久久久久综合| 你懂的网址亚洲精品在线观看| 多毛熟女@视频| 欧美日韩亚洲高清精品| 国产精品久久久久久人妻精品电影 | 在线天堂中文资源库| av片东京热男人的天堂| 国产亚洲一区二区精品| 大型av网站在线播放| 欧美大码av| 晚上一个人看的免费电影| 色94色欧美一区二区| 两人在一起打扑克的视频| 夜夜骑夜夜射夜夜干| av视频免费观看在线观看| 91国产中文字幕| xxxhd国产人妻xxx| 乱人伦中国视频| 男女之事视频高清在线观看 | 久久人人爽人人片av| 深夜精品福利| 亚洲av男天堂| 日韩中文字幕欧美一区二区 | 久久国产亚洲av麻豆专区| 精品人妻一区二区三区麻豆| 免费观看人在逋| 亚洲精品乱久久久久久| 午夜两性在线视频| 欧美日韩亚洲综合一区二区三区_| 欧美乱码精品一区二区三区| 涩涩av久久男人的天堂| 青草久久国产| 精品国产国语对白av| 精品一区二区三卡| 中国美女看黄片| 国产99久久九九免费精品| 99久久人妻综合| 欧美激情 高清一区二区三区| 亚洲激情五月婷婷啪啪| 久久久久国产一级毛片高清牌| 精品一区二区三区av网在线观看 | 五月开心婷婷网| a级毛片在线看网站| 爱豆传媒免费全集在线观看| 一区二区三区四区激情视频| 免费日韩欧美在线观看| 老司机靠b影院| 久久狼人影院| 亚洲视频免费观看视频| 赤兔流量卡办理| 亚洲精品美女久久久久99蜜臀 | 国产精品秋霞免费鲁丝片| 无限看片的www在线观看| 国产欧美日韩一区二区三 | 亚洲中文av在线| 亚洲国产中文字幕在线视频| 自线自在国产av| 婷婷色麻豆天堂久久| 王馨瑶露胸无遮挡在线观看| 看免费成人av毛片| 日本午夜av视频| 亚洲欧美清纯卡通| 黄色视频不卡| 中文字幕最新亚洲高清| 国产在线一区二区三区精| 欧美亚洲 丝袜 人妻 在线| 夜夜骑夜夜射夜夜干| 欧美黑人欧美精品刺激| 久久久国产一区二区| 一级a爱视频在线免费观看| 国产高清不卡午夜福利| 高清不卡的av网站| 久久精品熟女亚洲av麻豆精品| 久久久亚洲精品成人影院| 国产成人av激情在线播放| 亚洲欧美日韩另类电影网站| 日韩免费高清中文字幕av| 国产高清视频在线播放一区 | 国产精品国产av在线观看| 国产欧美日韩综合在线一区二区| 久久国产亚洲av麻豆专区| 黄色毛片三级朝国网站| 久久精品aⅴ一区二区三区四区| av福利片在线| 啦啦啦中文免费视频观看日本| 丁香六月天网| 免费黄频网站在线观看国产| 免费观看人在逋| 两个人看的免费小视频| 亚洲欧美日韩另类电影网站| 国产三级黄色录像| 满18在线观看网站| 亚洲成国产人片在线观看| 一本综合久久免费| 老司机深夜福利视频在线观看 | 欧美亚洲日本最大视频资源| 久久精品亚洲熟妇少妇任你| 一区二区日韩欧美中文字幕| 日韩av不卡免费在线播放| 91国产中文字幕| 日韩免费高清中文字幕av| 成年动漫av网址| 亚洲一码二码三码区别大吗| 欧美日韩亚洲国产一区二区在线观看 | 男人舔女人的私密视频| 在线天堂中文资源库| 1024视频免费在线观看| 视频区图区小说| cao死你这个sao货| 国产精品三级大全| 欧美日韩av久久| 亚洲伊人色综图| 精品人妻熟女毛片av久久网站| 亚洲伊人久久精品综合| 欧美日韩av久久| 国产精品久久久av美女十八| 香蕉丝袜av| 视频在线观看一区二区三区| 黄色怎么调成土黄色| 久久久久久久国产电影| 久久久欧美国产精品| 如日韩欧美国产精品一区二区三区| 伦理电影免费视频| 免费看av在线观看网站| 午夜福利,免费看| 一区二区三区精品91| 欧美国产精品一级二级三级| 我要看黄色一级片免费的| 日本91视频免费播放| 视频在线观看一区二区三区| 最新的欧美精品一区二区| 黄色毛片三级朝国网站| 亚洲久久久国产精品| √禁漫天堂资源中文www| 欧美精品高潮呻吟av久久| 最黄视频免费看| 亚洲成人免费电影在线观看 | 国产亚洲一区二区精品| 国产黄色免费在线视频| 久久久久国产精品人妻一区二区| 纵有疾风起免费观看全集完整版| 亚洲国产精品成人久久小说| 91精品三级在线观看| 狂野欧美激情性xxxx| 亚洲午夜精品一区,二区,三区| 日本av手机在线免费观看| 新久久久久国产一级毛片| 亚洲av国产av综合av卡| 91精品伊人久久大香线蕉| 午夜影院在线不卡| 制服诱惑二区| 丝袜脚勾引网站| 国产欧美亚洲国产| 精品第一国产精品| 国产亚洲一区二区精品| 欧美激情极品国产一区二区三区| 精品久久久久久电影网| a级毛片黄视频| 考比视频在线观看| 首页视频小说图片口味搜索 | 亚洲一码二码三码区别大吗| 丁香六月欧美| 精品亚洲成国产av| 日韩av不卡免费在线播放| 91精品伊人久久大香线蕉| tube8黄色片| 操出白浆在线播放| 亚洲色图综合在线观看| 日本欧美视频一区| 日本色播在线视频| 美国免费a级毛片| 国产av一区二区精品久久| 汤姆久久久久久久影院中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 国产激情久久老熟女| 久久久久久免费高清国产稀缺| 美女高潮到喷水免费观看| 国产成人一区二区三区免费视频网站 | 久久久亚洲精品成人影院| 亚洲中文日韩欧美视频| 在线观看免费视频网站a站| 一级毛片 在线播放| 国产精品熟女久久久久浪| 国产伦理片在线播放av一区| 18禁裸乳无遮挡动漫免费视频| 久久久精品区二区三区| 日本a在线网址| 久久久久久久精品精品| 国产野战对白在线观看| 亚洲色图 男人天堂 中文字幕| 国产精品熟女久久久久浪| 波野结衣二区三区在线| 日日爽夜夜爽网站| 国产91精品成人一区二区三区 | 亚洲久久久国产精品| 男女边吃奶边做爰视频| 日本五十路高清| 亚洲色图 男人天堂 中文字幕| 一边摸一边抽搐一进一出视频| 大香蕉久久网| 精品一区二区三区av网在线观看 | 成人18禁高潮啪啪吃奶动态图| 日韩 欧美 亚洲 中文字幕| 男人操女人黄网站| 久久99热这里只频精品6学生| av在线app专区| 免费看十八禁软件| 一本大道久久a久久精品| 国产精品.久久久| 久久久精品94久久精品| 精品人妻一区二区三区麻豆| 香蕉国产在线看| 亚洲精品一卡2卡三卡4卡5卡 | 欧美在线黄色| 亚洲熟女精品中文字幕| 老司机影院成人| 人体艺术视频欧美日本| 久久久久久久精品精品| 人妻一区二区av| 久久久久国产一级毛片高清牌| 久久久久久免费高清国产稀缺| 国产伦人伦偷精品视频| 热re99久久国产66热| 成年动漫av网址| 一本色道久久久久久精品综合| 欧美成人午夜精品| 国产有黄有色有爽视频| 免费一级毛片在线播放高清视频 | 在线精品无人区一区二区三| 成人国产av品久久久| 91精品三级在线观看| 日韩av免费高清视频| 亚洲精品国产区一区二| cao死你这个sao货| 9191精品国产免费久久| 国产淫语在线视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲成人免费电影在线观看 | www日本在线高清视频| 亚洲久久久国产精品| 91麻豆av在线| 好男人电影高清在线观看| 亚洲黑人精品在线| 欧美精品av麻豆av| 中文字幕亚洲精品专区| 亚洲一码二码三码区别大吗| 看免费av毛片| 国产又爽黄色视频| 亚洲av成人精品一二三区| 九色亚洲精品在线播放| 亚洲精品美女久久久久99蜜臀 | 亚洲少妇的诱惑av| 男的添女的下面高潮视频| 精品少妇久久久久久888优播| 日韩 欧美 亚洲 中文字幕| 国产成人av教育| 亚洲成人免费电影在线观看 | 欧美 亚洲 国产 日韩一| 国产精品.久久久| 亚洲欧洲精品一区二区精品久久久| 免费黄频网站在线观看国产| 亚洲天堂av无毛| 曰老女人黄片| 性高湖久久久久久久久免费观看| 亚洲精品国产色婷婷电影| 国产一区有黄有色的免费视频| 麻豆乱淫一区二区| 中文字幕人妻丝袜一区二区| 精品欧美一区二区三区在线| 97在线人人人人妻| 亚洲熟女精品中文字幕| 97精品久久久久久久久久精品| 精品国产一区二区三区久久久樱花| 亚洲 欧美一区二区三区| 久久毛片免费看一区二区三区| 国产精品人妻久久久影院| 婷婷成人精品国产| 国产伦人伦偷精品视频| 中文字幕av电影在线播放| 少妇被粗大的猛进出69影院| 各种免费的搞黄视频| av线在线观看网站| 香蕉国产在线看| 99热全是精品| 在线看a的网站| 久久 成人 亚洲| 一本一本久久a久久精品综合妖精| 操美女的视频在线观看| 免费在线观看视频国产中文字幕亚洲 | 黄色a级毛片大全视频| 精品一品国产午夜福利视频| 极品人妻少妇av视频| 精品福利观看| 久久人人爽人人片av| 国产精品九九99| 午夜免费鲁丝| 国产精品国产av在线观看| 大片电影免费在线观看免费| 少妇粗大呻吟视频| 男女边摸边吃奶| 国产片特级美女逼逼视频| 久久人妻熟女aⅴ| 亚洲精品国产av成人精品| 亚洲av日韩精品久久久久久密 | 国产成人a∨麻豆精品| 亚洲男人天堂网一区| 国产亚洲午夜精品一区二区久久| 欧美精品高潮呻吟av久久| 精品久久久久久电影网| 中文字幕制服av| 欧美日韩视频精品一区| 又紧又爽又黄一区二区| 亚洲欧美精品自产自拍| 咕卡用的链子| 欧美日韩综合久久久久久| 91精品国产国语对白视频| 亚洲精品美女久久久久99蜜臀 | 1024香蕉在线观看| 久久久国产一区二区| 久久狼人影院| 一本色道久久久久久精品综合| 妹子高潮喷水视频| 又大又黄又爽视频免费| 亚洲av成人不卡在线观看播放网 | bbb黄色大片| 亚洲av成人不卡在线观看播放网 | 少妇猛男粗大的猛烈进出视频| 最近中文字幕2019免费版| 色婷婷久久久亚洲欧美| 考比视频在线观看| 国产1区2区3区精品| 性色av乱码一区二区三区2| 亚洲欧美日韩另类电影网站| 亚洲国产av影院在线观看| 搡老岳熟女国产| 国产一区二区 视频在线| 19禁男女啪啪无遮挡网站| 免费日韩欧美在线观看| 亚洲一区二区三区欧美精品| 亚洲av成人不卡在线观看播放网 | 久久青草综合色| 国产一区二区三区综合在线观看| 日韩制服骚丝袜av| 久久热在线av| 久久国产精品男人的天堂亚洲| 极品人妻少妇av视频| 亚洲av成人不卡在线观看播放网 | 国产精品成人在线| 水蜜桃什么品种好| 国产精品av久久久久免费| 丝袜美足系列| 亚洲av日韩精品久久久久久密 | 亚洲国产欧美日韩在线播放| 欧美日韩视频精品一区| 亚洲精品国产区一区二| 一区在线观看完整版| 欧美 日韩 精品 国产| 美女主播在线视频| 午夜精品国产一区二区电影| 最新的欧美精品一区二区| 18禁观看日本| 日本猛色少妇xxxxx猛交久久| www.av在线官网国产| 国产精品熟女久久久久浪| 久久性视频一级片| av在线app专区| 久久国产精品人妻蜜桃| 国产一区二区激情短视频 | 赤兔流量卡办理| 黄色毛片三级朝国网站| 欧美精品一区二区大全| 另类精品久久| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品一二三| 国产成人精品久久二区二区免费| 国产黄色免费在线视频| 精品人妻1区二区| 成人三级做爰电影| 女人久久www免费人成看片| 国产片特级美女逼逼视频| 一区在线观看完整版| 男人添女人高潮全过程视频| 亚洲精品久久午夜乱码| 久久99一区二区三区| 欧美日韩亚洲高清精品| 成人三级做爰电影| 女性生殖器流出的白浆| 国产精品99久久99久久久不卡| 中文欧美无线码| 国产xxxxx性猛交| 久久青草综合色| 国产av精品麻豆| 在线观看www视频免费| e午夜精品久久久久久久| 日韩电影二区| av在线播放精品| 人成视频在线观看免费观看| 三上悠亚av全集在线观看| 欧美另类一区| 亚洲成色77777| 亚洲三区欧美一区| 亚洲国产av新网站| 少妇被粗大的猛进出69影院| 一本久久精品| 亚洲精品久久久久久婷婷小说| 国产深夜福利视频在线观看| 国产高清国产精品国产三级| 国产精品一区二区在线观看99| 日韩中文字幕欧美一区二区 | 午夜免费观看性视频| 免费一级毛片在线播放高清视频 | 90打野战视频偷拍视频| 亚洲专区国产一区二区| 日韩大片免费观看网站| 欧美黑人欧美精品刺激| 久久精品亚洲熟妇少妇任你| 欧美精品人与动牲交sv欧美| 欧美精品亚洲一区二区| 亚洲精品美女久久久久99蜜臀 | svipshipincom国产片| 一区福利在线观看| 国产片内射在线| 热99久久久久精品小说推荐| 久久综合国产亚洲精品| 久久鲁丝午夜福利片| 欧美少妇被猛烈插入视频| 大香蕉久久网| 国产精品熟女久久久久浪| 欧美亚洲 丝袜 人妻 在线| 一区二区三区激情视频| 少妇人妻 视频| 午夜福利一区二区在线看| 精品少妇久久久久久888优播| 黄色一级大片看看| 成年动漫av网址| 国产高清国产精品国产三级| 久久女婷五月综合色啪小说| 久久久久网色| 欧美日韩一级在线毛片| 丝袜喷水一区|