Abstract Generated abstract
This paper derives formulas connecting the coefficients of a generalized elastic energy density, including dependence on strain gradients, with microscopic force constants of a crystal lattice. Starting from the Born model of a simple lattice in the harmonic approximation, the discrete lattice description is interpolated to a continuum representation and expressed in Fourier space through tensors relating particle displacements, strains, and elastic moduli. By expanding the resulting elastic tensor in powers of the wave vector, the work obtains explicit expressions for the continuum elastic constants of all gradient orders in terms of lattice sums over the force constants. The zeroth order term is identified with the ordinary elastic modulus tensor in the long wavelength approximation.
Full Text
UDC 539.31
THEORY OF ELASTICITY
S. S. KUKUSHKINA
COMPUTATION OF ELASTIC CONSTANTS FROM MICROSCOPIC QUANTITIES FOR AN ELASTIC MEDIUM WHOSE ENERGY DEPENDS ON STRAIN GRADIENTS
(Presented by Academician L. I. Sedov, 13 XI 1967)
In a number of works (see, for example, \((^{1-4})\)) a generalized theory of elasticity is considered in which the elastic energy is represented in the form of a quadratic form in the strains and the strain gradients with respect to the coordinates. In the classical theory of elasticity the density of elastic energy \(U\) is represented in the form
\[ U = {}^{1}/_{2} c^{iklm}\varepsilon_{ik}\varepsilon_{lm}, \tag{1} \]
where the elastic constants \(c^{iklm}\) can be expressed in terms of the force constants of the crystal lattice. In the present work formulas are obtained which express the coefficients of the quadratic form of the elastic-energy density in the generalized theory of elasticity in terms of the force constants of the crystal lattice.
In passing from the micromodel to the continuum we shall use the method developed in \((^{5,6})\). As the initial micromodel we take the Born model of a simple crystal lattice in the harmonic approximation \((^7)\). The potential energy of the lattice is taken in the form
\[ \Phi = {}^{1}/_{2}\sum_{n,n'} u_{\alpha}(n)\Phi^{\alpha\beta}(n-n')u_{\beta}(n'), \tag{2} \]
where \(n=n^{\alpha}e_{\alpha}\) is a vector in an oblique lattice coordinate system with basis \(e_{\alpha}\) (the lattice sites correspond to integral components \(n^{\alpha}\)), \(u(n)\) are the displacements of the particles, and \(\Phi^{\alpha\beta}(n-n')\) are the force constants.
It was shown in \((^{5,6})\) that, by means of a certain interpolation, one can pass from functions of the discrete argument \(n\) to analytic functions of \(x\), with \(\Phi^{\alpha\beta}(n)\) replaced by the function \(\Phi^{\alpha\beta}(x)\), whose Fourier transform is
\[ \Phi^{\alpha\beta}(k)=\frac{1}{V_A}\sum_n \Phi^{\alpha\beta}(n)e^{in\cdot k}; \tag{3} \]
here \(V_A\) is the volume of the elementary cell of the lattice. Then the potential energy \(\Phi\) is written in the form
\[ \begin{aligned} 2\Phi &= \iint u_{\alpha}(x)\Phi^{\alpha\beta}(x-x')u_{\beta}(x')\,dx\,dx' = \frac{1}{(2\pi)^3}\int \overline{u_{\alpha}}(k)\Phi^{\alpha\beta}(k)u_{\beta}(k)\,dk \\ &= \iint \varepsilon_{\nu\alpha}(x)c^{\nu\alpha\mu\beta}(x-x')\varepsilon_{\mu\beta}(x')\,dx\,dx' = \frac{1}{(2\pi)^3}\int \overline{\varepsilon_{\nu\alpha}}(k)c^{\nu\alpha\mu\beta}(k)\varepsilon_{\mu\beta}(k)\,dk, \end{aligned} \tag{4} \]
where \(\varepsilon_{\nu\alpha}(x)=\partial_{(\nu}u_{\alpha)}(x)\) is the strain. Here a Fourier transformation with respect to \(x\) has been carried out.
The tensor \(c^{\nu\alpha\mu\beta}\) is symmetric in the pairs \(\nu\alpha\) and \(\mu\beta\) and Hermitian with respect to interchange of the pairs.
It was shown in \((^5)\) that \(\Phi^{\alpha\beta}(k)\) can be represented in the form
\[ \Phi^{\alpha\beta}(k)=\psi^{\alpha\beta\nu\mu}(k)k_{\nu}k_{\mu}, \tag{5} \]
where the tensor \(\psi^{\alpha\beta\nu\mu}\) is symmetric in \(\nu\mu\) and Hermitian in \(\alpha\beta\). The tensor \(c^{\nu\alpha\mu\beta}(k)\) is constructed as follows
\[ c^{\nu\alpha\mu\beta} = {}^{1}/_{2}\bigl(a^{\nu\alpha(\mu\beta)}+\overline{a^{\mu\beta(\nu\alpha)}}\bigr) + k_\lambda\bigl(a^{\nu\alpha\lambda(\mu\beta)} +\overline{a^{\mu\beta\lambda(\nu\alpha)}}\bigr), \tag{6} \]
here
\[ a^{\nu\alpha\mu\beta}(k) = \psi^{\alpha\beta\nu\mu}(k) + \psi^{\nu\beta\alpha\mu}(k) - \psi^{\mu\beta\alpha\nu}(k), \tag{7} \]
\[ a^{\nu\alpha[\mu\beta]}(k) = k_\lambda a^{\nu\alpha[\mu\beta]\lambda}(k), \]
the tensor \(a^{\nu\alpha\mu\beta\lambda}\) being symmetric in \(\nu\lambda\).
Represent \(\psi^{\alpha\beta\nu\mu}(k)\) in the form
\[ \psi^{\alpha\beta\nu\mu}(k) = \psi_{0}^{\alpha\beta\nu\mu} + k_\lambda \psi^{\alpha\beta\nu\mu\lambda}(k) \tag{8} \]
(here the subscript zero denotes the value of the function at \(k=0\)). Then \(a^{\nu\alpha\mu\beta}(k)\) can be represented in the form
\[ a^{\nu\alpha\mu\beta}(k) = \psi_{0}^{\alpha\beta\nu\mu} + \psi_{0}^{\nu\beta\alpha\mu} - \psi_{0}^{\mu\beta\alpha\nu} + k_\lambda\bigl(\psi^{\alpha\beta\nu\mu\lambda}(k) + \psi^{\nu\beta\alpha\mu\lambda}(k) - \psi^{\mu\beta\alpha\nu\lambda}(k)\bigr) = a_{0}^{\nu\alpha\mu\beta} + k_\lambda a^{\nu\alpha\mu\beta\lambda}. \tag{9} \]
Substituting (7) and (9) into (6) and taking into account the symmetry properties of \(\psi^{\alpha\beta\nu\mu}\), we obtain an expression for \(c^{\nu\alpha\mu\beta}\) in terms of \(\psi_{0}^{\alpha\beta\nu\mu}\) and \(\psi^{\alpha\beta\nu\mu\lambda}(k)\):
\[ c^{\nu\alpha\mu\beta}(k) = S_{0\chi\sigma\rho\tau}^{\nu\alpha\mu\beta} \psi_{0}^{\chi\sigma\rho\tau} + {}^{1}/_{4}k_\lambda \bigl[ S_{1\chi\sigma\rho\tau}^{\nu\alpha\mu\beta} \psi^{\lambda\chi\sigma\rho\tau}(k) + S_{2\chi\sigma\rho\tau}^{\nu\alpha\mu\beta} \psi^{\chi\lambda\sigma\rho\tau}(k) + S_{3\chi\sigma\rho\tau}^{\nu\alpha\mu\beta} \psi^{\chi\sigma\rho\lambda\tau}(k) \bigr]. \tag{10} \]
Here the following symmetrization operators have been introduced:
\[ S_{0\chi\sigma\rho\tau}^{\nu\alpha\mu\beta} = \delta_{\chi\sigma\rho\tau}^{\alpha\beta\nu\mu} + \delta_{\chi\sigma\rho\tau}^{\nu\beta\alpha\mu} - \delta_{\chi\sigma\rho\tau}^{\mu\beta\alpha\nu}, \]
\[ S_{1\chi\sigma\rho\tau}^{\nu\alpha\mu\beta} = \delta_{\chi\sigma\rho\tau}^{\nu\mu\beta\alpha} + \delta_{\chi\sigma\rho\tau}^{\alpha\beta\mu\nu} - 3\delta_{\chi\sigma\rho\tau}^{\beta\alpha\nu\mu} - 3\delta_{\chi\sigma\rho\tau}^{\mu\alpha\nu\beta}, \]
\[ S_{2\chi\sigma\rho\tau}^{\nu\alpha\mu\beta} = \delta_{\chi\sigma\rho\tau}^{\beta\mu\nu\alpha} + \delta_{\chi\sigma\rho\tau}^{\lambda\alpha\nu\beta} - 3\delta_{\chi\sigma\rho\tau}^{\alpha\nu\beta\mu} - 3\delta_{\chi\sigma\rho\tau}^{\nu\alpha\beta\mu}, \]
\[ S_{3\chi\sigma\rho\tau}^{\nu\alpha\mu\beta} = 4\delta_{\chi\sigma\rho\tau}^{\alpha\beta\nu\mu} + 4\delta_{\chi\sigma\rho\tau}^{\nu\beta\alpha\mu} + 4\delta_{\chi\sigma\rho\tau}^{\alpha\mu\nu\beta} + 4\delta_{\chi\sigma\rho\tau}^{\nu\mu\alpha\beta} - \]
\[ - \delta_{\chi\sigma\rho\tau}^{\mu\beta\alpha\nu} - \delta_{\chi\sigma\rho\tau}^{\mu\alpha\alpha\nu} - \delta_{\chi\sigma\rho\tau}^{\alpha\nu\beta\mu} - \delta_{\chi\sigma\rho\tau}^{\nu\alpha\beta\mu}, \]
where
\[ \delta_{\chi\sigma\rho\tau}^{\alpha\beta\nu\mu} = \delta_{\chi}^{\alpha}\delta_{\sigma}^{\beta}\delta_{\rho}^{\nu}\delta_{\tau}^{\mu}. \]
Expand the analytic function \(c^{\nu\alpha\mu\beta}(k)\) in a series in \(k\):
\[ c^{\nu\alpha\mu\beta}(k) = \sum_{p=0}^{\infty} c_{p}^{\nu\alpha\mu\beta\lambda_{1}\ldots\lambda_{p}} (-ik_{\lambda_{1}})\cdots(-ik_{\lambda_{p}}) \tag{11} \]
and carry out the same expansion under the integral sign in (4). Then the expression for the energy assumes the form:
\[ 2\Phi = \frac{1}{(2\pi)^3} \int \overline{\varepsilon_{\nu\alpha}(k)} c^{\nu\alpha\mu\beta}(k) \varepsilon_{\mu\beta}(k)\,dk = \]
\[ = \frac{1}{(2\pi)^3} \int \overline{\varepsilon_{\nu\alpha}(k)} \bigl( c_{0}^{\nu\alpha\mu\beta} + c_{1}^{\nu\alpha\mu\beta\lambda_{1}}(-ik_{\lambda_{1}}) + c_{2}^{\nu\alpha\mu\beta\lambda_{1}\lambda_{2}}(-ik_{\lambda_{1}})(-ik_{\lambda_{2}}) +\ldots \bigr) \times \]
\[ \times \varepsilon_{\mu\beta}(k)\,dk = \int \bigl( c_{0}^{\nu\alpha\mu\beta} \varepsilon_{\nu\alpha}(x)\varepsilon_{\mu\beta}(x) + c_{1}^{\nu\alpha\mu\beta\lambda_{1}} \varepsilon_{\nu\alpha}(x)\partial_{\lambda_{1}}\varepsilon_{\mu\beta}(x) - \]
\[ - c_{2}^{\nu\alpha\mu\beta\lambda_{1}\lambda_{2}} \partial_{\lambda_{1}}\varepsilon_{\nu\alpha}(x) \partial_{\lambda_{2}}\varepsilon_{\mu\beta}(x) +\ldots \bigr]\,dx. \tag{12} \]
For the density of elastic energy \(\varphi(x)\) we obtain the expression
\[ \varphi(x) = {}^{1}/_{2} \bigl( c_{0}^{\nu\alpha\mu\beta}\varepsilon_{\nu\alpha}\varepsilon_{\mu\beta} + c_{1}^{\nu\alpha\mu\beta\lambda_{1}} \varepsilon_{\nu\alpha}\partial_{\lambda_{1}}\varepsilon_{\mu\beta} +\ldots \]
\[ \ldots + (-1)^{m} c_{2m}^{\nu\alpha\mu\beta\lambda_{1}\ldots\lambda_{2m}} \partial_{\lambda_{1}}\cdots\partial_{\lambda_{m}}\varepsilon_{\nu\alpha} \partial_{\lambda_{m+1}}\cdots\partial_{\lambda_{2m}}\varepsilon_{\mu\beta} + \]
\[ + (-1)^{m} c_{2m+1}^{\nu\alpha\mu\beta\lambda_{1}\ldots\lambda_{2m+1}} \partial_{\lambda_{1}}\cdots\partial_{\lambda_{m}}\varepsilon_{\nu\alpha} \partial_{\lambda_{m+1}}\cdots\partial_{\lambda_{2m+1}}\varepsilon_{\mu\beta} +\ldots \bigr). \tag{13} \]
Let us expand the functions \(\psi^{\nu\alpha\mu\beta\lambda}\) in (10) also in a series in \(k\). Taking (3), (5), (8), and (10) into account, we obtain an expression for the elastic constants in (13) in terms of the force constants of the lattice:
\[ c_0^{\nu\alpha\mu\beta} = -\frac{1}{2V_A}\, S_{0\chi\sigma\rho\tau}^{\nu\alpha\mu\beta} \sum_n \Phi^{\chi\sigma}(n)n^\rho n^\tau, \]
\[ \begin{aligned} c_p^{\nu\alpha\mu\beta\lambda_1\cdots\lambda_p} &= \frac{(-1)^{p+1}}{4V_A(p+2)!} \Bigg\{ S_{1\chi\sigma\rho\tau}^{\nu\alpha\mu\beta} \sum_n \Phi^{\lambda_p\chi}(n) n^\sigma n^\rho n^\tau n^{\lambda_1}\cdots n^{\lambda_{p-1}} \\ &\quad+ S_{2\chi\sigma\rho\tau}^{\nu\alpha\mu\beta} \sum_n \Phi^{\chi\lambda_p}(n) n^\sigma n^\rho n^\tau n^{\lambda_1}\cdots n^{\lambda_{p-1}} \\ &\quad+ S_{3\chi\sigma\rho\tau}^{\nu\alpha\mu\beta} \sum_n \Phi^{\chi\sigma}(n) n^\rho n^\tau n^{\lambda_1}\cdots n^{\lambda_p} \Bigg\}. \end{aligned} \tag{14} \]
In the zero long-wavelength approximation, \(c_0^{\nu\alpha\mu\beta}\) is the ordinary tensor of elastic moduli.
Mathematical Institute named after V. A. Steklov
Academy of Sciences of the USSR
Received
15 X 1967
References Cited
\({}^{1}\) E. L. Aero, E. V. Kuvshinskii, FTT, 5, No. 9, 2591 (1963).
\({}^{2}\) R. D. Mindlin, Arch. Rat. Mech., Anal., 16, 51 (1964).
\({}^{3}\) E. Kröner, Int. J. Eng. Sci., 1, 261 (1963).
\({}^{4}\) M. A. Indin, PMM, 30, issue 3, 531 (1966).
\({}^{5}\) I. A. Kunin, PMM, 30, issue 3, 542 (1966).
\({}^{6}\) V. E. Vdovin, I. A. Kunin, PMM, 30, issue 6, 1071 (1966).
\({}^{7}\) M. Born, Huan Kun, Dynamical Theory of Crystal Lattices, IL, 1958.