• <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| 一进一出抽搐gif免费好疼 | 免费看a级黄色片| 韩国精品一区二区三区| 色哟哟哟哟哟哟| 成人特级黄色片久久久久久久| 黄色 视频免费看| 久久人人97超碰香蕉20202| 国产av在哪里看| 亚洲成人精品中文字幕电影 | 中文字幕人妻丝袜制服| а√天堂www在线а√下载| 午夜福利在线免费观看网站| 欧美日韩黄片免| 亚洲成国产人片在线观看| 久久久久久亚洲精品国产蜜桃av| 国产国语露脸激情在线看| 视频区图区小说| 国产97色在线日韩免费| 国产精品久久电影中文字幕| 欧美最黄视频在线播放免费 | 久久精品aⅴ一区二区三区四区| 国产av精品麻豆| 真人一进一出gif抽搐免费| svipshipincom国产片| 91九色精品人成在线观看| 18禁观看日本| 女生性感内裤真人,穿戴方法视频| 国产精品偷伦视频观看了| 校园春色视频在线观看| 欧美不卡视频在线免费观看 | 校园春色视频在线观看| a级毛片黄视频| 国产成人系列免费观看| 日本a在线网址| e午夜精品久久久久久久| 在线观看免费高清a一片| 少妇粗大呻吟视频| 久久精品人人爽人人爽视色| 国产麻豆69| 免费在线观看日本一区| 韩国精品一区二区三区| 五月开心婷婷网| av福利片在线| 日韩大尺度精品在线看网址 | 亚洲第一青青草原| 午夜精品在线福利| 国产激情欧美一区二区| 韩国精品一区二区三区| 99国产综合亚洲精品| x7x7x7水蜜桃| 9热在线视频观看99| 涩涩av久久男人的天堂| 最近最新中文字幕大全免费视频| 男女高潮啪啪啪动态图| 悠悠久久av| 亚洲av成人av| 热re99久久精品国产66热6| 在线观看一区二区三区| 高清毛片免费观看视频网站 | 亚洲国产精品sss在线观看 | 成人永久免费在线观看视频| 久久欧美精品欧美久久欧美| 999久久久国产精品视频| 婷婷丁香在线五月| √禁漫天堂资源中文www| 日日爽夜夜爽网站| 亚洲午夜精品一区,二区,三区| 12—13女人毛片做爰片一| av视频免费观看在线观看| 亚洲人成77777在线视频| 色精品久久人妻99蜜桃| 热99re8久久精品国产| 少妇的丰满在线观看| 国产av一区在线观看免费| 国产乱人伦免费视频| 久久精品aⅴ一区二区三区四区| 人成视频在线观看免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 99国产精品一区二区三区| 巨乳人妻的诱惑在线观看| 黄频高清免费视频| svipshipincom国产片| av免费在线观看网站| 国产一区在线观看成人免费| 久久久久九九精品影院| 亚洲精品av麻豆狂野| 国产高清视频在线播放一区| 亚洲人成77777在线视频| 性少妇av在线| 午夜日韩欧美国产| 久久国产精品男人的天堂亚洲| 久久人妻熟女aⅴ| 成人国语在线视频| 777久久人妻少妇嫩草av网站| 一区二区日韩欧美中文字幕| 欧美精品一区二区免费开放| 9191精品国产免费久久| 午夜a级毛片| 夜夜爽天天搞| 淫妇啪啪啪对白视频| 亚洲精品中文字幕在线视频| 国产一区二区三区在线臀色熟女 | 久久久久久人人人人人| 亚洲性夜色夜夜综合| av片东京热男人的天堂| 在线十欧美十亚洲十日本专区| avwww免费| 亚洲精华国产精华精| 国产高清国产精品国产三级| 午夜视频精品福利| 极品教师在线免费播放| 黑人猛操日本美女一级片| 午夜两性在线视频| 国产精品久久久av美女十八| 男人舔女人下体高潮全视频| 色播在线永久视频| 又紧又爽又黄一区二区| 精品国产超薄肉色丝袜足j| 免费少妇av软件| 啦啦啦 在线观看视频| 久久精品91无色码中文字幕| 69精品国产乱码久久久| 首页视频小说图片口味搜索| 丁香六月欧美| 在线观看一区二区三区激情| 69av精品久久久久久| 多毛熟女@视频| 老司机在亚洲福利影院| 一级黄色大片毛片| 夫妻午夜视频| x7x7x7水蜜桃| 国产高清国产精品国产三级| 日日爽夜夜爽网站| 欧美日韩av久久| 国产精品影院久久| 国产精品免费一区二区三区在线| 男女之事视频高清在线观看| 久久香蕉激情| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产看品久久| 一级,二级,三级黄色视频| 制服人妻中文乱码| 亚洲av成人不卡在线观看播放网| 少妇裸体淫交视频免费看高清 | 亚洲一区中文字幕在线| 亚洲成人国产一区在线观看| 女警被强在线播放| 一边摸一边抽搐一进一出视频| 国产欧美日韩一区二区三| 中文欧美无线码| 一级毛片精品| 国产又色又爽无遮挡免费看| 男人操女人黄网站| 久久青草综合色| 国产99白浆流出| 国产熟女xx| 国产极品粉嫩免费观看在线| 狂野欧美激情性xxxx| 脱女人内裤的视频| 757午夜福利合集在线观看| 欧美日韩av久久| 午夜两性在线视频| 国产午夜精品久久久久久| 男人的好看免费观看在线视频 | 国产日韩一区二区三区精品不卡| 欧美日韩av久久| 亚洲欧美日韩无卡精品| a级毛片在线看网站| 狂野欧美激情性xxxx| 欧美性长视频在线观看| 男女午夜视频在线观看| 国产精品九九99| 久99久视频精品免费| 人人澡人人妻人| 午夜两性在线视频| 老汉色av国产亚洲站长工具| 中文欧美无线码| 9色porny在线观看| 国产精品98久久久久久宅男小说| 97人妻天天添夜夜摸| 亚洲国产精品合色在线| 一区福利在线观看| 电影成人av| e午夜精品久久久久久久| 亚洲av成人一区二区三| av中文乱码字幕在线| 中文字幕人妻丝袜一区二区| 国产精品av久久久久免费| 这个男人来自地球电影免费观看| 欧美日韩国产mv在线观看视频| 国产男靠女视频免费网站| 亚洲中文av在线| 亚洲avbb在线观看| 亚洲熟女毛片儿| 在线永久观看黄色视频| 首页视频小说图片口味搜索| 18禁观看日本| 成人手机av| av电影中文网址| 黄色女人牲交| 久久草成人影院| 国产成人欧美| 久久久久国内视频| 亚洲国产精品合色在线| 成在线人永久免费视频| 精品国产乱子伦一区二区三区| 久久人人精品亚洲av| 国产成人影院久久av| 热re99久久国产66热| 国产主播在线观看一区二区| 最新在线观看一区二区三区| 99久久久亚洲精品蜜臀av| 日本a在线网址| 欧美最黄视频在线播放免费 | 国产欧美日韩综合在线一区二区| 一区二区三区激情视频| 老鸭窝网址在线观看| 久久国产亚洲av麻豆专区| av欧美777| 国产极品粉嫩免费观看在线| 久久久久久免费高清国产稀缺| 免费观看人在逋| 午夜免费鲁丝| 国产高清激情床上av| 日本wwww免费看| 九色亚洲精品在线播放| 性色av乱码一区二区三区2| 人妻丰满熟妇av一区二区三区| 国产一区在线观看成人免费| 精品熟女少妇八av免费久了| 黄网站色视频无遮挡免费观看| 日韩视频一区二区在线观看| 精品久久久久久久久久免费视频 | 女警被强在线播放| 国产亚洲精品久久久久久毛片| 成人免费观看视频高清| 日韩大尺度精品在线看网址 | www.自偷自拍.com| 女人爽到高潮嗷嗷叫在线视频| 亚洲专区中文字幕在线| 90打野战视频偷拍视频| 一级片'在线观看视频| 757午夜福利合集在线观看| 日本免费一区二区三区高清不卡 | 久久久国产精品麻豆| 极品人妻少妇av视频| www日本在线高清视频| 一a级毛片在线观看| 久久精品国产亚洲av高清一级| av片东京热男人的天堂| 国产99白浆流出| 久久久国产精品麻豆| 首页视频小说图片口味搜索| 亚洲欧美日韩另类电影网站| 黄片小视频在线播放| 在线观看日韩欧美| 亚洲自偷自拍图片 自拍| 天天影视国产精品| 午夜福利在线观看吧| 精品国内亚洲2022精品成人| 啦啦啦免费观看视频1| 视频在线观看一区二区三区| 国产成人免费无遮挡视频| 在线视频色国产色| 99riav亚洲国产免费| 啦啦啦在线免费观看视频4| 久久人人97超碰香蕉20202| 99国产综合亚洲精品| 欧洲精品卡2卡3卡4卡5卡区| 99久久久亚洲精品蜜臀av| 国产av精品麻豆| 高清av免费在线| 精品国产超薄肉色丝袜足j| 亚洲七黄色美女视频| 国产日韩一区二区三区精品不卡| 亚洲精品一区av在线观看| 51午夜福利影视在线观看| 一进一出抽搐动态| 久久精品亚洲熟妇少妇任你| 中文字幕人妻丝袜一区二区| 久久久国产成人精品二区 | 性色av乱码一区二区三区2| 高清在线国产一区| 午夜免费成人在线视频| 亚洲人成电影观看| 亚洲欧美精品综合一区二区三区| 女生性感内裤真人,穿戴方法视频| 少妇的丰满在线观看| 中文亚洲av片在线观看爽| 一夜夜www| 久久热在线av| 日本三级黄在线观看| 女人被躁到高潮嗷嗷叫费观| 久久久久久久精品吃奶| 窝窝影院91人妻| 久热这里只有精品99| 久久香蕉激情| 黄色成人免费大全| 中文字幕另类日韩欧美亚洲嫩草| 怎么达到女性高潮| 美女高潮喷水抽搐中文字幕| 校园春色视频在线观看| ponron亚洲| 高清毛片免费观看视频网站 | 搡老岳熟女国产| 午夜免费激情av| 精品国产乱子伦一区二区三区| 国产高清激情床上av| 欧美人与性动交α欧美精品济南到| 亚洲午夜理论影院| 欧美黄色片欧美黄色片| 成人国语在线视频| 亚洲精品一区av在线观看| 这个男人来自地球电影免费观看| 国产亚洲精品一区二区www| 色哟哟哟哟哟哟| 欧美日韩亚洲高清精品| 午夜久久久在线观看| 最近最新中文字幕大全电影3 | 首页视频小说图片口味搜索| 中文字幕高清在线视频| 欧美激情久久久久久爽电影 | 激情在线观看视频在线高清| 亚洲欧美日韩无卡精品| 在线观看一区二区三区| 人人妻人人爽人人添夜夜欢视频| 高清毛片免费观看视频网站 | 两性午夜刺激爽爽歪歪视频在线观看 | 看片在线看免费视频| 身体一侧抽搐| 女性生殖器流出的白浆| www国产在线视频色| 又紧又爽又黄一区二区| 亚洲av成人av| 国产欧美日韩综合在线一区二区| 淫妇啪啪啪对白视频| 日韩中文字幕欧美一区二区| 久久精品国产清高在天天线| 国产精品久久久久久人妻精品电影| 成人亚洲精品av一区二区 | 日韩 欧美 亚洲 中文字幕| 黄频高清免费视频| 亚洲男人的天堂狠狠| 中文亚洲av片在线观看爽| 亚洲精品粉嫩美女一区| 欧美日韩一级在线毛片| 黑人操中国人逼视频| 神马国产精品三级电影在线观看 | 亚洲伊人色综图| 国内毛片毛片毛片毛片毛片| 亚洲精品久久午夜乱码| 久久精品影院6| 亚洲国产中文字幕在线视频| 精品一品国产午夜福利视频| 可以免费在线观看a视频的电影网站| 亚洲美女黄片视频| 一进一出抽搐动态| 亚洲少妇的诱惑av| 99在线视频只有这里精品首页| 精品一区二区三区视频在线观看免费 | 亚洲一码二码三码区别大吗| 黑人操中国人逼视频| 亚洲人成77777在线视频| 国产精品一区二区在线不卡| 国产av一区二区精品久久| 国产97色在线日韩免费| 在线观看一区二区三区激情| 久久天躁狠狠躁夜夜2o2o| 欧美人与性动交α欧美软件| 成人免费观看视频高清| 色播在线永久视频| 99国产综合亚洲精品| 脱女人内裤的视频| 啦啦啦 在线观看视频| 久久久久精品国产欧美久久久| 国产精品国产av在线观看| 国产精品亚洲av一区麻豆| 无遮挡黄片免费观看| 99精品久久久久人妻精品| ponron亚洲| 男人的好看免费观看在线视频 | 免费在线观看完整版高清| 搡老乐熟女国产| 亚洲国产精品999在线| 首页视频小说图片口味搜索| 美女大奶头视频| 国产成人av教育| 成人国产一区最新在线观看| 午夜精品久久久久久毛片777| ponron亚洲| 久久午夜亚洲精品久久| 成人永久免费在线观看视频| 国产成人精品在线电影| 一级毛片高清免费大全| 亚洲精品粉嫩美女一区| 在线国产一区二区在线| 美女高潮到喷水免费观看| 久久国产亚洲av麻豆专区| 宅男免费午夜| xxxhd国产人妻xxx| 国产精品二区激情视频| 欧美最黄视频在线播放免费 | 亚洲精品中文字幕一二三四区| av电影中文网址| 欧美日本中文国产一区发布| 男男h啪啪无遮挡| 午夜福利在线免费观看网站| 麻豆一二三区av精品| 极品教师在线免费播放| 色播在线永久视频| 亚洲av第一区精品v没综合| 色老头精品视频在线观看| 少妇裸体淫交视频免费看高清 | 亚洲国产精品999在线| 成人三级做爰电影| 一区二区日韩欧美中文字幕| 美女扒开内裤让男人捅视频| 国产av一区二区精品久久| 久久久久国产一级毛片高清牌| 精品一区二区三卡| 精品午夜福利视频在线观看一区| 男女做爰动态图高潮gif福利片 | 国产不卡一卡二| 法律面前人人平等表现在哪些方面| 午夜免费观看网址| 亚洲色图综合在线观看| 露出奶头的视频| 亚洲色图综合在线观看| 欧美激情高清一区二区三区| 丝袜美足系列| 51午夜福利影视在线观看| av免费在线观看网站| 夜夜爽天天搞| 国产一区二区三区视频了| 久久精品成人免费网站| 国产高清视频在线播放一区| 丝袜美腿诱惑在线| 在线播放国产精品三级| 女人被躁到高潮嗷嗷叫费观| 香蕉丝袜av| 国产成人av激情在线播放| 亚洲自拍偷在线| 97超级碰碰碰精品色视频在线观看| 日韩免费av在线播放| 成人国语在线视频| 亚洲自拍偷在线| 校园春色视频在线观看| 人人妻人人爽人人添夜夜欢视频| 看免费av毛片| 丝袜美腿诱惑在线| 不卡av一区二区三区| 中亚洲国语对白在线视频| 国产精品久久久av美女十八| 18禁黄网站禁片午夜丰满| 老汉色∧v一级毛片| 国产麻豆69| ponron亚洲| 国产区一区二久久| 亚洲成人免费av在线播放| 精品久久久久久电影网| 国产一区在线观看成人免费| 婷婷丁香在线五月| 久久国产乱子伦精品免费另类| 国产激情久久老熟女| 国产不卡一卡二| 两个人看的免费小视频| 十八禁网站免费在线| 黑人欧美特级aaaaaa片| 国产亚洲欧美在线一区二区| 成人精品一区二区免费| 国产亚洲欧美在线一区二区| 日韩国内少妇激情av| 欧美日韩乱码在线| 欧美性长视频在线观看| 又大又爽又粗| 亚洲中文字幕日韩| 国产成人影院久久av| 好男人电影高清在线观看| 操出白浆在线播放| 日本免费一区二区三区高清不卡 | 亚洲欧美精品综合久久99| 丁香欧美五月| 成人18禁高潮啪啪吃奶动态图| 99国产精品免费福利视频| 精品国产乱子伦一区二区三区| 俄罗斯特黄特色一大片| 欧美日韩瑟瑟在线播放| 男人操女人黄网站| 国产成人系列免费观看| 久久久久国产一级毛片高清牌| 久久中文看片网| 男男h啪啪无遮挡| 成人亚洲精品一区在线观看| 99久久99久久久精品蜜桃| 嫩草影院精品99| 激情在线观看视频在线高清| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲 国产 在线| 熟女少妇亚洲综合色aaa.| 88av欧美| 日日摸夜夜添夜夜添小说| 国产亚洲精品久久久久5区| 侵犯人妻中文字幕一二三四区| 国产精品偷伦视频观看了| av天堂在线播放| 夫妻午夜视频| 99精品欧美一区二区三区四区| 亚洲国产欧美网| av超薄肉色丝袜交足视频| 成人特级黄色片久久久久久久| 精品无人区乱码1区二区| 亚洲欧美日韩无卡精品| 人人澡人人妻人| 一个人免费在线观看的高清视频| 亚洲va日本ⅴa欧美va伊人久久| 国产色视频综合| 女人高潮潮喷娇喘18禁视频| 成年女人毛片免费观看观看9| 国产一卡二卡三卡精品| 久久久久久久久免费视频了| 欧美激情高清一区二区三区| 激情视频va一区二区三区| 日本撒尿小便嘘嘘汇集6| 一进一出好大好爽视频| 午夜91福利影院| 欧美激情高清一区二区三区| 亚洲精品在线美女| 欧美丝袜亚洲另类 | 欧美性长视频在线观看| 18美女黄网站色大片免费观看| ponron亚洲| 欧美丝袜亚洲另类 | 亚洲 欧美 日韩 在线 免费| 制服诱惑二区| 每晚都被弄得嗷嗷叫到高潮| 91麻豆av在线| 国产蜜桃级精品一区二区三区| 午夜日韩欧美国产| 欧美日韩乱码在线| 国产乱人伦免费视频| 精品熟女少妇八av免费久了| www.自偷自拍.com| 欧美中文综合在线视频| 国产91精品成人一区二区三区| 两个人看的免费小视频| 精品免费久久久久久久清纯| 777久久人妻少妇嫩草av网站| 午夜老司机福利片| 亚洲自拍偷在线| av天堂久久9| 在线天堂中文资源库| 在线观看舔阴道视频| 丝袜美腿诱惑在线| 在线观看66精品国产| 丝袜人妻中文字幕| 国产亚洲精品第一综合不卡| 精品日产1卡2卡| 亚洲精品久久成人aⅴ小说| 法律面前人人平等表现在哪些方面| 美女高潮喷水抽搐中文字幕| 欧美日韩瑟瑟在线播放| av超薄肉色丝袜交足视频| 欧美日韩精品网址| 久久这里只有精品19| 精品国产一区二区三区四区第35| 久久 成人 亚洲| 亚洲人成77777在线视频| 亚洲 欧美 日韩 在线 免费| 精品无人区乱码1区二区| 午夜影院日韩av| 久久人人爽av亚洲精品天堂| 久久香蕉国产精品| 亚洲欧美精品综合一区二区三区| 日韩一卡2卡3卡4卡2021年| 久热这里只有精品99| 欧美乱码精品一区二区三区| av超薄肉色丝袜交足视频| 男女午夜视频在线观看| 久久久久久久久免费视频了| 中文字幕人妻丝袜一区二区| 成年版毛片免费区| 精品久久久久久电影网| 欧美黑人精品巨大| 久久 成人 亚洲| 日韩高清综合在线| 亚洲av成人不卡在线观看播放网| 两个人看的免费小视频| 男女下面插进去视频免费观看| 国产亚洲av高清不卡| 亚洲成国产人片在线观看| 久久香蕉国产精品| 欧美国产精品va在线观看不卡| 在线观看日韩欧美| 91麻豆精品激情在线观看国产 | 大陆偷拍与自拍| 999精品在线视频| 久久这里只有精品19| 亚洲成人国产一区在线观看| 欧美人与性动交α欧美软件| 国产极品粉嫩免费观看在线| 欧美在线黄色| 一级毛片高清免费大全| 国产精品偷伦视频观看了| 国产不卡一卡二| 啦啦啦在线免费观看视频4| a在线观看视频网站| 国产不卡一卡二| 国产99久久九九免费精品|