傳統損傷模型對于單元的尺寸十分敏感，不同單元尺寸會導致有限元模型精度出現明顯偏差。針對該問題，梯度損傷（Gradient-damage）模型的概念被提了出來。 本文詳細介紹了如何將梯度損傷模型應用于4節點平面單元，并在有限元模型中進行模擬。 ABAQUS提供了UEL（user defined element）給使用者進行開發。筆者利用UEL開發4節點平面單元，其邊界條件如下圖所示。其中，節點

利用ABAQUS自定義單元子程序，既可以開發新的單元，同時也可以定義新的材料本構模型。本文以損傷模型簡單應用于4節點平面單元為案例，介紹ABAQUS UEL的開發和使用。 如上圖所示，該單元包含4個節點，每個節點有兩個自由度，分別在水平（X）和垂直（Y）方向運動。節點1的兩個自由度被固定，節點4的水平自由度被固定，節點2的垂直自由度被固定。節點3和節點4在垂直方向上向上運動，位移為0.1mm

基于matlab的實現四節點板單元剛度矩陣求解，振動模態分析。可自主輸入材料參數，板單元長寬厚尺寸。輸出振型結果，輸出多少階可自主設置。最后以可視化的形式展示，程序已調通，可直接運行。

作者模擬得到單調拉伸以及循環加載下材料的宏觀應力應變響應為 微觀響應結果為 基于兩類疲勞指示因子，作者通過線性外推得到了基于模擬的壽命預測結果: 基于作者提供的思路和參數，對黃永剛原始程序進行修改，考慮背應力效應，并進行簡單的數值驗證 1，建立包含200晶粒的二維多晶模型（0.1*0.03mm），并使用四節點平面應變單元進行網格劃分

為什么要導出單元剛度矩陣 在學習有限元方法時，我們會需要編寫程序計算結構的單元剛度矩陣。此外，當我們需要做有限元軟件二次開發時，我們也需要驗證所做的開發是否正確。為了驗證程序正確性，我們可以從商業有限元軟件中導出單元剛度矩陣來驗證程序的計算結果。下面簡單介紹從ansys軟件中導出平面四邊形四節點單元的單元剛度矩陣。 平面四邊形四節點單元示例 如圖所示

然而，晶粒取向對宏觀應力-應變曲線的影響很小 基于作者提供的思路和參數，對黃永剛原始程序進行修改，考慮背應力效應，并進行簡單的數值驗證 1，建立包含200晶粒的二維多晶模型（0.1*0.03mm），并使用四節點平面應變單元進行網格劃分，如下圖 2，施加正弦形式的循環拉壓的位移載荷（1%），引力比為-1 3，模擬結果如下： 第一個滑移系統的背應力：

密實橡膠材料用不可壓的四節點平面應變 Herrmann 單元 （MSC. Marc 中單元類型 80），海綿橡膠一般采用Foam 模型，可選用任意四邊形平面應變單元 （MSC.Marc 中單元類型11），金屬骨架采用四節點平面應變全積分單元（MSC. Marc 中單元類型27）。

2階8節點減縮積分平面應變單元子程序UELMAT源代碼及計算算例

導讀：大家好，我是SimPC博士，主要從事工程結構抗震及減隔震研究，玻璃成型熱工設備流動及傳熱研究，玻璃材料力學性能研究。精通有限元等數值算法的實現，有限元軟件二次開發，數據處理，偏微分方程求解，優化算法，GUI界面開發等。有多項科研成果，其中SCI論文4篇，EI3篇，專利2篇。 近日我注冊并認證了技術鄰專欄，將在技術鄰官網和App給平臺用戶帶來Matlab有限元編程、復雜函數擬合和偏/常微分方程求解

（6）劃分網格 ①土體：四邊形網格；結構劃分技術；CPE4（四節點平面應變單元）；隧道內部網格進行局部加密。 ②管片：周長方向120個網格，CPE4I（平面應變四節點非協調單元）。 ③鋼管：周長方向120個網格，CPE4I。 （7）修改輸入文件，設置場變量 執行模型、編輯關鍵字、model-1；在相應位置添加下列語句，以此實現隧道內土體在reduce分析步模量進行軟化。