| 查看: 3576 | 回复: 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） 体积计算部分是不是有问题呢？希望知道的朋友能够详细说说，先谢谢啦 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有14人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有16人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有15人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有17人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有21人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有21人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
Nater1ver
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 3501.2
- 红花: 14
- 帖子: 772
- 在线: 48小时
- 虫号: 1526557
- 注册: 2011-12-06
- 专业: 水文、水资源
| 印象中compute voronoi/atom是唯一计算体系部分原子体积的方法，这个算法有些思路清奇，严格上说可能不太准确，可我用过一次，觉得精度还行。但在用之前你总要知道这个命令输出了什么东西吧。刚才简单看了一下手册，我有些记不清了，只能大致说一下当时是怎么做的。1)voronoi体积不是真正的体积，所以你是不是要V(总体积)*voronoi(A原子)/voronoi(所有原子)来得到A原子的体积？ 2)我记得这个命令输出的不是单个的数，而是数组，好像还是二维数组，所以你是不是应该输出到一个文本看看它的输出格式，把其中一个数组所有的结果加起来，甚至为了结果准确，还要做一个系综平均 |
9楼2019-03-12 14:36:24
2009云飞扬
铁杆木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 6404.2
- 散金: 99
- 红花: 3
- 帖子: 1797
- 在线: 365.1小时
- 虫号: 1745038
- 注册: 2012-04-09
- 专业: 固体力学
2楼2019-01-24 14:00:32
2009云飞扬
铁杆木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 6404.2
- 散金: 99
- 红花: 3
- 帖子: 1797
- 在线: 365.1小时
- 虫号: 1745038
- 注册: 2012-04-09
- 专业: 固体力学
3楼2019-01-24 18:47:04
2009云飞扬
铁杆木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 6404.2
- 散金: 99
- 红花: 3
- 帖子: 1797
- 在线: 365.1小时
- 虫号: 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 (幼儿园)
- 金币: 6404.2
- 散金: 99
- 红花: 3
- 帖子: 1797
- 在线: 365.1小时
- 虫号: 1745038
- 注册: 2012-04-09
- 专业: 固体力学
7楼2019-03-12 13:51:26
Nater1ver
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 3501.2
- 红花: 14
- 帖子: 772
- 在线: 48小时
- 虫号: 1526557
- 注册: 2011-12-06
- 专业: 水文、水资源
8楼2019-03-12 14:35:55
2009云飞扬
铁杆木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 6404.2
- 散金: 99
- 红花: 3
- 帖子: 1797
- 在线: 365.1小时
- 虫号: 1745038
- 注册: 2012-04-09
- 专业: 固体力学
10楼2019-03-12 23:33:06
20