• <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玫瑰| 老司机福利观看| 日韩精品青青久久久久久| 亚洲精品自拍成人| 亚洲国产欧洲综合997久久,| 只有这里有精品99| 少妇丰满av| 中文精品一卡2卡3卡4更新| 人妻夜夜爽99麻豆av| 级片在线观看| 国产在视频线精品| 日本免费在线观看一区| 99热精品在线国产| 亚洲av不卡在线观看| 欧美一区二区国产精品久久精品| 少妇高潮的动态图| 国产精品国产三级国产av玫瑰| 日本一本二区三区精品| 久久久精品欧美日韩精品| 日韩中字成人| 床上黄色一级片| 亚洲国产精品成人综合色| 1024手机看黄色片| 精品国内亚洲2022精品成人| 1000部很黄的大片| 日韩三级伦理在线观看| 国内精品宾馆在线| 国产精品麻豆人妻色哟哟久久 | 真实男女啪啪啪动态图| 欧美日韩一区二区视频在线观看视频在线 | 女人十人毛片免费观看3o分钟| 亚洲欧美中文字幕日韩二区| 亚洲av中文字字幕乱码综合| 精品国产露脸久久av麻豆 | 亚洲精品乱久久久久久| 欧美xxxx黑人xx丫x性爽| 亚洲精品国产成人久久av| 1024手机看黄色片| 亚洲av中文字字幕乱码综合| 日本五十路高清| 国产黄a三级三级三级人| 国模一区二区三区四区视频| 国产大屁股一区二区在线视频| 亚洲国产精品合色在线| 精品国产三级普通话版| av国产久精品久网站免费入址| 特大巨黑吊av在线直播| 欧美不卡视频在线免费观看| 草草在线视频免费看| 最近中文字幕高清免费大全6| 午夜激情欧美在线| 国产伦理片在线播放av一区| 久久国产乱子免费精品| 久久6这里有精品| 嘟嘟电影网在线观看| 亚洲性久久影院| 国内少妇人妻偷人精品xxx网站| 亚州av有码| 国模一区二区三区四区视频| 国产成人精品一,二区| 人体艺术视频欧美日本| 国产精品人妻久久久影院| 日韩av在线大香蕉| 亚洲国产精品国产精品| 日本黄色视频三级网站网址| 天堂√8在线中文| 99久国产av精品国产电影| 亚洲国产精品国产精品| 亚洲中文字幕一区二区三区有码在线看| 99久久中文字幕三级久久日本| 亚洲五月天丁香| 男人狂女人下面高潮的视频| 国产黄片美女视频| 在线免费观看的www视频| 国产在视频线精品| 久久这里有精品视频免费| 两个人的视频大全免费| 狂野欧美激情性xxxx在线观看| 精品人妻熟女av久视频| 国产伦在线观看视频一区| 日本五十路高清| 听说在线观看完整版免费高清| 午夜激情欧美在线| 亚洲欧美日韩东京热| av播播在线观看一区| 我要搜黄色片| 九九爱精品视频在线观看| 三级经典国产精品| 久久婷婷人人爽人人干人人爱| 免费搜索国产男女视频| 久久久a久久爽久久v久久| 高清午夜精品一区二区三区| 一边亲一边摸免费视频| 一级黄片播放器| 久久久久久久久中文| 久久精品人妻少妇| 国产精品女同一区二区软件| 免费在线观看成人毛片| 国产成人一区二区在线| 婷婷六月久久综合丁香| 久久国内精品自在自线图片| 国产午夜精品一二区理论片| 深夜a级毛片| 一个人免费在线观看电影| 桃色一区二区三区在线观看| www.av在线官网国产| 久久精品国产自在天天线| 日日摸夜夜添夜夜爱| 99热这里只有精品一区| 九九在线视频观看精品| 国产黄a三级三级三级人| 欧美xxxx性猛交bbbb| 欧美高清成人免费视频www| 国产精品国产三级专区第一集| 国产亚洲午夜精品一区二区久久 | 亚洲最大成人手机在线| 国产精品人妻久久久久久| 国产精品福利在线免费观看| 99国产精品一区二区蜜桃av| 久久久久久久亚洲中文字幕| 汤姆久久久久久久影院中文字幕 | 国产一区二区在线av高清观看| 黄片无遮挡物在线观看| 久久亚洲精品不卡| 国产成人freesex在线| 中文字幕人妻熟人妻熟丝袜美| 青春草国产在线视频| 成人性生交大片免费视频hd| 99久久精品国产国产毛片| 国语对白做爰xxxⅹ性视频网站| 成年免费大片在线观看| 午夜久久久久精精品| av在线老鸭窝| 久久久久久久久中文| av在线亚洲专区| 最近最新中文字幕大全电影3| 欧美+日韩+精品| 人人妻人人看人人澡| 亚洲成av人片在线播放无| 亚洲图色成人| 精品人妻一区二区三区麻豆| 婷婷色av中文字幕| 免费观看性生交大片5| 精品久久国产蜜桃| 亚洲一级一片aⅴ在线观看| 99久久中文字幕三级久久日本| 久久综合国产亚洲精品| 一边摸一边抽搐一进一小说| 久久午夜福利片| 亚洲精品乱久久久久久| 噜噜噜噜噜久久久久久91| 日本午夜av视频| 成年女人看的毛片在线观看| 欧美区成人在线视频| 人妻夜夜爽99麻豆av| 两个人的视频大全免费| 天天躁夜夜躁狠狠久久av| 免费观看a级毛片全部| 亚洲内射少妇av| 亚洲国产精品成人综合色| 三级毛片av免费| 老司机福利观看| 欧美成人免费av一区二区三区| 亚洲,欧美,日韩| 国产av一区在线观看免费| 久久久久久久久久成人| 天天躁夜夜躁狠狠久久av| 国产精品久久久久久精品电影| 日韩av在线大香蕉| 国产黄色小视频在线观看| 日日摸夜夜添夜夜爱| 超碰av人人做人人爽久久| 一本久久精品| 成人鲁丝片一二三区免费| 国产亚洲一区二区精品| 久久这里有精品视频免费| 国产老妇女一区| 亚洲欧美精品专区久久| 亚洲精品久久久久久婷婷小说 | 国产精品女同一区二区软件| 国产成年人精品一区二区| 亚洲自偷自拍三级| 亚洲av电影在线观看一区二区三区 | 国产淫语在线视频| 国产av码专区亚洲av| 国产老妇女一区| 亚洲综合色惰| 丝袜美腿在线中文| 毛片女人毛片| 1024手机看黄色片| 成人一区二区视频在线观看| 久久久久国产网址| 免费看光身美女| 在线观看av片永久免费下载| 中文天堂在线官网| 色吧在线观看| 国产亚洲91精品色在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 啦啦啦韩国在线观看视频| 大香蕉久久网| 波野结衣二区三区在线| 亚洲精品aⅴ在线观看| 哪个播放器可以免费观看大片| 免费电影在线观看免费观看| 丰满乱子伦码专区| 亚洲av免费高清在线观看| 亚洲国产成人一精品久久久| 久久久欧美国产精品| 看十八女毛片水多多多| 国产探花极品一区二区| 精品一区二区三区人妻视频| 中文字幕av成人在线电影| 亚洲av.av天堂| 狂野欧美激情性xxxx在线观看| 午夜老司机福利剧场| 亚洲精品影视一区二区三区av| 一级黄片播放器| 国产又色又爽无遮挡免| av在线亚洲专区| 久久精品国产亚洲av涩爱| 免费看光身美女| 内射极品少妇av片p| 久久精品国产自在天天线| 18禁裸乳无遮挡免费网站照片| 一级爰片在线观看| 免费看光身美女| 国产一区二区三区av在线| 爱豆传媒免费全集在线观看| 两个人视频免费观看高清| 欧美3d第一页| 亚洲无线观看免费| 久久久亚洲精品成人影院| 色哟哟·www| 国产精品久久久久久精品电影小说 | 亚洲久久久久久中文字幕| 亚洲在久久综合| 91av网一区二区| 亚洲av电影在线观看一区二区三区 | av天堂中文字幕网| 亚洲av男天堂| 一本一本综合久久| www.av在线官网国产| 久久久久久国产a免费观看| 国产久久久一区二区三区| 波多野结衣高清无吗| av国产久精品久网站免费入址| 中文资源天堂在线| 九色成人免费人妻av| 亚洲不卡免费看| 自拍偷自拍亚洲精品老妇| 久久国内精品自在自线图片| 一级黄色大片毛片| 女人十人毛片免费观看3o分钟| 波多野结衣巨乳人妻| 高清日韩中文字幕在线| 又爽又黄无遮挡网站| 欧美精品一区二区大全| 亚洲综合精品二区| 伦理电影大哥的女人| 男女边吃奶边做爰视频| 午夜精品一区二区三区免费看| eeuss影院久久| 水蜜桃什么品种好| 亚洲电影在线观看av| 日韩人妻高清精品专区| 一级毛片久久久久久久久女| 国产成人午夜福利电影在线观看| 精品久久久久久久久av| 午夜精品一区二区三区免费看| 亚洲欧洲国产日韩| 如何舔出高潮| 亚洲精品乱码久久久久久按摩| 亚洲精品国产av成人精品| 少妇熟女aⅴ在线视频| 亚洲国产欧美人成| 成人一区二区视频在线观看| 国产精品一区二区三区四区久久| 中文在线观看免费www的网站| av免费观看日本| 最近的中文字幕免费完整| 国产精品久久久久久久电影| 国产免费视频播放在线视频 | 日韩欧美精品v在线| 能在线免费观看的黄片| 人妻夜夜爽99麻豆av| 国模一区二区三区四区视频| 亚洲成色77777| 免费观看人在逋| 午夜视频国产福利| 精品久久国产蜜桃| 亚洲丝袜综合中文字幕| 欧美精品国产亚洲| 国产精品麻豆人妻色哟哟久久 | 两个人视频免费观看高清| 国产精品人妻久久久久久| 国产国拍精品亚洲av在线观看| 久久精品国产亚洲av天美| 黄片无遮挡物在线观看| 久久精品久久久久久久性| 成人毛片60女人毛片免费| 亚洲国产最新在线播放| 亚洲av二区三区四区| 视频中文字幕在线观看| 亚洲,欧美,日韩| 日韩强制内射视频| 日韩在线高清观看一区二区三区| 久久久精品94久久精品| 看片在线看免费视频| 免费不卡的大黄色大毛片视频在线观看 | 中文在线观看免费www的网站| 26uuu在线亚洲综合色| 久久久国产成人精品二区| 91午夜精品亚洲一区二区三区| 国产精品熟女久久久久浪| av天堂中文字幕网| 一级毛片我不卡| 国产毛片a区久久久久| 超碰av人人做人人爽久久| 人人妻人人看人人澡| 亚洲内射少妇av| 日本一本二区三区精品| 精品午夜福利在线看| 天堂影院成人在线观看| 久久精品夜夜夜夜夜久久蜜豆| 久久草成人影院| 黄片无遮挡物在线观看| 一级av片app| 2022亚洲国产成人精品| 韩国高清视频一区二区三区| 亚洲精品,欧美精品| 一本久久精品| 特大巨黑吊av在线直播| 亚洲欧美中文字幕日韩二区| 91精品一卡2卡3卡4卡| 亚洲成人中文字幕在线播放| 国产淫片久久久久久久久| 九草在线视频观看| 国产精品爽爽va在线观看网站| 又爽又黄无遮挡网站| 日本黄色视频三级网站网址| 精品久久久久久久久av| 非洲黑人性xxxx精品又粗又长| 伊人久久精品亚洲午夜| 免费av毛片视频| 熟女人妻精品中文字幕| 伦精品一区二区三区| 九九久久精品国产亚洲av麻豆| 老司机影院成人| 特级一级黄色大片| 最近最新中文字幕免费大全7| 熟女人妻精品中文字幕| 精品酒店卫生间| 国产成人福利小说| 亚洲精品乱码久久久v下载方式| 男的添女的下面高潮视频| 卡戴珊不雅视频在线播放| 国产又色又爽无遮挡免| 在线免费观看的www视频| 蜜桃亚洲精品一区二区三区| 两个人的视频大全免费| 天天躁日日操中文字幕| 免费看光身美女| 日本爱情动作片www.在线观看| 国产伦一二天堂av在线观看| 岛国在线免费视频观看| 三级经典国产精品| 国语对白做爰xxxⅹ性视频网站| 美女脱内裤让男人舔精品视频| 亚洲国产最新在线播放| 精品熟女少妇av免费看| 久久久国产成人免费| 午夜福利网站1000一区二区三区| 最近的中文字幕免费完整| 国产精品一区二区在线观看99 | 欧美激情国产日韩精品一区| 日韩人妻高清精品专区| 一区二区三区四区激情视频| 免费观看精品视频网站| 国产一区二区三区av在线| 日韩欧美精品免费久久| 久久这里有精品视频免费| 国产伦精品一区二区三区视频9| 少妇熟女欧美另类| 精品国产露脸久久av麻豆 | 床上黄色一级片| 亚洲国产精品久久男人天堂| 亚洲精品久久久久久婷婷小说 | 国产激情偷乱视频一区二区| 日韩制服骚丝袜av| 久久久国产成人精品二区| ponron亚洲| 18禁裸乳无遮挡免费网站照片| 搞女人的毛片| 色综合站精品国产| videos熟女内射| 日韩欧美精品v在线| 精品人妻熟女av久视频| 国内精品美女久久久久久| 丝袜美腿在线中文| 九草在线视频观看| 最后的刺客免费高清国语| av天堂中文字幕网| 精品欧美国产一区二区三| 成人毛片60女人毛片免费| 水蜜桃什么品种好| 黄片无遮挡物在线观看| 欧美高清性xxxxhd video| 亚洲怡红院男人天堂| 中文在线观看免费www的网站| 黑人高潮一二区| 夜夜爽夜夜爽视频| 亚洲最大成人av| 伦精品一区二区三区| av视频在线观看入口| 三级经典国产精品| 大香蕉97超碰在线| 久久热精品热| 亚洲精品456在线播放app| 男女下面进入的视频免费午夜| 国产久久久一区二区三区| 99久久中文字幕三级久久日本| 国产免费又黄又爽又色| 91久久精品国产一区二区三区| 91狼人影院| 日日干狠狠操夜夜爽| 国产精品一区二区在线观看99 | 成年免费大片在线观看| 免费播放大片免费观看视频在线观看 | 伊人久久精品亚洲午夜| 成人无遮挡网站| 又黄又爽又刺激的免费视频.| 搡女人真爽免费视频火全软件| 久久久久久久久久久丰满| 日韩大片免费观看网站 | 久久精品夜夜夜夜夜久久蜜豆| 午夜久久久久精精品| 国产亚洲91精品色在线| 国产成人午夜福利电影在线观看| 狠狠狠狠99中文字幕| 一级黄色大片毛片| 国产精品一区www在线观看| 国产精品国产三级专区第一集| 久热久热在线精品观看| av免费在线看不卡| 99久久无色码亚洲精品果冻| 国产精品1区2区在线观看.| 性色avwww在线观看| 边亲边吃奶的免费视频| 51国产日韩欧美| 日韩欧美精品v在线| 亚洲不卡免费看| 五月玫瑰六月丁香| 久久久久免费精品人妻一区二区| 男女视频在线观看网站免费| 精品熟女少妇av免费看| 人妻系列 视频| 精华霜和精华液先用哪个| 精品久久久久久久久久久久久| 97超视频在线观看视频| 日韩在线高清观看一区二区三区| 亚洲av熟女| 国产熟女欧美一区二区| 日本欧美国产在线视频| 亚洲真实伦在线观看| 欧美另类亚洲清纯唯美| 我的女老师完整版在线观看| 男女那种视频在线观看| 亚洲国产欧美人成| 国产毛片a区久久久久| 禁无遮挡网站| 2022亚洲国产成人精品| 国产v大片淫在线免费观看| 久久久成人免费电影| 啦啦啦韩国在线观看视频| 欧美一区二区亚洲| 国产av码专区亚洲av| 最近的中文字幕免费完整| 3wmmmm亚洲av在线观看| 级片在线观看| 欧美一区二区国产精品久久精品| 国产成人freesex在线| 国产在视频线精品| 国产av在哪里看| 女人被狂操c到高潮| 看非洲黑人一级黄片| 久久人妻av系列| 久久精品综合一区二区三区| 国产精品.久久久| 亚洲欧美日韩高清专用| 欧美日韩在线观看h| 中文字幕人妻熟人妻熟丝袜美| 久久精品综合一区二区三区| 成年版毛片免费区| 国产真实乱freesex| 97人妻精品一区二区三区麻豆| videos熟女内射| 淫秽高清视频在线观看| 黄色日韩在线| 国产美女午夜福利| 成人三级黄色视频| 亚洲人成网站高清观看| 亚洲av一区综合| 人体艺术视频欧美日本| 人妻系列 视频| 热99在线观看视频| 久久精品久久久久久久性| 亚洲av男天堂| 人人妻人人看人人澡| 国产69精品久久久久777片| 国产 一区精品| 色尼玛亚洲综合影院| 国国产精品蜜臀av免费| 看十八女毛片水多多多| 简卡轻食公司| 亚洲av中文av极速乱| 99热6这里只有精品| 男女啪啪激烈高潮av片| 国产乱来视频区| 国产 一区精品| 精品久久久久久成人av| 久久久成人免费电影| 别揉我奶头 嗯啊视频| 看非洲黑人一级黄片| 中文字幕制服av| 亚洲欧洲国产日韩| 精品人妻一区二区三区麻豆| 最近最新中文字幕免费大全7| 国产一区有黄有色的免费视频 | 亚洲乱码一区二区免费版| 久久久久性生活片| 蜜臀久久99精品久久宅男| 久久精品久久久久久久性| 中文欧美无线码| 免费看日本二区| 热99在线观看视频| 午夜久久久久精精品| 美女黄网站色视频| 国产日韩欧美在线精品| 亚洲欧洲日产国产| 久久久欧美国产精品| 激情 狠狠 欧美| 日韩国内少妇激情av| 性插视频无遮挡在线免费观看| 中文亚洲av片在线观看爽| 国产亚洲av嫩草精品影院| 两个人的视频大全免费| 精品无人区乱码1区二区| 人人妻人人澡欧美一区二区| 特大巨黑吊av在线直播| 亚洲av成人精品一二三区| 国产在线男女| 亚洲国产最新在线播放| 村上凉子中文字幕在线| 又黄又爽又刺激的免费视频.| 在线观看美女被高潮喷水网站| 免费观看在线日韩| 亚洲精品456在线播放app| 久久人人爽人人爽人人片va| 女人久久www免费人成看片 | 97人妻精品一区二区三区麻豆| 大又大粗又爽又黄少妇毛片口| 女人十人毛片免费观看3o分钟| 22中文网久久字幕| 欧美成人午夜免费资源| 激情 狠狠 欧美| 亚洲电影在线观看av| 国产精品蜜桃在线观看| 国产av不卡久久| 黄色一级大片看看| 亚洲欧美一区二区三区国产| 午夜a级毛片| 国产伦精品一区二区三区四那| 亚洲婷婷狠狠爱综合网| 久久99精品国语久久久| 国产日韩欧美在线精品| 丝袜美腿在线中文| 久久欧美精品欧美久久欧美| av在线老鸭窝| 日日摸夜夜添夜夜爱| 波多野结衣高清无吗| 国产色爽女视频免费观看| 狠狠狠狠99中文字幕| 91午夜精品亚洲一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 日本免费a在线| 身体一侧抽搐| 亚洲av不卡在线观看| a级毛片免费高清观看在线播放| 精品国产露脸久久av麻豆 | 国产又黄又爽又无遮挡在线| 69人妻影院| 22中文网久久字幕| 老师上课跳d突然被开到最大视频| 视频中文字幕在线观看| 久久久久久九九精品二区国产| 国产乱人偷精品视频| 久久综合国产亚洲精品| 91久久精品电影网| 国产在视频线在精品| 亚洲av一区综合| 成人毛片60女人毛片免费| 91精品国产九色| 久久国产乱子免费精品| 久久精品久久久久久久性| 少妇被粗大猛烈的视频| av在线老鸭窝| 日韩大片免费观看网站 | www.av在线官网国产| 最后的刺客免费高清国语| 一区二区三区高清视频在线| 国产精品国产三级国产专区5o | 午夜福利成人在线免费观看| 99久久中文字幕三级久久日本| 国产在线男女| 欧美激情在线99|