晶體塑性耦合
關注創建者：匿名 創建時間：2026-01-04
晶體塑性耦合的視頻教程
4.DAMASK晶體塑性有限元平臺案例實戰教程——多晶體晶體塑性分析——Mg(HCP)
DAMASK是德國馬普所推出的功能強大的開源晶體塑性有限元計算平臺，其可以實現多尺度、多場耦合的晶體塑性分析。本次課程旨在讓大家快速入門軟件使用，尤其是其譜求解器的使用和前后處理。
￥99 5分鐘 121播放
查看
基于Abaqus軟件的晶體塑性有限元分析(2)-基于UMAT的晶體塑性有限元程序
為了幫助大家在學習晶體塑性有限元分析過程中少犯錯和少走彎路，系列課程基于Abaqus軟件進行晶體塑性有限元分析(2)-基于UMAT的晶體塑性有限元程序。
￥125 26分鐘 800播放
查看
2.DAMASK晶體塑性有限元平臺案例實戰教程——2D多晶體鋁合金晶體塑性分析(FCC)
DAMASK是德國馬普所推出的功能強大的開源晶體塑性有限元計算平臺，其可以實現多尺度、多場耦合的晶體塑性分析。本次課程旨在讓大家快速入門軟件使用，尤其是其譜求解器的使用和前后處理。
￥99 18分鐘 89播放
查看
晶體塑性耦合的實例教程
通過和黃永剛晶體塑性模型進行耦合可以實現介觀尺度下，多晶材料的完整彈-塑-損傷力學行為分析，并且相比與其他損傷模型耦合方式而言，耦合相場法物理含義更加清晰，數值實現格式簡介，處理雅可比矩陣方便且易于收斂。因此逐漸受到介觀尺度分析材料損傷分析學者的青睞。
這里通過耦合常用的晶體塑性模型（黃-umat（修改取向到狀態變量））和斷裂相場方法，剛度和應力退化使用二次退化函數形式。模擬包含200個晶粒的多晶模型，使用平面應變簡化，拉伸變形30%，模擬拉伸過程中裂紋的產生和發展，其中斷裂總能量包含彈性變形能和塑性耗散功兩部分，模擬結果如下
初始多晶模型：
網格劃分（CPE4網格）：
相場分布（0：材料完好，1：材料完全失效）：
退化程度分布：
歐拉角（phi）分布：
可以看到耦合相場的晶體塑性模型具有潛在的預測裂紋萌生和發展的能力，其準確程度取決于斷裂能參數的選擇，與更精細的實驗對比，如原位的ebsd拉伸將成為良好的校核手段。這可能成為介觀尺度下斷裂力學的應用提供良好的參考
展開 經典文章推薦
《Necking behavior of AA 6022-T4 based on the crystal
plasticity and damage models
是最經典的耦合晶體塑性理論和連續損傷的文章之一，損傷力學有兩種主要方法。第一種是Gurson提出的基于微觀力學的損傷模型。在基于微觀力學的方法中，損傷演化通過孔隙成核、生長和聚結來描述。對空穴成核和生長進行了建模，必須使用實驗數據確定相關系數。另一種方法是連續損傷力學（CDM）。在CDM框架中，使用應力、壓力、溫度和應力三軸性確定斷裂應變。在這些研究之后，提出了許多改進的模型，以包括洛德角和各向異性損傷的影響，
作者在研究中使用的損傷模型基于連續損傷力學（CDM）。然而，通過結合CPFEM可以預測孔隙的萌生、生長和聚結行為。此外，材料因損傷而弱化用于描述頸縮后承載能力的突然下降，通過顯式時間積分方案進行了分析，這為通過CPFEM預測頸縮行為提供了可能性。然而，沒有預測頸縮形狀和載荷位移曲線。為了準確預測頸縮和載荷位移曲線，使用隱式時間積分方案進行了分析，可以獲得更合理的載荷位移曲線。此外，還進行了實驗，并與分析結果進行了比較。最后，新提出了四種不同的帶系數校準的損傷模型，并提出了一種最能描述頸縮行為的模型。
作者使用的四類連續損傷模型理論如下
（1）最大塑性應變損傷模型：該模型將損傷定義為當主塑性應變大于某一臨界值時開始和累積的損傷。此模型寫為：
ε1f.ini是損傷萌生塑性應變值，ε1f.ini是最大塑性應變值，D是損傷因子，M是損傷指數（通常取值大于1.0有利于流動應力平滑過渡）
（2）等效塑性應變損傷模型：該模型將損傷定義為當等效塑性應變大于某一臨界值時開始和累積的損傷。
展開 參考應變率：ε0
當滿足下列條件時，損傷初始化準則得以滿足：
等效塑性應變認為與應力三軸度和應變率相關聯。
θ^是無量綱溫度，表示為：
其中，θ是當前溫度，θ-melt是熔化溫度，θ_transition是指轉變溫度，在該溫度或低于該溫度時，損傷應變εpl_D的表達式不存在溫度依賴性。材料參數必須在轉變溫度或低于轉變溫度時測量。
損傷的發展可以公式化為:
公式中分母表示單元失效對應的Johnson-cook等效塑性應變，公式為：
分子表示為等效塑性應變增量，公式為：
公式中可以看到，損傷隨著塑性應變的增大不斷累積，直至材料的失效，通過損傷變量進一步與晶體材料的屈服面或者彈性性能的退化可以實現材料彈-塑-損傷的耦合模擬，當不對其進行耦合時，可以用來判斷材料的失效狀態與相關參數的關系。
參考文獻：《Crystal plasticity finite element modeling and simulation of diamond cutting of polycrystalline copper》編寫對應的材料子程序。在huang晶體塑性程序的基礎上，調用johnson-cookd損傷函數，編寫過程中，需要自定義響應的狀態變量，如等效塑性應變，等效塑性應變率，損傷變量，以及是否進行損傷單元的刪除分析。其中等效塑性應變增量的計算，通過滑移系統的分切應力與對應滑移系統剪切應變的乘積絕對值之后與等效應力的比值獲得。并最終實現損傷的表征，采用umat子程序進行編寫。
展開 DAMASK解決損傷與晶體塑性耦合問題 ￥100
基于DAMASK譜求解器算法的FFT解決金屬材料的損傷演化及裂紋起源及擴展問題
溫度場通過初始溫度以及塑性產熱計算，同時忽視局部的熱傳導，準靜態加載速率下的泰勒-昆尼系數η為0.0，1000 s加載速率下為0.95?1及以上(塑性功轉化為熱的比例)
通過經典的熱激活模型，將溫度效應引入流動方程，并考慮溫度對剛度的退化
位錯密度模型演化遵循經典的KM模型，同時考慮位錯之間的相互作用，即考慮了位錯的產生和湮滅，以及湮滅半徑與溫度的關系。因此有利于由實驗進行對照分析。
損傷基于經典的JC損傷，并等效的對應力進行退化
拉伸模型
網格劃分（每個單元表示一個單獨取向的晶粒，即初始的取向不同）
局部斷裂時溫度場分布（初始293K，假設taylor-Q系數為0.95）
局部斷裂時局部位錯密度分布（僅考慮統計儲存位錯密度）
局部斷裂時損傷分布
局部斷裂時等效塑性應變分布
局部斷裂時mises等效應力分布
展開 
晶體塑性耦合的相關專題、標簽、搜索
晶體塑性耦合的最新內容
原始文獻：《An elasto-viscoplastic formulation based on fast Fourier transforms for the prediction of micromechanical fields in polycrystalline materials》
DOI：10.1016/j.ijplas.2011.12.005
在計算微觀力學領域，如何高效預測多晶體內部的異質應力場量一直是核心難題
文章名稱《A three dimensional (3D) thermo-elasto-viscoplastic constitutive model for FCC polycrystals》
DOI：10.1016/j.ijplas.2015.04.001
在鋁合金、鎂合金等輕質材料成形過程中，溫度往往不是一個可以忽略的因素。尤其是在溫成形條件下，材料的流動應力、硬化能力、延性、應變率敏感性以及彈性回復都會發生明顯變化
文章名稱《Concurrent multi-scale crush simulations with a crystal plasticity model》
DOI：10.1016/j.tws.2011.12.019
在汽車防撞梁、吸能盒和薄壁管結構中，壓潰吸能能力直接影響結構安全性。傳統有限元分析通常采用各向同性塑性模型，通過宏觀應力–應變曲線描述材料響應。但實際金屬材料并不是“均勻黑箱
文章名稱《Simulation of polycrystal deformation with grain and grain boundary effects》
DOI：10.1016/j.ijplas.2011.03.001
做多晶材料模擬時，我們經常會遇到一個很現實的問題：晶粒尺寸明明會顯著影響強度，但在普通晶體塑性有限元模型里，這個效應并不會自然出現。
傳統 CP-FEM 可以很好地描述晶粒取向
文章推薦：《Reduced-order representations of crystallographic texture for application to surrogate modelling of austenitic stainless steel》
晶體塑性有限元（CPFE）模型在預測多晶材料宏觀性能與微觀晶體學織構的相互作用中扮演著核心角色 。然而，極其龐大的計算成本成為了將其推廣至宏觀工程部件
晶體塑性：構建Dream3D pipeline用于將EBSD模型制作成Abaqus可執行文件
案例實操
用于生成模型的Dream3D pipeline文件，只需要你設置EBSD數據的路徑和導出路徑即可，可以直接生成abaqus的晶體塑性模型，提供原始文件！
包含老版本Dream3D 6.5的管道文件，并且根據官方的使用說明文件已經成功移植到最新版Dream3D 7.4版本了。
晶體塑性模擬中的大變形網格重劃分4個月前
參考文獻《Large-deformation crystal plasticity simulation of microstructure and microtexture evolution through adaptive remeshing》
在我們進行大變形晶體塑性時，做到后期，最常見的“翻車點”不是本構收斂性問題，而是網格畸變：單元被壓扁/拉長后，數值誤差會明顯放大，輕則結果不準，
參考文獻：《A straightforward 3D polycrystal plasticity finite element method for dynamic/static recrystallization simulation》
文章doi：10.1016/j.jmst.2024.09.005
在這個文章中，作者提出了一種直接在 CPFEM 中實現 DRX/SRX 的方法，以位錯密度為核心變量
文獻一：《Gradient plasticity in gradient nano-grained metals》
文獻二：《Grain rotations during uniaxial deformation of gradient nano-grained metals using crystal plasticity finite element simulations》
推薦理由
初識彈塑性與晶體塑性——面向新手8個月前
文檔以ppt形式展示，從簡單走向深入，講解了彈塑性與晶體塑性在實現時的基本邏輯，內容展示：
