多種流形覆蓋方式耦合的數(shù)值流形法及其應(yīng)用

武鑫 郭璇 康歌 甯尤軍

摘要：為了更好地利用數(shù)值流形法（NMM）模擬連續(xù)和非連續(xù)問(wèn)題，根據(jù)NMM數(shù)學(xué)覆蓋選擇的靈活性，將有限元網(wǎng)格作為一種數(shù)學(xué)覆蓋方式生成流形覆蓋系統(tǒng)以避免細(xì)小流形單元的出現(xiàn)，并提出了將單個(gè)材料體描述為一個(gè)流形單元的獨(dú)立覆蓋方法。運(yùn)用多種覆蓋方式相耦合的方法模擬了混凝土的細(xì)觀拉壓模型，以及一個(gè)復(fù)雜結(jié)構(gòu)巖質(zhì)邊坡模型，模擬結(jié)果較好地再現(xiàn)了相應(yīng)的力學(xué)變形破壞過(guò)程。多種覆蓋方式的耦合使得構(gòu)建的NMM模擬模型更加合理；獨(dú)立覆蓋方式的使用降低了構(gòu)建離散體系NMM模型的復(fù)雜程度，并減少了模型中的單元數(shù)量以及不連續(xù)面上的接觸數(shù)量，提高了模擬計(jì)算效率。

關(guān)鍵詞：數(shù)值流形法（NMM）；獨(dú)立覆蓋；有限覆蓋；耦合使用；離散體系

中圖分類號(hào)：O347.7 文獻(xiàn)標(biāo)志碼：A文章編號(hào)：1000-582X（2018）05-060-08

Abstract： To better use the numerical manifold method （NMM） to simulate complex continuous and discontinuous deformation problems， for the flexibility of the choice of mathematical covers in the NMM， we adopt the finite element mesh to generate NMM mathematical covers so as to avoid the appearance of too-small manifold elements， and put forward an independent cover method， in which a single material domain is described as one manifold element. Different cover types are coupled in the simulations of the tensile and compressive concrete meso-scale models， as well as a rock slope model with complex structures. The corresponding mechanical deformation and failure processes are well reproduced by NMM simulations. The coupling use of different cover types makes the generated NMM simulation models more reasonable. The independent cover method decreases the complexity to construct NMM simulation models for discrete systems， and reduces the numbers of manifold elements and the contact numbers along discontinuity surfaces， which helps to improve the NMM computational efficiency.

Keywords： numerical manifold method （NMM）；independent cover；finite cover；coupling use；discrete system

混凝土、巖石是典型的非均勻、非連續(xù)工程及地質(zhì)材料，其變形、破壞是一個(gè)復(fù)雜的連續(xù)、非連續(xù)固體介質(zhì)力學(xué)行為問(wèn)題。數(shù)值模擬是研究這類復(fù)雜力學(xué)問(wèn)題的重要手段?；谶B續(xù)理論的有限元法（FEM），在模擬材料不連續(xù)變形和運(yùn)動(dòng)力學(xué)行為時(shí)能力存在不足。在基于非連續(xù)理論的方法中，離散元法（DEM）運(yùn)用顯式有限差分法模擬多體不連續(xù)系統(tǒng)的大位移和旋轉(zhuǎn)過(guò)程[1]；非連續(xù)變形分析（DDA）方法則基于最小勢(shì)能原理，較好地解決了多材料體的大變形和大位移問(wèn)題。對(duì)于巖石力學(xué)行為的模擬，DEM和DDA方法適用于對(duì)離散巖體結(jié)構(gòu)的穩(wěn)定狀態(tài)和變形模式及過(guò)程的模擬，但對(duì)于連續(xù)體內(nèi)裂紋的產(chǎn)生和擴(kuò)展模擬則存在一定不足[2]。顯然，2種數(shù)值方法對(duì)于混凝土這種顆?；w材料的細(xì)觀變形破壞模擬也同樣存在難度。連續(xù)非連續(xù)數(shù)值方法的耦合，如FEM-DEM[3]，以及連續(xù)非連續(xù)理論在數(shù)值方法中的統(tǒng)一，如數(shù)值流形法（NMM）[4]，則是更好的解決途徑。

數(shù)值流形法（NMM， numerical manifold method）通過(guò)2套獨(dú)立的覆蓋系統(tǒng)很好地將連續(xù)與非連續(xù)問(wèn)題統(tǒng)一在一個(gè)理論框架內(nèi)進(jìn)行描述，將FEM、DDA和解析法統(tǒng)一起來(lái)。覆蓋系統(tǒng)是NMM的重要組成部分，對(duì)此已有了較多的研究。蔡永昌等[5]論述了二維彈性問(wèn)題里使用矩形覆蓋流形方法的理論和實(shí)施技術(shù)，并給出了其覆蓋系統(tǒng)的全自動(dòng)生成算法；陳剛等[6]將有限元網(wǎng)格作為生成流形覆蓋系統(tǒng)的數(shù)學(xué)網(wǎng)格，利用有向圖的關(guān)聯(lián)矩陣結(jié)合塊體搜索算法，形成最終的流形法覆蓋系統(tǒng)；李海楓等[7]研究了三維流形單元的生成方法；蔡永昌等[8]提出了一種完全獨(dú)立覆蓋的高階NMM，在獨(dú)立覆蓋之間假設(shè)幾何厚度可以忽略的連接彈簧，以解決高階流形法的線性相關(guān)性問(wèn)題和物理覆蓋系統(tǒng)的生成復(fù)雜性問(wèn)題；蘇海東等[9]研究了矩形獨(dú)立覆蓋及其覆蓋加密，在獨(dú)立覆蓋中使用解析級(jí)數(shù)，并提出了任意形狀的獨(dú)立覆蓋形式。NMM在類混凝土材料的細(xì)觀變形破壞，以及巖質(zhì)邊坡失穩(wěn)等巖石力學(xué)與工程問(wèn)題的研究中都已得到了一定的應(yīng)用。黃濤[10]和趙妍[11]分別采用NMM模擬研究了混凝土以及高聚物粘結(jié)炸藥這兩類顆?；w材料的細(xì)觀變形破壞問(wèn)題；張國(guó)新等[2]驗(yàn)證了采用NMM可以正確計(jì)算邊坡傾倒安全系數(shù)和模擬邊坡傾倒破壞過(guò)程；黃濤等[12]應(yīng)用NMM對(duì)巖石雙孔爆破過(guò)程進(jìn)行了模擬；Ning和An等[13-14]采用基于莫爾庫(kù)侖開(kāi)裂準(zhǔn)則的NMM模擬了巖質(zhì)邊坡中的裂紋擴(kuò)展和邊坡破壞全過(guò)程。

根據(jù)NMM中數(shù)學(xué)覆蓋和物理覆蓋的定義，提出了將單個(gè)材料體描述為單個(gè)流形單元的獨(dú)立流形覆蓋方法，并與傳統(tǒng)的規(guī)則網(wǎng)格有限覆蓋以及由有限元網(wǎng)格生成的有限覆蓋相結(jié)合，以便對(duì)包含離散體系的復(fù)雜材料結(jié)構(gòu)的力學(xué)行為進(jìn)行模擬。通過(guò)混凝土細(xì)觀變形破壞行為模擬和復(fù)雜結(jié)構(gòu)巖質(zhì)邊坡失穩(wěn)破壞模擬算例，驗(yàn)證和討論了這一數(shù)值模擬途徑的可行性及優(yōu)勢(shì)。

1 數(shù)值流形法的有限覆蓋與獨(dú)立覆蓋

1.1 數(shù)值流形法的基本思想

采用2套獨(dú)立的覆蓋系統(tǒng)：數(shù)學(xué)覆蓋系統(tǒng)和物理覆蓋系統(tǒng)[3]。數(shù)學(xué)覆蓋（MC， mathematical cover）系統(tǒng)由許多部分重疊的覆蓋片組成，它覆蓋住整個(gè)物理區(qū)域。在二維NMM中，一般多采用規(guī)則三角形網(wǎng)格生成數(shù)學(xué)覆蓋系統(tǒng)。如圖1（a）所示，當(dāng)采用三角形網(wǎng)格生成數(shù)學(xué)覆蓋系統(tǒng)時(shí)，對(duì)于網(wǎng)格中的任一節(jié)點(diǎn)，含有這一節(jié)點(diǎn)的所有三角形組成1個(gè)數(shù)學(xué)覆蓋。在圖1（a）中，1-6號(hào)三角形組成的正六邊形就是該節(jié)點(diǎn)處的數(shù)學(xué)覆蓋。

物理覆蓋（PC， physical cover）是數(shù)學(xué)覆蓋和物理區(qū)域的交集，當(dāng)將數(shù)學(xué)覆蓋置于物理區(qū)域上時(shí)，數(shù)學(xué)覆蓋與物理區(qū)域重疊的部分就形成物理覆蓋系統(tǒng)。如果材料邊界或內(nèi)部不連續(xù)面將1個(gè)數(shù)學(xué)覆蓋切割成多個(gè)完全獨(dú)立的區(qū)域，則位于材料域內(nèi)的每一獨(dú)立區(qū)域就是1個(gè)物理覆蓋。在圖1（a）中，以節(jié)點(diǎn)a、b為中心的六邊形數(shù)學(xué)覆蓋Ma、Mb完全位于物理區(qū)域內(nèi)，并且均沒(méi)有被不連續(xù)面切割成完全獨(dú)立的區(qū)域，從而使得Ma、Mb分別成為1個(gè)物理覆蓋，記為P1a、P1b；以節(jié)點(diǎn)c為中心的數(shù)學(xué)覆蓋Mc只有一部分與物理區(qū)域相重疊，只有位于物理區(qū)域內(nèi)的這一部分才形成1個(gè)物理覆蓋，記為P1c，同理，以節(jié)點(diǎn)e為中心的數(shù)學(xué)覆蓋Me處形成1個(gè)物理覆蓋，記為P1e；而對(duì)于以節(jié)點(diǎn)d為中心的數(shù)學(xué)覆蓋Md而言，其位于物理區(qū)域內(nèi)的部分被物理區(qū)域中的不連續(xù)面切割成2個(gè)完全獨(dú)立的區(qū)域，所以，在節(jié)點(diǎn)d處形成2個(gè)物理覆蓋，記為P1d、P2d。

對(duì)于由以上方法產(chǎn)生的物理覆蓋，他們也是相互部分重疊的，這些相互重疊的物理覆蓋所形成的交集則構(gòu)成流形單元（ME， manifold element）。如圖1所示，圖1（a）中的物理覆蓋P1a、P1b和P1c相交形成流形單元e1，物理覆蓋P1b、P1d、P1e相交形成流形單元e2，物理覆蓋P1b、P2d、P1e相交形成流形單元e3。顯然，每個(gè)流形單元是3個(gè)物理覆蓋的交集。

1.2 有限覆蓋與獨(dú)立覆蓋

如圖1所示的由規(guī)則三角形網(wǎng)格生成數(shù)學(xué)覆蓋，是NMM廣泛采用的一種覆蓋形式。由于不需要網(wǎng)格邊與物理邊界相重合，特別是對(duì)于復(fù)雜形態(tài)的物理區(qū)域，前處理變得易于實(shí)現(xiàn)，但同時(shí)也會(huì)在物理區(qū)域的邊界產(chǎn)生一些面積極小或細(xì)長(zhǎng)的流形單元，如圖2（a）所示，從而對(duì)計(jì)算性能造成一定影響。生成NMM覆蓋系統(tǒng)的另外一種方式是采用傳統(tǒng)有限元網(wǎng)格作為數(shù)學(xué)覆蓋，如圖2（b）所示，當(dāng)采用三角形有限元網(wǎng)格時(shí)，對(duì)于其中任一節(jié)點(diǎn)，含有這一節(jié)點(diǎn)的所有三角形單元形成1個(gè)數(shù)學(xué)覆蓋，進(jìn)而生成相應(yīng)的物理覆蓋。對(duì)于任意一個(gè)有限單元，其3個(gè)頂點(diǎn)處的物理覆蓋相互重疊形成的流形單元就是該有限單元本身。這時(shí)NMM方法實(shí)質(zhì)上退化為了有限單元法，只是在NMM理論框架下，材料開(kāi)裂破壞的描述較有限元更加易于實(shí)現(xiàn)。無(wú)論是采用規(guī)則三角形網(wǎng)格還是采用有限元非規(guī)則網(wǎng)格生成的流形覆蓋，都可以稱為NMM的有限覆蓋。采用有限覆蓋可以很方便地描述連續(xù)材料體的變形和開(kāi)裂破壞行為。

對(duì)于離散體系而言，主要關(guān)心離散體的位移和相互之間的作用，而離散體自身的變形不重要時(shí)，采用如上所述的有限覆蓋則不具必要性和計(jì)算經(jīng)濟(jì)性。這里提出獨(dú)立覆蓋的概念，即每一個(gè)離散體物理區(qū)域采用獨(dú)立的覆蓋進(jìn)行描述，每個(gè)塊體形成一個(gè)流形單元，并擁有屬于自身的獨(dú)立的位移函數(shù)，實(shí)現(xiàn)離散體間位移不連續(xù)的描述。這時(shí)，NMM實(shí)質(zhì)上退化為了DDA方法，只是采用了不同的具體未知量和位移函數(shù)。值得注意的是，這里提出的獨(dú)立覆蓋的概念和文獻(xiàn)[8-9]中獨(dú)立覆蓋的概念是不一樣的。

當(dāng)采用三角形網(wǎng)格時(shí)，獨(dú)立覆蓋的構(gòu)造方法是生成1個(gè)包絡(luò)單個(gè)離散體的三角形，并以該三角形的3個(gè)頂點(diǎn)為數(shù)學(xué)覆蓋的中心，產(chǎn)生3個(gè)數(shù)學(xué)覆蓋，進(jìn)而生成相應(yīng)3個(gè)物理覆蓋。如圖3（a）所示，1個(gè)四邊形離散體物理區(qū)域被覆蓋其上的Δijk完全包含。根據(jù)數(shù)學(xué)覆蓋的概念，節(jié)點(diǎn)i、j、k處的數(shù)學(xué)覆蓋分別為以節(jié)點(diǎn)i、j、k為中心的3個(gè)六邊形， 記為Mi、Mj和Mk。它們與物理區(qū)域的交集分別構(gòu)成3個(gè)四邊形物理覆蓋Pi、Pj、Pk（區(qū)域均為物理區(qū)域本身），這3個(gè)四邊形物理覆蓋相互重疊形成1個(gè)四邊形流形單元，即該物理區(qū)域本身。這樣以來(lái)，單個(gè)離散體物理區(qū)域自身形成1個(gè)流形單元。對(duì)于如圖3（b）所示的離散體系統(tǒng)，其中每一個(gè)四邊形離散體均擁有屬于自身的獨(dú)立覆蓋，各自形成獨(dú)立的流形單元。

以上獨(dú)立覆蓋生成方法略去了數(shù)學(xué)覆蓋被物理區(qū)域再分割的過(guò)程，簡(jiǎn)化了物理覆蓋與流形單元的生成過(guò)程。在保證相互對(duì)應(yīng)的前提下，覆蓋系統(tǒng)和流形單元可以分開(kāi)建立。具體可以先建立離散體系統(tǒng)模型，每個(gè)離散體形成1個(gè)流形單元，進(jìn)而根據(jù)離散體的幾何和坐標(biāo)信息，對(duì)每個(gè)離散體引入1個(gè)三角形并形成相應(yīng)的數(shù)學(xué)與物理覆蓋。如圖3（b）所示，由此生成的覆蓋系統(tǒng)，保證了屬于任何2個(gè)離散體的物理覆蓋間的獨(dú)立性。

對(duì)于具體問(wèn)題的模擬，在覆蓋系統(tǒng)的生成過(guò)程中，可以根據(jù)問(wèn)題的幾何及物理特征，將由規(guī)則三角形或有限元網(wǎng)格生成的有限覆蓋，與獨(dú)立覆蓋進(jìn)行耦合使用，得到更加合理的NMM模擬模型。在文中，有限覆蓋和獨(dú)立覆蓋采用相同的覆蓋位移函數(shù)，自然地將各種覆蓋方式耦合到了統(tǒng)一的NMM數(shù)學(xué)描述框架下，具有簡(jiǎn)潔性和實(shí)用性。

2 不同流形覆蓋方式的耦合應(yīng)用

2.1 混凝土細(xì)觀模型的單軸拉壓模擬

混凝土在細(xì)觀尺度上可看成是由骨料、水泥砂漿及其交界面組成的顆?；w復(fù)合材料。其中，骨料與水泥砂漿的交界面是混凝土材料內(nèi)部最薄弱的相，在外荷載作用下，微裂紋往往率先在交界面上開(kāi)始萌生，并逐漸演化、發(fā)展和貫穿，最終形成宏觀裂紋導(dǎo)致混凝土材料的破壞。混凝土拉伸破壞后的斷面形態(tài)表明，絕大部分破壞面是粗骨料與水泥砂漿的界面粘結(jié)破壞（拉脫），而小部分是骨料間的水泥砂漿被拉斷。即使是在壓縮載荷作用下的破壞，一般也是以水泥砂漿的開(kāi)裂和界面開(kāi)裂為主，較少發(fā)生骨料破壞的情況。因此，對(duì)于混凝土這種顆粒基體復(fù)合材料的細(xì)觀模擬，骨料的變形和破壞一般可以認(rèn)為并不重要。

在混凝土細(xì)觀模型的NMM模擬中，可以采用獨(dú)立覆蓋描述骨料，采用有限覆蓋描述砂漿，界面則通過(guò)粘接強(qiáng)度加以考慮。如圖4（a）所示，其中的骨料、加載板以及測(cè)量板采用了獨(dú)立覆蓋，砂漿則采用了有限元網(wǎng)格來(lái)生成覆蓋系統(tǒng)，以避免采用規(guī)則三角形網(wǎng)格生成覆蓋時(shí)在砂漿骨料界面處產(chǎn)生細(xì)長(zhǎng)或極小尺寸流形單元。測(cè)量板位于加載板和細(xì)觀模型之間，因采用獨(dú)立覆蓋而具有常應(yīng)力，它在豎直方向上的應(yīng)力反映了混凝土細(xì)觀模型在豎向方向上的平均應(yīng)力。為加以對(duì)比，在圖4（b）的模型中，對(duì)骨料也采用有限元網(wǎng)格進(jìn)行了覆蓋。

在方件上下兩端均通過(guò)加載板進(jìn)行準(zhǔn)靜態(tài)線性拉伸或壓縮位移加載，加載速率為2 mm/s。假定骨料與水泥砂漿為線彈性材料，各類材料參數(shù)如表1所示。NMM模擬中的最大步位移比、時(shí)間步長(zhǎng)和接觸彈簧剛度分別取為5×10-4，1×10-5 s和600 GPa。

圖5（a）與圖5（b）所示為單軸拉伸的細(xì)觀模擬結(jié)果。骨料采用2種不同覆蓋方式，模型中主貫通裂紋均產(chǎn)生于模型中下部，并伴有多處分散裂紋，其中，砂漿骨料界面開(kāi)裂占整個(gè)開(kāi)裂面的較大部分。2種覆蓋條件下模型的細(xì)觀破壞形態(tài)極為相似。圖6（a）給出了由測(cè)量板得到的模型在豎直方向上的平均應(yīng)力變化曲線。曲線在初期表現(xiàn)為線性，隨著加載位移的增大，界面的逐步脫粘和砂漿的逐步開(kāi)裂導(dǎo)致軟化現(xiàn)象的出現(xiàn)。在主裂紋逐漸貫通的過(guò)程中，曲線開(kāi)始出現(xiàn)明顯下降。由于材料內(nèi)部結(jié)構(gòu)的不均勻性，曲線在下降過(guò)程中出現(xiàn)了一定的起伏現(xiàn)象。可以看出，骨料分別采用有限覆蓋和獨(dú)立覆蓋時(shí)，模型豎直方向的應(yīng)力位移曲線在前期極為近似，只是在后期的起伏階段才產(chǎn)生了較明顯差異。

圖5（c）與圖5（d）所示為單軸壓縮的模擬結(jié)果。在壓縮載荷作用下，細(xì)觀結(jié)構(gòu)的非均勻性導(dǎo)致大量分散裂紋的產(chǎn)生，沒(méi)有形成明顯優(yōu)勢(shì)裂紋帶，但單條分散裂紋基本都為傾斜方向的剪切裂紋。可以看出，在骨料采用2種不同覆蓋的條件下，模型中細(xì)觀破壞形態(tài)仍極為相似。圖6（b）所示的模型豎直方向平均應(yīng)力變化曲線表明，初期模型整體響應(yīng)呈線性，隨著分散裂紋的逐步產(chǎn)生，軟化現(xiàn)象逐步顯現(xiàn)。對(duì)比骨料采用2種不同覆蓋方式得到的應(yīng)力曲線，可以發(fā)現(xiàn)它們?nèi)苑浅＝咏?/p>

綜上所述，對(duì)于混凝土這樣的顆?；w材料，當(dāng)骨料顆粒分別采用有限覆蓋和獨(dú)立覆蓋時(shí)，得到了基本一致的模擬結(jié)果。但當(dāng)骨料采用獨(dú)立覆蓋時(shí)，覆蓋系統(tǒng)的生成更加簡(jiǎn)單，并且由于界面處接觸數(shù)量的減少，以及整個(gè)模型中單元數(shù)目的減少，模擬計(jì)算效率也得到了一定的提高。

2.2 復(fù)雜巖質(zhì)邊坡的失穩(wěn)破壞模擬

巖體是由連續(xù)巖石介質(zhì)和不連續(xù)面構(gòu)成的天然地質(zhì)材料，巖質(zhì)邊坡復(fù)雜的組成結(jié)構(gòu)對(duì)其破壞形式有著直接重要影響。根據(jù)Goodman[15-16]關(guān)于巖質(zhì)結(jié)構(gòu)邊坡失穩(wěn)的研究，建立如圖7（a）所示的由連續(xù)巖石介質(zhì)與離散巖塊構(gòu)成的邊坡模型。在由連續(xù)巖石介質(zhì)和少量不連續(xù)面組成的節(jié)理欠發(fā)育巖體中，連續(xù)巖石介質(zhì)的開(kāi)裂破壞是非常重要的破壞形式。而在包含大量不連續(xù)面的節(jié)理較充分的發(fā)育巖體中，巖體可看做是由獨(dú)立巖塊構(gòu)成的離散巖塊系統(tǒng)，巖體的變形破壞主要沿著已有的不連續(xù)面進(jìn)行。下面通過(guò)NMM方法模擬該邊坡在重力作用下的失穩(wěn)破壞過(guò)程?；谇拔乃峒暗挠邢薷采w和獨(dú)立覆蓋2種不同覆蓋方式，邊坡中節(jié)理欠發(fā)育巖體部分采用規(guī)則三角形網(wǎng)格進(jìn)行有限覆蓋以描述其變形破壞，邊坡中節(jié)理充分發(fā)育巖體部分采用獨(dú)立覆蓋，以便于描述離散巖塊系統(tǒng)的變形破壞及運(yùn)動(dòng)。同時(shí)，用作邊界約束條件的模型邊框也采用獨(dú)立覆蓋。NMM網(wǎng)格模型如圖7（b）所示。

模擬中，節(jié)理充分發(fā)育部分巖體（記為巖層）的材料參數(shù)為：密度2 500 kg/m3，彈性模量30 GPa，泊松比0.25。節(jié)理欠發(fā)育部分巖體（記為巖塊）的材料參數(shù)為：密度2 700 kg/m3，彈性模量30 GPa，泊松比0.25。巖層的強(qiáng)度參數(shù)取為：內(nèi)摩擦角27°，內(nèi)聚力0.5 MPa，抗拉強(qiáng)度0.5 MPa，巖塊部分不考慮開(kāi)裂破壞。模型中不連續(xù)面的強(qiáng)度參數(shù)均取為：摩擦角10°，內(nèi)聚力0.05 MPa，抗拉強(qiáng)度0.05 MPa。這里取相對(duì)一般實(shí)際更小的強(qiáng)度參數(shù)，是為了使得邊坡在重力作用下即發(fā)生破壞失穩(wěn)。

圖7（c）和圖7（d）給出了邊坡破壞過(guò)程的NMM模擬圖像。在自重和上部巖體的壓載作用下，節(jié)理發(fā)育良好的邊坡下部巖體沿既有不連續(xù)面發(fā)生破壞，巖塊系統(tǒng)發(fā)生松動(dòng)。下部離散巖塊系統(tǒng)承載能力的降低，導(dǎo)致上部巖層的不穩(wěn)定并沿豎向貫通節(jié)理面發(fā)生下錯(cuò)，巖層中既有的非貫通節(jié)理面進(jìn)而發(fā)生擴(kuò)展，將巖層逐步折斷。隨后，下部離散巖塊被向右擠出并向前方傾倒，上部破壞巖層下端與下部離散巖塊一起向右側(cè)運(yùn)動(dòng)，整個(gè)上部巖層則發(fā)生向后傾倒并下滑。這種破壞稱為滑動(dòng)傾倒的“二次傾倒模式”[15]，當(dāng)上部巖層滑動(dòng)趨勢(shì)被下端巖塊阻止時(shí)發(fā)生。以上模擬結(jié)果表明，采用有限覆蓋和獨(dú)立覆蓋耦合的方式能很好地模擬復(fù)雜巖質(zhì)邊坡的失穩(wěn)破壞問(wèn)題。

3 結(jié) 論

基于NMM有限覆蓋理論，采用有限元網(wǎng)格作為一種數(shù)學(xué)覆蓋方式生成NMM流形覆蓋系統(tǒng)，避免采用規(guī)則網(wǎng)格生成流形覆蓋系統(tǒng)時(shí)可能產(chǎn)生的細(xì)長(zhǎng)或極小尺寸流形單元，并提出了將單個(gè)材料體描述為單個(gè)流形單元的獨(dú)立覆蓋方法。通過(guò)多種流形覆蓋方式相耦合，模擬了混凝土方件的細(xì)觀變形破壞，以及復(fù)雜結(jié)構(gòu)巖石邊坡的失穩(wěn)破壞過(guò)程。結(jié)果表明，獨(dú)立覆蓋方式的采用，降低了模型的生成過(guò)程以及最終生成模型的復(fù)雜程度，減少了計(jì)算過(guò)程中的接觸數(shù)量以及流形單元數(shù)量，提高了計(jì)算效率。獨(dú)立覆蓋與有限覆蓋的耦合使用，合理地描述了混凝土的細(xì)觀結(jié)構(gòu)和巖質(zhì)邊坡的復(fù)雜構(gòu)造，NMM模擬較好地再現(xiàn)了相應(yīng)的物理力學(xué)過(guò)程。文中提出的獨(dú)立覆蓋方法，適用于單體變形不重要的離散體的模擬，也可用于模型中加載板、約束邊框等單體的描述。

參考文獻(xiàn)：

[1] 劉凱欣，高凌天.離散元法研究的評(píng)述[J].力學(xué)進(jìn)展，2003，33（4）：483-490.

LIU Kaixin， GAO Lingtian. A review on the discrete element method[J]. Advances in Mechanics， 2003， 33（4）： 483-490.（in Chinese）

[2] 張國(guó)新，趙妍，石根華，等.模擬巖石邊坡傾倒破壞的數(shù)值流形法[J].巖土工程學(xué)報(bào)，2007，29（6）：800-805.

ZHANG Guoxin， ZHAO Yan， SHI Genhua， et al. Toppling failure simulation of rock slopes by numerical manifold method[J]. Chinese Journal of Geotechnical Engineering， 2007， 29（6）： 800-805. （in Chinese）

[3] Mahabadi O K， Lisjak A， Munjiza A， et al. Y-Geo： New combined finite-discrete element numerical code for geomechanical applications[J]. International Journal of Geomechanics， 2012， 12（6）： 676-688.

[4] Ma G W， An X M， He L. The numerical manifold method： a review[J]. International Journal of Computational Methods， 2010， 7（1）： 1-32.

[5] 蔡永昌，張湘?zhèn)?流形方法的矩形覆蓋系統(tǒng)及其全自動(dòng)生成算法[J].重慶大學(xué)學(xué)報(bào)（自然科學(xué)版），2001，24（1）：42-46.

CAI Yongchang， ZHANG Xiangwei. Rectangular cover system of manifold method and its auto mesh algorithm[J]. Journal of Chongqing University （Natural Science Edition）， 2001，24（1）： 42-46. （in Chinese）

[6] 陳剛，劉佑榮.流形元覆蓋系統(tǒng)的有向圖遍歷生成算法研究[J].巖石力學(xué)與工程學(xué)報(bào)，2003，22（5）：711-716.

CHEN Gang， LIU Yourong. Generation of cover system for numerical manifold in travel theory of the oriented graph[J]. Chinese Journal of Rock Mechanics and Engineering， 2003， 22（5）： 711-716. （in Chinese）

[7] 李海楓，張國(guó)新，石根華，等.流形切割及有限元網(wǎng)格覆蓋下的三維流形單元生成[J].巖石力學(xué)與工程學(xué)報(bào)，2010，29（4）：731-742.

LI Haifeng， ZHANG Guoxin， SHI Genhua， et al. Manifold cut and generation of three-dimensional manifold element under FE mesh cover[J]. Chinese Journal of Rock Mechanics and Engineering， 2010， 29（4）： 731-742.（in Chinese）

[8] 蔡永昌，劉高揚(yáng).基于獨(dú)立覆蓋的高階流形方法[J].同濟(jì)大學(xué)學(xué)報(bào)（自然科學(xué)版），2015，43（12）：1794-1800.

CAI Yongchang， LIU Gaoyang. High-order manifold method with independent covers[J]. Journal of Tongji University （Natural Science）， 2015， 43（12）： 1794-1800.（in Chinese）

[9] 蘇海東，頡志強(qiáng)，龔亞琦，等.基于獨(dú)立覆蓋的流形法的收斂性及覆蓋網(wǎng)格特性[J].長(zhǎng)江科學(xué)院院報(bào)，2016，33（2）：131-136.

SU Haidong， XIE Zhiqiang， GONG Yaqi， et al. Characteristics of convergence and cover mesh in numerical manifold method based on independent covers[J]. Journal of Yangtze River Scientific Research Institute， 2016， 33（2）： 131-136.（in Chinese）

[10] 黃濤.聚合物粘結(jié)炸藥損傷與斷裂的流形元數(shù)值模擬[D].北京：北京理工大學(xué)，2006.

HUANG Tao. Numerical simulation of damage and fracture of polymer explosive by manifold method[D]. Beijing： Beijing Institute of Technology， 2006. （in Chinese）

[11] 趙妍.基于流形元法的裂縫追蹤技術(shù)在結(jié)構(gòu)破壞模擬中的應(yīng)用研究[D].北京：中國(guó)水利水電科學(xué)研究院，2008.

ZHAO Yan. Crack tracking techniques based on the manifold method and its applications in structure failure simula-tions[D]. Beijing： China Institute of Water Resources and Hydropower Research， 2008. （in Chinese）

[12] 黃濤，陳鵬萬(wàn)，張國(guó)新，等.巖石雙孔爆破過(guò)程的流形元法模擬[J].爆炸與沖擊，2006，26（5）：434-440.

HUANG Tao， CHEN Pengwan， ZHANG Guoxin， et al. Numerical simulation of two hole blasting using numerical manifold method[J]. Explosion and Shock Waves， 2006， 26（5）： 434-440. （in Chinese）

[13] Ning Y J， An X M， Ma G W. Footwall slope stability analysis with the numerical manifold method[J]. International Journal of Rock Mechanics and Mining Sciences， 2011， 48（6）： 964-975.

[14] An X M， Ning Y J， Ma G W， et al. Modeling progressive failures in rock slopes with non-persistent joints using the numerical manifold method[J]. International Journal for Numerical and Analytical Methods in Geomechanics， 2014， 38（7）： 679-701.

[15] Goodman R E， Kieffer D S. Behavior of rock in slopes[J]. Journal of Geotechnical and Geoenvironmental Engineering， 2000， 126（8）： 675-684.

[16] Goodman R E . Toppling-A fundamental failure mode in discontinuous materials—description and analysis[C]∥[S.l.]Geo-Congress 2013： Stability and Performance of Slopes and Embankments III. 2013： 2338-2368.

（編輯 陳移峰）