| 查看: 1725 | 回复: 11 | |||
lchunmei001木虫 (正式写手)
|
[求助]
冲击结果异常
|
|
请教大家，我沿x方向给粒子2KM/s的速度冲击加载，具体命令如下，但得到的结果很奇怪，压强很低的情况下，温度就已达到几万K，显然很不合理，但是在不知道是哪里出问题了，测试了几种物质，均出现这种情况，应该是命令哪里没对，还请高手指点一二，感激不尽！ velocity all set 20 NULL NULL sum yes units box fix 1 all nve fix 2 all wall/reflect xlo xhi |
回复此楼» 收录本帖的淘帖专辑推荐
 condensed matter physics |
» 猜你喜欢
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有13人回复
-
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性
已经有11人回复
-
表哥与省会女结婚，父母去帮带孩子被省会女气回家生重病了
已经有9人回复
-
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有10人回复
-
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有11人回复
ddll_slsl
木虫 (正式写手)
- 应助: 101 (高中生)
- 金币: 2939.2
- 红花: 6
- 帖子: 670
- 在线: 255.3小时
- 虫号: 2121784
- 注册: 2012-11-12
- 专业: 动力学与控制
2楼2013-01-16 11:11:25
lchunmei001
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1738.5
- 散金: 55
- 红花: 1
- 帖子: 797
- 在线: 151.3小时
- 虫号: 863507
- 注册: 2009-10-05
- 专业: 凝聚态物性I:结构、力学和
|
输入和结果如下，请予以指点，谢谢！ 输入： units metal boundary s p p lattice fcc 4.07 atom_style atomic region box block 0 100 0 10 0 10 create_box 1 box create_atoms 1 box pair_style eam pair_coeff * * Pt_u3.eam timestep 1.0e-3 neighbor 1.0 bin neigh_modify every 1 delay 10 check yes thermo 100 thermo_style custom step press vol temp pxx pyy pzz lx ly lz velocity all set 20 NULL NULL sum yes units box fix 1 all nve fix 2 all wall/reflect xlo xhi dump id all atom 10000 dump.melt run 50000 结果log文件： Step Press Volume Temp Pxx Pyy Pzz Lx Ly Lz 0 19146.265 674326.27 31285.986 535234.08 -238897.64 -238897.64 407.0814 40.7 40.7 100 23424.423 673363.58 31311.738 540451.5 -235089.11 -235089.11 406.50024 40.7 40.7 200 27653.106 671960.25 31331.714 546177.73 -231609.21 -231609.21 405.65307 40.7 40.7 300 31059.835 671084.06 31347.923 550838.79 -228829.64 -228829.64 405.12413 40.7 40.7 400 36428.574 669579.59 31362.412 557393.91 -224054.09 -224054.09 404.21589 40.7 40.7 500 39253.628 668750.35 31377.061 561471.39 -221855.25 -221855.25 403.7153 40.7 40.7 600 44925.431 667277.7 31394.038 568443.51 -216833.61 -216833.61 402.82628 40.7 40.7 700 47685.276 666310.63 31413.841 572717.7 -214830.93 -214830.93 402.24247 40.7 40.7 。。。。。。 |
3楼2013-01-16 11:36:01
ddll_slsl
木虫 (正式写手)
- 应助: 101 (高中生)
- 金币: 2939.2
- 红花: 6
- 帖子: 670
- 在线: 255.3小时
- 虫号: 2121784
- 注册: 2012-11-12
- 专业: 动力学与控制
4楼2013-01-16 14:55:28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
lchunmei001: 金币+5 2013-01-16 16:23:13
zh1987hs: 金币+3, 谢谢 2013-01-19 14:04:12
感谢参与，应助指数 +1
lchunmei001: 金币+5 2013-01-16 16:23:13
zh1987hs: 金币+3, 谢谢 2013-01-19 14:04:12
|
楼主，我有点儿不懂，请教你哈： 1、你是要模拟一个什么样的冲击呢？我把你的in文件运行了一下，4万多个原子都一起向x方向运动吗？安装我的理解，所谓的“冲击”应该是一个物体以较高的速度撞击另一个物体，比如炮弹攻击掩体的时候。我看过使用lammps模拟纳米颗粒冲击石墨烯的。等等这些。我没弄明白你的模拟要干啥。 2、在你的in文件中，你能告诉我哪句命令是thermostat的？你设置的温度是多少？ PS 把dump命令改成 dump id all atom 10000 dump.lammpstrj 吧，这样方便模拟完了使用VMD观看轨迹动画。 |
5楼2013-01-16 15:20:55
lchunmei001
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1738.5
- 散金: 55
- 红花: 1
- 帖子: 797
- 在线: 151.3小时
- 虫号: 863507
- 注册: 2009-10-05
- 专业: 凝聚态物性I:结构、力学和
6楼2013-01-16 15:30:24
lchunmei001
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1738.5
- 散金: 55
- 红花: 1
- 帖子: 797
- 在线: 151.3小时
- 虫号: 863507
- 注册: 2009-10-05
- 专业: 凝聚态物性I:结构、力学和
7楼2013-01-16 15:34:46
ddll_slsl
木虫 (正式写手)
- 应助: 101 (高中生)
- 金币: 2939.2
- 红花: 6
- 帖子: 670
- 在线: 255.3小时
- 虫号: 2121784
- 注册: 2012-11-12
- 专业: 动力学与控制
8楼2013-01-16 15:56:59
lchunmei001
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1738.5
- 散金: 55
- 红花: 1
- 帖子: 797
- 在线: 151.3小时
- 虫号: 863507
- 注册: 2009-10-05
- 专业: 凝聚态物性I:结构、力学和
9楼2013-01-16 16:12:53
  |
10楼2013-01-16 16:41:01