ABAQUS uel子程序的案例
基于ABAQUS的UEL子程序定義4節點平面應變等參單元的剛度問題
摘要：
采用基于ABAQUS的UEL子程序開發4節點平面應變等參單元，采用雙線性形函數，4點高斯積分，本構關系為線彈性各向同性材料，得到的單元剛度矩陣和ABABUS自帶的CPE4單元的單元剛度矩陣（剛度矩陣輸出方式為*element matrix output, elset= ALLE, stiffness=yes, OUTPUT FILE=USER DEFINED）不同；對比ANSYS的單元剛度矩陣，結果顯示兩者也不相同。問題出在哪里呢？本文檔將對此問題進行回答。
本文可以作為ABAQUS高級子程序UEL的入門級教程，做UEL的應該關注下！
基于ABAQUS的UEL子程序定義4節點平面應變等參單元的剛度問題（技術鄰 藍牙）.pdf
展開 Abaqus子程序代碼分享
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.
展開 abaqus最簡單的uel子程序自定義單元例子 ￥3
通過abaqus的uel子程序定義了2節點單元的剛度矩陣，直接運行分析可以得到節點應力、位移等參數，可以快速了解uel的組織架構。
運行方法：abaqus創建job——來源——點擊輸入文件——加載inp文件——編輯作業——通用模塊——加載用戶子程序for文件——運行
【JY】 ABAQUS子程序UEL的有限元原理與應用
三、計算結果對比
在相同邊界條件、載荷條件下，二維三角形單元的UEL子程序計算結果與有限元ABAQUS計算結果如下表3及圖4-5所示，從圖中可看出，子程序的位移計算結果與abaqus中CPS3單元位移計算結果一致。
ABAQUS UMAT UEL子程序咨詢解疑 ￥400
ABAQUS UMAT UEL子程序咨詢解疑
【全源碼】MATLAB相場脆性斷裂模擬代碼（AT1/AT2）【附對應文獻公式說明】 ￥1000
現有資源的門檻：網上的開源代碼多為Fortran編寫的Abaqus UEL/UMAT子程序，調試極其困難，且相當于“黑盒”，難以直觀理解算法邏輯。
驗證的難題：寫出了代碼，但不知道結果對不對，缺乏權威的Benchmark（基準）進行對比。
現在以帶偏心孔的缺口板為例，說明我編寫的MATLAB代碼準確性。幾何和邊界條件如下圖所示：
2. 驗證準確性
本案例提供了一套基于 MATLAB 編寫的相場斷裂有限元代碼，完整實現了 AT1 和 AT2 兩種經典的相場損傷模型。
本代碼的核心價值在于“精準驗證”： 代碼邏輯清晰、注釋詳細，更重要的是，選取了較為復雜的的斷裂力學算例（帶偏心孔的缺口板），將本代碼的計算結果與 吳建營教授（相場領域權威）發布的Abaqus UEL子程序計算結果 進行了逐點對比，驗證了代碼準確性。
位移-反力曲線：兩條曲線幾乎完全重合。
裂紋路徑（Crack Path）：裂紋擴展形態結果高度一致。
這是吳建營教授發布的Abaqus UEL子程序計算的裂紋路徑：
這是MATLAB代碼計算的裂紋路徑：
3. 代碼特點
雙模型支持：代碼內通過參數設置，可切換 AT1 和 AT2 模型。
交替迭代算法：采用了魯棒性較好的位移場與相場交替求解策略，收斂性好，適合初學者學習。
展開 發布ABAQUS 的 UEL算例 3， 4 ，9節點INP和uel子程序文件，提供入門級別的幫助 ￥50
發布ABAQUS 的 UEL算例 3， 4 ，9節點INP和用戶子程序
ABAQUS最強大的地方有兩點（1）強大的非線性求解能力；（2）友好的二次開發子程序接口。
對于用戶子程序而言，終極模式就是UEL或者VUEL（幾乎除了求解器，其它的都自己DIY），下面提供幾個UEL入門級別的算例，最早發布于2010年 simwe論壇，遷移至此。
3節點梁（非本人原創，特此聲明） UEL of 3node beam.rar
以下為付費內容
四節點平面單元（本人原創，下同）
3節點平面單元
9節點彈性程序
for 和inp文件如下
展開 Abaqus用戶子程序UEL資料及源碼整理 (暫停0903) ￥1
資料已刪，此貼暫停，請勿購買
木網殼雙節點考慮轉動受軸力影響的自定義連接單元-UCEL ￥200
（1）基于ABAQUS UEL子程序接口創建了此單元；
（2）基于WE-P理論分析模型編寫的子程序關系。
有限元理論基礎及Abaqus內部實現方式研究系列20： UEL用戶子程序開發步驟
iSolver子程序的接口完全按照Abaqus的標準實現，而Abaqus的子程序接口在近幾年內已經基本不再變化了，同樣的，雖然iSolver在不斷發展，但iSolver子程序接口將維持不變，所有在iSolver上編寫的算法子程序都只要維護自己的算法部分就行，而不是維護整個有限元求解的整個過程。
前面第八、九篇介紹了UMAT用戶自定義材料的開發，這里將介紹UEL用戶自定義單元的開發，本文首先簡單的討論了UEL的一般含義，并詳細的介紹了基于Fortran和Matlab兩種方式的UEL的開發步驟，對比發現開發步驟基本相同，同時采用Matlab更加高效和靈活。
展開 UEL子程序開發——梁單元
圖 3 變形協調條件
子程序講解
編制矩陣操作函數
為了程序的易讀性，我們可以將經常使用的的計算部分編寫成一個子程序，使用時直接call一下即可，比如在UMAT中，我們會常常用到Mises應力公式，所以我們可以將Mises的求解編寫為一個子程序，在UMAT中調用即可，大大減少了代碼量，更能提升易讀性。
本程序中，我編寫了矩陣置零函數，用于數組的初始化，大家在自己的程序中編制矩陣相乘函數和矩陣求逆函數，大大提升編寫代碼的效率。
【7月2日項目懸賞】
【單號6447】
預算范圍：1000
使用軟件：abaqus
需求描述：邊坡樁板墻地震動力響應建模
立即搶單
【單號6446】
預算范圍：2000
使用軟件：ABAQUS+VS+IVF(FORTRAN語言)
需求描述：鋼材在反復荷載作用下的損傷累積力學模型，需要用Fortran語言二次開發編譯子程序，有參考文獻，有一段循環塑性的程序和相應的驗證結果。
立即搶單
【單號6418】
預算范圍：1000
試用軟件：abaqus
需求描述：要應用自定義單元（abaqus UEL子程序）解決巖土工程中的問題，可是如何在UEL單元里施加重力？其地應力平衡應該如何實現？
立即搶單
【單號6416】
預算范圍：2000
使用軟件：COMSOL
需求描述：重復論文里的模型，COMSOL應力-滲流-損傷耦合，定義損傷變量，做出巖石壓縮-滲流過程的損傷變化，QQ554769224價格可商量
立即搶單
下載“技術鄰”APP，或點擊鏈接查看所有派單： http://www.yqgqt.org.cn/requirement/more
注：目前手機不支持搶單，如需搶單請用電腦瀏覽器打開該網址，或打開技術鄰首頁，在首頁右側點擊“懸賞”進行搶單。
客服微信：jishulink999
展開 材料脆性斷裂有限元模擬的UEL子程序實現方法
UEL (User Element Subroutine) 子程序允許用戶自定義單元的切線剛度矩陣及節點力向量，在實現相場斷裂模型時具有靈活性與便利性。UMAT (User Material Subroutine)子程序可以供用戶自定義材料的本構模型，同時可彌補UEL子程序無法可視化的缺陷。本文中的計算結合UEL子程序與UMAT子程序，采用雙層模型進行交互計算材料的斷裂過程。
【6月28日項目懸賞】
【單號6416】
預算范圍：2000
使用軟件：COMSOL
需求描述：重復論文里的模型，COMSOL應力-滲流-損傷耦合，定義損傷變量，做出巖石壓縮-滲流過程的損傷變化
立即搶單
【單號6341】
預算范圍：1500
使用軟件：abaqus
需求描述：八分之一石墨烯包裹碳化硅顆粒增強鋁基單胞模型，加對稱邊界條件，石墨烯和鋁之間加一層cohesive單元觀察脫粘，加周期性邊界條件，加載荷得到應力應變云圖和曲線，基體jc損傷，石墨烯橫觀各項同性，碳化硅當做線彈性的（八分之一模型不能加周期性邊界條件就一個正方體包一個球）
立即搶單
【單號6418】
預算范圍：1000
試用軟件：abaqus
需求描述：要應用自定義單元（abaqus UEL子程序）解決巖土工程中的問題，可是如何在UEL單元里施加重力？其地應力平衡應該如何實現？
立即搶單
【單號6428】
預算范圍：500
使用軟件：workbench
需求描述：求取一對斜齒輪嚙合的時變嚙和剛度，包括仿真和最后提取結果計算剛度的過程 工況：恒定轉速，但是負載是不規則的變化曲線
立即搶單
下載“技術鄰”APP，或點擊鏈接查看所有派單： http://www.yqgqt.org.cn/requirement/more
注：目前手機不支持搶單，如需搶單請用電腦瀏覽器打開該網址，或打開技術鄰首頁，在首頁右側點擊“懸賞”進行搶單。
客服微信：jishulink999
展開 二維三維彈塑性斷裂相場------uel子程序 ￥69.9
包含二維三維彈塑性斷裂相場程序，原始代碼，以及對應的程序公式，job文件，計算收斂性較好，運算穩定，可使用二維的三角形，四邊形單元，以及三維四面體，六面體單元的程序計算案例的結果：
二維：
二維斷裂相場分布：
三維斷裂相場模型：
