單元編號
關注創建者：320科技工作室 創建時間：2022-02-06
單元編號的視頻教程
使用Python批量提取多個odb的集合里面的Mises應力并寫入excel里面
本視頻詳細介紹了如何使用python代碼批量提取多個odb中的集合里面的Mises應力隨時間的變化關系，一方面可以在集合中實現提取，另一方面可以在單元編號里面實現提取，對批量處理操作有極大的幫助，可以大量節約人力物力財力。
￥99 8分鐘 40播放
查看
Abaqus土木仿真簡單教程（持續更新）
課程目錄： 章節1：Abaqus輸出骨架曲線和滯回曲線； 章節2：滯回曲線中提取骨架曲線數據的2個方法（Matlab代碼和exe程序）； 章節3：Abaqus與VS、IVF關聯；pqfiber和cs使用問題總結； 章節4：Abaqus多版本共存問題； 章節5：Abaqus與PythonReader； 章節6：Visualization中對比2個odb結果； 章節7：查看節點單元編號
免費 43分鐘 12081播放
查看
ABAQUS基坑非飽和滲流詳解
通過交互導出inp文件的方法對單元進行分組，確定單元面的編號，并從單元拓撲關系的角度介紹了確定單元面編號的方法。 講解了滲流分析中存在的四種滲流邊界，著重講解了自由溢出邊界設置的兩種方法，并通過基坑管井降水模型介紹操作過程。 講解了管井井點定流量抽水的實現方法，分析了管井井點定流量抽水過程中土體孔壓、位移變化的時間特征。
￥40 39分鐘 1275播放
查看
單元編號的實例教程
在編寫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勾選上，可以視覺上顯示出來單元刪除的效果。結果如下圖所示：
需要注意的是如果一個單元有多個積分點，則循環中需要跳過一部分積分點。該方法需要單核計算才能使用。
最后有子程序開發等相關需求，歡迎聯系我們。
展開 解決的方法：在command line輸入 mdb.models['modelName'].setValues(noPartsInputFile=ON) 重新生成inp文件，這樣在寫出的inp文件里面就沒有instance信息，而且單元節點編號也不會有重復了。
將單元i的面積存入數組
k = k+1
*endif
*enddo
但是此方法存在一個問題，命令流將提取單元編號從EMIN到EMAX的所有單元面積，而不能保證單元編號對應的單元屬于當前選擇的單元集合。例如，本例中選擇單元1、2和12，執行上述命令流時，程序仍然會提取單元編號從1到12一共12個單元的面積。
方法二
針對方法一存在的不足，筆者進行了改進，整體思路如下：
1、采用單元個數作為循環控制變量；
1
*do, i, 1, ENUM
2、采用 elnext() 函數獲取下一個單元編號；
3、將單元編號和單元面積存入數組后，采用*vwrite命令將數組寫入文本文件；
本例寫出的單元編號和單元面積結果如下
特別說明：*vwrite命令不能在命令流窗口中運行，而必須在ANSYS Batch模式下運行。
小結
經筆者改進的方法二，可在ANSYS中一次性提取任意數量、編號不連續平面單元的面積，并將結果保存到文本文件，方便后續的繪圖或是統計工作。
展開 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.
展開 **
**節點集合和單元集合的表示方法：
**1.如果集合中的節點或單元編號是連續的，可以表示為
**節點集合： *Nset， Nset=節點集合名稱， Generate
** 起始節點編號，結束節點編號，節點編號增量
**單元集合： *Elset， Elset=單元集合名稱，Generate
** 起始單元編號，結束單元編號，單元編號增量
**注意：節點集合和單元集合的名稱不得超過80個字符，必須以字母或下劃線開始。
**2.如果集合中的節點或單元編號是不連續的，表示方法為：依次列出集合中的所有節點或單元（每個數據行接點或單元編號不得超過16個），格式：
**節點集合： *Nset， Nset=節點集合名稱
節點編號1，節點編號2，...節點編號16
**單元集合： Elset，Elset=單元集合名稱
單元編號1，單元編號2，...單元編號16
*Nset, nset=_PickedSet2, internal, generate
1, 281, 1
*Elset, elset=_PickedSet2, internal, generate
1, 80, 1
** Region: (Section-1,Picked)
*Elset, elset=_PickedSet2, internal, generate
1, 80, 1
**截面屬性的基本表示方法為：*Solid Section， Elset=單元集合名稱， Material=材料名稱
** 界面參數
**界面參數可以是二維模型的厚度或一維模型的截面面積等。
展開 
單元編號的相關專題、標簽、搜索
單元編號的最新內容
考慮熱源的瞬態熱傳導有限元求解器22小時前
那么前端數據在生成求解器輸入的時候，就要告知求解器所有單元的編號和其對應的體熱功率。
當求解器拿到單元編號以后，就需要索引或者計算其面積，并根據單元三個節點編號，將功率加到載荷列陣對應的位置。
驗證
設計案例如下，區域外部為20℃空氣，對流換熱系數取5W/(m2K)，時間總長18000s，每步時間間隔60s。
由于單元是按層依次生成的，每層的單元編號是連續的，因此可以通過單元編號的范圍來定義每層的單元集合。
ABAQUS支持使用generate關鍵字來定義連續編號的單元集合，只需要指定起始編號、終止編號和步長即可。
對于N層的層合板，會生成N個單元集合，每個集合對應一層材料。這種按層分組的方式為后續的材料屬性指定提供了清晰的邏輯結構。
6.
與求解器耦合時，確保網格拓撲、單元類型、節點編號在內部和外部求解器間一致。</p><p><br></p><p><strong>四、材料與物理性質模塊（Materials & Physics）</strong></p><p>1.內置材料模型庫（線性/非線性彈性、塑性、粘彈性、粘塑性、損傷、疲勞、斷裂等），以及溫度、速率、熱-結構耦合效應。
：若你只想查看特定編號范圍的單元信息，可使用ELIST,IMIN,IMAX命令，IMIN和IMAX分別為起始和結束的單元編號。
image_process=/format,webp" data-initial-src="https://img.jishulink.com/202502/attachment/ba2f6debee4641a8b9730513ee59e15c.png">
</figure>
</figure><p>具體每根桿單元的節點編號以及材料屬性，在input文件中詳細介紹。
首先對此結構進行單元剖分，確定單元按照有限元的分析流程，要先對此結構進行單元剖分，確定單元與結點編號、以及單元的自由度編號。因為這里是平面應力問題，所以可以采用常應變三角形單元進行網格劃分，并且采用的是非結構化的網格。
</p><h2>3．桁架單元和節點的編號以及結點坐標的定義</h2><p>l 黑色為節點坐標</p><p>l 紅色為節點編號</p><p>l 紫色為單元編號</p><p><br></p><p><br></p><div contenteditable="false" width="100%">
<figure class="figure-image" data-img
生成的嵌入cohesive單元的新part中節點和單元的編號從1開始并連續，方便后處理， 和一些其他方法中分裂節點時加一個大數的編號方法不同。
</p><p>建一個沙漏控制卡，選4號或6號，附給泡沫單元的part6</p><p>接觸中將soft改為1，將sfs和sfm改為0.1</p><p> </p><p>負體積的原因是由于單元畸變引起的，單元節點編號有一個順序，當變形過大，或者不合理時， 某個或某些節點穿透所屬單元的面，造成負體積。
本案例是基于tcl語言實現每個單元的id號及對應單元的節點號，創建一個列表。具體實現過程見本案例的程序部分。
以下是輸出的格式：
1339780 {1394271 1394272 1394273 1394274}
詳情見收費的程序部分，凡購買本案例的朋友針對該案例有疑問，可私信，謝謝！
