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

    曲拐彎扭耦合振動影響規(guī)律研究

    2016-10-21 11:41:34陳宏威向建華吳斌輝北京理工大學(xué)機(jī)械與車輛學(xué)院北京100081
    關(guān)鍵詞:成正比內(nèi)燃機(jī)圓盤

    陳宏威 向建華 劉 帥 吳斌輝(北京理工大學(xué)機(jī)械與車輛學(xué)院北京100081)

    ·振動·噪聲·

    曲拐彎扭耦合振動影響規(guī)律研究

    陳宏威向建華劉帥吳斌輝
    (北京理工大學(xué)機(jī)械與車輛學(xué)院北京100081)

    基于圓盤非對稱轉(zhuǎn)子簡化模型和達(dá)朗貝爾原理,建立了曲軸單拐彎-扭耦合非線性振動數(shù)學(xué)模型,并利用解析推導(dǎo)法和四階龍格庫塔的迭代數(shù)值求解方法研究了彎曲振動和扭轉(zhuǎn)振動間的相互影響規(guī)律。研究結(jié)果表明,彎曲振動(或扭轉(zhuǎn)振動)會耦合產(chǎn)生多種頻率成分的扭轉(zhuǎn)振動(或彎曲振動),耦合振動幅值隨著激勵頻率和偏心距的增加而加大,且彎曲振動對扭轉(zhuǎn)振動的耦合影響相對較大。

    曲拐彎扭耦合振動非線性振動影響規(guī)律

    引言

    內(nèi)燃機(jī)曲軸應(yīng)用范圍廣泛,結(jié)構(gòu)復(fù)雜,功能繁多,是內(nèi)燃機(jī)系統(tǒng)關(guān)鍵部件[1]。曲軸承受周期性變化的切向力和法向力作用,常常產(chǎn)生彎曲振動、扭轉(zhuǎn)振動及二者的耦合振動形式,很容易引發(fā)多種振動問題,導(dǎo)致軸段的疲勞損壞、軸承磨損以及工作噪聲等[2]。以往人們關(guān)于內(nèi)燃機(jī)曲軸的研究主要集中于軸的扭轉(zhuǎn)振動,但近年來隨著認(rèn)識的深入,一些學(xué)者開始著重關(guān)注扭轉(zhuǎn)—彎曲耦合振動,形成了比較完備的理論分析、數(shù)值模擬與實驗研究相結(jié)合的體系[3-5]。由于曲軸高速回轉(zhuǎn)時的質(zhì)量不平衡或系統(tǒng)不對稱等非線性因素,使得數(shù)值仿真中曲軸彎曲振動的數(shù)學(xué)建模和求解比較復(fù)雜[6-8],其相應(yīng)的曲軸彎扭耦合振動分析的相關(guān)文獻(xiàn)目前尚不多見。李慧珍[9]在有限元計算曲軸振動中將曲軸視為圓截面梁和變截面矩形梁的組合;梁興雨[10]等建立了由多個自由度組成的發(fā)動機(jī)剛?cè)狁詈隙囿w動力學(xué)系統(tǒng)模型;何成兵[11]推導(dǎo)了三類微元軸段的彎扭耦合振動微分方程組,提出數(shù)值計算解法。Nataraj C[12]利用攝動的方法得到了非線性彎扭耦合振動微分方程組的初步解;梁明軒[13]分別通過諧波平衡法與數(shù)值積分法得到系統(tǒng)幅頻特性,分析了質(zhì)量偏心和轉(zhuǎn)速對軸系非線性動力學(xué)響應(yīng)的影響;劉占生[14]得到了共振區(qū)彎曲振動和扭轉(zhuǎn)振動的頻率特征,分析了不平衡量對彎扭耦合共振區(qū)的影響。目前,對于內(nèi)燃機(jī)曲軸的振動耦合分析,大多是基于計算精度較差的離散質(zhì)量或分布參數(shù)模型,且彎扭耦合振動影響等理論研究不夠深入。

    本文針對某內(nèi)燃機(jī)曲軸單拐,并基于圓盤非對稱轉(zhuǎn)子簡化模型和達(dá)朗貝爾原理,建立了曲拐彎-扭耦合非線性振動數(shù)學(xué)模型,并利用解析推導(dǎo)和龍格庫塔的迭代數(shù)值求解方法深入研究了彎曲振動和扭轉(zhuǎn)振動間的相互影響規(guī)律,從而為內(nèi)燃機(jī)曲軸彎扭耦合研究提供理論指導(dǎo)。

    1 非線性彎扭耦合模型建立和求解

    對于截取的曲拐,在考慮扭轉(zhuǎn)振動的前提下,采用離散化模型分析彎曲振動運動狀況。將曲軸單拐簡化為一個剛性圓盤并安裝在具有彎曲和扭轉(zhuǎn)彈性而無質(zhì)量的軸上。假定軸是各向同性,兩端軸承為剛性支撐,忽略了圓盤的陀螺效應(yīng),進(jìn)而將整個曲拐簡化為一個單質(zhì)量圓盤轉(zhuǎn)子模型,如圖1所示。

    圖1 曲拐和簡化后的單質(zhì)量圓盤轉(zhuǎn)子模型

    建立考慮扭轉(zhuǎn)振動耦合影響和偏心距的單質(zhì)量圓盤轉(zhuǎn)子示意圖如圖2所示。圓盤旋轉(zhuǎn)中心為O,質(zhì)心為O1,偏心距為e,φ為質(zhì)心旋轉(zhuǎn)的角度,ω為轉(zhuǎn)動角頻率,即φ=ωt+θ,θ為扭轉(zhuǎn)角。

    圖2 圓盤轉(zhuǎn)子示意圖

    由于偏心距e的存在,圓盤具有回轉(zhuǎn)速度φ˙和回轉(zhuǎn)加速度,使得質(zhì)心上有切向加速度eφ¨和法向加速度eφ˙2φ¨,如圖3所示。由此產(chǎn)生的慣性力為:

    圖3 慣性力分析示意圖

    同理,質(zhì)心具有彎曲加速度y¨和z¨,產(chǎn)生的慣性力為-my¨和-mz¨,由此施加在回轉(zhuǎn)中心的額外扭矩為:

    在假設(shè)圓盤為剛體且左右支撐軸彎扭各向同性的前提下,根據(jù)達(dá)朗伯原理得到彎扭耦合振動平衡方程組:式中:T為作用在曲柄銷處的切向力產(chǎn)生的扭矩;Fy、Fz為作用在曲柄銷處的徑向力在y、z方向的投影;c和k為彎曲阻尼和剛度,此處不再考慮y和z方向的各向異性;ct和kt為扭轉(zhuǎn)阻尼和剛度。

    從式(3)定性分析可以看出,耦合振動的產(chǎn)生是由于偏心距e的存在,e越大耦合作用就越大。而當(dāng)e=0時,圓盤的彎扭耦合振動平衡方程組變成:

    即當(dāng)偏心距不存在時,彎曲扭轉(zhuǎn)振動間不存在明顯的耦合關(guān)系。

    式(3)所示的彎扭耦合振動平衡方程組中變量z、y、θ均為時間的函數(shù),呈現(xiàn)高度非線性。本文利用Matlab/Simulink搭建式(3)所對應(yīng)的彎扭耦合振動平衡方程組仿真模型(如圖4所示),并基于四階-龍格庫塔法(Runge-Kutta)迭代求解,獲得相應(yīng)的數(shù)值解。

    圖4中,左右兩邊是模型輸入和輸出模塊,淺灰色的是兩個彎曲振動和1個扭轉(zhuǎn)振動計算模塊,深灰色的是彎扭和扭彎耦合計算模塊。仿真求解時,將初始值(作用力或力矩)賦給對應(yīng)的振動計算模塊,得到的值經(jīng)由各耦合模塊計算獲得新的耦合力或是力矩,再一并輸入作為總激勵來計算下一輪的值。以此反復(fù)迭代求解,直至滿足收斂條件,從而實現(xiàn)彎扭耦合振動平衡方程組的數(shù)值求解。

    由于彎扭耦合影響是相互的,本文在分析彎扭耦合振動作用時,采用單獨控制一種振動而激發(fā)另一種振動的辦法,即通過單獨施加力或是力矩,觀察一種振動產(chǎn)生時另一種振動頻域下的幅值變化,從而獲得彎扭振動之間的耦合作用規(guī)律。

    圖4 彎扭耦合振動求解模型

    2 彎曲振動對扭轉(zhuǎn)振動的影響研究

    由式(3)可以看出,由于偏心距e的存在,彎曲振動在圓盤上產(chǎn)生了額外扭矩M,即M=me(y¨sinφz¨cosφ),而重力作用產(chǎn)生的額外扭矩大小為megsin(tω+θ)。對M進(jìn)行傅立葉變換,將y、z方向的彎曲振動分解為若干簡諧振動的疊加,得到:

    式中:ωi為第i個彎曲振動分量的振動角頻率;Yi、Zi分別為y、z兩個方向上振動的振幅;αi、βi為各自的相位。

    將第i分量代入M=me(y¨sinφ-z¨cosφ),得到頻率為ωi的額外扭矩:

    從式(6)可以看出,角頻率為ωi的彎曲振動可以引起角頻率為ωi±ω的扭轉(zhuǎn)振動;彎曲振動產(chǎn)生的額外扭矩的幅值與不平衡量me成正比,與彎曲振幅(Yi,Zi)成正比,與該諧分量下角頻率的平方成正比,彎曲振動相位角(αi,βi)對扭矩的幅值影響也較大。此外,由于重力作用也會產(chǎn)生額外扭矩,其大小為meg sin(ωt+θ),故重力會引起角頻率為ω的扭轉(zhuǎn)振動。

    取單位曲拐模型仿真計算的基本參數(shù)為:轉(zhuǎn)速為3 000 r/min,即轉(zhuǎn)動頻率為50 Hz,k=1.736×1010N/m,J=0.045 kg·m2,kt=1.285×106(N·m)/rad,e= 0.02 m,m=9.07 kg,彎曲激勵頻率為35 Hz,幅值大小為Fy=1.0×107N,F(xiàn)z=5×105N,c=80(N·s)/m,Ct=100(N·m)/rad,通過圖4的數(shù)值迭代計算,結(jié)果如圖5所示。

    從圖5可以看出,在激勵力35 Hz作用下,彎曲振動在35 Hz產(chǎn)生了明顯的波峰。同時,在扭轉(zhuǎn)振動的15 Hz、50 Hz、85 Hz附近產(chǎn)生較為明顯的峰值;耦合作用產(chǎn)生的扭轉(zhuǎn)振動幅值和彎曲振動幅值比(無量綱)為1∶40。扭轉(zhuǎn)振動在15 Hz和85 Hz下產(chǎn)生的峰值符合│ω±ωi│的規(guī)律,即為激勵力頻率與轉(zhuǎn)動頻率之差或之和的絕對值。50 Hz處的峰值產(chǎn)生是因為重力所導(dǎo)致的不平衡引起的,重力可以產(chǎn)生meg sin(ωt+θ)的額外扭矩進(jìn)而引起頻率為ω的扭轉(zhuǎn)振動。值得注意的是,雖然由彎曲振動產(chǎn)生的扭矩M通過傅里葉變換可以得出N個不同角頻率的簡諧振動,進(jìn)而引起相應(yīng)的一系列的扭轉(zhuǎn)振動;但只有峰值較為明顯(通常是低階)的彎曲振動,耦合產(chǎn)生的扭轉(zhuǎn)振動才較為明顯。其余情況下產(chǎn)生的扭轉(zhuǎn)振動可以忽略。

    此外,在轉(zhuǎn)速一定的情況下,隨著彎曲振動激勵力頻率從55 Hz到90 Hz的變化,扭轉(zhuǎn)振動幅值也呈一定的變化。其中圖6顯示了│ω-ωi│扭轉(zhuǎn)振動隨著激勵力ω變化關(guān)系。

    從圖6中可以看出,隨著彎曲振動激勵力頻率從55 Hz到90 Hz的變化,頻率│ω-ωi│扭轉(zhuǎn)振動幅值逐漸增加。因為根據(jù)公式(6)可知,彎曲振動耦合產(chǎn)生的額外扭矩大小的決定因素中有激勵力頻率,且與其平方成正比,同時受重力產(chǎn)生的旋轉(zhuǎn)扭矩的影響,但總體上扭轉(zhuǎn)振動會隨著激勵力頻率的增加而相應(yīng)的增加。特別指出,當(dāng)扭轉(zhuǎn)固有頻率和彎曲固有頻率都不大時,會出現(xiàn)由于彎曲振動耦合作用所導(dǎo)致的扭轉(zhuǎn)共振的情況發(fā)生。

    圖5 彎曲振動和耦合產(chǎn)生的扭轉(zhuǎn)振動頻域圖

    圖6 │ω-ωi│扭轉(zhuǎn)振動隨著激勵力ω變化關(guān)系

    3 扭轉(zhuǎn)振動對彎曲振動的影響研究

    由式(3)可以看出,由于存在偏心距和扭轉(zhuǎn)角的原因,在圓盤部產(chǎn)生了切向力和法向力,這些激勵力作用在彎曲振動方向,使其產(chǎn)生相應(yīng)的彎曲振動。這里,以Z方向彎曲振動為研究對象;為進(jìn)一步分析微分方程,可以把扭轉(zhuǎn)角θ分解為若干諧振動的和,即

    扭對彎耦合項Z方向力me(φ˙2sinφ-φ¨cosφ)進(jìn)行分析,即:

    式中:第1項是由于旋轉(zhuǎn)不平衡引起的,為傳統(tǒng)的不平衡彎曲振動響應(yīng)激勵力,它將引起同頻率的彎曲振動。后邊各項均是由于扭轉(zhuǎn)振動引起的離心激勵力,頻率為:

    1)ω±ωi,幅值與、不平衡me、轉(zhuǎn)速ω成正比,與ωi或是成正比。

    3)ω±ωi±ωj,幅值與乘積、不平衡me、ωiωj乘積成正比。

    4)2ωi±ω,幅值與、不平衡me成正比。

    通過對比可知,這些激勵力由于幅值引入了扭振角,扭轉(zhuǎn)振動引起的慣性力相比不平衡直接引起的同頻率彎曲振動幅值要小很多;而其中ω±ωi頻率下對應(yīng)的振幅相對大一些。

    同樣,取單位曲拐模型仿真計算的基本參數(shù)為:轉(zhuǎn)速為3 000 r/min,即轉(zhuǎn)動頻率為50 Hz,k=1.736× 1010N/m,J=0.045 kg·m2,kt=1.285×106(N·m)/rad,e=0.02 m,m=9.07 kg,扭矩頻率為35 Hz,幅值大小為T=2e4N·m,c=80(N·s)·m-1,ct=100(N·m·s)·rad-1。通過圖4的數(shù)值迭代計算,結(jié)果如圖7~圖9所示。

    圖7 扭轉(zhuǎn)振動頻域圖

    圖8 耦合引起的同頻率彎曲振動幅值

    圖9 耦合引起的ω±ωi頻率彎曲振動幅值

    從圖7可以看出,在35 Hz扭矩的作用下,產(chǎn)生了比較明顯的波峰;圖9中彎曲振動在15 Hz、50 Hz和85Hz附近產(chǎn)生了明顯的波峰;耦合作用產(chǎn)生的彎曲振動幅值和扭轉(zhuǎn)振動幅值比(無量綱)為1∶12 000。圖8中在彎曲振動50 Hz處產(chǎn)生明顯峰值的原因是傳統(tǒng)的同頻彎曲振動作用的結(jié)果,其幅值與不平衡成正比,與轉(zhuǎn)速ω平方成正比。圖9中彎曲振動在15 Hz和85 Hz處產(chǎn)生峰值的原因與│ω±ωi│相符合,即是轉(zhuǎn)動角頻率與扭轉(zhuǎn)振動頻率之差或之和,其幅值與不平衡me成正比,與轉(zhuǎn)速ω成正比,與ωi或是ω2i成正比,與θi成正比。由于這里的扭轉(zhuǎn)剛度很大,故θi一般比較小,使得幅值要遠(yuǎn)小于同頻率引起的彎曲振動;而其它頻率成分對應(yīng)的幅值更小,耦合作用不明顯,可以忽略。

    圖10和圖11描述了在軸系轉(zhuǎn)速頻率為50 Hz,外扭矩頻率從55 Hz逐漸增加到90 Hz時,扭轉(zhuǎn)振動引起的│ω±ωi│頻率的彎曲振動隨著外扭矩頻率的變化趨勢。隨著扭轉(zhuǎn)激勵力頻率的增加,彎曲振動逐漸增加,這與彎曲振動對扭轉(zhuǎn)振動的結(jié)論相似。

    圖10 Y向彎曲振動隨扭振頻率的變化

    圖11 Z向彎曲振動隨扭振頻率的變化

    4 結(jié)論

    本文建立了曲軸單拐彎-扭耦合振動的數(shù)學(xué)模型,并針對彎-扭耦合振動和扭-彎耦合振動的影響規(guī)律進(jìn)行了解析推導(dǎo)和數(shù)值迭代求解,計算結(jié)果具有很好的一致性。具體結(jié)論如下:

    1)彎曲振動(或扭轉(zhuǎn)振動)會耦合產(chǎn)生多種頻率成分的扭轉(zhuǎn)振動(或彎曲振動),都包括頻率為ω±ωi的振動成分。

    2)耦合產(chǎn)生的振動幅值隨著激勵頻率和偏心距的增加而加大。

    3)無量綱下的彎扭耦合幅值比遠(yuǎn)大于扭彎耦合幅值比,故本曲軸單拐的彎曲振動對扭轉(zhuǎn)振動的耦合影響相對較大。

    1Fuchun Ren,Kun Yang,Zhong Guang Fu,et al.Study on the coupling of the bending and vibrations of the unbalanced rotor[C]:Xi'an:Proceedings of the 2nd International Conference on Hydrodynamic Bearing-rotor System Dynamics,March 20-22,1997

    2李震,桂長林,孫軍.內(nèi)燃機(jī)曲軸軸系振動分析研究的現(xiàn)狀、討論與展望[J].內(nèi)燃機(jī)學(xué)報,2002,20(05):469-474

    3Nevsan S.Performance increase in turbo molecular pumps with curved type blades[J].Vacuum,2012,86(11):1764-1769

    4韓建鑫,王煒,張琪昌.內(nèi)燃機(jī)曲軸系統(tǒng)變慣量公式的修正[J].內(nèi)燃機(jī)學(xué)報,2014,32(03):271-275

    5Munawar I,Abdul W,Dimitri B,et al.Design modification in rotor blade of turbo molecular pump[J].Nuclear Instruments and Methods in Physics Research A,2012,678:88-90

    6張勇,蔣滋康.軸系彎扭耦合振動的數(shù)學(xué)模型[J].清華大學(xué)學(xué)報(自然科學(xué)版),1998(08):115-118

    7周桐,徐健學(xué),傅衛(wèi)平.彎扭耦合振動對次同步諧振響應(yīng)的影響[J].應(yīng)用力學(xué)學(xué)報,2000,17(1):12-17

    8王雪嬌.內(nèi)燃機(jī)曲軸彎扭耦合非線性振動分析[D].天津:天津大學(xué),2009

    9李慧珍,張德平.用有限元進(jìn)行曲軸扭轉(zhuǎn)振動計算[J].內(nèi)燃機(jī)學(xué)報,1991,9(2):157-162

    10梁興雨,舒歌群,李東海,等.基于柔體曲軸多體動力學(xué)的軸系扭振響應(yīng)分析[J].內(nèi)燃機(jī)工程,2007(04):46-49

    11何成兵,楊昆,顧煜炯.質(zhì)量偏心對碰磨轉(zhuǎn)子彎曲振動和扭轉(zhuǎn)振動特性的影響[J].中國電機(jī)工程學(xué)報,2002,22(7):105-110

    12 Nataraj C.On the Interaction of Torsion and Bending in Rotating Shafts[J].Journal of Applied Mechanics,1993,60(3):239-241

    13梁明軒,袁惠群,趙天宇,等.內(nèi)燃機(jī)曲軸軸系彎扭耦合非線性振動響應(yīng)[J].東北大學(xué)學(xué)報(自然科學(xué)版),2014(03):402-405,414

    14劉占生,崔穎,黃文虎,等.轉(zhuǎn)子彎扭耦合非線性動力學(xué)特性研究[J].中國機(jī)械工程,2003,14(7):603-605

    The Nonlinear Study on Coupled Bending and Torsion Vibration of Crank Throw

    Chen Hongwei,Xiang Jianhua,Liu Shuai,Wu Binhui
    School of Mechanical Engineering,Beijing Institute of Technology(Beijing,100081,China)

    Based on the simplified model of asymmetric disc rotor and D'alembert's principle,the nonlinear mathematical model of coupled bending and torsion vibration is established,and then the mutual influence between the bending vibration and torsion vibration is studied through the theoretical derivation and the solving methods of iterative numerical of fourth-order Runge-Kutta.The results show that the bending vibration(or torsion vibration)will produce various frequency of torsion vibration(or bending vibration)if coupled.The coupling vibration amplitude will raise with the increasing of the excitation frequency and eccentricity,the bending vibration has a relatively large coupling influence on torsion vibration.

    Crank throw,Coupled bending and torsion vibration,Nonlinear vibration,Influence rule

    TK422

    A

    2095-8234(2016)04-0065-06

    2016-02-03)

    陳宏威(1991-),男,碩士研究生,主要研究方向為內(nèi)燃機(jī)結(jié)構(gòu)振動。

    向建華(1976-),男,博士,副教授。

    猜你喜歡
    成正比內(nèi)燃機(jī)圓盤
    圓盤鋸刀頭的一種改進(jìn)工藝
    石材(2020年6期)2020-08-24 08:27:00
    內(nèi)燃機(jī)的興與衰
    最高效和超低排放
    ——內(nèi)燃機(jī)4.0 Highest Efficiency and Ultra Low Emission–Internal Combustion Engine 4.0
    汽車文摘(2018年11期)2018-10-30 02:32:34
    單位圓盤上全純映照模的精細(xì)Schwarz引理
    奇怪的大圓盤
    一種證明引力大小與質(zhì)量乘積成正比的數(shù)學(xué)方法*
    物理通報(2016年7期)2016-10-25 02:31:32
    基于Profibus-DP的圓盤澆鑄控制系統(tǒng)的應(yīng)用
    內(nèi)燃機(jī)再制造產(chǎn)業(yè)已初具規(guī)模
    語絲
    讀者(2014年9期)2014-05-14 11:40:52
    《國外內(nèi)燃機(jī)》2014年總索引
    一进一出抽搐gif免费好疼| 国内少妇人妻偷人精品xxx网站 | 免费在线观看日本一区| 黑人操中国人逼视频| 久久精品亚洲精品国产色婷小说| 欧美日韩精品网址| 亚洲全国av大片| 精品人妻1区二区| 久久久国产欧美日韩av| 国产精品二区激情视频| 午夜亚洲福利在线播放| 成年人黄色毛片网站| 免费在线观看黄色视频的| 性欧美人与动物交配| 久久久久久亚洲精品国产蜜桃av| 国产精品电影一区二区三区| 国语自产精品视频在线第100页| 在线观看66精品国产| 动漫黄色视频在线观看| 午夜福利高清视频| www.999成人在线观看| 两个人视频免费观看高清| 久久性视频一级片| 色综合欧美亚洲国产小说| 精品久久久久久久久久久久久 | 亚洲欧美日韩无卡精品| 久久久国产精品麻豆| 麻豆av在线久日| 色精品久久人妻99蜜桃| 免费在线观看成人毛片| 欧洲精品卡2卡3卡4卡5卡区| 成人一区二区视频在线观看| 女同久久另类99精品国产91| 亚洲 国产 在线| 国产精品,欧美在线| 国产黄色小视频在线观看| 亚洲色图 男人天堂 中文字幕| 法律面前人人平等表现在哪些方面| 午夜久久久久精精品| 亚洲中文日韩欧美视频| 精品不卡国产一区二区三区| 国产在线观看jvid| 亚洲狠狠婷婷综合久久图片| 亚洲人成伊人成综合网2020| 夜夜躁狠狠躁天天躁| 日本熟妇午夜| 18禁观看日本| 亚洲精品色激情综合| 国产一区二区在线av高清观看| 叶爱在线成人免费视频播放| 在线观看日韩欧美| 国产av一区在线观看免费| 成人欧美大片| 黄网站色视频无遮挡免费观看| 久久久国产欧美日韩av| 成人免费观看视频高清| 99在线视频只有这里精品首页| 精品久久蜜臀av无| 精品高清国产在线一区| 老汉色∧v一级毛片| 国产又色又爽无遮挡免费看| 亚洲黑人精品在线| 在线看三级毛片| 国产精品一区二区免费欧美| 亚洲色图av天堂| 精品人妻1区二区| 亚洲成av片中文字幕在线观看| 免费搜索国产男女视频| 久久精品国产99精品国产亚洲性色| 最近最新中文字幕大全电影3 | 天堂√8在线中文| 麻豆av在线久日| 91麻豆av在线| 欧美乱色亚洲激情| 亚洲专区中文字幕在线| 1024香蕉在线观看| 日韩大码丰满熟妇| 国产aⅴ精品一区二区三区波| 亚洲色图av天堂| 亚洲欧美日韩无卡精品| 精品国产乱码久久久久久男人| 在线观看日韩欧美| 亚洲天堂国产精品一区在线| 亚洲精品一卡2卡三卡4卡5卡| 99国产极品粉嫩在线观看| АⅤ资源中文在线天堂| 又大又爽又粗| 国产久久久一区二区三区| netflix在线观看网站| 欧美一区二区精品小视频在线| 宅男免费午夜| 亚洲一卡2卡3卡4卡5卡精品中文| 久久人妻av系列| 白带黄色成豆腐渣| 国内揄拍国产精品人妻在线 | 国产aⅴ精品一区二区三区波| 一本综合久久免费| av天堂在线播放| 亚洲国产欧洲综合997久久, | 欧美成人午夜精品| 日本 av在线| 日韩视频一区二区在线观看| 久久这里只有精品19| 精品第一国产精品| 麻豆成人av在线观看| 一区二区日韩欧美中文字幕| 啪啪无遮挡十八禁网站| 在线av久久热| www日本在线高清视频| 欧美av亚洲av综合av国产av| 99国产精品一区二区蜜桃av| 熟女少妇亚洲综合色aaa.| 男男h啪啪无遮挡| 在线观看www视频免费| 少妇被粗大的猛进出69影院| 老鸭窝网址在线观看| 国产亚洲欧美98| 国产精品电影一区二区三区| 久久精品人妻少妇| 白带黄色成豆腐渣| 久久久久亚洲av毛片大全| 国产精品二区激情视频| a级毛片在线看网站| 狂野欧美激情性xxxx| 男男h啪啪无遮挡| 亚洲中文字幕一区二区三区有码在线看 | 欧美成人免费av一区二区三区| 啦啦啦韩国在线观看视频| 亚洲熟女毛片儿| 黄色女人牲交| 久久亚洲真实| 日本三级黄在线观看| 亚洲av中文字字幕乱码综合 | 亚洲五月婷婷丁香| 在线播放国产精品三级| 在线观看www视频免费| 国产av一区二区精品久久| 亚洲五月天丁香| 欧美一级毛片孕妇| 国产成+人综合+亚洲专区| 国产野战对白在线观看| 人成视频在线观看免费观看| 校园春色视频在线观看| 久久香蕉精品热| 18禁国产床啪视频网站| x7x7x7水蜜桃| 亚洲欧美精品综合一区二区三区| 90打野战视频偷拍视频| 熟女少妇亚洲综合色aaa.| 日本黄色视频三级网站网址| 久久久久久大精品| 日本三级黄在线观看| 亚洲一区二区三区色噜噜| 欧美性猛交黑人性爽| x7x7x7水蜜桃| 精品国产亚洲在线| 欧美成人性av电影在线观看| 精品久久久久久久久久久久久 | 日韩大码丰满熟妇| 亚洲精品av麻豆狂野| 91字幕亚洲| av在线天堂中文字幕| 精品欧美一区二区三区在线| 国产视频一区二区在线看| 一二三四社区在线视频社区8| 日本撒尿小便嘘嘘汇集6| 婷婷亚洲欧美| 免费高清在线观看日韩| 欧美国产日韩亚洲一区| 亚洲中文字幕一区二区三区有码在线看 | 国产一区二区在线av高清观看| 亚洲全国av大片| 亚洲国产欧美日韩在线播放| 999久久久精品免费观看国产| 午夜免费观看网址| 国产精华一区二区三区| 久久久国产成人精品二区| 两个人看的免费小视频| 免费av毛片视频| 神马国产精品三级电影在线观看 | 人妻丰满熟妇av一区二区三区| 真人做人爱边吃奶动态| 国产亚洲精品久久久久久毛片| avwww免费| 亚洲av熟女| 狂野欧美激情性xxxx| 日韩国内少妇激情av| 村上凉子中文字幕在线| 女生性感内裤真人,穿戴方法视频| www.www免费av| 国产久久久一区二区三区| 免费无遮挡裸体视频| 国产成人系列免费观看| av在线天堂中文字幕| 国产97色在线日韩免费| 99国产精品一区二区三区| 97人妻精品一区二区三区麻豆 | 成年女人毛片免费观看观看9| 国产97色在线日韩免费| 亚洲国产精品合色在线| 淫秽高清视频在线观看| 两个人看的免费小视频| 国内揄拍国产精品人妻在线 | 日本黄色视频三级网站网址| 久久久久久久精品吃奶| 亚洲中文日韩欧美视频| 亚洲国产欧洲综合997久久, | 深夜精品福利| 欧美中文综合在线视频| 91字幕亚洲| 日本 欧美在线| 窝窝影院91人妻| 一a级毛片在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 又紧又爽又黄一区二区| 成人午夜高清在线视频 | 51午夜福利影视在线观看| 欧美久久黑人一区二区| avwww免费| 人人妻人人澡人人看| 国产精品1区2区在线观看.| 亚洲av日韩精品久久久久久密| 午夜成年电影在线免费观看| 日韩成人在线观看一区二区三区| 欧美日韩一级在线毛片| 精品一区二区三区视频在线观看免费| 免费人成视频x8x8入口观看| 91九色精品人成在线观看| 香蕉av资源在线| 免费在线观看成人毛片| 久久久久久人人人人人| 大香蕉久久成人网| 国产精品99久久99久久久不卡| 变态另类丝袜制服| 国产成人啪精品午夜网站| 国产精品99久久99久久久不卡| 国产精品免费一区二区三区在线| 女生性感内裤真人,穿戴方法视频| 中文字幕av电影在线播放| 日韩精品中文字幕看吧| 国产在线观看jvid| 99国产综合亚洲精品| 国产视频内射| 看黄色毛片网站| 黄色 视频免费看| 禁无遮挡网站| 一卡2卡三卡四卡精品乱码亚洲| 少妇 在线观看| 在线免费观看的www视频| 1024手机看黄色片| 在线观看舔阴道视频| 亚洲国产精品sss在线观看| 国产成人精品无人区| 日本三级黄在线观看| 香蕉久久夜色| 日韩 欧美 亚洲 中文字幕| 亚洲色图 男人天堂 中文字幕| 成人欧美大片| 亚洲成人久久性| 亚洲狠狠婷婷综合久久图片| 免费在线观看成人毛片| 国产激情偷乱视频一区二区| 欧美+亚洲+日韩+国产| 国产精品一区二区三区四区久久 | 成人亚洲精品av一区二区| 国产视频一区二区在线看| 国产精品,欧美在线| a在线观看视频网站| 国产人伦9x9x在线观看| 亚洲熟女毛片儿| 88av欧美| 操出白浆在线播放| 一本精品99久久精品77| 天天躁夜夜躁狠狠躁躁| 日韩欧美免费精品| 丰满人妻熟妇乱又伦精品不卡| 久久国产精品影院| 88av欧美| e午夜精品久久久久久久| 久久精品亚洲精品国产色婷小说| 男女午夜视频在线观看| a在线观看视频网站| 亚洲在线自拍视频| 日本免费一区二区三区高清不卡| 久久中文字幕人妻熟女| 人妻久久中文字幕网| 精品少妇一区二区三区视频日本电影| 久久99热这里只有精品18| 中国美女看黄片| 怎么达到女性高潮| 午夜两性在线视频| 亚洲精品一区av在线观看| avwww免费| 欧美成人性av电影在线观看| 一个人观看的视频www高清免费观看 | 欧美黑人巨大hd| 久久精品国产99精品国产亚洲性色| 亚洲午夜理论影院| 欧美黑人精品巨大| 97超级碰碰碰精品色视频在线观看| 无限看片的www在线观看| 欧美在线黄色| 啪啪无遮挡十八禁网站| 午夜日韩欧美国产| 18禁观看日本| 一夜夜www| 在线播放国产精品三级| 国产视频一区二区在线看| 色综合站精品国产| 国产精品 欧美亚洲| 国产精品亚洲一级av第二区| 国产真人三级小视频在线观看| 性色av乱码一区二区三区2| 久久久久久亚洲精品国产蜜桃av| 亚洲午夜理论影院| 亚洲精品中文字幕在线视频| 免费在线观看完整版高清| 久久天躁狠狠躁夜夜2o2o| 成人精品一区二区免费| 在线天堂中文资源库| 人人澡人人妻人| 宅男免费午夜| 久久热在线av| 99精品欧美一区二区三区四区| 精品久久久久久久久久久久久 | 在线观看日韩欧美| 欧美一区二区精品小视频在线| 两性夫妻黄色片| 亚洲人成网站在线播放欧美日韩| 一边摸一边抽搐一进一小说| 侵犯人妻中文字幕一二三四区| 亚洲aⅴ乱码一区二区在线播放 | 老汉色∧v一级毛片| 一夜夜www| 国产伦一二天堂av在线观看| 免费av毛片视频| 黄色片一级片一级黄色片| 免费在线观看影片大全网站| 国产真实乱freesex| 亚洲一区二区三区色噜噜| 亚洲五月婷婷丁香| 国产视频一区二区在线看| 久久精品91无色码中文字幕| 久久久久国产精品人妻aⅴ院| 制服丝袜大香蕉在线| 亚洲av片天天在线观看| 黑丝袜美女国产一区| 久久久久久久久中文| 色综合站精品国产| 韩国av一区二区三区四区| 露出奶头的视频| 久久国产亚洲av麻豆专区| 19禁男女啪啪无遮挡网站| 人妻丰满熟妇av一区二区三区| 亚洲男人天堂网一区| 中文字幕最新亚洲高清| 99精品欧美一区二区三区四区| 两个人看的免费小视频| 天堂√8在线中文| 啦啦啦观看免费观看视频高清| 欧美黄色片欧美黄色片| 免费在线观看日本一区| 19禁男女啪啪无遮挡网站| 国产激情欧美一区二区| 色综合站精品国产| 欧美乱码精品一区二区三区| 男女之事视频高清在线观看| 美女国产高潮福利片在线看| a级毛片在线看网站| 婷婷精品国产亚洲av| 人成视频在线观看免费观看| 一级黄色大片毛片| 大香蕉久久成人网| 一区二区三区高清视频在线| 精品一区二区三区四区五区乱码| 亚洲成人精品中文字幕电影| 成人亚洲精品av一区二区| 一a级毛片在线观看| 亚洲男人的天堂狠狠| 日韩精品免费视频一区二区三区| 精品久久久久久久久久久久久 | 成人亚洲精品av一区二区| 午夜福利欧美成人| 国内久久婷婷六月综合欲色啪| 色综合亚洲欧美另类图片| 中文字幕最新亚洲高清| 亚洲av成人av| 亚洲中文日韩欧美视频| 国产精品精品国产色婷婷| 国产精品亚洲美女久久久| 国产亚洲精品久久久久5区| 成人永久免费在线观看视频| av视频在线观看入口| 亚洲色图av天堂| 18美女黄网站色大片免费观看| 精品不卡国产一区二区三区| 1024手机看黄色片| 久久久久亚洲av毛片大全| 亚洲免费av在线视频| 精品久久久久久,| 特大巨黑吊av在线直播 | 变态另类丝袜制服| 成人精品一区二区免费| 制服丝袜大香蕉在线| 欧美日韩瑟瑟在线播放| 18美女黄网站色大片免费观看| 中文亚洲av片在线观看爽| 制服诱惑二区| 久久精品91蜜桃| 99在线视频只有这里精品首页| 午夜福利高清视频| 亚洲一区高清亚洲精品| 久久精品aⅴ一区二区三区四区| 成人18禁在线播放| 亚洲精品美女久久久久99蜜臀| 国产精品电影一区二区三区| 亚洲性夜色夜夜综合| 久久精品亚洲精品国产色婷小说| 一个人观看的视频www高清免费观看 | 在线观看午夜福利视频| 欧美亚洲日本最大视频资源| 无人区码免费观看不卡| 波多野结衣av一区二区av| 欧美一级毛片孕妇| 1024香蕉在线观看| 日日爽夜夜爽网站| 免费高清视频大片| netflix在线观看网站| 搡老妇女老女人老熟妇| 欧美国产精品va在线观看不卡| 久久亚洲真实| 欧美亚洲日本最大视频资源| 在线看三级毛片| 18禁黄网站禁片免费观看直播| 精品日产1卡2卡| 国产蜜桃级精品一区二区三区| 国产精品一区二区免费欧美| 露出奶头的视频| 超碰成人久久| 99riav亚洲国产免费| 丝袜在线中文字幕| 最新在线观看一区二区三区| 久久精品国产亚洲av香蕉五月| 最新美女视频免费是黄的| 变态另类丝袜制服| 9191精品国产免费久久| 久久中文字幕人妻熟女| 亚洲精品一卡2卡三卡4卡5卡| 国产亚洲av嫩草精品影院| 久久久久久人人人人人| 久久精品国产99精品国产亚洲性色| 啦啦啦免费观看视频1| 国内久久婷婷六月综合欲色啪| 99国产极品粉嫩在线观看| 国产极品粉嫩免费观看在线| 精品久久久久久,| 亚洲国产精品久久男人天堂| 99精品欧美一区二区三区四区| 麻豆一二三区av精品| 欧美日韩乱码在线| 亚洲人成网站在线播放欧美日韩| 国产国语露脸激情在线看| 亚洲欧美日韩无卡精品| 天堂影院成人在线观看| av有码第一页| 岛国在线观看网站| tocl精华| 久久久国产欧美日韩av| 日韩欧美国产在线观看| 视频区欧美日本亚洲| 超碰成人久久| 亚洲五月婷婷丁香| 国产男靠女视频免费网站| 久久久久国产一级毛片高清牌| 成人三级黄色视频| 香蕉国产在线看| 很黄的视频免费| 欧美日韩黄片免| a在线观看视频网站| 手机成人av网站| 久久中文字幕一级| 可以在线观看毛片的网站| 久久天堂一区二区三区四区| 国产单亲对白刺激| av中文乱码字幕在线| 国产精品美女特级片免费视频播放器 | 男女午夜视频在线观看| x7x7x7水蜜桃| 99国产极品粉嫩在线观看| 亚洲人成伊人成综合网2020| 久久久久久九九精品二区国产 | 国产精品美女特级片免费视频播放器 | 18禁美女被吸乳视频| 熟妇人妻久久中文字幕3abv| 香蕉丝袜av| 日韩国内少妇激情av| 男女视频在线观看网站免费 | 成人国产一区最新在线观看| 一卡2卡三卡四卡精品乱码亚洲| 99精品欧美一区二区三区四区| 91在线观看av| www.www免费av| 神马国产精品三级电影在线观看 | 午夜激情av网站| 一级作爱视频免费观看| 搡老岳熟女国产| 亚洲天堂国产精品一区在线| 天天添夜夜摸| 亚洲人成网站在线播放欧美日韩| 日韩一卡2卡3卡4卡2021年| 久久婷婷人人爽人人干人人爱| 日日摸夜夜添夜夜添小说| 亚洲一码二码三码区别大吗| 12—13女人毛片做爰片一| 成年女人毛片免费观看观看9| 在线免费观看的www视频| 中文资源天堂在线| 久久精品成人免费网站| 人人妻,人人澡人人爽秒播| 日韩欧美在线二视频| 国产黄a三级三级三级人| 国产精品乱码一区二三区的特点| 91国产中文字幕| 色在线成人网| 日本一本二区三区精品| 欧美激情高清一区二区三区| 国产午夜福利久久久久久| 女人被狂操c到高潮| 国产欧美日韩一区二区三| 精品卡一卡二卡四卡免费| 色精品久久人妻99蜜桃| 久久久久亚洲av毛片大全| 精品久久久久久久久久久久久 | 久久青草综合色| 亚洲av第一区精品v没综合| 不卡一级毛片| 成熟少妇高潮喷水视频| 久久香蕉国产精品| av免费在线观看网站| 亚洲成人久久爱视频| 久久亚洲精品不卡| 中文字幕高清在线视频| 丝袜人妻中文字幕| 国产三级在线视频| 男女那种视频在线观看| 免费观看人在逋| 亚洲一区二区三区色噜噜| 99久久综合精品五月天人人| АⅤ资源中文在线天堂| 亚洲色图av天堂| 热99re8久久精品国产| 他把我摸到了高潮在线观看| 女人爽到高潮嗷嗷叫在线视频| 大香蕉久久成人网| 亚洲第一av免费看| 精品国内亚洲2022精品成人| 一个人观看的视频www高清免费观看 | 可以免费在线观看a视频的电影网站| 国产精品二区激情视频| www日本在线高清视频| av福利片在线| 亚洲精品久久成人aⅴ小说| 中文字幕高清在线视频| 18禁裸乳无遮挡免费网站照片 | 99久久国产精品久久久| 国产欧美日韩一区二区三| 国产伦在线观看视频一区| 禁无遮挡网站| www日本黄色视频网| 可以在线观看毛片的网站| 亚洲精华国产精华精| 国产男靠女视频免费网站| 国产在线精品亚洲第一网站| 日日干狠狠操夜夜爽| 18禁观看日本| 国产人伦9x9x在线观看| 国产伦人伦偷精品视频| 国产精品综合久久久久久久免费| svipshipincom国产片| av在线天堂中文字幕| 999精品在线视频| 亚洲国产中文字幕在线视频| 午夜福利免费观看在线| 亚洲人成网站在线播放欧美日韩| 老司机午夜福利在线观看视频| 国产黄a三级三级三级人| 亚洲欧美精品综合一区二区三区| 亚洲免费av在线视频| АⅤ资源中文在线天堂| 亚洲专区国产一区二区| 成年免费大片在线观看| 亚洲久久久国产精品| 婷婷亚洲欧美| 国产日本99.免费观看| 男女床上黄色一级片免费看| 白带黄色成豆腐渣| 免费在线观看亚洲国产| 亚洲久久久国产精品| 国产乱人伦免费视频| 日本撒尿小便嘘嘘汇集6| 波多野结衣高清作品| 亚洲国产精品合色在线| 免费在线观看亚洲国产| 天堂动漫精品| 亚洲专区国产一区二区| 一卡2卡三卡四卡精品乱码亚洲| 亚洲人成77777在线视频| 亚洲中文字幕一区二区三区有码在线看 | 亚洲av熟女| 亚洲精品在线观看二区| 精品国产乱子伦一区二区三区| 亚洲国产精品成人综合色| 神马国产精品三级电影在线观看 | 日本成人三级电影网站| 亚洲成a人片在线一区二区|