24小时热门版块排行榜    

查看: 3490  |  回复: 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的回帖

ddll_slsl

木虫 (正式写手)

有初速度, 下一步骤的速度呢?
看看原子坐标, 有小小移动吗?
受力情况如何?
4楼2013-02-14 09:02:30
已阅   回复此楼   关注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的回帖

paladinlee

铁杆木虫 (正式写手)

★ ★
zh1987hs: 金币+2, 鼓励交流 2013-02-16 10:45:26
引用回帖:
4楼: Originally posted by ddll_slsl at 2013-02-14 09:02:30
有初速度, 下一步骤的速度呢?
看看原子坐标, 有小小移动吗?
受力情况如何?

下一步还是保持原速度,粒子位置没动。
估计是粒子运动积分的group定义有问题。即(fix 4 al  nve)。但是我看说明没有要用到命令积分运动方程啊?
有没有什么方法让喷涂粒子的运动速度保持初速度,然后喷到基底表面的阶段半径内相互作用?
PS:如果对所有粒子all进行NVE积分的话,喷涂粒子有运动,但是运动太剧烈了,并且大部分在insert 的空间内运动,体系温度很高10000K...
我是胖人,不是粗人
5楼2013-02-14 23:41:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 286分人工智能专业请求调剂愿意跨考! +3 lemonzzn 2026-03-17 5/250 2026-03-21 11:28 by lemonzzn
[考研] 336求调剂 +3 rmc8866 2026-03-21 3/150 2026-03-21 11:16 by 求调剂zz
[考研] 286求调剂 +3 Faune 2026-03-21 3/150 2026-03-21 11:14 by ms629
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +7 枫桥ZL 2026-03-18 9/450 2026-03-21 09:56 by JineShine
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 化学求调剂 +4 临泽境llllll 2026-03-17 5/250 2026-03-21 02:23 by JourneyLucky
[考研] 328求调剂,英语六级551,有科研经历 +4 生物工程调剂 2026-03-17 8/400 2026-03-21 02:12 by JourneyLucky
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 281求调剂(0805) +14 烟汐忆海 2026-03-16 25/1250 2026-03-20 15:47 by yuncha
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 085601求调剂 +4 Du.11 2026-03-16 4/200 2026-03-17 17:08 by ruiyingmiao
[考研] 一志愿,福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
信息提示
请填处理意见