• <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)的影響
    欧美老熟妇乱子伦牲交| 丰满少妇做爰视频| 午夜免费男女啪啪视频观看| 午夜老司机福利片| 亚洲熟女毛片儿| 一级片免费观看大全| 我要看黄色一级片免费的| 成人亚洲精品一区在线观看| 黑人猛操日本美女一级片| h视频一区二区三区| 女人精品久久久久毛片| 精品国产超薄肉色丝袜足j| 视频在线观看一区二区三区| 水蜜桃什么品种好| 老鸭窝网址在线观看| 亚洲国产精品成人久久小说| 性色av乱码一区二区三区2| 91九色精品人成在线观看| videosex国产| 9色porny在线观看| 尾随美女入室| 搡老乐熟女国产| 99久久精品国产亚洲精品| www.熟女人妻精品国产| 国产成人精品久久二区二区免费| 91九色精品人成在线观看| 成在线人永久免费视频| 亚洲色图综合在线观看| 19禁男女啪啪无遮挡网站| 亚洲,一卡二卡三卡| 成年人免费黄色播放视频| xxxhd国产人妻xxx| 国产精品.久久久| 精品视频人人做人人爽| 日韩av免费高清视频| 亚洲精品久久成人aⅴ小说| 午夜久久久在线观看| 赤兔流量卡办理| 高清av免费在线| 久久毛片免费看一区二区三区| 国产精品一区二区在线不卡| 免费看不卡的av| 人人澡人人妻人| www.av在线官网国产| 午夜日韩欧美国产| 美女大奶头黄色视频| 久久久亚洲精品成人影院| 亚洲,欧美,日韩| 丁香六月天网| 水蜜桃什么品种好| 亚洲精品av麻豆狂野| 首页视频小说图片口味搜索 | 老司机靠b影院| 日韩制服骚丝袜av| 免费少妇av软件| 国产高清国产精品国产三级| 国产精品二区激情视频| 女人爽到高潮嗷嗷叫在线视频| 一边摸一边做爽爽视频免费| 国产日韩欧美亚洲二区| 一本一本久久a久久精品综合妖精| www.av在线官网国产| 国产av精品麻豆| 最新的欧美精品一区二区| 最新在线观看一区二区三区 | 亚洲欧美激情在线| 精品一品国产午夜福利视频| 中文字幕人妻丝袜制服| 亚洲人成电影免费在线| 国产成人啪精品午夜网站| 啦啦啦视频在线资源免费观看| av国产精品久久久久影院| 热99国产精品久久久久久7| 国产亚洲一区二区精品| 婷婷色综合大香蕉| 黑丝袜美女国产一区| 日韩av在线免费看完整版不卡| 一本—道久久a久久精品蜜桃钙片| 大片电影免费在线观看免费| 女警被强在线播放| e午夜精品久久久久久久| 欧美日韩亚洲国产一区二区在线观看 | 巨乳人妻的诱惑在线观看| 视频区欧美日本亚洲| 青草久久国产| 丰满人妻熟妇乱又伦精品不卡| 久久精品久久精品一区二区三区| 亚洲中文字幕日韩| 国产xxxxx性猛交| 午夜福利在线免费观看网站| 亚洲精品美女久久久久99蜜臀 | 手机成人av网站| 亚洲伊人色综图| 日韩,欧美,国产一区二区三区| 汤姆久久久久久久影院中文字幕| 香蕉国产在线看| 狂野欧美激情性xxxx| 国产成人免费无遮挡视频| 亚洲五月色婷婷综合| 欧美97在线视频| 午夜福利在线免费观看网站| 后天国语完整版免费观看| 久久ye,这里只有精品| 亚洲国产精品成人久久小说| 一级a爱视频在线免费观看| 建设人人有责人人尽责人人享有的| 亚洲精品国产av蜜桃| 在线观看免费视频网站a站| 国产欧美日韩一区二区三 | 嫩草影视91久久| 少妇精品久久久久久久| 婷婷色综合大香蕉| 日本色播在线视频| 大片电影免费在线观看免费| 丝袜美足系列| 精品久久久精品久久久| 丝袜脚勾引网站| 下体分泌物呈黄色| av线在线观看网站| 免费看不卡的av| 亚洲,欧美精品.| 亚洲人成77777在线视频| 日韩一本色道免费dvd| 日韩中文字幕视频在线看片| 亚洲av男天堂| 日韩中文字幕视频在线看片| 免费av中文字幕在线| 丁香六月欧美| 国产精品三级大全| 国产片内射在线| 欧美人与性动交α欧美软件| 超碰97精品在线观看| 建设人人有责人人尽责人人享有的| 美女午夜性视频免费| 在线观看人妻少妇| 丝瓜视频免费看黄片| 久久天堂一区二区三区四区| 1024香蕉在线观看| 日韩一本色道免费dvd| 日韩一卡2卡3卡4卡2021年| 欧美日韩亚洲国产一区二区在线观看 | 国产欧美日韩一区二区三区在线| 国产成人系列免费观看| 久久久久精品人妻al黑| 久久99一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 欧美精品av麻豆av| av不卡在线播放| 久久久久久久国产电影| 蜜桃在线观看..| 免费在线观看完整版高清| 久久久精品免费免费高清| 女性生殖器流出的白浆| 丝袜人妻中文字幕| 七月丁香在线播放| 欧美日韩视频高清一区二区三区二| www日本在线高清视频| 丝袜在线中文字幕| 国产精品国产av在线观看| 国产一区二区激情短视频 | 久久九九热精品免费| 2018国产大陆天天弄谢| 男男h啪啪无遮挡| 精品少妇一区二区三区视频日本电影| 国产视频一区二区在线看| 涩涩av久久男人的天堂| 欧美日本中文国产一区发布| 婷婷色麻豆天堂久久| 精品欧美一区二区三区在线| 丁香六月欧美| 国产无遮挡羞羞视频在线观看| 久久狼人影院| 中文字幕另类日韩欧美亚洲嫩草| 99热国产这里只有精品6| 久久热在线av| 欧美日韩av久久| 涩涩av久久男人的天堂| 性高湖久久久久久久久免费观看| 欧美成狂野欧美在线观看| 久久久精品94久久精品| av不卡在线播放| 涩涩av久久男人的天堂| 1024视频免费在线观看| 精品高清国产在线一区| 在线观看免费视频网站a站| 成年人黄色毛片网站| 久久精品久久精品一区二区三区| 亚洲第一av免费看| 国产精品偷伦视频观看了| 国产av国产精品国产| 一个人免费看片子| 亚洲天堂av无毛| 美国免费a级毛片| 午夜免费观看性视频| 91国产中文字幕| 男男h啪啪无遮挡| 国产真人三级小视频在线观看| 汤姆久久久久久久影院中文字幕| 亚洲国产欧美网| 日韩电影二区| 欧美日韩一级在线毛片| a 毛片基地| 精品少妇黑人巨大在线播放| 国产精品人妻久久久影院| 国产精品国产三级国产专区5o| 成人亚洲精品一区在线观看| 精品一区在线观看国产| 人人妻,人人澡人人爽秒播 | 最近手机中文字幕大全| 亚洲精品一二三| 视频区图区小说| 国产老妇伦熟女老妇高清| 亚洲av美国av| 婷婷丁香在线五月| 久久久精品区二区三区| 国产亚洲精品久久久久5区| 国产日韩一区二区三区精品不卡| 国产亚洲精品第一综合不卡| 国产精品国产三级专区第一集| 色播在线永久视频| 久久av网站| 精品一区二区三区四区五区乱码 | 波野结衣二区三区在线| 国产日韩一区二区三区精品不卡| 国产精品麻豆人妻色哟哟久久| 超色免费av| 午夜91福利影院| 久久久精品免费免费高清| 日本一区二区免费在线视频| 亚洲精品第二区| 久久女婷五月综合色啪小说| 黑人欧美特级aaaaaa片| 搡老岳熟女国产| 后天国语完整版免费观看| 91老司机精品| 一区二区日韩欧美中文字幕| 操美女的视频在线观看| 日韩中文字幕欧美一区二区 | 亚洲一区中文字幕在线| www.999成人在线观看| 亚洲精品乱久久久久久| 9热在线视频观看99| 国产主播在线观看一区二区 | 日本五十路高清| 久久精品成人免费网站| 好男人视频免费观看在线| 久久久欧美国产精品| av片东京热男人的天堂| 免费在线观看影片大全网站 | 天天躁日日躁夜夜躁夜夜| 18禁国产床啪视频网站| 在线av久久热| 亚洲国产成人一精品久久久| 欧美在线一区亚洲| 亚洲欧美一区二区三区久久| 成年人免费黄色播放视频| 亚洲国产中文字幕在线视频| 国产在线一区二区三区精| 汤姆久久久久久久影院中文字幕| 女性被躁到高潮视频| 九草在线视频观看| 男女边吃奶边做爰视频| 丝袜喷水一区| 精品一区二区三区四区五区乱码 | 大型av网站在线播放| 日日摸夜夜添夜夜爱| 精品一品国产午夜福利视频| 成人国产av品久久久| 久久热在线av| 亚洲精品久久成人aⅴ小说| 久久这里只有精品19| 在线观看www视频免费| 狠狠婷婷综合久久久久久88av| 中国国产av一级| 国产成人系列免费观看| 极品人妻少妇av视频| 热re99久久精品国产66热6| 你懂的网址亚洲精品在线观看| 亚洲精品久久午夜乱码| 久久精品aⅴ一区二区三区四区| 久久这里只有精品19| 性色av一级| 欧美在线一区亚洲| 2021少妇久久久久久久久久久| 啦啦啦 在线观看视频| 男女高潮啪啪啪动态图| 国产亚洲精品第一综合不卡| 少妇人妻 视频| 国产视频一区二区在线看| 亚洲国产精品一区二区三区在线| 欧美日韩国产mv在线观看视频| 欧美黄色片欧美黄色片| 夫妻午夜视频| 国产成人系列免费观看| 一本色道久久久久久精品综合| 在线亚洲精品国产二区图片欧美| 在线观看免费日韩欧美大片| 日韩 亚洲 欧美在线| 国产有黄有色有爽视频| av国产久精品久网站免费入址| 十八禁网站网址无遮挡| 精品第一国产精品| av视频免费观看在线观看| 国产97色在线日韩免费| 97在线人人人人妻| 国产精品秋霞免费鲁丝片| 国产亚洲一区二区精品| 超色免费av| 一边摸一边抽搐一进一出视频| 亚洲av电影在线进入| 色视频在线一区二区三区| 午夜免费男女啪啪视频观看| 婷婷丁香在线五月| 久久久国产欧美日韩av| 亚洲欧洲国产日韩| 每晚都被弄得嗷嗷叫到高潮| 男人操女人黄网站| 国产日韩欧美视频二区| 亚洲av国产av综合av卡| 国产精品一区二区精品视频观看| 国产伦人伦偷精品视频| 亚洲国产最新在线播放| 国产成人精品在线电影| 热re99久久国产66热| 99re6热这里在线精品视频| 国产精品 欧美亚洲| 丝袜喷水一区| 国产一区亚洲一区在线观看| 波多野结衣av一区二区av| 一区二区av电影网| 久久人人97超碰香蕉20202| 亚洲美女黄色视频免费看| 国产黄色视频一区二区在线观看| 美女大奶头黄色视频| 亚洲成人手机| 国产成人精品在线电影| 王馨瑶露胸无遮挡在线观看| 欧美老熟妇乱子伦牲交| 久久ye,这里只有精品| 精品一品国产午夜福利视频| 啦啦啦在线观看免费高清www| 国产日韩欧美亚洲二区| 十分钟在线观看高清视频www| 好男人电影高清在线观看| 免费观看a级毛片全部| 久久久精品区二区三区| 欧美成狂野欧美在线观看| 秋霞在线观看毛片| 国产黄色免费在线视频| 97精品久久久久久久久久精品| 亚洲欧美一区二区三区黑人| 色婷婷久久久亚洲欧美| 伊人亚洲综合成人网| 又黄又粗又硬又大视频| 免费观看人在逋| 一本久久精品| 狠狠婷婷综合久久久久久88av| 91九色精品人成在线观看| 看免费成人av毛片| 啦啦啦视频在线资源免费观看| 亚洲久久久国产精品| 两个人看的免费小视频| 亚洲,欧美,日韩| 国产人伦9x9x在线观看| 国产不卡av网站在线观看| 一级毛片电影观看| 亚洲精品久久成人aⅴ小说| 国产亚洲欧美精品永久| 久久狼人影院| 国产视频首页在线观看| 亚洲视频免费观看视频| 国产片内射在线| 人妻 亚洲 视频| 高清av免费在线| 日韩 欧美 亚洲 中文字幕| 操美女的视频在线观看| 国产在线观看jvid| 一边亲一边摸免费视频| 日韩欧美一区视频在线观看| 一本—道久久a久久精品蜜桃钙片| av在线app专区| 99久久精品国产亚洲精品| 黑人巨大精品欧美一区二区蜜桃| 少妇人妻久久综合中文| 无遮挡黄片免费观看| 菩萨蛮人人尽说江南好唐韦庄| 欧美+亚洲+日韩+国产| 久久这里只有精品19| 亚洲av成人不卡在线观看播放网 | 日日摸夜夜添夜夜爱| 国产三级黄色录像| 精品久久久精品久久久| 满18在线观看网站| 老司机亚洲免费影院| 亚洲欧美一区二区三区黑人| 丝袜脚勾引网站| 一级,二级,三级黄色视频| 久久久久国产精品人妻一区二区| 亚洲欧美一区二区三区黑人| 老司机在亚洲福利影院| 亚洲激情五月婷婷啪啪| 国产视频一区二区在线看| 亚洲人成电影观看| 久久久久久久大尺度免费视频| 如日韩欧美国产精品一区二区三区| 久久久久久免费高清国产稀缺| 国产精品成人在线| 女人被躁到高潮嗷嗷叫费观| 国产极品粉嫩免费观看在线| 精品亚洲乱码少妇综合久久| 成年av动漫网址| 一区二区av电影网| 国产精品麻豆人妻色哟哟久久| xxx大片免费视频| 免费少妇av软件| 午夜激情久久久久久久| 国产主播在线观看一区二区 | 国产欧美亚洲国产| 亚洲成人手机| 精品国产一区二区三区久久久樱花| 亚洲自偷自拍图片 自拍| 国产成人精品久久二区二区91| 国产片特级美女逼逼视频| 亚洲av欧美aⅴ国产| 一级黄片播放器| 又大又爽又粗| 妹子高潮喷水视频| 婷婷色综合大香蕉| 精品亚洲成国产av| 午夜久久久在线观看| 久热这里只有精品99| 国产视频首页在线观看| 国产欧美日韩精品亚洲av| 国产成人免费无遮挡视频| 大香蕉久久网| 高清欧美精品videossex| 日本a在线网址| 女性被躁到高潮视频| 亚洲国产欧美一区二区综合| 99热全是精品| 欧美精品一区二区免费开放| 伊人久久大香线蕉亚洲五| 后天国语完整版免费观看| 老司机靠b影院| 国产免费福利视频在线观看| 男人爽女人下面视频在线观看| 黑丝袜美女国产一区| 热99国产精品久久久久久7| 久久99一区二区三区| 老司机影院成人| 伊人亚洲综合成人网| 一边摸一边做爽爽视频免费| 亚洲av国产av综合av卡| 久久久久久久大尺度免费视频| 三上悠亚av全集在线观看| 国产av一区二区精品久久| 日本猛色少妇xxxxx猛交久久| 久久精品aⅴ一区二区三区四区| 桃花免费在线播放| 欧美精品高潮呻吟av久久| 操美女的视频在线观看| 久久狼人影院| xxx大片免费视频| 两个人免费观看高清视频| 在线观看国产h片| 国产精品二区激情视频| 91精品伊人久久大香线蕉| 黑丝袜美女国产一区| 2021少妇久久久久久久久久久| 电影成人av| 亚洲欧美激情在线| 欧美日韩一级在线毛片| 极品少妇高潮喷水抽搐| 日本猛色少妇xxxxx猛交久久| 久久毛片免费看一区二区三区| 欧美大码av| 日本黄色日本黄色录像| 国精品久久久久久国模美| 老司机影院毛片| 老司机深夜福利视频在线观看 | 日本av免费视频播放| 成年av动漫网址| 日韩伦理黄色片| 精品人妻一区二区三区麻豆| 又粗又硬又长又爽又黄的视频| 在线观看免费视频网站a站| 黑人欧美特级aaaaaa片| 人人妻人人澡人人爽人人夜夜| 久久精品国产综合久久久| 99国产综合亚洲精品| 老司机深夜福利视频在线观看 | 少妇人妻久久综合中文| 男人添女人高潮全过程视频| 午夜久久久在线观看| 亚洲视频免费观看视频| 美女扒开内裤让男人捅视频| 18禁国产床啪视频网站| 真人做人爱边吃奶动态| 国产伦理片在线播放av一区| 国产亚洲午夜精品一区二区久久| 99久久人妻综合| 国产精品98久久久久久宅男小说| 在线观看一区二区三区| 在线观看午夜福利视频| 亚洲avbb在线观看| 伊人久久大香线蕉亚洲五| 两性夫妻黄色片| 熟女电影av网| 一边摸一边抽搐一进一小说| 国产成人欧美在线观看| ponron亚洲| 久久人妻福利社区极品人妻图片| 可以在线观看毛片的网站| 欧美av亚洲av综合av国产av| 99热6这里只有精品| 亚洲av片天天在线观看| 国产精品,欧美在线| 免费看日本二区| 色哟哟哟哟哟哟| 久久人妻福利社区极品人妻图片| 黑人欧美特级aaaaaa片| 精品国产美女av久久久久小说| 老司机福利观看| 日韩成人在线观看一区二区三区| 人成视频在线观看免费观看| 欧美成人免费av一区二区三区| 悠悠久久av| 精品国产乱码久久久久久男人| 欧美大码av| 国产精品,欧美在线| 国产av一区在线观看免费| 男女做爰动态图高潮gif福利片| 国产激情久久老熟女| 婷婷六月久久综合丁香| 日本三级黄在线观看| 一区福利在线观看| 精品久久久久久成人av| 欧美在线黄色| 色综合欧美亚洲国产小说| 国产不卡一卡二| 老司机午夜福利在线观看视频| 久久国产精品男人的天堂亚洲| 午夜福利成人在线免费观看| avwww免费| 亚洲国产精品久久男人天堂| 国产精品久久久久久精品电影 | 日本熟妇午夜| 性色av乱码一区二区三区2| 一本大道久久a久久精品| 黑丝袜美女国产一区| 国产伦在线观看视频一区| 午夜亚洲福利在线播放| 美国免费a级毛片| 成人午夜高清在线视频 | 欧美成人性av电影在线观看| 一本大道久久a久久精品| 久久久国产精品麻豆| 国产精品野战在线观看| 欧美午夜高清在线| 国产av一区在线观看免费| 国产一卡二卡三卡精品| 国产极品粉嫩免费观看在线| 精品国产亚洲在线| 久久午夜亚洲精品久久| 黄色成人免费大全| 嫁个100分男人电影在线观看| 国产精品99久久99久久久不卡| 亚洲aⅴ乱码一区二区在线播放 | 欧美激情高清一区二区三区| 无限看片的www在线观看| av视频在线观看入口| 亚洲自偷自拍图片 自拍| 国产亚洲精品久久久久5区| 天天一区二区日本电影三级| 69av精品久久久久久| 色婷婷久久久亚洲欧美| av超薄肉色丝袜交足视频| av免费在线观看网站| 日本 欧美在线| 欧美三级亚洲精品| 成人手机av| ponron亚洲| 在线观看免费午夜福利视频| 97碰自拍视频| 免费在线观看成人毛片| 亚洲国产精品999在线| 久久草成人影院| 午夜福利一区二区在线看| 久久伊人香网站| 美女免费视频网站| 成人三级做爰电影| 1024香蕉在线观看| 欧美久久黑人一区二区| 久久久久久国产a免费观看| 亚洲av片天天在线观看| 国产精品自产拍在线观看55亚洲| 国产伦在线观看视频一区| а√天堂www在线а√下载| 婷婷亚洲欧美| 国产精品 国内视频| 国产精品免费视频内射| 国产精品久久视频播放| 欧美精品亚洲一区二区| 欧美黑人巨大hd| 欧美精品啪啪一区二区三区| 黑丝袜美女国产一区| 一本综合久久免费| 国产黄色小视频在线观看| 婷婷亚洲欧美| 少妇裸体淫交视频免费看高清 | www日本黄色视频网| 亚洲一卡2卡3卡4卡5卡精品中文| 熟女电影av网| 男女视频在线观看网站免费 | 亚洲天堂国产精品一区在线|