lammps计算应力

2024-04-14 21:02:45

lammps计算应力

lammps计算的应力有两种：

1.体系整体的应力状态

通过在thermo_style custom里加上pxx pyy pzz pxy pxz pyz字段可将给定时间步（由thermo N命令所指定）的体系应力值输出，再求时间平均即可（实际上求出的是压强张量，即负的应力值）；该应力的计算用的是统计力学里的Virial定理（参见<> by Allen & Tildesley），所算出来的应力与宏观应力是一致的（强调一下，用于平衡态）；

compute 1 all pressure thermo_temp virial

2.单个原子的应力

通过compute stress/atom命令所计算出来的六个值；这个应力的计算则是通过对上述Virial定理所定义的体系应力按原子分解，即将公式中的按原子求和的算符拿掉（同时应将体系的体积换为单个原子的体积），不过，正如楼主所指出的，由于单个原子的体积计算太麻烦，lammps在计算时干脆去掉了体积项，这就是为什么用compute stress/atom命令所算出来的“应力”具有能量的单位的原因。
可以看出，在lammps里，如果要计算体系中某个区域（由region定义，可以是整个模拟盒）所围成的“块”的应力，只需将该区域里的所有原子的单原子应力值加起来，再除以这个区域的体积即可，无须进行单个原子体积的计算。

对原子应力进行Voronoi体积加权平均即可得到系统瞬时应力；系统瞬时应力的系综平均值为宏观测量的系统应力值。

compute s all stress/atom
compute p all reduce sum c_s[1] c_s[2] c_s[3]
variable Press equal -(c_p[1]+c_p[2]+c_p[3])/(3*vol)
thermo_style custom step temp etotal press v_Press

如果要求X方向的应力,就把所有原子用compute stress/atom得到的sigma(xx)加起來,去除掉"所有原子占有的体积",就可以得到系统宏观的sigma(xx).得到的值与別人發表的計算paper相差无几(与实验值相差大约10%)。

转自：http://blog.sina.com.cn/s/blog_c0468c8f0101tfea.html#opennewwindow

lammps计算应力相关推荐

Ubuntu-多用户SSH登陆及LAMMPS计算
前言最近给实验室的Ubuntu服务器进行软件方面的配置,需要实现的功能如下: 实现多用户SSH登陆: 实现外网直接访问服务器(内网穿透): 服务器多用户创建.切换.权限授予: Ubuntu之间.Ub ...
lammps计算聚合物例子_LAMMPS中的系综（NPT/NVT）命令
官方手册的翻译及补充,关于系综命令可见: fix NVT/NPT/NPHlammps.sandia.gov 补充部分在段落结尾以"添注"形式补充. 语法: fix ID grou ...
lammps计算聚合物例子_（从examples的简单例子，到完成自己的科研课题）LAMMPS分子动力学技术与应用专题...
1 分子动力学模拟入门理论--掌握LAMMPS的in文件中各命令的意义 1.1系综理论1.2主要算法介绍1.3积分步长的选取1.4温度和压力控制1.5周期性边界条件1.6分子动力学模拟流程 2 LAM ...
lammps计算聚合物例子_LAMMPS模拟聚合物结构,非晶态聚合物变形行为的模拟,纳米线变形模拟,单轴张力模拟,晶格参数计算...
推荐一个网站,上面有LAMMPS模拟聚合物结构,非晶态聚合物变形行为的模拟,纳米线变形模拟,单轴张力模拟,晶格参数计算的lammps脚本,如下面是晶格参数计算的lammps脚本,具体网址是: LAMM ...
lammps计算聚合物例子_MD原理与Lammps实例.ppt
MD原理与Lammps实例.ppt * Laboratory of Space Materials Science and Technology 分子动力学模拟与 LAMMPS 软件的使用李留辉 2 ...
lammps计算团簇的方法
大家好,我是小马老师. 本文介绍在lammps中如何进行团簇的分析. 在一个group内,原子之间的距离小于一定距离时,可以认为这组原子时一个团簇. 在lammps中,提供了团簇计算命令: compu ...
lammps计算聚合物例子_有了它，实验人员入门计算超容易！材料设计与性质预测平台（MedeA）...
MedeA是全功能材料设计与性质预测平台,对材料及相关领域问题提供快速且可靠的解答,其应用领域非常之广,如电力发电.汽车.储能.合金设计.微电子及石油化工.航天航空.电池等行业. 在MedeA平台中, ...
matlab计算应力位移,2012年-有限元作业-matlab编程实现有限元求解简单结构位移及应力.doc...
<2012年-有限元作业-matlab编程实现有限元求解简单结构位移及应力.doc>由会员分享,可在线阅读,更多相关<2012年-有限元作业-matlab编程实现有限元求解简单结构位 ...
Lammps计算表面粘附力
关注 M r . m a t e r i a l , \color{Violet} \rm Mr.material\ , Mr.material

最新文章

热门文章