24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

Ballendtoo

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 341.3
散金: 3
帖子: 36
在线: 100.3小时
虫号: 2750563
注册: 2013-10-24
专业: 凝聚态物性I:结构、力学和

[求助] 怎么样在Lammps里面固定输出间隔，然后自动输出？已有1人参与

比如400000步的程序，我希望每跑间隔50000步自动输出一个当时的盒子大小x, y, z 该怎么实现呢？

有做过的同学吗
谢谢大家了

回复此楼

» 猜你喜欢

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

lammps中原子固定已经有4人回复
lammps 离子液体，离子中相邻隔的挺远的原子成键了，不知什么原因已经有8人回复
用lammps模拟填充水的碳管的热导率，能设置碳管为固定边界，水为周期性边界条件吗？已经有6人回复
FLUENT瞬态仿真中如何输出固定时间步长下某一个面上所有点的温度已经有19人回复

1楼 2014-04-03 17:39:15

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

Ballendtoo

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 341.3
散金: 3
帖子: 36
在线: 100.3小时
虫号: 2750563
注册: 2013-10-24
专业: 凝聚态物性I:结构、力学和

引用回帖:

4楼: Originally posted by ysshj at 2014-04-11 16:08:30
我只能说这种原因很复杂，导致出现除数很小的情况。
势函数如果按参考文献上的取，感觉一般不会有大错。
要不把in文件，data文件发上来大家看看？...

in 文件应该没有问题，因为跑一个结构可以正确输出结果，但是我把结构从单层改成双层，就输出nan了，具体的如下：

Ewald initialization ...
  G vector (1/distnace) = -nan
  estimated absolute RMS force accuracy = -nan
  estimated relative force accuracy = -nan
  KSpace vectors: actual max1d max3d = 0 1 13
Memory usage per processor = 5.82992 Mbytes
Step Lx Ly Lz Xlo Xhi Ylo Yhi Zlo Zhi
   0       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan
   10       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan
   20       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan
   30       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan
   40       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan
   50       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan
   60       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan
   70       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan
   80       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan
   90       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan
   100       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan       -nan
Loop time of 0.03101 on 16 procs for 100 steps with 11508 atoms

Ewald initialization ...
  G vector (1/distance) = 0.374164
  estimated absolute RMS force accuracy = 1.64594e-05
  estimated relative force accuracy = 1.14304e-06
  KSpace vectors: actual max1d max3d = 38056 28 92596
Setting up run ...
Memory usage per processor = 43.2059 Mbytes
Step Temp E_pair E_mol TotEng Press Volume
   0       1373       -nan          0       -nan       -nan 182551.45
   50       -nan       -nan          0       -nan       -nan       -nan
Loop time of 0.211481 on 16 procs for 50 steps with 11508 atoms

就是这两段出现了问题，但是跑9000个原子相同的条件下就可以输出

赞一下

回复此楼

5楼2014-04-11 19:29:24

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

查看全部 24 个回答

ysshj

金虫 (小有名气)

应助: 3 (幼儿园)
金币: 669.6
红花: 4
帖子: 111
在线: 43.3小时
虫号: 2237709
注册: 2013-01-12
性别: GG
专业: 机械摩擦学与表面技术

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
Ballendtoo: 金币+2, ★★★很有帮助, 十分感谢~~~~~我先试试看，确实是需要输出lx ly lz 2014-04-03 19:05:55
月只蓝: 金币+2, 感谢指导！ 2014-04-07 09:54:16
Ballendtoo: 金币+3, ★★★★★最佳答案 2014-04-11 19:32:53

thermo 50000
thermo_style custom step lx ly lz xlo xhi ylo yhi zlo zhi

赞一下

回复此楼

一个人要像一支队伍。