24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 2735  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

水-天一线

至尊木虫 (职业作家)

[求助] 不知道问题出在哪里,LAMMPS请大家帮我一下(附图片) 已有1人参与

我的目的是:
在一个半径为6.5,高为18.5的圆柱容器中pour 2300个球状颗粒(半径大小为0.5),
这些颗粒只受到重力影响。

我想要的结果是这样的:(可以忽略网格,只看颗粒)
不知道问题出在哪里,LAMMPS请大家帮我一下(附图片)

但是我得到的结果却是这样的:(可以忽略网格,只看颗粒)
不知道问题出在哪里,LAMMPS请大家帮我一下(附图片)-1

问题就出现了,为什么堆积效果不是我想要的结果呢?是不是颗粒之间的作用力需要做一定的限制呢?

这是生成颗粒堆积的代码:
atom_style        sphere
atom_modify     sort 0 0.0
boundary        ff ff ff   
newton                off        
communicate        single vel yes         

region          reg cylinder z 0 0 6.5 0.0 18.5 side in units box
create_box        1 reg

neighbor        0.2  bin   
neigh_modify         delay  0

timestep        0.001

fix                2 all gravity 10.0 spherical 0.0 -180.0  
fix                ins all pour 2300 1 230222 region reg diam one 0.5   

dump            ins all custom 1000 cylinder id x y z  
run             1000

我花在LAMMPS上的时间很长,但是一直都是一个人在瞎摸索,现在也不知道问题出在哪里,希望大家能够帮忙解答一下!
回复此楼

» 猜你喜欢

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

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

水-天一线

至尊木虫 (职业作家)

引用回帖:
2楼: Originally posted by abinitio at 2015-01-26 19:52:33
看堆积后结果,是不是重力方向反了?另外的原子变得无序,是不是加上了温度?还有,请教一下,你的这个voronoi是用哪个程序算的?voro++吗?

重力方向应该没错,虽然上面一部分非常紧密,但是中间却很松散。

原子虽然无序,但是我也没有添加温度。

我只是觉得原子之间的受力是不是要考虑进去?

voronoi是用voro++算的。
3楼2015-01-26 19:58:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

abinitio

版主 (著名写手)


水-天一线(月只蓝代发): 金币+1, 感谢指导 2015-01-27 21:39:25
看堆积后结果,是不是重力方向反了?另外的原子变得无序,是不是加上了温度?还有,请教一下,你的这个voronoi是用哪个程序算的?voro++吗?

[ 发自小木虫客户端 ]
2楼2015-01-26 19:52:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abinitio

版主 (著名写手)

引用回帖:
3楼: Originally posted by 水-天一线 at 2015-01-26 19:58:43
重力方向应该没错,虽然上面一部分非常紧密,但是中间却很松散。

原子虽然无序,但是我也没有添加温度。

我只是觉得原子之间的受力是不是要考虑进去?

voronoi是用voro++算的。...

没看到你的势函数是如何加上的

[ 发自小木虫客户端 ]
4楼2015-01-26 21:19:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

水-天一线

至尊木虫 (职业作家)

引用回帖:
4楼: Originally posted by abinitio at 2015-01-26 21:19:29
没看到你的势函数是如何加上的
...

的确,我没有添加势函数,我只是用了一句“fix    2 all gravity 10.0 spherical 0.0 -180.0  ”。
其实我想表达的意思是:
我想用lammps模拟现实的颗粒堆积过程,好比把N个大小一样的球放到固定尺寸的圆柱容器中,每个球体与周围球体都是相互接触的。

由于原本不是学这块内容的,所以不清楚怎么设置这方面的参数。
5楼2015-01-26 21:38:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 297求调剂 +22 ORCHID1 2026-04-10 25/1250 2026-04-13 23:11 by pies112
[考研] 085600材料与化工,求调剂 +14 won_qii 2026-04-07 14/700 2026-04-13 22:21 by pies112
[考研] 化学070300 求调剂 +21 哈哈哈^_^ 2026-04-12 21/1050 2026-04-13 21:36 by nxybio2007
[考研] 一志愿华工085600 331分 +7 天下ww 2026-04-09 7/350 2026-04-13 09:01 by lhj2009
[考研] 2本,初试303,0860求调剂 +6 floriea 2026-04-12 8/400 2026-04-12 18:13 by zhouxiaoyu
[考研] 291求调剂 +11 关忆北. 2026-04-09 12/600 2026-04-12 10:32 by 逆水乘风
[考研] 求调剂,262机械专硕 +8 嗯yyl 2026-04-08 8/400 2026-04-12 02:31 by 秋豆菜芽
[考研] 22408调剂315分 +3 zhuangyan123 2026-04-09 3/150 2026-04-12 00:25 by 蓝云思雨
[考研] 070300化学279求调剂 +19 哈哈哈^_^ 2026-04-08 20/1000 2026-04-11 20:43 by stoner78
[考研] 085410 273分调剂 +4 X1999 2026-04-09 4/200 2026-04-11 13:05 by pies112
[考研] 化学308分求调剂 +22 你好明天你好 2026-04-07 24/1200 2026-04-11 11:14 by ChemPharm
[考研] 农业管理302分求调剂 +3 xuening1 2026-04-10 3/150 2026-04-11 10:18 by zhq0425
[考研] 346,工科0854求调剂,专硕 +7 moser233 2026-04-10 8/400 2026-04-11 08:52 by 猪会飞
[考研] 287求调剂 +15 Fnhc 2026-04-07 21/1050 2026-04-10 19:09 by chemisry
[考研] 071000生物学调剂求助 +17 zzzzwww 2026-04-09 20/1000 2026-04-10 15:55 by 求调剂zz
[考研] 301求调剂 +5 149. 2026-04-10 5/250 2026-04-10 15:45 by 柴小白
[考研] 材料专硕调剂 +16 哈哈哈吼吼吼哈 2026-04-07 17/850 2026-04-09 21:16 by wutongshun
[考研] 348求调剂 +3 candyyyi 2026-04-09 3/150 2026-04-09 17:20 by 段伟艳
[考研] 337求调剂 +4 Gky09300550, 2026-04-09 4/200 2026-04-09 17:18 by 帕尔马拉特
[考研] 085501机械英二77总分294求调剂,接受跨专业学习 +6 守法公民亓纪 2026-04-08 6/300 2026-04-09 15:55 by wp06
信息提示
请填处理意见