這里使用之前帖子中的單軸算例，進(jìn)行了一點(diǎn)修改，前三個(gè)文件保持不變，在單軸前的處理中，刪除墻平衡的時(shí)候，趁機(jī)挖了一個(gè)裂縫。然后進(jìn)行單軸。

圖中便是加上的部分，kuandu是豎向裂縫的額長(zhǎng)度，houdu是裂縫的法向厚度的一半。

zhongxinX和zhongxinY是裂縫的中心位置。

輸入上面四個(gè)參數(shù)便可以生成裂縫。

---

配圖1-6：裂紋發(fā)展圖

可以看到裂紋不像傳統(tǒng)單軸是從兩端生成裂紋，而是在裂縫附近先生成裂紋。以一個(gè)X的方向向四周擴(kuò)展。

---

配圖7-11：塊體圖

這里可以看到宏觀裂紋的擴(kuò)展情況，宏觀上裂紋在裂縫的兩端向受壓端延伸。

形成貫穿的裂縫后，受壓端較窄的那部分又發(fā)展新的宏觀裂紋，這個(gè)和微觀裂紋圖是一致的。

---

配圖12： 應(yīng)力應(yīng)變曲線

這里可以看到在應(yīng)力下降后有一個(gè)階梯，懷疑是兩端窄的部分產(chǎn)生裂縫后形成之后劇烈的衰減，沒有比對(duì)過，需要再驗(yàn)證一下。

---

配圖13： 裂紋演變曲線

這里的巖體為高摩擦角的性質(zhì)，以剪壞為主，拉壞只有很少的一部分。

這里的代碼是在單軸的基礎(chǔ)上修改的，唯一修改的部分就是一開始說(shuō)明的，這里就不貼上代碼了。

后面會(huì)寫直剪、單剪、滾刀、邊坡等，有時(shí)間就做一做，比較慢。