ansys 四面體單元的案例
插件案例：任意單元面之間插入0厚度cohesive單元（支持四面體or六面體單元） ￥450
1.插件長什么樣子
2.插件有什么功能
在四面體或者六面體單元之間插入 0厚度的cohesive單元，
并對原基體材料和cohesive材料賦予材料屬性及網格屬性。
支持對整個part操作或進對part中的部分單元集合操作
3.插件效果
以下是六面體單元組成的球體，右側是在兩個單元面之間插入的cohesive效果圖
以下是某脆性材料沖擊破碎的效果圖
以下是在六面體之間插入COH3D8界面單元
（1） 無cohesive六面體單元 （2）插入0厚度cohesive單元
以下是在四面體單元之間插入COH3D6界面單元
（1） 無cohesive四面體單元 （2）插入0厚度cohesive單元
4.該插件的安裝
下載相應版本的插件程序，解壓縮到ABAQUS安裝目錄或者工作目錄下的abaqus_plugins文件夾內，工作目錄下默認是沒abaqus_plugins文件夾的，可以自己創建一個，然后重新啟動ABAQUS，在plug_ins菜單中就可以看到該工具。
以下是付費內容，主要包括
(1)該插件的編譯代碼，適合于abaqus6.14，請將你的Abaqus自行更新到此版本，
(2)附帶一個小球沖擊破碎的案例CAE文件+inp文件，模型如下圖所示
附件內容如下圖所示：
看清介紹，附件是執行代碼，不要盲目亂拍，源代碼價格另議，需要者站內聯系或添加QQ3539358512
展開 基于tcl語言創建四面體/五面體實體單元 ￥20
本案例是基于tcl語言實現用戶自定義的單元，并獲取單元的中心點，并依據單元中心點及單元節的最短邊去移動單元中心點，實現四面體/五面體實體單元的創建。具體實現過程見本案例的程序部分。
詳情見收費的程序部分，凡購買本案例的朋友針對該案例有疑問，可私信，謝謝！
ANSA中四面體單元質量的檢查與提高
本文將介紹在前處理軟件ANSA中，進行四面體網格劃分時，如何進行單元質量的檢查及提高。
1、單元質量的檢查
在ANSA界面中，點擊F11功能鍵彈出Presentation
Parameter窗口，如圖1所示，在這一窗口的Solid標簽中，用戶可以根據具體的要求進行質量標準的設定。
設定的內容主要包括：
然后，將模型顯示切換到HIDDEN顯示模式，違反質量準則的單元根據右側圖標的顏色顯示出來，如圖2所示。視圖中違反準則的單元數量顯示在窗口左側的OFF后面。鼠標放置在OFF上，按右鍵選擇show
only可將不合格的單元單獨顯示。
2、單元質量的提高
四面體網格的生成是基于三角形面網格的，所以，在提高四面體網格質量時，可采用兩種方法，一種是對體網格質量較差區域的面網格進行調整，然后再重新生成此處的體網格；另一種是直接采用體網格功能里Improve中的Reconstruct、Fix
Quality、Smooth、Move幾種功能來進行質量提高。
在使用Fix
Quality功能時，需設定網格調整時相應參數的變化范圍以使網格在允許的范圍內調整，如圖3所示。
調整后的網格模型如圖4所示，沒有off單元的存在，說明網格都符合要求，完成調整。
內容來源:有限元在線 版權歸作者所有
展開 案例實操：四面體單元懸臂梁的Matlab有限元編程過程講解
之前的課程我們學習了一維梁單元，二維平面單元，三維板殼單元的matlab有限元編程，本次案例主要講解如何用matlab實現針對四面體單元劃分的三維結構進行有限元編程，具體案例是一個懸臂梁受集中荷載的問題。圖1為本案例Matlab編程計算得到的結果。主要內容涉及四面體單元的有限元基本理論的推導，主要是單元剛度矩陣的推導，此外還包括等參單元和Hammer數值積分以及三維問題的后處理計算。
圖1 懸臂梁受集中荷載的應力云圖
一個完整的有限元程序基本組成部分包括前處理模塊、分析主程序模塊和后處理模塊。在前處理模塊中，實現節點坐標輸入、單元節點編號、網絡劃分以及邊界條件輸入等工作；在分析主程序模塊中，求解整體剛度方程；在后處理模塊中，實現結果顯示、數據輸出等工作。對應的有限元法的基本步驟：（1）幾何域離散，獲得標準化的單元；（2）通過能量原理（虛功原理或最小勢能原理，獲得單元剛度方程；（3）單元的集成（裝配）；（4）處理位移邊界條件；（5）計算位移場；（6）計算單元的其他物理量（應力應變）。這幾步中，最核心的內容是單元研究，具體包括：（1）節點描述（不同坐標系節點坐標的變化）；（2）場描述（位移場，應變場，應力場，形函數）；（3）單元剛度方程（基于能量原理推導）。需要說明的是后文的四面體單元有限元方程的推導過程是基于等參單元的基本理論從局部坐標（自然坐標、體積坐標）出發來推導四面體單元的剛度矩陣，因為這樣做比較規范自然，推導過程也適用于其他類型單元。但是因為四面體單元相對簡單也可以直接從直角坐標（全局坐標）進行推導，具體推導過程可參考清華大學曾攀老師的課程，直接從直角坐標（全局坐標）進行推導的過程省去了等參單元雅各比矩陣呀等坐標系映射的各種概念，理解起來相對容易。
展開 
四面體網格，六面體網格，低高階單元，對比研究
00 網格怎么選
四面體網格適應性強，自動化高。六面體網格雖然質量高，但劃分起來更麻煩。到底該怎么選擇？本文用一個例子進行對比研究。
01 幾何模型
02 部分網格展示
04 用低階六面體單元進行仿真計算
某兩點的位移隨節點數的變化趨勢：
某應力梯度較小位置的應力隨節點數的變化趨勢：
某應力梯度較大位置的應力隨節點數的變化趨勢：
05 用高階六面體單元進行仿真計算
某兩點的位移隨節點數的變化趨勢：
某應力梯度較小位置的應力隨節點數的變化趨勢：
某應力梯度較大位置的應力隨節點數的變化趨勢：
06 六面體單元的相關結論
01 位移結果可靠，節點數和單元階數的影響較??；
02 應力梯度較小位置的應力結果可靠，節點數和單元階數的影響較小；
03 應力梯度較大位置的應力結果不可靠，節點數和單元階數的影響較大；
07 四面體單元仿真計算與相關結論
01 高階四面體單元的位移結果可靠，節點數的影響較小；
02 低階四面體單元的位移結果不可靠，建議不要使用；
03 高階單元在應力梯度較小位置的應力結果可靠，節點數的影響較小；
04 低階單元在應力梯度較小位置的應力結果不可靠，建議不要使用；
05 應力梯度較大位置的應力結果不可靠，節點數和單元階數的影響較大；
08 總結論
01 在結構有限元分析中，建議不要使用低階四面體單元；
02 對于位移結果來說，六面體單元，高階四面體單元的求解都是可靠的，并且節點數影響較小。
展開 將單元zone轉化為四面體tetrahedron---block zone list poly (BBM)
1 引言
list命令一般用于屏幕顯示塊體或單元的屬性值，如果不附加寫入文件的命令，那么這些顯示的數據僅保存在內存中，不能直接保存到文件里面。但有兩個例外，一個是 block list poly <s>命令，這個命令產生出一個包含block create polyhedron face命令的文件，這個文件可以直接作為一個數據文件來調用。這個命令的主要用途是提取一個模型的可見塊內容，可以用來減少在建造過程中刪除了許多塊后模型所占的內存量。如果不指定文件名，那么系統會命名該文件為spoly.txt。有效位數等于12位或全局模型精度，以較大者為準。如果要增加精度到12位以上，使用model precision命令來設置。
另一個與之類似的命令是block zone list poly <s>，這個命令產生一個包含一系列block create tetrahedron命令的文件，其作用是將每個單元(zone)轉化成變成一個四面體塊(block)。原有的塊組、區域、材料號和本構號將應用于新的四面體塊。如果沒有指定文件名，那么系統會命名為 "poly.dat"。有效位數等于12或全局模型精度，以較大者為準。如果要增加精度到12位以上，使用model precision命令設置。
使用block zone list poly的主要用途是在FLAC3D或3DEC中產生粘結塊體模型[bonded block model (BBM)], BBM的生成由兩個步驟組成，第一步產生離散的塊體，第二步插入隨機的DFN或者輸入的節理組。下面分別演示這兩種途徑。
展開 基于ansys的梁單元、實體單元徐變精細化分析（含各參數解釋） ￥25
2、改網格模型，改成自己對應的網格模型，網格用ansys，hypermesh，ansa等前處理軟件都沒問題。
3、改材料參數，改成你想要的徐變模型，對著規范或者是你做出來的試驗擬合曲線。
以上即可實際應用。
ANSYS單元類型選擇方法 附ansys結構單元與材料應用手冊下載
下面是有關ANSYS分析中的單元選擇方法：
一、單元類型選擇概述：
ANSYS的單元庫提供了100多種單元類型，單元類型選擇的工作就是將單元的選擇范圍縮小到少數幾個單元上；
單元類型選擇方法：
1.設定物理場過濾菜單，將單元全集縮小到該物理場涉及的單元；
二、單元類型選擇方法
2.根據模型的幾何形狀選定單元的大類，如線性結構則只能用“Plane、Shell”這種單元去模擬；
3.根據模型結構的空間維數細化單元的類別，如確定為“Beam”單元大類之后，在對話框的右欄中，有2D和3D的單元分類，則根據結構的維數繼續縮小單元類型選擇的范圍；
三、單元類型選擇方法
4.確定單元的大類之后，又是也可以根據單元的階次來細分單元的小類，如確定為“Solid-Quad”，此時有四種單元類型：Quad 4node 42 Quad4node 183 Quad 8node 82 Quad 8node 183 前兩組即為低階單元，后兩組為高階單元；
四、單元類型選擇方法
5.根據單元的形狀細分單元的小類，如對三維實體，此時則可以根據單元形狀是“六面體”還是“四面體”，確定單元類型為“Brick”還是“Tet”；
五、單元類型選擇方法
6.根據分析問題的性質選擇單元類型，如確定為2D的Beam單元后，此時有三種單元類型可供選擇，如下：2D elastic 3 2Dplastic 23 2D tapered 54，根據分析問題是彈性還是塑性確定為“Beam3”或“Beam4”，若是變截面的非對稱的問題則用“Beam54”。
展開 ANSYS中桿單元和殼單元的單元耦合問題
在比較復雜的結構的有限元分析中，不同的結構部件通常使用不同類型的單元來模擬。
通常情況下，不同類型的單元的各個節點的自由度數目是不同的，不同類型單元的連接節點處的自由度的耦合問題，是一個比較令人頭疼的問題。
在ANSYS中通?？梢杂民詈厦頒P來耦合不同類型單元在連接節點處的自由度(DOF)。
也可以用CE命令來認為添加自由度之間的約束方程來達到耦合的目的。
下面是一個簡單的算例，使用了CE命令來耦合連接節點處的自由度。
模型是航天器的機翼的一個Section的某一個隔框。上下表皮是薄殼結構，用Shell63單元來模擬，在上下表皮之間有起支撐作用的桿件，用link8單元來模擬。
建模的時候，link8單元和shell63單元在連接有各自獨立的節點。即：link8單元和shell63單元的節點在連接處是重合的，但是，節點編號是各自獨立的。
link8單元在每個節點有 ux,uy,uz3個平動自由度；
shell63在每個節點有ux,uy,uz這3個平動自由度和rotx,roty,rotz這3個轉個自由，共6個自由度。
在耦合節點處，兩個耦合節點的ux,uy,uz自由度應該是相等的。
這個等式可以用CE命令來描述。
完整的命令流如下：
finish
/clear,start
/prep7
!定義第一種材料屬性；
mp,ex,1,30e6
mp,prxy,1,0.3
!定義shell63單元和實常數；
et,1,shell63
r,1,1e-3
!建立幾何模型；
rectng,31.8,33.2,0,0.3556
agen,2,1,1,1,0,0,1
a,1,4,8,5
a,6,7,3,2
KL,7,0.5, ,
KL,3,0.5, ,
在關鍵點處生成節點；
nkpt,100,4 !與編號為117的節點耦合
nkpt,101,9 !
展開 ANSYS各類型單元連接專題講解(五)之3D梁單元與殼單元剛接
例如采用ANSYS模擬一個多層混凝土框架結構，一般除計算整體指標外，我們在計算具體荷載作用時（如風荷載、地震作用、恒載、活載等），樓板一般采用彈性版，此時可用殼單元模擬，主梁、次梁采用梁單元模擬，此時變為梁單元包含在殼面內的情況，當然此類情況是否需要考慮截面偏置，可根據具體工程而定。
對這中梁單元包含在殼單元面內的情況，只需要將梁單元與殼單元共用節點即可，而無須格外建立約束方程。
三、梁單元在殼單元內但不包含
此種情況為梁與殼位于同一面內，但其中面不包含梁線，適用于多尺度建模分析（如下圖）。梁單元與殼單元的連接在端部可以通過剛性梁和剛性區域兩種方式連接。剛性梁采用MPC184單元，剛性區域采用Cerig命令，具體使用方法下期文章討論。
展開 ANSYS APDL實體單元和殼單元（不共節點）之間的連接 ￥100
實體單元和殼單元之間的連接是ANSYS中常見的問題。即使兩種單元之間共節點，但單元之間不連續（實體單元每個節點有3個平動自由度，而殼單元每個節點有3個平動自由度和3個轉動自由度），對于兩種單元之間面面接觸，可直接定義剛域，本文主要采用MPC法對實體-殼單元的連接方法進行說明。
1 單元類型
算例模型中，實體單元采用SOLID45，殼單元采用SHELL63，接觸位置不共節點。對于兩種單元之間的連接，通過目標單元TARGE170和接觸單元CONTA175實現，定義約束為實體-殼約束，接觸單元為MPC算法，接觸類型為綁定接觸。
2 有限元模型和綁定接觸
圖1 底部固定約束，殼單元施加均布荷載
圖2 目標單元和接觸單元
3 計算結果
圖3 von Mises stress
圖4 X-Component of displacement
付費內容為相關命令流。
展開 
ANSYS梁單元與實體單元的耦合與約束方程
ANSYS梁單元與實體單元的耦合與約束方程
By長安CAE
1 概述
在ANSYS計算過程中，有時候會遇到不同單元之間進行連接，由于不同的單元自由度不同，連接時通常需要通過耦合和約束方程建立節點自由度的關系，保證結果的準確性。
耦合可以理解成是將耦合的對象某個自由度作相等處理，而約束方程則不局限于相等這個關系，其可以描述具有某種關系的自由度。如圖1所示，為梁單元與平面單元的連接。如果不采用約束方程，力矩的傳遞無法完成，因為平面單元沒有轉動自由度。
圖1 梁單元與平面單元連接
為使節點2具有力矩傳遞的能力，要求1、2、3節點之間的自由度滿足以下關系：
ROTZ2 = (UY3 - UY1)/10
再通過CE命令，即可將此關系通過約束方程的形式施加給1、2、3節點。
2 命令
查看ANSYS的幫助文檔，查詢CE命令的解釋，如圖2所示。
圖2 ANSYS的CE命令解釋
CE, NEQN, CONST, NODE1, Lab1, C1, NODE2, Lab2, C2, NODE3, Lab3, C3
其中，NEQT表示常數，用于區別約束方程，一般可以用數字1、2、3表示即可，表示第幾個約束方程；
CONST表示方程的常數項，一般為0；
NODE1，表示第一個節點；
Lab1，表示自由度標簽，對于結構而言，就是三個平移和三個轉動自由度；
C1，表示該自由度的系數；
同理，后面的也一樣。
展開 在ANSYS 中3維坐標下的 shell structure 使用2D 平面單元劃分，應該使用哪個單元型號的單元
在ANSYS 中3維坐標下的 shell structure 使用2D 平面單元（僅考慮平面內的位移）劃分，應該使用哪個單元型號的單元？
ANSYS特殊單元——Follw201（隨動荷載）單元
ANSYS的Follw201單元是ANSYS的幾個特殊單元（比如mesh200）之一，稱為隨動荷載單元。都知道在ansys里面施加壓力載荷pressure時，其實載荷是可以隨動的，也就是能夠一直保持著面的法線方向，而施加集中了或者力矩時則不能保證。Follw201單元便能解決這個問題。
Follw201單元是一個單節點的3D單元，具有六個自由度，只能夠覆蓋在既有單元節點上，而且節點必須具有3個平移自由度和3個轉動自由度，也即是只能用在梁單元和殼單元上，實體單元僅有三個自由度。
Follw201單元主要用于幾何非線性分析問題中，在這類問題的分析過程中幾何會發生比較大的變形，面或線的法線方向可能發生比較大的變化，施加的載荷的方向是否隨動對結果的影響非常大。也就是用到此單元時會配對使用Nlgeon,on命令以打開大變形開關。如下圖所示為單元示意圖。
圖1
每個單元有兩個面，面1用于設定集中力的大小，面2用于設置力矩的大小，面的方向在應用時是通過單元的實常數進行定義的。
另外還需要注意，有限元求解的時候大部分是求解對稱矩陣，但隨動荷載單元的應用則包括了隨動荷載剛化效應，使剛度矩陣為非對稱的，因此需要采用非對稱求解器進行計算。
下面是具體應用，建一根梁單元，在梁的端部施加隨動集中力。
/prep7
!定義參數
EE=207E3
B=10
LCD=300
AA=B*B
IZ=B**4/12
PHZ=EE*IZ/LCD/LCD
!定義單元和材料
!201單元不需要定義材料
et,1,beam4
et,2,follw201
mp,ex,1,ee
mp,prxy,1,0.3
!定義實常數，實常數1設置梁單元的參數
r,1,aa,iz,iz,b,b
r,2,,1.0
!
展開 ANSYS不同單元類型連接專題（一）Solid-Beam單元的連接
通過對比兩次計算的結果發現：
1）全部使用Solid單元進行分析和使用Solid單元和Beam單元連接起來進行分析，
計算結果幾乎完全一致；
2）使用Solid單元和Beam單元建模和全部使用solid單元進行建模相比，節點數量大大減少，
顯著
降低了計算量。
三、連接原理。
在Workbench中，我們很容易就建立了solid-beam的連接，那么，軟件究竟是根據什么原理建立的呢？我們去ANSYS經典中一探究竟。
通過查看單元類型我們發現，ANSYS生成了計算用的5種單元類型。而我們沒有定義接觸，怎么會有接觸單元174和目標單元170呢？
通過查看接觸向導我們發現，ANSYS生成了一個
單點控制接觸，控制節點為173184?？吹竭@我們就大概明白了，在梁模型和實體模型接觸的位置，軟件建立了一個170點目標單元，在實體模型的端面上，軟件建立了174單元。使用170單元的173184節點控制174單元上的節點。
那么端面上的實體單元又是怎么和梁單元連接的呢？我們發現，還有一個
MPC184單元沒派上用場。我們單獨顯示MPC184單元，發現它連接了173183和173184節點，173184就是我們剛才提到的控制節點，而173183為軟件在梁模型的端點上建立的170單元上的節點。
至此，本文完結。
歡迎大家點擊在看和轉發支持！掃描二維碼關注公眾號，一起聊聊力學和有限元那點兒事。
展開 