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

    基于BIM的鐵路車輛-無砟軌道護(hù)軌性能分析

    2021-10-14 05:06:14劉沛陳衛(wèi)國張耀
    鐵道建筑 2021年9期
    關(guān)鍵詞:軌枕輪軌鋼軌

    劉沛 陳衛(wèi)國 張耀

    中鐵二院工程集團(tuán)有限責(zé)任公司,成都610031

    建筑信息模型(Building Information Modeling,BIM)技術(shù)在鐵路工程設(shè)計(jì)領(lǐng)域發(fā)展迅速,各專業(yè)正逐步實(shí)現(xiàn)完全或部分BIM正向設(shè)計(jì)。伴隨設(shè)計(jì)環(huán)境與設(shè)計(jì)內(nèi)容的三維化進(jìn)程,對設(shè)計(jì)方案性能的校核分析工作也應(yīng)實(shí)現(xiàn)空間可視化?;贐IM平臺和BIM模型在三維環(huán)境中直接對設(shè)計(jì)模型開展三維仿真計(jì)算分析是當(dāng)代軌道交通領(lǐng)域數(shù)字化設(shè)計(jì)階段的關(guān)鍵手段之一[1]。

    目前國內(nèi)外多用剛體動力學(xué)理論分析脫軌后無砟軌道的護(hù)軌性能,忽略了車輛與鋼軌的柔性變形和軌道板自身的撓度,導(dǎo)致仿真結(jié)果與實(shí)際運(yùn)動狀態(tài)相差較大[2-3]。從有限元角度出發(fā),將BIM模型離散為全尺寸有限元模型,可以降低仿真模型與實(shí)際狀態(tài)的差異。利用動態(tài)非線性有限元方法解決運(yùn)動和變形問題,可以使列車在脫軌后撞擊軌枕塊及護(hù)軌過程中的傳力過程和運(yùn)動狀態(tài)具有真實(shí)的柔性體本質(zhì),從而提高鐵路工程BIM三維設(shè)計(jì)的效率和準(zhǔn)確性[4-5]。

    本文選取某高速鐵路雙塊式雙線直線無砟軌道護(hù)軌段BIM設(shè)計(jì)模型為研究對象,以某市域列車頭車車輛作為脫軌車輛來分析護(hù)軌的防護(hù)性能,通過創(chuàng)建基于BIM模型的無砟軌道-車輛有限元脫軌分析模型,設(shè)置基于歐洲標(biāo)準(zhǔn)EN 15227—2020《Railway Applications—Crashworthiness Requirements for Rail Vehicles》[4]和國內(nèi)外列車脫軌事故[5]的脫軌場景和邊界條件,利用LS?DYNA軟件對該護(hù)軌結(jié)構(gòu)的防護(hù)性能和設(shè)置方案進(jìn)行仿真,研究軌道的護(hù)軌結(jié)構(gòu)對車輛的保護(hù)作用和自身的力學(xué)響應(yīng)。

    1 脫軌分析理論

    車輛在脫軌后與軌道結(jié)構(gòu)發(fā)生碰撞,通過鋼軌和車輛自身的變形及兩者之間的摩擦來消耗動能。為了更真實(shí)地模擬車輛與軌道碰撞時的柔性特征,以動態(tài)非線性大變形有限元理論作為脫軌分析的理論基礎(chǔ)。

    1.1 動態(tài)非線性有限元基本方程

    基于宏觀的運(yùn)動學(xué)和微觀的彈塑性變形原理,車輛的脫軌與護(hù)軌碰撞過程的任意質(zhì)點(diǎn)須滿足運(yùn)動方程、連續(xù)性方程、動量守恒方程、能量守恒方程以及各類邊界條件(面力邊界、位移邊界、接觸邊界),并最終滿足以下平衡方程[6]:

    式中:Ω為空間域;ρ為質(zhì)量密度;x?為質(zhì)點(diǎn)加速度;σij為應(yīng)力張量,向j求導(dǎo)即為σij,j;f為單位質(zhì)量所受體積力;xi為質(zhì)點(diǎn)在任意時刻ti的位置坐標(biāo),δxi表示位移微量;S1、S0分別為面力邊界、接觸邊界;nj為表面外法線方向單位向量;分別為接觸邊界兩側(cè)對邊界上產(chǎn)生的應(yīng)力。

    1.2 積分算法

    求解運(yùn)動方程多采用顯式中心差分法,其為預(yù)測性質(zhì)的遞推方法。LS?DYNA軟件采用的是一種計(jì)算精度較高的中心差分法的變異表達(dá),求解時位移、速度、加速度在不同的時刻交錯,因此被稱為“蛙跳式中心差分法”[7],其表達(dá)式為

    式中:M為總體質(zhì)量矩陣;P為由節(jié)點(diǎn)載荷、面力和體積力組成的總體節(jié)點(diǎn)載荷向量;F為單元應(yīng)力場集成的等效節(jié)點(diǎn)載荷向量;H為每個單元節(jié)點(diǎn)處的沙漏阻尼力集成向量;C為總體阻尼矩陣;x?為質(zhì)點(diǎn)速度,x為質(zhì)點(diǎn)位移;t的下標(biāo)表示不同的時刻;Δtn、Δtn+1分別為tn、tn+1時刻的時間步長。

    根據(jù)式(2)—式(4)可得到tn+1時刻的位移x(tn+1),加上初始構(gòu)型便可得當(dāng)前構(gòu)型。

    2 無砟軌道-車輛模型

    2.1 無砟軌道模型

    2.1.1 無砟軌道BIM正向設(shè)計(jì)流程(圖1)[8]

    圖1 軌道BIM正向設(shè)計(jì)流程

    設(shè)計(jì)前先依照IFC/IFD編碼標(biāo)準(zhǔn)創(chuàng)建EC Schema文件,用建庫工具創(chuàng)建軌道標(biāo)準(zhǔn)構(gòu)件及模板,整理好鋼軌、軌枕、道床等部件的布置和規(guī)格信息文件。開始軌道設(shè)計(jì)時,通過協(xié)同設(shè)計(jì)平臺參考路線模型,分析曲線段超高,導(dǎo)入鋼軌、軌枕、道床等部件的相關(guān)數(shù)據(jù)文件,交互編輯超高信息,然后生成軌道設(shè)計(jì)方案?;诖嗽O(shè)計(jì)方案,生成軌道三維BIM模型,然后檢查模型參數(shù)的正確性和信息的完備性直至符合設(shè)計(jì)技術(shù)要求。

    2.1.2 創(chuàng)建無砟軌道“BIM+有限元”模型

    選取某高速鐵路雙塊式雙線直線無砟軌道護(hù)軌段為研究對象。該路段全長45 m,左右線線路中心距4.4 m。護(hù)軌與外側(cè)鋼軌安裝在同一軌枕塊上,與運(yùn)行軌間距0.330 m;內(nèi)側(cè)鋼軌不設(shè)置護(hù)軌,軌枕塊縱向間距0.625 m;護(hù)軌與運(yùn)行軌頂面高程相等。護(hù)輪軌道設(shè)計(jì)如圖2所示。

    圖2 護(hù)輪軌道設(shè)計(jì)方案示意(單位:mm)

    根據(jù)上述設(shè)計(jì)方案,利用基于Bentley?OpenRail Designer軟件二次開發(fā)的鐵路軌道三維正向設(shè)計(jì)工具,通過超高分析、計(jì)算軌道配板及布設(shè)參數(shù),設(shè)計(jì)并創(chuàng)建出該路段LOD3.0等級的軌道BIM模型,包括軌道板、軌枕塊、鋼軌、扣件、鋼軌墊板、扣件螺栓等,如圖3所示[9-10]。

    圖3 雙塊式無砟軌道護(hù)軌段BIM模型及主要構(gòu)件細(xì)節(jié)

    根據(jù)有限元理論的網(wǎng)格劃分原則,采用HyperMesh14.0軟件把整個雙線軌道的BIM模型離散成具有真實(shí)的構(gòu)件間連接和接觸關(guān)系的有限元模型[11]。其中,將軌道板、軌枕塊、鋼軌、鋼軌墊板、扣件螺栓離散成8節(jié)點(diǎn)的實(shí)體單元,扣件彈簧用3節(jié)點(diǎn)梁單元模擬,并設(shè)置預(yù)緊力。雙塊式雙線無砟軌道有限元模型共有641 253個單元。

    2.2 車輛模型

    選取某市域列車頭車車輛作為脫軌車輛來分析護(hù)軌的防護(hù)性能。為使脫軌仿真結(jié)果具備更真實(shí)的脫軌形態(tài)和力學(xué)變化過程,須創(chuàng)建準(zhǔn)確而合理的車輛有限元模型,主要包括車輛轉(zhuǎn)向架與車體結(jié)構(gòu)。

    轉(zhuǎn)向架作為車輛的獨(dú)立部件,是脫軌分析過程中傳遞輪軌作用力的唯一部件。利用HyperMesh14.0軟件將車輛主體結(jié)構(gòu)離散為24 561個4節(jié)點(diǎn)薄殼單元和1 897個8節(jié)點(diǎn)實(shí)體單元,一系、二系空氣彈簧分別用36、72個6自由度的離散梁單元模擬,枕梁與車體之間的連接關(guān)系用點(diǎn)焊模擬。轉(zhuǎn)向架有限元模型如圖4所示。根據(jù)LS?DYNA的關(guān)鍵字手冊[12],將模型薄殼單元的積分算法設(shè)置為Belytschko?Tsay法,實(shí)體單元為CSSE法,離散梁單元為discrete beam法。

    圖4 市域列車頭車轉(zhuǎn)向架有限元模型

    該車輛車體為鋁合金結(jié)構(gòu)。首先將車輛三維模型中面化,并確保正確的組件連接關(guān)系,然后將車體三維中面模型離散成由875 903個薄殼單元組成的有限元模型,用Mass21單元搭配Beam單元模擬車內(nèi)安裝的各類設(shè)備,通過節(jié)點(diǎn)耦合的方式將二系彈簧與車體進(jìn)行連接。

    2.3 無砟軌道-車輛模型

    車輛在脫軌時對處于運(yùn)行軌道旁邊的護(hù)軌有一定的沖擊作用,護(hù)軌結(jié)構(gòu)對脫軌列車起保護(hù)作用,綜合護(hù)軌軌道和車輛的有限元模型形成脫軌系統(tǒng)。脫軌系統(tǒng)質(zhì)量為54.165 t,車輪與軌道之間設(shè)置防止互相穿透的面面接觸,車輛自身設(shè)置防止自身構(gòu)件互相穿透的單面接觸。最終建立的脫軌系統(tǒng)有限元分析模型如圖5所示。

    圖5 脫軌系統(tǒng)有限元分析模型

    3 車輛脫軌仿真計(jì)算及結(jié)果分析

    根據(jù)列車的運(yùn)行特征,對列車模型施加橫向和縱向平動初始速度作為動能輸入條件,依據(jù)文獻(xiàn)[4-5]確定計(jì)算工況的邊界條件,見表1。

    表1 車輛脫軌邊界條件

    利用LS?DYNA求解器的顯式求解算法對脫軌系統(tǒng)進(jìn)行計(jì)算。為了保證模型系統(tǒng)在計(jì)算時順利進(jìn)行,除了設(shè)置上述邊界條件外,還須將積分時間步長設(shè)為10-6s,分析終止時間設(shè)為200 s,文件輸出時間間隔設(shè)為10-3s;為了得到更好的輸出結(jié)果內(nèi)容與形式,設(shè)置二進(jìn)制輸出文件控制卡片(DATABASE_OPTION),主要的輸出控制參數(shù)包括DEFORC、ELOUT、GLSTAT、MATSUM、NODFOR、NODOUT、RCFORC、SLEOUT、SPCFORC、SWFORC、BINARY_D3PLOT。針 對LS?DYNA計(jì)算結(jié)果和所設(shè)置的輸出數(shù)據(jù)進(jìn)行相應(yīng)的數(shù)據(jù)濾波處理,得到各主要指標(biāo)的分析結(jié)果。

    3.1 脫軌過程形態(tài)變化

    車輛在脫軌過程中與軌道的主要形態(tài)變化如圖6所示。可知:脫軌前,車輪與鋼軌之間保持正常的輪軌接觸關(guān)系[13];車輛脫軌時,護(hù)軌側(cè)車輪首先從運(yùn)行軌掉落,同時遠(yuǎn)離護(hù)軌側(cè)的車輪爬上運(yùn)行軌道,被爬鋼軌發(fā)生橫向彎曲;爬軌輪越過鋼軌后,兩側(cè)車輪同時向下掉落至軌道板及軌枕塊,車輪在縱向運(yùn)行的同時繼續(xù)橫向靠近護(hù)軌;護(hù)軌側(cè)車輪與護(hù)軌發(fā)生接觸后,護(hù)軌阻止車輪繼續(xù)橫向移動,護(hù)軌出現(xiàn)橫向彎曲,車輪從軌枕塊掉至軌道板,致使車輛發(fā)生傾斜但未傾覆,最終車輛以傾斜但未傾覆的穩(wěn)態(tài)與護(hù)軌摩擦向前滑行。

    圖6 車輛脫軌過程形態(tài)變化

    3.2 能量變化

    車輛運(yùn)行的動能決定了脫軌后車輛對護(hù)軌的沖擊破壞程度[14]。脫軌過程中車輛的動能變化見圖7??芍涸?~0.05 s系統(tǒng)發(fā)生了第一次動能穩(wěn)定下降,此過程為爬軌側(cè)車輪與鋼軌碰撞摩擦吸能;在0.12 s之后動能穩(wěn)定下降,此時為護(hù)軌側(cè)車輪與護(hù)軌碰撞摩擦吸能,之后的動能均以摩擦熱的形式耗散。

    圖7 脫軌動能隨時間變化曲線

    3.3 輪軌接觸力

    車輛脫軌后輪軌接觸力大小決定了鋼軌給予車輛的橫向沖擊加速度大小,從而決定了護(hù)軌對車內(nèi)乘客與設(shè)備的橫向沖擊程度[15]。輪軌接觸力隨時間變化曲線見圖8??芍涸?.008~0.052 s,爬軌側(cè)車輪與鋼軌接觸造成了輪軌間第一次橫向沖擊,x、y、z方向的沖擊力峰值分別為540、1 000、2 750 kN;0.052 s后爬軌側(cè)車輪從鋼軌上脫落,整個轉(zhuǎn)向架在重力作用下橫向空行;在0.105 s時護(hù)軌側(cè)車輪與護(hù)軌接觸造成第二次輪軌沖擊,x、y、z方向的沖擊力峰值分別為352、1 500、300 kN。

    圖8 輪軌接觸力隨時間變化曲線

    3.4 輪對抬升量

    輪對抬升量決定了車輛脫軌后的危害程度。輪對抬升量(垂向位移)越大,脫軌后的危害程度越大[16]。輪對抬升量隨時間變化曲線見圖9。可知:脫軌過程中,爬軌側(cè)車輪與鋼軌發(fā)生碰撞后,在巨大的橫向壓力與動能作用下,輪緣爬過鋼軌,發(fā)生較大的垂向位移,爬至軌頂時車輪的垂向位移達(dá)到最大值,為17.622 mm;同時護(hù)軌側(cè)車輪垂向位移也達(dá)到最大值,為5.540 mm。爬軌側(cè)和護(hù)軌側(cè)車輪垂向位移均滿足車輛運(yùn)行安全性要求。

    圖9 輪對抬升量隨時間變化曲線

    3.5 護(hù)軌橫向變形

    護(hù)軌的橫向變形決定了其對車輪的保護(hù)能力。接觸點(diǎn)處護(hù)軌橫向側(cè)彎位移隨時間變化曲線及最終形態(tài)見圖10??芍涸?.105 s時車輪開始與護(hù)軌發(fā)生碰撞接觸,護(hù)軌在車輛的橫向擠壓下發(fā)生了側(cè)彎,繼而發(fā)生橫向變形;在0.2 s時護(hù)軌停止橫向變形,橫向側(cè)彎位移保持在31 mm,能夠滿足鐵路軌道設(shè)計(jì)規(guī)范要求。

    圖10 接觸點(diǎn)處護(hù)軌橫向變形情況

    3.6 道床受力狀態(tài)

    脫軌后軌道板的von Mises應(yīng)力分布如圖11所示??芍?,車輪與鋼軌碰撞的整個過程中,軌道板的最大von Mises應(yīng)力為35.17 MPa,發(fā)生位置為軌道板護(hù)軌安裝處靠近線路中心側(cè),能夠滿足C40混凝土的強(qiáng)度要求。

    圖11 軌道板von Mises應(yīng)力分布(單位:MPa)

    3.7 車輛偏移車輛限界情況

    脫軌過程結(jié)束后,車輛以圖12所示的穩(wěn)定狀態(tài)繼續(xù)沿軌道向前滑行,未發(fā)生傾覆,說明護(hù)軌對車輛起到了重要的保護(hù)作用。該狀態(tài)下,車體最外輪廓線橫向偏移529 mm,車輛重心偏移軌道中心線396 mm。

    圖12 車輛偏移形態(tài)(單位:mm)

    4 結(jié)論

    本文以非線性大變形有限元法為理論基礎(chǔ),以三維空間為設(shè)計(jì)與分析計(jì)算環(huán)境,利用BIM軟件設(shè)計(jì)并建立了某軌道及其護(hù)軌結(jié)構(gòu)的BIM模型,并基于BIM模型建立了該軌道與某車輛組成的脫軌系統(tǒng)有限元模型,通過LS?DYNA程序可視化研究了某高速鐵路軌道護(hù)軌裝置對車輛脫軌后的保護(hù)作用和自身的力學(xué)響應(yīng)。主要結(jié)論如下:

    1)車輛脫軌全過程中,輪軌之間的接觸力上升及摩擦吸能集中在爬軌側(cè)車輪與鋼軌接觸和護(hù)軌側(cè)車輪與護(hù)軌接觸兩個階段,考慮了柔性變形的有限元法使得列車在脫軌后撞擊軌枕塊和護(hù)軌的過程具有合理、可靠的傳力過程和運(yùn)動狀態(tài)。

    2)脫軌過程中,爬軌側(cè)車輪爬至軌頂時車輪的垂向位移達(dá)到最大值,為17.622 mm,同時護(hù)軌側(cè)車輪垂向位移也達(dá)到最大值,為5.540 mm,滿足車輛運(yùn)行安全性要求。

    3)車輪與護(hù)軌碰撞接觸過程中,護(hù)軌在車輛的橫向擠壓下發(fā)生側(cè)彎,繼而發(fā)生橫向變形,其最大值為31 mm,滿足鐵路軌道設(shè)計(jì)規(guī)范要求。

    4)脫軌后車輪與鋼軌的碰撞過程中,軌道板的最大von Mises應(yīng)力為35.17 MPa,發(fā)生在軌道板護(hù)軌安裝處靠近線路中心側(cè),滿足C40混凝土強(qiáng)度要求。

    5)脫軌過程結(jié)束后,車體最外輪廓線橫向偏移529 mm,車輛重心偏移軌道中心線396 mm,車輛以穩(wěn)定的傾斜形態(tài)繼續(xù)向前滑行,未發(fā)生傾覆,護(hù)軌對車輛起到了保護(hù)作用。

    猜你喜歡
    軌枕輪軌鋼軌
    裝配式可調(diào)間距的X形抗滑動軌枕研究
    中低速磁浮道岔與輪軌道岔的差異
    鋼軌焊后雙中頻感應(yīng)加熱工藝研究
    軌枕外觀及修補(bǔ)工藝研究
    中低速磁浮與輪軌交通信號系統(tǒng)的差異
    非線性穩(wěn)態(tài)曲線通過時輪軌滾動接觸的數(shù)值求解方法
    高速鐵路鋼軌疲勞過程的超聲非線性系數(shù)表征
    國內(nèi)外高速鐵路鋼軌性能對比研究
    一種新的混凝土軌枕裝卸機(jī)設(shè)計(jì)
    減緩高速鐵路鋼軌波磨的仿真分析
    亚洲自拍偷在线| 一区二区三区四区激情视频 | 国产精华一区二区三区| 亚洲欧美清纯卡通| 国产精品三级大全| 国产黄色小视频在线观看| 搡老妇女老女人老熟妇| 午夜免费激情av| 中文字幕精品亚洲无线码一区| 简卡轻食公司| 内地一区二区视频在线| 免费在线观看日本一区| 内射极品少妇av片p| 精品不卡国产一区二区三区| 国产男靠女视频免费网站| 亚洲美女视频黄频| 日韩强制内射视频| 日本一二三区视频观看| 亚洲avbb在线观看| 亚洲美女搞黄在线观看 | 男女边吃奶边做爰视频| 国产精品一及| 亚洲av第一区精品v没综合| 国产日本99.免费观看| 日日夜夜操网爽| 一区二区三区激情视频| 波多野结衣高清作品| 俺也久久电影网| 日本在线视频免费播放| 国产精品亚洲一级av第二区| 深夜精品福利| 成人av在线播放网站| 亚洲人成网站高清观看| av.在线天堂| 亚洲综合色惰| 日本欧美国产在线视频| 99久久精品一区二区三区| 免费在线观看影片大全网站| 欧美一区二区国产精品久久精品| 很黄的视频免费| a级一级毛片免费在线观看| 麻豆久久精品国产亚洲av| 国产精品国产三级国产av玫瑰| 精品久久久久久久末码| 十八禁国产超污无遮挡网站| 国产午夜福利久久久久久| 国产私拍福利视频在线观看| 真人一进一出gif抽搐免费| 国产又黄又爽又无遮挡在线| 1000部很黄的大片| 淫秽高清视频在线观看| 成人二区视频| 不卡视频在线观看欧美| 国产一区二区在线观看日韩| 中亚洲国语对白在线视频| 精品久久久久久成人av| 精品不卡国产一区二区三区| 国产精品久久久久久精品电影| 国模一区二区三区四区视频| 亚洲在线自拍视频| 亚洲专区中文字幕在线| 色综合站精品国产| 美女高潮的动态| 韩国av一区二区三区四区| 亚洲av中文av极速乱 | 在线观看舔阴道视频| 午夜福利在线在线| 最近中文字幕高清免费大全6 | 久久精品人妻少妇| 免费大片18禁| 亚洲avbb在线观看| 精品一区二区三区视频在线观看免费| 国产欧美日韩精品亚洲av| 国产精品一区www在线观看 | 一边摸一边抽搐一进一小说| 午夜精品一区二区三区免费看| 国产一级毛片七仙女欲春2| 亚洲av美国av| bbb黄色大片| 成人精品一区二区免费| 亚洲国产精品sss在线观看| 俄罗斯特黄特色一大片| x7x7x7水蜜桃| 特级一级黄色大片| av国产免费在线观看| 午夜亚洲福利在线播放| 最新在线观看一区二区三区| 免费看光身美女| 最好的美女福利视频网| 一夜夜www| 国产单亲对白刺激| 国产精品国产高清国产av| 在线观看一区二区三区| 在线观看午夜福利视频| 成人美女网站在线观看视频| 亚洲电影在线观看av| 免费电影在线观看免费观看| 99精品久久久久人妻精品| 欧美潮喷喷水| 日韩一本色道免费dvd| 日本精品一区二区三区蜜桃| 简卡轻食公司| 欧美+日韩+精品| 亚洲美女黄片视频| 免费人成视频x8x8入口观看| 老师上课跳d突然被开到最大视频| 国产高清不卡午夜福利| 十八禁网站免费在线| 久久午夜亚洲精品久久| 国产精品电影一区二区三区| 色在线成人网| 亚洲国产精品成人综合色| 欧美另类亚洲清纯唯美| 在线天堂最新版资源| 五月玫瑰六月丁香| 国产麻豆成人av免费视频| 少妇的逼水好多| 波多野结衣高清无吗| 俄罗斯特黄特色一大片| 校园春色视频在线观看| 在线天堂最新版资源| 欧美黑人欧美精品刺激| bbb黄色大片| 中文资源天堂在线| 在线观看一区二区三区| 99久久精品国产国产毛片| 啦啦啦观看免费观看视频高清| 日韩欧美 国产精品| 欧美日韩中文字幕国产精品一区二区三区| 色视频www国产| 国内精品一区二区在线观看| 国产高清视频在线播放一区| 亚洲av美国av| 两人在一起打扑克的视频| 日韩,欧美,国产一区二区三区 | 国产伦精品一区二区三区四那| 久久欧美精品欧美久久欧美| 男女视频在线观看网站免费| 窝窝影院91人妻| 国产av在哪里看| 欧美日韩乱码在线| 国内精品久久久久精免费| 日本免费一区二区三区高清不卡| 日本爱情动作片www.在线观看 | 欧美三级亚洲精品| 丰满乱子伦码专区| 97超级碰碰碰精品色视频在线观看| 日韩国内少妇激情av| 亚洲美女视频黄频| 精品一区二区三区视频在线| 天天一区二区日本电影三级| 国产精品人妻久久久影院| 一个人观看的视频www高清免费观看| 老师上课跳d突然被开到最大视频| 1000部很黄的大片| 亚洲国产精品成人综合色| 国产在线精品亚洲第一网站| 色在线成人网| 午夜老司机福利剧场| 欧美最黄视频在线播放免费| 五月伊人婷婷丁香| 亚洲成人久久爱视频| 亚洲最大成人中文| 变态另类丝袜制服| 国产久久久一区二区三区| 999久久久精品免费观看国产| 九九热线精品视视频播放| 成年女人毛片免费观看观看9| 波野结衣二区三区在线| 欧美日韩中文字幕国产精品一区二区三区| 国内久久婷婷六月综合欲色啪| 国产精品一区二区性色av| 国产精品精品国产色婷婷| 国产久久久一区二区三区| 成人av一区二区三区在线看| 欧美在线一区亚洲| 欧美潮喷喷水| 一区福利在线观看| 国产日本99.免费观看| 黄色日韩在线| 亚洲无线在线观看| 亚洲久久久久久中文字幕| 国产女主播在线喷水免费视频网站 | 国产精品人妻久久久久久| 五月伊人婷婷丁香| 51国产日韩欧美| 国产精品乱码一区二三区的特点| 亚洲av免费在线观看| 国内精品久久久久久久电影| 精品一区二区三区视频在线观看免费| 欧美zozozo另类| 淫妇啪啪啪对白视频| 国产av在哪里看| 午夜老司机福利剧场| 在线国产一区二区在线| 高清日韩中文字幕在线| 国产精品伦人一区二区| ponron亚洲| 亚洲黑人精品在线| 日韩在线高清观看一区二区三区 | 成人二区视频| 欧美日韩黄片免| 国内毛片毛片毛片毛片毛片| 国产精品电影一区二区三区| 91麻豆精品激情在线观看国产| 在线看三级毛片| 在线观看66精品国产| 国产一区二区亚洲精品在线观看| 九色国产91popny在线| 婷婷精品国产亚洲av在线| 亚洲精品影视一区二区三区av| 国产精品电影一区二区三区| 一区福利在线观看| 狂野欧美白嫩少妇大欣赏| 99久久九九国产精品国产免费| 亚洲精品粉嫩美女一区| 欧美成人免费av一区二区三区| av在线蜜桃| 欧美zozozo另类| 国产激情偷乱视频一区二区| 黄片wwwwww| 在线免费十八禁| 少妇熟女aⅴ在线视频| 国产精品人妻久久久影院| 麻豆av噜噜一区二区三区| 国产亚洲精品久久久久久毛片| 国产一区二区三区在线臀色熟女| 久久久久精品国产欧美久久久| 国内少妇人妻偷人精品xxx网站| 欧美激情在线99| 可以在线观看毛片的网站| 黄片wwwwww| 亚洲国产色片| 亚洲专区国产一区二区| 亚洲最大成人av| 99久久成人亚洲精品观看| 男人舔奶头视频| 亚洲va在线va天堂va国产| 尾随美女入室| 九九热线精品视视频播放| 亚洲天堂国产精品一区在线| 日本与韩国留学比较| 春色校园在线视频观看| 日韩国内少妇激情av| 亚洲精品粉嫩美女一区| 天堂网av新在线| 3wmmmm亚洲av在线观看| 热99re8久久精品国产| 波野结衣二区三区在线| 久久6这里有精品| 可以在线观看毛片的网站| 九九久久精品国产亚洲av麻豆| 国产中年淑女户外野战色| 国产老妇女一区| 偷拍熟女少妇极品色| 精品久久久久久久久久免费视频| 婷婷六月久久综合丁香| 成熟少妇高潮喷水视频| 久久人人精品亚洲av| 亚洲av免费高清在线观看| 人妻制服诱惑在线中文字幕| 国产一区二区三区在线臀色熟女| 啦啦啦韩国在线观看视频| 99久久精品热视频| 亚洲av不卡在线观看| 久久热精品热| 别揉我奶头~嗯~啊~动态视频| 网址你懂的国产日韩在线| 联通29元200g的流量卡| 欧美日韩精品成人综合77777| 人妻少妇偷人精品九色| 中文亚洲av片在线观看爽| 日韩欧美免费精品| 搡老熟女国产l中国老女人| 99在线人妻在线中文字幕| 在线播放国产精品三级| 一区福利在线观看| 亚洲av电影不卡..在线观看| 亚洲自拍偷在线| 18禁黄网站禁片午夜丰满| 精品不卡国产一区二区三区| 嫩草影院新地址| 成人一区二区视频在线观看| 99久久精品热视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲欧美日韩无卡精品| aaaaa片日本免费| 三级国产精品欧美在线观看| 国产亚洲精品久久久com| 国语自产精品视频在线第100页| 黄色视频,在线免费观看| 久久久精品大字幕| 久久久久久久久大av| 女人被狂操c到高潮| 男女啪啪激烈高潮av片| 国产精品美女特级片免费视频播放器| 亚洲一区二区三区色噜噜| 中亚洲国语对白在线视频| 日本熟妇午夜| 欧美日韩乱码在线| 内射极品少妇av片p| a级毛片a级免费在线| 淫妇啪啪啪对白视频| 久久久久久伊人网av| 久久草成人影院| 婷婷精品国产亚洲av在线| 简卡轻食公司| 身体一侧抽搐| 久久国内精品自在自线图片| 亚洲av二区三区四区| 欧美日韩乱码在线| 尾随美女入室| 丰满乱子伦码专区| 九九爱精品视频在线观看| 一进一出抽搐gif免费好疼| 高清日韩中文字幕在线| 91麻豆av在线| 免费看a级黄色片| 人妻少妇偷人精品九色| 日韩人妻高清精品专区| 99精品久久久久人妻精品| 日本-黄色视频高清免费观看| 男女视频在线观看网站免费| 黄色日韩在线| 在线观看午夜福利视频| 99视频精品全部免费 在线| 亚洲av熟女| 国产伦在线观看视频一区| 国产精品日韩av在线免费观看| 欧美一级a爱片免费观看看| 色综合亚洲欧美另类图片| 欧美高清性xxxxhd video| 深夜a级毛片| 桃色一区二区三区在线观看| 免费看a级黄色片| 午夜激情福利司机影院| 嫩草影视91久久| 18+在线观看网站| 久久精品国产亚洲av香蕉五月| a级毛片a级免费在线| 亚洲美女黄片视频| 男人和女人高潮做爰伦理| 免费看光身美女| 欧美日韩中文字幕国产精品一区二区三区| 国产在视频线在精品| 精品久久久久久久久亚洲 | 99精品在免费线老司机午夜| 亚洲乱码一区二区免费版| 两个人视频免费观看高清| 久久国产精品人妻蜜桃| 91av网一区二区| 在线看三级毛片| 12—13女人毛片做爰片一| 国产精品久久久久久av不卡| 国产av一区在线观看免费| 无遮挡黄片免费观看| 人妻少妇偷人精品九色| 熟女电影av网| 午夜亚洲福利在线播放| 真人做人爱边吃奶动态| 免费不卡的大黄色大毛片视频在线观看 | 97人妻精品一区二区三区麻豆| 69人妻影院| 亚洲成a人片在线一区二区| 国产高清不卡午夜福利| 精品免费久久久久久久清纯| 99久久精品一区二区三区| 久久久成人免费电影| 99在线视频只有这里精品首页| 日韩国内少妇激情av| 色综合亚洲欧美另类图片| 午夜爱爱视频在线播放| 色综合色国产| 久久人人爽人人爽人人片va| 国产69精品久久久久777片| av天堂中文字幕网| 色综合站精品国产| 黄色丝袜av网址大全| 免费看a级黄色片| 给我免费播放毛片高清在线观看| 看黄色毛片网站| 香蕉av资源在线| 久久久国产成人精品二区| 亚洲色图av天堂| www日本黄色视频网| 在线免费十八禁| 精品福利观看| 亚洲av免费在线观看| 免费黄网站久久成人精品| ponron亚洲| 一边摸一边抽搐一进一小说| 欧美激情国产日韩精品一区| 久久99热这里只有精品18| 成年女人毛片免费观看观看9| 色在线成人网| 国产精品久久电影中文字幕| 国产精品免费一区二区三区在线| 亚洲在线自拍视频| 欧美一区二区国产精品久久精品| 国产成人一区二区在线| av视频在线观看入口| 网址你懂的国产日韩在线| 国模一区二区三区四区视频| 国产精品久久久久久久电影| 1024手机看黄色片| 91午夜精品亚洲一区二区三区 | 少妇被粗大猛烈的视频| 九九爱精品视频在线观看| 免费av毛片视频| 国产色爽女视频免费观看| 亚洲av成人av| 中文字幕高清在线视频| 舔av片在线| 国产探花极品一区二区| 一个人观看的视频www高清免费观看| av天堂中文字幕网| 麻豆国产av国片精品| 欧美一区二区精品小视频在线| 能在线免费观看的黄片| 精品久久久久久,| 校园春色视频在线观看| 97超视频在线观看视频| 国产欧美日韩精品一区二区| 噜噜噜噜噜久久久久久91| 日日摸夜夜添夜夜添av毛片 | 久久精品91蜜桃| 午夜精品久久久久久毛片777| 久久热精品热| 色视频www国产| 亚洲成人中文字幕在线播放| 精品人妻一区二区三区麻豆 | 免费观看人在逋| 99久久九九国产精品国产免费| 亚洲精品久久国产高清桃花| 91久久精品电影网| 亚洲成a人片在线一区二区| 欧美日韩精品成人综合77777| 色尼玛亚洲综合影院| 国产精品人妻久久久影院| 亚洲av中文字字幕乱码综合| 免费看美女性在线毛片视频| 免费av不卡在线播放| 男人的好看免费观看在线视频| 国产高潮美女av| 精品人妻熟女av久视频| 色在线成人网| 国内久久婷婷六月综合欲色啪| 亚洲国产精品sss在线观看| 精品无人区乱码1区二区| 少妇人妻一区二区三区视频| 欧美成人免费av一区二区三区| 国产av在哪里看| 国产成人av教育| 色在线成人网| av黄色大香蕉| 国产一区二区在线av高清观看| 尾随美女入室| 午夜福利在线观看吧| 成年女人看的毛片在线观看| aaaaa片日本免费| 日日摸夜夜添夜夜添av毛片 | 特级一级黄色大片| 亚洲五月天丁香| 全区人妻精品视频| 18禁黄网站禁片午夜丰满| 亚洲国产精品sss在线观看| 久久精品国产鲁丝片午夜精品 | 国产精品亚洲美女久久久| 男女视频在线观看网站免费| 我的老师免费观看完整版| 午夜爱爱视频在线播放| 97人妻精品一区二区三区麻豆| 国产av一区在线观看免费| 亚洲内射少妇av| 国产免费男女视频| 男人的好看免费观看在线视频| 中文字幕av成人在线电影| 亚洲精品456在线播放app | 午夜日韩欧美国产| av视频在线观看入口| 欧美三级亚洲精品| 桃色一区二区三区在线观看| 亚洲综合色惰| 日韩一区二区视频免费看| 我要看日韩黄色一级片| 国产白丝娇喘喷水9色精品| 日韩亚洲欧美综合| 亚洲国产精品sss在线观看| 国产在视频线在精品| 久久精品综合一区二区三区| 老熟妇乱子伦视频在线观看| 国产淫片久久久久久久久| 自拍偷自拍亚洲精品老妇| 色尼玛亚洲综合影院| 亚洲精品粉嫩美女一区| 欧美激情久久久久久爽电影| 亚洲最大成人手机在线| 国产一区二区三区av在线 | 中文字幕免费在线视频6| 国产精品国产三级国产av玫瑰| 午夜日韩欧美国产| 啪啪无遮挡十八禁网站| 国内少妇人妻偷人精品xxx网站| 91麻豆av在线| 国产亚洲91精品色在线| 天天一区二区日本电影三级| 亚洲欧美清纯卡通| 午夜影院日韩av| 中文在线观看免费www的网站| 亚洲无线在线观看| .国产精品久久| 18禁黄网站禁片午夜丰满| 午夜a级毛片| 国产欧美日韩精品一区二区| 亚洲aⅴ乱码一区二区在线播放| 国产精品一区二区性色av| 久久精品国产99精品国产亚洲性色| 国产精品久久久久久精品电影| 成人精品一区二区免费| 亚洲精华国产精华精| 亚洲专区中文字幕在线| 亚洲不卡免费看| 亚洲国产欧洲综合997久久,| 伦精品一区二区三区| 日本一二三区视频观看| 国产精品无大码| 精品久久久久久久久久免费视频| 极品教师在线免费播放| 成人特级黄色片久久久久久久| 99热这里只有是精品在线观看| 久久久久久伊人网av| 亚洲成av人片在线播放无| 性插视频无遮挡在线免费观看| 国产高潮美女av| 亚洲av电影不卡..在线观看| 麻豆国产97在线/欧美| 天天一区二区日本电影三级| 成年人黄色毛片网站| 日本色播在线视频| 淫妇啪啪啪对白视频| 亚洲av第一区精品v没综合| 麻豆久久精品国产亚洲av| 一级av片app| 欧美黑人欧美精品刺激| 日韩,欧美,国产一区二区三区 | 日韩中字成人| 大又大粗又爽又黄少妇毛片口| 亚洲av熟女| 12—13女人毛片做爰片一| 中文字幕免费在线视频6| 国产一区二区激情短视频| 国产精品一区二区三区四区免费观看 | 美女 人体艺术 gogo| 国产高清激情床上av| 国产激情偷乱视频一区二区| 色精品久久人妻99蜜桃| 亚洲真实伦在线观看| 国产高清三级在线| av在线天堂中文字幕| 日本-黄色视频高清免费观看| 久久草成人影院| 九色成人免费人妻av| 国语自产精品视频在线第100页| 熟女人妻精品中文字幕| 亚洲综合色惰| 国产精品99久久久久久久久| 午夜福利视频1000在线观看| 无遮挡黄片免费观看| 大又大粗又爽又黄少妇毛片口| 大型黄色视频在线免费观看| 亚洲不卡免费看| 最近视频中文字幕2019在线8| 小说图片视频综合网站| 成人鲁丝片一二三区免费| 亚洲精品色激情综合| 看免费成人av毛片| 国产麻豆成人av免费视频| 亚洲av.av天堂| 国产午夜精品论理片| 1024手机看黄色片| 搡女人真爽免费视频火全软件 | 成年版毛片免费区| 国产成人av教育| 久久中文看片网| 久久久久久久久久成人| 99精品久久久久人妻精品| 日本 av在线| 久久久久久久久中文| 桃红色精品国产亚洲av| eeuss影院久久| 午夜日韩欧美国产| 国内精品宾馆在线| 午夜福利欧美成人| 亚洲精华国产精华液的使用体验 | 精品久久久久久久末码| 国产精品福利在线免费观看| 亚洲在线自拍视频| 少妇猛男粗大的猛烈进出视频 | 成人毛片a级毛片在线播放| 久久精品国产鲁丝片午夜精品 | 一区二区三区激情视频| 亚洲av不卡在线观看| 国产精品一区二区三区四区久久| 国产乱人伦免费视频| aaaaa片日本免费| 国产欧美日韩精品亚洲av| av福利片在线观看| 免费看日本二区| 亚洲av日韩精品久久久久久密| 日本黄色视频三级网站网址| 亚洲av不卡在线观看| 最近视频中文字幕2019在线8| 人妻夜夜爽99麻豆av| av黄色大香蕉| 欧美成人性av电影在线观看|