[论文解读] The role of the patch test in 2D atomistic-to-continuum coupling methods

Christoph Ortner|arXiv (Cornell University)|Jan 27, 2011

Microstructure and mechanical properties参考文献 37被引用 8

一句话总结

该论文证明，在四个关键技术条件下——能量一致性、局部化界面校正、校正的体积缩放以及原子区域的连通性——二维原子-连续（a/c）耦合方法中消除鬼力的片测试一致性，可保证一阶精度。该分析通过新颖的应力张量构造和 W^{1,p} 范数下的变分罪行技术，为基于能量的 a/c 方法提供了严格的理论基础，将一致性与收敛速率联系起来。

ABSTRACT

For a general class of atomistic-to-continuum coupling methods, coupling multi-body interatomic potentials with a P1-finite element discretisation of Cauchy--Born nonlinear elasticity, this paper adresses the question whether patch test consistency (or, absence of ghost forces) implies a first-order error estimate. In two dimensions it is shown that this is indeed true under the following additional technical assumptions: (i) an energy consistency condition, (ii) locality of the interface correction, (iii) volumetric scaling of the interface correction, and (iv) connectedness of the atomistic region. The extent to which these assumptions are necessary is discussed in detail.

研究动机与目标

• 确定二维 a/c 方法中的片测试一致性是否意味着一阶误差估计。
• 弥合能量基 a/c 方法理论理解中的关键空白，此前一致性尚未与收敛性建立联系。
• 确立在何种充分条件下，无鬼力方法可实现最优收敛速率。
• 基于 W^{1,p} 范数与变分罪行技术，建立 a/c 方法误差估计的一般分析框架。

提出的方法

• 引入一类通用的 a/c 方法，将多体原子间势能与 P1 有限元方法耦合，用于柯西-伯恩弹性模型。
• 开发一种专为二维和三维 a/c 耦合设计的新符号与分析框架，包括用于衡量局部平滑度的离散振荡算子。
• 应用 Shapeev 的键密度引理，将键积分与体积分关联，实现应力张量的比较。
• 构建类似于维里应力的原子与 a/c 应力函数，通过应力一致性实现误差分解。
• 采用改进的 a/c 应力函数与校正函数，控制在存在界面校正时的建模误差。
• 应用插值误差估计与对偶范数技术，以 W^{-1,p} 范数界定向建模误差，从而隔离关键的稳定性挑战。

实验结果

研究问题

• RQ1 在二维 a/c 方法中，片测试一致性是否意味着基于能量的公式具有第一阶收敛性？
• RQ2 为确保鬼力消除可导致最优误差估计，还需哪些额外的技术假设？
• RQ3 如何构建原子与 a/c 模型的应力张量，以实现严谨的误差分析？
• RQ4 能量一致性、局部化校正、体积缩放与原子区域连通性等假设，在实现第一阶收敛性方面，其必要程度如何？
• RQ5 建模误差能否在不依赖网格粗化的情况下被界定向？在何种条件下可以实现？

主要发现

• 仅靠片测试一致性不足以保证一阶收敛性；必须结合四项额外技术假设。
• 若方法满足能量一致性、局部化界面校正、校正的体积缩放以及原子区域的连通性，则在二维中可实现一阶一致性。
• 原子与 a/c 应力函数的构建，使建模误差可被精确分解为应力一致性项与校正项。
• 采用改进的 a/c 应力函数可控制界面校正存在时的误差，其中校正函数确保了界面处的一致性。
• 分析表明，只要满足稳定性假设，建模误差可独立于网格粗化被界定向。
• 主要结果，定理 6，证明在所述假设下，a/c 方法在 W^{1,p} 范数下实现了第一阶一致性。