利用ABAQUS自定義單元子程序，既可以開發新的單元，同時也可以定義新的材料本構模型。本文以損傷模型簡單應用于4節點平面單元為案例，介紹ABAQUS UEL的開發和使用。

如上圖所示，該單元包含4個節點，每個節點有兩個自由度，分別在水平（X）和垂直（Y）方向運動。節點1的兩個自由度被固定，節點4的水平自由度被固定，節點2的垂直自由度被固定。節點3和節點4在垂直方向上向上運動，位移為0.1mm。該正方形單元的邊長為100mm。在input文件里，坐標表示為，

定義節點組合與邊界條件為，

為了讓模型收斂性更好，采用quasi-newton 求解器。時間步設置為，

在文件夾中通過Powershell提交job和子程序，

單個單元的變形為，

采用不同的 ??

 ，在后處理中得到損傷因子的變化，

相對應的力-時間關系為，

對于多個單元的情況，比如9單元組成的模型，

具體介紹見知乎：ABAQUS UEL - 損傷材料本構簡單應用于4節點平面單元 - 知乎 (zhihu.com)

相對應的UEL代碼和input文件在付費內容中，