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

    小半徑曲線鋼軌波磨預(yù)測(cè)模型及在波磨抑制中的應(yīng)用*

    2019-07-02 11:49:06
    潤(rùn)滑與密封 2019年6期
    關(guān)鍵詞:波磨阻尼比輪軌

    (西南交通大學(xué)摩擦學(xué)研究所,四川成都 610031)

    隨著我國(guó)科學(xué)技術(shù)的發(fā)展,近幾年來(lái),我國(guó)在軌道交通運(yùn)輸領(lǐng)域取得了卓越的成就。但在鐵路運(yùn)行過(guò)程中,出現(xiàn)的一系列問(wèn)題也困擾著鐵路工作人員和相關(guān)學(xué)者,其中鋼軌波磨就是最突出問(wèn)題之一。特別是在城軌交通蓬勃發(fā)展之際,鋼軌波磨問(wèn)題越來(lái)越嚴(yán)重,嚴(yán)重的時(shí)候,一條新地鐵線投入使用不到2~3個(gè)月就出現(xiàn)鋼軌波磨。自從出現(xiàn)鋼軌波磨以來(lái),各國(guó)學(xué)者進(jìn)行了大量研究,提出的鋼軌波磨形成機(jī)制主要分為兩大類[1-2]:(1)輪軌表面粗糙度導(dǎo)致的輪軌瞬態(tài)動(dòng)力學(xué)相互作用引起摩擦功波動(dòng)產(chǎn)生鋼軌波磨[3];(2)由于輪軌間黏-滑導(dǎo)致自激振動(dòng)產(chǎn)生鋼軌波磨[4-6]。目前大部分研究者接受第一種理論[7-8]。

    盡管對(duì)鋼軌波磨的研究已有長(zhǎng)達(dá)100多年的歷史,但時(shí)至今日,占全世界鋼軌波磨總量80%以上的小半徑曲線內(nèi)軌波磨還是近乎百分之百發(fā)生,沒(méi)有預(yù)期的減少。這個(gè)尷尬的局面可能與目前大家認(rèn)知的鋼軌波磨機(jī)制存在明顯的缺陷有關(guān)。如現(xiàn)有的鋼軌波磨理論不能解釋為什么直線或者大半徑曲線上的左右2根鋼軌,以及小半徑曲線外軌難以發(fā)生波磨。為此,陳光雄教授及所在的團(tuán)隊(duì)提出了摩擦耦合自激振動(dòng)引發(fā)波磨的新觀點(diǎn)[9-12],該觀點(diǎn)認(rèn)為當(dāng)輪軌間蠕滑力達(dá)到飽和狀態(tài)時(shí),這個(gè)飽和蠕滑力容易引起輪軌系統(tǒng)的摩擦耦合自激振動(dòng),從而導(dǎo)致鋼軌波磨。本文作者根據(jù)摩擦自激振動(dòng)引起鋼軌波磨的觀點(diǎn),建立地鐵小半徑曲線軌道鋼軌波磨預(yù)測(cè)模型,運(yùn)用復(fù)特征值分析方法探討輪對(duì)采用過(guò)盈配合模擬方式和整體輪對(duì)模擬方式時(shí)輪軌系統(tǒng)運(yùn)動(dòng)穩(wěn)定性,并研究單側(cè)輪軌摩擦因數(shù)對(duì)鋼軌波磨的影響。

    1 輪軌系統(tǒng)摩擦耦合自激振動(dòng)模型

    1.1 輪軌系統(tǒng)模型

    通過(guò)對(duì)地鐵線路現(xiàn)場(chǎng)調(diào)查發(fā)現(xiàn),在曲線半徑R≤350 m的線路上,低軌都會(huì)出現(xiàn)波磨,而在高軌上幾乎沒(méi)有波磨產(chǎn)生。車輛在小半徑曲線線路上通過(guò)時(shí),導(dǎo)向輪對(duì)和鋼軌間的蠕滑力一般趨于飽和,輪軌間蠕滑力等于法向力乘以動(dòng)摩擦因數(shù),從而引起輪軌摩擦自激振動(dòng),導(dǎo)致鋼軌波磨的產(chǎn)生[10]。所以文中建立的是地鐵線路小半徑曲線軌道上只包含導(dǎo)向輪對(duì)輪軌系統(tǒng)的摩擦耦合自激振動(dòng)模型。其接觸模型如圖1所示。

    圖1 輪軌系統(tǒng)接觸模型

    通過(guò)調(diào)查實(shí)際線路上車輛運(yùn)行情況,并結(jié)合Simpack動(dòng)力學(xué)仿真計(jì)算可得,當(dāng)?shù)罔F車輛以65 km/h速度通過(guò)半徑為350 m曲線線路時(shí),軸箱垂向懸掛力FSVL=44.5 kN,F(xiàn)SVR=38.7 kN;橫向懸掛力FSLL=FSLR=5.9 kN,低軌與內(nèi)輪接觸角δR=2.1°,高軌和外輪接觸角δL=33°。鋼軌彈簧支撐的垂向剛度KRV=40.73 MN/m,橫向剛度KRL=8.79 MN/m;垂向阻尼和橫向阻尼分別為CRV=9 898.70 N·s/m和CRL=1 927.96 N·s/m。軌枕底部聯(lián)合支撐垂向和橫向剛度分別為KSV=89 MN/m和KSL=50 MN/m,垂向和橫向阻尼分別為CSV=89 800 N·s/m和CSL=40 000 N·s/m[14]。根據(jù)輪軌接觸模型,運(yùn)用ABAQUS軟件建立有限元模型如圖2所示,其中鋼軌選用長(zhǎng)為36 m的60 kg/m鋼軌;軌枕間距為625 mm;軌底坡為1/40;車輪與鋼軌間摩擦因數(shù)為0.45;輪對(duì)及鋼軌材料密度為7 800 kg/m3,彈性模量為210 GPa,泊松比為0.3。

    圖2 輪軌系統(tǒng)有限元模型

    1.2 輪軌系統(tǒng)摩擦自激振動(dòng)方程

    輪軌接觸時(shí)存在著輪軌摩擦相對(duì)滑動(dòng)作用,運(yùn)用ABAQUS軟件進(jìn)行復(fù)特征值分析時(shí),在穩(wěn)態(tài)滑動(dòng)平衡位置建立的運(yùn)動(dòng)微分方程[10,14]如下:

    Mx″+Cx′+Kx=0

    (1)

    式中:x為節(jié)點(diǎn)位移;M、C、K分別為質(zhì)量矩陣、阻尼矩陣和剛度矩陣。

    沒(méi)有摩擦?xí)r,方程(1)中的M、C、K都是對(duì)稱矩陣,所以方程(1)的特征值的實(shí)部Rm?0,系統(tǒng)處于穩(wěn)定運(yùn)動(dòng)狀態(tài)。當(dāng)考慮輪軌間的摩擦力之后會(huì)導(dǎo)致M、C和K變?yōu)榉菍?duì)稱矩陣,此時(shí)方程(1)的特征方程為

    (λ2M+λC+K)φ=0

    (2)

    式中:λ為系統(tǒng)運(yùn)動(dòng)方程特征值;φ為特征值對(duì)應(yīng)的特征向量。

    根據(jù)特征方程,求得其通解為

    (3)

    式中:t為時(shí)間,αi+jwi=λi為復(fù)特征值。

    運(yùn)用復(fù)特征值法分析系統(tǒng)摩擦自激振動(dòng)時(shí),根據(jù)等效阻尼比來(lái)評(píng)價(jià)不穩(wěn)定振動(dòng)發(fā)生的趨勢(shì)。等效阻尼比的數(shù)值小于0時(shí),系統(tǒng)才會(huì)發(fā)生摩擦自激振動(dòng),而且數(shù)值越小,系統(tǒng)發(fā)生摩擦自激振動(dòng)的趨勢(shì)越大。其計(jì)算公式為

    ξ=-α/(π|ω|)

    (4)

    式中:α為特征值實(shí)部;ω為特征值虛部。

    2 計(jì)算結(jié)果及分析

    2.1 網(wǎng)絡(luò)劃分

    考慮到模型中單元大小對(duì)模擬結(jié)果的影響,文中分析了輪軌接觸區(qū)單元尺寸變化對(duì)模型模擬結(jié)果的影響。根據(jù)輪軌系統(tǒng)模型規(guī)模大小,在輪軌系統(tǒng)有限元模型中,車輪踏面設(shè)為主面,輪軌接觸區(qū)單元平均長(zhǎng)度約為12 mm,鋼軌頂面設(shè)為從面,接觸區(qū)單元平均長(zhǎng)度約為8 mm。輪軌接觸區(qū)有限元網(wǎng)格模型如圖3所示。文中輪軌接觸區(qū)單元尺寸分別取上述單元尺寸平均長(zhǎng)度的1/2、1/3和1/4進(jìn)行計(jì)算,研究了單元細(xì)化過(guò)程中,輪軌系統(tǒng)不穩(wěn)定振動(dòng)主頻率的變化情況。

    圖3 輪軌接觸區(qū)網(wǎng)格模型

    圖4所示為不同單元尺寸下輪軌系統(tǒng)對(duì)應(yīng)不穩(wěn)定振動(dòng)主頻率變化情況,可知,當(dāng)輪軌接觸區(qū)單元尺寸為原有尺寸的1/2時(shí),計(jì)算得出輪軌系統(tǒng)不穩(wěn)定振動(dòng)主頻率為499.71 Hz;當(dāng)單元尺寸減小為原尺寸的1/3時(shí),不穩(wěn)定振動(dòng)主頻率為499.66 Hz;當(dāng)單元尺寸減小為原尺寸的1/4時(shí),不穩(wěn)定振動(dòng)主頻率為499.5 Hz。

    圖4 單元尺寸對(duì)不穩(wěn)定振動(dòng)主頻率的影響

    通過(guò)對(duì)比發(fā)現(xiàn),細(xì)化單元尺寸前后,模型計(jì)算出的不穩(wěn)定振動(dòng)主頻率結(jié)果很相近,相對(duì)誤差均小于0.05%,頻率大小保持不變趨勢(shì)。由此說(shuō)明輪軌接觸區(qū)單元尺寸變化對(duì)文中有限元模型模擬結(jié)果影響較小,因此為減少計(jì)算成本,有限元模型輪軌接觸面采用初始的網(wǎng)格尺寸就可以了。

    2.2 整體輪的輪軌摩擦自激振動(dòng)分析

    通過(guò)對(duì)整體輪的輪軌系統(tǒng)進(jìn)行復(fù)特征值分析,可以得到輪軌系統(tǒng)負(fù)等效阻尼比分布情況和主要不穩(wěn)定振動(dòng)模態(tài)。如圖5所示,輪軌系統(tǒng)存在2個(gè)負(fù)等效阻尼比,其分別為-0.029 37和-0.008 36。2個(gè)負(fù)等效阻尼比對(duì)應(yīng)的不穩(wěn)定振動(dòng)模態(tài)圖如圖6所示。

    圖5 整體式輪對(duì)模式下輪軌系統(tǒng)等效阻尼比分布

    圖6 整體式輪對(duì)模式下輪軌系統(tǒng)自激振動(dòng)模態(tài)

    由于等效阻尼比越小,發(fā)生不穩(wěn)定振動(dòng)趨勢(shì)越大,故等效阻尼比為-0.029 37時(shí),越容易發(fā)生不穩(wěn)定自激振動(dòng),其對(duì)應(yīng)的頻率為499.50 Hz,即輪對(duì)采用整體式模擬方式預(yù)測(cè)得到的波磨頻率為499.5 Hz左右。

    由圖6(a)可知,在最容易發(fā)生不穩(wěn)定振動(dòng)頻率下,輪軌系統(tǒng)摩擦自激振動(dòng)主要發(fā)生在低軌和內(nèi)側(cè)車輪上,從而在低軌上產(chǎn)生波磨,而高軌上沒(méi)有波磨。

    2.3 輪軸過(guò)盈配合時(shí)輪軌摩擦自激振動(dòng)分析

    在實(shí)際設(shè)計(jì)生產(chǎn)中,車軸和車輪通常通過(guò)壓裝法過(guò)盈裝配在一起,輪軸之間存在一定的過(guò)盈量?;谳嗆壞Σ榴詈献约ふ駝?dòng)引起鋼軌波磨的觀點(diǎn)建立的小半徑曲線上的輪軌有限元模型,是將車輪和車軸簡(jiǎn)化為一個(gè)整體,因而忽略了輪軸之間力的作用和摩擦接觸屬性,模擬方法有別于實(shí)際運(yùn)行工況。為了更好地研究鋼軌波磨的形成機(jī)制,在原有的輪軌系統(tǒng)有限元模型之上將車軸和車輪分離開(kāi)來(lái),建立了帶有摩擦接觸屬性且設(shè)置過(guò)盈配合的輪軸接觸關(guān)系。輪軸接觸模型如圖7所示。

    圖7 輪軸接觸模型

    車輪和車軸接觸參數(shù)設(shè)置如下,其中輪座直徑為180 mm,一般來(lái)說(shuō)輪軸配合過(guò)盈量為輪座直徑的0.08%~0.15%[15],即輪軸過(guò)盈量取值為0.144~0.27 mm。取輪軸過(guò)盈量為0.24 mm;輪座和輪轂接觸面摩擦因數(shù)為0.15,輪軌系統(tǒng)其他參數(shù)與上文參數(shù)設(shè)置相同。

    通過(guò)設(shè)置的模型參數(shù),計(jì)算得到的負(fù)等效阻尼比的分布情況如圖8所示。

    從圖8可以看出,當(dāng)輪對(duì)采用過(guò)盈配合進(jìn)行模擬時(shí),輪軌系統(tǒng)可能產(chǎn)生三階不穩(wěn)定振動(dòng)。3個(gè)小于0的等效阻尼比分別為-0.031 97、-0.027 68和-0.011 60,其對(duì)應(yīng)的頻率分別為322.99、502.06和536.22 Hz。對(duì)比等效阻尼比數(shù)值可以發(fā)現(xiàn),其中等效阻尼比-0.031 97最小,其不穩(wěn)定振動(dòng)頻率為322.99 Hz,該頻率即為輪軌系統(tǒng)主要不穩(wěn)定振動(dòng)頻率,所以該輪軌系統(tǒng)模型預(yù)測(cè)出的鋼軌波磨的頻率為322.99 Hz。

    圖8 輪軸過(guò)盈配合時(shí)輪軌系統(tǒng)自激振動(dòng)頻率分布

    為了探究輪軸采用過(guò)盈配合時(shí)有限元模型預(yù)測(cè)鋼軌波磨出現(xiàn)特點(diǎn)是否與實(shí)際線路波磨發(fā)生特點(diǎn)一樣,即在小半徑曲線線路上,外側(cè)鋼軌幾乎無(wú)波磨,內(nèi)側(cè)鋼軌幾乎都有波磨,文中提取了3個(gè)不穩(wěn)定振動(dòng)頻率下的振動(dòng)模態(tài),輪軌系統(tǒng)自激振動(dòng)模態(tài)如圖9所示。

    圖9 輪軸過(guò)盈配合時(shí)輪軌系統(tǒng)自激振動(dòng)模態(tài)

    由圖9(a)可知,在最容易發(fā)生不穩(wěn)定振動(dòng)頻率下,輪軌系統(tǒng)摩擦自激振動(dòng)主要發(fā)生在低軌和內(nèi)側(cè)車輪上,從而在低軌上產(chǎn)生波磨,而高軌上沒(méi)有波磨。這與實(shí)際線路上波磨發(fā)生的特點(diǎn)相吻合。

    2.4 兩種模擬方式的結(jié)果對(duì)比分析

    在實(shí)際地鐵線路中,在小半徑(R≤350 m)曲線普通軌道上,鋼軌波磨波長(zhǎng)普遍為40~50 mm?,F(xiàn)場(chǎng)調(diào)研顯示車輛通過(guò)該區(qū)段速度為55~65 km/h,由速度為頻率與波長(zhǎng)的積計(jì)算可知,普通軌道線路上波磨頻率范圍為305.56~451.38 Hz。整體輪對(duì)模型計(jì)算得出輪軌系統(tǒng)主要不穩(wěn)定振動(dòng)頻率為499.5 Hz,即該模型預(yù)測(cè)波磨頻率為499.5 Hz左右,相比實(shí)測(cè)數(shù)據(jù)推算結(jié)果略偏大。當(dāng)輪對(duì)采用過(guò)盈配合方式模擬時(shí),該模型預(yù)測(cè)波磨頻率為322.99 Hz左右,完全符合現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)推算的波磨頻率,并且振動(dòng)模態(tài)圖和線路上波磨發(fā)生特點(diǎn)完全吻合。

    綜合對(duì)比輪對(duì)2種不同模擬方法可以知道,輪軸采用過(guò)盈配合模擬更加符合實(shí)際情況,輪軌系統(tǒng)主要不穩(wěn)定振動(dòng)頻率與現(xiàn)場(chǎng)實(shí)測(cè)波磨頻率更加吻合;相比整體輪對(duì)模擬方式,其計(jì)算結(jié)果準(zhǔn)確性更高,模擬結(jié)果更好。

    3 抑制鋼軌波磨發(fā)生的方法

    EADIE等[16-17]提出利用摩擦調(diào)節(jié)劑控制輪軌摩擦因數(shù)可以消除鋼軌波磨,但分別討論兩側(cè)輪軌摩擦因數(shù)影響的研究較少。為了更好反映摩擦因數(shù)與鋼軌波磨關(guān)系,文中利用上文建立的高精度鋼軌波磨預(yù)測(cè)模型,分別研究了內(nèi)側(cè)輪軌和外側(cè)輪軌不同摩擦因數(shù)對(duì)鋼軌波磨的影響。在小半徑曲線軌道上輪軌系統(tǒng)主要不穩(wěn)定振動(dòng)頻率約為322 Hz,實(shí)際線路中輪軌摩擦因數(shù)一般為0.2~0.6。文中兩側(cè)輪軌摩擦因數(shù)取值0.2~0.6進(jìn)行討論,結(jié)果如圖10所示。

    圖10 摩擦因數(shù)對(duì)摩擦自激振動(dòng)的影響

    如圖10所示,當(dāng)改變外側(cè)輪軌摩擦因數(shù)時(shí),系統(tǒng)等效阻尼比變化較小,即外側(cè)輪軌摩擦因數(shù)對(duì)鋼軌波磨影響較小;當(dāng)逐漸增大內(nèi)側(cè)輪軌摩擦因數(shù)時(shí),系統(tǒng)等效阻尼比變化較為明顯且呈現(xiàn)減小趨勢(shì),因此,內(nèi)側(cè)輪軌摩擦因數(shù)對(duì)小半徑曲線上輪軌系統(tǒng)不穩(wěn)定振動(dòng)影響較大,適當(dāng)減小內(nèi)側(cè)輪軌摩擦因數(shù)有助于抑制鋼軌波磨。

    4 結(jié)論

    (1)輪軌系統(tǒng)中輪對(duì)采用過(guò)盈配合模擬方式時(shí),模擬得到主要不穩(wěn)定振動(dòng)頻率為322.99 Hz,該主要不穩(wěn)定發(fā)生在低軌和內(nèi)輪上,這與現(xiàn)場(chǎng)實(shí)測(cè)鋼軌波磨情況更為一致。相比整體輪對(duì)模擬方式,輪對(duì)采用分開(kāi)式模擬更加符合實(shí)際設(shè)計(jì)生產(chǎn)情況,模擬結(jié)果精度更高,該模擬方式更加有利于進(jìn)行鋼軌波磨相關(guān)研究工作。

    (2) 在小半徑曲線軌道上,外側(cè)輪軌摩擦因數(shù)對(duì)波磨影響較小,內(nèi)側(cè)輪軌摩擦因數(shù)比外輪軌摩擦因數(shù)對(duì)波磨的影響更為明顯,適當(dāng)減小內(nèi)側(cè)輪軌摩擦因數(shù)有利于抑制鋼軌波磨。

    猜你喜歡
    波磨阻尼比輪軌
    城軌小半徑曲線鋼軌波磨的成因與整治研究綜述
    北京地鐵鋼軌波磨現(xiàn)狀及防治措施
    基于細(xì)觀結(jié)構(gòu)的原狀黃土動(dòng)彈性模量和阻尼比試驗(yàn)研究
    地震研究(2021年1期)2021-04-13 01:05:24
    中低速磁浮道岔與輪軌道岔的差異
    黏滯阻尼器在時(shí)程分析下的附加有效阻尼比研究
    波形分析法求解公路橋梁阻尼比的探討
    上海公路(2018年3期)2018-03-21 05:55:48
    結(jié)構(gòu)構(gòu)件阻尼比對(duì)大跨度懸索橋地震響應(yīng)的影響
    中低速磁浮與輪軌交通信號(hào)系統(tǒng)的差異
    非線性穩(wěn)態(tài)曲線通過(guò)時(shí)輪軌滾動(dòng)接觸的數(shù)值求解方法
    高速鐵路鋼軌波磨對(duì)車輛—軌道動(dòng)態(tài)響應(yīng)的影響
    视频在线观看一区二区三区| 日本a在线网址| 最新的欧美精品一区二区| 女人精品久久久久毛片| 日日夜夜操网爽| a级毛片黄视频| 日韩制服丝袜自拍偷拍| 啦啦啦 在线观看视频| 啪啪无遮挡十八禁网站| 精品熟女少妇八av免费久了| 久久久久精品国产欧美久久久 | 两个人免费观看高清视频| 欧美激情极品国产一区二区三区| 欧美日韩成人在线一区二区| 欧美另类亚洲清纯唯美| 99九九在线精品视频| 啦啦啦啦在线视频资源| 午夜福利一区二区在线看| 国产av精品麻豆| 免费观看a级毛片全部| 国产精品欧美亚洲77777| 咕卡用的链子| 日韩制服丝袜自拍偷拍| videosex国产| 久久久国产成人免费| 多毛熟女@视频| 中文字幕精品免费在线观看视频| 天天躁夜夜躁狠狠躁躁| 日本精品一区二区三区蜜桃| 亚洲av国产av综合av卡| 夫妻午夜视频| 又黄又粗又硬又大视频| 日本五十路高清| 欧美黄色片欧美黄色片| 18在线观看网站| 自拍欧美九色日韩亚洲蝌蚪91| 午夜91福利影院| 午夜视频精品福利| 国产三级黄色录像| 69精品国产乱码久久久| 免费在线观看日本一区| 老汉色∧v一级毛片| 午夜免费成人在线视频| 婷婷色av中文字幕| 精品福利观看| 亚洲欧美精品综合一区二区三区| av天堂在线播放| 国产欧美日韩精品亚洲av| 男女边摸边吃奶| 亚洲国产欧美在线一区| 精品国产国语对白av| 黄色 视频免费看| 国产av一区二区精品久久| 免费观看人在逋| 国产国语露脸激情在线看| 精品国产乱子伦一区二区三区 | 如日韩欧美国产精品一区二区三区| 国产亚洲欧美精品永久| 亚洲欧美一区二区三区久久| 一级毛片女人18水好多| 国产亚洲av高清不卡| avwww免费| 黄频高清免费视频| 日韩中文字幕视频在线看片| 女人爽到高潮嗷嗷叫在线视频| 久久人人97超碰香蕉20202| 亚洲第一青青草原| 国产黄色免费在线视频| 久久影院123| a 毛片基地| 日本欧美视频一区| 精品国产一区二区三区久久久樱花| 欧美日韩一级在线毛片| xxxhd国产人妻xxx| www.自偷自拍.com| 欧美激情极品国产一区二区三区| 一区在线观看完整版| 中文字幕最新亚洲高清| 久久精品亚洲熟妇少妇任你| 十八禁高潮呻吟视频| 欧美日韩中文字幕国产精品一区二区三区 | 午夜激情av网站| 啦啦啦免费观看视频1| 黑丝袜美女国产一区| 操出白浆在线播放| 啦啦啦中文免费视频观看日本| 久久久精品区二区三区| 欧美激情 高清一区二区三区| 汤姆久久久久久久影院中文字幕| 淫妇啪啪啪对白视频 | 欧美在线一区亚洲| 多毛熟女@视频| 水蜜桃什么品种好| 久久亚洲精品不卡| a在线观看视频网站| av天堂久久9| 一级a爱视频在线免费观看| 日韩电影二区| 电影成人av| 亚洲国产中文字幕在线视频| 岛国在线观看网站| 国产xxxxx性猛交| 国产亚洲精品久久久久5区| 日韩一卡2卡3卡4卡2021年| 国产1区2区3区精品| 激情视频va一区二区三区| 国产又色又爽无遮挡免| 免费人妻精品一区二区三区视频| 精品一区在线观看国产| 丝袜在线中文字幕| 捣出白浆h1v1| 99香蕉大伊视频| 人人澡人人妻人| 黄网站色视频无遮挡免费观看| 国产精品久久久人人做人人爽| 电影成人av| 高清在线国产一区| 亚洲国产欧美日韩在线播放| 国产精品.久久久| 中文字幕人妻熟女乱码| 国产在视频线精品| 午夜免费观看性视频| 国产精品 国内视频| √禁漫天堂资源中文www| a 毛片基地| 少妇被粗大的猛进出69影院| 超色免费av| 久久国产精品人妻蜜桃| 久久久久久免费高清国产稀缺| 亚洲av美国av| 中文字幕人妻熟女乱码| 午夜福利影视在线免费观看| 久久国产精品男人的天堂亚洲| 免费女性裸体啪啪无遮挡网站| 色综合欧美亚洲国产小说| 老司机亚洲免费影院| 亚洲国产毛片av蜜桃av| 91av网站免费观看| 美女扒开内裤让男人捅视频| 韩国高清视频一区二区三区| 亚洲成人国产一区在线观看| 在线观看免费午夜福利视频| 亚洲专区字幕在线| 日韩 欧美 亚洲 中文字幕| 俄罗斯特黄特色一大片| 亚洲精品自拍成人| 日韩欧美免费精品| 美女高潮到喷水免费观看| 青春草视频在线免费观看| 国产淫语在线视频| 国产精品国产av在线观看| 啦啦啦 在线观看视频| 俄罗斯特黄特色一大片| 超色免费av| 爱豆传媒免费全集在线观看| 多毛熟女@视频| 丁香六月欧美| 91大片在线观看| 色播在线永久视频| 老熟妇乱子伦视频在线观看 | 国产一区二区激情短视频 | 99久久人妻综合| 亚洲性夜色夜夜综合| 亚洲第一欧美日韩一区二区三区 | 99热全是精品| 国产精品一区二区在线观看99| 性少妇av在线| 中文字幕精品免费在线观看视频| 亚洲第一青青草原| 成人黄色视频免费在线看| 精品少妇内射三级| 久久久欧美国产精品| 日本五十路高清| 久久精品国产a三级三级三级| 中亚洲国语对白在线视频| 国产精品 国内视频| 电影成人av| 日本欧美视频一区| 在线观看舔阴道视频| 91大片在线观看| 国产av国产精品国产| 久久精品亚洲熟妇少妇任你| 黄色 视频免费看| 久久久精品免费免费高清| 成人亚洲精品一区在线观看| 亚洲国产欧美一区二区综合| 欧美人与性动交α欧美精品济南到| 欧美少妇被猛烈插入视频| 咕卡用的链子| 在线观看www视频免费| 人人妻,人人澡人人爽秒播| 久久久欧美国产精品| 肉色欧美久久久久久久蜜桃| 99香蕉大伊视频| 可以免费在线观看a视频的电影网站| 国产精品一区二区在线不卡| 丰满少妇做爰视频| 午夜精品国产一区二区电影| 日韩欧美一区二区三区在线观看 | 母亲3免费完整高清在线观看| 国产男人的电影天堂91| 国产亚洲精品一区二区www | a在线观看视频网站| 久久久国产一区二区| 王馨瑶露胸无遮挡在线观看| 激情视频va一区二区三区| 亚洲第一青青草原| 1024视频免费在线观看| 亚洲精品国产一区二区精华液| 亚洲国产欧美网| 国产视频一区二区在线看| 亚洲国产看品久久| 亚洲国产精品一区二区三区在线| 一区二区av电影网| 一级毛片精品| 欧美成人午夜精品| 少妇被粗大的猛进出69影院| 久久免费观看电影| 欧美精品av麻豆av| 老司机深夜福利视频在线观看 | 在线av久久热| 亚洲欧美一区二区三区久久| 日韩欧美国产一区二区入口| 国产欧美亚洲国产| 99热国产这里只有精品6| 亚洲精品久久午夜乱码| 又紧又爽又黄一区二区| 美女中出高潮动态图| 狠狠婷婷综合久久久久久88av| 首页视频小说图片口味搜索| 国产在视频线精品| 天天影视国产精品| 国产伦人伦偷精品视频| 十分钟在线观看高清视频www| 无限看片的www在线观看| 久久精品国产a三级三级三级| 性色av乱码一区二区三区2| 99热网站在线观看| 欧美日韩国产mv在线观看视频| 麻豆乱淫一区二区| 国产91精品成人一区二区三区 | 啦啦啦 在线观看视频| av有码第一页| 免费人妻精品一区二区三区视频| 十分钟在线观看高清视频www| 热99国产精品久久久久久7| 伊人久久大香线蕉亚洲五| 深夜精品福利| 久久免费观看电影| 午夜久久久在线观看| 女性被躁到高潮视频| 波多野结衣av一区二区av| 大香蕉久久网| 精品熟女少妇八av免费久了| 日日摸夜夜添夜夜添小说| 999精品在线视频| 欧美激情 高清一区二区三区| 久久天躁狠狠躁夜夜2o2o| 美女扒开内裤让男人捅视频| 精品少妇内射三级| 国产伦理片在线播放av一区| 国产区一区二久久| 超碰97精品在线观看| 亚洲七黄色美女视频| 精品国内亚洲2022精品成人 | 人妻人人澡人人爽人人| www.av在线官网国产| 午夜激情久久久久久久| 午夜两性在线视频| 亚洲国产精品一区二区三区在线| 91精品伊人久久大香线蕉| avwww免费| 免费在线观看黄色视频的| 老司机影院成人| 一二三四社区在线视频社区8| 国产精品 国内视频| 精品高清国产在线一区| 久久精品国产亚洲av香蕉五月 | 最近最新中文字幕大全免费视频| 真人做人爱边吃奶动态| 啦啦啦免费观看视频1| 男人操女人黄网站| 黑人操中国人逼视频| 一区二区日韩欧美中文字幕| 亚洲欧美成人综合另类久久久| a级毛片黄视频| 亚洲精品国产精品久久久不卡| 亚洲成国产人片在线观看| 日韩一区二区三区影片| 国产欧美日韩精品亚洲av| 精品国产乱码久久久久久小说| 人人澡人人妻人| 国产在线观看jvid| 中亚洲国语对白在线视频| av一本久久久久| 新久久久久国产一级毛片| 欧美精品一区二区大全| 两个人看的免费小视频| 亚洲精品第二区| 性色av一级| 深夜精品福利| 热99re8久久精品国产| 无限看片的www在线观看| av欧美777| 国产伦人伦偷精品视频| 久久久国产一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 一级黄色大片毛片| 在线观看一区二区三区激情| 国产免费福利视频在线观看| 老司机午夜福利在线观看视频 | 久久久久视频综合| 亚洲va日本ⅴa欧美va伊人久久 | 中文字幕另类日韩欧美亚洲嫩草| 国产精品亚洲av一区麻豆| 91九色精品人成在线观看| 国产亚洲av片在线观看秒播厂| 男女下面插进去视频免费观看| 91老司机精品| 久久亚洲国产成人精品v| 成人国语在线视频| 97在线人人人人妻| 精品久久久久久电影网| 欧美+亚洲+日韩+国产| 汤姆久久久久久久影院中文字幕| 国产精品1区2区在线观看. | 在线观看一区二区三区激情| 成人国产av品久久久| 国产精品久久久av美女十八| 热99re8久久精品国产| 国产精品久久久久久精品古装| 欧美精品一区二区大全| 女性被躁到高潮视频| 午夜福利在线免费观看网站| 久久久久久久精品精品| 国产精品二区激情视频| 久久午夜综合久久蜜桃| 亚洲精品久久成人aⅴ小说| 少妇人妻久久综合中文| 欧美亚洲日本最大视频资源| 视频在线观看一区二区三区| 黑人巨大精品欧美一区二区mp4| 天天操日日干夜夜撸| 一本综合久久免费| 国产精品av久久久久免费| 国产成人a∨麻豆精品| 在线天堂中文资源库| 久久人人爽av亚洲精品天堂| kizo精华| 成人国语在线视频| 丝袜在线中文字幕| 51午夜福利影视在线观看| 欧美另类一区| 老熟妇乱子伦视频在线观看 | videosex国产| 麻豆av在线久日| 久久这里只有精品19| 欧美日韩av久久| 性高湖久久久久久久久免费观看| 亚洲伊人久久精品综合| 色老头精品视频在线观看| 亚洲精品美女久久av网站| 国产男人的电影天堂91| 悠悠久久av| 亚洲av日韩在线播放| 人人澡人人妻人| 亚洲第一av免费看| 在线观看免费日韩欧美大片| 无限看片的www在线观看| 十八禁网站免费在线| 一本—道久久a久久精品蜜桃钙片| 丝袜脚勾引网站| 视频区欧美日本亚洲| 黑人猛操日本美女一级片| 丝瓜视频免费看黄片| 亚洲精品自拍成人| 欧美乱码精品一区二区三区| 精品亚洲成国产av| 亚洲精品乱久久久久久| 午夜福利视频在线观看免费| 亚洲精品乱久久久久久| 男女无遮挡免费网站观看| av网站在线播放免费| tube8黄色片| 爱豆传媒免费全集在线观看| 亚洲欧美色中文字幕在线| 成人国产av品久久久| 一本—道久久a久久精品蜜桃钙片| 亚洲av成人一区二区三| 一区福利在线观看| 黄色片一级片一级黄色片| 国产精品99久久99久久久不卡| 午夜免费成人在线视频| 岛国毛片在线播放| 777米奇影视久久| 精品免费久久久久久久清纯 | 久久精品人人爽人人爽视色| 欧美变态另类bdsm刘玥| 国产亚洲欧美在线一区二区| 成年美女黄网站色视频大全免费| 国产麻豆69| 久久青草综合色| 狠狠精品人妻久久久久久综合| 国产精品九九99| 欧美久久黑人一区二区| 一本大道久久a久久精品| 国产精品免费视频内射| 韩国精品一区二区三区| 青春草视频在线免费观看| 国产一区二区激情短视频 | 国产亚洲精品久久久久5区| 亚洲精品久久久久久婷婷小说| 狂野欧美激情性xxxx| 国产成人av激情在线播放| 久久综合国产亚洲精品| 精品福利观看| 亚洲精品国产av蜜桃| 在线永久观看黄色视频| cao死你这个sao货| 爱豆传媒免费全集在线观看| 国产成人a∨麻豆精品| 岛国在线观看网站| 精品亚洲成a人片在线观看| 国产成人免费无遮挡视频| 精品国产一区二区三区四区第35| 满18在线观看网站| 伦理电影免费视频| 日本91视频免费播放| 岛国毛片在线播放| 少妇的丰满在线观看| 色婷婷久久久亚洲欧美| 久久精品人人爽人人爽视色| 亚洲欧美日韩高清在线视频 | 亚洲专区字幕在线| 国产精品久久久久久人妻精品电影 | 亚洲欧美成人综合另类久久久| 日韩 亚洲 欧美在线| 免费人妻精品一区二区三区视频| 9热在线视频观看99| 亚洲一区二区三区欧美精品| 欧美 亚洲 国产 日韩一| 国产淫语在线视频| 国产亚洲午夜精品一区二区久久| 涩涩av久久男人的天堂| www.精华液| 国产亚洲av片在线观看秒播厂| 欧美精品一区二区免费开放| 亚洲av成人不卡在线观看播放网 | 老汉色∧v一级毛片| 亚洲九九香蕉| 精品一区二区三区av网在线观看 | 热99国产精品久久久久久7| 国产深夜福利视频在线观看| 国产精品av久久久久免费| 老司机福利观看| 亚洲av成人一区二区三| 国产精品久久久av美女十八| 伊人亚洲综合成人网| 午夜成年电影在线免费观看| 秋霞在线观看毛片| 久久人妻熟女aⅴ| 欧美国产精品va在线观看不卡| 在线观看www视频免费| 女性生殖器流出的白浆| 国产成人影院久久av| 搡老熟女国产l中国老女人| 免费少妇av软件| 精品一区二区三区四区五区乱码| 亚洲一区中文字幕在线| 久久国产精品大桥未久av| www.精华液| 1024香蕉在线观看| 爱豆传媒免费全集在线观看| 考比视频在线观看| 日韩中文字幕视频在线看片| 18禁国产床啪视频网站| 日韩中文字幕视频在线看片| 宅男免费午夜| 色婷婷av一区二区三区视频| 男人添女人高潮全过程视频| 国产麻豆69| 日韩视频在线欧美| 精品国产一区二区久久| 久久精品熟女亚洲av麻豆精品| 日韩大片免费观看网站| 在线永久观看黄色视频| 亚洲国产中文字幕在线视频| 男女床上黄色一级片免费看| 99热国产这里只有精品6| 亚洲一区中文字幕在线| 中文字幕人妻丝袜制服| 男男h啪啪无遮挡| 最近最新中文字幕大全免费视频| 极品人妻少妇av视频| 悠悠久久av| 欧美黄色淫秽网站| 麻豆乱淫一区二区| 这个男人来自地球电影免费观看| 建设人人有责人人尽责人人享有的| 人人妻人人澡人人看| 侵犯人妻中文字幕一二三四区| 精品国产一区二区三区四区第35| 丰满迷人的少妇在线观看| 亚洲国产看品久久| 中国国产av一级| 日韩欧美一区二区三区在线观看 | 丝袜喷水一区| 亚洲欧洲精品一区二区精品久久久| 久久久久精品人妻al黑| 一区在线观看完整版| 蜜桃国产av成人99| 欧美日韩成人在线一区二区| 人人妻,人人澡人人爽秒播| www.精华液| 99久久人妻综合| 亚洲国产中文字幕在线视频| 成年人免费黄色播放视频| 人人澡人人妻人| 国产av精品麻豆| 国产成人免费无遮挡视频| 91成年电影在线观看| 国产亚洲一区二区精品| 美女中出高潮动态图| 婷婷丁香在线五月| 亚洲精品一卡2卡三卡4卡5卡 | 他把我摸到了高潮在线观看 | 午夜福利在线免费观看网站| 国精品久久久久久国模美| 久久国产精品男人的天堂亚洲| 性少妇av在线| 午夜两性在线视频| 免费高清在线观看视频在线观看| 国产成人a∨麻豆精品| 亚洲精品国产av成人精品| av视频免费观看在线观看| 两个人免费观看高清视频| 午夜精品久久久久久毛片777| 国产精品熟女久久久久浪| 日韩欧美免费精品| 我的亚洲天堂| 精品人妻在线不人妻| 欧美日韩福利视频一区二区| 免费日韩欧美在线观看| 国产精品自产拍在线观看55亚洲 | 美女福利国产在线| 老司机在亚洲福利影院| 老汉色av国产亚洲站长工具| 美女大奶头黄色视频| 777久久人妻少妇嫩草av网站| 巨乳人妻的诱惑在线观看| 中亚洲国语对白在线视频| 亚洲精品国产一区二区精华液| 女人爽到高潮嗷嗷叫在线视频| 一个人免费看片子| 国产国语露脸激情在线看| 久久青草综合色| 免费黄频网站在线观看国产| 日韩视频一区二区在线观看| 精品欧美一区二区三区在线| 婷婷成人精品国产| 欧美日韩亚洲高清精品| 成人黄色视频免费在线看| 国产区一区二久久| 下体分泌物呈黄色| 精品少妇黑人巨大在线播放| 青春草亚洲视频在线观看| 日本五十路高清| 亚洲欧美成人综合另类久久久| 精品福利观看| 亚洲一码二码三码区别大吗| 一本综合久久免费| 欧美日韩视频精品一区| 少妇粗大呻吟视频| 中文字幕人妻熟女乱码| 九色亚洲精品在线播放| 最近最新免费中文字幕在线| 动漫黄色视频在线观看| 久久九九热精品免费| 国产精品免费视频内射| 99热网站在线观看| 99久久人妻综合| 伦理电影免费视频| 一本一本久久a久久精品综合妖精| 日日爽夜夜爽网站| 欧美激情极品国产一区二区三区| av网站在线播放免费| 在线天堂中文资源库| 亚洲中文日韩欧美视频| 十八禁人妻一区二区| 亚洲三区欧美一区| 精品免费久久久久久久清纯 | 久久午夜综合久久蜜桃| 叶爱在线成人免费视频播放| 少妇 在线观看| 如日韩欧美国产精品一区二区三区| 欧美97在线视频| 精品高清国产在线一区| 精品久久久久久久毛片微露脸 | 少妇精品久久久久久久| 日本av手机在线免费观看| 伊人亚洲综合成人网| 久热爱精品视频在线9| 国产精品偷伦视频观看了| 我要看黄色一级片免费的| videosex国产| 精品国内亚洲2022精品成人 | 国精品久久久久久国模美| 亚洲一码二码三码区别大吗| 免费日韩欧美在线观看| 69精品国产乱码久久久| 国产主播在线观看一区二区| 男女高潮啪啪啪动态图| 肉色欧美久久久久久久蜜桃| 国产亚洲一区二区精品| 久久性视频一级片|