在進行多晶粒材料力學數值模擬時，voronoi模型被廣泛應用，目前算法也較多，有興趣的同學可參考計算機圖形學相關教材。
       就筆者個人經驗，voronoi晶粒模型的網格質量往往會對計算過程和結果產生較大影響。原因就在于：常規算法得到的voronoi圖形出現的短邊和小平面會導致模型整體單元數量增加以及不良單元出現概率增大。解決這類問題的方法就是：voronoi圖形優化，去除短邊和小平面。
       本期案例為一個立方體模型，晶粒數量均為20個，對比voronoi多面體經過圖形優化和沒有優化的網格差異。

        圖1、voronoi晶粒形狀優化與網格質量對比，左圖為Cubic_1、右圖為Cubic_2，從上至下依次為幾何結構圖、網格剖分圖，零厚度內聚力單元晶界圖
        從上圖可知，經過形狀優化后的voronoi晶粒模型網格質量得到了有效提升，網格數量明顯下降，不良單元得到了消除。不難看出，經過優化的模型，計算時長和結果精度將會優于未優化的模型。
       注： 本項目目前不接受答疑，僅提供工程協作，協作范圍：各類立方體狀、柱狀、球狀voronoi晶粒模型的構建，包含但不限于：常規模型，優化模型，晶粒長大模型、包含亞晶粒的多尺度晶粒模型等。
如需協作請提前將個人需求整理成word，私信留言，我會及時回復。

更多文章請關注微信公眾號：ABAQUS二次開發