晶體塑性有限元子程序
關注創建者：匿名 創建時間：2022-01-18
晶體塑性有限元子程序的視頻教程
基于Abaqus軟件的晶體塑性有限元分析(2)-基于UMAT的晶體塑性有限元程序
晶體塑性有限元程序教學分為以下3個方面： 晶體塑性理論 | 晶體塑性用戶子程序(UMAT) | 用戶子程序(UMAT)材料參數的確定 Abaqus軟件進行晶體塑性有限元分析(2)-基于UMAT的晶體塑性有限元程序 關鍵字：Abaqus軟件；晶體塑性有限元；多晶體模型；周期邊界條件 Finite element analysis of crystal plasticity
￥125 26分鐘 800播放
查看
基于Huang的晶體晶體塑性有限元程序入門介紹
為了幫助晶體塑性有限元同學入門，以Huang的晶體塑性有限元程序為例，入門級介紹微觀組織模型建立、材料屬性賦予、結果導出與分析等相關內容。本課程是系列視頻，各部分均有自己的重點縱容，本次內容僅為簡單介紹。后續課程主要包含以下幾個方面： 1 晶體塑性基本介紹 2 微觀組織模型的建立 3 材料屬性的賦予 4. 子程序輸出的控制 5. 后續結果的導出與分析
￥20 20分鐘 1266播放
查看
基于Huang的晶體塑性有限元程序入門-3-晶體取向的輸入
本節將介紹如何將以歐拉角形式表達的晶體取向轉化為Prof.Huang晶體塑性有限元程序能識別的Miller表達形式。
￥30 12分鐘 1152播放
查看
晶體塑性有限元子程序的實例教程
圖4.4 邊界條件設置參數
網格單元類型
將單元類型設置成incompatibles modes(默認的可能出現錯誤zero hour glass stiffness)：
圖4.5 網格設置參數
提交運算
選擇子程序所在目錄，提交計算：
圖4.6 選擇子程序所在目錄
后處理界面
應力應變分布：
圖4.7 后處理應力應變分布(相同晶體取向參數)
圖4.8 后處理應力應變分布(隨機晶體取向參數)
晶粒隨機取向的結果與相同取向的差異較大。
inp文件源文件：
鏈接：https://pan.baidu.com/s/11UPvZHl26QpDRupopi8LXQ
提取碼：voro
以上是對晶體塑性有限元仿真的簡單講解，如果對voronoi幾何模型，晶體塑性力學，周期邊界條件等感興趣歡迎繼續關注。
展開 晶體塑性有限元仿真入門(5)—歐拉角與晶體取向.pdf
晶體塑性有限元仿真入門(5)—歐拉角與晶體取向
備注：網頁排版有亂碼，建議下載附件pdf查看
晶體取向是材料學科中的重要分支，當晶粒發生擇優取向時，則導致材料性能（力學，物理和化學性能）的各向異性。各向異性會造成材料實際應用中的各種問題，如鋁合金典型的制耳現象，再如取向硅鋼中存在Goss織構時，有利于其磁學性能。在基礎研究領域，織構的形成與演變是基本的科學問題。在工業應用領域，通過織構的設計和控制可以提高材料的性能。隨著近年來EBSD和XRD等表征技術的發展，各種SCI期刊的發文都已離不開對材料晶體學取向的分析。這篇文章介紹晶體塑性有限元仿真過程中的歐拉角與晶體取向。
圖1 塑性變形過程導致的材料各向異性
全文包括以下幾個部分：
1) 材料晶體結構
2) EBSD工作原理
3) 晶體取向分析
4) 晶體塑性材料模型
5) 織構演變結果
6) 參考資料
7) 附錄
材料晶體結構
在晶體學中，晶體結構是對晶體材料中原子、離子或分子有序排列的描述。有序結構由組成粒子的內在性質產生，形成沿物質三維空間的主要方向重復的對稱模式，如圖2所示。
圖2 高分辨率透射電子顯微鏡圖片的鐵晶體，完美單晶的二維示意圖
構成這種重復圖案的材料中最小的一組粒子是結構的晶胞。晶胞完全反映了整個晶體的對稱性和結構，這是通過晶胞沿其主軸重復平移而建立的。平移向量定義布拉維點陣的節點，不同的晶體內部原子排列稱為具有不同的晶格結構。各種晶格結構可以歸納為七大晶系，各種晶系分別與十四種空間格（稱為Bravais晶格）相對應，如圖3所示。
展開 而研究微細觀尺度的變形不均勻性是新材料開發及優選的重要準則，晶體塑性有限元方法將晶體塑性理論和有限元軟件進行了恰當的融合，成為研究細觀層次塑性變形行為的強有力工具。來自華東理工大學機械與動力工程學院的艾鑫團隊，基于Voronoi方法建立了Custom450 鋼拉伸的二維晶體塑性模型，分析了初始硬化模量、參考剪切應變率、應變率敏感系數、初始屈服應力以及飽和流動應力對材料應力——應變曲線的影響，并對晶體塑性參數進行了標定。
在文獻中，作者所建立的單晶本構模型參考了HUANGY的單晶體模型的子程序UMAT，此率相關硬化晶體塑性模型需要確定的參數包括初始硬化模量h0、初始屈服應力τ0、參考剪切應變率γ，應變率敏感系數n和飽和流動應力τs，其他參數通過計算和查找文獻獲得。基于Voronoi方法，作者在有限元軟件Abaqus中建立了Custom450材料的多晶體二維幾何模型并將本構關系嵌入軟件中，進行拉伸過程的模擬。
圖1所示是微結構模型及其網格劃分，幾何模型尺寸長度為0. 2 mm，寬度為0. 5mm，共包含100個晶粒，大小和形狀隨機，且晶粒取向隨機分布。
圖1包含100個晶粒的微結構模型及其網格劃分
圖2是邊界條件的約束情況，模型的上端面和下端面的所有節點在y方向上具有均勻的位移，左側所有節點在x方向上設置約束，使其不能橫向移動，y方向自由，在右邊界施加載荷，右側的所有節點x方向上經受同等應變載荷，而在y方向上是自由的。
圖2邊界條件示意圖
對于體心立方晶體來說，3個滑移系包括1個主滑移系和2個次滑移系。分別對包含1、2、3組滑移系開動的情形進行模擬，結果如圖3所示，只有主滑移系 { 110} < 111 >啟動時，應力——應變曲線在彈塑性區間過渡的位置存在明顯拐點，并與試驗曲線吻合良好。
展開 晶體塑性有限元有做疲勞壽命預測的嗎？
有沒有關于研究晶體塑性分析的
晶體塑性有限元子程序的相關專題、標簽、搜索
晶體塑性有限元子程序的最新內容
原始文獻：《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
在計算微觀力學領域，如何高效預測多晶體內部的異質應力場量一直是核心難題
文獻一：《Gradient plasticity in gradient nano-grained metals》
文獻二：《Grain rotations during uniaxial deformation of gradient nano-grained metals using crystal plasticity finite element simulations》
推薦理由
Damask2.0.3版本是最后一個和Abaqus有接口的版本，在Damask2.0.3的官網中已經進行了聲明，目前最新的版本已經不支持和Abaqus聯合使用。但是Damask和abaqus聯合使用仍然是學習晶體塑性有限元方法不錯的工具。我曾使用過Abaqus聯合damask平臺，這里我對使用過程中個人的一點經驗進行簡單的介紹，希望各位讀者在研究過程中少一些技術上的障礙。
1.軟件安裝
作者：辭殤
關鍵詞：CPFEM；鈦合金；單軸拉伸；織構極圖；孿晶
晶體塑性有限元是一種結合了晶體塑性理論和有限元方法的數值模擬技術?。這種方法考慮了晶體材料的各向異性、滑移系統的開動和相互作用、以及變形過程中的硬化效應。它主要用于分析和預測晶體材料的塑性變形行為，特別是在微觀尺度上的變形機制。
晶體塑性有限元在材料科學和工程領域有著廣泛的應用，特別是在金屬加工、航空航天、汽車制造和生物醫學等領域
8.0版本介紹：
晶體塑性有限元 Abaqus 三維泰森多邊形（voronoi模型）插件 V8.0
9.0相比于8.0增添三項功能如下：
1. 新增功能概覽
1.1 三維光順晶體模塊
圖1.1 三維光順晶體模塊
1.2 顯示晶體晶界后處理模塊
圖1.2 顯示晶體晶界后處理模塊
1.3 光順泡沫幾何模型生成
圖1.3 三維光順泡沫模型
abaqus后處理中顯示晶界可以是多晶塑性分析更加直觀，但abaqus未內置此功能，需要通過二次開發實現，這里分享一個插件用于實現該功能，插件源于一位法國讀博士老哥的分享，將該插件放入到abaqus plug-in中即可輕松的實現后處理晶界的顯示問題
軟件用戶界面：
得到的效果圖如下：
如果您在文章中使用了該插件，請引用該作者對應的兩篇文獻：
1，A physically-based
黃永剛院士編寫的單晶晶體塑性UMAT，主要用于在Abaqus有限元軟件中進行單晶及多晶晶體塑性變形的計算。許多科研人員通過改寫其晶體塑性UMATs，建立自己的本構子程序并發表論文。利用該本構子程序進行材料模擬的研究方向很廣泛，包括但不限于修改材料滑移、孿晶系，修改硬化方程，加入損傷，將隱式分析的UMATs修改為顯式分析的VUMATs等，理解其基本理論并將公式與代碼對應是改寫的關鍵。
晶體塑性的理論和
umat與uexternal之間的變量調用（用于晶體塑性有限元中的析出物定義或再結晶）
參考文獻：《Low-cycle fatigue life prediction of a polycrystalline nickel-base superalloy using crystal plasticity modelling approach》
在原始程序中修改流動方程，加入背應力項，引入運動硬化項，從而可以描述多晶金屬循環加載中的包辛格效應
背應力的演化遵循
擴展黃永剛原始晶體塑性程序加入AF背應力模擬金屬疲勞問題
參考文獻：《Low-cycle fatigue life prediction of a polycrystalline nickel-base superalloy using crystal plasticity modelling approach》
在原始程序中修改流動方程，加入背應力項，引入運動硬化項，從而可以描述多晶金屬循環加載中的包辛格效應
