纖維骨料細觀混凝土的案例
三維隨機纖維-球體骨料細觀混凝土模型/細觀混凝土/纖維混凝土 ￥369
在前面的視頻中已經有了三維隨機球體骨料的細觀模型建立方法和靜力抗壓仿真分析的方法（課程封面如圖）。
簡單三維球體細觀模型的課程已更新一年半有余，后臺私信和留言中有不少同學咨詢“怎么添加纖維/有沒有纖維混凝土的模型”，由于之前一直在忙著干其他的事情，最近心血來潮，那么今天就寫了這個帖子，分享自己關于在三維細觀混凝土模型中添加纖維的一些思路和腳本建立方法。
三維隨機纖維-球體骨料細觀混凝土模型——四相（砂漿、骨料、ITZ和纖維）組分的復合混凝土模型組成效果和仿真效果如下。
在前邊的三維隨機球體模型中，投放骨料的思路十分簡單，只需要使用python代碼隨機生成一個球體特征數據，再與已存儲的球體數據進行判斷（判斷球心距與兩球體半徑之和）即可。若球心距大于半徑之和則存儲在骨料數據庫中，進行下一個球體的生成與判斷；若小于，則不存儲，進行下一個球體的生成判斷。最后使用python與abaqus之間的接口，把數據轉化為圖形即可。
把這個思路放到纖維與骨料之間的判斷中來，似乎也能進行相應的判斷。只需要生成隨機的纖維，用纖維端點坐標與骨料球心坐標，計算出球心到直線的距離就可以了，如下圖所示。
使用點到直線的距離公式判斷球體與纖維的相交，這樣看著好像沒啥問題，但其實纖維能在混凝土中分布的區域已經大大縮小了。如下圖所示，當纖維的方向指向骨料時，雖然纖維與骨料并沒有相交，但簡單地使用點到直線的距離公式，會被判斷為相交狀態，這根纖維就將被認為不能放在混凝土中。
展開 Abaqus纖維混凝土3D 泡沫混凝土 三維隨機幾何 三維混凝土細觀 多面體骨料建模
模型實例
以下是Abaqus內纖維混凝土的模型，纖維是采用三維圓柱體模擬的，混凝土內的骨料采用的是實體的球體。纖維及骨料均可設置不同的尺寸，并且各類型的數目不受限制，即可設置多種纖維及球體骨料大小。
研究進展
在Abaqus內建立混凝土細觀模型，如鋼纖維混凝土、不干涉球體骨料、多面體骨料模型等，是進行混凝土性能研究的主流方法之一。而在進行Abaqus混凝土細觀模擬時，隨機骨料及隨機纖維等幾何模型的構件是主要的難點所在。
為了在Abaqus內建立混凝土模型，有學者采用Abaqus命令的方式，但這需要有一定的程序設計基礎，并且需要反復改參、調試，極為不便。也有采用Abaqus混凝土建模插件實現的方式，這極大的節省了模型建立的耗時，如Abaqus混凝土多邊形或Abaqus混凝土三維球體骨料插件等，但其實現的模型較為簡單，幾何模型單一。
建模方案
這里介紹一種通過AutoCAD軟件建立纖維混凝土三維模型后導入到Abaqus內的方式。可實現多種混凝土模型的快速構建。CAD導入Abaqus的方法簡單，將CAD文件輸出為.sat格式，然后在Abaqus內選擇導入部件，選擇對應的.sat文件即可。
下面是通過該方法建立的Abaqus隨機幾何模型。
插件介紹
本插件可以生成多種形式的隨機三維幾何，用于Abaqus混凝土模型的建立，也可用于再生骨料混凝土、泡沫混凝土、加氣混凝土等方面。理論上講，只要幾何存在相似性，可進行模型簡化的，均可采用這種方式進行建模。
插件的詳細介紹及下載見下方鏈接：
CAD隨機幾何3D插件
展開 comsol聯合Matlab生成纖維、骨料細觀混凝土模型（附球形骨料代碼、纖維代碼） ￥99
本課程旨在介紹如何利用matlab與comsol連接，并利用matlab語言批量對comsol進行幾何建模，生成復雜、隨機的模型，如纖維、骨料等。可根據需要進行開裂分析等，效果圖如下：
寫在前面：[首先確定自己已安裝COMSOL Multiphysics 5.6 with MATLAB，
如果電腦上先安裝comsol，再安裝matlab的話一般不會出現這個程序。
解決方法：卸載已安裝的comsol，先安裝matlab，再安裝comsol，在安裝過程中會提示關聯matlab，安裝完成后即可出現該程序。]
1
comsol與matlab連接之隨機球形骨料生成腳本教學
（1）運行COMSOL Multiphysics 5.6，并以此點擊模型向導→三維→完成。此時在模型開發器中右鍵幾何，選擇球體
此時我們可以定義球體半徑為2，坐標[x,y,z]為[3，4，5]并構建選定對象，如下圖所示
至此為止，我們已在comsol中生成了1個球體，那么接下來介紹如何利用Matlab生成一定數量和半徑的球體。
展開 ABAQUS纖維骨料ITZ孔隙細觀混凝土模型
本案例介紹在ABAQUS內建立包含骨料、界面過渡區、纖維、孔隙在內的多相材料二維纖維混凝土細觀模型。
采用CAD纖維混凝土2D插件在AutoCAD內建立二維纖維混凝土模型，纖維混凝土模型的不同組分已分圖層進行繪制，需要對不同圖層內容分別另存為dxf文件。
在Abaqus內將建立的模型文件以草圖的形式分別導入，這里只展示骨料草圖內容。
通過草圖建立二維部件模型，這里展示的是纖維部件。
其中界面過渡區部件需要將直接導入的過渡區與骨料進行一次切割幾何完成。
將各個部件裝配為整體，將模型中的孔隙部分進行挖除。
可對模型中的不同部分分別設置材料。
展開 
ANSYS\ABAQUS纖維混凝土細觀骨料模型建立及網格劃分 ￥1.1
《基于三維隨機細觀模型的珊瑚混凝土力學性能模擬》一文中建立了考慮界面層（ITZ）、骨料、砂漿的三相混凝土模型，并采用“背景投影法（網格映射法）”建立了六面體非均質混凝土有限元模型。</p><p class="ql-align-justify">相比均質有限元模型，非均質有限元模型的仿真結果可信度更高，仿真效果更好，與實際破壞情況更為吻合，該方法具有廣泛的運用前景，可用于靜態力學試驗、動態力學試驗、爆破領域、建筑結構領域等。</p><p class="ql-align-justify">在有限元分析中，網格質量的好壞極大程度影響模擬的收斂性，尤其對于顯式動態分析案例中，為了避免網格畸變導致計算時間長、計算結果不收斂等問題，大多采用六面體網格進行計算。因此，本文對非均質纖維混凝土模型分別進行了四面體網格、六面體網格劃分的對比，并對該類網格劃分問題的步驟進行闡述。</p><p class="ql-align-justify">步驟一：采用Python、Fortran、APDL等編程語言生成隨機骨料及纖維，判定骨料與骨料之間，纖維與纖維之間，纖維與骨料之間互不侵入。基于此，生成骨料半徑、中心坐標，纖維起始點和終止點的坐標。</p><p class="ql-align-justify">步驟二：將坐標信息導入ANSYS或ABAQUS中，結合軟件自帶建模語言進行建模及網格劃分，四面體網格可通過hypermesh進行精細網格劃分，也可采用自編網格投影法進行六面體網格劃分，不同方法均存在利弊。六面體網格計算時間大量縮短，但骨料形狀為類球體，是否能投影為球體與單元網格尺寸大小有關，四面體網格計算時間較長，劃分形狀與球體基本一致。</p><p class="ql-align-justify">步驟三：進行材料、單元幅值，開展不同有限元分析。
展開 ANSYS混凝土三維隨機骨料 混凝土細觀 隨機球體 顆粒增強復合材料建模
研究進展
通過ANSYS進行混凝土細觀模型的構建是進行混凝土性能分析的有效方法，在ANSYS內構建混凝土細觀模型是分析的前提。現階段在ANSYS內進行隨機混凝土模型構建的主流方法是通過APDL命令流等形式，這要求研究者應具有一定的程序設計能力。
為了方便快捷的構建出混凝土細觀幾何模型，這里提出另一種建模方案，通過AutoCAD模型導入的方式，實現無編程構建混凝土隨機骨料。
模型構建
1、CAD模型生成
首先采用CAD隨機球體顆粒插件在AutoCAD內構建三維球體幾何模型：
插件可指定生成隨機分布的不相交的球體顆粒，同時生成與球體顆粒裝配的帶有孔洞的長方體基體。同時對顆粒的粒徑大小、比例等都能進行控制。
將生成的三維球體幾何模型導出為.sat格式文件備用。
2、ANSYS Workbench 導入
打開ANSYS Workbench，在幾何內進行導入預先保存的.sat文件：
后續進行網格劃分等操作，在ANSYS Workbench內進行即可：
插件下載
建模用到的CAD插件下載：
CAD隨機球體顆粒插件
展開 ANSYS隨機骨料 纖維混凝土 三維隨機纖維骨料 隨機纖維 隨機裂縫 隨機幾何模型
1、ANSYS三維纖維骨料混凝土：
2、ANSYS球形試件隨機模型：
3、ANSYS隨機裂縫巖石節理裂隙
建模插件：
CAD隨機幾何3D插件
COMSOL隨機多面體骨料 三維凸多面體骨料 無規則孔隙 三維混凝土細觀 三維骨料模型
混凝土模型
三維混凝土細觀模型的建立是進行混凝土性能模擬的有效方法，而在comsol建模過程中隨機凸多面體骨料的生成是幾何模型的難點。這里提供一種快速高效的三維凸多面體骨料建模的方案，以實現不同集配的混凝土模型。
建模教程
首先采用CAD隨機多面體3D插件在AutoCAD內生成所需要的三維混凝土細觀模型。
將該模型分圖層導出為.iges格式文件，這里分圖層導出是為了可以分部件導入到comsol軟件內，更方便材料賦值等操作。
本模型共導出四個iges文件，分別是帶有多面體孔洞的基體材料以及三種不同粒徑的多面體。
然后將iges文件分別導入到comsol內，這里建議每導入一部分后緊接著進行材料賦值操作，材料賦值完成并將該部分隱藏，然后再導入另一部分，否則可能會出現材料賦值難以選取的問題。
最后進行網格劃分、邊界條件、模擬計算等操作即可。
這里再放一張賦值不同材料后的模型：
插件下載
CAD隨機多面體3D插件
模型樣圖
隨機多面體骨料_AbyssFish.rar
展開 ANSYS隨機多面體骨料 三維多面體投放 隨機骨料混凝土細觀模型
ANSYS隨機多面體骨料模型，采取精確的干涉判斷，采用多面體相交判別程序，不同于常見的球體干涉，本程序可達到更好的隨機度，以實現大粒徑與小粒徑的匹配度。
輕骨料混凝土細觀損傷演化分析
(a)輕骨料混凝土骨料內聚力單元 (b)界面過渡區內聚力單元 (c)砂漿內聚力單元
圖2 輕骨料混凝土細觀模型中內聚力單元示意圖
2.1.3 輕骨料混凝土細觀模型的建立
采用1.5mm的單元尺寸網格對生成的多邊形輕骨料混凝土模型進行單元劃分，同時在輕骨料混凝土模型的全局區域內嵌入內聚力單元，這樣將會形成足夠精細的裂紋擴展區，使得輕骨料混凝土的裂紋擴展與真實混凝土裂紋擴展更加類似，由于內聚力單元非線性的損傷演化，為提高輕骨料混凝土細觀模型分析的收斂性，采用Abaqus動力顯式模塊進行準靜態分析，且加載過程中采用速度平滑加載來實現準靜態加載過程。
輕骨料混凝土細觀模型尺寸為100mm×300mm，輕骨料混凝土細觀數值模型單軸受壓如圖3所示。
圖3 輕骨料混凝土細觀數值模型單軸受壓示意圖
2.2 輕骨料混凝土細觀組分本構關系
砂漿與骨料之間的界面過渡區（ITZ）、輕骨料內界面和砂漿內界面（MII）均采用內聚力單元模擬，內聚力模型采用了適用于模擬裂縫的牽引力-分離本構關系，本文采用的是沒有幾何厚度的內聚力單元，因此在進行數值分析時，需要假設一個初始厚度，在ABAQUS軟件中默認零厚度內聚力單元的初始厚度為1。由于開裂位移很難測定，本文通過斷裂能參數 來表示位移參數，斷裂能與開裂位移的關系如下：
2.3 輕骨料混凝土單軸受壓數值模擬
輕骨料混凝土受壓模擬破壞結果如圖4所示，本文模擬得到的輕骨料混凝土應力應變曲線與王振宇等人開展的輕骨料混凝土試驗結果進行對比，如圖5所示。可知通過細觀模型得到輕骨料混凝土應力應變曲線與試驗吻合較好，表明模型能夠有效預測輕骨料混凝土的力學響應以及參數選取的合理性。
展開 ABAQUS三維多面體骨料密堆積混凝土細觀建模
ABAQUS三維多面體骨料密堆積建模通過重力堆積算法構建混凝土細觀結構，克服了傳統隨機分布模型與實際骨料沉降行為的偏差，更精準反映骨料在混凝土中的分布特征，可實現高骨料占比下的力學響應模擬，為混凝土損傷機理研究、材料參數標定及多尺度耦合分析提供可靠依據。本案例介紹在ABAQUS內建立三維混凝土多面體骨料重力密堆積模型。
混凝土細觀骨料堆積模型采用CAD多面體密堆積_圓柱體試件3D插件生成，在AutoCAD內建立模型后將骨料導出為iges格式文件。
在ABAQUS CAE中將骨料模型導入建立部件，如需對不同骨料分別設置材料屬性也可將骨料分為三個部件分別導入。
對混凝土細觀模型中的各組分設置不同的材料屬性，在混凝土細觀研究中可選用EasyCDP_Mortar&ITZ插件生成混凝土損傷塑性模型材料。
采用四面體單元進行混凝土骨料的網格劃分。
建立不同骨料部件的裝配。
或是通過AutoCAD軟件對骨料進行處理后導出卵石狀骨料模型，并導入到ABAQUS內。
如需進行三維混凝土細觀受壓損傷開裂方面的仿真模擬也可參考ABAQUS三維混凝土細觀受壓損傷斷裂數值模擬視頻教程。
展開 
ABAQUS隨機多邊形骨料及ITZ細觀混凝土CDP模型
混凝土的強度很大程度上取決于粗骨料與水泥砂漿之間的界面過渡區（ITZ）。本案例在ABAQUS內建立隨機多邊形骨料模型，并設置界面過渡區部件，采用CDP材料建立骨料、砂漿、ITZ三相混凝土細觀模型，并研究模型的軸壓破壞情況。
混凝土細觀模型采用CAD隨機多邊形顆粒插件建模生成，將插件生成的CAD文件按照不同圖層內容分三份以草圖的形式導入到ABAQUS內。
在ABAQUS內利用導入的草圖分別建立砂漿、界面過渡區、多邊形骨料三部分部件并進行裝配。
在本案例中采用EasyCDP插件快速生成不同強度的混凝土塑性損傷（CDP）材料，并將其分別指派到砂漿及ITZ部件，多邊形骨料不考慮損傷破壞。
選擇動力，顯式分析步，設置時間長度為5，將模型下部設置為固定約束，上部添加向下的位移作為載荷。
進行網格劃分，單元尺寸的設置推薦小于ITZ厚度及插件中設置的顆粒最小間距的較小值。
提交分析查看結果。
CAD隨機多邊形顆粒
https://www.yqgqt.org.cn/post/1787116
展開 ABAQUS三維隨機多面體骨料及ITZ細觀混凝土建模
在混凝土結構中，粗骨料與水泥砂漿之間的界面過渡區（ITZ）是決定混凝土承受荷載能力的關鍵因素之一,研究表明，混凝土在承受外力時的破壞主要發生在界面過渡區的損傷。本案例介紹通過CAD隨機多面體&過渡區3D插件建立隨機分布的三維多面體骨料及界面過渡區細觀混凝土模型，并將模型導入ABAQUS內進行多相材料的指定。
在AutoCAD軟件內，采用CAD隨機多面體&過渡區3D V1.0插件建立隨機投放的多面體骨料、界面過渡區(ITZ)部件及水泥砂漿基體三維模型。插件在建模時已將不同組分分圖層繪制，方便批量導出及后續的材料指定。這里將砂漿、ITZ、不同粒徑的骨料分別導出為iges格式文件。
將導出的iges模型文件以部件的形式導入到ABAQUS內。
可對多面體骨料、ITZ、砂漿分別手動設置材料，也可采用EasyCDP插件快速生成混凝土損傷塑性材料后，再指派到砂漿及ITZ部件。
在裝配中將兩組多面體骨料、水泥砂漿、過渡區部件進行裝配，完成三維細觀混凝土模型的建模。
CAD_隨機多面體&過渡區插件
https://www.yqgqt.org.cn/post/1909149
展開 Abaqus隨機圓形骨料ITZ細觀混凝土CDP軸壓模型
混凝土在細觀層面上由水泥砂漿、粗骨料和界面過渡區（ITZ）組成，在Abaqus內基于粗骨料隨機投放建立混凝土細觀模型，是研究混凝土軸壓下的本構關系及損傷演化的有效方法。本案例建立隨機圓形粗骨料及實體界面過渡區，對二維細觀混凝土在單軸壓縮下的力學行為進行有限元模擬，展示混凝土的破壞形態。
在Abaqus CAE軟件內，采用AbyssFish RandomAggregate V3.2插件建立圓形粗骨料、實體界面過渡區、水泥砂漿三部件混凝土細觀模型。
對各部件進行材料截面的指派，其中水泥砂漿部分采用CDP模型；界面過渡區采用弱化的砂漿模型；骨料部分不考慮其損傷破壞。
為模擬模型的軸壓狀態，新建長方形部件作為壓力機的鋼制壓板，設置壓板與混凝土試件間的接觸，并將下側壓板固定，上側壓板添加豎直向下的位移。
對模型進行網格劃分。
創建作業提交分析并查看結果。
展開 ANSYS/LS-DYNA三相細觀骨料混凝土SHPB沖擊壓縮模擬
關于SHPB數值模擬的研究已較為深入，模擬優勢主要在于可通過修正參數使模擬結果與實際一致，以此為基礎對材料的動態破壞過程及更為復雜的工況進行模擬研究，主要研究對象主要分為混凝土、巖石、金屬、陶瓷等材料，并通過LS-DYNA中的RHT、HJC、JC、K&C、CSC等材料模型來模擬上述材料在中高、高應變率荷載作用下裂紋擴展及損傷規律，試件往往采用的是均質模型。
近年來，關于非均質模型的研究已取得一些進展：
1.《Study of concrete damage mechanism under hydrostatic pressure by numerical simulations》一文中建立了考慮骨料、砂漿的兩相混凝土模型，并采用“背景投影法（網格映射法）”建立了六面體非均質混凝土有限元模型。
2.《3D mesoscopic investigation of the specimen aspect-ratio effect on the compressive behavior of coral aggregate concrete》一文中建立了考慮界面層（ITZ）、骨料、砂漿的三相混凝土模型，并采用“背景投影法（網格映射法）”建立了六面體非均質混凝土有限元模型。
3.《基于三維隨機細觀模型的珊瑚混凝土力學性能模擬》一文中建立了考慮界面層（ITZ）、骨料、砂漿的三相混凝土模型，并采用“背景投影法（網格映射法）”建立了六面體非均質混凝土有限元模型。
相比均質有限元模型，非均質有限元模型的仿真結果可信度更高，仿真效果更好，與實際破壞情況更為吻合，該方法具有廣泛的運用前景，可用于靜態力學試驗、動態力學試驗、爆破領域、建筑結構領域等。
展開 