多相介質(zhì)仿真的案例
COMSOL顆粒夾雜多孔介質(zhì)多相材料達(dá)西滲流模擬
在實(shí)際工程中滲流路徑往往不是單一材料，如滲流發(fā)生在夾雜碎石的土體中，這就造成滲流的復(fù)雜性。這里采用兩項(xiàng)材料通過COMSOL達(dá)西定律模塊對(duì)滲流進(jìn)行模擬。
模型采用CAD隨機(jī)球體顆粒&過渡區(qū)插件建立后導(dǎo)入到COMSOL軟件內(nèi)。
模型包括滲流發(fā)生的外側(cè)基體、內(nèi)部顆粒、顆粒及基體過渡區(qū)（ITZ）三部分組成，由于內(nèi)部顆粒的滲透系數(shù)遠(yuǎn)小于基體，因此可將其省略，邊界置為無流動(dòng)。設(shè)置過渡區(qū)的目的是在實(shí)際情況中，土體及內(nèi)部碎石顆粒間往往會(huì)有孔隙，這就造成了接觸面的實(shí)際滲透率遠(yuǎn)高于土體，模型剖切面如下。
模型設(shè)置左右兩側(cè)的水頭差，最終壓力及流速模擬結(jié)果如下。
土石混合體多相多物理場(chǎng)耦合數(shù)值仿真 ￥5000
基于COMSOL軟件對(duì)土石混合體進(jìn)行了數(shù)值仿真，考慮了土石混合體孔隙變化，細(xì)顆粒侵蝕，骨架結(jié)構(gòu)變形，此問題是一個(gè)多場(chǎng)（滲流場(chǎng)、變形場(chǎng)、應(yīng)力場(chǎng)、損傷場(chǎng)）多相介質(zhì)（土顆粒集合體，塊石，空隙，孔隙）耦合的復(fù)雜問題。仿真結(jié)果如圖2所示。
圖1 幾何模型
顆粒運(yùn)動(dòng)分布
應(yīng)力分布
孔隙滲流下的細(xì)顆粒遷移運(yùn)動(dòng)
圖2 數(shù)值仿真結(jié)果
感興趣的朋友可下載模型源文件，歡迎合作交流
LS-DYNA 基于CT掃描的多相材料動(dòng)力學(xué)仿真
本貼的內(nèi)容是關(guān)于LS-DYNA的高級(jí)應(yīng)用，基于真實(shí)多相材料的物質(zhì)分布進(jìn)行建模，從而真實(shí)模擬多相物質(zhì)的動(dòng)力學(xué)特性的技術(shù)。
以混凝土材料為例。
首先，獲取CT數(shù)據(jù)。這邊有一個(gè)關(guān)于CT數(shù)據(jù)的開源網(wǎng)站https://www.digitalrocksportal.org/
獲取到raw文件后，用avizo軟件打開，并生成二值化的tif圖。
接下來，就是實(shí)現(xiàn)多圖向網(wǎng)格生成的關(guān)鍵步驟了，可以參考 淵魚 大佬的帖子，或者自己編程。
一般無法直接生成lsdyna需要的k文件，需要借助其他建模或仿真軟件生成。
然后按照lsdyna的一般建模思路計(jì)算。
最后的損傷結(jié)果如圖。
展開 晶體塑性有限元仿真入門(2)--BCC、FCC、HCP晶格材料以及多相材料的有限元模擬
晶體塑性有限元仿真入門(2)--BCC、FCC、HCP晶格材料以及多相材料的有限元模擬
這篇文章講解如何使用晶體塑性有限元方法(CPFEM)進(jìn)行不同晶格材料以及多相材料的變形模擬，CPFEM是基于商業(yè)有限元軟件ABAQUS完成的建模，晶體塑性本構(gòu)模型是使用的開源的UMAT用戶子程序(源碼和inp文件見附件)。采用CPFEM模擬了面心立方結(jié)構(gòu)(FCC)、體心立方結(jié)構(gòu)(BCC)和密排六方結(jié)構(gòu)(HCP)的單晶、多晶及多相材料受到外部載荷時(shí)的力學(xué)響應(yīng)。基于滑移原理的晶體變形理論，隨著變形的進(jìn)行各滑移系統(tǒng)的臨界剪應(yīng)力都會(huì)增大，CPFEM將捕捉到材料的力學(xué)響應(yīng)(應(yīng)力-應(yīng)變曲線)。這些應(yīng)力-應(yīng)變數(shù)據(jù)有助于從根本上理解晶粒尺度下金屬變形的性質(zhì)。
首先我們從一個(gè)簡(jiǎn)單的FCC晶格材料的例子入手，講解如何進(jìn)行有限元模型的創(chuàng)建，從完全新手的角度出發(fā)，一步步講解如何建模，賦予材料和處理仿真結(jié)果。
本文章包括以下八個(gè)部分：
1) FCC晶格材料的變形模擬-單晶體
2) FCC晶格材料的變形模擬-多晶體
3) BCC晶格材料的變形模擬-單晶體
4) BCC晶格材料的變形模擬-多晶體
5) HCP晶格材料的變形模擬-單晶體
6) HCP晶格材料的變形模擬-多晶體
7) 多相材料的變形模擬
8) 參考資料
1. FCC晶格材料的變形模擬-單晶體
幾何模型
如圖1.1a在草圖里繪制R0.015mm的圓形，拉伸0.05mm，最后得到如圖1.1b所示的圓柱體(R0.015mm&H0.05mm)。
展開 
