多晶體模型采用三維Voronoi算法生成，試件尺寸為150×150×300mm棱柱模型，對晶格指定五種不同材料，實現晶格間的差異性。

對試件進行力學模擬，下側為固定邊界，限制z方向的位移，上表面通過給定位移的方式實現軸壓模擬。

不同晶格的楊氏模量如下圖所示。

計算所得的應力及位移云圖。

COMSOL晶體建模可采用CAD Voronoi 3D插件進行，插件下載鏈接：

http://www.yqgqt.org.cn/post/1915603