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

<p>本資源包含一份 PDF 文檔和可直接編譯運行的 Fortran UMAT 代碼，具體內容為：</p><p>Chaboche硬化本構模型 + 隱式積分 + 徑向返回</p><p>完整公式推導 + Fortran 源碼直接編譯</p><p>任意個數背應力分量 + 解析一致切線模量</p><p>PDF 包含規范化的本構方程、隱式積分、徑向返回與一致切線模量推導，可供初學者學習。

<p>因為要仿真混凝土破壞實驗，考慮用abaqus里面的CDP模型，查閱了相關資料進行了理論總結，并根據理論編寫計算程序。</p><p>ABAQUS中CDP 模型中采用的是混凝土在單軸受力狀態下的應力和非彈性應變，非彈性應變根據混凝土的單軸應力-應變曲線換算。

DTS自動連接到他專屬的虛擬桌面，里面已經<strong>預裝了CATIA、Abaqus、HyperMesh等工具——全部在云端，本地無需安裝任何大型軟件。</strong></p><p><strong>上午9:00，協同設計。</strong>項目經理在DTS共享會話中發起緊急評審：供應商剛更新了翼肋的細節設計，需要結構團隊確認。小李加入會話，與遠在成都的同事同時查看模型。

相關做法完美的集中到damask3.0版本里面，然而需要指出的是：DAMASK/譜方法更偏向規則網格與RVE范式，而工程里經常需要：任意幾何與復雜邊界（非周期、接觸、局部細化等），以及不同工藝路徑（多道次、換向、局部約束），Abaqus CPFEM（UMAT/VUMAT）在這些方面更“通用”，所以把“remesh + 狀態變量映射”做成一套工作流，就能把大變形晶體塑性更穩地推進到更高壓縮/更大應變階段

對于損傷的判斷，纖維采用Hashin準則： 樹脂采用最大應力： 仿真效果 應力： 損傷分布： 編輯 跳轉 應力應變曲線結果： （《基于嵌入式約束的機織復合材料細觀建模與分析

問題是在UAMT/VUMAT里面是很難做這種偏導的，包括迭代數值計算，不是完全不能，而是寫出來大概率各種報錯，還不好調試找原因。在子程序里面，最穩妥的就是寫加減乘除。 那時候寫彈塑性本構，對我理解子程序以及ABAQUS邏輯，起到了非常重要的作用。我的體會是，學寫子程序，應該先寫彈性，接著就寫彈塑性，這樣才能打好基礎。像我當時屬于是回頭補課。

（2）ABAQUS可以直接導出節點上的應力、位移、應變等數據，我們只要在表格中，利用基礎的EXCEL公式，就可以嵌入自己想分析的參數，比如嵌入一個強度準則，計算失效因子，進行渲染。 （3）方便不同工況的比較。對于同一個模型，以往我們計算了不同的工況，只能截取兩個圖在那對比。有了這個工具，我們可以通表格操作，直接計算出兩個工況的參數差異，然后渲染出來。結果量化了，放到論文里面更具說服力。

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

由此可以反映出該處的應力集中分布情況，然后根據疲勞實驗，驗正計算結果是否合理。結果表明，輪轂中輪輞和輪輻所承受的應力也許不是對稱應力循環[26]。 重慶大學周渝慶[27]利用CAE軟件ABAQUS，通過計算機模擬，精確地控制了輪輻與輪輞的厚度，并且通過彎曲疲勞實驗，確定了輪轂的最大載荷，以及整個輪轂的重量，從而實現了輪轂的輕量化[28][30]。