24小时热门版块排行榜

返回列表

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

2009云飞扬

铁杆木虫 (著名写手)

应助: 1 (幼儿园)
金币: 6327.2
散金: 99
红花: 2
帖子: 1752
在线: 363小时
虫号: 1745038
注册: 2012-04-09
专业: 固体力学

[求助] 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）
体积计算部分是不是有问题呢？希望知道的朋友能够详细说说，先谢谢啦

回复此楼

» 猜你喜欢

请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
2025冷门绝学什么时候出结果已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

1楼2019-01-23 15:46:17

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

Nater1ver

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 3258.6
红花: 14
帖子: 719
在线: 46.8小时
虫号: 1526557
注册: 2011-12-06
专业: 遥感机理与方法

印象中compute voronoi/atom是唯一计算体系部分原子体积的方法，这个算法有些思路清奇，严格上说可能不太准确，可我用过一次，觉得精度还行。但在用之前你总要知道这个命令输出了什么东西吧。刚才简单看了一下手册，我有些记不清了，只能大致说一下当时是怎么做的。1)voronoi体积不是真正的体积，所以你是不是要V(总体积)*voronoi(A原子)/voronoi(所有原子)来得到A原子的体积？ 2)我记得这个命令输出的不是单个的数，而是数组，好像还是二维数组，所以你是不是应该输出到一个文本看看它的输出格式，把其中一个数组所有的结果加起来，甚至为了结果准确，还要做一个系综平均