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

    基于EHL軸承支撐的曲軸多體系統(tǒng)動(dòng)力學(xué)仿真

    2015-03-09 03:33:04馬星國(guó)單楊智姝尤小梅
    機(jī)床與液壓 2015年21期
    關(guān)鍵詞:軸心油膜曲軸

    馬星國(guó),單楊智姝,尤小梅

    (沈陽(yáng)理工大學(xué)機(jī)械工程學(xué)院,遼寧沈陽(yáng) 110159)

    基于EHL軸承支撐的曲軸多體系統(tǒng)動(dòng)力學(xué)仿真

    馬星國(guó),單楊智姝,尤小梅

    (沈陽(yáng)理工大學(xué)機(jī)械工程學(xué)院,遼寧沈陽(yáng) 110159)

    為了更精確地分析無(wú)連桿發(fā)動(dòng)機(jī)曲軸系統(tǒng)動(dòng)力學(xué)特性,基于RecurDyn建立了無(wú)連桿發(fā)動(dòng)機(jī)曲軸軸系與EHL軸承耦合多體動(dòng)力學(xué)模型,進(jìn)行了發(fā)動(dòng)機(jī)曲軸系的動(dòng)力學(xué)特性分析與曲軸應(yīng)力分析,得到主軸承載荷、最大油膜壓力及最小油膜厚度等數(shù)據(jù)。與不計(jì)入油膜作用時(shí)的曲軸系多體動(dòng)力學(xué)仿真結(jié)果進(jìn)行比較,無(wú)連桿發(fā)動(dòng)機(jī)曲軸系與EHL動(dòng)力潤(rùn)滑耦合仿真分析結(jié)果得到的軸心軌跡更接近實(shí)際。

    曲軸系統(tǒng);多體動(dòng)力學(xué)仿真;無(wú)連桿發(fā)動(dòng)機(jī);動(dòng)力潤(rùn)滑

    0 前言

    目前,在進(jìn)行發(fā)動(dòng)機(jī)曲軸系統(tǒng)動(dòng)力學(xué)與油膜動(dòng)力潤(rùn)滑耦合分析時(shí),均是采用油膜軸承數(shù)學(xué)模型即直接求解Reynolds方程得到油膜壓力后利用數(shù)值積分求出軸承反力,不能很好地模擬復(fù)雜的軸系載荷變化。何芝仙等在建立曲軸軸承系統(tǒng)動(dòng)力學(xué)仿真模型時(shí)的基本思路就是去掉主軸承約束,用軸承反力替代。另外,曲軸軸承系統(tǒng)是一個(gè)不可分割的完整系統(tǒng),以往對(duì)于軸承的潤(rùn)滑分析并沒(méi)有考慮動(dòng)力學(xué)效應(yīng)對(duì)軸承潤(rùn)滑性能的影響,其實(shí)軸系動(dòng)力學(xué)與軸承摩擦學(xué)之間是相互影響的。所以要使計(jì)算結(jié)果符合工程實(shí)際,不僅要考慮滑動(dòng)軸承的各項(xiàng)參數(shù),還有軸系與軸承的耦合因素。

    本文作者基于RecurDyn多體動(dòng)力學(xué)軟件以及借助于RecurDyn中自帶的EHL流體潤(rùn)滑軸承子模塊建立的無(wú)連桿發(fā)動(dòng)機(jī)曲軸系統(tǒng)與EHL流體動(dòng)力潤(rùn)滑耦合仿真模型,進(jìn)行快速、可靠的軸承液體動(dòng)力學(xué)計(jì)算,完成了基于EHL流體潤(rùn)滑軸承支撐的無(wú)連桿發(fā)動(dòng)機(jī)曲軸系的動(dòng)力學(xué)仿真分析。文中所采用的方法可以更加真實(shí)地模擬曲軸系統(tǒng)的工況,更加準(zhǔn)確地分析曲軸的受力和變形。

    1 曲軸系統(tǒng)多體動(dòng)力與動(dòng)力潤(rùn)滑耦合仿真建模

    1.1 發(fā)動(dòng)機(jī)剛?cè)狁詈蟿?dòng)力學(xué)模型

    將在Solidworks中裝配后的無(wú)連桿發(fā)動(dòng)機(jī)曲軸系三維模型導(dǎo)入RecurDyn中,按實(shí)際運(yùn)動(dòng)情況確定各部件相互間的運(yùn)動(dòng)副。由于曲軸的彈性振動(dòng)變形對(duì)油膜承載及運(yùn)動(dòng)都會(huì)產(chǎn)生直接的影響,所以需要將曲軸建為柔性體。

    柔性體的運(yùn)動(dòng)方程是建立在廣義坐標(biāo)的基礎(chǔ)上的基于拉格朗日方程的控制微分方程的形式為

    曲軸剛體模型在RecurDyn中可直接通過(guò)Mesh Interface傳送至有限元前后處理器Femap中,將曲軸在Femap中劃分網(wǎng)格之后,通過(guò)Mesh Interface將Femap的網(wǎng)格傳送至RecurDyn中作為曲軸FFlex柔性體模型,無(wú)須再通過(guò)FFlex工具包中的Import按鈕導(dǎo)入在其他有限元軟件中處理的柔性體。中間省去數(shù)據(jù)的傳遞,減少了模型在不同軟件之間相互傳遞過(guò)程中數(shù)據(jù)的丟失。發(fā)動(dòng)機(jī)剛?cè)狁詈蟿?dòng)力學(xué)模型建立如圖1。

    圖1 發(fā)動(dòng)機(jī)剛?cè)狁詈隙囿w動(dòng)力學(xué)模型

    1.2 發(fā)動(dòng)機(jī)EHLBearing滑動(dòng)軸承單元的建立

    RecurDyn中帶有EHLBearing流體潤(rùn)滑軸承子模塊,可確定軸承負(fù)載,分析載荷分布并且通過(guò)對(duì)軸承進(jìn)行網(wǎng)格劃分,采用流體的方法分析油膜壓力變化。EHL軸承基于Reynolds方程。對(duì)潤(rùn)滑油作為不可壓縮牛頓流體、等黏度等假設(shè)情況下的Reynolds方程為:

    式中:Uw為軸頸表面線速度;Us為軸承表面線速度;η為潤(rùn)滑油動(dòng)力黏度;p為流體動(dòng)壓力;h為油膜厚度;t為時(shí)間。

    取軸承直徑為80 mm,軸承寬度為38 mm,軸瓦材料為錫青銅ZCuSn10Pb1,潤(rùn)滑間隙為0.02 mm,潤(rùn)滑油動(dòng)力黏度為0.01 Pa·s,油膜網(wǎng)格設(shè)置為周向30 mm,軸向7 mm,如圖2所示。定義完成后就可以建立發(fā)動(dòng)機(jī)多體動(dòng)力學(xué)與油膜動(dòng)力潤(rùn)滑耦合仿真模型。

    圖2 EHLBearing滑動(dòng)軸承油膜節(jié)點(diǎn)分布圖

    2 仿真及結(jié)果分析

    2.1 發(fā)動(dòng)機(jī)運(yùn)動(dòng)學(xué)及動(dòng)力學(xué)仿真結(jié)果

    在RecurDyn中采用Akami曲線擬合技術(shù)對(duì)示功圖的燃?xì)獗l(fā)壓力曲線進(jìn)行擬合,按照發(fā)火次序?qū)⑷細(xì)獗l(fā)壓力作用在活塞端部,施加轉(zhuǎn)速為6 000 r/min時(shí)各缸的燃?xì)獗l(fā)壓力如圖3所示。讓發(fā)動(dòng)機(jī)穩(wěn)定在6 000 r/min,對(duì)發(fā)動(dòng)機(jī)進(jìn)行仿真,仿真時(shí)間為一個(gè)周期0.02 s。

    圖3 燃?xì)獗l(fā)壓力圖 (按照發(fā)火順序)

    活塞軸向?yàn)閤軸,如圖4所示。

    圖4 曲軸系運(yùn)動(dòng)件質(zhì)心位移

    活塞和曲軸的質(zhì)心x向位移與動(dòng)平衡滑塊的質(zhì)心y向位移在一個(gè)工作周期內(nèi)都以正弦規(guī)律運(yùn)動(dòng),這個(gè)規(guī)律與理論相符合,說(shuō)明所建立的多剛體曲軸模型是正確的,可以模擬實(shí)際無(wú)連桿發(fā)動(dòng)機(jī)的運(yùn)動(dòng)規(guī)律。

    曲軸軸系模型進(jìn)行多體動(dòng)力學(xué)仿真,得到曲軸系在一個(gè)工作周期內(nèi)的仿真數(shù)據(jù)。圖5所示為活塞與動(dòng)平衡滑塊的加速度曲線圖,活塞與動(dòng)平衡滑塊的加速度變化規(guī)律基本相似但是方向相反,通過(guò)對(duì)曲軸軸系系統(tǒng)進(jìn)行動(dòng)平衡計(jì)算,結(jié)果發(fā)現(xiàn)活塞產(chǎn)生的不平衡慣性力與慣性力矩恰好與動(dòng)平衡滑塊產(chǎn)生慣性力與慣性力矩達(dá)到平衡。

    圖5 動(dòng)平衡滑塊與活塞加速度

    2.2 曲軸應(yīng)力分析

    通過(guò)RecurDyn中自帶的EHLBearing子模塊建立油膜模型,并對(duì)曲軸軸系進(jìn)行多體動(dòng)力學(xué)仿真。仿真得到曲軸柔性體的應(yīng)力分布云圖是模擬實(shí)際變化的油膜力施加在曲軸上。圖6為第一缸燃?xì)鈮毫Ρl(fā)時(shí)刻0.003 s曲軸應(yīng)力云圖,從圖中可以看出主軸頸處最大應(yīng)力值為10 MPa發(fā)生部位為過(guò)渡圓角處,圖7為0.003 2 s時(shí)刻曲軸應(yīng)力云圖,主軸頸處最大應(yīng)力值為22 MPa發(fā)生部位為過(guò)渡圓角處,兩圖對(duì)比看出曲軸受油膜壓力作用影響出現(xiàn)應(yīng)力值增大趨勢(shì)。這說(shuō)明燃?xì)獗l(fā)壓力影響油膜壓力的同時(shí)也影響主軸軸頸處的應(yīng)力。并且應(yīng)力集中在主軸頸過(guò)渡圓角處是由于它們局部尺寸比較小。圖8(a)所示為0.014 s時(shí)刻剛性支撐作用下曲軸應(yīng)力云圖,主軸頸處應(yīng)力值為56 MPa,圖8(b)所示為主軸頸過(guò)渡圓角處應(yīng)力值為10 MPa,對(duì)比兩圖看出油膜壓力對(duì)曲軸應(yīng)力影響比較平緩無(wú)沖擊。

    圖6 0.003 s時(shí)刻曲軸應(yīng)力云圖

    圖7 0.003 2 s時(shí)刻 曲軸應(yīng)力云圖

    圖8 0.014 s時(shí)刻曲軸應(yīng)力云圖

    如圖9所示曲軸的應(yīng)力最大點(diǎn)發(fā)生在0.014 s時(shí)刻,發(fā)生位置在曲軸與平衡圓滑塊接觸處如圖8(b)所示,最大應(yīng)力值為467.50 MPa,小于曲軸材料38CrMoAl的屈服強(qiáng)度690 MPa,曲軸強(qiáng)度滿足要求。

    圖9 曲軸各時(shí)間最大應(yīng)力圖

    曲軸發(fā)生最大位移時(shí)刻為0.014 s,發(fā)生部位在曲軸與平衡圓滑塊接觸處如圖10所示,最大位移值為0.031 312 mm,小于許用撓度0.133 mm,曲軸的剛度滿足要求。

    圖10 0.014 s時(shí)刻曲軸位移云圖

    2.3 EHL潤(rùn)滑動(dòng)力分析

    2.3.1 主軸承載荷分析

    曲軸輸入端的主軸承的負(fù)荷較大,故對(duì)輸入端主軸承的載荷進(jìn)行分析。全局坐標(biāo)系下曲軸軸向?yàn)閦軸,氣缸軸線方向?yàn)閤軸,圖11與圖12分別為主軸承x向載荷和y向載荷。從圖中可以看出,由于一個(gè)工作循環(huán)內(nèi)氣缸爆發(fā)壓力的間歇作用,主軸承x方向和y方向上載荷會(huì)相應(yīng)的表現(xiàn)出峰值響應(yīng)。由于曲軸系統(tǒng)設(shè)置恒定轉(zhuǎn)速,系統(tǒng)運(yùn)動(dòng)仿真開(kāi)始主軸承載荷x方向和y方向會(huì)出現(xiàn)沖擊。該主軸承與第二、三缸相鄰,則y方向載荷最大波峰出現(xiàn)在0.015 s左右,而第三缸的發(fā)火時(shí)間也在該時(shí)刻,此外在0.005 s左右y方向載荷出現(xiàn)沖擊由于燃?xì)獗l(fā)力開(kāi)始作用于第一缸發(fā)火在0.003 s左右。x方向載荷的最大波峰由于油膜阻尼作用,滯后出現(xiàn)在0.018 s左右。這種交變載荷作用在主軸承上使其受力變得復(fù)雜,并使機(jī)體產(chǎn)生振動(dòng)。

    圖11 主軸承x向載荷

    圖12 主軸承y向載荷

    2.3.2 主軸承最小油膜厚度分析

    圖13所示為主軸承最小油膜厚度隨時(shí)間變化關(guān)系圖,從圖中可以看出,主軸承在x、y向載荷最大峰值時(shí)刻油膜厚度達(dá)到最小峰值,說(shuō)明主軸承載荷變化對(duì)于軸承最小油膜厚度有著很大影響。油膜厚度最小峰值為1.3μm,大于奧地利李斯特研究所的推薦值 (主軸承 [hmin]:1.2~1.4μm),因而可以推斷出軸承工作可靠。

    圖13 主軸承最小油膜厚度

    2.3.3 主軸承最大油膜壓力分析

    圖14所示為主軸承周向油膜壓力隨曲軸轉(zhuǎn)角變化曲線。

    圖14 主軸承周向最大油膜壓力

    最大油膜壓力為43.7 MPa,出現(xiàn)在280°轉(zhuǎn)角處對(duì)應(yīng)于0.015 s左右,可以看出主軸承最大油膜壓力值出現(xiàn)的時(shí)間與最小油膜厚度最小值出現(xiàn)的時(shí)間是一致的,說(shuō)明主軸承受載荷作用及軸瓦變形對(duì)最大油膜壓力和最小油膜厚度有明顯影響,并且此時(shí)刻軸承工作條件最?lèi)毫?。最大油膜壓力峰值小?50 MPa的一般車(chē)用發(fā)動(dòng)機(jī)油膜壓力極限值,表明該軸承潤(rùn)滑程度良好。

    2.3.4 潤(rùn)滑油膜動(dòng)態(tài)平均壓力分布分析

    隨著載荷量的變化,軸變形導(dǎo)致的軸頸傾斜使油膜壓力產(chǎn)生偏布,變化規(guī)律如圖15所示。由于油膜厚度不足,最大油膜壓力位置向軸承端部移動(dòng)同時(shí)最大油膜壓力增大。在0.015 s時(shí)刻出現(xiàn)最大油膜壓力43.7 MPa,容易造成曲軸軸頸表面點(diǎn)蝕和軸承破壞。

    圖15 各時(shí)刻油膜壓力分布變化

    2.3.5 主軸承軸心軌跡對(duì)比分析

    圖16和圖17分別表示不計(jì)入油膜作用和采用EHLBearing軸承兩周期下的主軸承軸心軌跡。從圖上可以看出,兩種軸心軌跡所表現(xiàn)的軸心軌跡有明顯的差別。EHLBearing軸承在兩周期內(nèi)的軸心軌跡表現(xiàn)的比較一致,而不計(jì)入油膜作用的軸心軌跡在兩周期內(nèi)則無(wú)規(guī)律可遵循。這表明EHLBearing軸承可以模擬實(shí)際軸承模型并在動(dòng)力學(xué)仿真中考慮軸承的動(dòng)力潤(rùn)滑摩擦行為將更加真實(shí)的反應(yīng)實(shí)際工作條件下的發(fā)動(dòng)機(jī)曲軸系統(tǒng)。

    圖16 不計(jì)入油膜作用的軸心軌跡

    圖17 EHLBearing軸承軸心軌跡

    曲軸的振動(dòng)會(huì)對(duì)軸承上作用的載荷以及油膜厚度產(chǎn)生影響,因此通過(guò)分析軸心運(yùn)動(dòng)軌跡可以了解潤(rùn)滑油膜的分布情況及變化規(guī)律。由圖17可看出第二象限區(qū)為軸承承受最大負(fù)荷部位既第三缸燃?xì)鈮毫Ρl(fā)時(shí)刻,油膜壓力峰值劇增,軸承工作條件惡劣,要針對(duì)此處進(jìn)行對(duì)于軸承的最優(yōu)設(shè)計(jì),保證軸承處于液體潤(rùn)滑條件下工作;第四象限區(qū)負(fù)荷輕,軸承應(yīng)在此處開(kāi)油孔或油槽,使供油舒暢。

    3 結(jié)論

    通過(guò)建立無(wú)連桿發(fā)動(dòng)機(jī)曲軸軸系與EHL流體動(dòng)力潤(rùn)滑耦合多體動(dòng)力學(xué)模型,進(jìn)行了曲軸的應(yīng)力分析和油膜動(dòng)力潤(rùn)滑耦合分析,結(jié)果表明:

    (1)通過(guò)對(duì)活塞與動(dòng)平衡滑塊加速度仿真曲線的對(duì)比分析進(jìn)行動(dòng)平衡計(jì)算,結(jié)果表明發(fā)動(dòng)機(jī)整個(gè)運(yùn)動(dòng)過(guò)程基本上處于一個(gè)動(dòng)平衡狀態(tài)。

    (2)曲軸的最大應(yīng)力點(diǎn)在0.014 s曲軸與平衡圓滑塊的接觸點(diǎn),油膜壓力對(duì)曲軸應(yīng)力影響最大處為主軸頸過(guò)渡圓角處。

    (3)第三缸燃?xì)獗l(fā)時(shí)刻主軸承載荷沖擊最大,位于油膜厚度最薄處,是危險(xiǎn)部位,容易發(fā)生軸頸與軸瓦直接接觸現(xiàn)象。

    (4)與不計(jì)入油膜作用的軸心軌跡對(duì)比,EHLBearing軸承的軸心軌跡更接近實(shí)際系統(tǒng)。因此在做無(wú)連桿發(fā)動(dòng)機(jī)多體動(dòng)力學(xué)仿真分析的時(shí)候考慮流體動(dòng)力潤(rùn)滑耦合是非常必要的。

    [1]袁兆成.內(nèi)燃機(jī)設(shè)計(jì)[M].北京:機(jī)械工業(yè)出版社,2008.

    [2]林瓊,郝志勇,郭磊.曲軸系統(tǒng)多體動(dòng)力與油膜動(dòng)力潤(rùn)滑耦合的數(shù)字化仿真研究[J].內(nèi)燃機(jī)工程,2007,28(3):45-48.

    [3]程穎,宋瀟,孫善超.曲軸系柔性多體動(dòng)力學(xué)與動(dòng)力潤(rùn)滑耦合仿真[J].北京理工大學(xué)學(xué)報(bào),2006,26(4):314-317.

    [4]馬星國(guó),李想,尤小梅.曲軸系多柔體動(dòng)力學(xué)與動(dòng)力潤(rùn)滑耦合仿真分析[J].潤(rùn)滑與密封,2012,37(3):37 -41.

    [5]焦曉娟,張湝渭,彭曉斌.RecurDyn多體系統(tǒng)優(yōu)化仿真技術(shù)[M].北京:清華大學(xué)出版社,2010.

    [6]尹正兵,馬星國(guó),邵威.無(wú)連桿發(fā)動(dòng)機(jī)動(dòng)力學(xué)仿真研究[J].沈陽(yáng)理工大學(xué)學(xué)報(bào),2006,25(4):55 -58.

    [7]郝志勇,段秀兵,程金林.柴油機(jī)曲軸軸系的柔性多體動(dòng)力學(xué)仿真分析[J].鐵道機(jī)車(chē)車(chē)輛,2003(3):86-89.

    [8]何芝仙,桂長(zhǎng)林,李震,等.曲軸軸承系統(tǒng)動(dòng)力學(xué)和摩擦學(xué)行為耦合分析[J].車(chē)用發(fā)動(dòng)機(jī),2007(4):86-91.

    Crankshaft System Multibody Dynam ics Simulation Based on EHL Bearing Support

    MA Xingguo,SHAN Yangzhishu,YOU Xiaomei
    (School of Mechanical Engineering,Shenyang Ligong University,Shenyang Liaoning 110159,China)

    In order to analyze the dynamic characteristics of the no-link engine crankshaft system,amulti-body dynamics couplingmodel which consists of the no-link engine crankshaft system and EHL bearingwas builtbased on RecurDyn.In addition,the dynamic characteristics analysis of engine crankshaft system and crankshaft stress analysiswere analyzed,the EHL bearingwas studied to get themain bearing load,theminimum oil film thickness and themaximum oil film pressure.Afterwards,compared with the results of crankshaft system multi-body dynamics simulation without the oil film,it is concluded that the axis path ofmulti-body dynamics couplingmodel,which consists of the no-link engine crankshaft system and EHL hydrodynamic lubrication has been closer to the actual system.

    Crankshaft system;Multi-body dynamics simulation;No-link engine;EHL hydrodynamic lubrication

    TK422

    A

    1001-3881(2015)21-188-5

    10.3969/j.issn.1001 -3881.2015.21.047

    2014-08-20

    馬星國(guó) (1963—),男,工學(xué)博士,教授,主要研究方向?yàn)闄C(jī)械設(shè)計(jì)及理論、車(chē)輛動(dòng)力學(xué)及仿真。

    單楊智姝,E-mail:313344578@qq.com。

    猜你喜歡
    軸心油膜曲軸
    長(zhǎng)城油膜軸承油在高速棒材生產(chǎn)線的應(yīng)用
    石油商技(2021年1期)2021-03-29 02:36:04
    淺析曲軸磨床的可靠性改進(jìn)
    鋼結(jié)構(gòu)軸心受壓構(gòu)件穩(wěn)定性分析
    一種曲軸豎置熱處理裝置設(shè)計(jì)
    基于支持向量機(jī)的曲軸動(dòng)平衡控制方法
    大型數(shù)控立式磨床靜壓轉(zhuǎn)臺(tái)油膜熱特性仿真及其實(shí)驗(yàn)分析
    CFRP和角鋼復(fù)合加固混凝土矩形柱軸心受壓承載力
    以門(mén)靜脈-腸系膜上靜脈為軸心的腹腔鏡胰十二指腸切除術(shù)16例報(bào)道
    冷軋軋機(jī)油膜軸承系統(tǒng)故障分析與對(duì)策
    新疆鋼鐵(2016年4期)2016-05-25 00:37:38
    以文化為軸心,睿智教文言
    在线观看免费高清a一片| 国产成人精品久久久久久| 精品免费久久久久久久清纯 | 一区二区三区四区激情视频| 欧美97在线视频| 狂野欧美激情性xxxx| 欧美 日韩 精品 国产| 国产91精品成人一区二区三区 | 国产一级毛片在线| 欧美日韩成人在线一区二区| 自线自在国产av| 操美女的视频在线观看| 在线观看一区二区三区激情| 一区二区av电影网| 50天的宝宝边吃奶边哭怎么回事| 男女午夜视频在线观看| 成人影院久久| 99re6热这里在线精品视频| 50天的宝宝边吃奶边哭怎么回事| 大香蕉久久网| 男人爽女人下面视频在线观看| 1024视频免费在线观看| 中文乱码字字幕精品一区二区三区| 日韩熟女老妇一区二区性免费视频| 欧美精品一区二区大全| 男人爽女人下面视频在线观看| 美女大奶头黄色视频| 国产精品一二三区在线看| 中文字幕色久视频| 久久国产精品男人的天堂亚洲| 欧美人与性动交α欧美精品济南到| 亚洲国产最新在线播放| 国产亚洲av高清不卡| 国产成人精品久久久久久| 国产精品99久久99久久久不卡| 一区在线观看完整版| 涩涩av久久男人的天堂| av在线播放精品| 久久久国产精品麻豆| 久久久久久久精品精品| 精品久久久久久久毛片微露脸 | 女警被强在线播放| 日日夜夜操网爽| 亚洲av日韩精品久久久久久密 | 亚洲国产精品一区三区| 香蕉丝袜av| 亚洲欧美一区二区三区国产| 人人妻人人澡人人看| 国产真人三级小视频在线观看| 日本欧美国产在线视频| 国产黄频视频在线观看| 国产高清视频在线播放一区 | 亚洲久久久国产精品| 又紧又爽又黄一区二区| 青青草视频在线视频观看| 国产精品欧美亚洲77777| a级片在线免费高清观看视频| 亚洲国产最新在线播放| 夫妻性生交免费视频一级片| 国产成人欧美| 欧美日本中文国产一区发布| 免费在线观看日本一区| 免费高清在线观看日韩| 久久久久久久久免费视频了| 一区二区三区激情视频| 狂野欧美激情性xxxx| 国产伦理片在线播放av一区| 久久国产亚洲av麻豆专区| 欧美日韩福利视频一区二区| av在线老鸭窝| 精品欧美一区二区三区在线| 嫩草影视91久久| 青草久久国产| 晚上一个人看的免费电影| 欧美国产精品一级二级三级| 婷婷色av中文字幕| 少妇裸体淫交视频免费看高清 | 伊人久久大香线蕉亚洲五| 日韩熟女老妇一区二区性免费视频| 精品一区在线观看国产| 熟女av电影| 99国产精品99久久久久| 亚洲七黄色美女视频| 欧美日韩精品网址| 18禁裸乳无遮挡动漫免费视频| 自线自在国产av| 色综合欧美亚洲国产小说| 男人添女人高潮全过程视频| av一本久久久久| 欧美成人午夜精品| 午夜福利免费观看在线| 老司机靠b影院| 一区二区三区四区激情视频| 最黄视频免费看| 美女中出高潮动态图| 日日夜夜操网爽| 99热网站在线观看| 亚洲欧美精品自产自拍| 国产一卡二卡三卡精品| 黄网站色视频无遮挡免费观看| 欧美日韩福利视频一区二区| 国产午夜精品一二区理论片| 纯流量卡能插随身wifi吗| 欧美日韩一级在线毛片| 丁香六月天网| 欧美精品av麻豆av| 男女免费视频国产| 黄网站色视频无遮挡免费观看| 日韩一卡2卡3卡4卡2021年| 十八禁高潮呻吟视频| 亚洲色图综合在线观看| 国产色视频综合| 午夜福利一区二区在线看| 真人做人爱边吃奶动态| 国产成人精品久久久久久| 精品人妻一区二区三区麻豆| 精品卡一卡二卡四卡免费| 婷婷色综合大香蕉| 久久久精品免费免费高清| 免费在线观看影片大全网站 | 亚洲精品久久久久久婷婷小说| 亚洲国产毛片av蜜桃av| 欧美日韩av久久| 少妇人妻 视频| 欧美日本中文国产一区发布| 丝瓜视频免费看黄片| 夜夜骑夜夜射夜夜干| 99热网站在线观看| 亚洲av欧美aⅴ国产| 另类亚洲欧美激情| 一本久久精品| 国产国语露脸激情在线看| 看免费成人av毛片| 欧美人与善性xxx| h视频一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 国产精品一二三区在线看| 欧美日韩黄片免| 蜜桃国产av成人99| 亚洲精品第二区| 男人爽女人下面视频在线观看| 亚洲 欧美一区二区三区| 日韩制服骚丝袜av| 日本a在线网址| 国产高清不卡午夜福利| 赤兔流量卡办理| 新久久久久国产一级毛片| 少妇 在线观看| 日韩一卡2卡3卡4卡2021年| 亚洲,欧美,日韩| 国产精品国产三级国产专区5o| 久久ye,这里只有精品| 天天躁夜夜躁狠狠久久av| 国产一卡二卡三卡精品| 日本五十路高清| 国产黄色视频一区二区在线观看| 亚洲欧洲精品一区二区精品久久久| 女性被躁到高潮视频| 中文字幕人妻丝袜制服| 亚洲欧美色中文字幕在线| 日韩中文字幕欧美一区二区 | 久久中文字幕一级| 晚上一个人看的免费电影| 老司机午夜十八禁免费视频| 亚洲中文av在线| av天堂在线播放| 午夜激情av网站| 黄色毛片三级朝国网站| 日韩一卡2卡3卡4卡2021年| 19禁男女啪啪无遮挡网站| 女人高潮潮喷娇喘18禁视频| 一边摸一边做爽爽视频免费| 国产成人a∨麻豆精品| 日韩av免费高清视频| 男人操女人黄网站| 纵有疾风起免费观看全集完整版| 欧美激情极品国产一区二区三区| 老熟女久久久| 亚洲av片天天在线观看| 欧美日韩亚洲国产一区二区在线观看 | 亚洲熟女毛片儿| 校园人妻丝袜中文字幕| 青草久久国产| 久久人妻熟女aⅴ| av视频免费观看在线观看| 国产爽快片一区二区三区| 热99国产精品久久久久久7| av国产精品久久久久影院| av福利片在线| 免费看不卡的av| 午夜视频精品福利| 久久久久久久精品精品| 亚洲精品久久午夜乱码| 少妇裸体淫交视频免费看高清 | 在线精品无人区一区二区三| 丝袜在线中文字幕| 五月天丁香电影| 啦啦啦中文免费视频观看日本| 人人澡人人妻人| 日韩伦理黄色片| 久久人妻熟女aⅴ| 国产欧美日韩精品亚洲av| 国产精品久久久久久精品电影小说| 高潮久久久久久久久久久不卡| 伦理电影免费视频| 少妇人妻 视频| 999精品在线视频| 国产免费一区二区三区四区乱码| av网站在线播放免费| 色婷婷av一区二区三区视频| 中文字幕人妻丝袜制服| 热99久久久久精品小说推荐| 侵犯人妻中文字幕一二三四区| 免费一级毛片在线播放高清视频 | videos熟女内射| 免费一级毛片在线播放高清视频 | 欧美在线一区亚洲| 精品视频人人做人人爽| 丁香六月天网| 国产精品国产三级国产专区5o| 亚洲专区国产一区二区| 欧美日韩av久久| 国产视频一区二区在线看| 精品少妇一区二区三区视频日本电影| 免费看av在线观看网站| 丝瓜视频免费看黄片| 好男人电影高清在线观看| 欧美97在线视频| 国产一区二区三区综合在线观看| 夫妻性生交免费视频一级片| 免费女性裸体啪啪无遮挡网站| 日韩av不卡免费在线播放| 69精品国产乱码久久久| 午夜老司机福利片| 一边亲一边摸免费视频| 欧美日韩视频高清一区二区三区二| 美女福利国产在线| 男男h啪啪无遮挡| 成年人午夜在线观看视频| 精品视频人人做人人爽| videosex国产| 国产一区二区三区av在线| 国产精品一国产av| 黄色片一级片一级黄色片| 一级片'在线观看视频| 国产女主播在线喷水免费视频网站| 啦啦啦 在线观看视频| 精品免费久久久久久久清纯 | 国产精品av久久久久免费| 久久免费观看电影| 亚洲熟女精品中文字幕| 别揉我奶头~嗯~啊~动态视频 | 亚洲五月婷婷丁香| 麻豆国产av国片精品| 各种免费的搞黄视频| 黑人猛操日本美女一级片| 9色porny在线观看| 老司机深夜福利视频在线观看 | 97人妻天天添夜夜摸| 亚洲av成人不卡在线观看播放网 | 97在线人人人人妻| 久久久精品免费免费高清| 99国产精品免费福利视频| 最近最新中文字幕大全免费视频 | 亚洲精品第二区| 久久精品国产a三级三级三级| 国产午夜精品一二区理论片| 大香蕉久久网| 日韩一本色道免费dvd| 久久久精品国产亚洲av高清涩受| 搡老乐熟女国产| 91九色精品人成在线观看| 男人爽女人下面视频在线观看| 亚洲精品中文字幕在线视频| 新久久久久国产一级毛片| 欧美老熟妇乱子伦牲交| 欧美激情高清一区二区三区| 亚洲人成电影观看| 亚洲免费av在线视频| 在线观看国产h片| 我要看黄色一级片免费的| 晚上一个人看的免费电影| 午夜影院在线不卡| 看免费成人av毛片| 丝袜人妻中文字幕| 18禁黄网站禁片午夜丰满| 亚洲综合色网址| 国产人伦9x9x在线观看| 国产在视频线精品| 侵犯人妻中文字幕一二三四区| 涩涩av久久男人的天堂| 国产在线视频一区二区| 一级片免费观看大全| 日韩制服骚丝袜av| 欧美中文综合在线视频| 91精品国产国语对白视频| 国产日韩欧美视频二区| av线在线观看网站| 人人妻人人爽人人添夜夜欢视频| 精品福利永久在线观看| 国产精品秋霞免费鲁丝片| 国产一区二区 视频在线| 亚洲精品国产av成人精品| 久久人人97超碰香蕉20202| 亚洲成av片中文字幕在线观看| 亚洲av国产av综合av卡| 91麻豆av在线| 欧美 亚洲 国产 日韩一| 午夜免费成人在线视频| 日韩,欧美,国产一区二区三区| 在线观看免费日韩欧美大片| 多毛熟女@视频| 成年人午夜在线观看视频| 黑人巨大精品欧美一区二区蜜桃| 深夜精品福利| 大话2 男鬼变身卡| 成人国产av品久久久| 啦啦啦啦在线视频资源| 一区福利在线观看| 国产成人精品久久久久久| 十分钟在线观看高清视频www| 国产av国产精品国产| 国产无遮挡羞羞视频在线观看| 午夜免费观看性视频| 亚洲国产精品一区二区三区在线| 99热国产这里只有精品6| xxxhd国产人妻xxx| 婷婷成人精品国产| 18禁国产床啪视频网站| 免费女性裸体啪啪无遮挡网站| 国产成人91sexporn| av国产精品久久久久影院| 久久久久久免费高清国产稀缺| 欧美黑人欧美精品刺激| 国产一区二区三区av在线| 中文字幕最新亚洲高清| 久久国产精品影院| 国产午夜精品一二区理论片| 大型av网站在线播放| 久久久精品国产亚洲av高清涩受| 亚洲国产欧美网| 亚洲专区中文字幕在线| kizo精华| 大话2 男鬼变身卡| 2018国产大陆天天弄谢| 男人爽女人下面视频在线观看| 91成人精品电影| 国产精品.久久久| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲色图综合在线观看| 免费日韩欧美在线观看| 国产在线一区二区三区精| 人人妻人人澡人人爽人人夜夜| 亚洲国产av影院在线观看| 亚洲欧洲精品一区二区精品久久久| a级毛片黄视频| 18禁国产床啪视频网站| 国产有黄有色有爽视频| 久久久精品国产亚洲av高清涩受| 一边亲一边摸免费视频| 午夜日韩欧美国产| 免费在线观看影片大全网站 | 一个人免费看片子| 国产精品.久久久| 欧美日韩视频精品一区| 精品第一国产精品| 老汉色∧v一级毛片| 午夜福利,免费看| 波多野结衣一区麻豆| 人人澡人人妻人| 麻豆乱淫一区二区| 国产成人一区二区在线| 如日韩欧美国产精品一区二区三区| 在线观看国产h片| 韩国高清视频一区二区三区| 欧美黑人精品巨大| 国产女主播在线喷水免费视频网站| 最新在线观看一区二区三区 | 欧美国产精品va在线观看不卡| 久久久国产欧美日韩av| 成人亚洲欧美一区二区av| 一级毛片 在线播放| 王馨瑶露胸无遮挡在线观看| 狠狠精品人妻久久久久久综合| 亚洲成av片中文字幕在线观看| a级毛片在线看网站| 日韩av免费高清视频| 王馨瑶露胸无遮挡在线观看| 蜜桃国产av成人99| www.av在线官网国产| 国产成人一区二区三区免费视频网站 | 免费看十八禁软件| 国产男人的电影天堂91| 老司机亚洲免费影院| 国产又色又爽无遮挡免| 如日韩欧美国产精品一区二区三区| 91老司机精品| 国产欧美日韩一区二区三区在线| 亚洲av电影在线进入| 一区二区av电影网| 99国产精品一区二区三区| 国产主播在线观看一区二区 | 久久久久久久精品精品| 美女视频免费永久观看网站| 人人妻人人爽人人添夜夜欢视频| 免费人妻精品一区二区三区视频| 青青草视频在线视频观看| 妹子高潮喷水视频| 国产一区二区激情短视频 | 免费高清在线观看视频在线观看| 青青草视频在线视频观看| 国产欧美日韩一区二区三区在线| www.av在线官网国产| 国产成人欧美在线观看 | 国产亚洲欧美精品永久| 2021少妇久久久久久久久久久| 晚上一个人看的免费电影| 国产亚洲欧美在线一区二区| 欧美黄色淫秽网站| 男人爽女人下面视频在线观看| 国产午夜精品一二区理论片| 这个男人来自地球电影免费观看| 久久久精品区二区三区| 老司机亚洲免费影院| 久久亚洲国产成人精品v| 国产欧美日韩综合在线一区二区| 亚洲欧美清纯卡通| 午夜福利乱码中文字幕| 青草久久国产| 精品国产国语对白av| 欧美变态另类bdsm刘玥| 国产欧美日韩综合在线一区二区| 高清欧美精品videossex| 18禁国产床啪视频网站| 国产精品秋霞免费鲁丝片| 亚洲午夜精品一区,二区,三区| 亚洲,欧美精品.| 国产成人av教育| 国产老妇伦熟女老妇高清| 91麻豆精品激情在线观看国产 | 欧美成人精品欧美一级黄| 免费高清在线观看日韩| 多毛熟女@视频| 91字幕亚洲| 国产免费福利视频在线观看| 成人影院久久| 国产精品久久久久久精品电影小说| 欧美黑人欧美精品刺激| 国产精品亚洲av一区麻豆| 水蜜桃什么品种好| 九草在线视频观看| 国产熟女午夜一区二区三区| 日韩免费高清中文字幕av| 亚洲,欧美,日韩| 国产精品秋霞免费鲁丝片| av国产精品久久久久影院| 欧美在线黄色| 日本av免费视频播放| 最近手机中文字幕大全| 欧美大码av| 女人被躁到高潮嗷嗷叫费观| 高清黄色对白视频在线免费看| 天堂8中文在线网| 99久久人妻综合| 交换朋友夫妻互换小说| 91精品伊人久久大香线蕉| 80岁老熟妇乱子伦牲交| 亚洲天堂av无毛| 亚洲精品日本国产第一区| 捣出白浆h1v1| 国产激情久久老熟女| 欧美 日韩 精品 国产| 精品免费久久久久久久清纯 | 亚洲国产欧美在线一区| 亚洲熟女毛片儿| 久久国产精品大桥未久av| 伦理电影免费视频| 乱人伦中国视频| a级毛片在线看网站| 中文字幕高清在线视频| 另类亚洲欧美激情| 久久这里只有精品19| 日韩 亚洲 欧美在线| 成在线人永久免费视频| 欧美大码av| 在线观看人妻少妇| 永久免费av网站大全| 国产97色在线日韩免费| 秋霞在线观看毛片| 欧美老熟妇乱子伦牲交| 免费高清在线观看视频在线观看| 最近手机中文字幕大全| 黄片小视频在线播放| 亚洲成国产人片在线观看| 亚洲精品在线美女| 大型av网站在线播放| 99国产精品99久久久久| 欧美黑人精品巨大| 两人在一起打扑克的视频| 国产精品一二三区在线看| 美女大奶头黄色视频| 国语对白做爰xxxⅹ性视频网站| 久久天堂一区二区三区四区| 成人国语在线视频| 在线观看www视频免费| 大话2 男鬼变身卡| 日日爽夜夜爽网站| 午夜福利一区二区在线看| e午夜精品久久久久久久| 免费看av在线观看网站| 亚洲成色77777| 日韩av免费高清视频| 天堂俺去俺来也www色官网| 国产精品.久久久| 搡老乐熟女国产| 免费黄频网站在线观看国产| 久久性视频一级片| 好男人视频免费观看在线| 亚洲av美国av| 日本猛色少妇xxxxx猛交久久| 亚洲伊人色综图| 免费看十八禁软件| 亚洲欧美清纯卡通| 女人久久www免费人成看片| 亚洲成国产人片在线观看| 亚洲 国产 在线| 欧美成狂野欧美在线观看| 一级a爱视频在线免费观看| 亚洲九九香蕉| 亚洲成人手机| 国产精品二区激情视频| 国产精品欧美亚洲77777| 少妇精品久久久久久久| 国产欧美日韩综合在线一区二区| 韩国高清视频一区二区三区| 美女大奶头黄色视频| 久久精品亚洲熟妇少妇任你| 老司机影院毛片| 午夜福利视频精品| 成年人免费黄色播放视频| 国产一区二区 视频在线| 午夜免费男女啪啪视频观看| 日韩av免费高清视频| 亚洲国产毛片av蜜桃av| 亚洲欧美成人综合另类久久久| 婷婷成人精品国产| 欧美日韩成人在线一区二区| 黄网站色视频无遮挡免费观看| 亚洲,欧美,日韩| 人人澡人人妻人| 亚洲久久久国产精品| 色视频在线一区二区三区| 一本一本久久a久久精品综合妖精| 国产精品一区二区在线不卡| 国产av一区二区精品久久| 久久久国产一区二区| 美女午夜性视频免费| 可以免费在线观看a视频的电影网站| 99久久人妻综合| 亚洲av成人精品一二三区| 天天躁夜夜躁狠狠躁躁| 最黄视频免费看| 久久久久久久国产电影| 69精品国产乱码久久久| 欧美成人精品欧美一级黄| 成年人黄色毛片网站| 亚洲精品在线美女| 国产一区二区三区av在线| 久久久久久免费高清国产稀缺| 国产色视频综合| 另类精品久久| 午夜av观看不卡| 免费在线观看影片大全网站 | 伦理电影免费视频| 热99久久久久精品小说推荐| 久久久精品94久久精品| 精品视频人人做人人爽| 国产精品三级大全| 一区福利在线观看| 美女扒开内裤让男人捅视频| 日韩精品免费视频一区二区三区| 久久久国产一区二区| 午夜免费观看性视频| √禁漫天堂资源中文www| 国产午夜精品一二区理论片| 亚洲精品国产色婷婷电影| 久久天堂一区二区三区四区| av天堂久久9| 91精品三级在线观看| 国产成人av激情在线播放| 亚洲,欧美精品.| 黄色怎么调成土黄色| 国产三级黄色录像| 9191精品国产免费久久| 天天添夜夜摸| 免费在线观看视频国产中文字幕亚洲 | 久久精品久久久久久久性| 日本av免费视频播放| 久久久久久人人人人人| 操出白浆在线播放| 精品亚洲乱码少妇综合久久| 欧美黄色淫秽网站| 亚洲国产看品久久| 亚洲精品国产色婷婷电影| 久久av网站| 中文字幕av电影在线播放| 亚洲成色77777| 2021少妇久久久久久久久久久| 国产成人精品无人区| 97精品久久久久久久久久精品| 一区二区三区精品91| 欧美日韩亚洲高清精品| 女警被强在线播放| 精品人妻一区二区三区麻豆| 国产97色在线日韩免费|