單元編號處理的案例
基于tcl語言實現單元編號及節點編號的列表 ￥15
本案例是基于tcl語言實現每個單元的id號及對應單元的節點號，創建一個列表。具體實現過程見本案例的程序部分。
以下是輸出的格式：
1339780 {1394271 1394272 1394273 1394274}
詳情見收費的程序部分，凡購買本案例的朋友針對該案例有疑問，可私信，謝謝！
以C3D8R為例介紹如何在VUMAT中得到單元編號
VUMAT不像UMAT直接給出了單元的編號， VUMAT需要借助一個子程序vumatXtrArg 得到，這個方法最早發布在如下鏈接中
https://polymerfem.com/community/constitutive-models/nblock-in-vumat/
相信有不少朋友需要，摘錄如下
You are using C3D8R elements which means continuum 3D 8 noded reduced integration. A reduced integration of this kind has a single Gauss pt.
Thats why you have matching GP and element numbers.
展開 ABAQUS等有限元軟件的后處理中節點編號排序matlab程序 ￥5
以ABAQUS為例，在進行ABAQUS的節點信息后處理時，我們通常要分析，選取大量的節點，而我們在建模過程中節點的順序往往是不跟隨我們需求的，提取節點的速度、加速度、位移等數據并進行繪圖時，將節點編號與節點位置統一起來比較麻煩，在這里我會使用一個matlab小程序來調整節點編號與我們需要的空間位置進行對應。主要分為以下步驟
1.在ABAQUS中，選擇你要輸出的節點信息，通過report-xydate進行rpt文件的輸出。
2.對ABAQUS中的節點進行節點信息查詢，記錄節點編號信息。
3.使用文本文檔/notpad++將rpt文件打開，放到excel中
4.在excel中使用分列，將數據分開，并刪除第一行中沒有用的部分，以及第一列中的時間列，只保留節點編號與其對應的加速度/速度/位移時程等的變化。
5.使用matlab讀取文件位置，將你想要的正確的順序輸入matlab程序中，運行程序即可得到你想要的按順序編號的excel文件。
展開 Vusdfld\vumat等顯示動力學子程序獲取單元編號的一種方法
在編寫vusdfld\vumat等子程序的時候，有時需要知道單元編號，以便對單元直接賦值，不同于usdfld，umat子程序，他們的單元編號信息會傳入子程序中。但是顯示動力學對應的子程序傳入的是積分點編號，并沒有單元編號。可以有一種取巧的方法來實現積分點和單元編號的對應。具體思路如下：
以邊長為50mm立方體為例，網格種子密度為5，共1000個單元，單元類型C3R8R。而在顯示動力學中，使用單核運算時，abaqus會一次性傳入136個積分點，1000個單元會傳入8次，最后一次為48個積分點，可以通過write(*,*) ‘nblock’,nblock輸出來看傳入積分點的變化情況。因此可以在程序中設置一個全局變量，每8次循環后就歸0，重新循環。利用狀態變量statenew(k,1)來存儲單元。
在后處理中顯示單元編號，可以看出云圖從1到1000，正好對應所有單元的個數，利用查詢值功能，可以看到單元ID和狀態變量的值是相同的。
知道單元編號之后就可以對單元進行賦值操作，比如讀取一個隨機材料場文件，然后將數據賦值給單元。還可以利用該功能，實現類似“生死單元”的功能。一開始就根據單元編號抑制一部分，在之后的分析步在激活。比如第一個分析步殺死一些單元，然后下一個分析步在激活，比如刪除編號為501的單元。
可以看到編號501的狀態變量為值為0，單元應力為0，說明該單元被刪除，不參與運算過程。如果在場輸出中把STATUS勾選上，可以視覺上顯示出來單元刪除的效果。結果如下圖所示：
需要注意的是如果一個單元有多個積分點，則循環中需要跳過一部分積分點。該方法需要單核計算才能使用。
最后有子程序開發等相關需求，歡迎聯系我們。
展開 
輸出abaqus inp如何去除instance信息，而且單元節點編號也不會有重復
解決的方法：在command line輸入 mdb.models['modelName'].setValues(noPartsInputFile=ON) 重新生成inp文件，這樣在寫出的inp文件里面就沒有instance信息，而且單元節點編號也不會有重復了。
ABAQUS喵星人教你學會鋼筋混凝土殼單元的前處理與后處理
ABAQUS中的殼單元大家通常用于模擬鋼板等鋼結構，對于混凝土板殼，新手可能對內部的配筋方式，以及前后處理方法可能存在各種問題。實際上，ABAQUS提供了鋼筋混凝土板配筋的接口，這種“寫入式”而不進行直接建模的方法通常比較冷門且后處理相對不主流。今天喵星人就通過一個教程教你學會鋼筋混凝土殼單元的前處理與后處理。
0.前提
使用板殼單元的有限元模擬必須有兩個前提：
1、板殼力學及殼單元通常應用于一個方向尺寸遠小于另外兩個方向（通常不超過1/5）的結構。
喵星人點評：大家總有一個誤區，總覺得實體單元的精度最高，實則不然。對于板殼結構，由于其采用了Kirchhoff板假定，在此情況下相比實體單元，殼單元形函數更加逼近實際結構，其計算精度與計算代價均優于采用實體單元。
2、由于采用Kirchhoff板假定，即忽略混凝土板中鋼筋的粘結滑移行為，因此在精細化的鋼筋混凝土滯回模型中通常不再適用。
1、前處理
1.1 縱橫方向與局部坐標系
配筋的板殼單元，尤其是兩個平面方向差異配筋的板殼單元，必須指定坐標系，且喵星人建議使用局部坐標系。這是為了避免在裝配件中因旋轉導致整體坐標系的變換。本案例中的坐標系指派如圖所示。需要注意的是，鋼筋縱橫方向與局部坐標系方向直接掛鉤。
1.2 配筋面積/間距/方向
殼單元的配筋方法需在“編輯截面”中完成，不能直接建立線單元鋼筋。采用“寫入式”的建模方法，如下圖所示。
其實這種方法很像設計軟件中的操作，即通過加勁的方式考慮配筋混凝土。
展開 平面單元和3D單元在處理幾何非線性有何區別？
分別用C3D8R單元和CPS4R單元進行模擬，當關閉大變形時，得出的兩條曲線完全重合；打開大變形后，兩條曲線差別很明顯，應變越大，差距越大。請問版主，這是什么原因呀？
abaqus在處理幾何非線性的時候，兩種單元的處理方式有什么區別呢？
請版主幫我想想辦法，謝謝。
線單元后處理方法
顯示應力和彎矩
顯示梁單元應力和彎矩.doc
顯示線單元應力和彎矩.doc
ANSYS 有限元分析后處理結點解與單元解!
3.4求解
四、后處理
ANSYS 提供了兩個后外理器:通用后處理器(POST1)和 時間歷程后處理器(POST26)。通用后處理器(POST1):用來觀察整個模型在某一時刻的結果。時間歷程后處理器(POST26):用來觀察整個模型在不同時間段或荷載步上的結果，常用干處理瞬態分析和動力分析結果。本算例為靜力分析，因此，該模型的后處理主要用到 POST1 處理器。
4.1顯示變形形狀
4.2顯示位移云圖
PLNSOL 為用等值線或云圖的方式顯示結點處的計算結果;PLESOL為用等值線或云圖的方式顯示單元的計算結果。
4.3顯示應力云圖
4.3.1顯示連續應力云圖
4.3.2顯示非連續應力云圖
本文內容來源：Hulunbuir博客
展開 CAE前處理 | 高階單元在薄板網格劃分時的注意事項（2）
01
前言
在文章【CAE前處理 | 高階單元在薄板網格劃分時的注意事項（1）】中，筆者對比了不同長厚比下，厚度方向網格數量對薄板結構的剛度及強度影響
根據計算結果初步判斷，1層高階全積分單元是能夠滿足薄板結構常規計算需求
這里可能有伙伴會想，“高階單元既然精度這么高，豈不是網格隨便劃分下就能進行計算？”
這里暫且不討論其它，單就薄板結構網格劃分而言，還有很重要的一部分數據沒有進行對比，那就是“長度方向網格數量對計算精度有著怎樣的影響？”
Abaqus前處理插件-生成Voronoi單元集 ￥100
通過本款插件可方便生成類Voronoi的單元集，二維圖形和三維結構通用，適于任意幾何構型。
閉源程序保證支持2016及以下版本，高版本建議采用開源程序：
https://www.yqgqt.org.cn/content/post/f6b5ef4e-5878-43cb-bb9a-9cfd4dcb05b8
更新日志：
2020.7.19 ： 改進對早于2016的舊版本兼容性。
Abaqus前處理插件-生成Voronoi單元集(開源程序) ￥200
通過本款插件可方便生成類Voronoi的單元集，通用于二維圖形和三維結構，適于任意幾何構型。
包含程序的開源 .py 文件，可直接作為腳本文件執行，相信能夠幫你打開Abaqus二次開發的思路。
更新日志：
2020.7.19 ： 改進對早于2016的舊版本兼容性；
2020.9.15 ：增加隨機程度控制參數 (0~1)。
二維實例，隨機程度 = 0：
三維實例，隨機程度 = 0：
二維實例，隨機程度 = 1：
三維實例，隨機程度 = 1：
程序界面：
如有需求可私信詳聊。
Abaqus前處理插件-生成Voronoi單元集(開源程序) ￥150
通過本款插件可方便生成類Voronoi的單元集，通用于二維圖形和三維結構，適于任意幾何構型。
包含程序的開源 .py 文件，可直接作為腳本文件執行，相信能夠幫你打開Abaqus二次開發的思路。
更新日志：
2020.7.19 ： 改進對早于2016的舊版本兼容性。
如僅需閉源文件，請移步：
https://www.yqgqt.org.cn/content/post/19d81779-6d9a-4fd5-bb83-758694714c89
使用MPP并行計算處理ALE單元時
<p>并行計算最好采用MPP形式，以獲得良好的并行處理性能。使用MPP并行計算處理ALE單元時，需要設置如下關鍵字</p><p>*CONTROL_MPP_DECOMPOSITION_DISTRIBUTE_ALE_ELEMENTS</p><div contenteditable="false" width="100%"><img src="https://img.jishulink.com/upload/202210/43a0f8b10ebe4d07b840e80560730b10.png" title="mpp.png" alt="mpp.png" style="max-width:760px;" data-mobile-src="https://img.jishulink.com/upload/202210/43a0f8b10ebe4d07b840e80560730b10.png?image_process=/format,webp/resize,w_400" data-pc-src="https://img.jishulink.com/upload/202210/43a0f8b10ebe4d07b840e80560730b10.png?image_process=/format,webp/resize,w_760" data-initial-src="https://img.jishulink.com/upload/202210/43a0f8b10ebe4d07b840e80560730b10.png">
</div><p><br></p>
展開 CAE前處理 | 選擇合適的網格密度 | 實體單元（1）
進行有限元分析時使用者很容易陷入兩個極端，一是生怕網格數量不夠而影響計算精度，使用近乎變態的網格尺寸進行分析，二是輕視網格的重要性，直接盲目地使用軟件默認的網格進行處理。前者很容易帶來巨大的計算量從而加重分析負擔，后者經常拿出不合理的分析結果但卻不自知，其中實體單元的分析相對于梁，殼來說更容易出現上述問題，這也是該系列文章準備著重探討的地方。
問題的產生
如圖三種結構分別代表了狹長實體，常規實體，薄壁實體的特征，假設現在需要對這些結構進行有限元計算來得到其剛度和強度性質，那么請問整體和局部至少需要使用多少網格才能較好的捕捉到關注的問題？
說實話，在寫這系列文章時筆者并不知道具體多少合適，因為和大多數學習者一樣，網格數量的多少大部分時候憑感覺，反正就是：這個網格量應該夠了！顯然，這樣是非常不嚴謹的。
但是，要對這一問題進行探索并不容易，畢竟不像梁單元，實體結構從受力模式和網格維度方面都要復雜得多，因此，文章內容不可能遍歷所有結構特征，只能針對一些典型結構進行對比，得到一些基本的分析規律。
另外需要主要，文章基于的求解器為OptiStruct，不同求解器不同單元類型可能會得到不太一樣的結果，所以對于不同的求解器需要針對對應的問題進行各自規律的探索。
展開 