二维Comsol的Voronoi边界设置与多边形骨料、纤维骨料分析方法

二维comsol的Voronoi可设置方形边界圆形边界椭圆边界等等。 可选择条带过渡界面厚度。 需要ABAQUS2020及以上版本AUTOCAD2020及以上版本 以上两软件进行辅助生成。 另二维多边形骨料纤维骨料等均可采用此方法。在复合材料建模领域二维Voronoi结构生成是个挺有意思的活。最近折腾COMSOL的时候发现结合Python脚本可以玩出各种花样。比如搞个混凝土骨料模型边界形状随便换——方的、圆的、椭圆的想怎么整就怎么整。先上段核心代码看看怎么生成带边界的Voronoiimport numpy as np from scipy.spatial import Voronoi def generate_voronoi(boundary_typesquare, size10): points np.random.rand(50, 2) * size # 添加边界控制点 if boundary_type circle: theta np.linspace(0, 2*np.pi, 36) boundary size/2 * np.column_stack([np.cos(theta), np.sin(theta)]) else: # 默认方形边界 boundary [[0,0], [size,0], [size,size], [0,size]] all_points np.vstack([points, boundary]) return Voronoi(all_points)这段代码的骚操作在于动态添加边界控制点。比如处理圆形边界时先生成36个圆周上的点强制加入Voronoi生成过程这样生成的晶格会自动贴合圆形轮廓。参数size控制生成区域尺寸换成椭圆的话改改坐标变换就行。界面厚度控制是个技术活这里有个取巧的方法——给每个Voronoi边添加缓冲区def add_interface(vor, thickness0.1): new_ridges [] for ridge in vor.ridge_vertices: if -1 not in ridge: # 排除无限边 start vor.vertices[ridge[0]] end vor.vertices[ridge[1]] # 计算垂直方向偏移 normal np.array([end[1]-start[1], start[0]-end[0]]) normal thickness * normal / np.linalg.norm(normal) new_ridges.append([start normal, end normal]) new_ridges.append([start - normal, end - normal]) return new_ridges这招相当于给原始边线上下各偏移一个厚度值形成双线效果。参数thickness控制过渡层宽度调这个值就能改变界面带的尺寸。实际跑起来会发现当thickness超过相邻晶格间距的一半时界面带就会开始融合这个临界点需要注意。二维comsol的Voronoi可设置方形边界圆形边界椭圆边界等等。 可选择条带过渡界面厚度。 需要ABAQUS2020及以上版本AUTOCAD2020及以上版本 以上两软件进行辅助生成。 另二维多边形骨料纤维骨料等均可采用此方法。ABAQUS老司机们应该知道2020版开始支持直接导入SVG路径。把生成的Voronoi结构用matplotlib导出SVGimport matplotlib.pyplot as plt from svg.path import parse_path def export_to_svg(vor, filename): fig plt.figure() ax fig.add_subplot() voronoi_plot_2d(vor, axax, show_verticesFalse) ax.set_aspect(equal) plt.savefig(filename, formatsvg) plt.close() # 转换路径为ABAQUS可识别的格式 with open(filename) as f: paths parse_path(f.read()) return paths导出的路径在ABAQUS里可以直接当草图用比传统方法省事不少。AutoCAD 2020的Python API也挺给力处理复杂边界时可以用它来修修补补。比如遇到椭圆边界需要精细调整时调用acad.model.AddEllipse()方法直接创建参数化椭圆比手动描点精准多了。纤维骨料的生成其实可以看作是Voronoi的变种——把随机点换成线状分布就行。改改生成点的策略def generate_fiber_points(num20, length10): points [] for _ in range(num): angle np.random.rand() * np.pi x np.linspace(0, length, 50) y x * np.tan(angle) points.extend(np.column_stack([x, y])) return np.array(points)这种点阵生成的Voronoi结构会自动拉长形成纤维状分布。配合之前的边界控制方法轻松实现各向异性材料的建模。实际测试中发现当纤维角度分布范围控制在±30度时力学性能的各向异性表现最明显。整个流程跑下来从生成到仿线大概需要这么几步Python生成Voronoi基础结构AutoCAD处理特殊边界ABAQUS进行网格划分和力学分析COMSOL做多物理场耦合可选注意不同软件版本间的兼容性问题特别是AutoCAD 2020和ABAQUS 2020的API有较大改动。建议先在Jupyter里调试好生成算法再移植到完整脚本中。遇到边界点不闭合的情况可以试试在AutoCAD里用PEDIT命令做下合并比代码处理省心。