abaqus變量輸出變量的案例
ABAQUS能量平衡輸出變量
Total energy output quantities
ALLAE
“Artificial” strain energy associated with constraints used to remove singular modes (such as hourglass control), and with constraints used to make the drill rotation follow the in-plane rotation of the shell elements.
ALLCD
Energy dissipated by creep, swelling, and viscoelasticity.
ALLEE
Electrostatic energy.
ALLFD
Total energy dissipated through frictional effects. (Available only for the whole model.)
ALLIE
Total strain energy. (ALLIE = ALLSE + ALLPD + ALLCD + ALLAE + ALLQB + ALLEE + ALLDMD.)
ALLJD
Electrical energy dissipated due to flow of electrical current.
ALLKE
Kinetic energy.
ALLKL
Loss of kinetic energy at impact. (Available only for the whole model.)
ALLPD
Energy dissipated by rate-independent and rate-dependent
展開 Abaqus_Python 后處理結點數據提取、場變量云圖輸出 ￥1.99
今天體會了Python腳本操作ODB文件進行數據的輸出、處理等方面的快速高效。
首先介紹一下背景，同事的項目，計算case較多。每個case需要輸出幾個路徑上的部分節點值，而且需要針對不同的case、不同的幀、不同的視角截圖、保存云圖，如果人工出圖和提取節點數據，工作量會非常大，費時會遠遠超出計算所用時間；并且，同事還要求我根據case名稱、組件名稱（模型有多個組件組成，且網格是在part上畫的，所以同一個節點編號可能對應多個節點）等保存成不同的文件。
結點溫度數據沒有保存，因此沒有輸出。下圖是自動輸出的png格式云圖。
比較倉促，代碼也比價粗糙。回頭做一個詳細的介紹。
主要代碼見下。
展開 ABAQUS中隱式和顯式的節點和單元的輸出變量解析
FSLIP
Field: yes History: no .fil: no
Length of contact slip path at secondary nodes during contact (FSLIPEQ) and in some
cases (see About contact pairs in Abaqus/Explicit) components of net contact slip in local
tangent directions (FSLIP1 and FSLIP2). These variables remain constant while a
secondary node is not in contact.
FSLIPR
Field: yes History: no .fil: no
Magnitude of contact slip rate at secondary nodes during contact (FSLIPR) and in some
cases (see About contact pairs in Abaqus/Explicit) components of contact slip rate in local
tangent directions (FSLIPR1 and FSLIPR2). These variables are set to zero while a
secondary node is not in contact.
BONDSTAT
Field: no History: yes .fil: no
Spot weld bond status.
展開 ABAQUS歷史輸出中，各能量變量（ALLAE、ALLIE、ETOTAL等）意義
ABAQUS中，對于很多動態問題，尤其像高速沖擊模擬中，對結果評價很重要的一點就是要保證模型能量守恒，這就涉及到ABAQUS歷史輸出中各能量變量的意義，下面最各簡單整理：
ALLAE：人工偽應變能，六面體、殼網格中沙漏發生情況指標
ALLCD：蠕變、膨脹以及粘彈性消耗的能量
ALLFD：摩擦消耗的能量
ALLIE：總的內能，ALLIE=ALLSE+ALLPD+ALLCD+ALLAE+ALLQB+ALLEE+ALLDMD
ALLKE：動能
ALLKL：碰撞消耗的能量
ALLVD：粘性消耗的能量
ALLDMD：破壞消耗的能量
ETOTAL：所有能量的總和
NOTE：ALLAE（偽應變能）的理解：偽應變能就是控制沙漏變形所消耗的主要能量。如果偽應變能過高，說明過多的應變能被用來控制沙漏變形了。一般通過比較偽應變能和其他內部能量的值來判斷偽應變能是否過高，以及判斷過高的偽應變能的來源。
展開 
ABAQUS案例-場變量的應用及材料彈性模量隨場變量而變化 ￥3
ABAQUS中的場變量具有較高的應用價值，可以在一些復雜的工程應用中極大的減輕工作量。本實例即是展示一個場變量應用——材料彈性模量隨場變量而變化，其中它涉及到關鍵字的編輯（關鍵字的具體編輯也在附件中）。本實例在附件的inp文件中。
ParaView---FLAC3D和3DEC輸出變量到VTK文件
1 引言
在《免費二維有限元分析軟件ADONIS計算步驟》中曾經提及ADONIS的計算結果可以輸出為VTK文件，事實上3DEC和FLAC3D的計算結果也可以輸出為VTK文件。這個筆記簡要討論了輸出步驟。
2 VTK命令
3DEC和FLAC3D使用block vtk和zone vtk命令輸出vtk文件，可以輸出的變量有：
(1) displacements: 網格點的位移矢量
(2) zone-extra: 單元索引i的額外變量
(3) gridpoint-extra: 網格點索引i的額外變量
(4) forces: 網格點不平衡力的信息，包括最后的局部應力比【FLAC/Slope的強度降低過程(Strength Reduction Procedure in FLAC/Slope)】
(5) groups: slot中單元設置的組
(6) model-fluid: 流體屬性
(7) model-mechanical: 力學屬性
(8) model-thermal: 熱屬性
(9) pore-pressure: 網格點的孔隙壓力值
(10) range: 設置range, 用來控制哪些單元包括在vtk文件中
(11) saturation: 網格點飽和度的值
(11) stresses: 每個單元的平均應力。不輸出內部四面體的應力。
(12) temperatures: 網格點的溫度值
(13) velocities: 網格點的速度矢量
(14) filename: 設置vtk的文件名
默認情況下，只輸出應力(stresses)，組(groups)和位移(displacements). vtk的文件格式是XML語言。
展開 ABAQUS UMAT調用后為什么損傷變量出現負值？
ABAQUS UMAT調用后為什么損傷變量和單元刪除出現負值和大于1得值？
調用子程序后，計算復合材料損傷過程，損傷變量和單元刪除出現負值和大于1得值，真誠求助。附上子程序，請老師指教
20251203.txt
【公益帖】VUMAT 中間變量輸出到指定文件中的方法
背景： 在進行用戶材料本構用戶子程序調試過程中，我們需要輸出一些中間變量，跟蹤其演化規律，用于debug， 在umat中我們往往用 write（6，*） write（7，*）寫入到相應的msg或者dat文件中，但是在VUMAT中這種方法卻不可以，下面提供一種VUMAT輸出變量的方法
需要在vumat子程序中添加如下語句
open(1,file='絕對路徑\output.txt')
write(1,*) 變量
這樣會在指定路徑下生成一個名為output.txt，希望對大家有用
【LSDYNA操作小技巧七】LSDYNA中具體材料歷史變量的定義與輸出問題 ￥1.99
LSDYNA作為專門針對材料非線性、大變形問題而開發的大型顯示動力學仿真軟件，其擁有相當齊全的材料庫，針對其中材料的本構關系都有非常明確的軟件設置與輸出步驟，方便用戶對本構的快速設置甚至二次開發。本帖以典型的材料大變形問題磨粒切削加工引起的加工損傷為案例講解損傷歷史變量的定義與輸出步驟。
在我們沒有更好的文獻參考時，官方仿真軟件的幫助文檔是一個最佳的選擇。以下提供在LSDYNA中具體查閱歷史變量定義的鏈接如下，
History Variables for Certain Material Models — Welcome to the LS-DYNA support site.
基于此，必須首先定義目標歷史變量如圖1所示。具體通過，在關鍵字 *DATABASE_EXTENT_BINARY中定義NEIPH 或 NEIPS 實現。這里NEIPH 或 NEIPS 只能填入具體的數字加以定義材料具體哪些的歷史變量。
圖1
具體以此為例：為輸出磨粒切削加工引起的工件損傷，首先找到對應的工件材料編號110（工件為K9玻璃，采用JH-2本構），可知損傷因定義為2如圖1所示，故在NEIPH輸入數字2即可，若為了查看其他材料歷史變量，對應輸出1-4即可實現。
圖2
完成定義后，完成模型建模輸出為k文件，在LSPP中進行損傷云圖的輸出。圖3為定義的損傷在后處理中的查看方式，圖4為輸出的損傷云圖，紅色粒子表示損傷。
圖3
圖4
以上，基于此案例，可在LSDYNA軟件中定義并輸出所給材料的歷史變量。
展開 ABAQUS模擬多道次變形的變量繼承方法
一、引言
使用ABAQUS進行多道次加工時，往往牽扯道次之間變量的繼承（如晶粒尺寸、累積損傷等），這對多道次變形模擬結果的準確性有較大的影響。本文以VUHARD子程序及簡單的熱壓縮模型為例，分享雙道次壓縮之間的晶粒尺寸的繼承方法。
Abaqus后處理-云圖變量含義（部分）
結果中分量說明：
S11、S22、S33指各軸正應力；
S12指作用于XZ平面(與“2”，即Y軸垂直的平面)內，沿1方向剪應力；
S13指作用于XY平面內，沿1方向剪應力；
S23指作用于XY平面內，沿2方向剪應力。
若為柱坐標，S12、S13、S23分別指：由徑向向環向的剪應力、由徑向向軸向的剪應力、由環向向軸向的剪應力。
abaqus后處理插件—場變量結果標注 ￥45
abaqus后處理界面中探針功能附帶的標記樣式非常丑陋，基于abaqus的試圖注釋功能進行二次開發，形成了場變量標注插件，方便快速的標注關心區域的應力應變等結果。
插件介紹：
按鈕介紹
從左至右依次是：標記按鈕、隱藏標記按鈕、恢復顯示按鈕、刪除按鈕
示意動畫
使用方法：
1） Probe查看節點結果，并勾選需要標記的節點項；
2）點擊工具欄中的標記按鈕，進行標記。
特點
1） 標記速度快，即使在單元數目達到百萬級及以上的模型中，標記速度仍無明顯延遲；
2）所有標記注釋均在試圖注釋功能界面里，有利于對美觀度有更高要求者進一步修改美化。
展開 Abaqus python腳本開發 第三章 各類指令的方法對象變量
該部分主要介紹以Python scripting guidance為基礎的對象方法變量，在一些較難理解的模塊部分提供代碼示例，供讀者參考。
3.1 Mdb command
3.1.1 Mdb object
Mdb對象是高級Abaqus模型數據庫，其中儲存了模型和分析控制信息。
3.1.1.1 Mdb()
建立一個空的Mdb對象
可選變量
pathName
默認保存為.cae文件
返回值
Mdb對象
3.1.1.2 importDxf()
從AutoCAD .dxf 文件建立一個ConstraintedSketch對象
需要的參數
fileName
返回值
Mdb對象
importDxf(fileName=‘’）
3.1.1.3 openMdb()
打開現有的model databse file
必要參數：
pathName
一段用來指定打開路徑的字符串， 如果不提供擴展名的話，Abaqus/CAE會打開帶 .cae后綴的文件。
展開 ABAQUS變量解讀：教你讀懂應力/應變/損傷
結語
ABAQUS中的變量雖多，但歸根結底都是為描述材料行為和判斷結構狀態服務的。理解每個變量的物理意義，結合實際分析目標選擇合適的輸出，才能讓仿真結果真正“說話”。
如果你對某個變量還有疑問，或者想了解更深入的應用場景，歡迎在評論區留言！
getMask()討論：Abaqus Part對象的成員變量常用函數 ￥1
這樣可以有效的減少輸出結果文件的大小。
本以為這應該是一個極其簡單的小問題，不會浪費我太多的精力。但是當開始寫程序的時候才發現一個巨大的坑點：Part對象的Skin函數不接受列表類型的參數或者數組類型的參數，它只接受ABAQUS特定的幾何數組類型或者序列類型，按理說list類型屬于序列類型，但是該函數并不接受。如果通過getSequenceFromMask生成可用faces，mask究竟如何獲得又是一個問題。
mask為什么難以獲得呢？因為getMask()函數是FaceArray對象的成員函數，而該對象的建立是通過part對象建立的（p.faces），它會生成所有faces的合集（比如[0,1,2,3,4,5,6]），而我要進行skin的cell只包含這些faces中的一部分，并且編號是隨機的不可預測的，但是組成該cell的face編號可以通過c.getFaces()（(2,3.4,5)）。按照我最初的想法，我有了這些面的ID號，把它們從p.faces中讀取出來重新組合就可以了，但是組合完成的類型是list，Skin函數不接受它作為參數。我嘗試了多種方法均不能將最終的結果轉化成Skin可以使用的格式。并且由于組合后的不是abaqus的幾何數組類型，所以它并沒有getMask()方法，也就無法通過該方法獲得可以使用的mask。
無奈之下只好使用getSequenceFromMask函數，這就要解決一個麻煩的問題：getSequenceFromMask的mask值究竟怎么來？目前看來唯一的辦法就是解密getMask()函數。（未知原因，公司網絡連接不了ABAQUS幫助文件。）
展開 