多晶体模型采用三维Voronoi算法生成，试件尺寸为150×150×300mm棱柱模型，对晶格指定五种不同材料，实现晶格间的差异性。

对试件进行力学模拟，下侧为固定边界，限制z方向的位移，上表面通过给定位移的方式实现轴压模拟。

不同晶格的杨氏模量如下图所示。

计算所得的应力及位移云图。

COMSOL晶体建模可采用CAD Voronoi 3D插件进行。