三維無限單元的案例
Abaqus無限單元的建立方法
什么是無限單元
無限單元是有限元的一種單元，它的單元形式如下圖所示。
2. 無限單元的作用
使用無限單元作為反射邊界，應力波將無反射，防止在邊界上產生的應力波反射重新進入模型，從而確保結果的正確性。
3. 怎么建立無限單元
（1）先建立有限元模型，然后將需要設置為無限單元的部分分割，在mesh模塊中設置該部分單元類型，這樣在inp文件中需要設置為無限單元的部分就會集中在一起，方便修改其節點的編號。
（2）在inp文件里對無限單元進行單元屬性及編號的改變，需要注意無限單元的方向，在二維無限單元中，前兩個節點所組成直線中點A與后兩個節點所組成直線中點B，無限單元的方向就是A指向B的方向，如圖1，二維無限元的方向是朝下；在三維無限單元中，前四個節點所組成平面的中心點C與后四個節點所組成平面的中心點D，無限單元的方向就是C指向D的方向，如圖1，三維無限單元的方向朝右。
（3）將修改后的inp導入，建立job提交就可以了。
實例講解-鋼丸撞擊金屬板
機械噴丸的模擬，其中設置金屬板邊界部分為無限單元。
（1）建立有限元模型，將金屬板的邊界partition切割，全部設置為C3D8R單元類型，修改邊界部分為C3D8I（邊界部分的單元類型和其他部分不一樣，方便下步統一更改為無限單元），如圖，建立job-write input 生成inp文件。
（2）修改inp文件的單元屬性及編號
inp文件以記事本的方式打開，將Element, type（單元類型）C3D8I修改為CIN3D8（無限單元），
下一步修改無限單元節點編號的排列順序，在inp的data lines中第一列是單元編號，后面的是節點編號，如圖7。在本文中，無限單元的方向是從金屬板的中心向外。
展開 基于LS-DYNA的柱狀藥包在無限水域中爆炸動態響應三維模擬 ￥12.9
為了更好地模擬柱狀藥包在水域中爆炸后沖擊波的傳播過程，模型采用ALE（任意拉格朗日歐拉算法），為了使模擬達到無線水域的效果，在模型邊界處施加無反射邊界條件，有限元模型及計算結果如下：
圖1 柱狀藥包在無限水域中爆炸動態響應有限元模型
圖2 柱狀藥包在無限水域中爆炸動態響應等效應力
圖3 柱狀藥包在無限水域中爆炸過程應力波傳播過程
本案例適用于研究爆炸、沖擊、侵徹動力學的朋友，下面附上該模擬的K文件，大家有疑問可以在私信我，歡迎交流！
展開 求含有板單元、梁單元、質量單元、彈簧單元、三維實體單元的實例
有哪位大神能給小弟提供上述實例哦，非常感謝！！！！
ABAQUS中無限單元的使用（一）——2D操作案例
ABAQUS中無限單元的使用（一）
——2D操作案例
作者：Nick_Liu
使用無限單元可以用消除應力波的反射，從而以較小的模型得到和整體模型一致的計算結果，在節省計算資源和提高計算精度方面具有很大優勢。一般局部加載對整體影響不大或者問題為無邊界時可以使用無限單元模型。
ABAQUS 6.14提供了17種無限單元。非聲學問題中常用的無限單元有：CINAX4、CIN3D8、CINPE4和CINPS4。使用時只需在input文件中修改單元類型的關鍵字即可。
使用無限單元的關鍵在于模型區域的劃分。本文通過隨時空變化壓力載荷作用于圓柱體的案例，演示二維無限單元模型的建立過程。
1）建模過程
Step-1：建立軸對稱模型，并將模型分成3個區域，如圖1所示：
圖1 建模和分塊
Step-2：設置邊界區域厚度方向的種子密度為1，寬度方向種子數量相同；內部區域單元類型為CAX4R，邊界區域單元類型為CAX4；劃分網格如圖2所示：
圖2 模型網格
Step-3：孤立網格，調整堆疊方向，如圖3所示：
圖3 調整堆疊方向
Step-4：賦予材料屬性，建立裝配體和動態分析步（Dynamic, Explicit）。
Step-5：在模型頂部通過VDLoad子程序施加大小隨時間和空間變化的壓力；同時對模型施加對稱和固定約束，如圖4所示：
圖4 載荷和約束
Step-5：新建Job，輸出Input文件，將Input文件中的單元類型CAX4改為CINAX4。
Step-6：通過修改后的Input文件新建Job并提交計算。
2）計算結果
計算結果如圖5所示：
圖5 計算結果
可以看到，應力傳遞至白色的無限單元邊界時沒有進行反射，其結果與實際情況是一致的。
3）注意事項
a. 厚度方向只能有一個單元；
b.
展開 
位移邊界條件：三維實體單元與梁單元的區別與聯系
在有限元中，三維實體單元的節點只有三個平移自由度；而梁單元，或者殼單元等還存在旋轉自由度。筆者當初留意到這點的時候，由于才疏學淺，還曾糾結了很長時間，筆者之前也寫過這方面的見解。那么當在實體單元下的約束關系和梁單元下的約束關系，有什么樣的區別和聯系呢。
梁構件 （ 梁單元）
建立一個矩形梁，長度為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)
建立一個矩形梁（實體）
長度為2米，截面是0.04m*0.06m。
對比（自由模態）
求解兩種單元下的自由模態。求解結果基本一致。
對比（固定約束模態）
梁單元模型約束一個端點六個自由度，實體單元模型約束一個端面的三個自由度。求解結果基本一致。
展開 Abaqus中選擇三維實體單元類型的基本原則 附abaqus三維筒體過渡網格劃分下載
來源：力學與Abaqus仿真
對于大多數Abaqus用戶，在選擇單元類型時都會有這樣的困惑，可選的單元類型很多，還有減縮積分、完全積分、線性單元、二次單元、非協調單元、雜交單元、沙漏控制等眾多選擇（圖1），在實際有限元分析時，究竟應該如何選擇合適的單元類型。從今天開始，陸續介紹單元類型的選取原則，供大家參考。
圖1 單元類型選擇對話框
選擇三維實體單元類型時應遵循以下原則：
● 對于三維區域，盡可能采用結構化網格劃分技術或掃掠網格劃分技術，從而得到Hex單元網格，減小計算代價，提高計算精度。當幾何形狀復雜時，也可以在不重要的區域使用少量楔形（Wedge）單元。
● 如果使用了自由網格劃分技術，Tet單元的類型應選擇二次單元。在Abaqus/Explicit中應選擇修正的Tet單元 C3D10M，在Abaqus/Standard中可以選擇C3D10，但如果有大的塑性變形，或模型中存在接觸，而且使用的是默認的“硬”接觸關系（“hard”contact relationship），則也應選擇修正的Tet單元 C3D10M。
● Abaqus的所有單元均可用于動態分析，選取單元的一般原則與靜力分析相同。但在使用Abaqus/Explicit模擬沖擊或爆炸載荷時，應選用線性單元，因為它們具有集中質量公式，模擬應力波的效果優于二次單元所采用的一致質量公式。
如果使用的求解器是Abaqus/Standard，在選擇單元類型時還應注意以下方面：
● 對于應力集中問題，盡量不要使用線性減縮積分單元，可使用二次單元來提高精度。如果在應力集中部位進行了網格細化，使用二次減縮積分單元與二次完全積分單元得到的應力結果相差不大，而二次減縮積分單元的計算時間相對較短。
展開 abaqus2020-三維-顯示分析-通用接觸或接觸對接觸-單元刪除法模擬裂紋，單元穿透問題！！
1 abaqus2020-三維-顯示分析-單元刪除法模擬裂紋，僅采用通用接觸時，模型中出現明顯穿透，結果不合理！
2 abaqus2020-三維-顯示分析-單元刪除法模擬裂紋，僅采用接觸對接觸時，模型中出現少許穿透，結果相對合理，但不是最理想狀態！
3 abaqus2020-三維-顯示分析-單元刪除法模擬裂紋，同時采用通用接觸+接觸對接觸時，模型中無明顯穿透，結果合理！
【JY】Abaqus 三維應力單元解析、選擇與應用指南
2.2 不同求解器下的單元選擇差異
Abaqus/Standard 和 Abaqus/Explicit 在單元選擇上存在一定差異：
Abaqus/Standard：實體單元多樣，有二維和三維的線性、二次單元，積分方式可選完全積分或縮減積分。還包含修正的二次 Tri、Tet 單元，以及非協調模式單元和雜交單元。
Abaqus/Explicit：實體單元包括二維和三維的線性縮減積分單元，以及修正的二次 Tri 和 Tet 單元，但不存在二次完全積分實體單元。
2.3 特定場景下的單元選擇
針對不同的工程場景，三維實體單元的選擇策略也有所不同：
應力集中問題：對于應力集中問題，盡量不要使用線性縮減積分單元，可使用二次單元來提高精度。如果在應力集中部位進行了網格細化，使用二次縮減積分單元與二次完全積分單元得到的應力結果相差不大，而二次縮減積分單元的計算時間相對較短。
彈塑性分析：對于彈塑性分析，如果材料是不可壓縮性的（例如金屬材料），則不能使用二次完全積分單元，否則會出現體積自鎖問題，也不要使用二次 Tri 單元或 Tet 單元。推薦使用的是修正的二次 Tri 單元或 Tet 單元、非協調單元，以及線性縮減積分單元。如果使用二次縮減積分單元，當應變超過 20%~40% 時要劃分足夠密的網格。
接觸分析或大扭曲變形：如果模型中存在接觸或大的扭曲變形，則應使用線性四邊形 (Quad) 或六面體 (Hex) 單元，以及修正的二次 Tri 單元或 Tet 單元，而不能使用其他的二次單元。
彎曲為主的問題：對于以彎曲為主的問題，如果能夠保證在所關心部位的單元扭曲較小，使用非協調單元（例如 C3D8I 單元）可以得到非常精確的結果。
展開 三維實體模型能不能用殼單元
為了容易裝配，模型全部建的實體模型，請問還能運用殼單元嗎？
如果想使用該如何操作？
Matlab給ABAQUS模型插入粘聚力單元（二維和三維）
前言：以常用的幾種實體單元為例，展示了在單一單元、兩種單元和三種單元混合的模型中插入黏聚力單元后的網格圖。二維模型中插入的黏聚力單元為coh2d4單元，三維模型中插入的黏聚力單元為coh3d6、coh3d8單元。
三維梁單元模型中，想輸出梁柱截面的內力
在三維梁單元模型中，想輸出梁柱截面的內力，比方說彎矩、剪力、軸力，好像說這些內力都是局部坐標系中的內力，但是在sf輸出中，有sf1,sf2,sf3,sm1,sm2,sm3
SF1 Axial force.
SF2 Transverse shear force in the local 2-direction (not available for B23, B23H, B33, B33H).
SF3 Transverse shear force in the local 1-direction (available only for beams in space, not available for B33, B33H).
SM1 Bending moment about the local 1-axis.
SM2 Bending moment about the local 2-axis (available only for beams in space).
SM3 Twisting moment about the beam axis (available only for beams in space).
梁單元中，sm3是繞梁軸線的彎矩，sm1是繞梁1方向的彎矩。sf1是3方向即梁軸線方向的力（軸力），sf3是梁1方向的力
是局部坐標，空間梁這六個量都可能有值，SF1是軸力，SF2、SF3是剪力，SM1、SM2是彎矩，SM3是扭矩。
展開 
應用三維有限單元法計算應力強度因子
應用三維有限單元法計算應力強度因子
來源：中國機械工程 作者：林曉斌
摘要 描述了兩種基于有限單元計算面形裂紋應力強度因子的方法，建議了一種創造三維有限單元網格的途徑。計算方法的精度通過和其它解析解或數值解的比較得到了說明。
關鍵詞 應力強度因子 有限元 損傷容限設計 斷裂評定
無論在損傷容限設計還是在缺陷評定階段，工程師們需要知道正在分析的構件中裂紋的應力強度因子，因為判斷含裂紋構件的斷裂，或者計算剩余疲勞壽命大多依賴于這一參量。因此，在斷裂力學發展中，如何求取應力強度因子一直是一個重要的課題。當前已有許多方法可用來計算應力強度因子，較為典型的有解析法、邊界配位法、有限單元法、邊界元素法、體力法、權函數法和線彈簧模型。利用這些方法，大量的應力強度因子解已經獲得，已出版的應力強度手冊［1］中收編了許多典型的解。盡管如此，工程師們仍然會感到自己所需要的應力強度因子解很難找到，這是因為要解決的工程問題往往是一些受復雜載荷的構件，包含的裂紋也往往是一些不規則裂紋。
本文簡單介紹了兩種基于三維有限單元法計算面形裂紋應力強度因子的方法。有限單元法已經成為工程設計分析領域中一個強有力的計算工具，它能模擬非常復雜的構件。基于有限元的應力強度因子計算方法，自然也將具有卓越的工程能力。除了計算方法的介紹以外，還將簡單描述一種簡化網格的生成方法。最后提供了一些所得到的典型應力強度因子解，并和大家熟知的解進行了比較，以說明本文所描述的方法的可靠性。
展開 基于cohesive的三維脆性材料斷裂（abaqus cohesive單元本構） ￥999
1、根據論文《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 編寫的cohesive單元本構
2、適用于三維模型
3、包含umat以及vumat
4、umat適用范圍小，多個cohesive單元一般采用vumat進行計算
5、軟化曲線為Hordijk和bilinear
展開 ANSYS提取具體三維單元的體積，面元的面積和線的長度
在ANSYS中，能提取具體三維單元的體積，面元的面積和線的長度。
如：*GET,E_VOLUME,ELEM,10,VOLU 為提取編號為10的單元的體積
*GET,a_area,AREA,50,AREA 為提取編號為50的面元的面積
*GET,l_length,LINE,100,LENG 為提取編號為100的線的長度
以上對應的GUI操作： Utility Menu>Parameters>Get Scalar Data 如果要一次性提取多個元素的相關參數，可以用命令 *VGET, ParR, Entity, ENTNUM, Item1, IT1NUM, Item2, IT2NUM, KLOOP 對應GUI操作：Utility Menu>Parameters>Get Array Data
輸入命令 alist,p 出選項框，選你要看的那個面積，提取選中的單元面積。
*cfopen,'area','txt',
*GET,MaxEleNum,ELEM,,NUM,MAX
*GET,MinEleNum,ELEM,,NUM,MIN
*do,i,MinEleNum,MaxEleNum,1
*if,esel(i),eq,1,then
*get,volu,elem,i,volu
*vwrite,i,volu
(f5.0,f15.12)
*end if
*enddo
*cfcols
展開 cohesive單元分層開膠斷裂模擬-雙懸臂梁剝離DCB（三維模型） ￥2
cohesive單元分層開膠斷裂模擬-雙懸臂梁剝離DCB（三維模型）
