晶體有限元
關注創建者：淵魚 創建時間：2023-06-28
晶體有限元的視頻教程
基于Abaqus軟件的晶體塑性有限元分析(2)-基于UMAT的晶體塑性有限元程序
為了幫助大家在學習晶體塑性有限元分析過程中少犯錯和少走彎路，系列課程基于Abaqus軟件進行晶體塑性有限元分析(2)-基于UMAT的晶體塑性有限元程序。
￥125 26分鐘 800播放
查看
基于Abaqus軟件的晶體塑性有限元分析(0)-前言
為了幫助大家在學習晶體塑性有限元分析過程中少犯錯和少走彎路，本次課程基于Abaqus軟件進行晶體塑性有限元分析。
免費 7分鐘 6016播放
查看
基于Abaqus軟件的晶體塑性有限元分析(1)-Voronoi多晶體模型生成方法
為了幫助大家在學習晶體塑性有限元分析過程中少犯錯和少走彎路，系列課程基于Abaqus軟件進行晶體塑性有限元分析(1)-Voronoi多晶體模型生成方法。
￥199 32分鐘 763播放
查看
晶體有限元的實例教程
作者：辭殤
關鍵詞：CPFEM；鈦合金；單軸拉伸；織構極圖；孿晶
晶體塑性有限元是一種結合了晶體塑性理論和有限元方法的數值模擬技術?。這種方法考慮了晶體材料的各向異性、滑移系統的開動和相互作用、以及變形過程中的硬化效應。它主要用于分析和預測晶體材料的塑性變形行為，特別是在微觀尺度上的變形機制。
晶體塑性有限元在材料科學和工程領域有著廣泛的應用，特別是在金屬加工、航空航天、汽車制造和生物醫學等領域。通過這種技術，研究人員和工程師可以更好地理解材料的力學行為，從而開發出更輕、更強、更耐用的材料和產品。此外，晶體塑性有限元仿真還能夠考慮材料的微觀結構特征，如晶粒取向、晶界、相分布以及滑移系統的活動，從而能夠預測材料在細觀尺度上的織構演化。
利用CPFEM方法對鈦合金圓棒拉伸過程進行模擬，使用UMAT子程序以及Abaqus有限元軟件作為晶體塑性有限元分析的實現方式。并且，在一些復雜工藝條件下如切削、軋制、沖壓等，CPFEM方法同樣適用，能夠模擬材料變形過程中的非線性行為和動態響應。
在晶體塑性有限元中，首先在Abaqus中建立了單軸拉伸有限元模型如圖1所示，材料被建模為包含大量晶粒的集合體如圖2所示，每個晶粒都有其特定的晶體取向，并且每個晶粒的變形過程均考慮了滑移和孿晶的變形機制。
圖1 單軸拉伸有限元模型示意圖
圖2 單軸拉伸晶體塑性模型示意圖
通過有限元方法，可以計算出在給定拉伸載荷下，這些晶粒如何相互作用，以及它們如何隨時間變形。這種方法能夠提供關于晶體材料內部應力、應變和變形機制的詳細信息，有助于理解材料在受力時的響應，并優化材料的設計和加工過程。圖3所示為單軸拉伸過程應力云圖，圖4所示為單軸拉伸過程孿晶云圖。
展開 參考資料
Ti3Al單晶和雙相片層TiAl合金塑性行為的CPFEM模擬
Ti-6Al-4V合金納米壓痕變形與高周疲勞行為CPFEM研究
On Predicting the Channel Die Compression Behavior of HCP Magnesium AM30
兩相鈦合金拉伸力學行為的研究
密排六方金屬鎂的晶體塑性力學性能研究
HCP多晶體塑性的數值模擬
TA15鈦合金高溫變形多晶體塑性有限元模擬
γ-TiAl多晶體壓縮變形機制的晶體塑性有限元研究
純鈦單道次ECAP變形織構演化的細觀有限元模擬
純鈦晶體塑性力學性能研究
純鈦塑性變形行為的晶體塑性有限元模擬
純鈦壓縮變形下的晶體塑性有限元分析
考慮滑移與孿晶的鎂塑性本構研究
鈦合金雙態組織高溫拉伸行為的晶體塑性有限元研究
展開 晶體塑性有限元仿真入門(3)--開源代碼平臺EVOCD
晶體塑性有限元初學者較為熟知的兩個工具Huang's UMAT以及DAMASK平臺，這篇文章介紹另外一個晶體塑性有限元方法(CPFEM)的開源代碼平臺EVOCD，講解如何使用這些開源代碼進行材料的塑性變形模擬以及模擬變形過程中晶體取向的變化(織構)。
圖1 EVOCD的CPFEM流程圖
(E.B. Marin, Center for Advanced Vehicular Systems, Mississippi State University)
我們在網上搜索晶體塑性的關鍵字''CPFEM''時，會發現搜索引擎的網頁排名第一是馬普所(MPI, 大名鼎鼎的DAMASK就是他們團隊的成果)的研究成果，其次是密西西比州立大學先進車輛系統中心(Center for Advanced Vehicular Systems, Mississippi State University)的開源代碼平臺EVOCD，第三是基于Huang的晶體塑性有限元方法，由此可見EVOCD在晶體塑性有限元方法中的重要性。
圖2 CPFEM搜索結果
(從上到下分別是馬普所 (dierk-raabe.com) 、密西西比州立大學 (msstate.edu) 、哈佛大學 (columbia.edu) 的相關研究成果)
國內的晶體塑性有限元初學者，最主要的還是使用Huang's UMAT以及DAMASK平臺，而對密西西比州立大學的開源代碼平臺EVOCD不太常用。這篇文章將講解該平臺的使用方法以及如何使用該平臺進行晶體塑性有限元變形模擬。
展開 還有如果有同學需要damask和neper、ABAQUS等，請參考我的另一篇帖子——基于linux虛擬機的晶體塑性有限元計算平臺（已安裝damask和neper等）》
http://www.yqgqt.org.cn/content/post/1287159
【VMware】如何打開已有虛擬機：https://jingyan.baidu.com/article/6181c3e0c44ebd152ef15332.html
最后附上測試視頻：
neper測試.mp4
晶體塑性有限元仿真入門(5)—歐拉角與晶體取向.pdf
晶體塑性有限元仿真入門(5)—歐拉角與晶體取向
備注：網頁排版有亂碼，建議下載附件pdf查看
晶體取向是材料學科中的重要分支，當晶粒發生擇優取向時，則導致材料性能（力學，物理和化學性能）的各向異性。各向異性會造成材料實際應用中的各種問題，如鋁合金典型的制耳現象，再如取向硅鋼中存在Goss織構時，有利于其磁學性能。在基礎研究領域，織構的形成與演變是基本的科學問題。在工業應用領域，通過織構的設計和控制可以提高材料的性能。隨著近年來EBSD和XRD等表征技術的發展，各種SCI期刊的發文都已離不開對材料晶體學取向的分析。這篇文章介紹晶體塑性有限元仿真過程中的歐拉角與晶體取向。
圖1 塑性變形過程導致的材料各向異性
全文包括以下幾個部分：
1) 材料晶體結構
2) EBSD工作原理
3) 晶體取向分析
4) 晶體塑性材料模型
5) 織構演變結果
6) 參考資料
7) 附錄
材料晶體結構
在晶體學中，晶體結構是對晶體材料中原子、離子或分子有序排列的描述。有序結構由組成粒子的內在性質產生，形成沿物質三維空間的主要方向重復的對稱模式，如圖2所示。
圖2 高分辨率透射電子顯微鏡圖片的鐵晶體，完美單晶的二維示意圖
構成這種重復圖案的材料中最小的一組粒子是結構的晶胞。晶胞完全反映了整個晶體的對稱性和結構，這是通過晶胞沿其主軸重復平移而建立的。平移向量定義布拉維點陣的節點，不同的晶體內部原子排列稱為具有不同的晶格結構。各種晶格結構可以歸納為七大晶系，各種晶系分別與十四種空間格（稱為Bravais晶格）相對應，如圖3所示。
展開 
晶體有限元的相關專題、標簽、搜索
晶體有限元的最新內容
考慮熱源的瞬態熱傳導有限元求解器3小時前
關鍵詞：熱源，瞬態，熱傳導，有限元求解器，三角形單元，自研
在《瞬態熱傳導有限元求解器開發》一文中，我們介紹了自研的二維瞬態熱傳導求解器。
當時那個控制方程沒有考慮熱源，邊界條件中只涉及溫度、熱流、對流。然而在很多問題中，熱源才是最關鍵的邊界條件，比如電發熱、化學反應生熱。
熱源的處理
熱源是體熱，相對應的熱流是面熱。兩者處理方式類似，都是根據單位熱功率值和幾何尺寸計算熱功率，然后加到控制方程矩陣的右側
原始文獻：《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
在計算微觀力學領域，如何高效預測多晶體內部的異質應力場量一直是核心難題
對于高溫合金、鋁合金薄壁件、微尺度構件等問題，如果材料存在明顯織構或晶粒尺度效應，將晶體塑性與結構有限元耦合，能夠提供比傳統本構更豐富的物理信息。
我們可以將我之前推文提到的umat-taylor模型轉化為vumat子程序，進一步使用晶體塑性模型模擬大變形結構尺度材料變形行為。
第一層是我們熟悉的晶粒尺度模擬，也就是基于有限元的晶體塑性計算。它負責求解每個晶粒、每個單元中的應力、應變、滑移量和位錯密度演化。第二層是介觀尺度模擬，用來處理普通 CP-FEM 很難直接描述的部分：位錯在晶粒內部的重新分布、由位錯堆積產生的背應力，以及位錯穿過晶界時受到的阻礙。
這篇文章里，最值得關注的是它對晶界的處理。
文章推薦：《Reduced-order representations of crystallographic texture for application to surrogate modelling of austenitic stainless steel》
晶體塑性有限元（CPFE）模型在預測多晶材料宏觀性能與微觀晶體學織構的相互作用中扮演著核心角色 。
關鍵詞：CFD，有限元，對流項，繞流，迎風格式，湍流模型
在《流體有限元求解器開發-不可壓定常流動模型》一文中，我們介紹了考慮對流項的不可壓流動求解器的實現。
然而正如所預料的那樣，一旦流速高一些，或者粘性小一些，仿真結果就容易發散，收斂性成為一大難題。
為了解決這個問題，CFD大神們想出了各種手段，有的嚴格按照理論去處理盡力彌合。有的則主打靈感修正，問就是人工粘性、人工擴散、人工穩定
<h3><strong>【版權聲明與技術存證】關于某型“巷道超前支架”結構有限元分析報告的公開撤回聲明</strong></h3><p><strong>一、 成果歸屬與授權撤回</strong></p><p>本文發布內容為本人針對某型巷道超前支架所做的有限元分析（FEA）階段性成果。</p><p><strong>合作背景說明：</strong> > 合作方：<strong>西安某礦業學科背景高校相關研究團隊
隨后可進行網格劃分，設置邊界條件，并執行晶體結構的有限元仿真分析。
可劃分網格，并進行晶體結構有限元仿真分析。
[圖片]
