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

分子模拟论坛 Molecular Simulation Forums

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 1385|回复: 2

关于displace_atoms

[复制链接]

50

主题

150

帖子

356

积分

中级会员

Rank: 3Rank: 3

积分
356
发表于 2010-10-26 11:03:00 | 显示全部楼层 |阅读模式
这是我结合几个帖子写的
#file name: in.c11
units metal
boundary s s s
atom_style atomic
neighbor2.0 bin
neigh_modifydelay 1 check yes
lattice fcc 3.615
region box block 0 30 0 3 0 3
create_box 1 box
create_atoms 1 box
mass 1 63.546
#define group
region 1 block INF 1 INF INF INF INF
group left region 1
region 2 block 29 INF INF INF INF INF
group right region 2
group mobile subtract allleft
pair_style eam
pair_coeff * * Cu_u3.eam
compute p all pressure thermo_temp
variablepressx equal c_p[1]
variablepressy equal c_p[2]
variablepressz equal c_p[3]
thermo 1000
thermo_style custom step temp etotal press v_pressxv_pressy v_presszvol
velocityleft set 0.0 0.0 0.0
minimize 1.0e-6 1.0e-6 1000 10000
min_style cg
dump1 all atom 1000 dump.cupoly
fix 1 left setforce 0.0 NULL NULL
fix2mobile temp/rescale 100 1.0 1.0 5 1.0
fix 3 all nve
timestep 0.001
run 40000
computes all stress/atom
computestrea all reduce sum c_s[1]
variablestressa equal c_strea/vol
computestreb all reduce sum c_s[2]
variablestressb equal c_streb/vol
computestrec all reduce sum c_s[3]
variablestressc equal c_strec/vol
thermo_stylecustom step temp etotal press v_stressa v_pressx v_stressb v_stressc vol
label begin
variable t loop 100
displace_atoms mobileramp x 1 30 x 0 0.01
minimize 1.0e-6 1.0e-6 1000 10000
run 10000
next t
jump in.c11 begin
结果是从固定的left 原子处断裂,更没有necking现象。我想是displace_atoms 有问题吧。哪位告诉下正确的displace_atoms 用法,O(∩_∩)O谢谢。之所以选择displace_atoms,是因为deform 拉伸时stress-strain 不是很线性,不知道displace_atoms 线性符合的怎么样?谢谢指教
回复

使用道具 举报

6

主题

19

帖子

50

积分

注册会员

Rank: 2

积分
50
发表于 2011-12-20 13:35:00 | 显示全部楼层
你的mobile怎么定义的是all-left啊这不相当于把全部都拉么肯定是从left断裂了
回复

使用道具 举报

0

主题

1

帖子

24

积分

新手上路

Rank: 1

积分
24
发表于 2017-2-28 16:52:17 | 显示全部楼层
想问一下,这个问题应该怎么解决啊,求助各位大神
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2020-6-3 08:15 , Processed in 0.062229 second(s), 26 queries .

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

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