24小时热门版块排行榜    

查看: 989  |  回复: 2

Han19910102

银虫 (初入文坛)

[求助] 单晶ni加压 已有1人参与

如何在这个in文件中给Ni加应力
# 3d metal shear simulation

units                metal
boundary        s s p

atom_style        atomic
lattice                fcc 3.52  orient                x 1 0 0  orient                y 0 1 1 orient                z 0 -1 1 origin                0.5 0 0
region                box block 0 16.0 0 10.0 0 2.828427
create_box        3 box


create_atoms        1 box

pair_style        eam
pair_coeff        * * niu3

neighbor        0.3 bin
neigh_modify        delay 5

region                lower block INF INF INF 1.4 INF INF
region                upper block INF INF 8.6 INF INF INF
group                lower region lower
group                upper region upper
group                boundary union lower upper
group                mobile subtract all boundary

set                group lower type 2
set                group upper type 3

# void

#region                void cylinder z 8 5 2.5 INF INF
#delete_atoms        region void

# temp controllers

compute         new3d mobile temp
compute         new2d mobile temp/partial 0 1 1

# equilibrate

velocity        mobile create 300.0 5812775 temp new3d
fix                1 all nve
fix                2 boundary setforce 0.0 0.0 0.0

fix                3 mobile temp/rescale 10 300.0 300.0 10.0 1.0
fix_modify        3 temp new3d

thermo                25
thermo_modify        temp new3d

timestep        0.001
run                100

# shear

velocity        upper set 1.0 0 0
velocity        mobile ramp vx 0.0 1.0 y 1.4 8.6 sum yes

unfix                3
fix                3 mobile temp/rescale 10 300.0 300.0 10.0 1.0
fix_modify        3 temp new2d

dump                1 all atom 10 dump.shear
thermo                100
thermo_modify        temp new2d

reset_timestep        0
run                3000
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : dump.shear
  • 2015-09-02 20:46:31, 38.86 M
  • 附件 2 : niu3
  • 2015-09-02 20:46:40, 35.85 K

» 收录本帖的淘帖专辑推荐

分子动力学模拟

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alundilong

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
Han19910102: 金币+2, ★★★很有帮助 2015-09-04 18:30:40
try to add strain.
fix                2 R1 deform 1 x erate 1e-2 units box remap x
2楼2015-09-03 02:10:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qiaoshuai55

银虫 (初入文坛)

多看书,加油!
更努力,更幸运!
3楼2015-09-03 21:23:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Han19910102 的主题更新
信息提示
请填处理意见