24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1839  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 复试调剂,一志愿南农083200食品科学与工程 +4 XQTJZ 2026-03-26 4/200 2026-03-27 01:23 by dgnhs
[考研] 一志愿太原理工安全工程300分,求调剂 +3 0857求调剂. 2026-03-24 3/150 2026-03-27 01:18 by 求知若渴lz
[考研] 071000生物学求调剂,初试成绩343 +6 小小甜面团 2026-03-25 6/300 2026-03-26 23:01 by 不吃魚的貓
[考研] 一志愿211院校 344分 东北农业大学生物学学硕,求调剂 +3 丶风雪夜归人丶 2026-03-26 3/150 2026-03-26 19:35 by 求调剂zz
[考研] 0703化学/290求调剂/本科经历丰富/工科也可 +7 丹青奶盖 2026-03-26 7/350 2026-03-26 19:18 by macy2011
[考研] 290分调剂求助 +3 吉祥止止陈 2026-03-25 3/150 2026-03-25 19:58 by barlinike
[考研] 材料学硕,求调剂 6+4 糖葫芦888ll 2026-03-22 9/450 2026-03-25 11:19 by greychen00
[考研] 085600材料与化工调剂 +9 A-哆啦Z梦 2026-03-23 15/750 2026-03-25 11:18 by Ainin_
[考研] 289材料与化工(085600)B区求调剂 +4 这么名字咋样 2026-03-22 5/250 2026-03-25 08:20 by mx.yue
[考研] 一志愿北化315 求调剂 +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by 了了了了。。
[考研] 一志愿南航材料专317分求调剂 +5 炸呀炸呀炸薯条 2026-03-23 5/250 2026-03-24 16:52 by 星空星月
[考研] 材料专硕331求调剂 +4 鲜当牛 2026-03-24 4/200 2026-03-24 15:58 by JourneyLucky
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +7 vv迷 2026-03-22 7/350 2026-03-23 23:44 by Txy@872106
[考研] 384求调剂 +3 子系博 2026-03-22 6/300 2026-03-23 21:45 by 子系博
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
[考研] 293求调剂 +3 涛涛Wjt 2026-03-22 5/250 2026-03-22 22:21 by jiangpengfei
[考研] 319求调剂 +4 小力气珂珂 2026-03-20 4/200 2026-03-22 15:53 by ColorlessPI
[考研] 269专硕求调剂 +6 金恩贝 2026-03-21 6/300 2026-03-22 14:31 by ColorlessPI
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 一志愿南理工085701环境302求调剂院校 +3 葵梓卫队 2026-03-20 3/150 2026-03-20 19:28 by zhukairuo
信息提示
请填处理意见