24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3551  |  回复: 12
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

paladinlee

铁杆木虫 (正式写手)

[求助] lammps deposit 粒子初始速度为0?

大家新年好!
我最近在用lammps算金属表面喷涂,使用deposit命令,但是我按照说明设定了添加粒子的初速度,但是粒子速度一直为0,不运动。想请教下各位是什么原因?
我的in文件如下:
...
#Al为喷涂粒子type 2,Cu为基底 type 1,EAM势能,单位metal
...
region 4 block 0.0 72.3 0.0 72.3 210.0 225.0 side in units box  #设定添加粒子region,离基底Z方向距离20nm
...
group substrate region 1  #基底group
group al type 2    #喷涂的Al
...
compute 1 cuthermal temp
compute 2 cu temp
compute 3 al temp
compute_modify 3 dynamic yes
fix 1 substrate setforce        0        0        0  #基底fix
fix 2 cuthermal  nvt temp        300.0        300.0        10  #NVT恒温Cu
fix 3 cu  nve      #NVE恒定接触面
#fix 4 al  nve      #NVE恒定喷涂粒子,但是我尝试了有没有这条命令Al都没有运动
fix coat all deposit 10000 2 64 12345 region 4 near 5.0 vz -0.03 -0.03 units box    #deposit命令,是这有问题?

timestep   0.005  #fs
thermo     1000
thermo_style custom step temp ke pe etotal c_1 c_2 c_3
dump 1 all xyz 1000 metal.xyz
run 1000000
附图为计算结果,褐色粒子为Al,一直没有move,后续粒子一直添加到定义的region中。

谢谢大家!

1.jpg
回复此楼

» 猜你喜欢

我是胖人,不是粗人
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paladinlee

铁杆木虫 (正式写手)

送鲜花一朵
引用回帖:
6楼: Originally posted by ddll_slsl at 2013-02-15 08:42:55
具体的只有靠你自己去摸索, 看dump等输出,仔细琢磨。 你的模型不算复杂,应该容易实现的。

请问喷涂粒子的group的该怎么定义?我用type或者region之类的定义,产生初始喷涂粒子group的粒子数为0,最终导致粒子不能运动。
我是胖人,不是粗人
8楼2013-02-16 06:02:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

ddll_slsl

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
paladinlee: 金币+70, ★★★很有帮助, 10 2013-02-14 01:11:43
zh1987hs: 金币+2, 谢谢 2013-02-16 10:45:01
dump出原子速度看看?
2楼2013-02-13 17:47:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paladinlee

铁杆木虫 (正式写手)

引用回帖:
2楼: Originally posted by ddll_slsl at 2013-02-13 17:47:07
dump出原子速度看看?

dump 喷涂粒子有初速度,但是不运动?怎么回事?
我是胖人,不是粗人
3楼2013-02-14 01:11:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddll_slsl

木虫 (正式写手)

有初速度, 下一步骤的速度呢?
看看原子坐标, 有小小移动吗?
受力情况如何?
4楼2013-02-14 09:02:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 344 材料专业 求调剂211 无地域要求 +6 hualkop 2026-04-11 6/300 2026-04-12 15:00 by seattle40
[考研] 电气工程专硕320求调剂 +5 小麻子111 2026-04-10 5/250 2026-04-12 10:47 by zhouyuwinner
[考研] 297求调剂 +17 ORCHID1 2026-04-10 18/900 2026-04-12 09:59 by zhouwenxian
[考研] 268分085602化学工程调剂 +28 月照花林。 2026-04-09 28/1400 2026-04-12 02:38 by 秋豆菜芽
[考研] 332求调剂 +14 蕉蕉123 2026-04-10 14/700 2026-04-12 00:27 by 蓝云思雨
[考研] 22专硕求调剂 +6 haoyun上岸 2026-04-11 8/400 2026-04-11 23:21 by labixiaoqiao
[考研] 一志愿085502,267分求调剂 +11 再忙也要吃饭啊 2026-04-08 12/600 2026-04-11 11:08 by zhq0425
[考研] 266求调剂 +29 阳阳哇塞 2026-04-07 29/1450 2026-04-10 16:20 by 高维春
[考研] 293调剂 +25 yj1221 2026-04-08 26/1300 2026-04-10 15:02 by 柴小白
[考博] 博士自荐 +7 可可小胖 2026-04-08 7/350 2026-04-10 08:28 by kimhero
[考研] 本科211 工科085400 280分求调剂 可跨专业 +3 LZH(等待调剂中 2026-04-09 3/150 2026-04-09 21:29 by wutongshun
[考研] 材料299专硕求调剂 +10 +21 2026-04-09 10/500 2026-04-09 17:34 by 1753564080
[考研] 265求调剂 +4 风说她早忘了 2026-04-07 4/200 2026-04-09 13:59 by only周
[考研] 296求调剂 +3 汪!?! 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
[考研] 专硕0854初试考材科基,求调剂 +7 3220548044 2026-04-06 10/500 2026-04-08 21:59 by hypershenger
[考研] 材料调剂 +11 一样YWY 2026-04-07 11/550 2026-04-07 15:13 by shdgaomin
[考研] 调剂 +4 mcbbc 2026-04-06 5/250 2026-04-07 12:33 by upczlm1989
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 09:25 by 小黑不怕难
[考研] 22408 331分求调剂 +4 y__1 2026-04-06 4/200 2026-04-06 17:26 by 土木硕士招生
[考研] 求调剂到材料 +5 程9915 2026-04-06 5/250 2026-04-06 15:21 by yulian1987
信息提示
请填处理意见