24小时热门版块排行榜    

查看: 2335  |  回复: 3

naterriver

至尊木虫 (正式写手)

[求助] lammps中怎么固定一族原子 已有1人参与

现在做一个结构时core-shell模型,想要固定住core单独加热shell,小木虫的帖子搜出来的帖子和手册里的命令也试过,但一直没有成功过。

因为是刚刚学习lammps,请有经验的虫友指教下!

units metal
boundary p p p
atom_style atomic

read_restart ./restart/PdPt_relax.112.restart

pair_style eam/alloy
pair_coeff * * PdPt.set Pd Pt
neighbor 2.0 bin
neigh_modify delay 1 every 1

group core type 1
group shell type 2

compute pe all pe/atom
compute myRDF all rdf 100
fix 1 shell nvt temp 2000 2000 0.2
fix 2 core setforce 0.0 0.0 0.0

velocity core set 0 0 0 sum no units box

dump 1 all custom 2000 ./config/PdPt_2000.*.dump id type c_pe x y z

restart 10000 ./restart/PdPt_2000.*.restart
thermo 1000
thermo_style custom step temp lx ly lz press pxx pyy pzz
timestep 0.001

run 2000   every 1 "velocity core set 0 0 0 sum no units box"

图片1是原来的结构,图片2就是计算之后的结构,不知道为什么在x方向被压成了一层结构

lammps中怎么固定一族原子
1.jpg


lammps中怎么固定一族原子-1
2.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

chinaboy798

银虫 (小有名气)

【答案】应助回帖

velocity  sum    yes  呢?不要sum  no。
2楼2014-12-15 19:09:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luzihen

金虫 (正式写手)

为什么不用fix velocity 把core的速度设成0就行了,不用那么麻烦吧
3楼2014-12-18 17:58:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

没有烟的火柴

新虫 (初入文坛)

velocity low set 0.0 0.0 0.0 units box                                          
velocity up set 0.0 0.0 0.0 units box
4楼2014-12-22 16:38:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 naterriver 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-30 22:24 by xiangy672
[基金申请] 国社科又开始会评了,不知道这次命运如何 +5 雨打竹帘 2026-08-30 9/450 2026-08-30 22:23 by 余韵清
[基金申请] 能否申诉? +4 echo8914667 2026-08-30 5/250 2026-08-30 21:11 by echo8914667
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 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 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
信息提示
请填处理意见