https://www.yqgqt.org.cn/content/post/1266640 第二十七篇：Abaqus內部計算和顯示的應變。

以往廣泛應用的數值方案通常是：先把 積分點的數據外推到節點，再用線性形函數求梯度，然而這類方案只能用特定單元（如 C3D8），對自適應網格、復雜接觸不友好。

組建剛度矩陣K，abaqus自己處理 2. 載荷列陣F，abaqus自己處理 3. x=K^-1*F，所有節點的位移 abaqus自己處理 4. 根據x，求每個單元的應變增量值 abaqus自己處理 5. 根據單元應變和單元剛度矩陣，求應力。

求解熱場 提交 T-*.inp + DFLUX.for（或 user=UFLUX.for）作業，得到 T-*.odb。 檢查代表性點溫度–時間曲線、峰值溫度、等溫線是否合理。

1 說明 該本構完全從文檔《Writing User Subroutines with Abaqus》中摘抄而來，采用Fortran77格式編寫 2 理論文檔 需要考慮熱膨脹（熱應變）和材料參數隨溫度變化。

相比于 Subspace 子空間求解法，分塊子空間法對網格單元質量要求，計算機內存以及硬盤空間要求都適中，求解速度較快、精度高。

二、剪切自鎖 在小變形線彈性分析中，在求出節點位移向量的解后，需要進一步算出應變場；非線性分析中，在一個增量步迭代得到位移向量解后，也需要算出相關應變值，再代入本構數據中查詢本構點，進而構造下一個增量步迭代所需要的初始切線剛度矩陣。然而，與我們通常的印象不同，這里計算應力應變值，是在積分點上計算的，也就是是將積分點的坐標值代入應力應變的公式，而不是直接求節點的應力應變。

https://www.yqgqt.org.cn/content/post/1266640 第二十七篇：Abaqus內部計算和顯示的應變。

https://www.yqgqt.org.cn/content/post/1266640 第二十七篇：Abaqus內部計算和顯示的應變。

<p><strong>概述</strong>：結合HHT時程積分法，推導了ABAQUS <strong>靜/動力隱式算法</strong>中的UEL關鍵矩陣表達式，并將公式應用到自編CPS4/CPE4、C3D8 BBAR和C3D20用戶自定義單元中，計算結果均與ABAQUS自帶單元保持一致。其中，靜力計算中，關鍵矩陣AMARTX和RHS等可直接按照剛度矩陣和方程右端不平衡力輸出。