• <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)

    www日本在线高清视频| 性少妇av在线| 成人毛片a级毛片在线播放| 午夜福利,免费看| 韩国精品一区二区三区| 麻豆精品久久久久久蜜桃| 日韩 亚洲 欧美在线| 欧美日韩国产mv在线观看视频| 一级片免费观看大全| 在线天堂最新版资源| 叶爱在线成人免费视频播放| 日韩熟女老妇一区二区性免费视频| 一级爰片在线观看| 在线观看www视频免费| 人妻少妇偷人精品九色| 午夜福利影视在线免费观看| 色播在线永久视频| 色视频在线一区二区三区| 精品一区在线观看国产| 久久久国产一区二区| 亚洲内射少妇av| 成人国产av品久久久| 成人国语在线视频| 国产国语露脸激情在线看| 日本色播在线视频| 两性夫妻黄色片| a级毛片在线看网站| 丝袜美腿诱惑在线| 精品国产乱码久久久久久男人| 国产精品国产av在线观看| 亚洲三级黄色毛片| 一个人免费看片子| 日日啪夜夜爽| 18禁动态无遮挡网站| 国产一级毛片在线| 精品少妇一区二区三区视频日本电影 | 色播在线永久视频| 中文欧美无线码| 国产黄频视频在线观看| 成人毛片a级毛片在线播放| 日韩人妻精品一区2区三区| 男女啪啪激烈高潮av片| 大话2 男鬼变身卡| 免费女性裸体啪啪无遮挡网站| av有码第一页| videosex国产| 看十八女毛片水多多多| 亚洲精品美女久久av网站| 日本黄色日本黄色录像| 国产欧美亚洲国产| 成年人免费黄色播放视频| 看十八女毛片水多多多| 亚洲成色77777| 黄色一级大片看看| 免费大片黄手机在线观看| 99re6热这里在线精品视频| 成年女人在线观看亚洲视频| a级片在线免费高清观看视频| 美女午夜性视频免费| 美女午夜性视频免费| 久久婷婷青草| 国产xxxxx性猛交| 极品少妇高潮喷水抽搐| 丁香六月天网| 精品亚洲成a人片在线观看| 涩涩av久久男人的天堂| 久久久a久久爽久久v久久| 久久久a久久爽久久v久久| 精品视频人人做人人爽| 考比视频在线观看| 日韩一区二区视频免费看| 一本久久精品| 少妇猛男粗大的猛烈进出视频| 国产男女内射视频| 亚洲欧洲精品一区二区精品久久久 | 亚洲欧美中文字幕日韩二区| 曰老女人黄片| 国产老妇伦熟女老妇高清| 国产精品久久久久久精品电影小说| 亚洲综合精品二区| 欧美人与性动交α欧美软件| 欧美 亚洲 国产 日韩一| 一本大道久久a久久精品| 在线 av 中文字幕| 男人爽女人下面视频在线观看| 啦啦啦中文免费视频观看日本| 免费观看无遮挡的男女| 午夜91福利影院| av电影中文网址| 曰老女人黄片| 另类精品久久| 国产爽快片一区二区三区| 亚洲图色成人| tube8黄色片| 日韩精品免费视频一区二区三区| 韩国精品一区二区三区| 亚洲伊人色综图| 亚洲美女搞黄在线观看| www.自偷自拍.com| 久热这里只有精品99| 欧美变态另类bdsm刘玥| 久久久精品免费免费高清| 18禁国产床啪视频网站| 久久久久国产网址| 欧美日本中文国产一区发布| 亚洲人成77777在线视频| 亚洲精品一二三| 日韩,欧美,国产一区二区三区| 精品国产一区二区三区四区第35| 丁香六月天网| 国产午夜精品一二区理论片| 亚洲男人天堂网一区| 亚洲经典国产精华液单| 亚洲中文av在线| 国产有黄有色有爽视频| 亚洲激情五月婷婷啪啪| 香蕉国产在线看| 日本wwww免费看| 欧美精品国产亚洲| 国产亚洲午夜精品一区二区久久| 亚洲综合色网址| 校园人妻丝袜中文字幕| 天美传媒精品一区二区| 午夜91福利影院| 老司机亚洲免费影院| 麻豆乱淫一区二区| 免费黄网站久久成人精品| 丝瓜视频免费看黄片| 日韩精品免费视频一区二区三区| 男女无遮挡免费网站观看| 日韩在线高清观看一区二区三区| 男人添女人高潮全过程视频| 亚洲精品一二三| 老汉色av国产亚洲站长工具| 国产高清不卡午夜福利| 男女啪啪激烈高潮av片| 免费播放大片免费观看视频在线观看| 成人国产麻豆网| 日本vs欧美在线观看视频| 大码成人一级视频| 久久ye,这里只有精品| 日日摸夜夜添夜夜爱| 欧美人与性动交α欧美软件| 国产精品国产三级国产专区5o| 一区二区三区激情视频| 欧美成人午夜免费资源| 天天影视国产精品| 人人妻人人爽人人添夜夜欢视频| 国产毛片在线视频| 波野结衣二区三区在线| 一级,二级,三级黄色视频| 99香蕉大伊视频| 国产亚洲欧美精品永久| 国产视频首页在线观看| 国产av国产精品国产| a 毛片基地| 亚洲美女搞黄在线观看| 国产又色又爽无遮挡免| 岛国毛片在线播放| 永久网站在线| 国产一区二区三区av在线| 中国三级夫妇交换| 国产人伦9x9x在线观看 | 国产精品麻豆人妻色哟哟久久| 巨乳人妻的诱惑在线观看| 国产成人精品一,二区| 免费看不卡的av| 久久久国产精品麻豆| 黄频高清免费视频| 91aial.com中文字幕在线观看| 69精品国产乱码久久久| 可以免费在线观看a视频的电影网站 | 亚洲欧美色中文字幕在线| 熟女av电影| 久久女婷五月综合色啪小说| 99精国产麻豆久久婷婷| 亚洲精品自拍成人| 亚洲国产av影院在线观看| 久久久精品94久久精品| 亚洲欧美一区二区三区久久| 麻豆乱淫一区二区| 国产精品久久久av美女十八| 国产成人一区二区在线| 免费观看在线日韩| 看非洲黑人一级黄片| 久久ye,这里只有精品| 国产精品久久久久成人av| 韩国高清视频一区二区三区| 丁香六月天网| 亚洲精华国产精华液的使用体验| 亚洲成av片中文字幕在线观看 | 黄频高清免费视频| 在线观看免费日韩欧美大片| 中国三级夫妇交换| 在线观看美女被高潮喷水网站| 午夜老司机福利剧场| 黄网站色视频无遮挡免费观看| 国产精品蜜桃在线观看| 日韩精品免费视频一区二区三区| 久久精品国产亚洲av涩爱| 你懂的网址亚洲精品在线观看| 在线看a的网站| 又大又黄又爽视频免费| 中文字幕最新亚洲高清| 免费日韩欧美在线观看| 国产精品久久久久久久久免| 日韩精品免费视频一区二区三区| 精品酒店卫生间| 亚洲人成电影观看| 有码 亚洲区| av电影中文网址| 精品99又大又爽又粗少妇毛片| 亚洲欧美色中文字幕在线| 80岁老熟妇乱子伦牲交| 国产精品av久久久久免费| 国产精品国产三级专区第一集| 久久久久精品久久久久真实原创| 国产精品 欧美亚洲| 美国免费a级毛片| 黄网站色视频无遮挡免费观看| 自线自在国产av| 国产不卡av网站在线观看| 天天操日日干夜夜撸| 国产成人91sexporn| 精品第一国产精品| 下体分泌物呈黄色| 国产野战对白在线观看| 精品酒店卫生间| 丝袜人妻中文字幕| 熟女少妇亚洲综合色aaa.| 国产成人a∨麻豆精品| 另类精品久久| 久久精品熟女亚洲av麻豆精品| 国产毛片在线视频| 亚洲精品,欧美精品| 大片电影免费在线观看免费| 卡戴珊不雅视频在线播放| 少妇 在线观看| 精品一区二区三卡| 国产熟女欧美一区二区| 日本午夜av视频| 欧美精品av麻豆av| 欧美日韩av久久| 一级爰片在线观看| 国产男女超爽视频在线观看| 哪个播放器可以免费观看大片| 国产成人精品婷婷| 亚洲精品久久午夜乱码| 日韩 亚洲 欧美在线| av在线老鸭窝| 成人免费观看视频高清| 中文欧美无线码| 制服人妻中文乱码| 在线 av 中文字幕| tube8黄色片| 国产男女超爽视频在线观看| 国产欧美日韩综合在线一区二区| 丰满乱子伦码专区| 老司机亚洲免费影院| 乱人伦中国视频| 久久鲁丝午夜福利片| 香蕉国产在线看| 黑人欧美特级aaaaaa片| 亚洲欧美日韩另类电影网站| 国产淫语在线视频| 最近中文字幕2019免费版| 69精品国产乱码久久久| 成人18禁高潮啪啪吃奶动态图| 99久久综合免费| 午夜精品国产一区二区电影| 中文字幕人妻丝袜一区二区 | 黄片小视频在线播放| 男女啪啪激烈高潮av片| 精品午夜福利在线看| 久久韩国三级中文字幕| 激情视频va一区二区三区| 日韩一区二区视频免费看| 成年美女黄网站色视频大全免费| 日本午夜av视频| 色视频在线一区二区三区| 国产无遮挡羞羞视频在线观看| 国产免费又黄又爽又色| 亚洲欧美一区二区三区黑人 | 啦啦啦啦在线视频资源| 亚洲精品一二三| 国产成人精品久久二区二区91 | 亚洲,欧美精品.| 天美传媒精品一区二区| 久久久久精品性色| 久久这里有精品视频免费| 寂寞人妻少妇视频99o| 成人国产麻豆网| 国产av国产精品国产| 纯流量卡能插随身wifi吗| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 午夜日本视频在线| 一本色道久久久久久精品综合| 久久精品久久精品一区二区三区| 精品少妇一区二区三区视频日本电影 | 久久精品国产自在天天线| 久久精品久久精品一区二区三区| 天天躁夜夜躁狠狠躁躁| 波野结衣二区三区在线| 国产麻豆69| 亚洲国产最新在线播放| 精品亚洲成a人片在线观看| 色网站视频免费| 女人精品久久久久毛片| 人妻人人澡人人爽人人| 国产精品蜜桃在线观看| 亚洲精品美女久久av网站| 婷婷色综合大香蕉| 免费观看无遮挡的男女| 999久久久国产精品视频| 亚洲人成77777在线视频| 视频区图区小说| 精品少妇内射三级| 国产一级毛片在线| 看免费成人av毛片| 免费高清在线观看视频在线观看| 精品人妻偷拍中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 91精品国产国语对白视频| 亚洲精品日韩在线中文字幕| 久久鲁丝午夜福利片| 丝瓜视频免费看黄片| 1024香蕉在线观看| 亚洲av欧美aⅴ国产| 伦理电影免费视频| 亚洲成人一二三区av| 精品酒店卫生间| 亚洲国产av影院在线观看| 欧美人与性动交α欧美软件| 精品一区二区三卡| 亚洲精品自拍成人| 丰满迷人的少妇在线观看| 啦啦啦啦在线视频资源| 亚洲久久久国产精品| 搡女人真爽免费视频火全软件| 亚洲精品av麻豆狂野| www.自偷自拍.com| 99精国产麻豆久久婷婷| 99久久中文字幕三级久久日本| 一级毛片我不卡| 亚洲国产看品久久| 国产黄频视频在线观看| 亚洲精品av麻豆狂野| 国产亚洲av片在线观看秒播厂| 黑丝袜美女国产一区| 一本色道久久久久久精品综合| 国产精品熟女久久久久浪| 赤兔流量卡办理| 91aial.com中文字幕在线观看| 久久热在线av| 免费黄网站久久成人精品| 99国产综合亚洲精品| 日本猛色少妇xxxxx猛交久久| 下体分泌物呈黄色| 伦精品一区二区三区| 精品国产一区二区久久| 在现免费观看毛片| 尾随美女入室| 中文字幕人妻熟女乱码| 欧美人与性动交α欧美软件| 精品人妻在线不人妻| 午夜av观看不卡| 大片免费播放器 马上看| 美女午夜性视频免费| 亚洲精品久久久久久婷婷小说| 精品一区在线观看国产| 999精品在线视频| 99精国产麻豆久久婷婷| 大片电影免费在线观看免费| 麻豆精品久久久久久蜜桃| 美女高潮到喷水免费观看| 少妇的丰满在线观看| 欧美日韩精品网址| 在线精品无人区一区二区三| 久久久久精品久久久久真实原创| 久久 成人 亚洲| 午夜福利视频在线观看免费| 精品国产露脸久久av麻豆| √禁漫天堂资源中文www| 婷婷色av中文字幕| 国产高清不卡午夜福利| 2018国产大陆天天弄谢| 婷婷成人精品国产| av.在线天堂| 亚洲欧洲日产国产| 丰满饥渴人妻一区二区三| 五月开心婷婷网| 大码成人一级视频| 亚洲国产欧美日韩在线播放| 中文字幕人妻丝袜制服| www.精华液| 黄网站色视频无遮挡免费观看| 免费在线观看完整版高清| 国产精品无大码| av卡一久久| 亚洲精品久久成人aⅴ小说| 午夜av观看不卡| av电影中文网址| 少妇的逼水好多| 少妇人妻 视频| a 毛片基地| 99re6热这里在线精品视频| 七月丁香在线播放| 人妻少妇偷人精品九色| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 男人操女人黄网站| 久久综合国产亚洲精品| 99久久人妻综合| 永久免费av网站大全| 一区二区三区四区激情视频| 亚洲欧美色中文字幕在线| 日产精品乱码卡一卡2卡三| 亚洲成人手机| 久久久精品免费免费高清| 亚洲天堂av无毛| 女性被躁到高潮视频| 欧美bdsm另类| 春色校园在线视频观看| av视频免费观看在线观看| 亚洲欧美成人综合另类久久久| 毛片一级片免费看久久久久| 久久国产精品男人的天堂亚洲| 一级片免费观看大全| 国产av国产精品国产| 巨乳人妻的诱惑在线观看| 91aial.com中文字幕在线观看| 咕卡用的链子| 人人妻人人爽人人添夜夜欢视频| 久久这里只有精品19| 丝袜脚勾引网站| 成人亚洲欧美一区二区av| 中文字幕最新亚洲高清| www.自偷自拍.com| 这个男人来自地球电影免费观看 | 欧美中文综合在线视频| 狠狠精品人妻久久久久久综合| 成人亚洲欧美一区二区av| √禁漫天堂资源中文www| 久久久久久免费高清国产稀缺| av.在线天堂| 久久97久久精品| 精品一区二区三卡| 日韩不卡一区二区三区视频在线| 午夜福利影视在线免费观看| 尾随美女入室| 亚洲av福利一区| 女人被躁到高潮嗷嗷叫费观| 最近最新中文字幕大全免费视频 | 熟女av电影| 十分钟在线观看高清视频www| 国产日韩欧美视频二区| 一级毛片黄色毛片免费观看视频| 久久ye,这里只有精品| 建设人人有责人人尽责人人享有的| 国产乱来视频区| 亚洲欧美精品自产自拍| 中文字幕av电影在线播放| 考比视频在线观看| 少妇人妻久久综合中文| 久久久精品区二区三区| 国产成人精品久久久久久| 多毛熟女@视频| 亚洲精品一区蜜桃| 成人黄色视频免费在线看| 下体分泌物呈黄色| 人体艺术视频欧美日本| 久久久精品94久久精品| 国产乱人偷精品视频| 国产精品一国产av| 水蜜桃什么品种好| 99国产综合亚洲精品| 亚洲 欧美一区二区三区| 亚洲精品久久成人aⅴ小说| 美女国产视频在线观看| 亚洲国产欧美日韩在线播放| 1024视频免费在线观看| 最近手机中文字幕大全| 亚洲国产日韩一区二区| 中文字幕最新亚洲高清| av在线老鸭窝| 少妇人妻 视频| 久热久热在线精品观看| 日本免费在线观看一区| 看免费av毛片| kizo精华| 最近最新中文字幕免费大全7| 国产日韩欧美视频二区| av网站免费在线观看视频| 最近手机中文字幕大全| 岛国毛片在线播放| 国产亚洲最大av| 午夜影院在线不卡| 色婷婷av一区二区三区视频| 丝袜在线中文字幕| 最新的欧美精品一区二区| 亚洲精品日本国产第一区| 叶爱在线成人免费视频播放| 少妇猛男粗大的猛烈进出视频| av视频免费观看在线观看| 大片电影免费在线观看免费| 久久久a久久爽久久v久久| 亚洲综合色惰| 中文字幕制服av| 国产一区二区三区av在线| 少妇 在线观看| 欧美激情高清一区二区三区 | 美女xxoo啪啪120秒动态图| 亚洲精品久久成人aⅴ小说| 各种免费的搞黄视频| 99久久人妻综合| 国产国语露脸激情在线看| 亚洲av综合色区一区| 视频在线观看一区二区三区| 国产精品嫩草影院av在线观看| 99国产精品免费福利视频| 黑人猛操日本美女一级片| 亚洲国产看品久久| av免费在线看不卡| 日日啪夜夜爽| 亚洲av成人精品一二三区| 飞空精品影院首页| 亚洲av中文av极速乱| 777久久人妻少妇嫩草av网站| 99精国产麻豆久久婷婷| 啦啦啦在线免费观看视频4| 久久久久精品人妻al黑| 久久热在线av| 宅男免费午夜| 激情五月婷婷亚洲| 久久这里有精品视频免费| 久久久久国产一级毛片高清牌| 免费黄网站久久成人精品| 制服诱惑二区| 丰满饥渴人妻一区二区三| 18在线观看网站| 秋霞在线观看毛片| 亚洲av国产av综合av卡| 水蜜桃什么品种好| 母亲3免费完整高清在线观看 | av在线app专区| 桃花免费在线播放| 涩涩av久久男人的天堂| 精品一区二区三区四区五区乱码 | 国产免费现黄频在线看| 久久精品久久精品一区二区三区| 久久毛片免费看一区二区三区| av网站在线播放免费| 日韩视频在线欧美| 高清av免费在线| 天天影视国产精品| 一级片'在线观看视频| 青春草视频在线免费观看| 国产男女内射视频| 日韩精品免费视频一区二区三区| 女人精品久久久久毛片| 午夜福利视频在线观看免费| 中国三级夫妇交换| 在线观看一区二区三区激情| 中文字幕色久视频| 男人添女人高潮全过程视频| 午夜久久久在线观看| 18禁国产床啪视频网站| 久久影院123| 亚洲久久久国产精品| 综合色丁香网| 免费在线观看完整版高清| 亚洲精品日韩在线中文字幕| 秋霞在线观看毛片| 成人亚洲欧美一区二区av| 9191精品国产免费久久| www.av在线官网国产| 一本久久精品| 在线看a的网站| 一级黄片播放器| 国产一区二区三区av在线| 亚洲国产欧美日韩在线播放| 国产欧美日韩综合在线一区二区| 少妇熟女欧美另类| 久久99蜜桃精品久久| 久久久精品区二区三区| 2022亚洲国产成人精品| 激情视频va一区二区三区| 制服诱惑二区| 狠狠精品人妻久久久久久综合| 五月伊人婷婷丁香| 精品亚洲乱码少妇综合久久| 日韩三级伦理在线观看| 午夜福利视频在线观看免费| 在线观看三级黄色| 免费看不卡的av| 免费观看性生交大片5| 久久久久精品人妻al黑| 高清不卡的av网站| 亚洲在久久综合| 黄色毛片三级朝国网站| 老司机影院毛片| 亚洲国产色片| 亚洲av国产av综合av卡| 国产精品久久久久久精品古装| 亚洲美女黄色视频免费看| 欧美 日韩 精品 国产| 精品一区二区三区四区五区乱码 | 精品国产露脸久久av麻豆| 女性被躁到高潮视频| 日韩一本色道免费dvd| 啦啦啦啦在线视频资源| 国产国语露脸激情在线看| 波多野结衣av一区二区av| 国产男女内射视频| 国产毛片在线视频| 男男h啪啪无遮挡|