24小时热门版块排行榜    

查看: 1922  |  回复: 1

11B329

金虫 (正式写手)

[求助] 关于lammps里examples的问题(shear) 已有1人参与

各位大神,希望能指点一下。
在equilibrate之前,lower与upper的温度并没有设定,只是设定了mobile的温度,请问那其余两组原子的初速度怎样?还有就是原子的质量好像没有给啊,不用给定吗?
以下是shear的程序:
# 3d metal shear simulation


units                metal

boundary        s s p


atom_style        atomic

lattice                fcc 3.52

region                box block 0 16.0 0 10.0 0 2.828427

create_box        3 box


lattice                fcc 3.52 orient        x 1 0 0 orient y 0 1 1 orient z 0 -1 1
&
                origin 0.5 0 0

create_atoms        1 box


pair_style        eam

pair_coeff        * * Ni_u3.eam


neighbor        0.3 bin

neigh_modify        delay 5


region                lower block INF INF INF 0.9 INF INF

region                upper block INF INF 6.1 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 100 dump.shear


#dump                2 all image 100 image.*.jpg type type &

#                axes yes 0.8 0.02 view 0 0 zoom 1.5 up 0 1 0 adiam 2.0

#dump_modify        2 pad 4


#dump                3 all movie 100 movie.mpg type type &

#                axes yes 0.8 0.02 view 0 0 zoom 1.5 up 0 1 0 adiam 2.0

#dump_modify        3 pad 4


thermo                100

thermo_modify        temp new2d


reset_timestep        0

run                3000
回复此楼
往前走吧!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mgqqlwq

铜虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
11B329: 金币+10, ★★★很有帮助, 多谢!我一会试试! 2015-12-11 07:42:15
我记得eam势函数文件里有原子质量的信息,你可以看下pair_style eam的doc page。未设定初始温度的group的温度我不知道,不过是否可以试着把那两个group的原子速度输出出来看一下。希望对你有用。
2楼2015-12-10 18:59:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 11B329 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然评审意见 +12 wangmingqi 2026-08-28 17/850 2026-08-29 08:43 by appleren
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 8/400 2026-08-29 08:32 by foolishmani
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 07:59 by gy1nBQXYQJqL
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 07:55 by popo0625
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 07:46 by gy1nBQXYQJqL
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 05:45 by gy1nBQXYQJqL
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +8 苏东坡二世 2026-08-23 8/400 2026-08-28 14:41 by lkforward
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +4 gdd2018 2026-08-28 9/450 2026-08-28 14:22 by jklily
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见