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

    基于單元重構與節(jié)點分離的大型地下洞室軟弱結構面模擬方法

    2017-04-14 03:47:21趙健肖明陳俊濤李冬冬
    湖南大學學報·自然科學版 2017年3期

    趙健+肖明+陳俊濤+李冬冬

    摘 要:針對在大型地下洞室中軟弱結構面難以合理模擬的問題,通過引進“單元重構節(jié)點分離”思想,提出了新的非連續(xù)軟弱結構面模擬方法.將該方法用于實際地下廠房模型中并模擬出非連續(xù)結構面,并將該重構模型導入計算程序,與常規(guī)連續(xù)有限元模型的計算結果進行對比.結果表明,重構模型能夠用于數(shù)值計算,計算結果在非連續(xù)結構面附近呈現(xiàn)非連續(xù)性,且計算結果所反映的規(guī)律與所認識規(guī)律基本一致,能夠反映真實的結構面賦存狀態(tài).因此,該軟弱結構面建模方法實現(xiàn)了結構面的非連續(xù)化,為大型地下洞室的復雜地質非連續(xù)結構面模擬提供了新的實現(xiàn)途徑.

    關鍵詞:單元重構;節(jié)點分離;結構面;非連續(xù);建模方法

    中圖分類號:TU45 文獻標志碼:A

    文章編號:1674-2974(2017)03-0134-09DOI:10.16339/j.cnki.hdxbzkb.2017.03.017

    Abstract:As weak structural planes are difficult to be simulated reasonably in large underground chamber, “element reconstruction and node separation” was introduced, and the new simulation methodology of discontinuous weak structural planes was proposed in this study. The proposed methodology was used in the actual numerical model of underground powerhouse to simulate the discontinuous weak structural planes. The reconstructed model was then imported into calculation program, and the result was compared with the prediction by conventional continuous finite element model. The comparison shows that the reconstructed model can be used in the numerical calculation, and the discontinuity feature is presented near the discontinuous structural planes. Moreover, the rules in the calculation results are consistent with the general rules, which can describe the actual occurrence mode of the structural planes. Therefore, the simulation methodology of discontinuous weak structural planes achieves the transformation into the discontinuity of the structural planes, and it provides a new way for the simulation of complex geological discontinuous weak structural planes in large underground chamber.

    Key words:element reconstruction;node separation;structural planes;discontinuity;modelling

    軟弱結構面是力學強度明顯低于圍巖且一般充填有一定厚度軟弱物質的地質界面.在地下洞室開挖與施工中,軟弱結構面是控制圍巖穩(wěn)定的關鍵因素,它直接制約著工程巖體變形、破壞的發(fā)生和發(fā)展過程,巖體的失穩(wěn)也往往發(fā)生在這類巖體強度最薄弱的部位[1].同時大型地下洞室往往交錯縱橫,結構布置復雜,軟弱結構面與洞室任意切割,在建模前期直接建立結構面單元也較為復雜.因此軟弱結構面的正確模擬是分析地下洞室圍巖穩(wěn)定的關鍵,是確保圍巖穩(wěn)定計算結果合理的前提條件.

    目前針對軟弱結構面的模擬,主要有以下幾種常見的方法:①單元弱化.將結構面穿過的有限單元的力學參數(shù)降低而不建立單獨的斷層單元,對其進行單獨賦值處理,以此將結構面單元和周圍巖體區(qū)分開來.這種方法計算思路較為簡便,且容易操作,但是容易得到齒狀交錯的斷層單元,與實際工程狀態(tài)有差異.②長而薄的有限單元[2].直接在計算模型中依據(jù)軟弱結構面產(chǎn)狀建立有限單元模擬結構面,并用有限元方法進行分析.這種方法能夠得到較為準確的、符合實際的平面斷層單元,但是前期建模較為復雜,建模效率較低.③隱含斷層[3].肖明等采用隱含斷層單元的數(shù)值模擬計算方法將復雜斷層單元隱含在巖體單元中,用一個均質各向異性的巖體當量復合單元模擬斷層的結構效應.這種方法能夠將斷層效果耦合到巖體單元中,力學性質簡明,但是需要重新編制有限元程序,不能直接應用到現(xiàn)行的有限元計算程序和軟件中.總結方法①②③,可發(fā)現(xiàn)目前這幾種常見的結構面模擬方法仍然基于變形連續(xù)的有限元分析,不能反映由斷層所帶來的非連續(xù)特征,結構面兩側巖體仍然連續(xù),沒有斷開,與實際情況有出入[4].除了以上幾種方法,應用比較廣泛的還有Goodman單元[5].Goodman將其提出并應用于巖石力學中作為節(jié)理單元,后來被推廣應用于土與結構的共同作用、人工塊體結構的有限元計算中作為接觸面單元.這種方法能夠模擬結構面的非連續(xù)變形特性,但是單元無厚度,結構面兩側單元容易相互嵌入.

    在總結以往研究成果的基礎上,本文引進“單元重構節(jié)點分離”思想,改進了軟弱結構面模擬方法.在不考慮結構面的有限元模型的基礎上引入結構面,將結構面穿過的有限單元進行切割重構并生成新單元,同時將結構面兩側共用節(jié)點進行分離.這種方法不但使大型地下洞室群中復雜軟弱結構面建模簡單化,而且能夠實現(xiàn)結構面的非連續(xù)化,為含有復雜結構面的地下洞室圍巖穩(wěn)定分析提供了一種十分有效的新思路.

    1 軟弱結構面模擬過程

    大型地下洞室軟弱結構面的建??梢詺w納為以下5個步驟, 如圖1所示.

    1)不考慮結構面完成單元離散.根據(jù)地下結構的布置,考慮各管道和洞室之間的相對位置,建立不考慮結構面的有限元模型,完成單元離散.這一部分的發(fā)展已經(jīng)比較成熟,許多有限元軟件都可以完成模型的建立,并且不同軟件建立的模型也可以通過編制程序完成相互轉換.

    2)確定結構面的信息.根據(jù)不同的發(fā)育程度可將結構面劃分成Ⅰ~Ⅴ級,發(fā)育程度從區(qū)域性的斷裂帶逐步減小到微小的節(jié)理裂隙,由此帶來了結構面模擬方法的不同.根據(jù)結構面相對于地下廠房的延伸范圍,可以分為有限結構面和無限結構面;根據(jù)結構面的分布情況,可以分為確定結構面和隨機結構面.目前很多學者將結構面模擬為具有圓盤形狀的平面,同時通過Monte Carlo隨機方法[6]也可以生成具有裂隙特征的結構面信息.

    3)采用單元重構技術,完成結構面導入.由于地下洞室存在的結構面數(shù)量往往不止一條,因此可以采用分次導入的方式.在不考慮結構面的有限元模型的基礎上,導入第一條結構面切割單元并生成重構模型,之后依次在之前的重構模型的基礎上,導入剩余的結構面.同時值得注意的是,在第2)步的結構面信息確定中,當某條結構面發(fā)育程度較高,結構面間空隙較大且軟弱填充物較多時,可以確定結構面間隙的相對距離,并通過平移結構面,對有限元模型進行兩次切割,從而模擬出結構面之間空隙的效果,中間的狹長單元可視為軟弱填充物.

    4)采用節(jié)點分離技術,完成結構面非連續(xù)化.記錄在單元重構中新生成的節(jié)點,并通過添加重復節(jié)點完成結構面的非連續(xù)化.該步驟是區(qū)分重構模型與傳統(tǒng)有限元連續(xù)模型的最根本的技術.在有限元模型變形之前,結構面上重復節(jié)點的空間坐標相同,重復節(jié)點重合.當結構面兩側單元發(fā)生相對變形時,重復節(jié)點的空間位置發(fā)生相對錯動,結構面脫開分離.

    5)檢查重構單元的幾何信息.在完成單元重構和節(jié)點分離后,應重點檢查重構單元的幾何信息,防止出現(xiàn)一些幾何形態(tài)畸形的單元.可通過計算單元的雅克比矩陣或單元體積來判斷單元是否畸形.

    以下重點介紹在模擬過程中涉及到的“單元重構”和“節(jié)點分離”技術.

    2 單元重構與節(jié)點分離

    2.1 單元重構

    關于單元重構方法的介紹詳見文獻[7-8],其基本思想是對結構面穿過的單元進行二次劃分,生成新單元.首先,采用成熟的網(wǎng)格劃分技術,對研究區(qū)域進行單元離散,重點研究區(qū)域和開挖臨空面附近的單元劃分相對于其他區(qū)域應該較小,以滿足研究精度的需要,以此生成不考慮結構面的計算模型;其次,將結構面依次導入計算模型,對結構面穿過的單元進行二次劃分,生成滿足計算要求的新單元.

    為了形象具體地說明該算法,這里給出一個簡單易懂的二維實例.首先對該區(qū)域進行網(wǎng)格離散,生成4單元9節(jié)點有限元模型,如圖2所示;然后將結構面1和結構面2導入,對結構面穿過單元進行一次單元重構,單元E1被重新劃分為單元E1和E2,單元E2被重新劃分為單元E3和E4,單元E3被重新劃分為單元E5,E6,E7和E8,單元E4被重新劃分為單元E9和E10,如圖3所示.

    將結構面導入后,在結構面附近會生成許多形態(tài)不規(guī)則單元,例如圖3中的五邊形單元E1和E4.目前在二維有限元計算中,比較常見的單元形態(tài)為四面體單元和三角形單元,這2種單元形態(tài)簡單,適應性良好且能夠滿足工程計算精度需要,而其他形態(tài)的單元在有限元計算中往往比較復雜,其形函數(shù)很難構造[9-10],因此對這部分單元需要進行二次單元重構.

    將圖3中的單元E1和E4添加網(wǎng)格線進行二次單元重構,得到如圖4中所示結果.五邊形單元E1被重新劃分為三角形單元E1和四邊形單元E2,五邊形單元E4被重新劃分為四邊形單元E5和三角形單元E6.至此,圖2中所有單元被重新切割劃分為三角形和四邊形單元,重構生成的單元都能進行有限元分析計算.

    同理在三維分析中,比較常見的計算單元為八節(jié)點六面體單元和四面體單元,其他的復雜形態(tài)單元也可通過添加網(wǎng)格面等方式重新劃分為四面體單元[11],圖5和圖6顯示了五節(jié)點四棱錐和六節(jié)點三棱柱的重構過程,其他復雜單元也可參考下圖重構.

    結構面導入模型后,由于只在結構面附近進行單元重構,離結構面較遠的大部分圍巖單元仍然保持原單元形態(tài),因此需要進行重構的單元數(shù)量相對整體模型來說較少,計算量不大.且單元重構完成后模型中只存在八節(jié)點六面體和四面體2種單元形態(tài),這2種單元形函數(shù)表達方式都比較簡單[9],且能夠滿足工程計算精度要求.

    采用該單元重構方法,能夠減小含有結構面模型的建模難度,可以在不考慮結構面的有限元模型基礎上任意導入結構面.該重構模型能夠保留在關鍵研究部位和開挖臨空面附近的細小劃分單元,同時能夠獲得結構面附近的重構單元,這一部分單元往往劃分尺寸較小,能夠滿足研究精度要求.但是用此方法重新生成的計算模型仍然是連續(xù)的,結構面兩側巖體仍然共用節(jié)點,沒有實現(xiàn)節(jié)點分離,計算結果不能較好地呈現(xiàn)非連續(xù)特征.

    2.2 節(jié)點分離

    由于采用單元重構技術后的計算模型仍然是連續(xù)的,因此需要進行節(jié)點分離,其基本步驟如下:

    1)遍歷模型中的所有結構面,存儲每個結構面附近的所有節(jié)點.

    2)查找步驟1)中存儲節(jié)點所在的所有單元,并分析是否存在單元分布于結構面兩側.若存在,則記錄該共用節(jié)點和單元.

    3)計算共用節(jié)點所涉及的單元數(shù)量,在相應的單元增加重復節(jié)點,該重復節(jié)點的空間坐標在發(fā)生相對位移前相同.

    以下給出了一個二維節(jié)點分離的簡單示例,如圖7所示.結構面附近的節(jié)點編號為3, 5, 7, 10, 11, 12, 13和15,將其進行節(jié)點分離.除了節(jié)點11外,其余7個節(jié)點都分離為2個節(jié)點,由于節(jié)點11處在2個結構面的交界處,因此節(jié)點11分離為4個節(jié)點.

    通過增加共用節(jié)點,完成了節(jié)點分離,實現(xiàn)了結構面的非連續(xù)化.

    3 地下洞室斷層模擬

    本節(jié)以鎮(zhèn)安水電站地下廠房為例,將單元重構與節(jié)點分離技術應用到地下洞室斷層模擬中.

    鎮(zhèn)安水電站地下廠房模型高為622 m,寬為300 m,地下洞室群埋深約為210 m,根據(jù)實測地形資料模擬山體表面起伏.鎮(zhèn)安水電站地質勘探資料見表1,地下廠房區(qū)域斷層發(fā)育較為充分,主要存在表1所示2組斷層,而其他裂隙節(jié)理等發(fā)育不是非常明顯且數(shù)量不多,因此主要對以斷層為主的結構面進行模擬.

    以廠房1#機組為坐標原點建立局部坐標系,以廠房軸向指向4#機組段方向為y方向,以垂直y方向指向下游為x方向.

    將斷層模擬成具有一定范圍的圓盤面[12],由斷層產(chǎn)狀和斷層圓心點可唯一確定斷層的空間位置,同時由斷層半徑可確定斷層的延伸范圍.由于斷層間裂隙較明顯,厚度較大,因此對該模型中的斷層采用2個平行面進行切割,模擬出斷層厚度,中間切割重構而成的單元可看成是斷層的軟弱夾泥.

    根據(jù)表1中的斷層產(chǎn)狀,將斷層導入有限元模型中,經(jīng)過單元重構和節(jié)點分離,完成了斷層的模擬,如圖8,圖9和圖10所示.從圖中可看出,采用“單元重構節(jié)點分離”方法,能夠導入任意產(chǎn)狀結構面,減小了前期結構面建模的難度,且結構面附近單元相對密集,重構單元形態(tài)良好,能夠滿足工程需要.

    4 考慮軟弱結構面的重構模型數(shù)值計算

    為驗證所建重構模型能夠用于數(shù)值計算,將其采用NDDA方法(基于節(jié)點的非連續(xù)變形分析)[13-14]對第3節(jié)中考慮結構面的鎮(zhèn)安地下廠房模型進行分析,并與傳統(tǒng)連續(xù)斷層模型的有限元計算結果進行對比.傳統(tǒng)連續(xù)斷層采用真實存在的長而薄的薄層單元來模擬,地下洞室?guī)r體力學參數(shù)見表2.

    從圖9(b)可看出,2條斷層都穿過2#機組段,2#機組段受斷層切割作用明顯,因此選擇2#機組段作為典型機組段進行分析.計算結果如圖11和圖12所示,分別表示了2#機組段采用連續(xù)斷層和非連續(xù)斷層的洞周位移和最大主應力分布.規(guī)定壓應力為負,拉應力為正.

    當采用連續(xù)斷層時,斷層附近的應力值有所降低,斷層附近應力云圖沿斷層向下突出,意味著應力值由正常范圍向斷層帶降低[15].其計算結果往往存在一個“影響帶”,該“影響帶”區(qū)域往往大于斷層本身厚度.由于結構面的作用,該區(qū)域位移增大,應力強度降低,但仍然保持連續(xù)變化,該規(guī)律從圖11(b)和圖12(b)可以清晰地看出.而當采用非連續(xù)斷層后,斷層兩邊的位移和應力狀態(tài)往往不連續(xù),兩側差異明顯,出現(xiàn)齒狀和交錯的現(xiàn)象[16].

    為了更清晰、確切地觀察連續(xù)斷層和非連續(xù)斷層對洞室穩(wěn)定影響的區(qū)別,選取受斷層切割的主廠房上游邊墻底部作為研究對象.主廠房上游邊墻中下部高邊墻效應明顯,在洞室圍巖穩(wěn)定分析中是重點研究部位.

    從圖13和圖14可看出,采用非連續(xù)斷層時,斷層單元的位移和應力狀態(tài)與兩側巖體的位移和應力狀態(tài)完全不同,呈現(xiàn)出非連續(xù)變化[17],斷層與兩側巖體的位移和應力值相比相差較多,且斷層和兩側巖體出現(xiàn)錯開現(xiàn)象(如圖13(a)所示),而采用連續(xù)斷層時,斷層與兩側巖體的位移和應力值連續(xù)變化,沒有出現(xiàn)位移和應力值的不連續(xù)跳躍現(xiàn)象.從數(shù)值上看,非連續(xù)斷層和連續(xù)斷層相比,斷層位移增幅為0.3~0.5 cm,斷層最大主應力降幅在2.0 MPa左右.同時可觀察到,在采用2種不同模擬斷層情況下,“影響帶”范圍之外的巖體位移和應力狀態(tài)的變化規(guī)律基本一致,沒有較大差別.

    可以看出,考慮非連續(xù)斷層后,洞周圍巖變形呈現(xiàn)不連續(xù)性,計算結果所反映的規(guī)律與所認識的經(jīng)驗規(guī)律一致,表明含有非連續(xù)斷層的重構模型可以用于數(shù)值計算,能夠反映非連續(xù)斷層對于圍巖穩(wěn)定的影響.

    5 結 論

    1)基于單元重構提出了以增加共用節(jié)點為核心的節(jié)點分離技術,完成了軟弱結構面的模擬.該方法的突出優(yōu)勢在于能夠在傳統(tǒng)有限元模型的基礎上任意導入軟弱結構面,而非前期直接模擬結構面,很大程度上減小了結構面模擬的難度,并且實現(xiàn)了模型從連續(xù)到不連續(xù)的轉換.相比于傳統(tǒng)的結構面建模方法,該方法能夠合理地模擬模型中結構面位置,正確地反映結構面的實際空間分布,形成非連續(xù)結構面,更符合工程實際,是一種全新的、精確的結構面建模方法.

    2)以地下洞室為例,將提出的結構面建模方法應用于有限元模型中,并采用NDDA方法進行數(shù)值計算,與連續(xù)有限元模型的計算結果進行對比.結果表明,當采用非連續(xù)結構面時,結構面與周圍巖體的位移和應力狀態(tài)呈現(xiàn)非連續(xù)變化規(guī)律,其數(shù)值結果和連續(xù)結構面相比也有變化.計算結果所反映的規(guī)律與所認識規(guī)律基本一致,更接近工程實際,證明了該非連續(xù)模型用于實際工程的可靠性.

    由此可以看出,以“單元重構節(jié)點分離”為核心的地下廠房軟弱結構面建模方法是一種全新的結構面模擬方法,該方法實現(xiàn)了結構面的非連續(xù)化,并為結構面的非連續(xù)變形分析提供了新的實現(xiàn)途徑.

    參考文獻

    [1] 康石磊, 楊峰, 張箭,等. 基于強度折減和上限有限元的橢圓形毛洞隧道圍巖穩(wěn)定性分析[J]. 湖南大學學報:自然科學版, 2015, 42(9):104-109.

    KANG Shilei, YANG Feng, ZHANG Jian, et al. Finite element upper bound analysis of stability of unlined elliptical tunnel based on strength reduction method [J]. Journal of Hunan University: Natural Sciences, 2015, 42(9): 104-109.(In Chinese)

    [2] 邱穎, 孫云貴. 斷層力學模型及其數(shù)值模擬方法評述[J]. 水利與建筑工程學報, 2004, 2(1):14-16.

    QIU Ying, SUN Yungui. Study on fault mechanics model and numerical simulation methods[J]. Journal of Water Resources & Architectural Engineering, 2004, 2(1):14-16.(In Chinese)

    [3] 肖明, 陳俊濤. 大型地下洞室復雜地質斷層數(shù)值模擬分析方法[J]. 巖土力學, 2006, 27(6):880-884.

    XIAO Ming, CHEN Juntao. Numerical simulation and analysis method of complex geologic faults in large underground chamber[J]. Rock & Soil Mechanics, 2006, 27(6):880-884. (In Chinese)

    [4] 雷顯權, 陳運平. 應用非連續(xù)模型研究斷層對地殼應力的影響[J]. 中南大學學報:自然科學版, 2011, 42(8):2379-2386.

    LEI Xianquan, CHEN Yunping. Effects of faults on crustal stress using discontinuous model[J]. Journal of Central South University: Science and Technology, 2011, 42(8): 2379-2386.(In Chinese)

    [5] GOODMAN R E, TAYLOR R L, BREKKE T L. A model for the mechanics of jointed rock[J]. Journal of Soil Mechanics & Foundations Div, 2015, 94:637-660.

    [6] 方濤, 柴軍瑞, 胡海浪, 等. Monte Carlo 方法在巖體裂隙結構面模擬中的應用[J]. 露天采礦技術, 2007(1): 7-9.

    FANG Tao, CHAI Junrui, HU Hailang, et al. Application of Monte Carlo method to simulating the fracture interconnectivity of rock mass[J]. Opencast Mining Technology, 2007(1):7-9.(In Chinese)

    [7] 張雨霆, 肖明, 左雙英. 基于單元重構的巖土工程復雜地質斷層建模方法[J]. 巖石力學與工程學報, 2009, 28(9):1848-1855.

    ZHANG Yuting, XIAO Ming, ZUO Shuangying. Methodology for modeling of complex geological faults in geotechnical engineering based on element reconstruction[J]. Chinese Journal of Rock Mechanics & Engineering, 2009, 28(9):1848-1855.(In Chinese)

    [8] 張雨霆, 肖明, 丁秀麗,等. 復雜巖石塊體識別的單元重構-聚合方法[J]. 巖石力學與工程學報, 2012, 31(3):507-523.

    ZHANG Yuting, XIAO Ming, DING Xiuli, et al. Identification method of complex rock blocks using element reconstruction and aggregation technique[J]. Chinese Journal of Rock Mechanics & Engineering, 2012, 31(3):507-523. (In Chinese)

    [9] 朱伯芳. 有限單元法原理與應用[M]. 2版. 北京: 中國水利水電出版社, 1998:112-124.

    ZHU Bofang. Principles and applications of finite element method [M]. 2nd ed. Beijing: China Water & Power Press, 1998:112-124.(In Chinese)

    [10]龍述堯, 侯淑娟. 關于無網(wǎng)格方法中點插值形函數(shù)的研究[J]. 湖南大學學報:自然科學版, 2004, 31(3):94-97.

    LONG Shuyao, HOU Shujuan. A study on shape functions in the meshless point interpolation method[J]. Journal of Hunan University: Natural Sciences, 2004, 31(3):94-97. (In Chinese)

    [11]卞康, 肖明. 退化單元位移場插值計算方法研究[J]. 武漢大學學報:工學版, 2008, 41(6):16-20.

    BIAN Kang, XIAO Ming. Research on calculation methods for displacement field interpolation of degenerate elements[J]. Engineering Journal of Wuhan University, 2008, 41(6):16-20.(In Chinese)

    [12]李新強, 楊松青, 汪小剛. 巖體隨機結構面三維網(wǎng)絡的生成和可視化技術[J]. 巖石力學與工程學報, 2007, 26(12):2564-2569.

    LI Xinqiang, YANG Songqing, WANG Xiaogang. Generation and visualization technologies of three-dimensional network of rockmass stochastic structural plane[J]. Chinese Journal of Rock Mechanics & Engineering, 2007, 26(12):2564-2569.(In Chinese)

    [13]SHYU K. Nodal-based discontinuous deformation analysis[D]. Berkeley, USA: Department of Civil Engineering, University of California, Berkeley, 1993:10-20.

    [14]TIAN Q, ZHAO Z, BAO H. Block fracturing analysis using nodal-based discontinuous deformation analysis with the double minimization procedure[J]. International Journal for Numerical & Analytical Methods in Geomechanics, 2013, 81(4):409-415.

    [15]曹代勇, 吳李泉, 占文鋒,等. 斷層應力封閉性的FLAC模擬研究[J]. 煤田地質與勘探, 2005, 33(1):1-3.

    CAO Daiyong, WU Liquan, ZHAN Wenfeng, et al. Application of FLAC in simulation of fault stress seal[J]. Coal Geology & Exploration, 2005, 33(1):1-3. (In Chinese)

    [16]翟曉榮, 吳基文, 彭濤,等. 斷層對煤層開采應力分布影響的數(shù)值模擬分析[J]. 煤礦開采, 2013, 18(5):14-16.

    ZHAI Xiaorong, WU Jiwen, PENG Tao, et al. Numerical simulation for influence of faults on mining stress distribution[J]. Coal Mining Technology, 2013, 18(5):14-16. (In Chinese)

    [17]黃生文, 司鐵漢, 陳文勝,等. 斷層對大跨度隧道圍巖應力影響的有限元分析[J]. 巖石力學與工程學報, 2006, 25(z2):3788-3793.

    HUANG Shengwen, SI Tiehan, CHEN Wensheng, et al. Finite element analyses of influence of fault on large-span tunnel surrounding rock stress[J]. Chinese Journal of Rock Mechanics & Engineering, 2006, 25(z2):3788-3793. (In Chinese)

    国产又色又爽无遮挡免费看| 狂野欧美激情性xxxx| 一级毛片女人18水好多| 欧美三级亚洲精品| 在线永久观看黄色视频| 亚洲精华国产精华精| 日本撒尿小便嘘嘘汇集6| 黄色毛片三级朝国网站| 成人av一区二区三区在线看| 美女高潮到喷水免费观看| 欧美成狂野欧美在线观看| 久久草成人影院| 日日爽夜夜爽网站| 97碰自拍视频| or卡值多少钱| 国产97色在线日韩免费| 欧美乱色亚洲激情| 美女大奶头视频| 日韩欧美免费精品| 特大巨黑吊av在线直播 | 国产精品久久久人人做人人爽| 搡老妇女老女人老熟妇| 99久久综合精品五月天人人| 欧美激情久久久久久爽电影| 男女做爰动态图高潮gif福利片| 亚洲精品中文字幕一二三四区| 国产亚洲欧美98| 国产亚洲精品综合一区在线观看 | 国产精品九九99| 香蕉国产在线看| 午夜久久久久精精品| 精品福利观看| 亚洲第一青青草原| 欧美成狂野欧美在线观看| 精品第一国产精品| 亚洲精品国产区一区二| 久久久久久国产a免费观看| 久久九九热精品免费| 中国美女看黄片| 国产亚洲欧美精品永久| 日韩精品中文字幕看吧| 香蕉av资源在线| 国产精品 国内视频| 亚洲av电影不卡..在线观看| 一级黄色大片毛片| 国产一卡二卡三卡精品| 亚洲五月婷婷丁香| 国产99久久九九免费精品| 成人特级黄色片久久久久久久| 精品国产美女av久久久久小说| 日本 欧美在线| 亚洲国产精品久久男人天堂| 一边摸一边做爽爽视频免费| 波多野结衣高清作品| 桃红色精品国产亚洲av| 欧美最黄视频在线播放免费| 91字幕亚洲| 久久久久久国产a免费观看| 一进一出抽搐动态| 又黄又粗又硬又大视频| 亚洲精品久久成人aⅴ小说| 日本免费a在线| 成人国语在线视频| 老司机在亚洲福利影院| 岛国视频午夜一区免费看| 在线观看免费日韩欧美大片| 一区二区三区国产精品乱码| 成人三级黄色视频| 狠狠狠狠99中文字幕| 欧美一级毛片孕妇| 亚洲色图av天堂| 一a级毛片在线观看| 69av精品久久久久久| 亚洲av电影不卡..在线观看| 在线观看www视频免费| 日韩欧美三级三区| 午夜老司机福利片| 欧美绝顶高潮抽搐喷水| 男人舔女人下体高潮全视频| 桃红色精品国产亚洲av| 色综合婷婷激情| 男女午夜视频在线观看| 亚洲精品中文字幕一二三四区| 人妻久久中文字幕网| 曰老女人黄片| 免费在线观看完整版高清| 午夜成年电影在线免费观看| 无限看片的www在线观看| 热re99久久国产66热| 亚洲午夜理论影院| 国产精品久久久久久亚洲av鲁大| 欧美乱色亚洲激情| 亚洲国产精品sss在线观看| 露出奶头的视频| 一进一出抽搐动态| 色综合婷婷激情| 少妇粗大呻吟视频| 久久久久久久午夜电影| 首页视频小说图片口味搜索| 亚洲色图av天堂| 一a级毛片在线观看| 一夜夜www| 婷婷丁香在线五月| 国产精品99久久99久久久不卡| 午夜福利一区二区在线看| 人成视频在线观看免费观看| 久久热在线av| 国产精品久久久久久亚洲av鲁大| 日韩av在线大香蕉| 黄片播放在线免费| 亚洲精品一区av在线观看| 琪琪午夜伦伦电影理论片6080| 人人妻人人澡人人看| 色综合婷婷激情| 亚洲欧美日韩无卡精品| 99riav亚洲国产免费| 亚洲人成网站高清观看| av天堂在线播放| 免费看美女性在线毛片视频| 一本精品99久久精品77| 午夜视频精品福利| 国产亚洲精品久久久久5区| 精品午夜福利视频在线观看一区| 夜夜看夜夜爽夜夜摸| 亚洲av成人av| 十分钟在线观看高清视频www| 两性午夜刺激爽爽歪歪视频在线观看 | 十分钟在线观看高清视频www| 亚洲专区字幕在线| 精品国产美女av久久久久小说| 91国产中文字幕| 怎么达到女性高潮| 久久国产乱子伦精品免费另类| 国产高清有码在线观看视频 | 亚洲在线自拍视频| 777久久人妻少妇嫩草av网站| 熟妇人妻久久中文字幕3abv| 99久久国产精品久久久| 精品欧美一区二区三区在线| 国产av一区在线观看免费| 亚洲自偷自拍图片 自拍| 亚洲精品中文字幕在线视频| 亚洲精品中文字幕在线视频| 桃红色精品国产亚洲av| 日韩欧美一区视频在线观看| 天堂动漫精品| 黄色视频,在线免费观看| 男人操女人黄网站| 国产人伦9x9x在线观看| 丁香六月欧美| 在线国产一区二区在线| 久久国产亚洲av麻豆专区| 91老司机精品| 99久久99久久久精品蜜桃| 老司机深夜福利视频在线观看| 婷婷精品国产亚洲av| 黄色 视频免费看| 99国产极品粉嫩在线观看| 国内揄拍国产精品人妻在线 | 99在线人妻在线中文字幕| 90打野战视频偷拍视频| 男女做爰动态图高潮gif福利片| 国产精品久久久av美女十八| 久久亚洲精品不卡| 亚洲国产欧洲综合997久久, | 我的亚洲天堂| 成人三级黄色视频| 亚洲成av人片免费观看| 国产亚洲精品av在线| 91成人精品电影| 日韩av在线大香蕉| 在线国产一区二区在线| 国产不卡一卡二| 国产1区2区3区精品| 熟女少妇亚洲综合色aaa.| 精品日产1卡2卡| 我的亚洲天堂| 亚洲自拍偷在线| 中亚洲国语对白在线视频| 夜夜躁狠狠躁天天躁| 男人舔女人的私密视频| 中国美女看黄片| 脱女人内裤的视频| 中文字幕人成人乱码亚洲影| 99在线视频只有这里精品首页| 韩国av一区二区三区四区| 午夜a级毛片| 1024手机看黄色片| 女性被躁到高潮视频| 欧美激情 高清一区二区三区| 日韩有码中文字幕| 欧美精品啪啪一区二区三区| 国产一区二区在线av高清观看| 性欧美人与动物交配| 在线观看66精品国产| 99在线人妻在线中文字幕| 久久精品国产99精品国产亚洲性色| 久久国产精品影院| 亚洲性夜色夜夜综合| 亚洲精品国产一区二区精华液| 欧美日韩黄片免| 叶爱在线成人免费视频播放| 亚洲在线自拍视频| 亚洲片人在线观看| x7x7x7水蜜桃| 亚洲精品在线美女| 最近最新中文字幕大全免费视频| 一级毛片精品| 色综合亚洲欧美另类图片| 三级毛片av免费| 色婷婷久久久亚洲欧美| 亚洲男人天堂网一区| e午夜精品久久久久久久| 国产精品99久久99久久久不卡| 啦啦啦 在线观看视频| 亚洲第一欧美日韩一区二区三区| 日本免费一区二区三区高清不卡| 视频在线观看一区二区三区| 国产精品亚洲美女久久久| 婷婷六月久久综合丁香| 90打野战视频偷拍视频| 国产伦人伦偷精品视频| 一级片免费观看大全| 别揉我奶头~嗯~啊~动态视频| 熟女少妇亚洲综合色aaa.| 中文字幕久久专区| 中文字幕精品免费在线观看视频| avwww免费| 欧美人与性动交α欧美精品济南到| 熟妇人妻久久中文字幕3abv| 视频区欧美日本亚洲| 欧美中文综合在线视频| 色在线成人网| 热99re8久久精品国产| 精品福利观看| 国产亚洲精品久久久久5区| 我的亚洲天堂| 午夜免费观看网址| 亚洲国产中文字幕在线视频| 韩国精品一区二区三区| 久久久久九九精品影院| 婷婷亚洲欧美| 亚洲九九香蕉| 日日爽夜夜爽网站| 国产精品二区激情视频| 久久伊人香网站| 人人妻,人人澡人人爽秒播| 桃红色精品国产亚洲av| 白带黄色成豆腐渣| 精品国产乱码久久久久久男人| 我的亚洲天堂| 久久精品国产清高在天天线| 亚洲avbb在线观看| 亚洲美女黄片视频| 女警被强在线播放| 亚洲一区二区三区色噜噜| 欧美在线一区亚洲| 91在线观看av| 88av欧美| 黄色丝袜av网址大全| 亚洲成人精品中文字幕电影| 无人区码免费观看不卡| 亚洲精品中文字幕在线视频| a在线观看视频网站| 欧美 亚洲 国产 日韩一| 十八禁网站免费在线| 99久久综合精品五月天人人| 亚洲五月色婷婷综合| 亚洲精品色激情综合| 无人区码免费观看不卡| 999精品在线视频| 免费搜索国产男女视频| 99国产精品一区二区三区| 亚洲精品国产精品久久久不卡| 欧美国产日韩亚洲一区| 国内揄拍国产精品人妻在线 | 久久久久久免费高清国产稀缺| 两个人视频免费观看高清| 真人做人爱边吃奶动态| 亚洲熟妇中文字幕五十中出| 中文在线观看免费www的网站 | 成人午夜高清在线视频 | 亚洲av第一区精品v没综合| 黄色a级毛片大全视频| 久久精品aⅴ一区二区三区四区| 亚洲人成伊人成综合网2020| 妹子高潮喷水视频| 18禁裸乳无遮挡免费网站照片 | 亚洲av电影在线进入| 精品国产乱子伦一区二区三区| 午夜a级毛片| 亚洲一区高清亚洲精品| 人人妻,人人澡人人爽秒播| 麻豆国产av国片精品| 久久国产精品影院| 色精品久久人妻99蜜桃| 精品福利观看| 一个人观看的视频www高清免费观看 | av片东京热男人的天堂| 亚洲国产精品合色在线| 波多野结衣av一区二区av| 51午夜福利影视在线观看| 国产v大片淫在线免费观看| 欧美乱色亚洲激情| 国产精品香港三级国产av潘金莲| а√天堂www在线а√下载| 校园春色视频在线观看| 国产一级毛片七仙女欲春2 | 亚洲av片天天在线观看| 男女做爰动态图高潮gif福利片| 在线观看一区二区三区| 国产真人三级小视频在线观看| 69av精品久久久久久| 九色国产91popny在线| 午夜福利在线在线| 国产精品久久久久久人妻精品电影| 波多野结衣巨乳人妻| 欧美三级亚洲精品| 女生性感内裤真人,穿戴方法视频| 色综合欧美亚洲国产小说| 亚洲全国av大片| 在线观看午夜福利视频| 久久国产亚洲av麻豆专区| 成人国产一区最新在线观看| 色综合站精品国产| 国产视频内射| 国产精品99久久99久久久不卡| 听说在线观看完整版免费高清| 久久午夜综合久久蜜桃| 欧美亚洲日本最大视频资源| 国产精品免费一区二区三区在线| 日韩有码中文字幕| 在线观看舔阴道视频| www日本黄色视频网| 一边摸一边抽搐一进一小说| 欧美绝顶高潮抽搐喷水| 一夜夜www| 亚洲一区二区三区不卡视频| 国产精品爽爽va在线观看网站 | 999久久久国产精品视频| 精品欧美国产一区二区三| 成人特级黄色片久久久久久久| av超薄肉色丝袜交足视频| 欧美日韩瑟瑟在线播放| 搡老妇女老女人老熟妇| 国产免费男女视频| 免费女性裸体啪啪无遮挡网站| 色播在线永久视频| 天天添夜夜摸| 搡老熟女国产l中国老女人| 成人亚洲精品一区在线观看| 精品国产国语对白av| 美女高潮到喷水免费观看| 啦啦啦观看免费观看视频高清| 国产精品二区激情视频| 成人欧美大片| 露出奶头的视频| 中文字幕高清在线视频| or卡值多少钱| 无限看片的www在线观看| 香蕉丝袜av| 国产1区2区3区精品| 亚洲精品国产精品久久久不卡| 免费在线观看日本一区| 女生性感内裤真人,穿戴方法视频| 亚洲 欧美一区二区三区| 久久久水蜜桃国产精品网| 精品国内亚洲2022精品成人| 欧美黄色淫秽网站| 国产成人系列免费观看| 99热这里只有精品一区 | 怎么达到女性高潮| 国产精品免费一区二区三区在线| 精品一区二区三区av网在线观看| АⅤ资源中文在线天堂| 国产1区2区3区精品| 色综合站精品国产| 好看av亚洲va欧美ⅴa在| 极品教师在线免费播放| 老司机深夜福利视频在线观看| 搡老妇女老女人老熟妇| 在线十欧美十亚洲十日本专区| 性色av乱码一区二区三区2| 欧美性长视频在线观看| 亚洲黑人精品在线| 黄色成人免费大全| 两性夫妻黄色片| 欧美三级亚洲精品| 日韩有码中文字幕| 欧美日本亚洲视频在线播放| 国产精品亚洲av一区麻豆| 精品久久久久久,| 国产成人欧美在线观看| 久久狼人影院| 制服人妻中文乱码| 在线天堂中文资源库| 可以在线观看毛片的网站| www.999成人在线观看| 窝窝影院91人妻| 国产一区二区在线av高清观看| 精品久久久久久久人妻蜜臀av| 欧美性猛交黑人性爽| 在线观看舔阴道视频| 男女之事视频高清在线观看| 亚洲午夜精品一区,二区,三区| 亚洲五月天丁香| 熟女电影av网| 国产熟女午夜一区二区三区| 神马国产精品三级电影在线观看 | 亚洲自偷自拍图片 自拍| 最近最新中文字幕大全免费视频| 侵犯人妻中文字幕一二三四区| 这个男人来自地球电影免费观看| 亚洲精品在线美女| 久久久精品国产亚洲av高清涩受| 村上凉子中文字幕在线| 色av中文字幕| 久久精品国产清高在天天线| 国产精品电影一区二区三区| 俄罗斯特黄特色一大片| 老熟妇乱子伦视频在线观看| 久久性视频一级片| 久久精品人妻少妇| 久久精品国产99精品国产亚洲性色| 午夜福利欧美成人| 国产高清videossex| 欧美激情 高清一区二区三区| 脱女人内裤的视频| 久久中文字幕一级| 欧美不卡视频在线免费观看 | 天堂影院成人在线观看| 哪里可以看免费的av片| 9191精品国产免费久久| 免费女性裸体啪啪无遮挡网站| 欧美不卡视频在线免费观看 | 黄色a级毛片大全视频| av有码第一页| 国产伦在线观看视频一区| 日本a在线网址| 国产精品99久久99久久久不卡| 黄片大片在线免费观看| 老熟妇仑乱视频hdxx| 国内精品久久久久久久电影| 久久人人精品亚洲av| 首页视频小说图片口味搜索| 国产蜜桃级精品一区二区三区| 久久国产精品人妻蜜桃| 久久久久久久午夜电影| 色婷婷久久久亚洲欧美| 国产单亲对白刺激| 日本 av在线| 国产精品一区二区免费欧美| videosex国产| 日韩精品中文字幕看吧| 青草久久国产| 国产男靠女视频免费网站| 成人亚洲精品av一区二区| 久久天堂一区二区三区四区| 久久 成人 亚洲| 91av网站免费观看| 色婷婷久久久亚洲欧美| 国产精华一区二区三区| 午夜免费成人在线视频| 国产激情久久老熟女| 国产精品免费一区二区三区在线| 精品高清国产在线一区| 久久精品国产亚洲av高清一级| 婷婷亚洲欧美| 男女视频在线观看网站免费 | 19禁男女啪啪无遮挡网站| 最新美女视频免费是黄的| 欧美黄色片欧美黄色片| 亚洲第一电影网av| 亚洲欧美日韩无卡精品| 免费在线观看成人毛片| 一区二区三区高清视频在线| 久久久久九九精品影院| 女同久久另类99精品国产91| 69av精品久久久久久| 中文资源天堂在线| 免费高清视频大片| 午夜激情福利司机影院| 99久久精品国产亚洲精品| 波多野结衣巨乳人妻| 波多野结衣高清作品| 国产精品电影一区二区三区| 美女高潮到喷水免费观看| 可以在线观看毛片的网站| 十八禁网站免费在线| 久久青草综合色| 露出奶头的视频| 久久久久久久久中文| 亚洲精品美女久久久久99蜜臀| 国产成人啪精品午夜网站| 搡老岳熟女国产| 欧美中文日本在线观看视频| 每晚都被弄得嗷嗷叫到高潮| 成年人黄色毛片网站| 国产亚洲av高清不卡| 色精品久久人妻99蜜桃| 免费无遮挡裸体视频| 欧美在线一区亚洲| 成人国产综合亚洲| 亚洲国产精品成人综合色| 成年版毛片免费区| 老鸭窝网址在线观看| 婷婷亚洲欧美| 日韩av在线大香蕉| 曰老女人黄片| 欧美zozozo另类| 757午夜福利合集在线观看| 日韩欧美免费精品| 欧美乱色亚洲激情| 黄频高清免费视频| 国产精品乱码一区二三区的特点| 性色av乱码一区二区三区2| 在线十欧美十亚洲十日本专区| 两个人视频免费观看高清| 日韩精品中文字幕看吧| 久久亚洲真实| 在线观看66精品国产| 99久久国产精品久久久| 国产不卡一卡二| 成人一区二区视频在线观看| 最近最新中文字幕大全免费视频| 老司机靠b影院| 在线av久久热| 久久中文看片网| 香蕉久久夜色| 亚洲成av人片免费观看| 一级a爱片免费观看的视频| 男女那种视频在线观看| 国产黄色小视频在线观看| 精品国产超薄肉色丝袜足j| 亚洲成国产人片在线观看| 国产精品二区激情视频| 少妇熟女aⅴ在线视频| 久久九九热精品免费| 亚洲国产精品sss在线观看| 精品久久久久久久久久免费视频| 免费观看人在逋| 国产一区二区三区视频了| 91麻豆av在线| 一二三四社区在线视频社区8| 又黄又粗又硬又大视频| 欧美色视频一区免费| 中文字幕人成人乱码亚洲影| 免费一级毛片在线播放高清视频| 日本 欧美在线| 自线自在国产av| 成人18禁在线播放| 听说在线观看完整版免费高清| 亚洲熟女毛片儿| 十分钟在线观看高清视频www| 欧美黑人欧美精品刺激| 女性被躁到高潮视频| 日韩有码中文字幕| 在线永久观看黄色视频| 亚洲av成人一区二区三| 18美女黄网站色大片免费观看| 亚洲精品在线观看二区| 日本一本二区三区精品| 午夜亚洲福利在线播放| 少妇粗大呻吟视频| 日本三级黄在线观看| videosex国产| 日韩中文字幕欧美一区二区| 可以免费在线观看a视频的电影网站| 长腿黑丝高跟| 曰老女人黄片| 免费在线观看成人毛片| 国产成年人精品一区二区| 亚洲自偷自拍图片 自拍| 麻豆成人午夜福利视频| 男人舔女人下体高潮全视频| 热re99久久国产66热| 午夜免费激情av| 精品久久久久久,| 一级毛片女人18水好多| 久久久久久久精品吃奶| 大型黄色视频在线免费观看| 国产av一区在线观看免费| 久久人妻福利社区极品人妻图片| 国产一区二区三区视频了| 精品久久久久久成人av| 国产精品自产拍在线观看55亚洲| 亚洲熟女毛片儿| 欧美成人免费av一区二区三区| 国产精品日韩av在线免费观看| 国产精品一区二区精品视频观看| 日本一区二区免费在线视频| 国产91精品成人一区二区三区| 超碰成人久久| 黄色片一级片一级黄色片| 黄片小视频在线播放| 岛国在线观看网站| 桃红色精品国产亚洲av| 亚洲专区中文字幕在线| 午夜精品久久久久久毛片777| 欧美成人免费av一区二区三区| 久久精品夜夜夜夜夜久久蜜豆 | 久久人妻福利社区极品人妻图片| 777久久人妻少妇嫩草av网站| 亚洲成av人片免费观看| 亚洲av成人不卡在线观看播放网| xxxwww97欧美| 色综合亚洲欧美另类图片| 一级毛片女人18水好多| a级毛片a级免费在线| 欧美日本视频| 亚洲九九香蕉| 侵犯人妻中文字幕一二三四区| x7x7x7水蜜桃| 一二三四社区在线视频社区8| 久热爱精品视频在线9|