计算完成后打开输出的p.profile文件,文件中第一列是bin的编号,第二列是bin的位置,第三列是bin中的平均原子数,第四列到第六列是c_myp[1],c_myp[2] ,c_myp[3]。用第四列乘以第三列再除以bin的体积取相反数就是pxx,用第五列乘以第三列除以bin的体积取相反数就是pyy,用第六列乘以第三列除以bin的体积...
ERROR:Thermo keyword pxx in variable requires thermo to use /init press (../thermos.cpp:1171) 报错具体位置。 解决方法1:将p1、p2、p3去掉。使用thermo_style输出pxx pyy pzz,不推荐使用,画图时不方便。 …
一是体系整体的应力状态,通过在thermo_style custom里加上pxx pyy pzz pxy pxz pyz字段可将给定时间步(由thermo N命令所指定)的体系应力值输出,再求时间平均即可(实际上求出的是压强张量,即负的应力值);该应力的计算用的是统计力学里的Virial定理(参见<<Computer simulation of liquid>> by Allen & Tildesley...
2. 应力的输出 对于全部原子体系的应力,可通过在thermo_style custom里加上pxx pyy pzz pxy pxz pyz字段讲体系应力值输出。注意:pxx pyy pzz这些值为负值。 可以使用下面的公式进行转换,在units metal下,除以10000将应力单位转换为GPa: streesx=-pxx/10000 streesy=-pyy/10000 streesz=-pzz/10000 thermo 100 ...
一是体系整体的应力状态,通过在thermo_style custom里加上pxx pyy pzz pxy pxz pyz字段可将给定时间步(由thermo N命令所指定)的体系应力值输出,再求时间平均即可(实际上求出的是压强张量,即负的应力值);该应力的计算用的是统计力学里的Virial定理(参见<<Computer simulation of liquid>> by Allen & Tildesley...
解决方法1:将p1、p2、p3去掉。使用thermo_style输出pxx pyy pzz,不推荐使用,画图时不方便。 解决方法2:采用另一种计算应力方法,该法适用于在这个计算应力时报这种类型错误使用,适用性单一 解决方法3:适用于所有该类型报错,即由其他变量引起的该类型错误。可在力场后(minimize前)加上thermo命令并输出step press pxx...
Pxx、pyy、pzz、pxy、pyz、pxz、 Fmax、fnorm、 Cella、cellb、cellc、cellalpha、cellbeta、cellgamma、 c_ID、c_ID[I]、c_ID[I][J]、 f_ID、f_ID[I]、f_ID[I][J]、V_name Dump: Dump 命令主要是用来输出计算后单个原子的相关信息的。
1. 计算应变:假如在Z方向上进行拉伸,应变可以这样计算: variable tmp equal "lz" variable L0 equal $ {tmp } variable strain equal "(lz - v_L0)/v_L0"。 2. 输出应力:对于全部原子体系的应力,可通过在thermo_style custom里加上pxx pyy pzz pxy pxz pyz字段将体系应力值输出。注意:pxx pyy pzz...
Pxx、pyy、pzz、pxy、pyz、pxz、 Fmax、fnorm、 Cella、cellb、cellc、cellalpha、cellbeta、cellgamma、 c_ID、c_ID[I]、c_ID[I][J]、 f_ID、f_ID[I]、f_ID[I][J]、V_name Dump: Dump 命令主要是用来输出计算后单个原子的相关信息的。
一是体系整体的应力状态,通过在thermo_style custom里加上pxx pyy pzz pxy pxz pyz字段可将给定时间步(由thermo N命令所指定)的体系应力值输出,再求时间平均即可(实际上求出的是压强张量,即负的应力值); 该应力的计算用的是统计力学里的Virial定理(参见<<Computer simulation of liquid>> by Allen & Tildesle...