RVE的案例
IJP：從RVE到組件的跨尺度預測
為了克服多尺度模型在單個尺度上預測的不足，來自凝固加工國家重點實驗室的Xinxin Sun等人將元胞自動機晶體塑性有限元法和人工神經網絡巧妙地結合起來，建立了從代表體積元(RVE)到組分的跨尺度整體預測模型。CACPFEM模型充分耦合了非均勻變形和微觀結構演變，如動態再結晶(DRX)，用于解釋RVE的響應。為了反映應變率、溫度、微觀結構和變形模式對響應的依賴性，對RVE施加了大量恒定和變化的加載路徑。所有的響應(包括力學和微觀結構的響應)形成了一個巨大的數據庫，在此基礎上，通過訓練、驗證、測試和循環優化，建立了具有Marquardt-Levenberg (M-L)算法的反向傳播(BP) ANN模型。人工神經網絡模型的輸出設置為微觀結構演變(包括DRX體積分數和平均晶粒尺寸)和取決于加載路徑和微觀結構的J2-JBOY3樂隊本構模型的動態變化的宏觀尺度參數，然后應用于有限元模型以預測部件的響應。因此，建立了一座橋梁來連接RVE和組件的響應。反過來，部件局部區域的變形歷史也可以應用于RVE，以進一步研究微尺度變形機制和微結構演化。利用跨尺度模型，得到了反映各向異性、拉壓不對稱性、應變率、溫度、微觀結構和變形模式依賴性的結果。它得益于基于物理的CACPFEM、依賴于變形條件和微結構演化的J2-JBOY3樂隊本構模型、優化的ANN模型以及它們的創新組合。優化策略保證了跨尺度預測的準確性。跨尺度模型在旋轉坯料單軸壓縮中的應用以及在新形狀坯料的屈服面預測和鍛造過程中的推廣，表明了該模型的跨尺度預測能力。
圖1 跨尺度模型的框架
跨尺度模型的框架，如圖1所示，該框架基于四個部分和三種優化方法的組合。
CACPFEM用于數據樣本生成:CACPFEM模型充分耦合了非均勻變形、力學響應和微結構演化之間的相互作用。
展開 Ansys Workbench | 材料微觀結構：四種 RVE 的均質化分析
選擇隨機單向纖維作為代表性體積元（RVE）。設置纖維體積分數為0.4，纖維直徑為50μm。創建幾何模型（圖1），并使用默認設置生成網格。
4. 創建一個恒定材料，并求解工程常數。工程常數匯總如圖2所示。可以觀察到，纖維方向上的整體楊氏模量 E1 比 E2 和 E3 大100%以上。這是因為纖維的楊氏模量高于基體，從而增強了縱向剛度。這種微觀結構的典型例子是木材和一些復合材料。
圖1. 隨機單向纖維的 RVE
圖2. 隨機單向纖維結構材料的工程常數
案例2：體心立方結構（金屬）
5. 按照案例1的相同步驟操作。為顆粒定義各向同性材料屬性（E=25000MPa, ν=0.3），并為基體定義各向同性材料屬性（E=18000 MPa, ν=0.3）。
6. 定義體心立方結構 RVE（圖3）。顆粒尺寸設為1nm。生成網格。這種微觀結構是金屬的典型代表。
圖3. 體心立方結構的 RVE
7. 求解工程常數。工程常數概覽如圖 4 所示。由于顆粒在三個方向上的分布相同，因此得到的宏觀尺度材料是各向同性的，例如鋼和金。
圖4. 體心立方結構材料的工程常數
案例3：晶格結構（金剛石）
8. 按照案例1的相同步驟操作。為梁定義各向同性材料屬性（E=100MPa, ν=0.3）。注意，此材料屬性僅為示例，并不代表金剛石的真實值。
9. 定義金剛石形狀的晶格 RVE，屬性如圖 5 所示。生成網格。顧名思義，金剛石具有這種微觀結構。
圖5. 金剛石晶格結構的 RVE
10. 求解工程常數。工程常數概覽如圖 6 所示。這種金剛石晶格微觀結構也導致各向同性的材料行為。
圖6.
展開 周期性 RVE 幾何建模：復合材料單胞自動化生成插件
復合材料多尺度力學仿真中，代表性體積單元（RVE）的幾何建模與網格劃分是前處理階段的主要工作之一。受周期性邊界條件的約束，纖維在模型邊界處的切割精度直接影響后續網格匹配。當纖維端面與基體表面未能完全共面時，往往產生微小幾何階躍，導致節點投影誤差。這些問題在手動腳本處理時出錯的概率較高。
針對上述情況，基于Abaqus環境開發了Periodic RVE Generator插件，對纖維生成、布爾切削及空間排布算法進行了重新編寫，以提升建模穩定性與操作效率。以下就工具的主要算法邏輯和使用方式作簡要說明。
圖 1. 插件GUI界面（可輸入基體尺寸，纖維直徑，長度，纖維體積分數，短纖維/連續纖維，設置纖維方向）
一、纖維拓撲形態的定義
為適應不同分析層次的需求，插件將纖維的幾何拓撲與空間取向解耦。通過Fibre Form選項可切換短纖維與連續纖維兩種模式。
當選用連續纖維時，程序調用超限切削邏輯：先在計算基體尺寸后，使纖維初始生成時超出邊界，隨后通過全局布爾運算切除外部多余幾何體。這一處理方式使得所有纖維端面與基體表面具備一致的平齊度，避免了切割面階差對周期性網格對齊造成的影響。
圖 2. 連續纖維（左圖）和短纖維（右圖）周期性單胞
二、纖維空間分布算法
插件內置了兩種空間拓撲分布方式：
正交約束排布：控制纖維沿指定的X、Y或Z方向對齊，適用于單向板類RVE的構建；
三維隨機分布（Random 3D）：采用球面投影與隨機變量正弦變換生成取向向量，保證空間方向無統計偏置。通過干涉檢測算法，在較高體積分數條件下仍能保持一定的生成成功率。
圖 3.
展開 Abaqus混凝土周期性邊界代表體單元插件：Random Sphere RVE 3D (Mesh) - AbyssFish ￥698
插件介紹
Random Sphere RVE 3D (Mesh) - AbyssFish 插件可在Abaqus生成三維具備周期性邊界條件(Periodic Boundary Conditions, PBC)的隨機球體骨料及骨料-水泥界面過渡區(Interfacial Transition Zone, ITZ)模型。即采用周期性代表性體積單元法(Periodic Representative Volume Element，PRVE)，以代表體積單元(Representative Volume Element，RVE)或稱為表征單元體(Representative Elemental Volume, REV)微觀結構的計算來準確地模擬和預測混凝土材料的宏觀行為。插件采用體素網格方式，通過背景網格將砂漿、骨料、ITZ劃分為三個集（Set），并對單元映射三種空材料。
插件支持設置長方體部件的長度(Length)、寬度(Width)、高度(Height)，以及在網格劃分中單元的尺寸(Element size)。可設置生成球體的最小粒徑(D_min)及最大粒徑(D_max)，即球體尺寸的分布范圍，球體占整個長方體試件的比例(Ratio)，界面過渡區的厚度(ITZ)，以及超時終止參數(Time)。
模型可分為砂漿基體、界面層、球體骨料三相材料。
插件生成的模型均滿足周期性分布邊界條件。
可對每個集(Set) 批量插入嵌入0厚度cohesive粘結單元（注：需要自行添加，本插件不具備此功能）。
?
展開 
精沖鋼微觀組織對其力學性能和精沖性能影響的多尺度模擬研究
模型整體尺寸為100μm×100μm，網格劃分時將球狀碳化物等效為此RVE模型中的一個網格。
⑵基于實際微觀組織的RVE模型。
圖3 C15E熱軋碳鋼鐵素體-珠光體混合微觀組織的RVE模型建立過程
未球化退火熱軋鋼的微觀組織由片層狀珠光體與鐵素體基體混合而成，兩相形狀不規則，難以通過軟件直接生成RVE模型。因此需要基于材料的真實金相組織建立RVE模型，即利用軟件對真實的金相組織進行黑白兩值化處理，以分出鐵素體和片狀珠光體兩相，在建模軟件中按兩相的位置信息構建RVE模型。典型的未球化退火態RVE模型建模過程如圖3所示。考慮到兩相組織的尺度為10μm左右，RVE模型的整體尺寸為100μm×100μm。
結果與討論
微觀組織對拉伸、剪切力學性能的影響
通過對純微觀的RVE模型施加拉伸、剪切的邊界條件，模擬滲碳體顆粒的不同直徑、體積分數、分布特征對拉伸、剪切力學性能的影響，將模擬所得不同模型的名義應力—名義應變曲線示意圖如4、圖5所示。
圖4 滲碳體顆粒直徑、體積分數、分布狀態對拉伸性能的影響
圖5 滲碳體顆粒直徑、體積分數、分布狀態對剪切性能的影響
圖4a、4b分別反映了粒徑、體積分數的變化對拉伸性能的影響。圖5a表示在不同粒徑和體積分數下，剪切加載RVE模型至起始應變時的名義應力—應變曲線。根據模擬結果，當粒子直徑減小，粒子分數減小時，RVE模型內部發生斷裂的時間將被推遲，同時斷裂發生時對應的名義應力值將變大。
當球狀滲碳體呈帶狀分布時，即微觀組織中出現碳化物帶時，在帶中更容易發生應力和應變集中。描述碳化物帶分布特征的參數有兩個：碳化物帶中的碳化物顆粒占模型總顆粒的百分數和帶中碳化物顆粒聚集密度。圖4c、4d為含碳化物帶材料的拉伸名義應力—應變曲線，圖5b給出了碳化物帶狀分布對RVE剪切加載輸出曲線的影響。
展開 基于ansa的rve生成代碼 ￥50
<p>通過ansa軟件，寫對應的py代碼來自動生成隨機的3D的RVE模型，并自動進行均質化計算，得到剛度矩陣。代碼中已經對纖維長度，半徑，體積分數，還有基材和纖維的模量和泊松比進行了參數化，可以批量生成多種不同類型的RVE模型。</p><p><br></p><p><br></p><p><br></p><figure style="text-align: center;" class="ql-align-center">
<figure class="figure-image" contenteditable="false" data-img="https://img.jishulink.com/202504/attachment/2b7fecf82bed47cf8ae05a01df1ac618.png" style="display: inline-block;" data-regular="true">
<img src="https://img.jishulink.com/202504/attachment/2b7fecf82bed47cf8ae05a01df1ac618.png" data-mobile-src="https://img.jishulink.com/202504/attachment/2b7fecf82bed47cf8ae05a01df1ac618.png?image_process=/format,webp" data-pc-src="https://img.jishulink.com/202504/attachment/2b7fecf82bed47cf8ae05a01df1ac618.png?
展開 hypermesh怎么在復合材料RVE模型對稱面上畫分完全對稱網格問題
2、那么貌似只能不對體進行布爾運算，對每個面進行2D tetra , 然后對稱網格之后還需要不停的 equivalence, 生成一個enclosed 面網格，再用 tetra mesh 把封閉面網格生成四面體網格，這樣做的弊端是反復的equivalence真的非常麻煩，而且對于復合材料RVE模型，纖維束的關系是相互搭接連接關系是一個閉環，需要處理共節點的地方非常多。
3、嘗試分六面體，一樣的問題，在分好一個以后，分另一個掃略需要沿著已有面網格進行，掃略后又要檢查是否需要equivalence，但是整體反倒比畫四面體要來的方便。
4、由于模型具有對稱性，直接導入1/4模型，然后畫網格，對稱，再對稱。這是目前想到的可能更簡單的方法。
5、干脆纖維束和基體分別分網，之后設置個tie,這樣分網工作量小很多，也不要求共節點，后續運算可能效率低。
寫了幾條其實是提了一下問題，都沒有很好的解決。如果有懂復合材料RVE模型分網的，請指教。
微信截圖_20201125000219.png
微信截圖_20201125000312.png
微信截圖_20201125000254.png
展開 RVE，傳熱分析
有RVE建模或傳熱分析經驗或學習資源的朋友請聯系，可有償，謝謝
LS-DYNA人工智能多尺度計算技術及其在注塑成型復合材料領域的應用
在多尺度材料設計和分析方面，LS-DYNA軟件提供了RVE建模功能，其思路是數值化地重構材料樣本，這些數值化的材料樣本模型可以非常準確地代表真實材料的微觀幾何結構，我們將其稱為代表性體積單元，簡稱RVE。以纖維增強復合材料為例，如果我們知道纖維取向和體積分數的具體數值，并且可以分別測量出纖維和基體的材料屬性，那么就能為這種短纖維增強復合材料，構建對應的RVE數值模型，然后對該RVE模型開展有限元計算，以預測均質化的復合材料宏觀屬性。上圖展示了用于RVE分析的LS-DYNA關鍵字，RVE分析功能對復合材料，在材料層面上的虛擬設計和測試非常有用。
現在如果考慮更高層面，不僅是在材料樣本層面而是著眼于大尺度的復合材料部件，若要對復合材料部件開展非常準確的結構分析，那么可以考慮使用高精度的多尺度結構分析方法。多尺度結構分析方法的基本思路是，首先用有限元離散化全局的復合材料部件，然后將每個有限元分別耦合到與其局部材料微觀結構相對應的RVE模型，同時對這些RVE模型也用有限元方法進行離散化。因而在動態仿真中的每一個時間步，都可以根據纖維和基體的屬性對RVE模型開展局部有限元分析。
在RVE分析中，根據從全局有限元計算獲得的宏觀應變，預測均質化的復合材料宏觀應力和宏觀材料剛度，計算出復合材料的宏觀應力后，可以繼續完成復合材料部件的全局有限元計算，然后預測整體部件在下一個時間步內產生的結構變形，通過這種多尺度仿真能完全避免，因采用傳統的復合材料本構關系所引起的模型誤差。
因此理論來說，前文所描述的多尺度結構分析方法能實現非常高的預測精度。然而在實際中，對于大型結構而言這類高保真的多尺度結構分析仿真成本過高，對于這里顯示的模型即使采用并行計算，也需要花費一個多月的時間。
展開 LS-DYNA人工智能多尺度計算技術及其在注塑成型復合材料領域的應用
在多尺度材料設計和分析方面，LS-DYNA軟件提供了RVE建模功能，其思路是數值化地重構材料樣本，這些數值化的材料樣本模型可以非常準確地代表真實材料的微觀幾何結構，我們將其稱為代表性體積單元，簡稱RVE。以纖維增強復合材料為例，如果我們知道纖維取向和體積分數的具體數值，并且可以分別測量出纖維和基體的材料屬性，那么就能為這種短纖維增強復合材料，構建對應的RVE數值模型，然后對該RVE模型開展有限元計算，以預測均質化的復合材料宏觀屬性。上圖展示了用于RVE分析的LS-DYNA關鍵字，RVE分析功能對復合材料，在材料層面上的虛擬設計和測試非常有用。
現在如果考慮更高層面，不僅是在材料樣本層面而是著眼于大尺度的復合材料部件，若要對復合材料部件開展非常準確的結構分析，那么可以考慮使用高精度的多尺度結構分析方法。多尺度結構分析方法的基本思路是，首先用有限元離散化全局的復合材料部件，然后將每個有限元分別耦合到與其局部材料微觀結構相對應的RVE模型，同時對這些RVE模型也用有限元方法進行離散化。因而在動態仿真中的每一個時間步，都可以根據纖維和基體的屬性對RVE模型開展局部有限元分析。
在RVE分析中，根據從全局有限元計算獲得的宏觀應變，預測均質化的復合材料宏觀應力和宏觀材料剛度，計算出復合材料的宏觀應力后，可以繼續完成復合材料部件的全局有限元計算，然后預測整體部件在下一個時間步內產生的結構變形，通過這種多尺度仿真能完全避免，因采用傳統的復合材料本構關系所引起的模型誤差。
因此理論來說，前文所描述的多尺度結構分析方法能實現非常高的預測精度。然而在實際中，對于大型結構而言這類高保真的多尺度結構分析仿真成本過高，對于這里顯示的模型即使采用并行計算，也需要花費一個多月的時間。
展開 RVE使用均質化方法求平均應力應變
基于Python對二維rve計算提取EVOL 得到的總面積明顯大于實際
面積 可能是因為啥呢？
RVE模型添加載荷 ￥10
# 獲取RVE的最大邊界值
a = mdb.models[ModelName].rootAssembly
p = a.instances[PartName]
xPointList = []
yPointList = []
zPointList = []
verticleList = p.vertices
for v in verticleList:
Loc = v.pointOn[0]
xPointList.append(Loc[0])
yPointList.append(Loc[1])
zPointList.append(Loc[2])
# 獲得基體的尺寸信息
Xmax = max(xPointList)
Xmin = min(xPointList)
Ymax = max(yPointList)
Ymin = min(yPointList)
Zmax = max(zPointList)
Zmin = min(zPointList)
展開 多晶塑性任意復雜網格周期性邊界的施加
多晶材料的宏觀性能來自內部晶粒與晶界的復雜相互作用，而我們在計算中只能截取有限大小的 RVE。如果邊界處理不當，RVE 的響應會被“邊界效應”主導：例如邊界過度約束導致材料顯得過硬，或邊界過度自由導致材料顯得過軟，甚至出現非物理的應變局部化或旋轉模態。這種誤差會直接影響應力–應變曲線、各向異性參數（如 R 值）、晶粒內應變分布和損傷起裂位置等關鍵結論。
周期性邊界條件的目標是：讓 RVE 的對邊在變形上保持一致性，使得 RVE 在數值上可看作“無限周期平鋪”的材料內部單元，從而在有限計算域內盡可能逼近真實材料內部的連續性與統計代表性。尤其對于多晶結構，周期性邊界能顯著降低邊界引入的偏差，使 RVE 對晶粒數、晶粒形貌與取向統計的敏感性更可控。
周期性邊界的核心思想是“對邊協同變形”。直觀理解：如果把 RVE 沿三個方向無限重復拼接，那么相鄰兩個單元在拼接面上的位移與變形必須連續，否則材料會出現裂縫或重疊。因此，在周期性條件下，RVE 的一對相對邊界（或面）需要滿足兩點：（1）波動位移在對邊相同：也就是去掉宏觀均勻變形后，剩余的局部起伏在對邊要匹配；（2）宏觀應變通過對邊位移差來體現：對邊的位移差對應外加的平均變形（例如單軸拉伸、剪切等）。這樣，RVE 內部既允許存在晶粒尺度的不均勻變形與應力波動，又保證了邊界處與“相鄰單元”拼接時的兼容性，從而實現對無限材料內部行為的合理近似。
在很多工程實現里，周期性邊界常依賴“節點一一對應”：要求相對兩邊（兩面）具有嚴格相同的網格拓撲與節點分布。為了滿足這一點，建模時往往需要：在幾何上確保邊界嚴格對齊；在網格上強制生成周期匹配節點；甚至為滿足配對而犧牲局部網格質量。
對多晶模型而言，這帶來一個典型問題：邊界鋸齒形。
展開 ABAQUS中求解某部分單元的平均應力或平均應變 ￥10
1、參考模型：單向纖維的RVE模型；
2、腳本功能：針對指定的單元集合，在后處理中求解平均應力和平均應變。
3、應用的公式：一階均勻化計算方法。對于 RVE 模型的平均真應力和平均真應變，可通過對 RVE 內每一個單元的真應力 （真應變）取均值獲得。使用一階均勻化計算方法輸出的應力和應變適用于各種邊界條件，但需要對每個單元進行應力（應變）的輸出和計算。
案例 | 利用Digimat快速評估聚合物炸藥的材料性能
以材料微觀結構為基礎的熱力學數值預測模型，用代表體積元法（Representative Volume Element，RVE）可有效解決這個問題。采用Digimat-FE模塊中的RVE建模方法，建立了填充相、基體相和氣孔相三相有限元計算模型。
Digimat-FE 示意圖
具體操作流程
首先對材料特性進行設定，包括Binder和TATB材料特性以及TATB作為夾雜項的設定，和氣泡作為夾雜相的設定。在Digimat-FE中可以直接在geometry當中生成幾何模型PBX代表體單元模型，其中TATB體積分數設定為70%，孔隙度設定為5%，紅色為TATB，灰色為空穴，基材為粘結劑，可以從生成的下表中看出實際生成的TATB體積分數為0.692509，空穴的體積分數為0.0500096，與設定誤差非常小。
材料組成及RVE微結構設置生成
RVE單元各相體積分數比
然后利用內置前處理劃分網格。20秒鐘后，系統自動幾何網格，并劃分合適的網格。
接下來需要劃分設定載荷，FE當中支持內部設定載荷，可以設定軸向，雙軸，剪切等載荷類型。這里我們設置沿x軸拉伸。加載類型為應變加載，峰值為0.3。
然后就可以得到分析的結果，例如應力分布，位移分布云圖等必要信息：
如果需要計算它的工程常數，亦無需通過結果進行手算推導，在Digimat-FE下我們可以在載荷加載時做如下設定，選擇automatic properties evaluation，這里我們還可以快速計算計算E1，nu12，nu23，然后提交計算，我們可以獲得以下工程常數。
從建模到工程常數的提取，僅僅只需要10分鐘不到的時間就可以完成。
展開 