3D單元
關(guān)注創(chuàng)建者：CAE從業(yè)者 創(chuàng)建時間：2016-11-25
3D單元的視頻教程
復(fù)合材料3D實體單元建模
斯姆勒數(shù)值仿真研究院×技術(shù)鄰 復(fù)合材料3D實體單元建模，專業(yè)CAE數(shù)值仿真技術(shù)講解
免費 37分鐘 391播放
查看
ABAQUS復(fù)合材料3D殼單元建模
斯姆勒數(shù)值仿真研究院：從幾何尺寸、材料特性的問題描述來講解ABAQUS復(fù)合材料3D殼單元建模，詳細講解纖維增強復(fù)合材料的殼單元建模方法。
免費 22分鐘 1960播放
查看
3D單元的實例教程
前面文章主要講解了2d梁單元與2d實體單元的剛接問題，今日主要講解3d梁單元與殼單元的剛接問題。前面文章有講，梁單元除ROtZ外與殼單元有5個自由度物理意義相同，因而，當需要考慮梁單元與殼單元的剛接問題時，只需考慮該自由度與殼單元其他自由度的約束方程。具體處理方式可根據(jù)實際情況采用不同的處理方法。
3d梁單元與殼單元剛性連接按照位置關(guān)系的不同，可分為三類：
1）梁單元以一定角度與殼單元相交。
2）梁單元包含在殼單元內(nèi)。
3）梁單元在殼面但不包含的情況。
下面對這三類情況分別進行闡述。
一、梁單元以一定角度與殼單元相交
該類情況示意圖如下：
此種方式可以通過梁單元自由度ROTZ與殼單元其他自由度之間的數(shù)學關(guān)系，建立約束方程，命令流為CE。很多教材上面都有梁單元垂直穿過殼單元的經(jīng)典案例，例如一個典型的采用CE建立約束方程的命令流如下：
CE,1,0,142,ux,1,23,ux,-1,2,rotz,ny(142)-ny(58)
CE,2,0,92,uy,1,30,uy,-1,2,rotz,-(nx(92)-nx(30))
此種處理方式水哥個人不推薦，一則建立數(shù)學關(guān)系太煩，二則在一個工程中，如果此類情況較多，工作量實在太大，所以一般這類情況我們是通過建立剛性區(qū)域解決，這種方式會自動生成約束方程，雖不如之間建立約束方程合理，但能大大減少工作量。
關(guān)于剛性區(qū)域的討論，下期文章講解。
二、梁單元包含在殼面內(nèi)的情況
此種情況應(yīng)該是工程中最多的情況，典型情況便是我們經(jīng)常所見樓板與梁的連接。
展開 有限元基礎(chǔ)編程（終結(jié)篇）——C3D8單元程序編制
本篇推文是有限元基礎(chǔ)編程的終結(jié)篇，講述C3D8單元的程序編制及實現(xiàn)。主要內(nèi)容有：C3D8單元理論基礎(chǔ)、便于編程的“乘大數(shù)法”處理邊界條件、編制程序注意事項、云圖繪制函數(shù)、INP文件讀取函數(shù)、Abaqus仿真對比等，內(nèi)容量大，慢慢食用~
特別聲明：程序框架采用了吉林大學左文杰老師的腳本文件，計算單元剛度的核心計算程序仍延續(xù)我們以往編制程序的風格。代碼文件獲取方式詳見文末。
理論基礎(chǔ)
與Q4單元理論基礎(chǔ)相同，唯一的區(qū)別就是：每個節(jié)點的自由度由2變成了3，代碼具體變化看Ke函數(shù)和C3D8_cal_B函數(shù)的變化，理論部分可參考有限元基礎(chǔ)編程——Q4單元。
展開 我們今天介紹的通過python編程操作ABAQUS輸入文件生成PD3D單元顆粒，其可操作性更強，我們可以不采用粒子生成器內(nèi)部定義的隨機算法生成顆粒，用戶可以根據(jù)需求自定義顆粒分布算法，以契合實際工況。此外，可省去粒子生成顆粒的分析步，直接進行工況建模求解計算。
本貼只是個人興趣，只提供思路，不提供源碼，用戶需了解ABAQUS的inp文件的書寫規(guī)則、python操作文件語法和生成顆粒的底層邏輯(分布模型)。感興趣的可以私信，提供編寫思路。
下面我們采用這一方法生成直徑2mm、3mm、4mm和5mm的混合顆粒，數(shù)量為1000。具體生成結(jié)果如下圖所示。
分別用C3D8R單元和CPS4R單元進行模擬，當關(guān)閉大變形時，得出的兩條曲線完全重合；打開大變形后，兩條曲線差別很明顯，應(yīng)變越大，差距越大。請問版主，這是什么原因呀？
abaqus在處理幾何非線性的時候，兩種單元的處理方式有什么區(qū)別呢？
請版主幫我想想辦法，謝謝。
大家有沒有誰有3D奇異單元劃分網(wǎng)格的算例？簡單的3維結(jié)構(gòu)劃分的例子也行
3D單元的相關(guān)專題、標簽、搜索
3D單元的最新內(nèi)容
初始模型如下：
在step中使用熱力耦合分析步，在子程序中引入溫度相關(guān)的變形梯度
邊界條件設(shè)置：初始溫度場293K，同時設(shè)定Y+方向為393K，所有熱相關(guān)參數(shù)均使用文章的相關(guān)參數(shù)，左側(cè)固定，右側(cè)施加位移邊界條件，并使用C3D8T單元進行網(wǎng)格離散。
所有實體層采用 C3D8R 減縮積分單元并激活單元刪除，內(nèi)聚力層采用 COH3D8 單元，沖頭則使用離散剛體單元 R3D4。網(wǎng)格劃分基于掃掠技術(shù)（Advancing Front）生成。</p><p class="ql-align-justify"> 鋪層邏輯：支持非對稱鋪層序列輸入，用戶通過逗號分隔輸入各層角度。
3D陣列：天線單元呈立體排列，能夠在任何方向上控制一個或多個波束。
波束成形器的類型
無源電子掃描陣列（PESA）：無源相控陣列，是整個陣列只有一個收發(fā)器的天線。這是最常見的相控陣列配置類型。
有源電子掃描陣列（AESA）：在有源相控陣列天線中，每個天線單元或單元子集都有一個模擬收發(fā)器模塊，用于在每個單元中產(chǎn)生相移。這種更先進的方法通常用于軍事應(yīng)用。
考慮GND的大變形冷軋模擬4個月前
初始的計算模型如下所示：
軋制模型：
變形量為20%，整體包含500個晶粒，使用10萬C3D8R單元，整體計算時間為：34小時48分
變形后的結(jié)果如下圖所示
等效應(yīng)力分布：
等效塑性應(yīng)變分布：
幾何必須位錯密度分布：
統(tǒng)計儲存位錯密度分布：
可以看到和作者類似 的模擬趨勢，即GND分布于晶界相關(guān)，SSD分布主要是板材邊緣位置，
因科研需要，一直在研究一些單元算法，看著網(wǎng)上相關(guān)資料很多，但是和商軟對標的非線性單元技術(shù)相對較少。非線性這方面ABAQUS比較受人認可，所以打算用空余時間研究一下ABAQUS的單元技術(shù)，推導(dǎo)編寫一下相關(guān)程序供大家討論。本人水平十分有限，主要是學習ABAQUS的文檔，力學理論和代碼方面的問題請大家不吝賜教。
本文主要推導(dǎo)ABAQUS在幾何非線性（大變形）有限元分析中，用于計算單元切線剛度矩陣的算法
https://www.yqgqt.org.cn/content/post/403932
第七篇：C3D8六面體單元的剛度矩陣。
https://www.yqgqt.org.cn/content/post/430177
第八篇：UMAT用戶子程序開發(fā)步驟。
進行混凝土細觀三維模型的網(wǎng)格劃分，本案例中采用二次四面體單元（C3D10）,單元數(shù)量為362萬個。
在使用插件建立混凝土細觀三維幾何模型時，應(yīng)合理設(shè)置建模參數(shù)，以確保網(wǎng)格劃分的質(zhì)量。
基體選用金屬矩形板，彈性模量為70GPa，泊松比為0.3，尺寸為1×0.2×0.02(m)，選擇進行C3D8R單元進行網(wǎng)格劃分；壓電片材料選用PZT-5，采用上述壓電本構(gòu)模型，尺寸為0.1×0.1×0.01(m)。
3.邊界條件設(shè)置
邊界條件為基體板左側(cè)固定端約束，右端自由，壓電片上下表面施加5個周期的220V正弦交流電，如圖3所示。
這八個節(jié)點分為上下兩層，每層四個節(jié)點，節(jié)點的排列順序必須符合ABAQUS對C3D8R單元類型的要求。單元的生成采用了逐層處理的策略，對于每一層材料，將二維網(wǎng)格信息復(fù)制并擴展為三維單元信息。具體做法是，對于第i層，其下表面的四個節(jié)點來自第i層的二維網(wǎng)格節(jié)點，而上表面的四個節(jié)點來自第i加1層的二維網(wǎng)格節(jié)點。通過節(jié)點編號的系統(tǒng)性偏移，可以高效地生成所有層的單元連接關(guān)系。
