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

    含剪切節(jié)理面巖質(zhì)邊坡滑裂面位置及穩(wěn)定性研究

    2024-06-09 00:00:00陳東宇劉文連眭素剛許漢華
    關(guān)鍵詞:安全系數(shù)

    摘 要:同向雙平面滑動是存在單一地質(zhì)斷層面(剪切節(jié)理面)巖質(zhì)邊坡的常見破壞模式之一,但對該種類型的滑裂面計算方法并不充足。為了能夠更加高效準確地尋找邊坡的滑裂面位置,判斷邊坡的穩(wěn)定性,基于極限平衡理論和非線性數(shù)學(xué)規(guī)劃模型,假設(shè)滑體的滑動方式為同向雙平面滑動,再假設(shè)目標函數(shù)為該巖質(zhì)邊坡的安全系數(shù),運用MTALAB全局最優(yōu)搜索法,計算含剪切節(jié)理面工程邊坡在天然工況作用下的滑裂面位置及穩(wěn)定性,并與極限分析法、強度折減法和畢肖普法進行對比分析。研究結(jié)果表明:基于極限平衡理論和非線性數(shù)學(xué)規(guī)劃模型得出的滑裂面位置與安全系數(shù)基本一致,驗證此類方法的可行性,為存在單一地質(zhì)斷層面巖質(zhì)邊坡的滑裂面計算和穩(wěn)定性分析提供了新依據(jù)。

    關(guān)鍵詞:巖質(zhì)邊坡;滑裂面計算;安全系數(shù);極限平衡法;非線性理論;最優(yōu)化原理

    中圖分類號:TU452

    文獻標志碼:A

    巖質(zhì)邊坡的地層是分布在我國西南地區(qū)的一種工程性質(zhì)較差的軟巖,從上世紀60年代國家開始對西南地區(qū)進行開發(fā)建設(shè),巖質(zhì)邊坡就慢慢地出現(xiàn)在工程作業(yè)中。在工程作業(yè)中,滑坡等地質(zhì)災(zāi)害就出現(xiàn)在了我們的視線當(dāng)中,其不僅制約了工程建設(shè)的發(fā)展,而且還對人類的生產(chǎn)生活產(chǎn)生了許多影響[1-2]。針對這一特殊的軟質(zhì)巖土,預(yù)測邊坡穩(wěn)定性及滑坡失穩(wěn)面的位置,對土建、水利、交通等方面有非常重要的意義。

    在結(jié)構(gòu)面控制型邊坡中,主要的失穩(wěn)模式有平面滑動、楔型滑動、傾倒破壞等[3-4]。目前常見的邊坡穩(wěn)定性預(yù)測分析方法有:極限平衡法、強度折減法、系統(tǒng)分析法、數(shù)值模擬等[5-8]。李芬[9]、朱擎[10]、程小龍[11]等通過強度折減法和雙強度折減法解決了雙平面破壞模式下的臨界失穩(wěn)問題,討論了巖質(zhì)邊坡的破裂角和內(nèi)摩擦角之間的關(guān)系,進而提出安全系數(shù)的簡化估算公式。陳建宏等[12]使用了極限平衡分析的上下限法,并把模糊化處理參數(shù)和置信水平等概念加入到了平面滑動的巖質(zhì)邊坡工程中,此方法能夠有效地解決平面滑動巖質(zhì)邊坡的參數(shù)不確定問題。CHENG[13]、ARDESTANI[14]等采用極限平衡法,通過二維平面應(yīng)力分析求解邊坡的臨界滑動面及穩(wěn)定性系數(shù)。SCHLOTFELDT等[15]以反傾巖質(zhì)邊坡為研究對象,把數(shù)值分析與極限平衡法相結(jié)合,提出了這類巖質(zhì)邊坡的研究方法。張崇波[16]對香港秀茂坪巖質(zhì)邊坡進行了可靠性分析,利用拉丁超立方抽樣可靠性分析方法,分析了同時考慮參數(shù)不確定性和最危險滑面不確定性的平面滑動巖質(zhì)邊坡可靠性,優(yōu)化了可靠性分析方法的計算精度和效率,將研究成果應(yīng)用到工程實踐,并將該方法應(yīng)用到山東臨沂換流站站址邊坡的評價當(dāng)中。陳志強等[17]提出了如果巖質(zhì)邊坡的滑裂面位置和產(chǎn)狀無法確定,就可以用巖質(zhì)邊坡平面滑動滑面極限傾角來計算邊坡的最小安全系數(shù),并考慮拉裂縫深度對滑坡穩(wěn)定性的影響。高丙麗等[18]基于坐標投影法,提出了三維單滑面和雙滑面型塊體的穩(wěn)定系數(shù)計算方法,并基于Matlab開發(fā)出適于巖質(zhì)邊坡工程中平面多面體塊體和曲面塊體穩(wěn)定性分析的CPG程序,實現(xiàn)了結(jié)構(gòu)面、臨空面及不穩(wěn)定塊體的空間表示及可視化。但對存在剪切節(jié)理面巖質(zhì)邊坡的滑裂面位置的計算除了工程軟件外,并沒有其他更加簡便的方法了。

    針對以上研究的不足之處,本文以兩種不同地層的巖質(zhì)邊坡為研究對象,將極限平衡理論和非線性數(shù)學(xué)規(guī)劃理論結(jié)合起來,通過求解最優(yōu)化方法確定安全系數(shù)從而求解巖質(zhì)邊坡滑裂面,建立一種直接計算同向雙平面滑動邊坡滑裂面位置的計算方法。

    1 計算原理

    對于穩(wěn)定性受確定性結(jié)構(gòu)面控制的巖質(zhì)邊坡,應(yīng)根據(jù)結(jié)構(gòu)面的產(chǎn)狀與強度參數(shù),采用極限平衡方法計算邊坡的穩(wěn)定性。根據(jù)結(jié)構(gòu)面的空間展布,邊坡失穩(wěn)類型分別為平面滑動和空間滑動。嚴格意義上邊坡滑體都是空間塊體,但對于單一結(jié)構(gòu)面控制的滑體,或由兩個及兩個以上平面構(gòu)成的畫面,只要這些平面走向大致相同、與邊坡破面走向平行或接近平行,且滑體兩側(cè)不受約束或約束不大時,即可按平面滑動進行分析,即本文算法的基本原理。如圖1所示,本文研究對象主要是存在兩種不同地層的復(fù)合巖質(zhì)邊坡,該類邊坡主要是由于不同紀元產(chǎn)生的不同類型巖石共同組成的一種特殊的復(fù)合地層巖質(zhì)邊坡,地層Ⅱ的強度較大,其次是地層Ⅰ,該巖質(zhì)邊坡強度較弱的部分即最容易發(fā)生滑裂的部位是兩種地層的交界面(剪切節(jié)理面),通常的破壞形式為同向雙平面滑動。傳統(tǒng)計算方法常見的有有限元法、強度折減法、極限平衡法,可通過工程軟件大量的計算來實現(xiàn),但計算結(jié)果通常為圓弧型滑裂面或為不規(guī)則滑動面,不符合邊坡實際滑裂面形狀。

    本文算法基于塑性力學(xué)的基本理論,將邊坡安全系數(shù)作為目標函數(shù),將底滑面和后緣滑裂面的位置坐標作為決策變量,同時考慮滑體非線性約束,如極限平衡方程,底滑面、后緣滑裂面的屈服條件,即Mohr-Coulomb屈服條件,最后使用優(yōu)化算法可同時求解得到邊坡的安全系數(shù)及底滑面和后緣滑裂面控制點的位置坐標。

    2 邊坡滑塊計算模型

    將邊坡滑塊離散出來成為一塊狀單元如圖2所示,該塊體單元為一四面體,每條邊為一結(jié)構(gòu)面,除了該單元體的自重,底滑面和后緣滑裂面還受到法向力和剪力;由于已有極限平衡法里的剛性假設(shè),該滑體的內(nèi)部不會發(fā)生變形,故可以直接對該滑體進行受力分析。該滑塊單元的受力如圖2所示,滑塊在形心處受到自身重力G,θ為兩種地層交界線(底滑面)與水平方向的夾角,β為后緣滑裂面與水平方向的夾角,由于不施加外部荷載,故臨空面AD和頂面CD是不受力的,滑塊底滑面AB受到支持力NAB和剪切力SAB,后緣滑裂面也受到支持力NBC和剪切力SBC,此時該滑塊單元處于極限平衡狀態(tài)。

    3 巖質(zhì)邊坡的力學(xué)機制分析

    3.1 滑塊的力學(xué)模型及幾何關(guān)系

    根據(jù)含剪切節(jié)理面的巖質(zhì)邊坡變形破壞機制,可將發(fā)生破裂的面分為底滑面AB和后緣滑裂面BC兩個部分,如圖1所示,滑裂面、臨空面和頂面共同組成滑塊。

    假設(shè)滑塊為均質(zhì)剛體,則滑塊的非線性約束重力方程為

    式中:G為滑體的自重;ρ為第四系地層的密度;g為重力加速度;(xA,yA)、(xB,yB)、(xC,yC)、(xD,yD)分別為A、B、C、D 4點的坐標。

    3.2 底滑面和后緣滑裂面的幾何關(guān)系

    如圖1所示,以O(shè)點為原點來確定A、B、C、D各點的坐標。

    (1)底滑面AB的坐標及長度描述

    式中:θ是底滑面AB與水平方向的夾角,θ以逆時針方向為正;xA是A點的x坐標,yA是y點的y坐標,xB是B點的x坐標,yB是B點的y坐標,lAB是底滑面AB的長度。

    (2)后緣滑裂面BC的坐標及長度描述

    式中:β是后緣滑裂面與水平方向的夾角,β以逆時針為正;xC是C點的x坐標,yC是C點的y坐標,lBC是后緣滑裂面BC的長度。

    (3)后緣滑裂面BC兩點的附加約束

    式中:xD是D點的x坐標;L是巖質(zhì)邊坡的寬度;d是邊坡基底的高度;H是巖質(zhì)邊坡的高度。其中,xB在xD和L之間取間隔為1。

    3.3 含剪切節(jié)理面巖質(zhì)邊坡滑塊的非線性數(shù)學(xué)規(guī)劃模型

    根據(jù)滑塊的平衡方程、屈服條件和幾何約束條件,可以建立約束非線性數(shù)學(xué)規(guī)劃模型,由于約束條件較多,為了確保求解方程的可行性,本文采取最優(yōu)化方法中的全局最優(yōu)解。

    (1)建立目標函數(shù)

    Maximize K

    式中:K為巖質(zhì)邊坡的安全系數(shù);Maximize表示“使最大”。

    (2)建立滑體的非線性約束平衡方程

    式中:NAB是底滑面受到的法向力,NAB以受壓為正;SAB是底滑面AB的剪力,SAB以對滑體產(chǎn)生逆時針的轉(zhuǎn)動效果為正;NBC是后緣滑裂面BC的法向力,NBC以受壓為正;SBC是后緣滑裂面BC的剪力,SBC以對滑體產(chǎn)生逆時針的轉(zhuǎn)動效果為正。

    (3)建立底滑面AB的非線性約束Mohr-Coulomb屈服條件

    式中:K為巖質(zhì)邊坡的安全系數(shù);NAB是底滑面AB受到的法向力,NAB以受壓為正;SAB是底滑面AB的剪力,SAB以對滑體產(chǎn)生逆時針的轉(zhuǎn)動效果為正;φAB是底滑面AB的內(nèi)摩擦角;cAB是底滑面AB的凝聚力;lAB是底滑面AB的長度。

    (4)建立后緣滑裂面BC的非線性約束Mohr-Coulomb屈服條件

    式中:K為巖質(zhì)邊坡的安全系數(shù);NBC是后緣滑裂面BC受到的法向力,NBC以受壓為正;SBC是后緣滑裂面BC的剪力,SBC以對滑體產(chǎn)生逆時針的轉(zhuǎn)動效果為正;φBC是后緣滑裂面BC的內(nèi)摩擦角;cBC是后緣滑裂面BC的凝聚力;lBC是后緣滑裂面BC的長度。

    (5)建立非線性數(shù)學(xué)規(guī)化模型

    將上述目標函數(shù)安全系數(shù)K、滑塊受力的極限平衡方程式、底滑面和后緣滑裂面的屈服條件式以及各點坐標的附加幾何約束條件式聯(lián)立,即可得到復(fù)合地層巖質(zhì)邊坡滑體的非線性數(shù)學(xué)規(guī)劃模型如下:

    上式中,部分參數(shù)可以在邊坡中測量得知。

    4 非線性數(shù)學(xué)規(guī)劃模型的求解

    巖質(zhì)邊坡后緣滑裂面位置穩(wěn)定性計算方法的模型是一種比較典型的非線性數(shù)學(xué)規(guī)劃模型,本文的模型計算流程如下。

    (1)定義目標函數(shù):定義滑面的材料參數(shù),抗剪斷凝聚力和內(nèi)摩擦角度;定義幾何、荷載參數(shù),包括用幾何法確定各點坐標和角度等確定部分并給予賦值,將未知部分作為未知數(shù)后續(xù)來進行求解。

    (2)列出復(fù)合地層巖質(zhì)邊坡滑體的非線性數(shù)學(xué)規(guī)劃的函數(shù)模型:先建立底滑面和后緣滑裂面的方程和滑體的重力計算公式;再根據(jù)上一步驟點的坐標確定底滑面lAB和后緣滑裂面lBC的長度;最后列出滑體的平衡方程和底滑面、后緣滑裂面的屈服條件以及部分額外幾何平面約束。

    (3)利用MATLAB軟件中的功能函數(shù)進行求解,得到巖質(zhì)邊坡的安全系數(shù)K和后緣滑裂面BC的具體位置,即得到復(fù)合地層巖質(zhì)邊坡的后緣滑裂面位置的最優(yōu)解。

    5 基本算例分析

    5.1 邊坡算例參數(shù)

    本文選取的巖質(zhì)邊坡具體參數(shù)為坡底寬L=100 m,坡高H=35 m,邊坡基底高度d=15 m,邊坡頂面寬度為67.19 m;臨空面AD的傾斜角度為68°,底滑面AB的傾角為30°,邊坡巖體的密度ρ=2 500 kg/m3,后緣滑裂面BC的傾角β為決策變量;對于滑體ABCD,其中,A點和D點坐標已知xA=19.920 8 m、yA=18.457 1 m、xD=32.806 2 m、yD=50 m,B點和C點坐標為決策變量。

    在該算例中,由于滑塊長期只受重力荷載的影響,故底滑面AB和后緣滑裂面BC的凝聚力和內(nèi)摩擦角均為固定值,具體參數(shù)見表1。

    5.2 本文算法的邊坡算例計算結(jié)果

    巖質(zhì)邊坡滑體的非線性數(shù)學(xué)規(guī)劃模型是一個非線性數(shù)學(xué)規(guī)劃模型,其目標函數(shù)是安全系數(shù),使用全局最優(yōu)解求解巖質(zhì)邊坡滑體的非線性數(shù)學(xué)規(guī)劃模型,可求解得到以下決策變量:K、xC、β。其中,xB取坐標間距為1,根據(jù)xB的確定坐標求解得到的決策變量如表2所示。

    由表2可以看出,MTALAB最優(yōu)化算法,在以xB取間隔為1且能得到最優(yōu)解的情況下得到的安全系數(shù)最小為0.791 7,即文本文算法最終結(jié)果。

    5.3 本文算法結(jié)果驗證

    為了驗證最優(yōu)化算法求解最優(yōu)滑裂面和安全系數(shù)的可靠性與和合理性,采用OPTUM G2 2020軟件與Rocscience.Slide.v6007軟件進行對比驗算。本次驗算通過算例選取模型進行建模并施加標準邊界條件,選用OPTUM G2 2020軟件中極限分析的上下限法及強度折減的上下限法進行滑裂面搜索,但搜索的滑裂面并非平面雙折線型滑裂面,搜索得到的巖質(zhì)邊坡總位移如圖3—6,圖中紅色實線為極限平衡最優(yōu)搜索法結(jié)果,綠色多段線為軟件搜索滑裂面,為了方便對比參考,將軟件搜索的滑裂面取底滑面的右端點與滑裂面和頂面的交點的連線近似模擬同向雙平面型滑裂面,即藍色直線;選用Rocscience.Slide.v6007軟件中的畢肖普法可以完成非圓弧滑裂面(平面雙折線型滑裂面)的搜索,搜索得到巖質(zhì)邊坡的安全系數(shù)K,底滑面和后緣滑裂面的A點和B點坐標及后緣滑裂面傾角β如表3所示。

    圖3—6分別為G2軟件通過極限分析和強度折減法的上下限法搜索的滑裂面位置,可以看出,G2軟件搜索的滑裂面底滑面為直線形,但后緣滑裂面為近似圓弧形。在近似滑裂面與本文方法的結(jié)果滑裂面對比中,底滑面的長度誤差在0.5 m之內(nèi),后緣滑裂面誤差在1~2 m之間。從圖中的安全系數(shù)對比,基于極限平衡法的最優(yōu)化搜索的安全系數(shù)低于G2軟件的各個方法所求的安全系數(shù),最大差值為0.130 8,差距最小的是與極限平衡法的下限法的對比,誤差為8.36%,差距最大的是與強度折減法的上限法的對比,誤差為16.52%。通過對比近似滑裂面的傾角,極限平衡和強度折減法的上下限法與剛體平衡法的后緣滑裂面傾角相差3~4°,誤差為5.09%~6.78%。

    表3為畢肖普法與全局最優(yōu)搜索法結(jié)果,從表3可以看出,本文的全局最優(yōu)搜索法與Slide軟件的畢肖普法搜索結(jié)果非常接近,誤差最大的依然為所計算的安全系數(shù),誤差為11.04%。對于滑裂面位置,通過兩點坐標對比分析,誤差最大的為B點x坐標,誤差為5.88%,而C點坐標誤差僅為0.60%。后緣滑裂面傾角畢肖普法得到的結(jié)果為56.71°,本文算法得到的結(jié)果為58.91°,誤差為3.88%。

    圖7—9為畢肖普法與全局最優(yōu)搜索法的結(jié)果,從圖7—9可知,全局最優(yōu)搜索法與Slide軟件的畢肖普法搜索結(jié)果非常接近。兩種方法的安全系數(shù)都隨后緣滑裂面凝聚力c的增大而增大,在后緣滑裂面凝聚力c=120 kPa時差值最大,差值為0.176 3,誤差為16.45%。對于滑裂面位置,通過兩點坐標對比分析,B點x坐標在后緣滑裂面凝聚力c=85 kPa時差值最大,差值為2.84 m,但對一個100 m寬的邊坡來說,誤差僅為2.84%,而C點x坐標在后緣滑裂面凝聚力c=115 kPa時差值最大,差值為1.81 m,誤差僅為1.81%。

    5.4 結(jié)果分析與討論

    通過算例驗證,我們不難發(fā)現(xiàn)本文算法的安全系數(shù)普遍偏低,相差最多為0.1左右,但在邊坡防護治理工作中,可以起到預(yù)防邊坡失穩(wěn)的作用。由于計算方法和計算原理的不同,誤差必然是存在的,在圖3—7中我們可以發(fā)現(xiàn)極限平衡和強度折減的上下限法的滑裂面位置與本文算法搜索的滑裂面位置差距較大,這是因為OPTUM G2 軟件無法實現(xiàn)同向雙平面型滑裂面搜索,對于能夠?qū)崿F(xiàn)同向雙平面型滑裂面搜索的Slide軟件,畢肖普法與全局最優(yōu)搜索法(本文算法)結(jié)果非常接近。

    通過觀察后緣滑裂面傾角,我們可以發(fā)現(xiàn)復(fù)合地層巖質(zhì)邊坡的后緣滑裂面傾角接近60°,根據(jù)邊坡滑動破壞的普遍形式,平面滑動的邊坡巖體是沿單一地質(zhì)斷層面(剪切節(jié)理面)的剪切位移,此時邊坡傾角β、地層交界面傾角θ及其內(nèi)摩擦角φ之間的關(guān)系為βgt;θgt;φ。本文算法的實施例以及Slide軟件的計算結(jié)果符合規(guī)定,進一步驗證本文算法的可行性與正確性。

    6 結(jié)論

    (1)本文以復(fù)合地層巖質(zhì)邊坡為研究對象,基于非線性理論的非線性數(shù)學(xué)規(guī)劃模型計算復(fù)合地層巖質(zhì)邊坡的穩(wěn)定性及底滑面和后緣滑裂面位置,通過最優(yōu)化方法中的全局最優(yōu)搜索的計算結(jié)果,并與Slide軟件的結(jié)果進行對比,驗證了本文算法的可行性和準確性。

    (2)傳統(tǒng)的極限平衡法計算邊坡的穩(wěn)定性系數(shù)時,需要假定破壞面,然后再根據(jù)假定的每一個破壞面的安全系數(shù),來確定邊坡最易發(fā)生破壞的滑裂面,強度折減法對c、tan φ簡單折減后得到的應(yīng)力場已經(jīng)不是原巖質(zhì)邊坡的真實應(yīng)力場,故而得出的結(jié)果也存在偏差,而本文的方法可以省略滑裂面位置的假設(shè),直接通過全局優(yōu)化算出邊坡的滑裂面位置,相比于極限平衡法,本文的編制程序更加簡單、工程應(yīng)用簡便,可將其應(yīng)用于含剪切節(jié)理面地層巖質(zhì)邊坡的設(shè)計中。

    (3)計算邊坡穩(wěn)定性及邊坡滑裂面的工程軟件很多,主要以圓弧形滑裂面為主,但計算平面型滑裂面算法的軟件還需開發(fā),本文算法程序簡單且計算精度較高,為同向雙平面滑動的邊坡提供了一種簡單高效的計算方法。在本文算法中,邊坡的滑裂面位置與Slide軟件的計算結(jié)果非常接近,安全系數(shù)比軟件的計算結(jié)果要普遍偏低,在工程應(yīng)用中可以起到防范于未然的作用。

    (4)本文提出了一種用于特殊巖質(zhì)邊坡折線形滑裂面分析的計算方法,旨在解決現(xiàn)有軟件如OPTUM G2和Madis civil在處理此類問題時的不足。所開發(fā)的計算程序具有簡單直觀、原理明確且計算精度高的特點。本方法理論基礎(chǔ)扎實,且在工程實踐中操作簡便,適用于處理含有單個地質(zhì)斷裂面(如斷層面或剪切節(jié)理面)的巖質(zhì)邊坡的滑裂面位置分析。該算法不僅能夠準確預(yù)測滑裂面的位置,還能清晰地展示滑塊的自重、底部滑動面和后緣滑裂面的受力狀況,從而直觀地呈現(xiàn)邊坡的穩(wěn)定性狀況。此外,算法提供的詳細信息對于邊坡支護結(jié)構(gòu)的設(shè)計同樣具有重要價值。

    參考文獻:

    [1]黃潤秋. 20世紀以來中國的大型滑坡及其發(fā)生機制[J]. 巖石力學(xué)與工程學(xué)報, 2007(3): 433-454.

    [2] 夏元友, 李梅. 邊坡穩(wěn)定性評價方法研究及發(fā)展趨勢[J]. 巖石力學(xué)與工程學(xué)報, 2002(7): 1087-1091.

    [3] 楊肖鋒, 魯祖德, 陳從新, 等. 板裂結(jié)構(gòu)順層巖質(zhì)邊坡滑移-彎曲破壞機制的力學(xué)模型研究[J]. 巖土力學(xué), 2022, 43(增刊1): 258-266.

    [4] 黃少平, 晏鄂川, 尹曉萌, 等. 不同臨空條件的層狀反傾巖質(zhì)邊坡傾倒變形幾何特征參數(shù)影響規(guī)律[J]. 地質(zhì)科技通報, 2021, 40(1): 159-165.

    [5] LUO G, HU Q, TAN J. Process stability analysis for high slope based on limit equilibrium method and strength reduction finite element method[J]. Mining and Metallurgical Engineering, 2013, 33(2): 14-17.

    [6] YAN L, SUN Y. The analysis on rock slope stability based on strength reduction technique[J]. Journal of Guangxi University(Natural Science Edition), 2008, 33(3): 235-238.

    [7] 張文蓮, 孫曉云, 陳勇, 等. 基于巖體抗壓強度折減的邊坡穩(wěn)定性分析方法[J]. 巖土力學(xué), 2022, 43(增刊1): 607-615.

    [8] 李寧, 郭雙楓, 姚顯春. 再論巖質(zhì)高邊坡穩(wěn)定性分析方法[J]. 巖土力學(xué), 2018, 39(2): 397-406,416.

    [9] 李芬, 成濤. 基于雙折減系數(shù)法的巖質(zhì)邊坡平面滑動分析[J]. 公路工程, 2018, 43(5): 294-299.

    [10]朱擘, 何光春. 強度折減法在平面滑動型巖質(zhì)邊坡的應(yīng)用[J]. 水運工程, 2012(9): 65-69.

    [11]程小龍, 張緒進, 朱擘. 雙強度折減法在平面滑動巖質(zhì)邊坡中的應(yīng)用[J]. 水運工程, 2013(3): 72-76.

    [12]陳建宏, 鐘福生, 陳定坤. 平面滑動型巖質(zhì)邊坡極限平衡分析的上下限法[J]. 中南大學(xué)學(xué)報(自然科學(xué)版), 2013, 44(8): 3310-3315.

    [13]CHENG Y, LANSIVAARA T, WEI W. Two-dimensional slope stability analysis by limit equilibrium and strength reduction methods[J]. Computers and Geotechnics, 2007, 34(3): 137-150.

    [14]ARDESTANI A, AMINI M, ESNAEILI K. A two-dimensional limit equilibrium computer code for analysis of complex toppling slope failures[J]. Journal of Rock Mechanics and Geotechnical Engineering, 2021, 13(1): 114-130.

    [15]SCHLOTFELDT P, ELME D, PANTON B. Overhanging rock slope by design: an integrated approach using rock mass strength characterisation, large-scale numerical modelling and limit equilibrium methods-ScienceDirect[J]. Journal of Rock Mechanics and Geotechnical Engineering, 2018, 10(1): 72-90.

    [16]張崇波. 基于LHS方法的巖質(zhì)邊坡平面滑動可靠性分析[D]. 南京: 南京大學(xué), 2015.

    [17]陳志強, 王亮清, 劉順昌. 巖質(zhì)邊坡平面滑動滑面極限傾角的研究[J]. 人民長江, 2013, 44(13): 39-42.

    [18]高丙麗, 李鐸, 李朗, 等. 基于坐標投影法巖質(zhì)邊坡塊體穩(wěn)定性分析及其可視化研究[J]. 巖土力學(xué), 2022, 43(1): 181-194.

    Study on the Location and Stability of Slip Fracture Surface

    of Rocky Slope with Shear Joint

    Abstract:

    Codirectional biplane sliding is one of the common failure modes of rock slopes with a single geological fault plane (shear joint face), but the calculation methods for this type of slip surface are insufficient. In order to find the position of the slip surface of the slope more efficiently and accurately, and judge the stability of the slope, based on the limit equilibrium theory and nonlinear mathematical programming model, assuming that the slider has fallen off and the shape is an uncertain quadrilateral slide, and then assuming that the objective function is the safety factor of the rock slope, the position and stability of the slip fracture surface of the engineering slope with shear joint surface under natural working conditions are calculated by using MTALAB's global optimal search method. Then the results are analyzed and compared with those ofthe strength reduction method and the Bishop rigid body equilibrium method. The results show that the position of the slip fracture surface based on the limit equilibrium theory and the nonlinear mathematical programming model is basically consistent with the safety factor, which verifies the feasibility of such a method, providing a new basis for the calculation and stability analysis of the slip fracture surface of rock slope with a single geological fault layer.

    Key words:

    rock slope; calculation of slip surface; safety factor; limit equilibrium method; nonlinear theory; optimization principle

    猜你喜歡
    安全系數(shù)
    碎石土庫岸邊坡穩(wěn)定性及影響因素分析
    考慮材料性能分散性的航空發(fā)動機結(jié)構(gòu)安全系數(shù)確定方法
    不同因素對填筑路堤邊坡穩(wěn)定性影響分析
    反壓護道對填土路基的穩(wěn)定性影響研究
    邊坡穩(wěn)定性分析中傳遞系數(shù)法顯式解和隱式解安全系數(shù)的對應(yīng)關(guān)系探討
    基于滑動面傾角權(quán)重的邊坡穩(wěn)定性評價
    重力式擋土墻抗滑穩(wěn)定性安全系數(shù)的異性分析及經(jīng)驗安全系數(shù)方法
    閘室樁基處理后水平抗滑穩(wěn)定安全系數(shù)提高值的估算范圍研究
    電梯懸掛鋼絲繩安全系數(shù)方法的計算
    關(guān)于電梯懸掛鋼絲繩安全系數(shù)計算的一些探討
    久久这里只有精品19| 我的女老师完整版在线观看| 亚洲欧洲日产国产| 春色校园在线视频观看| 中文字幕亚洲精品专区| 日本欧美视频一区| 狠狠婷婷综合久久久久久88av| 午夜福利在线观看免费完整高清在| www.av在线官网国产| 欧美变态另类bdsm刘玥| 久久久欧美国产精品| 精品亚洲成国产av| 久久久精品区二区三区| 免费日韩欧美在线观看| 最新的欧美精品一区二区| 亚洲色图 男人天堂 中文字幕 | 国产精品免费大片| 欧美日韩国产mv在线观看视频| 国产成人aa在线观看| 日本黄色日本黄色录像| 秋霞在线观看毛片| av电影中文网址| 精品久久久精品久久久| 欧美人与性动交α欧美软件 | 日韩电影二区| 精品国产一区二区三区久久久樱花| 婷婷色综合大香蕉| 国产1区2区3区精品| 成人二区视频| 国产成人精品福利久久| 亚洲伊人久久精品综合| 一本大道久久a久久精品| 99久久人妻综合| 黄色一级大片看看| 精品国产乱码久久久久久小说| 午夜免费鲁丝| 自线自在国产av| 免费在线观看完整版高清| 亚洲精品自拍成人| 久久人人爽av亚洲精品天堂| 丰满饥渴人妻一区二区三| 国产毛片在线视频| 国产一区二区在线观看日韩| 大片免费播放器 马上看| 狠狠精品人妻久久久久久综合| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 永久免费av网站大全| 女人精品久久久久毛片| 曰老女人黄片| 精品国产一区二区三区四区第35| 天天躁夜夜躁狠狠久久av| 午夜激情久久久久久久| 国产成人aa在线观看| 午夜视频国产福利| 久久99蜜桃精品久久| 久久ye,这里只有精品| 亚洲欧美一区二区三区国产| 国产精品国产av在线观看| 国产精品 国内视频| 国产精品三级大全| 亚洲人与动物交配视频| 精品一品国产午夜福利视频| freevideosex欧美| 久久久久久久精品精品| 久久精品国产a三级三级三级| 成人亚洲精品一区在线观看| 国产成人免费观看mmmm| 免费日韩欧美在线观看| 亚洲,欧美精品.| 亚洲精品视频女| 91久久精品国产一区二区三区| 大香蕉97超碰在线| 日韩电影二区| 国产 一区精品| 日韩成人av中文字幕在线观看| 天堂俺去俺来也www色官网| 久久av网站| 国产精品成人在线| 一二三四中文在线观看免费高清| 在线观看人妻少妇| a级毛片在线看网站| 黑人猛操日本美女一级片| 少妇熟女欧美另类| 成年人免费黄色播放视频| 美女主播在线视频| 王馨瑶露胸无遮挡在线观看| 亚洲丝袜综合中文字幕| 国产一区二区三区av在线| 日本av手机在线免费观看| 在线观看www视频免费| 国产毛片在线视频| 亚洲综合色惰| 在线免费观看不下载黄p国产| 国产不卡av网站在线观看| 亚洲精品视频女| 成人18禁高潮啪啪吃奶动态图| 欧美成人午夜精品| 午夜视频国产福利| 国产一区二区在线观看日韩| 日韩欧美精品免费久久| 久久人人97超碰香蕉20202| 色94色欧美一区二区| www.色视频.com| 亚洲成人手机| 秋霞伦理黄片| 91aial.com中文字幕在线观看| 日韩成人av中文字幕在线观看| 狠狠精品人妻久久久久久综合| 成人亚洲欧美一区二区av| 成人国语在线视频| 99热6这里只有精品| 亚洲精品一二三| 超色免费av| 在线观看国产h片| 热re99久久精品国产66热6| 97超碰精品成人国产| 成人手机av| 国产av国产精品国产| 久久人人97超碰香蕉20202| 国产精品久久久久久久电影| av一本久久久久| 建设人人有责人人尽责人人享有的| 亚洲,欧美,日韩| 少妇被粗大猛烈的视频| 国产精品嫩草影院av在线观看| h视频一区二区三区| 婷婷色麻豆天堂久久| 久久精品久久精品一区二区三区| 涩涩av久久男人的天堂| 日韩av不卡免费在线播放| 成人18禁高潮啪啪吃奶动态图| 丝袜在线中文字幕| 久久人妻熟女aⅴ| 激情五月婷婷亚洲| 免费黄网站久久成人精品| 99视频精品全部免费 在线| 国产一区亚洲一区在线观看| 久久精品夜色国产| 久久久久久久久久久免费av| 高清av免费在线| 欧美成人午夜精品| 国产高清不卡午夜福利| 欧美少妇被猛烈插入视频| 美女视频免费永久观看网站| 亚洲中文av在线| 视频区图区小说| 国产在视频线精品| 久久精品熟女亚洲av麻豆精品| a级毛片在线看网站| 久久99热这里只频精品6学生| 丰满饥渴人妻一区二区三| 国产毛片在线视频| 街头女战士在线观看网站| av在线观看视频网站免费| 欧美日韩精品成人综合77777| 五月玫瑰六月丁香| 国产高清国产精品国产三级| 亚洲成人手机| 欧美人与性动交α欧美软件 | 中文欧美无线码| 国产亚洲午夜精品一区二区久久| 好男人视频免费观看在线| 日韩一区二区三区影片| 欧美精品高潮呻吟av久久| 十八禁网站网址无遮挡| 久久久久久久久久久久大奶| 人人妻人人澡人人看| 久久久久人妻精品一区果冻| 日本欧美国产在线视频| 黑人巨大精品欧美一区二区蜜桃 | 国产视频首页在线观看| 精品久久久精品久久久| 99re6热这里在线精品视频| 日韩不卡一区二区三区视频在线| 久久精品夜色国产| 色94色欧美一区二区| 欧美精品国产亚洲| 中文字幕亚洲精品专区| 在线看a的网站| 国产片内射在线| 婷婷成人精品国产| 午夜福利网站1000一区二区三区| 久久 成人 亚洲| 亚洲国产看品久久| 日本av免费视频播放| av卡一久久| 老女人水多毛片| 国产永久视频网站| av国产久精品久网站免费入址| 国产成人精品一,二区| 免费高清在线观看视频在线观看| 99久久精品国产国产毛片| 午夜视频国产福利| 18禁裸乳无遮挡动漫免费视频| 97精品久久久久久久久久精品| 不卡视频在线观看欧美| 一级片'在线观看视频| 亚洲成人一二三区av| 国产av精品麻豆| 久久久精品94久久精品| 涩涩av久久男人的天堂| 少妇被粗大猛烈的视频| 麻豆精品久久久久久蜜桃| 人人妻人人澡人人爽人人夜夜| 大片电影免费在线观看免费| av黄色大香蕉| 26uuu在线亚洲综合色| 国产精品久久久久成人av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日韩三级伦理在线观看| 女的被弄到高潮叫床怎么办| av.在线天堂| 亚洲国产精品一区三区| 日韩精品免费视频一区二区三区 | 亚洲精品成人av观看孕妇| 免费高清在线观看日韩| 色视频在线一区二区三区| 免费少妇av软件| 国产精品无大码| 中文字幕精品免费在线观看视频 | 久久精品国产综合久久久 | av卡一久久| 热99久久久久精品小说推荐| 各种免费的搞黄视频| 欧美成人午夜精品| 纵有疾风起免费观看全集完整版| 最后的刺客免费高清国语| 亚洲精华国产精华液的使用体验| 亚洲国产av新网站| 国产精品蜜桃在线观看| 亚洲精品,欧美精品| 欧美 亚洲 国产 日韩一| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲av.av天堂| 国产乱来视频区| 伦理电影大哥的女人| 色吧在线观看| 天天躁夜夜躁狠狠躁躁| 日韩av不卡免费在线播放| av国产久精品久网站免费入址| 亚洲,欧美精品.| 亚洲国产欧美在线一区| 韩国高清视频一区二区三区| 亚洲av欧美aⅴ国产| 免费大片18禁| 考比视频在线观看| 精品午夜福利在线看| 国产成人欧美| 日韩av不卡免费在线播放| 亚洲国产最新在线播放| 日韩成人伦理影院| 国产精品免费大片| videos熟女内射| 街头女战士在线观看网站| 国产片特级美女逼逼视频| 人人妻人人澡人人爽人人夜夜| 久久影院123| 精品亚洲成国产av| 菩萨蛮人人尽说江南好唐韦庄| 一边摸一边做爽爽视频免费| 午夜福利网站1000一区二区三区| 一区二区三区乱码不卡18| 午夜福利,免费看| 久久久国产一区二区| 日日啪夜夜爽| 日本wwww免费看| 日韩一区二区三区影片| 亚洲激情五月婷婷啪啪| 久久久久网色| 狠狠婷婷综合久久久久久88av| 亚洲,一卡二卡三卡| 亚洲精华国产精华液的使用体验| 欧美国产精品一级二级三级| 国产成人精品一,二区| 搡女人真爽免费视频火全软件| 久久精品久久久久久久性| 大话2 男鬼变身卡| 一本久久精品| 99视频精品全部免费 在线| 久久精品久久精品一区二区三区| 亚洲,欧美,日韩| 黑人猛操日本美女一级片| 少妇精品久久久久久久| 极品人妻少妇av视频| 国产精品久久久久久精品古装| 天天操日日干夜夜撸| 一级毛片我不卡| 亚洲,欧美,日韩| 久久久精品94久久精品| 两个人免费观看高清视频| 日日撸夜夜添| 国产在线一区二区三区精| 欧美人与性动交α欧美精品济南到 | 欧美精品亚洲一区二区| 日本-黄色视频高清免费观看| 看免费av毛片| 伦精品一区二区三区| 久久久精品94久久精品| 国产av一区二区精品久久| 一级毛片 在线播放| 久久久久久久国产电影| 99热国产这里只有精品6| 欧美3d第一页| 欧美日韩成人在线一区二区| 青春草国产在线视频| 男女下面插进去视频免费观看 | 国产精品 国内视频| 久久国内精品自在自线图片| 肉色欧美久久久久久久蜜桃| 9191精品国产免费久久| 黑丝袜美女国产一区| 中文字幕av电影在线播放| 久久久久久久久久成人| 观看美女的网站| 国产精品不卡视频一区二区| 青春草亚洲视频在线观看| 国产淫语在线视频| 午夜久久久在线观看| 亚洲国产精品国产精品| 亚洲av成人精品一二三区| 飞空精品影院首页| av福利片在线| 国产一级毛片在线| 99香蕉大伊视频| 建设人人有责人人尽责人人享有的| 天天躁夜夜躁狠狠久久av| 2021少妇久久久久久久久久久| 国产精品一区二区在线观看99| 青春草亚洲视频在线观看| 亚洲人成网站在线观看播放| 久久久精品94久久精品| 午夜视频国产福利| 男女高潮啪啪啪动态图| 在线观看免费视频网站a站| 十分钟在线观看高清视频www| 晚上一个人看的免费电影| 1024视频免费在线观看| 丝袜脚勾引网站| 一级毛片黄色毛片免费观看视频| 日日摸夜夜添夜夜爱| 尾随美女入室| 久久综合国产亚洲精品| 亚洲精华国产精华液的使用体验| 国产高清国产精品国产三级| 热99久久久久精品小说推荐| 欧美精品av麻豆av| 国产av一区二区精品久久| 久久人人爽人人爽人人片va| 久久久久久久亚洲中文字幕| 亚洲国产av影院在线观看| 午夜精品国产一区二区电影| 黄色毛片三级朝国网站| 国产黄色免费在线视频| 国内精品宾馆在线| 亚洲综合精品二区| 亚洲国产欧美日韩在线播放| 热99久久久久精品小说推荐| 中文欧美无线码| 有码 亚洲区| 午夜91福利影院| 久久久久久久亚洲中文字幕| 伊人亚洲综合成人网| 亚洲国产看品久久| 菩萨蛮人人尽说江南好唐韦庄| 色哟哟·www| 少妇熟女欧美另类| 91精品三级在线观看| 久久99精品国语久久久| 午夜久久久在线观看| 曰老女人黄片| 亚洲第一av免费看| 看十八女毛片水多多多| 精品一区二区三区四区五区乱码 | 国产成人精品在线电影| 制服诱惑二区| 免费观看性生交大片5| 高清毛片免费看| 欧美日韩综合久久久久久| 黑人高潮一二区| 在线 av 中文字幕| 国产色婷婷99| 18在线观看网站| 亚洲,一卡二卡三卡| 一个人免费看片子| 美女大奶头黄色视频| √禁漫天堂资源中文www| 99热这里只有是精品在线观看| 少妇高潮的动态图| 午夜久久久在线观看| 晚上一个人看的免费电影| 成人亚洲欧美一区二区av| 国产欧美另类精品又又久久亚洲欧美| 国产精品一区www在线观看| 亚洲在久久综合| 亚洲国产精品专区欧美| 亚洲精品日本国产第一区| 亚洲少妇的诱惑av| 免费女性裸体啪啪无遮挡网站| 多毛熟女@视频| 三上悠亚av全集在线观看| 中文字幕人妻丝袜制服| 97在线视频观看| 亚洲国产欧美日韩在线播放| 亚洲中文av在线| 日韩中文字幕视频在线看片| www日本在线高清视频| 97在线人人人人妻| 边亲边吃奶的免费视频| 久久 成人 亚洲| 青春草视频在线免费观看| 18在线观看网站| 丰满少妇做爰视频| 男人添女人高潮全过程视频| 18在线观看网站| 男男h啪啪无遮挡| 亚洲丝袜综合中文字幕| 超色免费av| 好男人视频免费观看在线| 夫妻性生交免费视频一级片| 国产精品三级大全| 亚洲精品国产av蜜桃| 大话2 男鬼变身卡| 久久ye,这里只有精品| 国产一级毛片在线| 国产精品久久久久久久电影| 大片电影免费在线观看免费| 婷婷成人精品国产| 国产有黄有色有爽视频| 国产激情久久老熟女| 建设人人有责人人尽责人人享有的| 色94色欧美一区二区| 高清av免费在线| 欧美日韩亚洲高清精品| 精品国产一区二区三区久久久樱花| xxx大片免费视频| 另类精品久久| 男女免费视频国产| 亚洲 欧美一区二区三区| 国产一区有黄有色的免费视频| 欧美日韩精品成人综合77777| 狂野欧美激情性bbbbbb| 久久久久久久久久久免费av| 亚洲国产欧美在线一区| 成人二区视频| 男女免费视频国产| 啦啦啦视频在线资源免费观看| 午夜免费观看性视频| 成人国产麻豆网| 国产精品嫩草影院av在线观看| 日韩成人伦理影院| 少妇的丰满在线观看| 黄网站色视频无遮挡免费观看| 国产精品一区二区在线观看99| 男女边吃奶边做爰视频| 性色av一级| 男女免费视频国产| 少妇猛男粗大的猛烈进出视频| 亚洲国产毛片av蜜桃av| 欧美精品一区二区免费开放| 久久久久久久大尺度免费视频| 亚洲精品456在线播放app| 免费大片黄手机在线观看| 一本—道久久a久久精品蜜桃钙片| 曰老女人黄片| 尾随美女入室| 久久久久视频综合| 欧美 日韩 精品 国产| 交换朋友夫妻互换小说| 国产深夜福利视频在线观看| 亚洲久久久国产精品| 精品国产一区二区久久| 亚洲人与动物交配视频| 国产欧美亚洲国产| 99热网站在线观看| 老司机亚洲免费影院| 久久久久久久久久成人| 日韩av免费高清视频| 纯流量卡能插随身wifi吗| 亚洲一区二区三区欧美精品| 卡戴珊不雅视频在线播放| 男人添女人高潮全过程视频| 91aial.com中文字幕在线观看| 高清不卡的av网站| 考比视频在线观看| 国产成人91sexporn| 秋霞在线观看毛片| 90打野战视频偷拍视频| 爱豆传媒免费全集在线观看| 青春草国产在线视频| 日本黄色日本黄色录像| 菩萨蛮人人尽说江南好唐韦庄| 婷婷色av中文字幕| 国产成人精品无人区| 91精品国产国语对白视频| 成年女人在线观看亚洲视频| 欧美人与性动交α欧美精品济南到 | 亚洲伊人久久精品综合| 九色亚洲精品在线播放| 熟女电影av网| 日韩中文字幕视频在线看片| 嫩草影视91久久| 免费看十八禁软件| 国产精品国产高清国产av | 黑人猛操日本美女一级片| 午夜免费成人在线视频| 黄色片一级片一级黄色片| 国产一区二区三区视频了| 免费看十八禁软件| 亚洲国产精品合色在线| 亚洲性夜色夜夜综合| 精品国内亚洲2022精品成人 | 男女下面插进去视频免费观看| 久久久久久亚洲精品国产蜜桃av| 变态另类成人亚洲欧美熟女 | 最近最新中文字幕大全电影3 | 大型黄色视频在线免费观看| 青草久久国产| 色老头精品视频在线观看| 在线看a的网站| 精品久久久精品久久久| 99精国产麻豆久久婷婷| 亚洲成人手机| 欧美日本中文国产一区发布| 久久久久久久久免费视频了| 久久精品aⅴ一区二区三区四区| 人妻久久中文字幕网| 一区二区三区激情视频| 国产成人精品无人区| 国产精品乱码一区二三区的特点 | 在线十欧美十亚洲十日本专区| 岛国毛片在线播放| 夜夜躁狠狠躁天天躁| 精品人妻在线不人妻| 老司机福利观看| 啦啦啦免费观看视频1| av网站免费在线观看视频| 欧美日韩瑟瑟在线播放| 精品国产一区二区三区久久久樱花| 色播在线永久视频| 精品人妻熟女毛片av久久网站| 欧洲精品卡2卡3卡4卡5卡区| 国产日韩欧美亚洲二区| 热99re8久久精品国产| 少妇粗大呻吟视频| 最新的欧美精品一区二区| 国产亚洲精品一区二区www | 久久人妻av系列| 热99久久久久精品小说推荐| 午夜久久久在线观看| av在线播放免费不卡| 婷婷丁香在线五月| 国产又色又爽无遮挡免费看| 日本vs欧美在线观看视频| 亚洲伊人色综图| e午夜精品久久久久久久| av天堂在线播放| 99久久精品国产亚洲精品| 乱人伦中国视频| 亚洲色图 男人天堂 中文字幕| 天堂俺去俺来也www色官网| 亚洲人成电影观看| 老司机在亚洲福利影院| 69av精品久久久久久| 久久精品亚洲熟妇少妇任你| 国产精品av久久久久免费| 咕卡用的链子| 亚洲精品美女久久久久99蜜臀| 最近最新中文字幕大全电影3 | 国产亚洲一区二区精品| 国产av又大| 水蜜桃什么品种好| 亚洲一区二区三区欧美精品| 国产视频一区二区在线看| 国产精品偷伦视频观看了| 不卡一级毛片| 国产国语露脸激情在线看| 老司机靠b影院| 视频区图区小说| 国产男靠女视频免费网站| 亚洲免费av在线视频| 久久久国产欧美日韩av| av片东京热男人的天堂| 黑人巨大精品欧美一区二区蜜桃| 99热网站在线观看| 国产又爽黄色视频| 亚洲国产中文字幕在线视频| 久久天堂一区二区三区四区| 人人妻人人添人人爽欧美一区卜| 99re6热这里在线精品视频| 欧美最黄视频在线播放免费 | 亚洲伊人色综图| 亚洲一区高清亚洲精品| 国产又爽黄色视频| 亚洲欧洲精品一区二区精品久久久| 美女高潮喷水抽搐中文字幕| 脱女人内裤的视频| 国产成人精品在线电影| 免费在线观看视频国产中文字幕亚洲| а√天堂www在线а√下载 | 麻豆av在线久日| 久久久精品区二区三区| 国产极品粉嫩免费观看在线| 天天躁狠狠躁夜夜躁狠狠躁| 国产男女内射视频| 国产蜜桃级精品一区二区三区 | 在线永久观看黄色视频| av天堂久久9| 欧美激情极品国产一区二区三区| 美国免费a级毛片| 老鸭窝网址在线观看| 一边摸一边做爽爽视频免费| 在线观看日韩欧美| 极品少妇高潮喷水抽搐| 极品人妻少妇av视频| 欧美黑人精品巨大|