| 查看: 2831 | 回复: 10 | ||
2009云飞扬铁杆木虫 (著名写手)
|
[求助]
lammps中某个区域压强的求解
|
|
|
各位朋友，还是前面提到过的一个问题。我建立了一个立方体，中心挖去原子形成一个球孔，在球孔中填入了另一种原子A，想求解球孔区域A原子所形成的压强是多大，目前知道应该局部压强采用compute stress/atom 来计算（单位是能量——压强.体积），但是这样就需要求解球孔区域A原子所形成的体积有多大，计算体积目前是采用的compute voronoi/atom 指令。但是在计算过程中总是觉得体积算的不对，相关的具体化指令如下，请知道的朋友帮忙指点指点，应该怎么计算才能正确得到。先谢谢啦。 compute stress voidA stress/atom He_bubble compute stress1 voidA reduce sum c_stress[1] compute stress2 voidA reduce sum c_stress[2] compute stress3 voidA reduce sum c_stress[3] compute volum voidA voronoi/atom only_group variable p equal -(${c_stress1}+${c_stress2}+${c_stress3})/（3.0*c_volum） 体积计算部分是不是有问题呢？希望知道的朋友能够详细说说，先谢谢啦 |
回复此楼» 猜你喜欢
全日制（定向）博士
已经有5人回复
-
假如你的研究生提出不合理要求
已经有10人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
Nater1ver
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 3339.1
- 红花: 14
- 帖子: 737
- 在线: 47.7小时
- 虫号: 1526557
- 注册: 2011-12-06
- 专业: 遥感机理与方法
| 印象中compute voronoi/atom是唯一计算体系部分原子体积的方法，这个算法有些思路清奇，严格上说可能不太准确，可我用过一次，觉得精度还行。但在用之前你总要知道这个命令输出了什么东西吧。刚才简单看了一下手册，我有些记不清了，只能大致说一下当时是怎么做的。1)voronoi体积不是真正的体积，所以你是不是要V(总体积)*voronoi(A原子)/voronoi(所有原子)来得到A原子的体积？ 2)我记得这个命令输出的不是单个的数，而是数组，好像还是二维数组，所以你是不是应该输出到一个文本看看它的输出格式，把其中一个数组所有的结果加起来，甚至为了结果准确，还要做一个系综平均 |
9楼2019-03-12 14:36:24
2009云飞扬
铁杆木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 6293.2
- 散金: 99
- 红花: 2
- 帖子: 1767
- 在线: 363.9小时
- 虫号: 1745038
- 注册: 2012-04-09
- 专业: 固体力学
2楼2019-01-24 14:00:32
2009云飞扬
铁杆木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 6293.2
- 散金: 99
- 红花: 2
- 帖子: 1767
- 在线: 363.9小时
- 虫号: 1745038
- 注册: 2012-04-09
- 专业: 固体力学
3楼2019-01-24 18:47:04
2009云飞扬
铁杆木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 6293.2
- 散金: 99
- 红花: 2
- 帖子: 1767
- 在线: 363.9小时
- 虫号: 1745038
- 注册: 2012-04-09
- 专业: 固体力学
4楼2019-01-26 11:22:38
ruibingqilin
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 682.7
- 帖子: 285
- 在线: 152.6小时
- 虫号: 1943683
- 注册: 2012-08-17
- 性别: GG
- 专业: 高分子材料结构与性能
5楼2019-02-12 11:49:08
吾艾高
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 368.9
- 散金: 10
- 帖子: 88
- 在线: 2.5小时
- 虫号: 13099613
- 注册: 2018-12-08
- 性别: GG
- 专业: 爆炸与冲击动力学
6楼2019-03-11 23:54:13
2009云飞扬
铁杆木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 6293.2
- 散金: 99
- 红花: 2
- 帖子: 1767
- 在线: 363.9小时
- 虫号: 1745038
- 注册: 2012-04-09
- 专业: 固体力学
7楼2019-03-12 13:51:26
Nater1ver
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 3339.1
- 红花: 14
- 帖子: 737
- 在线: 47.7小时
- 虫号: 1526557
- 注册: 2011-12-06
- 专业: 遥感机理与方法
8楼2019-03-12 14:35:55
2009云飞扬
铁杆木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 6293.2
- 散金: 99
- 红花: 2
- 帖子: 1767
- 在线: 363.9小时
- 虫号: 1745038
- 注册: 2012-04-09
- 专业: 固体力学
10楼2019-03-12 23:33:06