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

    聯(lián)合載荷下角接觸球軸承靜態(tài)剛度仿真研究

    2022-09-23 11:49:48王建祥俞成濤
    關(guān)鍵詞:理論值滾珠內(nèi)圈

    王建祥,俞成濤,蔣 鵬,陳 宇

    (1.江蘇理工學(xué)院 機(jī)械工程學(xué)院,江蘇 常州 213000;2常州技師學(xué)院 智能制造學(xué)院,江蘇 常州 213032)

    憑借在不同工況下能夠穩(wěn)定地以較高轉(zhuǎn)速工作的顯著優(yōu)點(diǎn),角接觸球軸承成為各類轉(zhuǎn)子-軸承系統(tǒng)中的關(guān)鍵支承件之一。在實(shí)際裝備設(shè)計(jì)中,由試驗(yàn)來(lái)得出軸承的性能參數(shù)往往需要付出高昂的代價(jià);因此,隨著技術(shù)的發(fā)展,對(duì)軸承進(jìn)行仿真研究已成為工程應(yīng)用中的通用方法。

    角接觸球軸承的靜態(tài)特性對(duì)轉(zhuǎn)子系統(tǒng)的性能起著關(guān)鍵作用,而靜態(tài)剛度又最能突出軸承的靜態(tài)特性;因此,對(duì)這方面的研究應(yīng)加以重視。目前,國(guó)內(nèi)外學(xué)者對(duì)角接觸球軸承靜態(tài)特性理論的研究基本集中在以下兩個(gè)方面。

    其一,是通過(guò)赫茲接觸理論建立軸承的滾珠受力接觸模型。如:Jones[1]提出了溝道控制理論并建立了三自由度計(jì)算模型,分析了聯(lián)合載荷下球軸承的載荷分布、變形以及應(yīng)力之間的關(guān)系;陳宗農(nóng)等人[2]基于Jones的溝道控制理論建立了聯(lián)合載荷下球軸承的五自由度剛度求解矩陣;張迅雷等人[3]提出一種角接觸球軸承靜態(tài)剛度的精確計(jì)算公式,并通過(guò)對(duì)比簡(jiǎn)化計(jì)算值和試驗(yàn)值,驗(yàn)證了計(jì)算公式的準(zhǔn)確性;劉艷華等人[4]針對(duì)定位與定壓兩種不同預(yù)緊方式,分析計(jì)算了聯(lián)合載荷下角接觸球軸承的靜態(tài)剛度;方兵等人[5]通過(guò)坐標(biāo)變換法分析軸承滾珠受力與變形的關(guān)系,并對(duì)剛度進(jìn)行了求解與試驗(yàn)驗(yàn)證;Liu等人[6]提出了一種在預(yù)加載下提高軸承載荷分布及剛度的精度和計(jì)算效率的方法;師浩浩等人[7]借助最小二乘法和線性回歸方法,對(duì)赫茲接觸理論進(jìn)行計(jì)算方法簡(jiǎn)化,并分析了軸承的靜態(tài)剛度;李震等人[8]基于Demul理論和赫茲接觸理論建立了五自由度的靜態(tài)分析模型,研究了外載荷對(duì)接觸角、接觸力及靜態(tài)剛度的影響。

    其二,是利用有限元分析法對(duì)球軸承的幾何模型進(jìn)行約束加載分析。如:Tang等人[9]利用ANSYS的APDL參數(shù)化命令流對(duì)球軸承進(jìn)行有限元建模,將加載后的接觸特性結(jié)果與理論結(jié)果進(jìn)行對(duì)比,并證明了方法的實(shí)用性;Guo等人[10]采用接觸力學(xué)和有限元相結(jié)合的方法求解兩種滾動(dòng)軸承的剛度矩陣,并提出了一種數(shù)值方法來(lái)確定軸承的剛度。張永強(qiáng)[11]使用了球-板接觸的有限元模型,研究了網(wǎng)格單元密度對(duì)接觸應(yīng)力結(jié)果的影響,并通過(guò)對(duì)比軸向變形的有限元值與理論值,證明了方法的可行性。

    然而,以上研究中角接觸球軸承有限元法的分析,對(duì)變形和剛度的變化規(guī)律并未與理論計(jì)算作詳細(xì)的對(duì)比;因此,本文以赫茲接觸理論為理論計(jì)算模型的基礎(chǔ),應(yīng)用有限元仿真方法研究聯(lián)合載荷對(duì)角接觸球軸承變形和剛度的影響,通過(guò)對(duì)比驗(yàn)證,為實(shí)際裝備設(shè)計(jì)中軸承的靜態(tài)剛度仿真分析提供一定的理論參考。

    1 角接觸球軸承理論計(jì)算模型

    1.1 滾珠受力變形分析

    如圖1所示,為角接觸球軸承第j個(gè)滾珠的受力接觸示意圖。設(shè)置坐標(biāo)OXr、OX與Or分別為軸向載荷與徑向載荷的施加方向;Oi與Oo分別為內(nèi)外圈滾道的溝曲率中心位置;a0為受載荷前滾珠與內(nèi)外圈的初始接觸角。

    圖1 滾珠受力接觸示意圖

    受載荷前內(nèi)外圈溝曲率中心Oi與Oo間的初始距離A0為:

    式中:ri、ro、Dw分別為內(nèi)、外圈溝曲率半徑和滾珠直徑;fo、fi分別為內(nèi)、外圈的溝道比。

    固定軸承外圈,對(duì)內(nèi)圈同時(shí)施加徑向載荷Fr和軸向載荷Fa,滾珠會(huì)受到由內(nèi)、外圈溝滾道對(duì)其的法向作用力?;诤掌澖佑|理論[12],第j個(gè)滾珠受到的接觸力Qj與其變形δj的關(guān)系式為:

    式中,Kn為內(nèi)外圈與滾珠間的總接觸剛度,由以下關(guān)系式求得:

    其中:Ki為內(nèi)圈與滾珠間的法向接觸剛度;Ko為外圈與滾珠間的法向接觸剛度。

    如圖2所示,為聯(lián)合載荷下軸承溝曲率中心軌跡。受載后,軸承外滾道曲率中心軌跡不變,而內(nèi)滾道發(fā)生變化并產(chǎn)生軸向位移δa、徑向位移δr和角位移θ,且滾珠與內(nèi)外圈的接觸角由a0變?yōu)閍 j。在滾珠位置角φ j處的內(nèi)外圈溝曲率中心距離由A0變?yōu)锳j,關(guān)系式為:

    圖2 聯(lián)合載荷下溝曲率中心軌跡

    式(4)中滾珠位置角φj可由下式確定:

    內(nèi)圈溝曲率中心軌跡半徑Ri表達(dá)式為:

    第j個(gè)滾珠的接觸變形總量δj為:

    受載后第j個(gè)滾珠實(shí)際接觸角αj為:

    1.2 軸承的整體載荷-位移關(guān)系及剛度求解

    在具有N個(gè)滾珠的角接觸球軸承內(nèi),可列出軸向、徑向以及力矩的平衡方程[13]為:

    式(9)是以位移δa、δr和θ為未知量的非線性方程組。編寫Matlab程序并設(shè)置加載的已知聯(lián)合載荷進(jìn)行理論計(jì)算,方程求解采用Newton-Raphson迭代法。對(duì)軸承位移未知量進(jìn)行迭代計(jì)算并對(duì)位移求導(dǎo),其剛度矩陣解析表達(dá)式為:

    2 角接觸球軸承有限元仿真模型

    2.1 軸承模型三維實(shí)體設(shè)置

    選用7003C角接觸球軸承作為仿真對(duì)象并建立有限元分析模型,材料為軸承鋼,彈性模量為2.07×105MPa,泊松比為0.3。如表1所示,為軸承具體的尺寸參數(shù)。

    表1 軸承基本參數(shù)

    為簡(jiǎn)化計(jì)算,不考慮保持架的建模,在SolidWorks中嚴(yán)格按照表中軸承實(shí)際尺寸建立三維模型,再導(dǎo)入ANSYS Workbench的靜力學(xué)模塊中,并設(shè)置軸承的相關(guān)材料屬性。首先,對(duì)內(nèi)外圈采用掃掠網(wǎng)格劃分的方法(單元尺寸為0.4 mm),滾珠設(shè)置為自動(dòng)四面體網(wǎng)格劃分(單元尺寸為0.1 mm),以保證滾珠接觸部分的網(wǎng)格細(xì)化;然后,對(duì)滾珠與內(nèi)外圈滾道的接觸部分進(jìn)行局部細(xì)化(細(xì)化網(wǎng)格單元尺寸為0.1 mm),局部細(xì)化可采用Contact Sizing、Face Sizing和Inflation等方法。

    如圖3所示,為模型的網(wǎng)格分布。網(wǎng)格單元類型默認(rèn)為solid186,共節(jié)點(diǎn)2 227 060個(gè),單元1 146 185個(gè)。

    圖3 7003C角接觸球軸承有限元網(wǎng)格劃分

    2.2 模型條件設(shè)定

    2.2.1 接觸設(shè)置

    (1)依據(jù)凸面或者易變形面為接觸面的準(zhǔn)則[14]來(lái)設(shè)置接觸對(duì),滾珠為接觸面,內(nèi)外圈滾道為目標(biāo)面。

    (2)接觸類型為摩擦接觸,設(shè)置接觸摩擦系數(shù)為0.1,打開Small Sliding并設(shè)置穿透容差系數(shù)為0.1,以提高收斂速度。

    (3)接觸算法可以采用純罰函數(shù)法或增廣拉格朗日法,這兩種方法的計(jì)算原理與理論計(jì)算中滾珠接觸面的接觸剛度均接近,考慮到接觸剛度的設(shè)置,此處采用純罰函數(shù)法。

    (4)法向接觸剛度設(shè)置可采用調(diào)節(jié)接觸剛度因子或定義絕對(duì)接觸剛度值的方法,理論計(jì)算模型已計(jì)算出內(nèi)外圈接觸剛度系數(shù)Ki、Ko,因此,此處直接定義絕對(duì)接觸剛度值。其中:滾珠與內(nèi)圈的接觸設(shè)置為7.037×105N/mm3,外圈設(shè)置為6.259×105N/mm3,在有限元軟件中該參數(shù)的單位意義是接觸剛度對(duì)域的比值,其參數(shù)值可參考文獻(xiàn)[12]中有關(guān)內(nèi)外圈接觸剛度系數(shù)的計(jì)算。

    2.2.2 約束設(shè)置

    (1)由于軸承外圈固定,選擇其側(cè)面和外圓柱面施加Fixed Support全約束。

    (2)為模擬保持架對(duì)滾珠的約束作用,首先,需要根據(jù)每個(gè)滾珠所在的方位角,在內(nèi)圈圓心處建立多個(gè)沿內(nèi)圈滾道切向重合的局部坐標(biāo)系;然后,所有滾珠通過(guò)在對(duì)應(yīng)坐標(biāo)系上施加沿內(nèi)圈的切向位移約束,使?jié)L珠在仿真時(shí)只會(huì)產(chǎn)生軸向和徑向的位移變形。

    (3)將軸承內(nèi)圈內(nèi)環(huán)面的一半單元節(jié)點(diǎn)通過(guò)Remote Points的方式耦合到一中心節(jié)點(diǎn),用于模擬對(duì)耦合點(diǎn)施加徑向力下,軸承內(nèi)圈同轉(zhuǎn)軸配合時(shí)其內(nèi)圈實(shí)際受到的徑向力分布情況。

    2.2.3 加載設(shè)置

    將軸向載荷施加到軸承內(nèi)圈面部,徑向載荷施加到耦合的節(jié)點(diǎn)上,其中,軸向力與徑向力的實(shí)際施加方向分別為三維坐標(biāo)的-X與Y。

    如圖4所示,為軸承施加約束和載荷后的模型。

    圖4 施加約束和載荷后的模型

    2.3 求解設(shè)置

    求解的分析類型為準(zhǔn)靜態(tài)結(jié)構(gòu)分析,將求解方式設(shè)置為幾何非線性,開啟大變形選項(xiàng),并設(shè)置計(jì)算子步進(jìn)行求解。

    2.4 結(jié)果處理

    如圖5(a)、圖5(b)所示,分別為軸向力1 500 N、徑向力500 N下軸承軸向與徑向變形,在軸承內(nèi)圈背部的網(wǎng)格單元或節(jié)點(diǎn)添加探針拾取位移結(jié)果。圖中標(biāo)簽處即模擬實(shí)際試驗(yàn)過(guò)程中傳感器的探測(cè)位置,通過(guò)力除以位移坐標(biāo)差的方法可得到不同載荷下的剛度值。

    圖5 聯(lián)合載荷下的變形結(jié)果

    3 聯(lián)合加載仿真結(jié)果對(duì)比分析

    以7003C軸承為例,對(duì)軸承的理論值與有限元值進(jìn)行仿真結(jié)果對(duì)比分析,以模擬軸向預(yù)緊力下對(duì)軸承進(jìn)行加載的靜態(tài)試驗(yàn)[5]。

    3.1 軸向載荷對(duì)軸向變形及剛度的影響

    如圖6所示,為施加軸向力與產(chǎn)生的軸向變形的關(guān)系曲線??梢钥闯觯狠S向載荷(0~2 000 N)的增大會(huì)使得軸承軸向變形的理論值與有限元值均呈非線性增長(zhǎng);總體上軸向變形的有限元值略大于理論值。

    圖6 軸向載荷與軸承軸向變形的關(guān)系

    如圖7所示,為軸向載荷與軸向剛度的關(guān)系曲線。可以看出:圖中軸向載荷的增大會(huì)使得軸承軸向剛度的理論值與有限元值均呈非線性增長(zhǎng);軸向剛度的理論值略大于有限元值。

    圖7 軸向載荷與軸承軸向剛度的關(guān)系

    3.2 徑向載荷對(duì)徑向變形及剛度的影響

    如圖8所示,為在兩種軸向載荷條件下徑向載荷與徑向變形的關(guān)系圖??梢钥闯觯涸谳S向載荷分別為1 000 N和1 500 N時(shí),徑向載荷(0~500 N)的增大會(huì)使軸承徑向變形的理論值與有限元值均接近于線性增長(zhǎng);徑向變形的有限元值大于理論值,且所施加徑向載荷的增大會(huì)使差別變得明顯;軸向載荷為1 500 N時(shí)比1 000 N時(shí)的徑向變形小,并且軸向載荷越大,徑向載荷對(duì)徑向變形影響的理論值與有限元值更加接近。

    圖8 徑向載荷與軸承徑向變形的關(guān)系

    如圖9所示,為在兩種軸向載荷條件下徑向載荷與徑向剛度的關(guān)系圖??梢钥闯觯涸谳S向載荷分別為1 000 N和1 500 N時(shí),徑向載荷(0~500 N)的增大會(huì)使得徑向剛度的理論值呈非線性增長(zhǎng),而徑向剛度的有限元值先增長(zhǎng)后緩慢減小,說(shuō)明徑向力對(duì)徑向剛度的弱化影響在有限元仿真中較為明顯;各徑向載荷下的有限元值均小于理論值;相同徑向力時(shí),軸向載荷1 500 N下的徑向剛度大于軸向載荷1 000 N下的徑向剛度。

    圖9 徑向載荷與軸承徑向剛度的關(guān)系

    徑向剛度的兩種仿真結(jié)果出現(xiàn)差別其原因在于:一方面,有限元仿真時(shí)探測(cè)點(diǎn)的結(jié)果不僅有內(nèi)外圈與滾珠的接觸變形,還包括了軸承內(nèi)外圈自身的變形、內(nèi)外圈與滾珠間的相對(duì)滑移和接觸算法產(chǎn)生的穿透量;另一方面,由于徑向力的增加會(huì)導(dǎo)致軸承內(nèi)圈產(chǎn)生偏轉(zhuǎn),進(jìn)而使得徑向剛度出現(xiàn)弱化現(xiàn)象,本文的理論計(jì)算模型未能體現(xiàn)出這一點(diǎn),而在有限元仿真中,接觸滑移和接觸穿透會(huì)使得這種情況較為明顯,這也會(huì)導(dǎo)致探測(cè)點(diǎn)的結(jié)果與理論值有些偏差。

    3.3 軸向載荷對(duì)徑向變形的影響

    如圖10所示,為相同徑向載荷500 N下軸向載荷(1 000~2 000 N)對(duì)徑向變形影響的理論解析與有限元仿真結(jié)果??梢钥闯觯狠S向載荷的增加使得軸承的徑向變形的理論值與有限元值均呈非線性減小,且減小趨勢(shì)會(huì)逐漸變得緩慢;徑向位移的有限元值大于理論值并且趨近于理論值。

    圖10 軸向載荷與軸承徑向變形的關(guān)系

    3.4 徑向載荷對(duì)軸向變形的影響

    如圖11所示,為相同軸向載荷1 500 N下徑向載荷(0~1 000 N)對(duì)軸向變形影響的理論解析與有限元仿真結(jié)果??梢钥闯觯簭较蜉d荷的增加使得軸承的軸向變形的理論值與有限元值均呈非線性減小,且減小趨勢(shì)較為緩慢;軸向位移的有限元值大于理論值。

    圖11 徑向載荷與軸承軸向變形的關(guān)系

    對(duì)比文獻(xiàn)[2]中相同接觸角的角接觸球軸承相關(guān)試驗(yàn)結(jié)果,可以明顯看出,軸向與徑向變形的有限元仿真結(jié)果及理論計(jì)算值的變化趨勢(shì)同試驗(yàn)結(jié)果基本一致。因此,在工程實(shí)際應(yīng)用中,本文的軸承仿真方法具有一定的參考價(jià)值,其中,有限元法仿真的精確性,可通過(guò)對(duì)比工程試驗(yàn)結(jié)果并優(yōu)化模型設(shè)置來(lái)進(jìn)一步提高。

    4 結(jié)論

    本文通過(guò)理論模型與有限元仿真模型的建立,計(jì)算了7003C軸承的軸向變形與剛度以及徑向變形與剛度,并對(duì)比分析了聯(lián)合載荷下角接觸球軸承的理論解析與有限元仿真結(jié)果。得到以下結(jié)論:

    (1)軸向載荷的增大使得軸承的軸向變形與軸向剛度均呈非線性增長(zhǎng)的趨勢(shì);軸向變形的有限元值大于理論值,軸向剛度的有限元值小于理論值。

    (2)當(dāng)軸向載荷一定時(shí),徑向載荷的增大使軸承徑向變形的理論值與有限元值均接近于線性增長(zhǎng),其中,有限元仿真結(jié)果的規(guī)律更趨近于試驗(yàn);徑向剛度的理論值呈非線性增長(zhǎng),而徑向剛度的有限元值先增長(zhǎng)后緩慢減??;相同徑向力下的軸向載荷越大,徑向剛度也會(huì)越大。

    (3)軸向載荷的增加會(huì)使得軸承的徑向變形呈非線性減小,并且徑向載荷的增加也會(huì)使得軸承的軸向變形呈非線性減小。

    (4)角接觸球軸承軸向與徑向變形的有限元值均大于理論值,軸向與徑向剛度的有限元值均小于理論值。

    猜你喜歡
    理論值滾珠內(nèi)圈
    特種復(fù)合軸承內(nèi)圈推力滾道磨削用工裝設(shè)計(jì)
    哈爾濱軸承(2021年4期)2021-03-08 01:00:48
    滾珠絲杠的熱力耦合作用仿真
    擴(kuò)招百萬(wàn)背景下各省區(qū)高職院校新增招生規(guī)模測(cè)度研究
    主軸軸承內(nèi)圈鎖緊用臺(tái)階套的裝配
    組合變形實(shí)驗(yàn)中主應(yīng)力方位角理論值的確定
    ASME規(guī)范與JB/T4730對(duì)接焊縫超聲檢測(cè)的靈敏度差異探討
    滾珠絲杠的四軸加工
    內(nèi)圈帶缺陷中介軸承的動(dòng)力學(xué)建模與振動(dòng)響應(yīng)分析
    精密滾珠絲杠副伺服加載試驗(yàn)臺(tái)設(shè)計(jì)
    高速滾珠絲杠副動(dòng)態(tài)接觸特性求解
    天天躁夜夜躁狠狠躁躁| 大香蕉久久成人网| 18禁裸乳无遮挡动漫免费视频| 综合色丁香网| 日韩一区二区三区影片| 人妻 亚洲 视频| 全区人妻精品视频| 99久久中文字幕三级久久日本| 黑人欧美特级aaaaaa片| 精品卡一卡二卡四卡免费| 老司机亚洲免费影院| 99热6这里只有精品| 欧美成人精品欧美一级黄| 另类亚洲欧美激情| 日产精品乱码卡一卡2卡三| 女人精品久久久久毛片| 黄色毛片三级朝国网站| 亚洲精品456在线播放app| 午夜免费男女啪啪视频观看| freevideosex欧美| 国产视频首页在线观看| 亚洲欧美日韩另类电影网站| 国产精品不卡视频一区二区| 在线观看免费日韩欧美大片| 国产精品女同一区二区软件| 久热久热在线精品观看| 99国产精品免费福利视频| 1024视频免费在线观看| 亚洲伊人久久精品综合| 最后的刺客免费高清国语| 精品99又大又爽又粗少妇毛片| 亚洲久久久国产精品| 亚洲精品久久久久久婷婷小说| 免费久久久久久久精品成人欧美视频 | 久久久久视频综合| 久久国内精品自在自线图片| 成人综合一区亚洲| 国产熟女欧美一区二区| 亚洲国产欧美在线一区| 国产精品99久久99久久久不卡 | 99久久人妻综合| 亚洲精品美女久久av网站| av播播在线观看一区| 成人漫画全彩无遮挡| 热99国产精品久久久久久7| 男女午夜视频在线观看 | 亚洲国产精品成人久久小说| 久久精品aⅴ一区二区三区四区 | 一区二区av电影网| 日日摸夜夜添夜夜爱| 欧美成人精品欧美一级黄| 欧美日本中文国产一区发布| 亚洲av成人精品一二三区| av在线播放精品| 欧美亚洲日本最大视频资源| av有码第一页| 亚洲成国产人片在线观看| 午夜91福利影院| 亚洲图色成人| 成年美女黄网站色视频大全免费| 国产精品三级大全| 久久99热6这里只有精品| 亚洲欧美成人精品一区二区| 午夜av观看不卡| 色视频在线一区二区三区| 激情视频va一区二区三区| 欧美老熟妇乱子伦牲交| 欧美日韩精品成人综合77777| 观看美女的网站| 久久精品久久久久久久性| 欧美+日韩+精品| 亚洲美女黄色视频免费看| 成人影院久久| 女的被弄到高潮叫床怎么办| 女性生殖器流出的白浆| 18禁裸乳无遮挡动漫免费视频| 我要看黄色一级片免费的| 日本黄大片高清| 最近最新中文字幕免费大全7| 久久精品国产综合久久久 | 一级毛片电影观看| 成人二区视频| 熟女电影av网| 国产高清国产精品国产三级| 国产爽快片一区二区三区| 五月开心婷婷网| 日韩不卡一区二区三区视频在线| 国产精品熟女久久久久浪| 免费少妇av软件| 边亲边吃奶的免费视频| 黑人欧美特级aaaaaa片| 黄色 视频免费看| 国产一区二区三区综合在线观看 | 国产爽快片一区二区三区| 免费高清在线观看日韩| 午夜激情久久久久久久| 久久99热6这里只有精品| 在线精品无人区一区二区三| www.色视频.com| 最近手机中文字幕大全| 成人午夜精彩视频在线观看| 波多野结衣一区麻豆| 人人妻人人添人人爽欧美一区卜| 亚洲av在线观看美女高潮| 国产av精品麻豆| 最近最新中文字幕免费大全7| 老司机影院成人| 久久精品熟女亚洲av麻豆精品| 男的添女的下面高潮视频| 久久久精品94久久精品| 丝袜在线中文字幕| 肉色欧美久久久久久久蜜桃| 精品一区在线观看国产| 亚洲国产成人一精品久久久| 欧美xxⅹ黑人| 两个人看的免费小视频| 国产在线一区二区三区精| 成年美女黄网站色视频大全免费| 免费大片18禁| 高清在线视频一区二区三区| 久久久久久久久久久久大奶| 久久影院123| 成人二区视频| 国产免费福利视频在线观看| 久久久精品94久久精品| 亚洲av在线观看美女高潮| 永久网站在线| 黄片播放在线免费| 国产亚洲精品第一综合不卡 | 婷婷成人精品国产| 亚洲欧美精品自产自拍| 国产精品女同一区二区软件| 少妇精品久久久久久久| 亚洲精品,欧美精品| 久久青草综合色| 内地一区二区视频在线| 欧美精品一区二区大全| 日本爱情动作片www.在线观看| 亚洲伊人色综图| 亚洲第一区二区三区不卡| 97在线视频观看| 建设人人有责人人尽责人人享有的| 国产探花极品一区二区| 国产一区亚洲一区在线观看| 人妻少妇偷人精品九色| 国产在线免费精品| 亚洲伊人久久精品综合| a级毛色黄片| 伦理电影大哥的女人| 日韩中文字幕视频在线看片| 热re99久久精品国产66热6| 九草在线视频观看| 九色成人免费人妻av| 国产乱来视频区| 亚洲经典国产精华液单| 777米奇影视久久| 免费黄色在线免费观看| 成年女人在线观看亚洲视频| 久久久久久久久久久久大奶| 精品少妇内射三级| 国产精品国产av在线观看| 最黄视频免费看| 性色avwww在线观看| 久久久久久久亚洲中文字幕| 黄片无遮挡物在线观看| 欧美人与性动交α欧美软件 | 少妇猛男粗大的猛烈进出视频| 国产在线一区二区三区精| 天天操日日干夜夜撸| 黑人猛操日本美女一级片| 欧美成人精品欧美一级黄| 欧美日韩综合久久久久久| 亚洲综合色惰| 水蜜桃什么品种好| 在线精品无人区一区二区三| 精品国产一区二区久久| 在现免费观看毛片| 久久鲁丝午夜福利片| 色婷婷av一区二区三区视频| 美国免费a级毛片| 久久久久久久亚洲中文字幕| 美女主播在线视频| 国产 精品1| av视频免费观看在线观看| 街头女战士在线观看网站| 99久久综合免费| 亚洲国产看品久久| 日本色播在线视频| 下体分泌物呈黄色| 国产成人精品在线电影| av天堂久久9| 国产爽快片一区二区三区| 满18在线观看网站| 卡戴珊不雅视频在线播放| 免费人成在线观看视频色| av免费观看日本| 国内精品宾馆在线| 成人亚洲精品一区在线观看| 国产精品一区二区在线不卡| 精品一区二区三卡| 又大又黄又爽视频免费| 欧美成人午夜精品| 老司机影院毛片| 亚洲欧洲日产国产| 亚洲国产精品国产精品| 国产黄频视频在线观看| 老司机亚洲免费影院| 亚洲精品美女久久久久99蜜臀 | 色5月婷婷丁香| 亚洲国产av影院在线观看| 9热在线视频观看99| 国产成人精品福利久久| 中国国产av一级| 久久精品国产鲁丝片午夜精品| 亚洲精品自拍成人| 最近手机中文字幕大全| 伦理电影大哥的女人| 亚洲欧洲日产国产| 国产色爽女视频免费观看| 亚洲一区二区三区欧美精品| 9191精品国产免费久久| 在线精品无人区一区二区三| www.熟女人妻精品国产 | 欧美精品av麻豆av| 91成人精品电影| 自拍欧美九色日韩亚洲蝌蚪91| 性色av一级| 国产有黄有色有爽视频| 欧美日韩视频高清一区二区三区二| 免费黄频网站在线观看国产| 99香蕉大伊视频| 久久久久精品性色| 中文字幕人妻丝袜制服| 18禁在线无遮挡免费观看视频| 97超碰精品成人国产| 少妇人妻久久综合中文| 久久久久久人人人人人| 亚洲精品乱码久久久久久按摩| 久久精品国产综合久久久 | 性色avwww在线观看| 国产成人一区二区在线| 一级a做视频免费观看| 99久久人妻综合| 大香蕉久久成人网| 日韩 亚洲 欧美在线| 岛国毛片在线播放| 视频中文字幕在线观看| 亚洲精华国产精华液的使用体验| 国产一区亚洲一区在线观看| 新久久久久国产一级毛片| 亚洲精品,欧美精品| 欧美国产精品va在线观看不卡| 九草在线视频观看| 日韩成人av中文字幕在线观看| 午夜影院在线不卡| 日日撸夜夜添| 纯流量卡能插随身wifi吗| 永久免费av网站大全| 国产免费视频播放在线视频| 亚洲精品美女久久久久99蜜臀 | 国产一区二区激情短视频 | 亚洲五月色婷婷综合| 亚洲成人av在线免费| 午夜免费鲁丝| 国产精品人妻久久久久久| 爱豆传媒免费全集在线观看| 韩国精品一区二区三区 | 国产在线一区二区三区精| 中文欧美无线码| av女优亚洲男人天堂| 性高湖久久久久久久久免费观看| xxxhd国产人妻xxx| 欧美激情 高清一区二区三区| 亚洲成国产人片在线观看| 亚洲国产成人一精品久久久| 日本色播在线视频| 另类精品久久| freevideosex欧美| 少妇 在线观看| 在线天堂最新版资源| 高清不卡的av网站| 久久亚洲国产成人精品v| 国产国拍精品亚洲av在线观看| 亚洲美女搞黄在线观看| 亚洲熟女精品中文字幕| 欧美3d第一页| 汤姆久久久久久久影院中文字幕| 少妇 在线观看| 大话2 男鬼变身卡| 亚洲欧美清纯卡通| 国产精品久久久久久久电影| 大片免费播放器 马上看| 亚洲性久久影院| 国国产精品蜜臀av免费| 亚洲精品久久成人aⅴ小说| 极品人妻少妇av视频| 国产无遮挡羞羞视频在线观看| 国产成人精品无人区| 最后的刺客免费高清国语| 国产一区二区激情短视频 | 黄色 视频免费看| 婷婷色麻豆天堂久久| av不卡在线播放| 看十八女毛片水多多多| 久久韩国三级中文字幕| 亚洲国产日韩一区二区| av在线播放精品| 亚洲国产精品999| 九色成人免费人妻av| tube8黄色片| 精品一区二区三区四区五区乱码 | 视频在线观看一区二区三区| 欧美日韩av久久| 美女脱内裤让男人舔精品视频| 免费观看性生交大片5| 日韩在线高清观看一区二区三区| 久久久久人妻精品一区果冻| 久久亚洲国产成人精品v| 国国产精品蜜臀av免费| av视频免费观看在线观看| 亚洲国产欧美在线一区| 成人综合一区亚洲| 少妇的逼好多水| 午夜影院在线不卡| 国产成人精品婷婷| 成年动漫av网址| xxx大片免费视频| 国产亚洲精品第一综合不卡 | 两个人免费观看高清视频| 欧美精品国产亚洲| 久久精品熟女亚洲av麻豆精品| 各种免费的搞黄视频| 韩国精品一区二区三区 | 男人添女人高潮全过程视频| 高清av免费在线| 国产日韩欧美在线精品| 国产欧美另类精品又又久久亚洲欧美| 黑人巨大精品欧美一区二区蜜桃 | 午夜福利视频在线观看免费| 少妇被粗大猛烈的视频| 18禁动态无遮挡网站| 在线观看免费日韩欧美大片| 成人二区视频| 国产一区二区在线观看日韩| 久久久久久人妻| 亚洲美女黄色视频免费看| 国产精品一二三区在线看| 亚洲一区二区三区欧美精品| a级毛色黄片| 人成视频在线观看免费观看| 久久国内精品自在自线图片| 国产免费一级a男人的天堂| 亚洲四区av| 成年美女黄网站色视频大全免费| 精品久久久精品久久久| 成年美女黄网站色视频大全免费| 精品一区二区三区视频在线| 久久免费观看电影| 国产精品99久久99久久久不卡 | 老司机影院毛片| 亚洲成人手机| 美女脱内裤让男人舔精品视频| 亚洲婷婷狠狠爱综合网| 黑丝袜美女国产一区| 91在线精品国自产拍蜜月| 国产伦理片在线播放av一区| 亚洲精品久久午夜乱码| 女性生殖器流出的白浆| 欧美日韩国产mv在线观看视频| 国产高清不卡午夜福利| 热re99久久国产66热| 久久久国产精品麻豆| 久久青草综合色| 在线天堂中文资源库| 啦啦啦中文免费视频观看日本| 亚洲欧美成人精品一区二区| 天美传媒精品一区二区| 国产一区二区三区综合在线观看 | 国产免费福利视频在线观看| 9热在线视频观看99| 久久久久久久久久久久大奶| 一边摸一边做爽爽视频免费| 亚洲色图 男人天堂 中文字幕 | 一级,二级,三级黄色视频| 九色亚洲精品在线播放| 日韩人妻精品一区2区三区| 亚洲欧美中文字幕日韩二区| 久久久亚洲精品成人影院| www日本在线高清视频| 狠狠婷婷综合久久久久久88av| 我要看黄色一级片免费的| 寂寞人妻少妇视频99o| 2021少妇久久久久久久久久久| 中文字幕人妻丝袜制服| 欧美日韩av久久| 男女边摸边吃奶| 精品一区二区三区视频在线| 亚洲精品日韩在线中文字幕| 久久久精品区二区三区| 免费观看a级毛片全部| 国产欧美日韩一区二区三区在线| 久久久久精品久久久久真实原创| 乱码一卡2卡4卡精品| 少妇的丰满在线观看| 欧美精品高潮呻吟av久久| 色吧在线观看| 青春草视频在线免费观看| 日韩av免费高清视频| 一区二区日韩欧美中文字幕 | 国产一区有黄有色的免费视频| 免费不卡的大黄色大毛片视频在线观看| 国产又爽黄色视频| 91久久精品国产一区二区三区| 日韩成人伦理影院| 青青草视频在线视频观看| 久久久久久久久久成人| 在线观看国产h片| 亚洲精品美女久久久久99蜜臀 | 午夜福利网站1000一区二区三区| a 毛片基地| 久久久亚洲精品成人影院| 国产乱来视频区| 波多野结衣一区麻豆| 国产爽快片一区二区三区| 美女国产视频在线观看| 免费人成在线观看视频色| 狂野欧美激情性xxxx在线观看| 18禁国产床啪视频网站| 看免费av毛片| 夫妻性生交免费视频一级片| 中文乱码字字幕精品一区二区三区| 成人毛片a级毛片在线播放| 99热6这里只有精品| 制服诱惑二区| 女的被弄到高潮叫床怎么办| 日本与韩国留学比较| 亚洲精品国产av蜜桃| 亚洲av.av天堂| 欧美精品高潮呻吟av久久| 亚洲欧洲日产国产| 国产探花极品一区二区| 日韩欧美精品免费久久| 亚洲av日韩在线播放| 在线免费观看不下载黄p国产| 只有这里有精品99| 日韩制服丝袜自拍偷拍| 777米奇影视久久| 热99久久久久精品小说推荐| 亚洲欧美日韩卡通动漫| 性色av一级| 欧美精品人与动牲交sv欧美| 亚洲精品成人av观看孕妇| 国产精品成人在线| 亚洲五月色婷婷综合| 欧美亚洲 丝袜 人妻 在线| 天堂俺去俺来也www色官网| 26uuu在线亚洲综合色| 久久精品国产鲁丝片午夜精品| av片东京热男人的天堂| 免费人成在线观看视频色| 精品第一国产精品| 性高湖久久久久久久久免费观看| 国产精品嫩草影院av在线观看| 这个男人来自地球电影免费观看 | av女优亚洲男人天堂| 欧美日韩成人在线一区二区| 狠狠精品人妻久久久久久综合| 日韩大片免费观看网站| 欧美97在线视频| 亚洲av电影在线观看一区二区三区| 九草在线视频观看| 欧美国产精品va在线观看不卡| 亚洲一区二区三区欧美精品| 久久久久久人妻| 亚洲国产欧美在线一区| 人人妻人人爽人人添夜夜欢视频| 王馨瑶露胸无遮挡在线观看| 母亲3免费完整高清在线观看 | 日日啪夜夜爽| 日本猛色少妇xxxxx猛交久久| 深夜精品福利| 一本久久精品| 最新的欧美精品一区二区| 亚洲精品国产av蜜桃| 精品久久蜜臀av无| 18禁国产床啪视频网站| 精品卡一卡二卡四卡免费| 午夜福利网站1000一区二区三区| 亚洲精品中文字幕在线视频| 99久久中文字幕三级久久日本| 80岁老熟妇乱子伦牲交| 一二三四中文在线观看免费高清| 高清黄色对白视频在线免费看| 秋霞在线观看毛片| 亚洲国产av影院在线观看| 成年动漫av网址| 精品少妇久久久久久888优播| 国产熟女欧美一区二区| 国产精品国产三级国产专区5o| 亚洲美女视频黄频| 精品人妻在线不人妻| 国产av国产精品国产| 全区人妻精品视频| 韩国高清视频一区二区三区| 日日撸夜夜添| 97超碰精品成人国产| 大片免费播放器 马上看| 中文天堂在线官网| 午夜精品国产一区二区电影| 高清av免费在线| 不卡视频在线观看欧美| 国产精品熟女久久久久浪| 91aial.com中文字幕在线观看| 亚洲精品一区蜜桃| 欧美精品av麻豆av| 成人亚洲欧美一区二区av| 免费黄色在线免费观看| 人人妻人人澡人人看| h视频一区二区三区| 日韩成人av中文字幕在线观看| 欧美日韩一区二区视频在线观看视频在线| 少妇熟女欧美另类| 国产一区二区在线观看日韩| 日韩精品免费视频一区二区三区 | 日韩一区二区视频免费看| 亚洲欧美日韩卡通动漫| 成人毛片a级毛片在线播放| 女的被弄到高潮叫床怎么办| 日韩av免费高清视频| 岛国毛片在线播放| 黄色视频在线播放观看不卡| 最黄视频免费看| 日本欧美视频一区| 人人妻人人添人人爽欧美一区卜| 欧美xxⅹ黑人| 亚洲国产欧美日韩在线播放| 国产一区亚洲一区在线观看| 欧美+日韩+精品| 9191精品国产免费久久| 蜜臀久久99精品久久宅男| 乱码一卡2卡4卡精品| 欧美性感艳星| 成年美女黄网站色视频大全免费| 国产成人欧美| 亚洲av福利一区| 久久久久国产网址| 国产片特级美女逼逼视频| 日本色播在线视频| 成年动漫av网址| 国产成人a∨麻豆精品| 亚洲av欧美aⅴ国产| 久久人人97超碰香蕉20202| 日韩av免费高清视频| av有码第一页| 少妇 在线观看| av片东京热男人的天堂| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 一级毛片我不卡| 在线观看免费日韩欧美大片| 国产精品久久久av美女十八| 免费观看性生交大片5| 精品福利永久在线观看| 国产深夜福利视频在线观看| 欧美另类一区| 99视频精品全部免费 在线| 欧美丝袜亚洲另类| 久久久久视频综合| 黄色配什么色好看| 亚洲国产精品成人久久小说| 精品国产一区二区久久| 熟女av电影| 亚洲av免费高清在线观看| 波多野结衣一区麻豆| 精品一品国产午夜福利视频| 啦啦啦视频在线资源免费观看| 精品福利永久在线观看| 国产在线一区二区三区精| 两个人看的免费小视频| 久久99一区二区三区| 熟妇人妻不卡中文字幕| 国产亚洲午夜精品一区二区久久| 成人亚洲欧美一区二区av| 热re99久久精品国产66热6| 一级片'在线观看视频| www.av在线官网国产| 水蜜桃什么品种好| 少妇人妻久久综合中文| 久久久久久久久久成人| 欧美日韩av久久| 美女国产视频在线观看| 亚洲久久久国产精品| 欧美日韩av久久| 美女国产视频在线观看| 十八禁网站网址无遮挡| 日韩 亚洲 欧美在线| 色吧在线观看| 黄色配什么色好看| 久久人人爽人人片av| 五月玫瑰六月丁香| 久久精品人人爽人人爽视色| 激情视频va一区二区三区| 一边亲一边摸免费视频| 99精国产麻豆久久婷婷| 飞空精品影院首页| 国产亚洲最大av| 日本免费在线观看一区| 国产欧美日韩一区二区三区在线| 中文字幕最新亚洲高清| 在线看a的网站| 内地一区二区视频在线| 9色porny在线观看| 亚洲av欧美aⅴ国产| 老女人水多毛片| 卡戴珊不雅视频在线播放| 久久热在线av|