石墨烯在LAMMPS中的弹性模量计算

大家好，我想在LAMMPS中通过MD模拟计算单层石墨烯片的弹性常数(C11, C22, C12和C66)。我使用下面的等式:

E = (C11*C22-C12*C21)/C22，泊松比= C12/C22

http://www.sciencedirect.com/science/article/pii/S0261306909003082

y和z方向上的应变都是零。我通过计算应力/原子命令从LAMMPS计算了x、y和xy方向的sigma11、sigma22和tau12。由应力-应变关系求解的方程为:

sigma11 = C11*epsilon11 + C22*epsilon22

sigma22 = C12*epsilon11 + C22*epsilon22

tau12 = C66*gamma12

当22 = 0时，简化为11 = C11*从这个方程我们可以计算出C11。

我的问题是，这里要用到的和的实际值是多少。如果它们是断裂应力和应变，那么我的C11 = 170 GPa。如果我取非常低的应变值，例如0.04，则C11 = 700 GPa。但应该是977 Gpa。

我的计算理解对吗?如果需要，我可以共享我的LAMMPS输入文件。