abaqus用戶變量的案例
在Linux(Centos7)下的用戶子程序調用的環境變量設置 ￥10
目前為止，網絡上依然沒有在linux系統下用戶子程序調用的環境變量設置與關聯，對Linux系統及其命令陌生的部分學生和科研人員而言，是個較大的挑戰。
為此，本文將以Centos 7 這一Linux操作系統為例，列出含用戶子程序調用的調試流程，以便給大家一些指導。
特別聲明，此文僅供學生、科研等人員本人參考、學習使用，請勿他用。
DEFORM二次開發后用戶單元變量初始值批量導入工具WriteUserVar ￥19.89
軟件功能：
方便Deform軟件二次開發過程后模擬設置時用戶變量的輸入問題，不必每次逐個輸入用戶變量的名稱和初始值，大大簡化了模擬工作。
使用說明：
1. 先在deform中進行模擬的參數設置，全部設置完以后，導出key文件：File – save as 保存到一個地方。
2. 打開軟件 input key 選擇剛剛保存的key文件，打開input USRV 選擇自己的用戶變量列表（見本文最后）。
3. 默認的object為1 ，可自行更改每次生成一個object的變量表，點擊RUN即可添加用戶變量，完成后會提示finish的提示。
4. 點擊 SA VE即可保存生成的key文件（注意，這里為了方便一些deform文本模式的用戶，所以直接SA VE得到的key文件中是只有生成的單元用戶變量表信息，直接在deform里import key是可以的。如果想在之前生成的key文件中加入該信息，保存時選擇原key即可，替換掉原來的key就可以了）。
5. 點擊deform中的File – import keyword，找到剛剛生成的key文件即添加成功。
6. 在deform中的Input – Object element中找到User可以看到剛剛導入的用戶變量（注意先選擇相對應的部件，不同部件都可以進行定義用戶變量需要用戶進行區分）
注：用戶變量表的設置的說明：
請自行創建一個txt文檔，在文檔中每行輸入一個用戶變量的名字（如 TEMPERATURE_C）和一個初始值，名稱中不能有空格等若需分隔請加下劃線_，變量名和初始值之間用空格進行分隔。
展開 ABAQUS案例-場變量的應用及材料彈性模量隨場變量而變化 ￥3
ABAQUS中的場變量具有較高的應用價值，可以在一些復雜的工程應用中極大的減輕工作量。本實例即是展示一個場變量應用——材料彈性模量隨場變量而變化，其中它涉及到關鍵字的編輯（關鍵字的具體編輯也在附件中）。本實例在附件的inp文件中。
ABAQUS UMAT調用后為什么損傷變量出現負值？
ABAQUS UMAT調用后為什么損傷變量和單元刪除出現負值和大于1得值？
調用子程序后，計算復合材料損傷過程，損傷變量和單元刪除出現負值和大于1得值，真誠求助。附上子程序，請老師指教
20251203.txt
一個abaqus用戶材料子程序（UMAT）問題解決說明 附Abaqus用戶子程序專題下載
在開展abaqus用戶子程序的項目時，容易出現很多bug問題，下面是一個比較典型的容易遇到的報錯。
The job input file "Job-1.inp" has been submitted for analysis.
Job Job-1: Analysis Input File Processor aborted due to errors.
Error in job Job-1: Analysis Input File Processor exited with an error.
檢查工作路徑下的log文件（Job-1.log），提示無法識別“ifort”，這就意味著這個命令在系統路徑中無法被識別。實際上，“ifort”適用于編譯Fortran文件的，如果無法使用則無法完成編譯。
log文件內容
為了解決這個問題，在系統文件中找到“ifort.exe”，并將其添加到環境變量path中，然后重啟軟件即可解決。
添加路徑
下載地址：Abaqus用戶子程序專題
展開 abaqus后處理插件—場變量結果標注 ￥45
abaqus后處理界面中探針功能附帶的標記樣式非常丑陋，基于abaqus的試圖注釋功能進行二次開發，形成了場變量標注插件，方便快速的標注關心區域的應力應變等結果。
插件介紹：
按鈕介紹
從左至右依次是：標記按鈕、隱藏標記按鈕、恢復顯示按鈕、刪除按鈕
示意動畫
使用方法：
1） Probe查看節點結果，并勾選需要標記的節點項；
2）點擊工具欄中的標記按鈕，進行標記。
特點
1） 標記速度快，即使在單元數目達到百萬級及以上的模型中，標記速度仍無明顯延遲；
2）所有標記注釋均在試圖注釋功能界面里，有利于對美觀度有更高要求者進一步修改美化。
展開 Abaqus后處理-云圖變量含義（部分）
結果中分量說明：
S11、S22、S33指各軸正應力；
S12指作用于XZ平面(與“2”，即Y軸垂直的平面)內，沿1方向剪應力；
S13指作用于XY平面內，沿1方向剪應力；
S23指作用于XY平面內，沿2方向剪應力。
若為柱坐標，S12、S13、S23分別指：由徑向向環向的剪應力、由徑向向軸向的剪應力、由環向向軸向的剪應力。
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模擬多道次變形的變量繼承方法
一、引言
使用ABAQUS進行多道次加工時，往往牽扯道次之間變量的繼承（如晶粒尺寸、累積損傷等），這對多道次變形模擬結果的準確性有較大的影響。本文以VUHARD子程序及簡單的熱壓縮模型為例，分享雙道次壓縮之間的晶粒尺寸的繼承方法。
ABAQUS變量解讀：教你讀懂應力/應變/損傷
有限元后處理直接與數據圖片處理、論文撰寫相關，除了典型的應力張量與應變張量外，ABAQUS還提供了大量可供使用者讀取的其他應力/應變/損傷參數，這都有助于結果的分析。今天喵星人就教你讀懂其中的應力、應變及損傷的后處理細節。
一、應力相關
根據用戶手冊及后處理分類，ABAQUS提供了三類典型的后處理變量：
1.不變量
不變量的定義是指張量在坐標旋轉下保持不變的量。這些量反映了材料內在的力學狀態，與觀察方向無關，因此常用于判斷材料的屈服、破壞或變形行為。
喵星人認為以下幾種應力相關不變量相對比較重要：
Mises：基于第四強度理論，用戶手冊定義如下：
Tresca：基于第三強度理論，用戶手冊定義如下：
Tresca equivalent stress, defined as the maximum difference between principal stresses.
Pressure：靜水壓力，注意正值為壓，負值為拉，用戶手冊定義如下：
編輯
跳轉
Third Invariant：第三應力不變量，用戶手冊定義如下：
編輯
跳轉
編輯
跳轉
這里，喵星人給出更加簡潔的定義：
TRIAX：應力三軸度，可用于評估斷裂行為，尤其適用于延性損傷模型，用戶手冊定義如下：
編輯
跳轉
2.主應力
主應力是指該點應力狀態中僅存在正應力、無剪應力的特殊方向上的應力值。
展開 ABAQUS顯式子程序調用規則及nblock變量解釋
ABAQUS在調用VUHARD子程序時，每次向子程序提供136個單元（單元數＜136則提供所有單元）進行計算，使用(do k=1,nblock → end do)計算每一個單元的相關變量。因此對于1000單元來說，一共調用8次子程序。
設置全局變量commom /globals/ kdtest，在ABAQUS每一次調用子程序之后，給其加一，統計模擬過程中的總循環數，在(do k=1,nblock → end do) 循環內部輸出變量kdtest的值。（全局變量可以不跟著k的循環而變化，用戶可以根據需求設置其在代碼中的功能），結果截圖如下：
共1000個數據，1~7各重復了136次，8重復了48次，與上述分析一致。
因此在nblock實際代表的是ABAQUS提供給子程序的材料點塊，這個塊區包含的單元數與模型單元數有關，而k則是對該材料點塊實現一個遍歷，確保每個單元都被考慮到。本研究僅針對于單核計算來講，多核模擬將在后續展開介紹。
展開 
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后綴的文件。
展開 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幫助文件。）
展開 Abaqus后處理插件——云圖場變量查詢器 ￥50
基于Python對Abaqus進行后處理的二次開發，形成了場變量查詢插件，寫論文寫報告的實用小工具，話不多說直接上圖。
其實Abaqus后處理模塊本身具有場變量查詢功能（Tools-Query-Probe values），可以滿足使用需求，但其顯示效果非常不友好，以節點（單元）編號+場變量值的形式顯示，理工科氣息十足。不得不說，Abaqus在操作體驗方面而言和Workbench還是有一些差距。
PS：常規的二次開發通過fieldvariable.getsubset方式查詢到的場變量為節點直接插值的計算結果，即節點上的值是從相鄰單元插值后未經過平均化的值，是離散的，而結果顯示云圖默認是按75%平均化處理后的結果，更接近實際。
本插件查詢到的場變量為平均化之后的結果，和云圖的顯示的值一致。
展開 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.
展開 