24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1838  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 343求调剂 +4 赠我一本书 2026-03-23 4/200 2026-03-27 00:40 by wxiongid
[考研] 349求调剂 +5 杰斯塔里斯 2026-03-21 5/250 2026-03-27 00:31 by wxiongid
[考研] 材料学硕333求调剂 +7 北道巷 2026-03-24 7/350 2026-03-27 00:30 by wxiongid
[考研] 336材料求调剂 +7 陈滢莹 2026-03-26 9/450 2026-03-27 00:20 by wxiongid
[考研] 071000生物学求调剂,初试成绩343 +6 小小甜面团 2026-03-25 6/300 2026-03-26 23:01 by 不吃魚的貓
[考研] 327求调剂 +7 prayer13 2026-03-23 7/350 2026-03-26 20:48 by 不吃魚的貓
[考研] 材料调剂 +7 匹克i 2026-03-23 7/350 2026-03-26 19:00 by 不吃魚的貓
[考研] 材料277求调剂 +5 min3 2026-03-24 5/250 2026-03-26 15:13 by zzll406
[考研] 0856求调剂 +8 zhn03 2026-03-25 9/450 2026-03-26 13:42 by zzll406
[考研] 总分293求调剂 +6 加一一九 2026-03-25 8/400 2026-03-26 13:30 by yujianx
[考研] 303求调剂 +6 蓝山月 2026-03-25 6/300 2026-03-25 22:47 by 418490947
[考研] 材料与化工328分调剂 +6 。,。,。,。i 2026-03-23 6/300 2026-03-25 22:30 by 418490947
[考研] 334分 一志愿武理-080500 材料求调剂 +4 李李不服输 2026-03-25 4/200 2026-03-25 21:26 by 星空星月
[考研] 材料专硕 335 分求调剂 +4 拒绝冷暴力 2026-03-25 4/200 2026-03-25 18:45 by haxia
[考研] 材料学硕,求调剂 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_
[考研] 0703化学求调剂 +6 奶油草莓. 2026-03-22 7/350 2026-03-25 10:00 by shangxh
[考研] 340求调剂 +5 话梅糖111 2026-03-24 5/250 2026-03-25 06:53 by ilovexiaobin
[考研] 一志愿北化315 求调剂 +3 akrrain 2026-03-24 3/150 2026-03-24 19:35 by 了了了了。。
[考研] 361求调剂 +3 Glack 2026-03-22 3/150 2026-03-23 22:03 by fuyu_
信息提示
请填处理意见