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

結合作者的理論（尤其是分段線性化和應力驅動的求解思路）我們可以把獨立的vpsc子程序編寫進abaqus里面，為了避免復雜的雅可比推導，以及適用各種復雜的變形工況，推薦使用abaqus的顯式求解器，即vumat程序以下展示一個使用vpsc-鎂合金本構模型，模擬包含1個單元，單元包好100個晶粒在RD方向壓縮20%的模擬效果（原始模型參數取自vpsc官方案例,為了減少計算時間使用高應變率進行計算，

這里以編寫VUMAT進行三點彎曲模擬為例子）值得注意的是當我們把單晶的本構模型與有限元求解器結合時，不同晶粒之間在有限元的意義上同時滿足應力平衡和應變協調。隨著計算器處理數據能力的提升，晶體塑性有限元方案CP-FEM逐漸替代了對應的均勻化方案，但在尺度跨越時，均勻化方案依然可以成為不同尺度下溝通的橋梁。

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

但實際金屬材料并不是“均勻黑箱”：晶粒取向、滑移系激活、織構演化都會影響局部塑性變形，尤其在薄壁管壓潰這類大變形、強局部化問題中，微觀結構可能對吸能行為產生重要影響。

重劃分 Remeshing + 狀態變量映射（最通用） 當網格畸變到閾值，換一張“干凈網格”，把舊網格的歷史狀態（取向、硬化、位錯密度等）映射到新網格繼續算——這是很多晶體塑性/微觀模擬里最常用的工程化路線。

單位統一使用m，pa，運行結束后的應力分布應變分布晶粒旋轉角度分布使用歐拉角作為輸出，變形結束后的極圖分布（初始隨機取向）：相應結果和linux平臺下的進行了詳細對比，結果保持一致。

考慮到這一點，VUMAT是用矢量化接口構建的。在每次增量開始時，應力和狀態變量數據以二維陣列的形式傳入。陣列中的每一列都包含與材料的積分點相關的信息。當使用多個處理器執行模擬時，可以將分析數據拆分為塊并獨立求解。因此，在子程序的編寫中保留了矢量化，以便可以實現最佳的處理器并行化。使用ABAQUS/explicit不需要相同的過程，因為時間增量足夠小。

</p><p>3 自定義材料笨狗或單元笨狗（umat或uel）</p><p>主要是通過自定義損傷材料笨狗或者單元笨狗模擬裂紋擴展，主要包括：usdfld、vusdfld、umat、vumat、uel等等。