|
大家好,我前一段时间发现这个论坛有些朋友对于如何在原子模拟中计算应力比较困惑。由于我个人担任原子模拟课程的讲授,所以对这个问题有一定的思考,希望在这里与大家进行交流。我对该问题的认识可以参考我最近的一篇相关的英文文章,发表在《Journal of Computational and Theoretical Nanoscience》2009, vol.6(5):1081-1089上的论文,见附件。下面是其摘要的中文版。另外为方便大家阅读,附件还给出一个简短的中文版文件,
原子模拟目前已成为力学中的有力研究工具之一,但是如何能合理地从中获取与柯西应力等效的应力在国际范围内还是很有争议,特别是在如何处理virial应力定义中的速度相关项。本文通过严格地遵循连续介质力学中柯西应力的定义,提出了拉格朗日描述下原子应力的基本定义,并可在任何情形下正确计算得到柯西应力。为了便于计算,又进一步提出拉格朗日virial应力,该定义仍然具有virial的形式,但是避免包含了有争议的速度项。我们还指出被广泛采用的经典的virial应力其实是欧拉描述下的virial应力,其中包含了原子的运动速度项,该定义只有在可用冲量-动量定理能正确计算内力时才适用。但在实际原子模拟中,很多时候这个条件并不能满足,如本文举出的一些例子,而我们提出拉格朗日原子应力可以有效避免这种速度相关的不客观性。
[ 本帖最后由 liubcn666 于 2009-3-11 14:54 编辑 ]
|
|