亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        城市軌道交通路基沉降與無(wú)砟軌道軌面不平順映射關(guān)系

        2020-02-26 09:45:10張克平和振興石廣田賈曉紅翟志浩
        關(guān)鍵詞:變形結(jié)構(gòu)

        張克平,和振興,石廣田,賈曉紅,翟志浩

        (1.蘭州交通大學(xué)機(jī)電工程學(xué)院,蘭州 730070; 2.中鐵第一勘察設(shè)計(jì)院集團(tuán)有限公司,西安 710043)

        由于城市軌道交通具有運(yùn)量大、效率高、運(yùn)行準(zhǔn)時(shí)及乘坐舒適等特性,現(xiàn)已成為我國(guó)大中型城市的主要交通工具之一[1]。其中,路基是城市軌道交通線路工程的基礎(chǔ),同時(shí)也是最薄弱和最不穩(wěn)定的結(jié)構(gòu)。在列車(chē)動(dòng)載和重力以及環(huán)境等復(fù)雜因素的綜合作用下,路基不可避免地產(chǎn)生不均勻沉降[2]。路基不均勻沉降發(fā)展到一定程度,上部軌面也將隨之下沉,使軌面產(chǎn)生不平順,軌面不平順又將加劇輪軌間動(dòng)荷載,進(jìn)而影響行車(chē)安全性和乘坐舒適性。

        針對(duì)路基沉降問(wèn)題,國(guó)內(nèi)外學(xué)者開(kāi)展了大量的研究工作。翟婉明等[3]針對(duì)高速鐵路基礎(chǔ)結(jié)構(gòu)性能演變對(duì)軌道結(jié)構(gòu)的損傷等關(guān)鍵問(wèn)題進(jìn)行了分析,提出了高速鐵路基礎(chǔ)結(jié)構(gòu)不均勻沉降與軌面幾何變形的映射關(guān)系;Dahlberg[4]在總結(jié)各國(guó)典型有砟軌道沉降模型的基礎(chǔ)上,利用有限元軟件對(duì)有砟軌道的長(zhǎng)期累積變形進(jìn)行了模擬分析;Li等[5]應(yīng)用車(chē)輛-軌道耦合動(dòng)力學(xué)理論和有限元方法,通過(guò)建立道岔區(qū)沉降預(yù)測(cè)模型,對(duì)道床沉降演化與軌道不平順間的耦合關(guān)系進(jìn)行研究;肖威等[6]基于有限元方法建立土質(zhì)路基上CRTSⅢ型板式無(wú)砟軌道系統(tǒng)空間耦合模型,針對(duì)路基不均勻沉降作用下板式無(wú)砟軌道結(jié)構(gòu)受力變形特性以及軌道脫空變化規(guī)律展開(kāi)研究;郭宇等[7-8]通過(guò)建立路基上考慮層間接觸方式的板式無(wú)砟軌道-路基有限元實(shí)體模型,分別就路基沉降引起縱連板式無(wú)砟軌道和單元板式無(wú)砟軌道軌面幾何變形情況展開(kāi)研究;陳兆偉等[9-10]分別就橋上鋪設(shè)單元板式無(wú)砟軌道和縱連板式無(wú)砟軌道結(jié)構(gòu)條件下就橋墩沉降與鋼軌軌面變形間的映射關(guān)系進(jìn)行了理論推導(dǎo)。周洺漢等[11]基于路基沉降的現(xiàn)場(chǎng)檢測(cè)和理論分析的方法,提出集成工程治理的措施。

        關(guān)于路基沉降方面的研究目前大多主要集中于鐵路,且研究方法往往只考慮靜平衡狀態(tài)下的軌道自重荷載,有關(guān)城市軌道交通在軌道自重荷載和列車(chē)動(dòng)荷載作用下路基沉降與軌面不平順間的映射關(guān)系研究相對(duì)較少?;诖?,本文基于溫克爾彈性地基耦合梁理論和有限元方法[12-14],建立考慮軌道層間連接特性的整體道床軌道梁-體空間有限元實(shí)體模型,就城市軌道交通線路路基沉降引起無(wú)砟軌道軌面不平順變化規(guī)律進(jìn)行分析,進(jìn)而為城市軌道交通線路路基沉降的控制和無(wú)砟軌道線路的養(yǎng)護(hù)維修提供參考,并為后期的動(dòng)力學(xué)研究工作提供理論依據(jù)。

        1 建立分析模型

        1.1 無(wú)砟軌道模型

        選用我國(guó)城市軌道交通線路廣泛采用的鋼筋混凝土整體道床典型軌道結(jié)構(gòu),其結(jié)構(gòu)組成主要包括鋼軌、高彈性扣件、預(yù)制混凝土支撐塊(軌枕)、混凝土道床板和混凝土支承層,如圖1所示。由于支撐塊式整體道床軌道結(jié)構(gòu)的預(yù)制混凝土支承塊(軌枕)與道床板通過(guò)混凝土澆筑在一起,因此可將其看作一個(gè)整體研究[15]。

        圖1 路基上整體道床軌道結(jié)構(gòu)

        無(wú)砟軌道三維實(shí)體單元模型能較好地反映無(wú)砟軌道的復(fù)雜空間力學(xué)特性[16]。因此,建立整體道床軌道結(jié)構(gòu)梁-體空間有限元實(shí)體模型,為了減小模型的計(jì)算量,建立有限元半尺寸模型,如圖2所示。模型中,將鋼軌視為彈性點(diǎn)支撐梁,扣件系統(tǒng)采用彈簧-阻尼單元模擬;考慮到道床板、支承層及路基在傳遞列車(chē)動(dòng)荷載過(guò)程中的重要性及其自身的彈塑性質(zhì),采用三維實(shí)體單元模擬,且只考慮其彈性行為。為消除邊界效應(yīng),模型長(zhǎng)度取60 m。整體道床軌道模型主要結(jié)構(gòu)參數(shù)如表1所示。

        圖2 城市軌道交通無(wú)砟軌道梁-體空間有限元模型

        結(jié)構(gòu)層厚度/m寬度mE/MPa泊松比密度/(kg·m-3)鋼軌——2.10×1050.307850道床板0.192.83.65×1040.207830底座板0.33.43.25×1040.202750路基0.46.01800.302300

        無(wú)砟軌道層間連接方式對(duì)模型計(jì)算的精確性起著重要作用。鋼軌與道床板間通過(guò)一系列離散支承的線性彈簧(模擬扣件)連接,為了考慮扣件的尺寸效應(yīng),將每個(gè)鋼軌節(jié)點(diǎn)與其對(duì)應(yīng)扣件尺寸范圍內(nèi)的道床板節(jié)點(diǎn)均相連[17];由于道床板與混凝土支承層之間連接緊密,且相對(duì)滑移較小,因此采用ABAQUS軟件中特有的Tie連接方式模擬接觸;考慮到路基表層出現(xiàn)不均勻沉降后可能產(chǎn)生的軌道離縫、脫空等缺陷,混凝土支承層與路基之間采用接觸單元模擬。由于列車(chē)動(dòng)荷載在無(wú)砟軌道各部件間形成的初始應(yīng)力對(duì)無(wú)砟軌道各部件間及無(wú)砟軌道與路基間接觸狀態(tài)非線性產(chǎn)生較大影響[18-19],因此在無(wú)砟軌道鋼軌上施加設(shè)計(jì)列車(chē)動(dòng)荷載,在路基上表面施加不均勻沉降荷載。

        1.2 路基沉降模型

        由于城市軌道交通線路路基的橫斷面一般都比較窄,且橫向不均勻沉降一般很少發(fā)生,因此,針對(duì)城市軌道交通線路路基不均勻沉降的研究中只考慮沿軌道縱向分布的垂向不均勻沉降,認(rèn)為沿橫向沉降分布均勻[20],據(jù)此,本文只考慮軌道結(jié)構(gòu)沿軌道縱向分布的垂向變形和縱向應(yīng)力。

        最常見(jiàn)的路基不均勻沉降模型主要包括3種典型形式。一種為余弦型不均勻沉降,一般最為常見(jiàn)且主要發(fā)生在路基上,如圖3所示。而在路橋過(guò)渡段或其他過(guò)渡段地區(qū),會(huì)發(fā)生錯(cuò)臺(tái)型或者折角型不均勻沉降[8],如圖4所示。而整體道床軌道結(jié)構(gòu)一般為分塊式結(jié)構(gòu),尤其是在結(jié)構(gòu)連接處必須留有一定的伸縮縫,因此在研究城市軌道交通線路路基不均勻沉降對(duì)軌道結(jié)構(gòu)受力變形的影響時(shí),將選用國(guó)內(nèi)使用最廣的下凹全波余弦型沉降曲線模型來(lái)模擬路基不均勻沉降,其描述函數(shù)如式(1)所示。

        (1)

        式中,y為沿線路縱向的位移,m;f0為沉降幅值,mm;l為沉降波長(zhǎng),m;Z為發(fā)生不均勻沉降的位置坐標(biāo)。

        圖3 余弦型不均勻沉降曲線

        圖4 錯(cuò)臺(tái)型和折角型沉降曲線

        1.3 荷載模擬

        當(dāng)路基發(fā)生不均勻沉降時(shí),引起無(wú)砟軌道結(jié)構(gòu)累積沉降變形的主要原因是無(wú)砟軌道自重荷載和列車(chē)動(dòng)荷載,因此在分析路基不均勻沉降引起無(wú)砟軌道軌面不平順的映射關(guān)系中同時(shí)考慮無(wú)砟軌道自重荷載和列車(chē)動(dòng)荷載的作用。其中,關(guān)于無(wú)砟軌道自重荷載和列車(chē)動(dòng)荷載的設(shè)計(jì)取值可參考文獻(xiàn)[21]。圖5為路基上發(fā)生余弦型不均勻沉降的示意,其中圖5(a)為不考慮軌道自重荷載和列車(chē)動(dòng)荷載,圖5(b)為考慮軌道自重荷載和列車(chē)動(dòng)荷載。

        2 路基沉降幅值影響

        2.1 沉降幅值對(duì)軌道結(jié)構(gòu)變形的影響

        表2列出了當(dāng)路基不均勻沉降波長(zhǎng)為20m時(shí),沉降幅值由5 mm等間距增加到45 mm時(shí)8種不同工況下軌道結(jié)構(gòu)最大跟隨性變形幅值。分析表2數(shù)據(jù)可知,在路基沉降波長(zhǎng)為20m條件下,鋼軌、道床板和支承層最大變形量均隨路基不均勻沉降幅值的增大呈增幅逐漸減小的增大趨勢(shì),且在同一路基沉降幅值條件下的最大變形量基本保持一致。以道床板為例,10~15 mm路基沉降幅值條件下道床板最大變形量增幅為4.07 mm,而當(dāng)路基沉降幅值為20~25 mm時(shí),道床板最大變形量增幅僅為1.83 mm。支承層與路基間的沉降差隨沉降幅值的增大而逐漸增大,進(jìn)而造成軌道脫空現(xiàn)象逐漸加劇。

        表2 不均勻沉降幅值對(duì)軌道結(jié)構(gòu)變形影響

        圖6為不同沉降幅值條件下鋼軌的跟隨性沉降變形曲線。由圖6可以看出,當(dāng)路基發(fā)生如圖3所示的余弦型不均勻沉降時(shí),鋼軌在軌道自重荷載和設(shè)計(jì)列車(chē)動(dòng)荷載的共同作用下發(fā)生跟隨性沉降變形,且在沉降的兩端邊緣區(qū)域出現(xiàn)了微小的局部上拱變形,使得鋼軌和道床板的沉降波長(zhǎng)相比路基沉降波長(zhǎng)出現(xiàn)了不同程度的擴(kuò)散,擴(kuò)散長(zhǎng)度隨路基沉降幅值增大而增大。此外,當(dāng)沉降幅值小于15 mm時(shí),鋼軌最大變形量與路基沉降幅值相差不大。當(dāng)沉降幅值超過(guò)15 mm后,隨著沉降幅值的進(jìn)一步增大,鋼軌的最大變形量增幅逐漸減小,與路基間沉降差值逐漸增大,此時(shí)出現(xiàn)圖5(b)所示的軌道脫空,且隨著路基沉降幅值的進(jìn)一步增大,軌道脫空將越來(lái)越嚴(yán)重。

        圖6 不同沉降幅值條件下鋼軌變形曲線

        2.2 沉降幅值對(duì)軌道結(jié)構(gòu)受力的影響

        表3列出了當(dāng)路基不均勻沉降波長(zhǎng)為20 m時(shí),沉降幅值由5 mm等間距增加到45 mm不同工況下鋼軌、道床板和支承層的縱向最大應(yīng)力值變化。圖7為20 m沉降波長(zhǎng)條件下鋼軌、道床板和支承層最大拉/壓應(yīng)力隨沉降幅值的變化曲線,其中圖7(a)為最大拉應(yīng)力,圖7(b)為最大壓應(yīng)力。

        表3 不均勻沉降幅值對(duì)軌道結(jié)構(gòu)應(yīng)力影響 MPa

        注:正值為拉應(yīng)力,負(fù)值為壓應(yīng)力。

        圖7 不同沉降幅值條件下軌道結(jié)構(gòu)應(yīng)力變化曲線

        由表3和圖7可知,當(dāng)路基沉降波長(zhǎng)一定時(shí),隨沉降幅值的增大,鋼軌、道床板和支承層的最大拉/壓應(yīng)力隨沉降幅值的增大呈增幅逐漸減小的增大趨勢(shì),當(dāng)沉降幅值小于20 mm時(shí),軌道各層結(jié)構(gòu)最大拉/壓應(yīng)力增幅較大,當(dāng)沉降幅值增大到20 mm之后時(shí),軌道各層結(jié)構(gòu)的最大拉/壓應(yīng)力增幅逐漸減小;相同沉降工況條件下,軌道各層結(jié)構(gòu)最大拉/壓應(yīng)力值從上往下依次減??;隨沉降幅值的進(jìn)一步增大,可能會(huì)導(dǎo)致軌道結(jié)構(gòu)的最大拉/壓壓力超出混凝土最大應(yīng)力極限,導(dǎo)致軌道結(jié)構(gòu)裂紋,進(jìn)而影響無(wú)砟軌道的服役壽命和安全性,所以在平時(shí)的養(yǎng)護(hù)維修作業(yè)中應(yīng)該嚴(yán)格控制路基不均勻沉降的幅值。

        3 路基沉降波長(zhǎng)影響

        3.1 沉降波長(zhǎng)對(duì)軌道結(jié)構(gòu)變形的影響

        表4列出了路基不均勻沉降幅值為20 mm條件下,沉降波長(zhǎng)由5 m等間距增加到45 m時(shí)8種不同工況下鋼軌、道床板和支承層的跟隨性沉降的最大變形量。由表4可知,在20 mm沉降幅值且相同沉降波長(zhǎng)條件下,鋼軌、道床板和支承層最大變形量依次增大,即鋼軌變形量最小,支承層變形量最大。由于鋼軌與道床板間通過(guò)彈簧阻尼元件連接,因此,道床板與鋼軌間存在微小的變形差異;由于道床板與支承層間是通過(guò)ABAQUS軟件中特有的Tie方式模擬接觸,因此,道床板與支承層變形量基本保持一致。

        表4 不均勻沉降波長(zhǎng)對(duì)軌道結(jié)構(gòu)變形影響

        圖8為20 mm沉降幅值不同沉降波長(zhǎng)條件下無(wú)砟軌道軌面變形曲線。由圖8可知,當(dāng)路基不均勻沉降幅值為20 mm時(shí),鋼軌最大變形量均隨沉降波長(zhǎng)的增大呈增幅逐漸減小的增大趨勢(shì),且最終趨于穩(wěn)定。當(dāng)路基不均勻沉降波長(zhǎng)小于25 m時(shí),此時(shí)有明顯的軌道脫空現(xiàn)象。如沉降波長(zhǎng)為15 m時(shí),軌道結(jié)構(gòu)最大變形量?jī)H為12.74 mm左右,遠(yuǎn)小于路基沉降幅值20 mm。當(dāng)沉降波長(zhǎng)增大到25 m之后,鋼軌最大變形量已接近路基沉降幅值,且基本保持穩(wěn)定,此時(shí)支承層與路基間脫空現(xiàn)象明顯緩解。因此,路基不均勻沉降幅值20 mm條件下沉降波長(zhǎng)小于25 m時(shí)對(duì)無(wú)砟軌道軌面不平順的影響顯著,應(yīng)引起足夠重視。

        圖8 不同沉降波長(zhǎng)條件下鋼軌變形曲線

        3.2 沉降波長(zhǎng)對(duì)軌道結(jié)構(gòu)受力的影響

        表5列出了當(dāng)路基不均勻沉降幅值為20 mm時(shí),沉降波長(zhǎng)由5 m等間距增加到45 m不同工況下鋼軌、道床板和支承層的縱向最大拉/壓應(yīng)力值變化,其中應(yīng)力為正表示軌道結(jié)構(gòu)在該位置處受拉力作用,應(yīng)力為負(fù)時(shí)表示軌道結(jié)構(gòu)在該位置處受壓力作用。圖9(a)和圖9(b)分別為20 mm沉降幅值不同沉降波長(zhǎng)工況下軌道結(jié)構(gòu)最大拉、壓應(yīng)力變化曲線。

        表5 不均勻沉降波長(zhǎng)對(duì)軌道結(jié)構(gòu)應(yīng)力影響 MPa

        注:正值為拉應(yīng)力,負(fù)值為壓應(yīng)力。

        從表5和圖9可以看出,鋼軌、道床板和支承層最大拉/壓應(yīng)力均隨沉降波長(zhǎng)的增大呈先增大后減小的變化規(guī)律。就鋼軌而言,當(dāng)路基不均勻沉降波長(zhǎng)達(dá)到25 m時(shí),其縱向最大拉/壓應(yīng)力達(dá)到最大值,路基不均勻沉降波長(zhǎng)大于25 m之后其拉/壓應(yīng)力開(kāi)始逐漸減小。對(duì)于道床板和支承層來(lái)說(shuō),當(dāng)沉降波長(zhǎng)達(dá)到25 m時(shí),縱向最大拉壓/應(yīng)力達(dá)到最大值,25 m波長(zhǎng)之后拉/壓應(yīng)力隨波長(zhǎng)的增大開(kāi)始逐漸減小。例如,以道床板為例,當(dāng)沉降波長(zhǎng)從5 m增加到25 m時(shí),最大拉、壓應(yīng)力增幅分別為4.4 MPa和3.82 MPa,而隨著路基不均勻沉降波長(zhǎng)繼續(xù)增大到45 m時(shí),當(dāng)波長(zhǎng)超過(guò)25 m之后道床板最大拉、壓應(yīng)力開(kāi)始逐漸減小。由此可見(jiàn),在路基沉降幅值20 mm條件下,沉降波長(zhǎng)為15~25 m時(shí)對(duì)軌道結(jié)構(gòu)應(yīng)力影響最為明顯,應(yīng)引起足夠重視。

        圖9 不同沉降波長(zhǎng)條件下軌道結(jié)構(gòu)應(yīng)力變化曲線

        4 軌面不平順與路基沉降之間的映射關(guān)系

        研究分析表明,當(dāng)路基發(fā)生圖3所示的余弦型不均勻沉降時(shí),靜平衡狀態(tài)下無(wú)砟軌道結(jié)構(gòu)在軌道自重和設(shè)計(jì)列車(chē)動(dòng)載的作用下會(huì)發(fā)生跟隨性沉降變形,并在沉降位置的兩端出現(xiàn)微小的局部上拱。由此可知,無(wú)砟軌道軌面不平順幾何變形曲線可由路基不均勻沉降的波長(zhǎng)、波深以及局部上拱幅值3項(xiàng)參數(shù)共同決定。因此,可以建立鋼軌軌面不平順與余弦型路基不均勻沉降之間的映射關(guān)系式[8]

        δr{λ,a,h}=δr{fλ(s,A),fa(s,A),fh(s,A)}

        (2)

        式中,δr為余弦型軌面不平順特征集;λ、a、h分別為軌面不平順變形曲線的波長(zhǎng)、幅值和局部上拱幅值;fλ、fa、fh分別為路基不均勻沉降對(duì)鋼軌變形波長(zhǎng)、幅值和上拱幅值的控制作用。

        5 結(jié)語(yǔ)

        針對(duì)城市軌道交通路基不均勻沉降問(wèn)題,通過(guò)建立考慮層間接觸非線性整體道床軌道梁-體空間有限元模型,分析了在軌道自重荷載和列車(chē)動(dòng)荷載作用下路基不均勻沉降對(duì)無(wú)砟軌道軌面不平順的影響規(guī)律,并在此基礎(chǔ)上,對(duì)有關(guān)城市軌道交通無(wú)砟軌道線路的路基不均勻沉降的安全限值進(jìn)行研究。主要得出如下結(jié)論。

        (1)無(wú)砟軌道軌面不平順受路基沉降幅值和波長(zhǎng)的綜合影響。一定沉降波長(zhǎng)條件下,沉降幅值越大軌面不平順越顯著,當(dāng)路基沉降幅值小于25 mm時(shí),軌道結(jié)構(gòu)變形與路基沉降的跟隨性較好,當(dāng)路基沉降幅值超過(guò)25 mm之后,軌道結(jié)構(gòu)與路基間形成明顯的脫空。

        (2)軌面不平順對(duì)路基沉降波長(zhǎng)也極為敏感。一定沉降幅值條件下,軌面不平順隨路基沉降波長(zhǎng)的增大呈先增大后穩(wěn)定,當(dāng)沉降波長(zhǎng)超過(guò)25 m時(shí)路基與軌道結(jié)構(gòu)間脫空現(xiàn)象明顯緩解,此時(shí)軌面不平順基本可與路基變形保持一致。

        (3)路基不均勻沉降導(dǎo)致軌道結(jié)構(gòu)在一定波長(zhǎng)范圍內(nèi)縱向最大應(yīng)力值總共出現(xiàn)3處,其中最大壓應(yīng)力出現(xiàn)在發(fā)生沉降的中心位置處,最大拉應(yīng)力出現(xiàn)在沉降位置的兩端;軌道結(jié)構(gòu)縱向拉/壓應(yīng)力隨沉降幅值的增大而增大,當(dāng)幅值超過(guò)25 mm之后拉/壓應(yīng)力增幅均逐漸減??;同一路基沉降條件下,軌道各層結(jié)構(gòu)最大拉/壓應(yīng)力值從上往下依次減小。

        (4)本文在研究路基沉降與軌面不平順間的映射關(guān)系中,著重考慮了無(wú)砟軌道結(jié)構(gòu)的自重荷載和設(shè)計(jì)列車(chē)動(dòng)荷載的影響,但未考慮列車(chē)循環(huán)荷載和周邊環(huán)境等復(fù)雜因素的綜合影響,還需在此基礎(chǔ)上結(jié)合現(xiàn)場(chǎng)觀測(cè)數(shù)據(jù)做進(jìn)一步的分析。研究結(jié)論將為下一步的系統(tǒng)動(dòng)力學(xué)特性研究提供理論基礎(chǔ)。

        猜你喜歡
        變形結(jié)構(gòu)
        《形而上學(xué)》△卷的結(jié)構(gòu)和位置
        談詩(shī)的變形
        論結(jié)構(gòu)
        新型平衡塊結(jié)構(gòu)的應(yīng)用
        模具制造(2019年3期)2019-06-06 02:10:54
        “我”的變形計(jì)
        變形巧算
        例談拼圖與整式變形
        會(huì)變形的餅
        論《日出》的結(jié)構(gòu)
        創(chuàng)新治理結(jié)構(gòu)促進(jìn)中小企業(yè)持續(xù)成長(zhǎng)
        国内免费AV网站在线观看| 国产精品永久久久久久久久久| 久久精品无码一区二区三区免费| 人妻去按摩店被黑人按中出| 亚洲无码性爱视频在线观看| 日本免费精品一区二区三区视频 | 亚洲男人的天堂在线播放| 韩国精品一区二区三区| 国产一区二区三区不卡在线播放 | 亚洲第一页在线免费观看| 男人天堂网2017| 成片免费观看视频大全| AV成人午夜无码一区二区| 中文字幕久久国产精品| 精品卡一卡二乱码新区| 99香蕉国产精品偷在线观看| 国产在线观看黄| 久久99国产精品久久99密桃| 免费人成在线观看视频高潮| 性一交一乱一伦a片| 国产亚洲精品hd网站| av免费在线国语对白| 日韩av无码精品一二三区| 国产欧美日产久久| 一本久久a久久精品综合| 国产网站一区二区三区| 久久综合九色综合欧美狠狠| 狠狠丁香激情久久综合| 日韩有码在线免费视频| a级毛片免费观看在线播放| 波多野结衣免费一区视频| 久久精品国产亚洲av麻豆四虎| 久久免费看的少妇一级特黄片| 日韩欧美亚洲综合久久影院ds| 中文字幕AⅤ人妻一区二区| 亚洲福利视频一区二区三区| 久久人人爽av亚洲精品| 乱码午夜-极国产极内射| 亚洲国产av一区二区三| 人妻少妇精品视频专区vr| 精品日韩欧美一区二区在线播放|