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

? 基于黃umat梯度結構晶粒變形模擬 案例實操 1，建立包含896個晶粒的梯度多晶模型 2，對多晶模型賦予對應的材料屬性 3，X0方向固定，施加X1方向50%工程應變的拉伸載荷 4，保留晶界形狀，使用CPE3單元 5，提交與后處理材料數據 梯度晶粒幾何模型 模型載荷示意圖 不同時刻材料的對數應變分布

并分析其他可能的影響，或者使用類似的研究思路，使用更加物理的本構模型如位錯密度模型等進行對比研究進行簡單修改兩個模型實現類似的效果：為構造典型梯度結構，使用了隨機尺寸的結構文獻一模型效果：晶粒尺寸分布：不同晶粒的初始強度分布：不同晶粒的飽和強度分布：不同晶粒的變形過程累計剪切分布：不同晶粒的變形過程應力分布：另一個模型效果一致

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

多尺度模型子程序考慮微觀晶粒與宏觀單元的相互作用：1) 宏觀單元由多個晶粒組成，每個晶粒的相變行為獨立計算；2) 單元整體響應為各晶粒響應的加權平均，可模擬晶粒取向對宏觀行為的影響（如案例中 Element 1 及其不同晶粒的應變差異）。3、 案例介紹和結果對比1.

采用CPFEM模擬了面心立方結構(FCC)、體心立方結構(BCC)和密排六方結構(HCP)的單晶、多晶及多相材料受到外部載荷時的力學響應。基于滑移原理的晶體變形理論，隨著變形的進行各滑移系統的臨界剪應力都會增大，CPFEM將捕捉到材料的力學響應(應力-應變曲線)。這些應力-應變數據有助于從根本上理解晶粒尺度下金屬變形的性質。

基于黃永剛umat的FCC多晶在簡單剪切變形下的有限元分析1，使用neper軟件生成包含500個晶粒的幾何模型2，賦予相應的材料參數（基于腳本完成材料的批量賦值）3，施加相應的邊界條件（X0方向完全固定，X1方向施加Z1方向0.15mm的位移）4，結果與后處理5，模型文件見附錄6，提供材料屬性賦予腳本

研究中使用的材料是316L不銹鋼，它是面心立方（FCC）；因此，滑移可能發生在12個滑移系統上（｛111｝＜110＞）只用一組滑移系開動的雙晶模型三維模型：模擬得到的尺寸效應效果需要注意的是當前模型不會考慮考慮變形過程中晶粒變形的影響（體素網格變形過程中與晶界距離保持不變）如果向結構施加大的載荷，則晶粒有可能發生顯著程度的變形。

（一）taylor均勻化方案 Taylor均勻化模型在abaqus實現時，通常一個積分點包含多個晶粒，其中積分點響應由晶粒響應的體積平均響應給出，積分內的各個晶粒的變形梯度均等于宏觀變形梯度。通過變形梯度計算各個晶粒的應力響應，晶粒彼此之間不存在應力平衡關系，但應變協調，計算得到各個晶粒的得到柯西應力，之后得到體平均柯西應力，更新積分點響應。

對于高溫合金、鋁合金薄壁件、微尺度構件等問題，如果材料存在明顯織構或晶粒尺度效應，將晶體塑性與結構有限元耦合，能夠提供比傳統本構更豐富的物理信息。我們可以將我之前推文提到的umat-taylor模型轉化為vumat子程序，進一步使用晶體塑性模型模擬大變形結構尺度材料變形行為。

它模擬了骨料在外部應變和應力作用下的塑性變形。VPSC 基于滑移和孿晶的物理剪切機制，并考慮了晶粒相互作用效應。除了提供宏觀應力-應變響應外，它還解釋了單個晶粒的硬化、重新定向和形狀變化。因此，它預測了與塑性成形相關的硬化和織構的演變。模擬程序可應用于金屬、金屬間化合物和地質聚集體的變形。

基于黃umat探究Mg-Cu雙相材料簡單拉伸下的變形行為 案例實操 1，建立包含500個晶粒的多晶模型，模型尺寸0.6*0.3*0.05（mm） 2，對晶粒編號1-250賦予Cu的屬性（參數來自于黃畢業論文）251-500賦予AZ31材料的屬性，考慮三組滑移系和一組拉伸孿晶系 3，X0方向固定，施加X1方向的25%工程應變的單向拉伸載荷 4，指定對應的單元類型

模擬包含200個晶粒的多晶模型，使用平面應變簡化，拉伸變形30%，模擬拉伸過程中裂紋的產生和發展，其中斷裂總能量包含彈性變形能和塑性耗散功兩部分，模擬結果如下初始多晶模型：網格劃分（CPE4網格）：相場分布（0：材料完好，1：材料完全失效）：退化程度分布：歐拉角（phi）分布：可以看到耦合相場的晶體塑性模型具有潛在的預測裂紋萌生和發展的能力

此外，晶體塑性有限元仿真還能夠考慮材料的微觀結構特征，如晶粒取向、晶界、相分布以及滑移系統的活動，從而能夠預測材料在細觀尺度上的織構演化。利用CPFEM方法對鈦合金圓棒拉伸過程進行模擬，使用UMAT子程序以及Abaqus有限元軟件作為晶體塑性有限元分析的實現方式。并且，在一些復雜工藝條件下如切削、軋制、沖壓等，CPFEM方法同樣適用，能夠模擬材料變形過程中的非線性行為和動態響應。

，基于原始的唯象晶體塑性模型進行修改，將初始屈服，硬化模量，飽和強度，以及率相關系數構造為晶粒尺寸的函數，實現建立具有尺寸效應的多晶本構模型，這對目前金屬梯度結構介觀尺度下力學性能的表征具有一定的啟發性</span></p><p><span style="background-color: rgb(193, 230, 198);">文獻一的研究使用Voronoi鑲嵌方法構建梯度納米晶結構，使用的本構模型如下

然而，如果模擬規模變小，例如在專注于納米壓痕（Zaafarani et al.，20082006）和微柱壓縮（Raabe，Ma和Roters，2007a）的研究中，則局部模型可能由于無法描述尺寸效應而不足，較小晶粒尺寸的強化效應是由于晶界附近非均勻塑性變形的體積分數較高。

使用文章提到的策略，嘗試進行數值顯示，首先在umat隱式中進行實現，并在后續中修改為vumat即可。實現策略驗證使用包含200個晶粒的二維模型拉伸驗證。

基于黃umat和cohesive實現疲勞過程中的沿晶斷裂分析 案例實操 1，建立包含50個晶粒的多晶模型，晶界處使用cohesive單元 2，分別對晶界cohesive單元和晶粒實體單元進行材料屬性賦值 3，X0方向固定，施加X1方向的疲勞載荷 4，指定對應的單元類型 5，提交與后處理材料數據 晶粒幾何模型 晶界幾何模型

同時需要注意ABAQUS/標準采用Jaumann應力率，而ABAQUS/explicit采用Green–Naghdi應力率Ω表示變形梯度極分解對應的旋轉部分，在ABAQUS/標準中，材料被視為基于固定的全局坐標系。增量旋轉在每次增量開始時傳遞給UMAT。這個數組dR是應力和應變數組在上一次增量結束和當前增量開始之間旋轉的量。

圖4 金屬增材制造殘余應力的形成機制 固有應變法固有應變法最早由日本學者上田（Ueda）提出，廣泛應用于大型焊接結構的扭曲與殘余應力預測。由于可以快速實現大型復雜零件的殘余應變與扭曲變形預測，固有應變法目前已成為零件級增材制造模擬的主流方法，并且已經被多種商用增材制造模擬軟件所集成。金屬增材制造模擬中固有應變的獲取主要有2種方法：微觀尺度模擬和標準件變形標定。

參數與文獻參數保持一致，初始位錯密度為1.2e16mm^-2幾何模型：包含0.3mm*0.3mm包含67個晶粒的二維模型，包含90000個四邊形網格：邊界條件：施加X方向變形的周期性邊界條件，應變為5%變形結束后滑移系統總位錯密度分布情況：變形結束后應力分布情況：變形結束后累計taylor剪切為：可以看到模型很好的符合了作者文獻的模擬效果