| 查看: 3116 | 回复: 6 | ||
[求助]
【求助】lammps模拟碰撞竟然丢失原子，难道是势函数的问题？？ 已有2人参与
|
 建了一个模型碰撞的时候一直丢失两个原子，都要疯了，为了检查错误，直接建立了一个简单盒子模拟碰撞，竟然还是报错丢失两个原子！请各路大神帮忙瞄瞄，很简单的小程序，咋会一直报错呢？请不要用什么截断半径过小，还是飞出盒子这样的话回复我。因为我觉得完全不可能飞出盒子，步长也调到过0.00001ps，但是不解决问题，而且是周期边界，丢失原子觉得很奇怪啊！ velocity 0 0 500的时候没问题，但是速度一提高就报错，比如velocity 0 0 1000。实在是不知道咋办了，所以请各路大神帮忙看看，小弟先谢过了。可能有些同学打开附件有问题，所以我直接把代码粘贴出来了，势能是lammps自带的势能 units metal boundary p p p atom_style atomic lattice bcc 2.855 #orient x 1 1 1 orient y -1 -1 2 orient z 1 -1 0 region box block -16 16 -15 15 -20 20 create_box 1 box mass 1 56.0 create_atoms 1 box pair_style eam/alloy pair_coeff * * Fe_mm.eam.fs Fe timestep 0.001 neighbor 1.5 bin neigh_modify every 30 delay 0 check no #displace_box all z delta -2 2 remap none compute pe all pe/atom compute dis all displace/atom dump 1 all custom 500 2.lammpstrj id type x y z vx vy vz c_pe c_dis[1] fix 3 all nvt temp 300 300 0.1 run 1000 unfix 3 write_restart dislocation.restart ——————————————————下面是restart文件 read_restart dislocation.restart lattice bcc 2.855 #orient x 1 1 1 orient y -1 -1 2 orient z 1 -1 0 pair_style eam/alloy pair_coeff * * Fe_mm.eam.fs Fe timestep 0.0001 neighbor 1.5 bin neigh_modify every 10 delay 0 check no compute pe all pe/atom compute dis all displace/atom dump 2 all custom 100 4.lammpstrj id type x y z vx vy vz c_pe c_dis[1] group kpa id 62410 velocity kpa set 0.0 0.0 1000.0 units box fix 5 all nve run 1000 unfix 5 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : 超简单碰撞.zip
2016-04-01 10:16:49, 328.17 K
» 猜你喜欢
各位老师好，求调剂，本科211，一志愿天津大学生物与医药学硕，差两名录取。
已经有3人回复
-
化工学硕294分，求导师收留
已经有20人回复
-
327求调剂
已经有8人回复
-
297求调剂
已经有22人回复
-
085801电气专硕272求调剂
已经有11人回复
-
271求调剂
已经有33人回复
-
调剂
已经有27人回复
-
一志愿2110，化学学硕310分，本科重点双非求调剂
已经有20人回复
-
085600材料与化工，求调剂
已经有13人回复
-
本科郑州大学，一志愿华东师范大学282求调剂
已经有31人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】lammps关于fix ave/spatial命令应用问题
已经有4人回复
- 【求助】lammps中是否有关于Bond-orientational order参数的
已经有3人回复
- 【求助】Lammps使用？
已经有0人回复
- 【求助】LAMMPS-运行报错
已经有9人回复
- 【求助】lammps势的选择
已经有3人回复
- 【求助】lammps模拟聚合物刷如何建模？
已经有1人回复
- 【求助】LAMMPS构造晶界
已经有0人回复
- 【求助】lammps或者dlpoly能否做粗粒化模拟吗？
已经有3人回复
- 【求助】Lammps能做布朗动力学模拟吗？
已经有3人回复
- 【求助】lammps怎么算声子谱？
已经有1人回复
- 【求助】lammps的使用
已经有9人回复
- 【求助】lammps 安装问题
已经有1人回复
- 【求助】Lammps 中加外力的问题
已经有1人回复
- 【求助】LAMMPS中添加Lees-Edwards边界条件
已经有2人回复
- 【求助】lammps这段代码的意思是？
已经有1人回复
- 【求助】lammps计算heat/flux出错
已经有1人回复
- 【求助】lammps如何测量模拟以后材料的晶格参数？
已经有7人回复
- 【求助】LAMMPS－REAX计算出错
已经有3人回复
- 【求助】lammps中多元化合物的晶格创建问题
已经有14人回复
- 【求助】lammps中消除平动对温度统计的影响问题
已经有2人回复
- 【求助】Lammps 多机并行的问题
已经有11人回复
hooge
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 1790.1
- 散金: 236
- 红花: 1
- 帖子: 975
- 在线: 135.1小时
- 虫号: 89636
- 注册: 2005-08-29
- 专业: 金属材料的力学行为
2楼2016-04-01 13:09:58
3楼2016-04-01 13:15:09
whyaquarius
新虫 (初入文坛)
- 应助: 4 (幼儿园)
- 金币: 253.8
- 散金: 70
- 帖子: 30
- 在线: 34.7小时
- 虫号: 2175398
- 注册: 2012-12-09
- 性别: GG
- 专业: 核技术及其应用
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
大懒猫王浩: 金币+300, ★★★★★最佳答案, 师兄竟然来赚金币了，，，，不好意思，只能给师兄了。。。 2016-04-01 17:29:35
月只蓝: IP检查，初步认定为马甲行为，涉嫌转移金币，违反论坛相关管理规定：http://emuch.net/bbs/viewthread.php?tid=3711903 2016-04-02 10:16:56
月只蓝: 经查，您和大懒猫王浩的IP地址十分相近，请于2天内通过站内短消息解释您和大懒猫王浩的关系！逾期按马甲处理，并可能追加额外金币处罚。 2016-04-02 10:18:04
月只蓝: 金币-50, 认定为金币转移行为，核减金币300×0.8=240个，现扣除50个，还需要扣除190个。 2016-04-07 13:38:49
月只蓝: 金币-50, 认定为金币转移行为，核减金币300×0.8=240个，现扣除50个，还需要扣除140个。 2016-04-07 13:39:08
月只蓝: 金币-50, 认定为金币转移行为，核减金币300×0.8=240个，现扣除50个，还需要扣除90个。 2016-04-07 13:39:34
月只蓝: 金币-50, 认定为金币转移行为，核减金币300×0.8=240个，现扣除50个，还需要扣除40个。 2016-04-07 13:40:20
月只蓝: 金币-40, 认定为金币转移行为，核减金币300×0.8=240个，现已扣除完毕。 2016-04-07 13:40:42
月只蓝: 考虑到用户可能不了解马甲行为规定且是初犯，免除额外处罚。 2016-04-07 13:42:00
感谢参与，应助指数 +1
大懒猫王浩: 金币+300, ★★★★★最佳答案, 师兄竟然来赚金币了，，，，不好意思，只能给师兄了。。。 2016-04-01 17:29:35
月只蓝: IP检查，初步认定为马甲行为，涉嫌转移金币，违反论坛相关管理规定：http://emuch.net/bbs/viewthread.php?tid=3711903 2016-04-02 10:16:56
月只蓝: 经查，您和大懒猫王浩的IP地址十分相近，请于2天内通过站内短消息解释您和大懒猫王浩的关系！逾期按马甲处理，并可能追加额外金币处罚。 2016-04-02 10:18:04
月只蓝: 金币-50, 认定为金币转移行为，核减金币300×0.8=240个，现扣除50个，还需要扣除190个。 2016-04-07 13:38:49
月只蓝: 金币-50, 认定为金币转移行为，核减金币300×0.8=240个，现扣除50个，还需要扣除140个。 2016-04-07 13:39:08
月只蓝: 金币-50, 认定为金币转移行为，核减金币300×0.8=240个，现扣除50个，还需要扣除90个。 2016-04-07 13:39:34
月只蓝: 金币-50, 认定为金币转移行为，核减金币300×0.8=240个，现扣除50个，还需要扣除40个。 2016-04-07 13:40:20
月只蓝: 金币-40, 认定为金币转移行为，核减金币300×0.8=240个，现已扣除完毕。 2016-04-07 13:40:42
月只蓝: 考虑到用户可能不了解马甲行为规定且是初犯，免除额外处罚。 2016-04-07 13:42:00
|
来来来，我来赚金币啦~ 在辐照损伤的模拟中，在离子间距离比较小时，已经不能用平衡态下的作用势来描述，一般用 ZBL 来描述这个碰撞过程。比如说常用到的 tersoff/zbl 作用势就势就是将 Tersoff 作用势与 ZBL 作用势通过一个费米函数连接起来，一般 认为小于 1Å 时需要用 ZBL 描述。 P.S. eam势木有官方的zbl衔接版本~ |
4楼2016-04-01 17:20:32
5楼2016-04-01 18:13:55
xxxxzzz
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1682
- 散金: 200
- 帖子: 303
- 在线: 83小时
- 虫号: 2018899
- 注册: 2012-09-21
- 性别: GG
- 专业: 核技术及其应用
6楼2016-04-03 08:16:51
7楼2016-12-02 15:44:09
