abaqus 三維單元的案例
Abaqus中選擇三維實(shí)體單元類型的基本原則 附abaqus三維筒體過(guò)渡網(wǎng)格劃分下載
來(lái)源：力學(xué)與Abaqus仿真
對(duì)于大多數(shù)Abaqus用戶，在選擇單元類型時(shí)都會(huì)有這樣的困惑，可選的單元類型很多，還有減縮積分、完全積分、線性單元、二次單元、非協(xié)調(diào)單元、雜交單元、沙漏控制等眾多選擇（圖1），在實(shí)際有限元分析時(shí)，究竟應(yīng)該如何選擇合適的單元類型。從今天開(kāi)始，陸續(xù)介紹單元類型的選取原則，供大家參考。
圖1 單元類型選擇對(duì)話框
選擇三維實(shí)體單元類型時(shí)應(yīng)遵循以下原則：
● 對(duì)于三維區(qū)域，盡可能采用結(jié)構(gòu)化網(wǎng)格劃分技術(shù)或掃掠網(wǎng)格劃分技術(shù)，從而得到Hex單元網(wǎng)格，減小計(jì)算代價(jià)，提高計(jì)算精度。當(dāng)幾何形狀復(fù)雜時(shí)，也可以在不重要的區(qū)域使用少量楔形（Wedge）單元。
● 如果使用了自由網(wǎng)格劃分技術(shù)，Tet單元的類型應(yīng)選擇二次單元。在Abaqus/Explicit中應(yīng)選擇修正的Tet單元 C3D10M，在Abaqus/Standard中可以選擇C3D10，但如果有大的塑性變形，或模型中存在接觸，而且使用的是默認(rèn)的“硬”接觸關(guān)系（“hard”contact relationship），則也應(yīng)選擇修正的Tet單元 C3D10M。
● Abaqus的所有單元均可用于動(dòng)態(tài)分析，選取單元的一般原則與靜力分析相同。但在使用Abaqus/Explicit模擬沖擊或爆炸載荷時(shí)，應(yīng)選用線性單元，因?yàn)樗鼈兙哂屑匈|(zhì)量公式，模擬應(yīng)力波的效果優(yōu)于二次單元所采用的一致質(zhì)量公式。
如果使用的求解器是Abaqus/Standard，在選擇單元類型時(shí)還應(yīng)注意以下方面：
● 對(duì)于應(yīng)力集中問(wèn)題，盡量不要使用線性減縮積分單元，可使用二次單元來(lái)提高精度。如果在應(yīng)力集中部位進(jìn)行了網(wǎng)格細(xì)化，使用二次減縮積分單元與二次完全積分單元得到的應(yīng)力結(jié)果相差不大，而二次減縮積分單元的計(jì)算時(shí)間相對(duì)較短。
展開(kāi) abaqus梁單元的熱分析？
各位大佬，請(qǐng)問(wèn)abaqus里三維梁單元無(wú)法進(jìn)行heat transfer分析嗎，在mesh類型里，選擇heat transfer后單元類型的后面沒(méi)有增加T，計(jì)算的時(shí)候也會(huì)報(bào)錯(cuò)，顯示單元缺少屬性定義。
abaqus2020-三維-顯示分析-通用接觸或接觸對(duì)接觸-單元刪除法模擬裂紋，單元穿透問(wèn)題！！
1 abaqus2020-三維-顯示分析-單元刪除法模擬裂紋，僅采用通用接觸時(shí)，模型中出現(xiàn)明顯穿透，結(jié)果不合理！
2 abaqus2020-三維-顯示分析-單元刪除法模擬裂紋，僅采用接觸對(duì)接觸時(shí)，模型中出現(xiàn)少許穿透，結(jié)果相對(duì)合理，但不是最理想狀態(tài)！
3 abaqus2020-三維-顯示分析-單元刪除法模擬裂紋，同時(shí)采用通用接觸+接觸對(duì)接觸時(shí)，模型中無(wú)明顯穿透，結(jié)果合理！
【JY】Abaqus 三維應(yīng)力單元解析、選擇與應(yīng)用指南
筆者近期在整理相關(guān)研究資料時(shí)，系統(tǒng)梳理了 Abaqus 中實(shí)體單元的分類邏輯、理論基礎(chǔ)及不同場(chǎng)景下的選擇策略，發(fā)現(xiàn)現(xiàn)有實(shí)踐中有粉絲仍存在單元類型誤用、特性理解不充分等問(wèn)題。鑒于此，本文將從單元分類、選擇原則、特定場(chǎng)景應(yīng)用及最佳實(shí)踐等方面展開(kāi)論述，旨在為從事 Abaqus 仿真分析的研究者與工程技術(shù)人員提供系統(tǒng)性參考，助力提升數(shù)值模擬的可靠性與科學(xué)性。
在現(xiàn)代工程仿真領(lǐng)域，Abaqus 作為一款功能強(qiáng)大的有限元分析軟件，提供了豐富的單元庫(kù)來(lái)滿足各種結(jié)構(gòu)力學(xué)分析需求。特別是在三維應(yīng)力分析場(chǎng)景中，正確選擇和使用單元類型對(duì)于獲得準(zhǔn)確、高效的計(jì)算結(jié)果至關(guān)重要。Abaqus 單元庫(kù)按照單元族分類，主要包括實(shí)體單元 (C)、殼單元 (S)、梁單元 (B)、桁架單元 (T)、剛體單元 (R) 等。每個(gè)單元族又包含多種具體單元類型，適用于不同的幾何特征、載荷條件和分析目標(biāo)。
本篇是三維應(yīng)力單元，即實(shí)體單元篇！
【相關(guān)閱讀】
【JY】有限元分析的單元類型分享一波~
【JY】有限單元分析的常見(jiàn)問(wèn)題及單元選擇
1 實(shí)體單元分類與理論基礎(chǔ)
實(shí)體單元是 Abaqus 中最基礎(chǔ)也最常用的單元類型，可在其任何表面與其他單元連接，能夠精確地反映幾何形狀的復(fù)雜性，適用于模擬具有復(fù)雜形狀的結(jié)構(gòu)。在 Abaqus 中，實(shí)體單元根據(jù)節(jié)點(diǎn)位移插值階數(shù)、積分方式和特殊功能可分為多種類型。
1.1 按節(jié)點(diǎn)位移插值階數(shù)分類
根據(jù)節(jié)點(diǎn)位移插值階數(shù)，實(shí)體單元主要分為三類：
線性單元（一階單元）：節(jié)點(diǎn)僅布置在單元角點(diǎn)，各方向采用線性插值，計(jì)算相對(duì)簡(jiǎn)單，適用于變形較簡(jiǎn)單的情況。線性單元的主要優(yōu)勢(shì)在于計(jì)算效率高，但精度相對(duì)較低，適用于對(duì)精度要求不高或初步分析階段。
展開(kāi) 
Matlab給ABAQUS模型插入粘聚力單元（二維和三維）
前言：以常用的幾種實(shí)體單元為例，展示了在單一單元、兩種單元和三種單元混合的模型中插入黏聚力單元后的網(wǎng)格圖。二維模型中插入的黏聚力單元為coh2d4單元，三維模型中插入的黏聚力單元為coh3d6、coh3d8單元。
基于cohesive的三維脆性材料斷裂（abaqus cohesive單元本構(gòu)） ￥999
1、根據(jù)論文《Three-dimensional modeling of fracture in quasi-brittle materials using plasticity and cohesive finite elements》DOI:https://doi.org/10.1007/s10704-021-00514-1 編寫(xiě)的cohesive單元本構(gòu)
2、適用于三維模型
3、包含umat以及vumat
4、umat適用范圍小，多個(gè)cohesive單元一般采用vumat進(jìn)行計(jì)算
5、軟化曲線為Hordijk和bilinear
展開(kāi) 采用ABAQUS連接單元等效建筑隔震支座，實(shí)現(xiàn)二維、三維隔震分析。
采用ABAQUS連接單元等效建筑基礎(chǔ)隔震支座，實(shí)現(xiàn)結(jié)構(gòu)二維、三維隔震分析。水平自由度可實(shí)現(xiàn)雙線性恢復(fù)力模型等，豎向自由度可實(shí)現(xiàn)彈簧恢復(fù)力模型、具有耗能能力的摩擦彈簧恢復(fù)力模型等。
求含有板單元、梁單元、質(zhì)量單元、彈簧單元、三維實(shí)體單元的實(shí)例
有哪位大神能給小弟提供上述實(shí)例哦，非常感謝！！！！
ABAQUS隔震分析—隔震單元雙線性模型的參數(shù)確定及設(shè)置方法（可實(shí)現(xiàn)三維隔震） ￥68
ABAQUS隔震分析—隔震單元雙線性模型的參數(shù)確定及設(shè)置方法
位移邊界條件：三維實(shí)體單元與梁單元的區(qū)別與聯(lián)系
在有限元中，三維實(shí)體單元的節(jié)點(diǎn)只有三個(gè)平移自由度；而梁單元，或者殼單元等還存在旋轉(zhuǎn)自由度。筆者當(dāng)初留意到這點(diǎn)的時(shí)候，由于才疏學(xué)淺，還曾糾結(jié)了很長(zhǎng)時(shí)間，筆者之前也寫(xiě)過(guò)這方面的見(jiàn)解。那么當(dāng)在實(shí)體單元下的約束關(guān)系和梁單元下的約束關(guān)系，有什么樣的區(qū)別和聯(lián)系呢。
梁構(gòu)件 （ 梁單元）
建立一個(gè)矩形梁，長(zhǎng)度為2米，截面是0.04m*0.06m。
FINISH$/CLEAR
/FILNAME,BEARIGN
/TITLE,COMBIN14 AND BEAM188
!UNITS,S-M-KG-N
/PREP7
ET,1,BEAM188$ET,2,COMBIN14,,,0
MP,EX,1,200E9$MP,PRXY,1,0.3$MP,DENS,1,7850
R,2,2E8$SECTYPE,1,BEAM,RECT$SECDATA,0.04,0.06
N,1,0,0$N,2,0.01,0.01,0$N,3,0.01,0,0.01
K,1,0.01,0$K,2,1,0$L,1,2
LATT,1,,1,,,,1$LESIZE,ALL,,,20
LMESH,ALL$BN=NODE(0.01,0,0)
建立一個(gè)矩形梁（實(shí)體）
長(zhǎng)度為2米，截面是0.04m*0.06m。
對(duì)比（自由模態(tài)）
求解兩種單元下的自由模態(tài)。求解結(jié)果基本一致。
對(duì)比（固定約束模態(tài)）
梁單元模型約束一個(gè)端點(diǎn)六個(gè)自由度，實(shí)體單元模型約束一個(gè)端面的三個(gè)自由度。求解結(jié)果基本一致。
展開(kāi) Abaqus單元的命名規(guī)則
概述
Abaqus/Standard殼單元庫(kù)包括：
? 用于三維殼體幾何的單元；
? 具有軸對(duì)稱變形的軸對(duì)稱幾何的單元；
? 具有關(guān)于一個(gè)平面對(duì)稱的一般變形的軸對(duì)稱幾何的單元；
? 用于應(yīng)力/位移、傳熱和完全耦合的溫度位移分析的單元；
? 通用單元，以及專門適用于分析“厚”或“薄”殼的單元；
? 使用簡(jiǎn)化或減縮積分的通用、三維、一階單元；
? 考慮有限膜應(yīng)變的單元；
? 在可能的情況下，每個(gè)節(jié)點(diǎn)使用五個(gè)自由度的單元，以及每個(gè)節(jié)點(diǎn)始終使用六個(gè)自由度的單元；
? 連續(xù)體殼單元。
Abaqus/Explicit 殼單元庫(kù)包括：
? 通用三維單元，用于模擬考慮有限膜應(yīng)變的“厚”或“薄”殼；
? 小應(yīng)變單元；
? 完全耦合的溫度位移分析單元；
? 用于具有軸對(duì)稱變形的軸對(duì)稱幾何形狀的單元；
? 連續(xù)體殼單元。
命名約定
殼單元的命名約定取決于單元的維度。
三維殼單元
Abaqus中的三維殼單元命名如下：
第一項(xiàng)：
常規(guī)應(yīng)力/位移單元（S）
連續(xù)應(yīng)力/位移單元（SC）
三角形應(yīng)力/位移薄殼（STRI）
熱傳輸殼（DS）
第二項(xiàng)：節(jié)點(diǎn)數(shù)
第三項(xiàng)：減縮積分（可選）
第四項(xiàng)：
可選五自由度
耦合溫度-位移（T）
ABAQUS/Explicit中的小應(yīng)變公式（S）
第五項(xiàng)：
在ABAQUS/Explicit中的小應(yīng)變公式中考慮翹曲（可選）
例如，S4R是一個(gè)四節(jié)點(diǎn)四邊形應(yīng)力/位移殼單元，具有減縮積分和大應(yīng)變公式；SC8R是一個(gè)八節(jié)點(diǎn)、四邊形、一階插值、應(yīng)力/位移連續(xù)體殼單元，具有減縮積分。
展開(kāi) 
【JY】Abaqus“殼”單元概述與應(yīng)用（三）——非線性擬協(xié)調(diào)固體連續(xù)殼單元CSS8
【相關(guān)閱讀】
【JY】Abaqus殼單元概述與應(yīng)用（一）
【JY】Abaqus 三維應(yīng)力單元解析、選擇與應(yīng)用指南
【JY】Abaqus“殼”單元概述與應(yīng)用（二）——固體殼單元
傳統(tǒng)固體殼單元在處理幾何非線性、材料非線性及復(fù)雜邊界條件時(shí)，存在諸多難以克服的缺陷，這促使研究者探索新的單元構(gòu)造方法。非線性擬協(xié)調(diào)固體殼單元的提出，正是為了突破這些局限，其研究動(dòng)因主要源于以下幾方面：
（一）傳統(tǒng)固體單元的固有缺陷
自鎖現(xiàn)象普遍存在
傳統(tǒng)固體單元（如C3D8R）在模擬薄板殼結(jié)構(gòu)時(shí)，易出現(xiàn)剪切自鎖、薄膜自鎖、體積自鎖等問(wèn)題。剪切自鎖源于單元位移插值無(wú)法準(zhǔn)確表征純彎曲狀態(tài)下的零剪切應(yīng)變，導(dǎo)致計(jì)算結(jié)果剛度偏高；薄膜自鎖則因低階形函數(shù)無(wú)法捕捉不可伸縮彎曲模式下的面內(nèi)應(yīng)變分布，使位移被低估；體積自鎖多見(jiàn)于近不可壓縮材料分析，由于單元無(wú)法準(zhǔn)確描述等體積運(yùn)動(dòng)，導(dǎo)致體積變化被過(guò)度約束。這些自鎖現(xiàn)象嚴(yán)重影響計(jì)算精度，尤其是在粗網(wǎng)格或大長(zhǎng)高比結(jié)構(gòu)中表現(xiàn)更為突出。
計(jì)算效率與精度的矛盾
為克服自鎖問(wèn)題，需要采用增強(qiáng)假設(shè)應(yīng)變法（EAS）、假設(shè)自然應(yīng)變法（ANS）或雜交應(yīng)力法等，這些方法往往需要引入額外的內(nèi)部參數(shù)或復(fù)雜的數(shù)值積分，使得單元列式復(fù)雜、相對(duì)殼單元計(jì)算成本增加。
幾何非線性處理的局限性
現(xiàn)有非線性固體殼單元多基于連續(xù)體變形梯度的極分解處理幾何非線性，該方法不僅計(jì)算量大，且在 Cartesian 坐標(biāo)系下難以保證旋轉(zhuǎn)描述的準(zhǔn)確性。在大變形、大轉(zhuǎn)動(dòng)問(wèn)題中，極分解可能導(dǎo)致切線剛度矩陣奇異，影響迭代收斂性。此外，傳統(tǒng)單元在處理不規(guī)則網(wǎng)格或畸變網(wǎng)格（如C3D8I）時(shí)，精度衰減明顯，難以滿足工程對(duì)復(fù)雜結(jié)構(gòu)分析的需求。
展開(kāi) (一維彈簧單元的直接剛度方法)Python編程和ABAQUS結(jié)果對(duì)比
在學(xué)習(xí)《有限元方法基礎(chǔ)教程》過(guò)程中，通過(guò)自己編程實(shí)現(xiàn)有限單元法的數(shù)值解答，加深對(duì)理論的理解。這個(gè)在我去年的時(shí)候就已經(jīng)發(fā)過(guò)帖子了，但是沒(méi)有講解代碼，沒(méi)有和ABAQUS有限元軟件對(duì)比。這次發(fā)帖子就是想講解代碼以及和ABAQUS結(jié)果對(duì)比，發(fā)現(xiàn)我又更深層次理解了有限元方法。
這次寫(xiě)得是最簡(jiǎn)單的模型：一維線彈簧單元。采用的是直接剛度法求解。
直接剛度法的求解思路如下所示，其中粉紅色的是輸入，淡藍(lán)色是輸出。主要是要區(qū)分齊次邊界條件和非齊次邊界條件，非齊次邊界條件的話就要修改【F】。
下面將貼出我用python寫(xiě)得一維彈簧單元的直接剛度法：
例子計(jì)算：
如圖是一個(gè)彈簧系統(tǒng)，單元節(jié)點(diǎn)信息如下，5節(jié)點(diǎn)受到一個(gè)強(qiáng)制位移20mm，明顯這是一個(gè)非齊次邊界條件問(wèn)題。
Python編程輸入信息如下：
結(jié)果如下：
可以看到，輸出結(jié)果和書(shū)上的答案一致。
下面進(jìn)行ABAQUS模擬：
添加彈簧單元
添加邊界條件
顯示單元編號(hào)、節(jié)點(diǎn)編號(hào)如下所示，紅色表示單元編號(hào)，黃色是節(jié)點(diǎn)編號(hào)
ABAQUS計(jì)算結(jié)果如下：
首先是變形圖前后對(duì)比
反力云圖如下所示，基本和直接剛度法計(jì)算的結(jié)果一致
位移云圖如下所示，基本和直接剛度法計(jì)算的結(jié)果一致
整體剛度矩陣如下所示，因?yàn)?em>ABAQUS彈簧單元是三維的，每一個(gè)節(jié)點(diǎn)有3個(gè)自由度，15x15，原味的剛度矩陣如下
我們把剛度矩陣轉(zhuǎn)化為一維的，方便和編程的結(jié)果對(duì)比
從結(jié)果可以看出，ABAQUS的整體剛度矩陣和直接剛度法計(jì)算出來(lái)的整體剛度矩陣有些差異，如圖標(biāo)紅的所示。
那么在整體剛度矩陣上為什么ABAQUS會(huì)和直接剛度法的整體剛度矩陣有差異呢，到底ABAQUS的整體剛度矩陣對(duì)不對(duì)呢？答案將在下一期揭曉。歡迎大家積極討論。
展開(kāi) 三維實(shí)體模型能不能用殼單元
為了容易裝配，模型全部建的實(shí)體模型，請(qǐng)問(wèn)還能運(yùn)用殼單元嗎？
如果想使用該如何操作？
Abaqus 三維鉆孔仿真案例教學(xué) ￥29.99
</p><h1><strong style="color: rgb(255, 0, 0);">附件：完整案例教學(xué)內(nèi)容和本案例中的abaqus模型文件（包括cae、odb和inp文件）</strong></h1><h2><br></h2><p><br></p>
