Diagonalization of the system of Lamé static equations of linear isotropic elasticity

We find the simplest representation of the general solution to the system of static Lamé equations of linear isotropic elasticity in the form of a linear combination of the first derivatives of three functions that satisfy three independent harmonic equations. The representation depends on 12 free parameters choosing which it is possible to obtain various representations of the general solution and simplify the boundary value conditions for the solution of boundary value problems as well as the representation of the general solution for dynamic Lamé equations. The system of Lamé equations diagonalizes, i.e., is reduced to the solution of three independent harmonic equations. The representation implies three conservation laws and a formula for producing new solutions making it possible, given a solution, to find new solutions to the Lamé static equations by derivations. In the two-dimensional case of a plane deformation, the so-found solution immediately implies the Kolosov–Muskhelishvili representation for shifts by means of two analytic functions of complex variable. Two examples are given of applications of the proposed method of diagonalization of two-dimensional elliptic systems.