Construction of the a/n/l problem
We first extend our domain with a region which we denote Then we consider the linearisation of the Cauchy Born model in this region. As what follows requires indices, we briefly employ indicial summation convention. Let be a reference strain, then
thus for our problem:
Where is known as the shear modulus of linear elasticity. If we denote the fourth order tensor , then for a small displacement , we obtain the linearised energy difference functional:
In particular, in our project we are concerned with 1D motion in an antiplane fashion, (i.e displacements occur perpendicular to the atomistic lattice). Under this assumption, which implies here, and if we use a potential uses nearest neighbour interactions, the linearised energy difference functional becomes
Due to the fundemental local definition of the continuum models, there are no long range interactions that are taken into account, and thus no force imbalances over the continuum interface. So coupling the energy is simply acheived by summing the energies in each domain, Hence:
then the A/N/L coupling problem seeks to find
Consistency errors and convergence rates
We now state the construction of the analytical error estimates in seminorm between this model and the ideal atomistic model. The consistency analysis is straightforward as we are able to split up the consistency difference into two parts:
The first of these is considered in the Atomistic/Nonlinear coupling section. The second term is then dealt with an application of the taylors theorem to yield.
From which we can obtain a rate of convergence based on the radius of the atomistic domain for our point defect (which satisfies ):
So, as in the nonlinear elasticity case the boundary is the biggest source of error. In fact, our numerical tests found there is no practical difference between the nonlinear and linear elasticity minimisers in the case of this defect.