请选择 进入手机版 | 继续访问电脑版

分子模拟论坛 Molecular Simulation Forums

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 1093|回复: 0

材料拉伸过程的剪切带分裂

[复制链接]

2

主题

2

帖子

83

积分

注册会员

Rank: 2

积分
83
发表于 2018-12-10 08:02:11 | 显示全部楼层 |阅读模式
请问有人在模拟非晶材料拉伸过程中发生剪切带分裂的问题吗?请教一下问题出在什么地方?谢谢。

我是在Ovito中观察的von mises strain。在Ovito是选择atomic strain 和 color coding,选择shear strain,并关闭 Use minimum image convention。


附上代码:
# ----- Initialization
units                        metal
dimension                3
boundary                s p p
atom_style                atomic

# ----- Variable
variable                t equal 10

# ----- Atom definition
read_data                data_LiSi.txt
mass                        1 6.941
mass                        2 28.09

# ----- Force field
pair_style                meam
pair_coeff                * * library.meam Li Si LiSi.meam Li Si
neighbor           3.0 bin
neigh_modify        every 5 delay 10 check yes

# ----- Minimization
reset_timestep        0
timestep                0.001
velocity                all create $t 589665 dist gaussian units box

thermo                        10000
thermo_style             custom step temp etotal ke pe press pxx pyy pzz pxy pxz pyz vol lx ly lz

# ----- Deformation
reset_timestep             0

variable                tmp equal "lz"
variable                L0 equal ${tmp}
print                        "Initial length is ${L0}"
variable                strain equal "(lz-v_L0)/v_L0"
variable                p1 equal "-pxx/10000"
variable                p2 equal "-pyy/10000"
variable                p3 equal "-pzz/10000"
variable                N equal 'step'
variable                V equal "vol"
variable                T equal "temp"
variable                PE equal "pe"
variable                ETOTAL equal "etotal"
variable                p equal "press"
variable                Lx equal "lx"
variable                Ly equal "ly"
variable                Lz equal "lz"

fix                        def1 all print 5000 "${N} ${T} ${strain} ${p1} ${p2} ${p3} ${p} ${V} ${Lx} ${Ly} ${Lz} ${PE} ${ETOTAL}" file tensile_data.txt screen no

dump                        1 all custom 10000 dump_LiSi.*.txt mass type xs ys zs fx fy fz

fix                        2 all npt temp $t $t 0.1 y 0.0 0.0 1.0 drag 1.0

variable                srate equal 1.0e10
variable                srate1 equal ${srate}/1.0e12

fix                        3 all deform 100 z erate ${srate1} units box remap x

run                        700000

write_data                data_LiSi_tensile.txt

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Archiver|手机版|小黑屋|分子模拟论坛

GMT+8, 2020-3-29 02:17 , Processed in 0.058643 second(s), 23 queries .

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表