首页/文章/ 详情

晶体塑性有限元 Abaqus 三维泰森多边形(voronoi模型)插件 V6.0

7天前浏览280

6.0版本相对于5.0版本做了一系列的功能升级和模块新增,首先升级了基础的三维voronoi模块以及梯度晶模块,然后增加了带权重的二维voronoi和周期性边界,当然三维周期性的voronoi也增加到当前的版本中,此外6.0版本新增了网格映射模块(Mapping VoronoiEG)、用户自定义模块(UserDefine VoronoiEG)、流动两相模块(Fluid VoronoiEG ) 、泡沫结构模块(Foam Structure )以及周期性网格模块(Periodic Mesh).

5.0版本介绍:

晶体塑性有限元 Abaqus 三维泰森多边形(voronoi模型)插件 V5.0

1. 插件启动方式

首先启动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-8

import math


defgetSize(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 Numx=4, y=4, z=4, ratio=0.5)

   

图(b) 流体模型示例2

(Period Numx=4, y=4, z=8, ratio=0.5)

图7.6 不同参数条件下的流体网格模型

8. 模型生成速度

图8.1 模型生成时间随晶体数量变化曲线


来源:320科技工作室
AbaqusUDF通用pythonUM控制
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2025-04-03
最近编辑:7天前
320科技工作室
硕士 | 结构工程师 lammps/ms/vasp/
获赞 225粉丝 384文章 341课程 0
点赞
收藏
作者推荐

基于LAMMPS模拟Cu单晶疲劳循环加载过程

关键词:循环载荷;Cu单晶,塑性变形,位错,lammps循环载荷是指在外力作用下,材料或结构经历周期性应力或应变变化的现象。这种周期性变化通常是由于机械振动、疲劳测试、交变工作环境等因素引起的。循环载荷的大小和方向随时间呈规律性变化,可以是正弦波、方波或其他形式的波形。循环载荷的影响一般采用应力-应变曲线或疲劳寿命实验来确定。通常根据材料在循环载荷下的应力幅值、应变幅值以及循环次数来定义其疲劳性能。应力幅值较小、应变幅值较小且循环次数较少的条件下,材料可能表现出良好的抗疲劳性能;而应力幅值较大、应变幅值较大且循环次数较多的条件下,材料则容易发生疲劳破坏。疲劳寿命是表征材料在循环载荷作用下抵抗破坏能力的参数,宏观上可以通过疲劳实验测量。但是,宏观与微观之间存在差异,例如,微观裂纹的萌生和扩展在宏观实验中可能难以直接观察。为获得循环载荷作用下材料行为的分子模拟和实验结果间的定量比较与普适性解释,通常以材料内部的微观结构变化为特征,极端情况是材料内部无损伤和完全断裂,介于两者之间的所有其他状态都可以认为是不同程度的损伤累积。分子动力学在研究材料的循环加载行为及其微观作用机理方面正逐渐展现出不可替代的价值。以铜(Cu)为例,作为广泛应用的工程材料,其力学性能和循环加载下的响应特性是科研和工业界关注的焦点。选取面心立方(FCC)结构的Cu作为研究对象,其晶格参数来源于标准的晶体学数据库,典型的Cu晶格参数为a=b=c=3.615Å,α=β=γ=90°,形成高度对称的立方晶胞结构。为了模拟实际材料中的加载情况,首先需要构建一个足够大的Cu单晶模型,确保模拟结果能够反映材料的宏观行为而不受模型尺寸的限制。 图1 (a)循环载荷加载曲线;(b)分子动力学模型模型采用第三章中的 (100) 取向立方结构模型,X、Y、Z 三个方向分别对应于 [100]、[010]、[001] 取向,三个方向均采用周期性边界条件以消除边界效应。通过控制应变,采用拉压循环的方式进行加载,应变比为 R =−1 ( R 为每次循环的最小应变与最大应变之比)。加载示意图如图1(a)所示。为了研究循环加载下温度和应变率对疲劳力学性能和变形机理的影响,分别在300K温度下和应变率为1×109s-1的条件下进行了模拟计算,此外,还考虑了正弦形波循环加载对力学性能和变形机制的影响。图2 循坏载荷下的应力应变曲线300 K 时的循环应力-应变曲线如图2所示,当高温合金受到循环加载时,最大应力随循环次数的增加而增大,即首先发生应力循环硬化,这主要是由于初始缺陷的积累,如位错、堆垛层错等。随着加载的进行,循环应力-应变曲线在最后几个循环中基本一致。随着循环次数的增加,最大应力趋于稳定并达到循环饱和状态,这符合金属的循环变形特征。图3 (a)循坏载荷下Cu模型;(b)剪切应变;(c)循环载荷下位错分析;(d)公共邻域分析图4 (a)循坏载荷下Cu模型;(b)剪切应变;(c)循环载荷下位错分析;(d)公共邻域分析图3和图4分别为不同应变下Cu单晶的循环载荷、剪切应变、位错分析、公共邻域分析的可视化图,通过ovito可视化后,可以发现循环载荷下Cu单晶存在明显的应力集中现象,同时发生均匀相变,在Cu单晶内部可以发现存在少量的bcc以及Other原子,这对Cu单晶的变形和力学性能有显著的影响。来源:320科技工作室

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈