不知火舞的被虐|伊人天伊人天天综合网|博洛尼亚天气|任你懆这里只有精品4|久久美日韩精品久久|掌中之物漫画免费阅读观看|0丨d老妇

巖土工程常用的Mohr-Coulomb模型和Drucker-Prager模型為理想彈塑性本構模型，MCC模型為硬化彈塑性模型，難以同時反映土體的剪切硬化和壓縮硬化，采用Mohr-Coulomb強度理論作為屈服準則，從Vermeer雙硬化模型發展而來的硬化土（HS）本構模型[2]作為一種雙屈服面硬化彈塑性本構模型，在描述軟土和較硬土的變形特性上有較好的表現[3]。

對于線性隨動硬化模型，可以選取三個數據點，保證三點處于同一直線上，對最后一組數據點進行一個特殊處理，可以選取一個很大的塑性應變值，以保證計算過程中的等效塑性應變都落在這三個數據點點，由此插值得到便滿足線性關系。

本構關系子程序采用彈塑性 - 相變耦合本構模型，總應變可分解為：1) 彈性行為：基于線彈性理論，由楊氏模量（E）和泊松比（ν）描述；2) 塑性行為：采用 J2 塑性理論，通過 Von Mises 屈服準則判斷屈服，關聯流動法則描述塑性流動；3) 相變耦合：馬氏體體積分數（tfv）通過硬化曲線調控屈服應力，塑性應變增量反哺 tfv 演化，形成 “塑性 - 相變” 雙向耦合

，基于原始的唯象晶體塑性模型進行修改，將初始屈服，硬化模量，飽和強度，以及率相關系數構造為晶粒尺寸的函數，實現建立具有尺寸效應的多晶本構模型，這對目前金屬梯度結構介觀尺度下力學性能的表征具有一定的啟發性文獻一的研究使用Voronoi鑲嵌方法構建梯度納米晶結構，使用的本構模型如下：流動方程：硬化方程為：通過假設：單晶水平上的所有抗滑移參數與局部晶粒尺寸D的平方根成反比

通過將等效塑性應變設置為零，將消除硬化記憶。如果為模型指定了背應力，這些背應力也將設置為零。如果在材料定義中包含退火行為，并且退火溫度定義為低于為金屬塑性模型指定的熔化溫度，則硬化記憶將在退火溫度下刪除，熔化溫度將嚴格用于定義硬化函數。否則，硬化記憶將在熔化溫度下自動移除。如果材料點的溫度在隨后的時間點低于退火溫度，則材料點可以再次加工硬化。

本構理論分成晶體塑性和再結晶兩部分，其中晶體塑性部分公式如下：流動方程（經典的唯象流動）：硬化方程使用的taylor位錯模型位錯密度的演化使用經典的KM方程：再結晶部分公式包含形核和晶界遷移兩部分，其中形核的理論公式是晶界遷移速度為：整體數值實現框架示意圖如下：作者以OFHC銅為研究對象，對775K和875K的熱壓縮進行了研究

文章doi：10.1016/j.ijplas.2019.09.006推薦理由：作者基于經典的K-M硬化的位錯密度模型，考慮晶界處更高的位錯塞積，提出了考慮晶界效應的位錯密度模型，探討滑移轉移對Al雙晶拉伸變形機制的影響。

黃永剛院士編寫的單晶晶體塑性UMAT，主要用于在Abaqus有限元仿真中進行單晶及多晶晶體塑性變形的計算，是許多科研工作者學習晶體塑性模擬的教學資源。可以在其基礎上對硬化模型進行修改，甚至引入損傷。UMAT主要應用于隱式分析，而對于大變形接觸問題，隱式分析往往計算效率較低。

可以根據學生的基礎適當增加UMAT和VUMAT編程的講解，Fortran基本語法等，亦可根據需求針對性的分析UMATs修改方法。課程的內容主要包括如下方面： 本構模型推導主要包含了應變分解、本構方程、硬化方程、本構時間積分方法及雅克比矩陣等推導。

（3）umat_documentation 翻譯主要對第一個pdf內容進行了翻譯（4）黃永剛程序迭代部分解釋該pdf文檔主要對umat子程序NR迭代部分進行解釋（5）晶體塑性硬化公式該pdf主要介紹常見的幾類唯象的硬化模型相關程序：（1）原始程序可參考文檔一（2）huang_umat_97對應于97版本修改后的子程序（3）考慮孿晶效應的本構模型參考

作者：辭殤關鍵詞：CPFEM；鈦合金；單軸拉伸；織構極圖；孿晶晶體塑性有限元是一種結合了晶體塑性理論和有限元方法的數值模擬技術?。這種方法考慮了晶體材料的各向異性、滑移系統的開動和相互作用、以及變形過程中的硬化效應。它主要用于分析和預測晶體材料的塑性變形行為，特別是在微觀尺度上的變形機制。

文章doi：10.1016/j.actamat.2023.119103文章通過原位EBSD實驗和晶體塑性有限元方法，研究了晶界遷移對局部變形的微觀影響機制，并進行了定量分析，作者的研究結果表明，晶界遷移不會改變激活的滑移系統，同時晶界遷移造成的應力下降與孔隙閉合造成應力下降的機理類似，并探討了GND造成的硬化，晶界遷移與材料拉伸行為的關聯機制。

.，20082006）和微柱壓縮（Raabe，Ma和Roters，2007a）的研究中，則局部模型可能由于無法描述尺寸效應而不足，較小晶粒尺寸的強化效應是由于晶界附近非均勻塑性變形的體積分數較高。文獻中有幾種基于位錯機制的解釋，如晶界前移動位錯的堆積，導致應力集中，從而增加晶界附近的滑移阻力或應變梯度，從而產生額外的位錯密度增量，從而增加滑移阻力（Evers等人，2002）。

該幾何模型由開源軟件Neper建立，約包含100個等軸晶粒。圖4 Cu多晶拉伸過程晶體塑性有限元模型示意圖 為了揭示多晶變形過程中孿晶對宏觀力學響應的影響，圖5給出了Cu多晶拉伸過程中宏觀塑性行為的演化結果。進一步分析Cu多晶和單晶變形過程中孿生機制對應變硬化行為的影響，圖6給出了Cu單晶和多晶拉伸變形過程的應變硬化率演化結果。

在文獻中，作者所建立的單晶本構模型參考了HUANGY的單晶體模型的子程序UMAT，此率相關硬化晶體塑性模型需要確定的參數包括初始硬化模量h0、初始屈服應力τ0、參考剪切應變率γ，應變率敏感系數n和飽和流動應力τs，其他參數通過計算和查找文獻獲得。

作者關于尺寸相關的晶體塑性模型的主要改進思路為：本構模型基于經典的亞彈性本構框架：其中流動方程為：硬化方程為：晶粒尺寸形態相關參數的引入：概念的引入：晶粒A滑移會導致相鄰晶粒晶界附近的應力上升，其上升程度與晶內區域距離晶界距離有關，由于晶粒之間的取向差，晶界處的應力分布受到影響。

并利用不同本構之間的建立跨尺度橋梁，為了驗證umat的能力，在同一次計算中使用兩類晶體本構（唯象與基于位錯密度的本構模型）建立包含30個晶粒的多晶模型，承受X方向5%的工程應變，其中晶粒1-15使用基于位錯密度的晶體塑性模型，16-30使用唯象的晶體塑性模型采用Voce硬化理論，變形后結果展示。

.，20082006）和微柱壓縮（Raabe，Ma和Roters，2007a）的研究中，則局部模型可能由于無法描述尺寸效應而不足，較小晶粒尺寸的強化效應是由于晶界附近非均勻塑性變形的體積分數較高。文獻中有幾種基于位錯機制的解釋，如晶界前移動位錯的堆積，導致應力集中，從而增加晶界附近的滑移阻力或應變梯度，從而產生額外的位錯密度增量，從而增加滑移阻力（Evers等人，2002）。

，在材料屈服面演化的過程中引入了MSG理論，其實現可以參考作者的另外的一篇文章《Failure and forming quality study of metallic foil blanking with different punch-die clearances》其實現方法，利用經典的taylor理論，與swift冪律硬化模型相關聯，使得材料的塑性應變硬化與材料微觀層面的位錯密度關聯。

為了拓展梯度效應引入了GND導致的應變硬化，滑移阻力的演化表示為為了更好的描述晶界，并對多晶進行離散，作者使用的的單元類型為C3D4，對應的雅可比和梯度矩陣分別為：作者模擬NiTi合金使用的材料參數為：作者的案例模型對應的數值結果根據作者的思路可以編寫對應的二維三角形單元和三維的四面體單元對應的應變梯度晶體塑性模型