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

    基于非線性M-C準(zhǔn)則的圓形基礎(chǔ)下伏土洞抗沖切破壞研究

    2019-04-16 07:31:08王志斌羅謝鑫夏浩祥
    關(guān)鍵詞:土洞剛性圓形

    王志斌,羅謝鑫,夏浩祥

    ?

    基于非線性M-C準(zhǔn)則的圓形基礎(chǔ)下伏土洞抗沖切破壞研究

    王志斌,羅謝鑫,夏浩祥

    (湖南科技大學(xué) 土木工程學(xué)院,湖南 湘潭 411201)

    針對(duì)土體材料特點(diǎn),引入非線性Mohr-Coulomb準(zhǔn)則及相關(guān)流動(dòng)法則,構(gòu)建圓形剛性基礎(chǔ)下伏土洞沖切破壞模式。采用極限分析上限法推導(dǎo)出其極限承載力的計(jì)算公式,得出參數(shù)任意組合條件下圓形剛性基礎(chǔ)下伏土洞沖切極限承載力。利用Plaxis 3D軟件建模分析,并與理論解進(jìn)行對(duì)比,驗(yàn)證本文方法的有效性。研究結(jié)果表明:非線性系數(shù)對(duì)極限承載力的影響較大,在其他參數(shù)條件一定的情況下,極限承載力與呈非線性負(fù)相關(guān),且趨于穩(wěn)定。極限承載力隨土體初始黏聚力0,頂板厚度以及剛性基礎(chǔ)寬度的增大而增大,隨土體重度、土體抗拉強(qiáng)度的增大而減小。

    非線性Mohr-Coulomb準(zhǔn)則;極限分析;極限承載力;圓形基礎(chǔ);沖切

    我國(guó)巖溶地質(zhì)分布范圍廣泛,土洞是發(fā)育于可溶巖層上覆土體內(nèi)的空腔,是巖溶區(qū)常見(jiàn)的一種巖溶作用產(chǎn)物[1]。在一些沖積區(qū),第四系土層覆蓋厚度較大,可達(dá)40余米,土洞深埋于第四系土層中[2]。當(dāng)樁基或獨(dú)立基礎(chǔ)位于土洞或軟弱下臥層上方時(shí),常產(chǎn)生沖切破壞[3]。目前國(guó)內(nèi)外研究重點(diǎn)在于巖溶區(qū)樁基抗沖切破壞模式,Hongseob等[4]得出混凝土和巖石等脆性材料在受沖切作用破壞時(shí)的沖切破壞體并非標(biāo)準(zhǔn)的圓錐臺(tái)形,而是以一定曲線為母線所形成的軸對(duì)稱旋轉(zhuǎn)體。趙明華等[5?6]通過(guò)將溶洞簡(jiǎn)化為梁板結(jié)構(gòu)模型,系統(tǒng)地研究了溶洞頂板的穩(wěn)定性與安全厚度;張慧樂(lè)等[7?8]采用室內(nèi)模型試驗(yàn)的方法,研究了巖溶區(qū)嵌巖樁的承載能力與破壞模式,得出嵌巖樁破壞模式和影響因素;劉善軍等[9]基于應(yīng)力平衡法建立了松散蓋層穩(wěn)定性計(jì)算公式。曾中林等[10]采用變分法對(duì)巖溶區(qū)條形剛性基礎(chǔ)下伏孔洞巖體沖切破壞機(jī)制進(jìn)行了研究。雷勇等[11]采用變分法得到巖溶區(qū)樁端巖層抗沖切安全厚度計(jì)算公式。這些研究指明了地下空洞變形破壞的方向。通過(guò)對(duì)以上研究現(xiàn)狀的總結(jié)發(fā)現(xiàn),以往研究對(duì)象多集中在巖溶,以巖質(zhì)材料為主,對(duì)土質(zhì)材料形成的孔洞承受地表荷載的研究甚少。本文采用極限分析上限法,構(gòu)建圓形剛性基礎(chǔ)下伏土洞破壞模式,根據(jù)非線性Mohr-Coulomb準(zhǔn)則及相關(guān)流動(dòng)法則,推導(dǎo)出其極限承載力的計(jì)算公式,在此基礎(chǔ)上采用數(shù)學(xué)軟件Matlab進(jìn)行編程求解,得出非線性Mohr-Coulomb屈服準(zhǔn)則[12?17]參數(shù)任意組合條件下圓形剛性基礎(chǔ)下伏土洞沖切極限承載力,利用Plaxis 3D軟件建模分析,并與理論解作了對(duì)比,兩者吻合較好,驗(yàn)證了本文方法的有效性,具有較強(qiáng)的理論與實(shí)踐意義。

    1 基本理論

    極限分析上限定理:對(duì)于任意一個(gè)有效破壞機(jī)構(gòu),假定存在1個(gè)同時(shí)滿足速度邊界條件和應(yīng)變?速度相容條件的速度場(chǎng),根據(jù)外荷載功率等于內(nèi)部能量損耗率的關(guān)系條件,能夠得到1個(gè)不小于實(shí)際破壞荷載的極限荷載。具體如式(1)所示。

    2 剛性圓形基礎(chǔ)下伏土洞沖切上限分析

    2.1 水平圓形剛性基礎(chǔ)下伏土洞沖切破壞模型計(jì)算假定

    根據(jù)極限分析上限定理,對(duì)于任意假定的沖切破壞機(jī)制,其速度場(chǎng)必須滿足機(jī)動(dòng)容許的要求,即滿足如下條件:1) 速度邊界條件;2) 幾何相容條件;3) 相關(guān)聯(lián)流動(dòng)法則。通過(guò)虛功原理的功率平衡方程,解得的承載力為實(shí)際承載力的一個(gè)上限解。

    由于極限分析上限定理是建立在將土體假想為理想的彈塑性體和忽略了其幾何變形的基礎(chǔ)上,故在構(gòu)建破壞模式時(shí)作如下假設(shè):1) 圓形剛性基礎(chǔ)受沖切作用時(shí),基礎(chǔ)下覆土體為理想的彈塑性介質(zhì),破壞時(shí)服從非線性M-C強(qiáng)度屈服準(zhǔn)則及相關(guān)流動(dòng)法則;2) 基礎(chǔ)下部的破壞土體為剛性體,且破壞過(guò)程中不與剛性基礎(chǔ)發(fā)生相對(duì)滑動(dòng);3) 圓形剛性基礎(chǔ)直徑為2,下覆土層的厚度為,沖切破壞半徑為,下伏土洞頂部支護(hù)力為。

    為求得土洞區(qū)圓形剛性基礎(chǔ)下覆土體受沖切作用時(shí)的極限承載力,必須先假定一種沖切破壞機(jī)制。為了使獲得的上限解盡可能接近真實(shí)解,假定的沖切破壞機(jī)制還必須與圓形剛性基礎(chǔ)下覆土體的實(shí)際破壞模式相同。Hongseob等[4]得出混凝土和巖石等脆性材料在受沖切作用破壞時(shí)的沖切破壞體是以一定曲線為母線所形成的軸對(duì)稱旋轉(zhuǎn)體。這是由于混凝土與巖石等材料的強(qiáng)度屈服準(zhǔn)則往往具有非線性的特點(diǎn),材料破壞面上應(yīng)力分布不均,從而導(dǎo)致破壞面母線為曲線形式。

    根據(jù)極限分析上限定理和沖切破壞體自身的特點(diǎn),構(gòu)建了圖1所示的沖切破壞機(jī)制,并對(duì)其進(jìn)行了相應(yīng)假定:

    1) 圓形剛性基礎(chǔ)頂部受力的作用時(shí)傳遞到基礎(chǔ)底部的作用為圓形均布荷載,破壞時(shí)呈三維軸對(duì)稱破壞,且破壞場(chǎng)與基礎(chǔ)對(duì)稱軸重合。

    2) 沖切破壞面方程開(kāi)始于剛性基礎(chǔ)端部邊緣,在土體內(nèi)部連續(xù),其破壞面的母線幾何方程為()。

    3) 圓形剛性基礎(chǔ)下覆土體為理想均質(zhì)彈塑性體,破壞時(shí)服從非線性MC強(qiáng)度屈服準(zhǔn)則。

    4) 圓形剛性基礎(chǔ)在達(dá)到極限狀態(tài)時(shí),下覆土體內(nèi)部產(chǎn)生的塑性變形必須滿足相關(guān)聯(lián)流動(dòng)法則的要求。

    圖1 圓形剛性基礎(chǔ)三維沖切破壞機(jī)制

    2.2 水平圓形剛性基礎(chǔ)下伏土洞沖切極限承載力變分分析

    2.2.1 水平圓形剛性基礎(chǔ)下伏土洞沖切三維破壞機(jī)制的構(gòu)建

    假設(shè)圓形剛性基礎(chǔ)下覆土體為理想均質(zhì)彈塑性體,破壞時(shí)服從非線性Mohr-Coulomb屈服準(zhǔn)則,剛性基礎(chǔ)與土體之間無(wú)相對(duì)位移,圓形剛性基礎(chǔ)極限狀態(tài)下對(duì)應(yīng)的三維曲線型破壞機(jī)制構(gòu)建過(guò)程如下:首先假設(shè)有一條關(guān)于軸對(duì)稱的速度間斷線(),()自基礎(chǔ)底部開(kāi)始向下,直至土洞頂部,形成一個(gè)上部窄而下部寬的破壞面;其次令()繞軸旋轉(zhuǎn)180°,形成一個(gè)“漏斗”型的沖切破壞體;該沖切破壞體由半徑分別為和的圓和以()為母線的旋轉(zhuǎn)曲面組成。該沖切破壞體的體積和側(cè)表面積均可用解析數(shù)學(xué)的方法求解,在極限分析上限定理的能耗計(jì)算中,可以根據(jù)沖切破壞機(jī)構(gòu)的體積和側(cè)表面積算出其破壞機(jī)制的內(nèi)能耗散功率于外力功率。

    2.2.2 水平圓形剛性基礎(chǔ)下伏土洞沖切能耗計(jì)算

    將剛性基礎(chǔ)下覆土體沖切破壞速度間斷面視為薄變形層,其厚度為,作用在變形層上的剪應(yīng)力為τ,正應(yīng)力為σ,滿足非線性M-C強(qiáng)度屈服準(zhǔn)則。屈服函數(shù)可表示為:

    在主應(yīng)力空間內(nèi),令塑性勢(shì)函數(shù)面與強(qiáng)度屈服面重合,根據(jù)相關(guān)流動(dòng)法則,可得塑性勢(shì)函數(shù)為[10]:

    根據(jù)塑性位勢(shì)理論[10]:

    則將式(3)代入式(4)中可得分離薄層上的塑性應(yīng)變率為:

    根據(jù)圖1中幾何運(yùn)動(dòng)關(guān)系,可得:

    沖切破壞體分離薄層上的塑性應(yīng)變率也可表示為:

    將式(6)代入式(7)中可得:

    為滿足變形的相容協(xié)調(diào),式(5)和式(8)中的塑性應(yīng)變率分量應(yīng)相等,故可得分離薄層上的正應(yīng)力分量和剪應(yīng)力分量:

    沖切破壞時(shí)內(nèi)部能耗僅發(fā)生在剛性基礎(chǔ)下覆土體的破壞分離面處,故薄變形層上單位體積的耗散功率為:

    沖切破壞時(shí),破壞土體的側(cè)面積可以按旋轉(zhuǎn)體的側(cè)面積公式得到:

    將式(10)沿沖切破壞體側(cè)面積積分可得速度間斷面上內(nèi)能耗散功率為:

    沖切破壞時(shí)土體重力的功率為:

    圓形剛性基礎(chǔ)極限承載力的功率為:

    圓形剛性基礎(chǔ)下伏土洞洞頂支護(hù)力的功 率為:

    根據(jù)虛功率原理建立平衡方程,內(nèi)能耗散功率等于外力功率,故有:

    將式(12)~式(15)代入式(16)中可得圓形剛性基礎(chǔ)沖切極限承載力的表達(dá)式為:

    式中:

    由極限分析上限法可知,求解結(jié)果大于真實(shí)的極限承載力,其最小值才更接近真實(shí)解。極限承載力的極值由式(17)所確定,通過(guò)變分運(yùn)算可以將這一泛函的極值問(wèn)題轉(zhuǎn)化為求解歐拉方程在滿足邊界條件下的定解問(wèn)題。對(duì)應(yīng)的歐拉方程為[10]:

    根據(jù)歐拉方程對(duì)()和()的偏導(dǎo)分別為:

    將式(20)代入式(19)并整理可得:

    式中:1為積分常數(shù),可由如下邊界條件求得。

    由于沖切破壞時(shí)沒(méi)有水平外力作用,且沖切破壞為對(duì)稱破壞,故根據(jù)剛性基礎(chǔ)端部微元體的平衡可知,水平剪應(yīng)力分量τ=0。其表達(dá)式為:

    其中:

    故有:

    將式(23)代入式(22)并由式(21)可得:

    即:

    故有:

    由式(27)可知,只有當(dāng)?1是整數(shù)時(shí),才能積分得到()的精確表達(dá)式,根據(jù)位移邊界條件求得積分常數(shù)以及沖切破壞寬度,進(jìn)而求得剛性基礎(chǔ)下覆單層土體沖切極限承載力u的表達(dá)式。根據(jù)圖1可知,位移邊界條件如下:

    式中:為下伏土洞沖切破壞寬度。

    當(dāng)其不是整數(shù)時(shí),無(wú)法求得其原函數(shù),但可根據(jù)差分原理并結(jié)合邊界條件分層迭代得到()的函數(shù)曲線形狀,同理,也可得到固定參數(shù)條件下極限承載力u的數(shù)值。

    2.2.3 差分原理

    為獲求任意參數(shù)條件下圓形剛性基礎(chǔ)下伏土洞沖切極限承載力,以下利用差分原理進(jìn)行沖切破裂面曲線函數(shù)()數(shù)值求解,得到破裂面曲線函數(shù)()和旋轉(zhuǎn)后的三維曲面形狀,進(jìn)一步分層積分迭加可以求得對(duì)應(yīng)剛性基礎(chǔ)極限承載力u。

    根據(jù)差分法的基本概念:

    給定式中參數(shù),,σ,0和,可得出任意點(diǎn)x的導(dǎo)數(shù)值。假設(shè)曲線的每一微段為x+1?x=Δ,則將每點(diǎn)的導(dǎo)數(shù)值可視為點(diǎn)(x,(x))和點(diǎn)(x+1,(x+1))連接直線的斜率。將微小直線段近似為實(shí)際解析函數(shù)的曲線段,根據(jù)式(29)可得:

    利用式(30)進(jìn)行迭加計(jì)算,可得出每點(diǎn)的(x+1)值,其表達(dá)式為:

    由于無(wú)法得到精確的解析曲線表達(dá)式,極限承載力u的求解采取同樣的思路。首先得出每一微段直線對(duì)應(yīng)的薄層極限承載力u()的積分表達(dá)式:

    將全部的薄層極限承載力積分進(jìn)行累加,可得整個(gè)破壞機(jī)構(gòu)的極限承載力u的表達(dá)式為:

    根據(jù)以上推導(dǎo)過(guò)程,利用數(shù)學(xué)軟件Matlab進(jìn)行編程求解,可得出非線性Mohr-Coulomb屈服準(zhǔn)則參數(shù)任意組合條件下圓形剛性基礎(chǔ)下伏巖溶土洞沖切極限承載力及其對(duì)應(yīng)三維破壞模式。

    2.3 圓形剛性基礎(chǔ)下覆單層土體特殊分析模型

    當(dāng)=1.0時(shí),非線性MC準(zhǔn)則退化為線性Mohr- Coulomb準(zhǔn)則,沖切破壞面退化為理想的圓臺(tái)面,如圖2所示。

    由圖2可知,破壞直線的斜率為:

    故其破壞直線的函數(shù)形式為:

    式中:為待定系數(shù),可由如下邊界條件確定。

    代入可得:

    破壞直線方程為:

    此時(shí)沖切破壞時(shí)內(nèi)能耗散功率為:

    圖2 圓形剛性基礎(chǔ)三維沖切直線破壞機(jī)制

    土體重力的功率為:

    式中:表示破壞圓臺(tái)的體積。

    根據(jù)虛功率原理建立平衡方程可得:

    3 對(duì)比分析

    采用Plaxis 3D建模進(jìn)行對(duì)比分析。在半無(wú)限空間中,沖切破壞為圓對(duì)稱破壞,故取1/4模型進(jìn)行分析,如圖3所示。模型邊界尺寸20 m×20 m×20 m,圓形剛性基礎(chǔ)半徑=1 m,其上作用有均布荷載,即未知的極限承載力,土層厚度=5 m,孔洞半徑=4.5 m,洞頂支護(hù)力=20 kN/m2,土體容重=18 kN/m3,土體黏聚力為=20 kPa,單軸抗拉強(qiáng)度t=30 kPa,非線性系數(shù)=1.0,內(nèi)摩擦角=33.7°,楊氏模量=10 000 kN/m2,泊松比=0.3,不考慮水的影響。

    圖3 1/4模型圖

    利用Plaxis 3D軟件進(jìn)行塑性分析和安全性分析,通過(guò)在剛性基礎(chǔ)上施加荷載,計(jì)算模型的安全系數(shù),當(dāng)計(jì)算得到的安全系數(shù)恰好為1時(shí),此時(shí)的荷載為極限荷載。通過(guò)計(jì)算,得到的沖切極限承載力約為177 kN,安全系數(shù)s=1.002,其塑性剪切帶三維圖如圖4所示。而本文基于極限分析上限定理所得沖切極限承載力為185.5 kN,孔洞破壞寬度約為4.33 m,其極限承載力和破壞寬度的誤差分別為4.58%和3.93%,均不超過(guò)5%。同時(shí),將剪切破壞帶與理論解的破壞曲線作了對(duì)比,并將其展示在同一張圖中,如圖5所示,兩者吻合較好,由此可說(shuō)明本文方法的可行性和有效性。

    圖4 塑性剪切帶

    圖5 三維破壞曲面剖面圖

    4 參數(shù)與破壞模式分析

    由式(17)可知,影響圓形剛性基礎(chǔ)土洞頂板的沖切極限承載力的因素眾多。且非線性M-C準(zhǔn)則的非線性系數(shù)對(duì)極限承載力的影響更為顯著,以下分析土洞頂板自由沖切極限承載力隨單一參數(shù)的影響變化情況以及破壞模式的變化。取剛性基礎(chǔ)半徑=0.5~1.5 m,土體抗拉強(qiáng)度t=40~60 kPa,土體初始黏聚力0=10~30 kPa,土體重度=16~24 kN/m3,頂板厚度=4~8 m,土洞頂板底部支護(hù)力=100 kPa,分別繪制出極限承載力隨=1.0~3.0之間變化時(shí)的影響曲線。如圖6~10所示。

    圖6 Pu隨C0變化曲線

    綜合圖6~10可知,非線性系數(shù)對(duì)極限承載力的影響較大,在其他參數(shù)條件一定的情況下,極限承載力與呈非線性負(fù)相關(guān),且趨于穩(wěn)定。極限承載力隨土體初始黏聚力0,頂板厚度以及剛性基礎(chǔ)寬度的增大而增大,隨土體重度、土體抗拉強(qiáng)度的增大而減小。

    圖7 Pu隨d變化曲線

    圖8 Pu隨γ變化曲線

    圖9 Pu隨σt變化曲線

    圖10 Pu隨H變化曲線

    5 結(jié)論

    1) 根據(jù)非線性Mohr-Coulomb準(zhǔn)則及相關(guān)流動(dòng)法則,構(gòu)造出圓形基礎(chǔ)下伏土洞土體沖切破壞的機(jī)動(dòng)速度場(chǎng),利用極限分析上限法結(jié)合變分原理,推導(dǎo)出圓形剛性基礎(chǔ)下伏土洞沖切破壞時(shí)的破壞曲面函數(shù),利用差分原理結(jié)合邊界條件迭代出破壞曲面函數(shù)形狀及其極限承載力,拓展分析了圓形剛性基礎(chǔ)下伏土洞沖切破壞時(shí)的特殊模型。利用Plaxis 3D軟件建模分析,并與理論解作了對(duì)比,兩者吻合較好,驗(yàn)證了本文方法的有效性。

    2) 對(duì)于圓形剛性基礎(chǔ)下伏土洞的沖切破壞,非線性系數(shù)對(duì)極限承載力的影響較大,在其他參數(shù)條件一定的情況下,極限承載力與呈非線性負(fù)相關(guān),且趨于穩(wěn)定。極限承載力隨土體初始黏聚力、頂板厚度以及剛性基礎(chǔ)寬度的增大而增大,隨土體重度、土體抗拉強(qiáng)度的增大而減小。

    [1] 劉鎮(zhèn), 黎杰明, 楊旭, 等. 考慮鄰近土洞影響的盾構(gòu)掘進(jìn)速度控制[J]. 工程地質(zhì)學(xué)報(bào), 2017, 25(6): 1633? 1639.LIU Zhen, LI Jieming, YANG Xu, et al. Optimal control of TBM excavation rate considering neighboring soil void[J]. Journal of Engineering Geology, 2017, 25(6): 1633?1639.

    [2] 丁春林, 甘百先, 鐘輝虹, 等. 含土洞、溶洞的機(jī)場(chǎng)滑行道路基穩(wěn)定性評(píng)估[J]. 巖石力學(xué)與工程學(xué)報(bào), 2003, 22(8): 1329?1333. DING Chunlin, GAN Baixian, ZHONG Huihong, et al. Stability evaluation of airfield runway subgrade containing earth caves and karst[J]. Chinese Journal of Rock Mechanics and Engineering, 2003, 22(8): 1329? 1333.

    [3] Wyllie D C. Foundations on rock[M]. London: Chapman and Hall, 1992.

    [4] Hongseob O H, Jongsung S I M. Punching shear strength of strengthened deck panels with externally bonded plates[J]. Composites part B: Engineering, 2004, 35(4): 313?321.

    [5] 趙明華, 張銳, 胡柏學(xué), 等. 巖溶區(qū)樁端下伏溶洞頂板穩(wěn)定性分析研究[J]. 公路交通科技, 2009, 26(9): 13?16. ZHAO Minghua, ZHANG Rui, HU Boxue, et al. Analysis of stability of cave roof under pile tip in karst area[J]. Journal of Highway and Transportation Research and Development, 2009, 26(9): 13?16.

    [6] 趙明華, 雷勇, 張銳.巖溶區(qū)樁基沖切破壞模式及安全厚度研究[J]. 巖土力學(xué), 2012, 33(2): 524?530.ZHAO Minghua, LEI Yong, ZHANG Rui. Study of punching failure mode and safe thickness of pile foundation in karst region[J]. Rock and Soil Mechanics, 2012, 33(2): 524?530.

    [7] 張惠樂(lè), 張智浩, 王述紅, 等. 巖溶區(qū)嵌巖樁的試驗(yàn)研究與分析[J]. 土木工程學(xué)報(bào), 2013, 46(1): 92?103. ZHANG Huile, ZHANG Zhihao, WANG Shuhong, et al. Experimental study and analysis on rock-socketed pile in karst area[J]. China Civil Engineering Journal, 2013, 46(1): 92?103.

    [8] 黃生根, 梅世龍. 南盤(pán)江特大橋巖溶樁基承載特性的試驗(yàn)研究[J]. 巖石力學(xué)與工程學(xué)報(bào), 2004, 23(5): 809? 813. HUANG Shenggen, MEI Shilong. Testing study on bearing behavior of piles for Nanpan River great bridge in karst area[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(5): 809?813.

    [9] 劉善軍, 宋長(zhǎng)斌. 隱伏碳酸鹽巖區(qū)第四系松散蓋層穩(wěn)定性分析[J]. 中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào), 2003, 14(1): 18?21.LIU Shanjun, SONG Changbin. Stability analysis of Quaternary loose cover in carbonate rock region[J]. The Chinese Journal of Geological Hazard and Control, 2003, 14(1): 18?21.

    [10] 曾中林, 曾大艇, 李韜宇, 等. 巖溶區(qū)條形剛性基礎(chǔ)下伏孔洞頂板沖切極限承載力變分分析[J]. 自然災(zāi)害學(xué)報(bào), 2018, 27(5): 27?38.ZENG Zhonglin, ZENG Dating, LI Taoyu, et al. Variation analysis of ultimate bearing punching capacity of roof underlying the rigid strip foundation in karst area[J]. Journal of Natural Disasters, 2018, 27(5): 27?38.

    [11] 雷勇, 陳秋南, 馬繽輝.基于極限分析的樁端巖層沖切分析[J]. 巖石力學(xué)與工程學(xué)報(bào), 2014, 33(3): 631?638.LEI Yong, CHEN Qiunan, MA Binhui. Punching analysis of rock at pile tip base on limit analysis[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(3): 631?638.

    [12] Maksimovic M. Nonlinear failure envelope for soils[J]. Journal of Geotechnical and Geoenvironmental Engineering, 1989, 115(4): 581?586.

    [13] Baker R. Nonlinear Mohr envelopes based on triaxial data[J]. Journal of Geotechnical and Geoenvironmental Engineering, ASCE, 2004, 130(5): 498?506.

    [14] Hoek E, Brown E T. Practical estimate the rock mass strength[J]. International Journal of Rock Mechanics and Mining Science, 1997, 34(8): 1165?1186.

    [15] YANG X L, YIN J H. Slope stability analysis with nonlinear failure criterion[J]. Journal of Engineering Mechanics, ASCE, 2004, 130(3): 267?273.

    [16] ZHAO L H, LI L, YANG F, et al. Upper bound analysis of slope stability with nonlinear failure criterion based on strength reduction technique[J]. Journal of Central South University of Technology, 2010, 17(4): 836?844.

    [17] ZHAO Lianheng, YANG Feng, ZHANG Yingbin, et al. Effects of shear strength reduction strategies on safety factor of homogeneous slope based on a general nonlinear failure criterion[J]. Computers & Geotechnics, 2015, 63: 215?228.

    [18] 譚亦高, 左仕, 胡世紅, 等. 三維圓形淺埋錨板極限抗拔力非線性上限解析[J]. 鐵道科學(xué)與工程學(xué)報(bào), 2017, 14(6): 1166?1173.TAN Yigao, ZUO Shi, HU Shihong, et al. Study on the ultimate pullout capacity of 3-D shallow circle plate anchors[J]. Journal of Railway Science and Engineering, 2017, 14(6): 1166?1173.

    [19] 張佳華, 張標(biāo). 非線性屈服準(zhǔn)則下淺埋隧道支護(hù)參數(shù)上限解[J]. 鐵道科學(xué)與工程學(xué)報(bào), 2018, 15(8): 2061? 2071.ZHANG Jiahua, ZHANG Biao. Upper bound solutions of support parameter for shallow tunnels in nonlinear soil masses[J]. Journal of Railway Science and Engineering, 2018, 15(8): 2061?2071.

    Punching analysis of soil under circular foundation based on nonlinear M-C criterion

    WANG Zhibin, LUO Xiexin, XIA Haoxiang

    (School of Civil Engineering, Hunan University of Science and Technology, Xiangtan 411201, China)

    The nonlinear Mohr-Coulomb failure criterion and the associated flow rule have been introduced according to the characteristics of soil mass, and then the punching failure mechanism of the soil with an underlying circular rigid foundation was constructed. Then the calculation formula of the ultimate bearing capacity was established using limit analysis method. Furthermore, the ultimate bearing capacity of the soil with an underlying circular rigid foundation was obtained in the condition of arbitrary combination of parameters. The Plaxis 3D software was used to analyze, and the numerical results were in agreement with the theoretical solutions, which verified the validity of the proposed method. The influence of parameters on ultimate bearing capacity was also analyzed, and the results shows that the influence of the nonlinear coefficienton the ultimate bearing capacity is obvious, and that the ultimate bearing capacity decreases with the increase of, and the decreasing trend tends to be gentle with given parameters,. The ultimate bearing capacity increases with the0of initial cohesive force of soil, the roof thickness and the rigidity foundation width, and decreases with the increase of soil weight and the increase of soil tensile strength.

    nonlinear Mohr-Coulomb failure criterion; limit analysis; ultimate bearing capacity; circular foundation; punching

    10.19713/j.cnki.43?1423/u.2019.03.011

    TU473

    A

    1672 ? 7029(2019)03 ? 0637 ? 09

    2018?04?01

    國(guó)家重點(diǎn)研發(fā)計(jì)劃資助項(xiàng)目(2017YFB1201204)

    王志斌(1975?),男,湖南桃江人,講師,博士,從事邊坡穩(wěn)定性及隧道圍巖穩(wěn)定性等研究;E?mail:71496933@qq.com

    (編輯 涂鵬)

    猜你喜歡
    土洞剛性圓形
    自我革命需要“剛性推進(jìn)”
    加權(quán)p-Laplace型方程的剛性
    下伏土洞加筋地基條形荷載下應(yīng)力擴(kuò)散計(jì)算
    剛性兌付的法律治理
    金融法苑(2018年2期)2018-12-07 00:59:52
    基于Schwarz交替法的巖溶區(qū)雙孔土洞地基穩(wěn)定性分析
    為什么窨井蓋大多都是圓形的
    基于數(shù)值模擬的土洞穩(wěn)定性分析
    肥皂泡為什么是圓形?
    土洞施工中含水段塌方處理方案探討
    圓形題
    国产精品熟女久久久久浪| 观看av在线不卡| 免费av不卡在线播放| 久久久久国产网址| 国产精品国产av在线观看| 91成人精品电影| 人妻系列 视频| av女优亚洲男人天堂| 精品国产国语对白av| 视频中文字幕在线观看| 国产男女内射视频| 人人妻人人澡人人看| 又大又黄又爽视频免费| 久久久久久久久久久久大奶| 美女脱内裤让男人舔精品视频| 精品少妇黑人巨大在线播放| 欧美激情国产日韩精品一区| 日韩欧美一区视频在线观看 | 久久久久网色| 美女福利国产在线| 99九九线精品视频在线观看视频| 各种免费的搞黄视频| 国产av精品麻豆| 九草在线视频观看| 亚洲欧美日韩东京热| 一个人免费看片子| 欧美日韩亚洲高清精品| 国产日韩欧美视频二区| 精品人妻熟女av久视频| 久久久久精品性色| 下体分泌物呈黄色| 免费观看无遮挡的男女| 18禁在线播放成人免费| 99久久精品热视频| 国产精品无大码| 天堂中文最新版在线下载| 精品人妻熟女毛片av久久网站| 国产免费又黄又爽又色| 在线亚洲精品国产二区图片欧美 | 少妇被粗大猛烈的视频| 两个人的视频大全免费| 人妻 亚洲 视频| av在线app专区| 亚洲av综合色区一区| 精品国产国语对白av| 国产黄片视频在线免费观看| 国产成人精品婷婷| 成人漫画全彩无遮挡| 18禁在线播放成人免费| 一区二区三区精品91| 人妻夜夜爽99麻豆av| 一个人看视频在线观看www免费| 97在线人人人人妻| 国产日韩欧美视频二区| 麻豆精品久久久久久蜜桃| 久久国内精品自在自线图片| 国产精品熟女久久久久浪| 国产老妇伦熟女老妇高清| 成人特级av手机在线观看| 又大又黄又爽视频免费| 又大又黄又爽视频免费| 日日啪夜夜撸| 国产成人精品无人区| 岛国毛片在线播放| 另类亚洲欧美激情| 国产精品久久久久久久久免| 久久99热6这里只有精品| 日韩一区二区三区影片| 国产 精品1| 国产成人精品婷婷| 亚州av有码| 久久免费观看电影| 内地一区二区视频在线| 丰满少妇做爰视频| 高清av免费在线| 伊人亚洲综合成人网| 国产成人一区二区在线| 少妇熟女欧美另类| 9色porny在线观看| 婷婷色综合大香蕉| 在线看a的网站| 人妻制服诱惑在线中文字幕| 欧美人与善性xxx| 欧美激情国产日韩精品一区| 亚洲欧美日韩东京热| 三级经典国产精品| 午夜福利在线观看免费完整高清在| 精品久久久噜噜| 国产精品久久久久久精品古装| 在线免费观看不下载黄p国产| 亚洲成人一二三区av| 亚洲美女搞黄在线观看| 亚洲图色成人| 熟妇人妻不卡中文字幕| av在线app专区| 中文字幕制服av| 99热6这里只有精品| 国产高清有码在线观看视频| 男男h啪啪无遮挡| 一区二区三区免费毛片| 下体分泌物呈黄色| 97精品久久久久久久久久精品| 国产美女午夜福利| 国产中年淑女户外野战色| 2021少妇久久久久久久久久久| 综合色丁香网| 精品国产国语对白av| a 毛片基地| 亚洲av男天堂| av网站免费在线观看视频| 全区人妻精品视频| 亚洲av日韩在线播放| 午夜av观看不卡| 色吧在线观看| a级毛片免费高清观看在线播放| 日韩一本色道免费dvd| 99国产精品免费福利视频| 国产爽快片一区二区三区| 国产成人午夜福利电影在线观看| 建设人人有责人人尽责人人享有的| 王馨瑶露胸无遮挡在线观看| 精品久久久精品久久久| 国产淫片久久久久久久久| 久久人人爽人人片av| 妹子高潮喷水视频| 国模一区二区三区四区视频| 人妻人人澡人人爽人人| 性色avwww在线观看| 我的女老师完整版在线观看| 国产极品粉嫩免费观看在线 | 综合色丁香网| 热re99久久国产66热| 亚洲四区av| 久久久久久久大尺度免费视频| 亚洲av成人精品一区久久| 日韩视频在线欧美| 18+在线观看网站| 久久狼人影院| 国产淫语在线视频| 黄色欧美视频在线观看| 纯流量卡能插随身wifi吗| 国产无遮挡羞羞视频在线观看| 美女内射精品一级片tv| 嫩草影院入口| 91在线精品国自产拍蜜月| 久久综合国产亚洲精品| 男女边吃奶边做爰视频| 七月丁香在线播放| 国产伦精品一区二区三区视频9| 亚洲精品日韩av片在线观看| 亚洲av.av天堂| 最新中文字幕久久久久| 91aial.com中文字幕在线观看| 这个男人来自地球电影免费观看 | 成年女人在线观看亚洲视频| 又黄又爽又刺激的免费视频.| 久久久久久久久久久丰满| 国产在线一区二区三区精| 久久久a久久爽久久v久久| 男女无遮挡免费网站观看| 日韩欧美一区视频在线观看 | 九九久久精品国产亚洲av麻豆| 精品久久久久久电影网| 十分钟在线观看高清视频www | 大片免费播放器 马上看| 午夜老司机福利剧场| 国产精品99久久99久久久不卡 | 久久久国产精品麻豆| 欧美+日韩+精品| 成人毛片60女人毛片免费| 国产精品一区二区在线观看99| av福利片在线观看| 国产熟女欧美一区二区| 亚洲国产日韩一区二区| 亚洲欧美成人精品一区二区| 亚洲性久久影院| 永久网站在线| 国产一区二区三区av在线| 亚洲综合色惰| 日本猛色少妇xxxxx猛交久久| 久久人人爽av亚洲精品天堂| 大码成人一级视频| 国内少妇人妻偷人精品xxx网站| 日本黄色日本黄色录像| 免费av不卡在线播放| 日韩熟女老妇一区二区性免费视频| 久久女婷五月综合色啪小说| 三级经典国产精品| 国产男女超爽视频在线观看| av网站免费在线观看视频| 午夜福利网站1000一区二区三区| 精品少妇黑人巨大在线播放| 王馨瑶露胸无遮挡在线观看| 国产精品女同一区二区软件| 丰满迷人的少妇在线观看| 国产男女内射视频| 国产伦精品一区二区三区视频9| 免费观看a级毛片全部| 国产视频首页在线观看| 精品卡一卡二卡四卡免费| 精品人妻熟女av久视频| 秋霞在线观看毛片| 日日啪夜夜爽| 日韩精品免费视频一区二区三区 | 青青草视频在线视频观看| 亚洲va在线va天堂va国产| 国产精品一二三区在线看| 国产日韩欧美亚洲二区| 亚洲精品一二三| 亚洲久久久国产精品| 波野结衣二区三区在线| 边亲边吃奶的免费视频| 国内揄拍国产精品人妻在线| 国产精品秋霞免费鲁丝片| 国产男人的电影天堂91| 青春草亚洲视频在线观看| 久久午夜福利片| 黄色毛片三级朝国网站 | 久久久精品免费免费高清| 十八禁网站网址无遮挡 | 亚洲一区二区三区欧美精品| 91午夜精品亚洲一区二区三区| 少妇人妻一区二区三区视频| 男女啪啪激烈高潮av片| 在线精品无人区一区二区三| 在线观看人妻少妇| 免费人成在线观看视频色| 97在线人人人人妻| 99久久精品一区二区三区| 在线精品无人区一区二区三| 夜夜看夜夜爽夜夜摸| 国内少妇人妻偷人精品xxx网站| 国产精品福利在线免费观看| 九九在线视频观看精品| 最近中文字幕2019免费版| 男女边摸边吃奶| 街头女战士在线观看网站| 丰满人妻一区二区三区视频av| 99热网站在线观看| 爱豆传媒免费全集在线观看| 大码成人一级视频| 精品一区在线观看国产| 久久99热这里只频精品6学生| 午夜福利视频精品| 精品少妇内射三级| 国产精品久久久久久精品电影小说| 特大巨黑吊av在线直播| 一级爰片在线观看| 蜜桃在线观看..| 成人亚洲精品一区在线观看| 精品国产一区二区三区久久久樱花| 久久99一区二区三区| 亚洲av中文av极速乱| 国产一区二区三区av在线| 人人妻人人澡人人看| 精品国产一区二区久久| 久久久久久久久大av| 成人漫画全彩无遮挡| 免费观看性生交大片5| 精品卡一卡二卡四卡免费| 久久综合国产亚洲精品| 激情五月婷婷亚洲| 国国产精品蜜臀av免费| 男女国产视频网站| 国产精品国产三级国产av玫瑰| 久久久久精品久久久久真实原创| 一区二区三区精品91| 亚洲经典国产精华液单| 国产高清不卡午夜福利| 国产精品一区二区在线观看99| 国产高清国产精品国产三级| 亚洲欧洲日产国产| 日日爽夜夜爽网站| 国产美女午夜福利| 久久久a久久爽久久v久久| 久久久久久久久久人人人人人人| 亚洲天堂av无毛| 日韩三级伦理在线观看| 能在线免费看毛片的网站| 国产视频内射| 国产淫语在线视频| 久久国产精品大桥未久av | 亚洲高清免费不卡视频| 国产69精品久久久久777片| 国产av一区二区精品久久| 日本黄大片高清| 亚洲国产精品专区欧美| 一本大道久久a久久精品| 国产日韩欧美亚洲二区| 国产av精品麻豆| 亚洲精品色激情综合| 中文字幕av电影在线播放| 欧美丝袜亚洲另类| 欧美97在线视频| 老司机影院毛片| 国产熟女午夜一区二区三区 | 精品久久久精品久久久| 免费观看无遮挡的男女| 亚洲国产精品国产精品| 国产成人精品久久久久久| 欧美另类一区| 久热久热在线精品观看| 成人亚洲欧美一区二区av| 婷婷色麻豆天堂久久| 亚洲欧洲精品一区二区精品久久久 | 国产欧美日韩综合在线一区二区 | 日韩视频在线欧美| 精品酒店卫生间| 91午夜精品亚洲一区二区三区| 啦啦啦在线观看免费高清www| 高清欧美精品videossex| 高清在线视频一区二区三区| 大香蕉久久网| 精品国产乱码久久久久久小说| 免费大片18禁| 久久99热6这里只有精品| 国产精品不卡视频一区二区| 亚洲精品亚洲一区二区| 成人特级av手机在线观看| 国产精品一区www在线观看| 久久国产精品男人的天堂亚洲 | 国产高清有码在线观看视频| 欧美丝袜亚洲另类| 国产成人午夜福利电影在线观看| 肉色欧美久久久久久久蜜桃| 亚洲欧洲国产日韩| 人人妻人人添人人爽欧美一区卜| 麻豆精品久久久久久蜜桃| 免费人成在线观看视频色| 久久99一区二区三区| 亚洲国产欧美在线一区| 在现免费观看毛片| 男人添女人高潮全过程视频| 在线看a的网站| 亚洲av福利一区| av女优亚洲男人天堂| 亚洲第一av免费看| 亚洲欧美精品自产自拍| 亚洲国产欧美在线一区| 国产片特级美女逼逼视频| 久久国产精品男人的天堂亚洲 | 曰老女人黄片| 免费少妇av软件| 久久人妻熟女aⅴ| 亚洲av男天堂| 春色校园在线视频观看| 亚洲国产毛片av蜜桃av| 国产精品一区二区性色av| 亚洲国产欧美在线一区| 老司机亚洲免费影院| 亚洲第一av免费看| 三级国产精品欧美在线观看| 少妇的逼好多水| a 毛片基地| 日韩电影二区| 午夜福利网站1000一区二区三区| 午夜激情久久久久久久| 曰老女人黄片| a级一级毛片免费在线观看| 国内少妇人妻偷人精品xxx网站| 国产精品三级大全| 51国产日韩欧美| 熟女av电影| 亚洲精品,欧美精品| 人人妻人人添人人爽欧美一区卜| 香蕉精品网在线| 久久av网站| 在现免费观看毛片| av免费观看日本| 日韩熟女老妇一区二区性免费视频| 国产亚洲5aaaaa淫片| 在线天堂最新版资源| 3wmmmm亚洲av在线观看| 亚洲欧洲日产国产| 日韩制服骚丝袜av| 视频中文字幕在线观看| 午夜精品国产一区二区电影| 欧美+日韩+精品| 汤姆久久久久久久影院中文字幕| 丰满乱子伦码专区| 一级av片app| 国产有黄有色有爽视频| 日韩精品有码人妻一区| 能在线免费看毛片的网站| a级一级毛片免费在线观看| 国产伦精品一区二区三区视频9| 人妻少妇偷人精品九色| 日韩,欧美,国产一区二区三区| 女人精品久久久久毛片| 国产 精品1| 精品酒店卫生间| 免费不卡的大黄色大毛片视频在线观看| 午夜激情福利司机影院| 三级经典国产精品| 一边亲一边摸免费视频| 色视频www国产| 中文字幕人妻丝袜制服| 狂野欧美激情性xxxx在线观看| 国产精品蜜桃在线观看| 欧美激情极品国产一区二区三区 | 欧美性感艳星| 少妇熟女欧美另类| 日本黄色日本黄色录像| 99久久精品热视频| 亚洲成人手机| 成人美女网站在线观看视频| 中文天堂在线官网| 99九九线精品视频在线观看视频| 如何舔出高潮| 国产精品99久久久久久久久| 视频区图区小说| 99九九在线精品视频 | 国产午夜精品久久久久久一区二区三区| 91在线精品国自产拍蜜月| 亚洲成人一二三区av| 大陆偷拍与自拍| 大片电影免费在线观看免费| 97精品久久久久久久久久精品| 精品国产一区二区三区久久久樱花| av国产精品久久久久影院| 在线播放无遮挡| 一级片'在线观看视频| 日本91视频免费播放| 最新中文字幕久久久久| 特大巨黑吊av在线直播| 国产亚洲5aaaaa淫片| 国产美女午夜福利| 国产成人精品无人区| 日韩免费高清中文字幕av| 成人亚洲精品一区在线观看| av在线播放精品| 国产伦在线观看视频一区| 国产视频首页在线观看| 中文字幕久久专区| 中国三级夫妇交换| 亚洲精品日本国产第一区| 国产精品99久久久久久久久| 一二三四中文在线观看免费高清| av天堂中文字幕网| 国产精品久久久久成人av| 欧美三级亚洲精品| 下体分泌物呈黄色| 久久久午夜欧美精品| 国产乱来视频区| 自拍偷自拍亚洲精品老妇| 国产无遮挡羞羞视频在线观看| 欧美日韩在线观看h| 嘟嘟电影网在线观看| 少妇被粗大猛烈的视频| 亚洲av在线观看美女高潮| 精品久久久久久电影网| 成人美女网站在线观看视频| 伦理电影免费视频| 亚洲在久久综合| 国产成人a∨麻豆精品| 在线天堂最新版资源| 国产亚洲精品久久久com| 国产免费一级a男人的天堂| 在线免费观看不下载黄p国产| 日本与韩国留学比较| 纯流量卡能插随身wifi吗| 亚洲综合色惰| 日韩一区二区视频免费看| 丰满迷人的少妇在线观看| 大香蕉久久网| 亚洲四区av| 亚洲一区二区三区欧美精品| 久久久久久久国产电影| 免费不卡的大黄色大毛片视频在线观看| h视频一区二区三区| 欧美性感艳星| 成人亚洲欧美一区二区av| 国产精品熟女久久久久浪| 亚洲精品乱码久久久久久按摩| 亚洲av综合色区一区| 亚洲av成人精品一二三区| 国产免费又黄又爽又色| 狂野欧美激情性xxxx在线观看| 九草在线视频观看| 久久国内精品自在自线图片| 久久热精品热| 91精品国产九色| 色吧在线观看| 另类亚洲欧美激情| 欧美性感艳星| 熟女电影av网| 久久精品夜色国产| 亚洲美女视频黄频| 成年女人在线观看亚洲视频| 一级毛片久久久久久久久女| 色网站视频免费| 国产精品一区www在线观看| 国产精品偷伦视频观看了| 高清毛片免费看| 中文精品一卡2卡3卡4更新| 蜜桃在线观看..| 啦啦啦视频在线资源免费观看| 国产有黄有色有爽视频| 免费黄色在线免费观看| 国产黄频视频在线观看| 亚洲精品乱码久久久v下载方式| 亚洲av福利一区| 久久av网站| 在线观看免费高清a一片| 欧美日本中文国产一区发布| 人人澡人人妻人| 免费黄频网站在线观看国产| 日韩在线高清观看一区二区三区| 高清av免费在线| 国产色爽女视频免费观看| 草草在线视频免费看| 久久久久久久大尺度免费视频| 伊人亚洲综合成人网| 老司机影院成人| 免费观看在线日韩| 午夜福利影视在线免费观看| 国产av码专区亚洲av| 嫩草影院入口| 看十八女毛片水多多多| 久久国产亚洲av麻豆专区| 成年人免费黄色播放视频 | 国产精品嫩草影院av在线观看| 夜夜看夜夜爽夜夜摸| 最后的刺客免费高清国语| 亚洲成色77777| 一区二区三区免费毛片| 成人免费观看视频高清| av一本久久久久| 另类亚洲欧美激情| 一边亲一边摸免费视频| 亚洲精品乱码久久久v下载方式| 免费大片黄手机在线观看| 校园人妻丝袜中文字幕| 免费观看在线日韩| 麻豆成人午夜福利视频| 观看av在线不卡| 亚洲人成网站在线播| 一边亲一边摸免费视频| 另类亚洲欧美激情| 成人无遮挡网站| 久久鲁丝午夜福利片| 黄色怎么调成土黄色| 尾随美女入室| 亚洲三级黄色毛片| .国产精品久久| 天美传媒精品一区二区| 亚洲,欧美,日韩| 国产成人免费无遮挡视频| 久久国产精品男人的天堂亚洲 | 国产黄频视频在线观看| 春色校园在线视频观看| 在线观看免费视频网站a站| 亚洲伊人久久精品综合| 婷婷色av中文字幕| 免费高清在线观看视频在线观看| 精品卡一卡二卡四卡免费| 日韩成人伦理影院| 天堂8中文在线网| 欧美一级a爱片免费观看看| 亚洲av.av天堂| 97在线人人人人妻| 精华霜和精华液先用哪个| 国产精品久久久久成人av| 国产伦精品一区二区三区四那| 一本色道久久久久久精品综合| 少妇精品久久久久久久| 美女主播在线视频| 丰满饥渴人妻一区二区三| av播播在线观看一区| 久久久久精品性色| 日韩电影二区| 欧美精品国产亚洲| 免费观看性生交大片5| 免费观看性生交大片5| 人人妻人人看人人澡| 18禁裸乳无遮挡动漫免费视频| 欧美成人午夜免费资源| 一本—道久久a久久精品蜜桃钙片| 亚洲精品乱码久久久v下载方式| 日韩制服骚丝袜av| 欧美 日韩 精品 国产| 中文欧美无线码| 特大巨黑吊av在线直播| 亚洲欧美精品专区久久| 免费黄网站久久成人精品| 亚洲精品,欧美精品| 久久人人爽人人片av| 国产一区二区三区av在线| 国产精品熟女久久久久浪| 久久av网站| 在线观看国产h片| 欧美老熟妇乱子伦牲交| 色视频在线一区二区三区| 欧美高清成人免费视频www| 国产高清有码在线观看视频| 日韩在线高清观看一区二区三区| 国产视频首页在线观看| 午夜激情福利司机影院| 人妻夜夜爽99麻豆av| 性色avwww在线观看| 大香蕉久久网| 男男h啪啪无遮挡| 一个人看视频在线观看www免费| 国产高清国产精品国产三级| 日韩电影二区| 亚洲美女黄色视频免费看| 欧美日韩亚洲高清精品| 日韩伦理黄色片| 新久久久久国产一级毛片| 草草在线视频免费看| 国产成人精品一,二区| xxx大片免费视频| 欧美变态另类bdsm刘玥| 久久精品久久精品一区二区三区| 欧美激情国产日韩精品一区|