3.1 第一步，在剛性板上施加-3.375mm 的位移以壓縮脊柱間隔器；第二步開始時，移除位移，使間隔器可以自由變形。 3.2 從第三步開始施加熱載荷，溫度從23.85℃ 升高到 37.85℃。在此期間，由于未發生相變，間隔器的形狀保持不變。第四步，溫度從 37.85℃ 升高到 50.85℃，由于此步中未發生主要的相變，計算再次快速收斂。

· 專業領域模塊： o Adams/Car：汽車行業專屬，快速搭建整車、懸架、轉向、制動系統模型，適配新能源汽車電池包、電驅動總成仿真。 Adams/Machinery：機械傳動專用，覆蓋齒輪、軸承、皮帶、鏈條等傳動部件的剛柔耦合仿真。 o o Adams/Flex：柔性體分析模塊，結合有限元法模擬部件彈性變形，適配精密機械、航空結構的振動與應力分析。

傳統的彈塑性模型無法準確模擬這種“又快、又熱、變形又大”的極端物理過程，而 JC 模型正是為了破解這些高能耗、高破壞性的力學難題而誕生的。 該模型的核心思想是將復雜的金屬材料行為進行“解耦”，認為材料的強度主要受到三個獨立因素的疊加影響：應變硬化、應變率（變形速度）強化和熱軟化。

結合作者的理論（尤其是分段線性化和應力驅動的求解思路）我們可以把獨立的vpsc子程序編寫進abaqus里面，為了避免復雜的雅可比推導，以及適用各種復雜的變形工況，推薦使用abaqus的顯式求解器，即vumat程序 以下展示一個使用vpsc-鎂合金本構模型，模擬包含1個單元，單元包好100個晶粒在RD方向壓縮20%的模擬效果（原始模型參數取自vpsc官方案例,為了減少計算時間使用高應變率進行計算，

結果表明，在 25–232 ℃ 范圍內，模型能夠較好描述溫度升高導致的流動應力降低和硬化行為變化；在 260 ℃ 時，模型在較小應變范圍內仍能較好預測，但高應變階段會出現偏差，這可能與動態回復、動態再結晶等高溫變形機制有關，而這些機制并未被該模型顯式考慮。 為了驗證模型的可遷移性，作者進一步進行了不同溫度下的簡單剪切模擬。

模擬的案例如下： 初始沖壓模型如下： 使用軸對稱單元可以減小模型的網格數量，顯著提高計算效率，因此模擬案例使用CAX4R單元，模型初始尺寸為R=0.015mm，H=0.0048mm，初始網格模型如下圖所示： 采用位移邊界條件加載，初始加載第一步ALE網格如下（網格會根據變形自動調整不同區域密度）： 第一步計算接觸時SSD分布： 第一步計算接觸時GND分布

后續很多孿晶模型基于此進行二次開發，因此實現該文章的數值模型對于孿晶的研究非常有幫助： 使用文章的公式，講整體算法集成到abaqus的vumat子程序相對容易，因為不需要推導一致性雅可比。但是率無關模型通常數值穩定性較差。

如果輸入的應變率曲線出現交叉（即高應變率下的應力低于低應變率下的應力），或者硬化曲線呈現負斜率（未激活損傷模塊時），求解器的材料剛度矩陣將出現非正定，導致不可控的網格畸變。此外，必須通過外推確保表格覆蓋到極高應變率（如10000 /s），以防求解器在局部高變形區發生錯誤的常數外推。

導入鞋底幾何模型（圖1）。 圖1. 鞋底幾何模型 3. 劃分網格。 使用六面體主導網格方法對整個部件進行網格劃分，設置全局網格尺寸為 3 mm。為內表面創建命名選擇，用于后續生成靜水壓流體單元。使用剖切視圖有助于選擇內表面。 4. 施加邊界條件并定義分析類型。 開啟大變形，并定義若干子步。固定底面，在頂面施加 600 N 的壓力載荷。

將插件文件夾復制至Abaqus插件目錄abaqus_plugins，例如： D:\ABAQUS2023\product\win_b64\code\python2.7\lib\abaqus_plugins 4.2. 啟動Abaqus，無需預先創建模型。