abaqus單元組
關(guān)注創(chuàng)建者：王靖雯 創(chuàng)建時(shí)間：2023-02-27
abaqus單元組的實(shí)例教程
在Abaqus中，用戶自定義元素子程序（UEL）的開發(fā)需要遵循一些特定的規(guī)則和約定。其中一個(gè)關(guān)鍵約定是關(guān)于子程序中的RHS（右手邊）向量的維數(shù)，我最近在嘗試用UEL做一些二次開發(fā)，也發(fā)現(xiàn)了RHS向量的維數(shù)比單元的總自由度數(shù)多了4個(gè)這一現(xiàn)象，結(jié)合在站內(nèi)一些同行的猜測，我認(rèn)為這可能是由于Abaqus的內(nèi)部工作方式所導(dǎo)致的。
首先，有同行懷疑是因?yàn)橛昧怂墓?jié)點(diǎn)單元，所以多了四個(gè)，我開發(fā)的單元是12個(gè)節(jié)點(diǎn)的，依然多了四個(gè)維度，因此排除是單元內(nèi)節(jié)點(diǎn)個(gè)數(shù)導(dǎo)致的。
在Abaqus中，RHS向量的維數(shù)實(shí)際上包括了除了單元的位移自由度外的其他項(xiàng)，這些項(xiàng)用于處理多種情況，例如：
體積力和表面力的計(jì)算：RHS向量可能包括用于計(jì)算體積力和表面力的額外自由度。這些自由度用于存儲(chǔ)單元內(nèi)的體積力和表面力的貢獻(xiàn)。
約束和邊界條件：Abaqus可能需要額外的自由度來處理約束條件和邊界條件，以確保數(shù)值穩(wěn)定性和正確的求解結(jié)果。
內(nèi)部狀態(tài)變量：某些材料模型和非線性分析可能需要存儲(chǔ)和更新一些內(nèi)部狀態(tài)變量，這些變量也可以占用RHS向量中的額外位置。
因此，RHS向量的維數(shù)不僅僅包括單元的位移自由度，還包括其他與分析和模型特性相關(guān)的項(xiàng)。這是Abaqus設(shè)計(jì)的一部分，旨在確保通用性和可擴(kuò)展性，以處理各種復(fù)雜的問題。其中我認(rèn)為可擴(kuò)展性是一項(xiàng)比較重要的應(yīng)用，我目前所做的工作可能會(huì)利用到這一點(diǎn)，等有結(jié)果了會(huì)繼續(xù)更新。
展開 
abaqus單元組的相關(guān)專題、標(biāo)簽、搜索
abaqus單元組的最新內(nèi)容
結(jié)合作者的理論（尤其是分段線性化和應(yīng)力驅(qū)動(dòng)的求解思路）我們可以把獨(dú)立的vpsc子程序編寫進(jìn)abaqus里面，為了避免復(fù)雜的雅可比推導(dǎo)，以及適用各種復(fù)雜的變形工況，推薦使用abaqus的顯式求解器，即vumat程序
以下展示一個(gè)使用vpsc-鎂合金本構(gòu)模型，模擬包含1個(gè)單元，單元包好100個(gè)晶粒在RD方向壓縮20%的模擬效果（原始模型參數(shù)取自vpsc官方案例,為了減少計(jì)算時(shí)間使用高應(yīng)變率進(jìn)行計(jì)算，
用戶可獨(dú)立設(shè)置中心加密區(qū)與外圍粗化區(qū)的單元尺寸，兼顧計(jì)算精度與效率。所有實(shí)體層采用 C3D8R 減縮積分單元并激活單元刪除，內(nèi)聚力層采用 COH3D8 單元，沖頭則使用離散剛體單元 R3D4。網(wǎng)格劃分基于掃掠技術(shù)（Advancing Front）生成。
考慮GND的大變形冷軋模擬4個(gè)月前
其中流動(dòng)方程使用經(jīng)典的唯象流動(dòng)方程：
硬化模型則使用了同時(shí)考慮SSD和GND的位錯(cuò)密度硬化模型：
作者構(gòu)建了包含 520 個(gè)晶粒的三維 RVE（NiTi 基體晶粒 442 個(gè)、β-Nb 晶粒 78 個(gè)），并在 ABAQUS 中進(jìn)行單道次軋制變形20%的模擬。
我只需要從均勻的基體網(wǎng)格中判斷，哪些單元被已有的纖維網(wǎng)格包裹了，就可以得到纖維的體素網(wǎng)格。</p><p>實(shí)際上，體素方法本身并不高端，它本質(zhì)是一直暴力搜索的方法。從技術(shù)上來說，沒有什么美感。</p><p>唯一需要注意的是，為了保證纖維足夠的特征，我們需要先創(chuàng)建大規(guī)模的基體單元。我目前測試下來，對于單胞級(jí)的，最好也要有100W網(wǎng)格。暴力搜索方法就會(huì)很慢。
ABAQUS的六面體網(wǎng)格，只要知道構(gòu)造網(wǎng)格的8個(gè)節(jié)點(diǎn)和排序規(guī)律，就可以用一行字符創(chuàng)建出網(wǎng)格。
類比到修真世界，就是口訣（代碼）+符咒（字符串）+陣法（數(shù)據(jù)排布）。
機(jī)織復(fù)合材料看似復(fù)雜、幾何參數(shù)多。但是由于周期性排布的特點(diǎn)，總能找到一個(gè)代表性單元。只要把握住這個(gè)代表性單元，就完成了一般的工作。
定義分析步，打開幾何非線性開關(guān)，設(shè)置步長為100s，每間隔1s輸出一組結(jié)果，采用動(dòng)力學(xué)隱式求解方法。
4.計(jì)算結(jié)果
通過ABAQUS有限元計(jì)算可以得到壓電復(fù)合結(jié)構(gòu)的正弦振動(dòng)響應(yīng)結(jié)果，如圖4所示，動(dòng)態(tài)圖展示了壓電復(fù)合結(jié)構(gòu)在交流電作用下動(dòng)力學(xué)響應(yīng)。圖5為基體板自由端某一節(jié)點(diǎn)位移時(shí)域曲線。
ABAQUS三維多面體骨料密堆積混凝土細(xì)觀建模5個(gè)月前
采用四面體單元進(jìn)行混凝土骨料的網(wǎng)格劃分。
建立不同骨料部件的裝配。
或是通過AutoCAD軟件對骨料進(jìn)行處理后導(dǎo)出卵石狀骨料模型，并導(dǎo)入到ABAQUS內(nèi)。
計(jì)算特點(diǎn):
- 隱式分析: 核心是求解大型稀疏線性方程組。計(jì)算量集中在矩陣的分解和迭代求解上，對內(nèi)存容量、內(nèi)存帶寬和CPU的單核性能（頻率和緩存）都比較敏感。
- 顯式分析: 核心是時(shí)間步進(jìn)。為了保證計(jì)算穩(wěn)定，時(shí)間步長極小，導(dǎo)致總計(jì)算步數(shù)巨大。但每一步中，每個(gè)單元的計(jì)算相對獨(dú)立，是典型的“ embarrassingly parallel”（高度并行）問題。
詳細(xì)的化學(xué)反應(yīng)機(jī)理可能包含數(shù)千個(gè)反應(yīng)，在每個(gè)網(wǎng)格單元、每個(gè)時(shí)間步都需要計(jì)算。強(qiáng)耦合性: 流場、溫度場、化學(xué)組分場相互影響，求解過程復(fù)雜且收斂困難。時(shí)間步長小: 為捕捉火焰鋒面，需要極小的時(shí)間步長，導(dǎo)致總計(jì)算步數(shù)巨大。
-計(jì)算平臺(tái):
CPU多核計(jì)算(傳統(tǒng)基石): 傳統(tǒng)上，這類問題運(yùn)行在大型CPU計(jì)算集群上，通過MPI并行。GPU計(jì)算(前沿方向): GPU為燃燒仿真帶來了革命性變化。
集合的定義
單元集合的定義是為了將屬于同一層的所有單元歸為一組，以便對這一組單元統(tǒng)一指定材料方向和截面屬性。
由于單元是按層依次生成的，每層的單元編號(hào)是連續(xù)的，因此可以通過單元編號(hào)的范圍來定義每層的單元集合。
ABAQUS支持使用generate關(guān)鍵字來定義連續(xù)編號(hào)的單元集合，只需要指定起始編號(hào)、終止編號(hào)和步長即可。
對于N層的層合板，會(huì)生成N個(gè)單元集合，每個(gè)集合對應(yīng)一層材料。
