Abaqus 泰森插件-voronoi 3D

6.0版本相对于5.0版本做了一系列的功能升级和模块新增,首先升级了基础的三维voronoi模块以及梯度晶模块,然后增加了带权重的二维voronoi和周期性边界,当然三维周期性的voronoi也增加到当前的版本中,此外6.0版本新增了网格映射模块(Mapping VoronoiEG)、用户自定义模块(UserDefine VoronoiEG)、流动两相模块(Fluid VoronoiEG ) 、泡沫结构模块(Foam Structure )以及周期性网格模块(Periodic Mesh).5.0版本介绍:晶体塑性有限元 Abaqus 三维泰森多边形（voronoi模型）插件 V5.01. 插件启动方式首先启动Abaqus界面，单击菜单栏AbqVoronoi，点击子菜单中的选项可启动不同的晶体模型生成模块，包括：Voronoi2D模块(生成二维Voronoi模型)、Voronoi3D模块(生成三维Voronoi模型)、VoronoiEG模块(生成网格型Voronoi模型)和VoronoiTool(其他辅助工具)，如图所示：图1.1 启动三维多晶模型生成插件2. Voronoi2D模块Voronoi2D模块包括：Basic Voronoi2D模块、BSpline Voronoi2D模块和Weighted Voronoi2D模块，其中Basic Voronoi2D模块可生成矩形和圆形边界的二维Voronoi晶体模型；BSpline Voronoi2D模块可生成在每个晶胞中填充B样条曲线的模型。2.1 Basic Voronoi2D模块BSpline Voronoi2D模块用户输入界面如下：图2.1 Basic Voronoi2D模块矩形边界模型用户界面图2.2 Basic Voronoi2D模块圆形边界模型用户界面2.2 BSpline Voronoi2D模块BSpline Voronoi2D模块用户输入界面如下：图2.3BSplineVoronoi2D模块用户输入界面2.3 Weighted Voronoi2D模块【新】Weighted Voronoi2D模块用户输入界面如下：图2.4WeightedVoronoi2D模块用户输入界面2.4 Periodic Voronoi2D模块【新】PeriodicVoronoi2D模块用户输入界面如下：图2.5PeriodicVoronoi2D模块用户输入界面3. Voronoi3D模块Voronoi3D模块包括：Basic Voronoi3D模块、Multiphase Voronoi3D模块、Weighted Voronoi3D模块、Gradient Voronoi3D模块、Periodic Voronoi3D和UserDefine Voronoi3D模块。3.1 Basic Voronoi3D模块【功能升级】(1) Basic Voronoi3D模块生成长方体边界模型的用户输入界面如下：图3.1Basic Voronoi3D模块长方体边界模型用户输入界面(2) Basic Voronoi3D模块生成圆柱体边界模型的用户输入界面如下：图3.2Basic Voronoi3D模块圆柱体边界模型用户输入界面(3) Basic Voronoi3D模块生成球体边界模型的用户输入界面如下：图3.3 Basic Voronoi3D模块球体边界模型用户输入界面(4) Basic Voronoi3D模块生成长方体拉伸模型的用户输入界面如下：图3.4Basic Voronoi3D模块长方体拉伸模型用户输入界面(5) Basic Voronoi3D模块生成圆柱体拉伸模型的用户输入界面如下：图3.5Basic Voronoi3D模块圆柱体拉伸模型用户输入界面3.2 Multiphase Voronoi3D模块Multiphase Voronoi3D模块的用户输入界面如下：图3.6Multiphase Voronoi3D模块用户输入界面Cell Size输入参数说明：(1) Size[radius]列指该相晶体的半径大小，记为; (2) Num列指该相晶体的数量，记为。和需保证以下关系：式中为多相晶体相数，为长方体总体积。该模块支持用户输入固定位置晶胞，点击“Set Fixed Grain”按钮可进行相应输入，输入界面如下：图3.7固定晶胞参数输入界面3.3 Weighted Voronoi3D模块Weighted Voronoi3D模块的用户输入界面如下：图3.8Weighted Voronoi3D模块的用户输入界面Cell Info输入参数说明：(1) Ratio列指该相晶体体积和占体积的比率，记为，必须保证; (2) Num列指该相晶体的数量。Ratio列中支持存在一个为“-1”的值，其Ratio值为。3.4 Gradient Voronoi3D模块【功能升级】Gradient Voronoi3D模块的用户输入界面如下：图3.9GradientVoronoi3D模块用户输入界面Cell Size输入参数说明：Size Limit代表Z方向上下边界的晶体大致半径，支持3种分布：Ramp、Pow2和UDF（用户自定义尺寸场），前两种分布晶体大致大小随Z的变化如下： 图(a) Ramp分布 图(b) Pow2分布图3.10GradientVoronoi3D模块晶体大小分布对于UDF分布，用户需自己使用Python编程语言编写一个尺寸场分布程序，即复写getSize(x, y, z)函数中的内容，函数返回值为一个大于0的浮点数，其意义为该点的晶体尺寸(半径)。例如：在1×1×1的方盒中，晶体在方盒中的尺寸分布为：将该尺寸分布函数转为Python程序，参考示例如下：#coding:utf-8import mathdefgetSize(x, y, z): ################################################ size =0.2*math.sqrt((x-0.5)**2+(y-0.5)**2)+0.025 ################################################ returnmax(size,1E-5) if __name__ =="__main__": print(getSize(0.,0.,1.)) 生成结果如下图所示：图3.11圆周分布的梯度晶体模型示例3.5 Periodic Voronoi3D模块【新】Periodic Voronoi3D模块生成用户输入界面如下：图3.12Periodic Voronoi3D模块用户输入界面3.6 UserDefine Voronoi3D模块UserDefineVoronoi3D模块生成用户输入界面如下：图3.13UserDefine Voronoi3D模块自定义Part用户输入界面图3.14UserDefine Voronoi3D模块自定义形核点用户输入界面4. VoronoiEG模块VoronoiEG模块包括：Basic VoronoiEG模块、Mapping VoronoiEG模块、UserDefineVoronoiEG模块、Smoothing VoronoiEG模块和Fluid VoronoiEG模块；4.1 Basic VoronoiEG模块Basic VoronoiEG模块用户输入界面如下：图4.1 Basic VoronoiEG模块用户输入界面4.2 Mapping VoronoiEG模块【新】Mapping VoronoiEG模块用户输入界面如下：图4.2 Mapping VoronoiEG模块用户输入界面4.3 UserDefine VoronoiEG模块【新】UserDefineVoronoiEG模块用户输入界面如下：图4.3UserDefineVoronoiEG模块用户输入界面4.4 Smoothing VoronoiEG模块Smoothing VoronoiEG模块用户输入界面如下：图4.4 Smoothing VoronoiEG模块用户输入界面其中参数Type有两个可选项，分别为Closed和Open。Close代表所有晶体均封闭；Open代表边界处的晶体为开口状态。4.5 Fluid VoronoiEG模块【新】Fluid VoronoiEG模块用户输入界面如下：图4.5Fluid VoronoiEG模块用户输入界面Fluid VoronoiEG模块Period Num可以控制某一方向的复杂程度，数值越大越复杂，取值范围为1-10；Volume Ratio参数用来控制两相的体积比；5. 其他工具5.1 Voronoi Orientation模块目前插件提供晶体随机取向设置模块，其用户界面如下图所示：图5.1 晶体随机取向设置模块用户界面其中Rotation Type参数有4个可选项，分别为Full Random、X Rotation、Y Rotation和Z Rotation。Full Random表示晶体取向完全随机；X Rotation表示所有晶体取向的X方向不变；其他两个选项亦同。5.2 Foam Structure模块【新】目前插件提供由Voronoi 3D模型转泡沫结构模型的一键转换工具，无需手动处理模型，如果选择的模型为0厚度晶界模型（不存在实体晶界部分），则会生成壳模型；否则仅保留实体晶界部分作为泡沫模型，其用户界面如下图所示：图5.2Voronoi 3D模型转泡沫结构模型模块用户界面5.3 Periodic Mesh模块【新】目前插件提供由PeriodicVoronoi模型周期性网格划分工具，可保证在对应的周期面上生成网格的节点一一对应，网格类型仅支持四面体网格，其用户界面如下图所示：图5.3Periodic Mesh模块用户界面6. 插件通用参数说明(1) 生成算法：Uniform和Random算法。Uniform生成的模型每个晶体大小比较均匀，Random算法生成的模型晶体大小比较随机，并受最小距离约束(该参数可一定程度控制晶体均匀度) 图(a) Random算法 图(b) Uniform算法图6.1不同生成算法生成的晶体模型结果(2) Gap参数：Gap参数用于设置晶界的厚度。若Gap值设置为0，代表晶体无实体晶界；若Gap值设置大于0，则会在所有晶体间生成指定厚度的实体晶界，厚度值为Gap值，如图所示：图6.2晶体实体晶界厚度示意图7. 其他可间接生成的模型7.1 线框模型生成的Voronoi模型经过处理后可生成如下图所示的线框模型，可用于梁/杆单元有限元分析：图7.1 线框模型示例7.2 支撑结构模型插件可对自定义模型生成Voronoi支撑结构，示例如下：图7.2 支撑结构模型示例7.3 泡沫结构模型插件生成的Voronoi模型经过处理后可得到泡沫结构模型，可用于泡沫结构的力学性能有限元分析，可生成的泡沫结构模型的类型如下： 图(a) 壳模型 图(b) 实体模型 图(c) 网格模型 图(d) 特殊分布模型(Weighted模型)图7.3泡沫结构模型示例注：任何可生成三维几何Voronoi模型的模块均可用于创建泡沫结构模型。7.4 多种梯度模型插件支持任意梯度分布模型，示例如下： 图(a) 自定义场分布曲线1 图(b) 自定义场分布曲线2 图(a) 自定义场分布曲线1结果 图(b) 自定义场分布曲线2结果图7.4 不同梯度分布模型7.5 几何模型到网格模型的快速映射插件支持Voronoi 3D几何模型到网格模型的映射，示例如下： 图(a) Voronoi 3D几何模型 图(b) 几何到网格的映射结果图7.5 几何模型到网格模型的映射结果7.6 流体网格模型不同参数条件下的流体网格模型，示例如下： 图(a) 流体模型示例1(Period Num：x=4, y=4, z=4, ratio=0.5) 图(b) 流体模型示例2(Period Num：x=4, y=4, z=8, ratio=0.5)图7.6 不同参数条件下的流体网格模型8. 模型生成速度图8.1 模型生成时间随晶体数量变化曲线来源：320科技工作室