ABAQUSUEL
關注創(chuàng)建者：關邇 創(chuàng)建時間：2021-03-25
ABAQUSUEL的視頻教程
ABAQUS UEL自定義單元子程序手把手實例研究（Fortran語言）
ABAQUS UEL子程序手把手實例研究，包括: (1) 經(jīng)典4節(jié)點平面應力單元的UEL子程序； (2) UEL子程序后處理顯示（比如單元應力云圖）； (3) UEL子程序傳參到其他子程序(比如UMAT)； (4) UEL子程序和ABAQUS自帶單元相互驗證； (5) UEL單元和ABAQUS自帶單元的混合模型；? (6) UEL單元單位轉(zhuǎn)換； (7) UEL一個單元與兩個單元的調(diào)試過程
￥200 11小時42分鐘 20517播放
查看
ABAQUS UEL/UMAT子程序綜合實例訓練營
本課程涉及的ABAQUS UEL/UMAT子程序內(nèi)容為中上難度，緊跟本人之前的ABAQUS UEL/UMAT子程序的基礎課程，點擊下面超鏈接（藍色文字）可看到該課程： ABAQUS UEL自定義單元子程序手把手實例研究（Fortran語言） 關于ABAQUS子程序的本人其他課程，可點擊下面超鏈接（藍色文字）看到： ABAQUS模擬混凝土水化熱溫度場、熱應力裂縫擴展（XFEM）
￥200 7小時26分鐘 5691播放
查看
淺談ABAQUS子程序UEL/UMAT開發(fā)及應用
適用人群：研究生、CAE仿真從業(yè)者和愛好者、ABAQUS子程序開發(fā)者 淺談ABAQUS子程序UEL/UMAT開發(fā)及應用（免費）【已結(jié)束】?直播時間：2022-09-08 19:00 ? ? ? ? ABAQUS作為一款開放度高、非線性能力強的工業(yè)軟件在行業(yè)內(nèi)深受追捧，特別是其強大的子程序模塊給眾多具有開發(fā)需求的人員提供了平臺，比如用戶自定義本構(gòu)UMAT、VUMAT，用戶自定義單元UEL、VUEL
￥10 1小時25分鐘 919播放
查看
ABAQUSUEL的實例教程
第1階模態(tài)云圖（分別為：ABAQUS、UEL_CPS4、UEL_SBFEM）
第2階模態(tài)云圖（分別為：ABAQUS、UEL_CPS4、UEL_SBFEM）
第3階模態(tài)云圖（分別為：ABAQUS、UEL_CPS4、UEL_SBFEM）
第4階模態(tài)云圖（分別為：ABAQUS、UEL_CPS4、UEL_SBFEM）
（三）Koyna混凝土壩地震動力響應分析
在壩體底部邊界輸入加速度時程（地表水平向、豎直向加速度地震動記錄見附件）。
（1）時程數(shù)據(jù)對比
壩頂-壩踵水平向相對位移時程
壩頂-壩踵豎直向相對位移時程
壩頂-壩踵水平向相對速度時程
壩頂-壩踵豎直向相對速度時程
（2）云圖對比
第4s水平向位移云圖對比（分別為：ABAQUS、UEL_CPS4、UEL_SBFEM）
第4s水平向位移云圖對比（分別為：ABAQUS、UEL_CPS4、UEL_SBFEM）
（3）峰值數(shù)據(jù)對比
統(tǒng)計三種情況的壩頂-壩踵相對數(shù)據(jù)峰值，比較表格如下。
(四)附件
附件包括包含兩個文件夾，分別為abaqus自帶單元計算文件和自編uel計算文件（for文件加密后的obj文件）。不包含sbfem的相關計算文件。
Koyna dam dynamic analysis.zip
?注：帖子不定時更新，也可能永遠不更新，慎重參考，如給您帶來誤導，深感抱歉。
展開 摘要：
采用基于ABAQUS的UEL子程序開發(fā)4節(jié)點平面應變等參單元，采用雙線性形函數(shù)，4點高斯積分，本構(gòu)關系為線彈性各向同性材料，得到的單元剛度矩陣和ABABUS自帶的CPE4單元的單元剛度矩陣（剛度矩陣輸出方式為*element matrix output, elset= ALLE, stiffness=yes, OUTPUT FILE=USER DEFINED）不同；對比ANSYS的單元剛度矩陣，結(jié)果顯示兩者也不相同。問題出在哪里呢？本文檔將對此問題進行回答。
本文可以作為ABAQUS高級子程序UEL的入門級教程，做UEL的應該關注下！
基于ABAQUS的UEL子程序定義4節(jié)點平面應變等參單元的剛度問題（技術鄰 藍牙）.pdf
展開 利用ABAQUS自定義單元子程序，既可以開發(fā)新的單元，同時也可以定義新的材料本構(gòu)模型。本文以損傷模型簡單應用于4節(jié)點平面單元為案例，介紹ABAQUS UEL的開發(fā)和使用。
如上圖所示，該單元包含4個節(jié)點，每個節(jié)點有兩個自由度，分別在水平（X）和垂直（Y）方向運動。節(jié)點1的兩個自由度被固定，節(jié)點4的水平自由度被固定，節(jié)點2的垂直自由度被固定。節(jié)點3和節(jié)點4在垂直方向上向上運動，位移為0.1mm。該正方形單元的邊長為100mm。在input文件里，坐標表示為，
定義節(jié)點組合與邊界條件為，
為了讓模型收斂性更好，采用quasi-newton 求解器。時間步設置為，
在文件夾中通過Powershell提交job和子程序，
單個單元的變形為，
采用不同的 ??
，在后處理中得到損傷因子的變化，
相對應的力-時間關系為，
對于多個單元的情況，比如9單元組成的模型，
具體介紹見知乎：ABAQUS UEL - 損傷材料本構(gòu)簡單應用于4節(jié)點平面單元 - 知乎 (zhihu.com)
相對應的UEL代碼和input文件在付費內(nèi)容中，
展開 ABAQUS提供了UEL（user defined element）給使用者進行開發(fā)。筆者利用UEL開發(fā)4節(jié)點平面單元，其邊界條件如下圖所示。其中，節(jié)點1的X、Y方向被限制住，節(jié)點2的Y方向被限制，節(jié)點4的X方向被限制，節(jié)點3、4的Y方向有豎向位移0.1mm。單元為100*100mm的二維正方形。
每個節(jié)點除了X和Y方向的位移，還帶有非局部應變（nonlocal strain）。
單個單元模型，
多個單元模型，
具體內(nèi)容可參見知乎文章：
ABAQUS UEL-梯度損傷模型應用于4節(jié)點平面單元 - 知乎 (zhihu.com)
相應的input文件和uel代碼付費可見，
展開 Journal of the Mechanics and Physics of Solids, 126: 136-150 (2019)
8、UELSGP.zip
9、ABAQUS UEL SUBROUTINE FOR PHASE FIELD FRACTURE AND HYDROGEN DIFFUSION
ABAQUS UEL子程序用于相場斷裂和氫擴散
ABAQUS user-element subroutine (UEL) for the coupled deformation – hydrogen transport – phase field fracture scheme presented in the associated paper.
Paper: E. Martínez-Pa?eda, A. Golahmar, C.F. Niordson. A phase field formulation for hydrogen assisted cracking. Computer Methods in Applied Mechanics and Engineering, 342: 742-761 (2018)
9、PhaseFieldH.zip
10、ABAQUS UEL SUBROUTINE FOR PHASE FIELD FRACTURE
ABAQUS UEL子程序用于相場斷裂
ABAQUS user-element subroutine (UEL) with a robust phase field formulation for fracture.
展開 
ABAQUSUEL的相關專題、標簽、搜索
ABAQUSUEL的最新內(nèi)容
現(xiàn)有資源的門檻：網(wǎng)上的開源代碼多為Fortran編寫的Abaqus UEL/UMAT子程序，調(diào)試極其困難，且相當于“黑盒”，難以直觀理解算法邏輯。
驗證的難題：寫出了代碼，但不知道結(jié)果對不對，缺乏權威的Benchmark（基準）進行對比。
現(xiàn)在以帶偏心孔的缺口板為例，說明我編寫的MATLAB代碼準確性。
（1）基于ABAQUS UEL子程序接口創(chuàng)建了此單元；
（2）基于WE-P理論分析模型編寫的子程序關系。
采用ABAQUS軟件通過UEL子程序進行了二維熱力耦合相場斷裂模型的求解，采用了能量分解（譜分解和球-偏分解），附件包括CAE模型（22版本）、INP文件和子程序
單個單元模型，
多個單元模型，
具體內(nèi)容可參見知乎文章：
ABAQUS UEL-梯度損傷模型應用于4節(jié)點平面單元 - 知乎 (zhihu.com)
相應的input文件和uel代碼付費可見，
本文以損傷模型簡單應用于4節(jié)點平面單元為案例，介紹ABAQUS UEL的開發(fā)和使用。
如上圖所示，該單元包含4個節(jié)點，每個節(jié)點有兩個自由度，分別在水平（X）和垂直（Y）方向運動。節(jié)點1的兩個自由度被固定，節(jié)點4的水平自由度被固定，節(jié)點2的垂直自由度被固定。節(jié)點3和節(jié)點4在垂直方向上向上運動，位移為0.1mm。該正方形單元的邊長為100mm。
概述：采用UEL接口二次開發(fā)實現(xiàn)八節(jié)點單元，考慮BBAR修正，避免體積自鎖，對標ABAQUS自帶的C3D8單元，計算的剛度矩陣、質(zhì)量矩陣和阻尼矩陣均與ABAQUS保持一致。并且采用UMAT子程序進行應力和應變數(shù)據(jù)的可視化，計算的應力應變數(shù)據(jù)同樣與ABAQUS保持一致，可視化效果同ABAQUS。以方塊的受動力簡諧荷載為例，采用上述程序，應用動力隱式計算分析步，最終計算的位移、應變等時程曲線均與ABAQUS
計算結(jié)果表明，自編UEL與ABAQUS自帶單元結(jié)果一致。
<p><strong>概述</strong>：結(jié)合HHT時程積分法，推導了ABAQUS <strong>靜/動力隱式算法</strong>中的UEL關鍵矩陣表達式，并將公式應用到自編CPS4/CPE4、C3D8 BBAR和C3D20用戶自定義單元中，計算結(jié)果均與ABAQUS自帶單元保持一致。其中，靜力計算中，關鍵矩陣AMARTX和RHS等可直接按照剛度矩陣和方程右端不平衡力輸出。動力隱式計算中，這兩者的輸出較為復雜
概述：開發(fā)了適用于靜力通用、頻率分析和動力隱式（固定增量步長和自適應增量步長均可）的三維八節(jié)點線性UEL，即ABAQUS自帶的C3D8單元，該UEL考慮了B-BAR修正，避免體積鎖死。
</p><h2>1.3 uel結(jié)果的可視化</h2><p>由于Abaqus不支持uel結(jié)果的可視化，因此我們可以利用umat來輔助進行可視化。使用Abaqus的標準單元覆蓋uel，即標準單元與uel共享節(jié)點，但是材料使用umat（設置一個很小的剛度），關鍵步驟在于使用<strong>common block</strong>將uel中的svars數(shù)據(jù)傳遞到umat中的statev中。
