Рассмотрен процесс построения численно-аналитического решения основной смешанной краевой задачи статической теории упругости (ТУ). Определяющим соотношениям линейной однородной изотропной упругой среды эквивалентна система трех разрешающих уравнений Ламе, каждое — линейное второго порядка относительно частных производных. Средствами энергетического метода граничных состояний (МГС) с помощью разложения искомого состояния в ряд Фурье по элементам сепарабельного базиса гильбертова пространства состояний краевая задача математической физики приводится к системе линейных алгебраических уравнений (СЛАУ). Построение базиса опирается на вариант общего решения системы уравнений эллиптического типа. Каждое базисное состояние создается посредством использования гармонических многочленов. Ортогонализация проводится алгоритмом Грама–Шмидта. Граничные условия (ГУ) основной смешанной задачи предполагают разбиение границы тела на два класса. На одной части границы заданы перемещения точек, на оставшейся — поверхностные усилия. СЛАУ содержит информацию как о содержимом базиса (матрица коэффициентов), так и о значениях характеристик ГУ (вектор правых частей). Точность решения оценивается двумя факторами: 1) неравенство Бесселя; 2) интегральная квадратичная невязка восстановленного граничного состояния с ГУ. Использование обоих факторов приводит к эффекту самодостаточности МГС: нет необходимости в сопоставлении пробных решений с эталонными, построенными иными методами. Выполнены конкретные расчеты для двух классов смешанных задач: 1) ограниченное двуполостное тело (шар с двумя симметрично расположенными сферическими полостями. Задано жесткое встречное смещение поверхностей полостей, внешняя граница тела свободна от нагрузки; 2) упругий круговой цилиндр закреплен с жестким штоком по части одного основания. Боковая поверхность нагружена касательными усилиями, направленными вдоль оси цилиндра. Во второй задаче испытывается сингулярный скачок характеристик граничного состояния по линии раздела ГУ, принадлежащей границе тела. Построены численно-аналитические решения обеих задач, проиллюстрированы поля характеристик НДС тела, выполнены комментарии состояний и сделаны выводы. 

The process of constructing a numerical-analytical solution to the basic mixed boundary value problem of the static theory of elasticity (TE) is considered. The defining relations of a linear homogeneous isotropic elastic medium are equivalent to a system of three resolving Lame equations, each linear of the second order with respect to partial derivatives. By means of the energetic method of boundary states (MBS), by decomposing the desired state into a Fourier series according to the elements of the separable basis of the Hilbert space of states, the boundary value problem of mathematical physics is reduced to a system of linear algebraic equations (SLAE). The construction of the basis is based on a variant of the general solution of a system of elliptic type equations. Each basic state is created by using harmonic polynomials. Orthogonalization is performed by the Gram-Schmidt algorithm. The main mixed task involves dividing the body boundary into two classes. Point movements are set on one part of the border, and surface movements are set on the remaining part. The accuracy of the solution is assessed by two factors: 1) Bessel's inequality; 2) integral quadratic residual of the restored boundary state with BC. The use of both factors leads to the effect of self-sufficiency of the MBS: there is no need to compare trial solutions with reference ones based on other methods. Specific calculations have been performed for two classes of mixed tasks: 1) a bounded bicavous body (a ball with two symmetrically arranged spherical cavities. A rigid counter-displacement of the surfaces of the cavities is set, the outer boundary of the body is free from load; 2) an elastic circular cylinder is fixed with a rigid rod along part of one base. The side surface is loaded with tangential forces directed along the axis of the cylinder. In the second problem, a singular jump in the characteristics of the boundary state along the interface line belonging to the boundary of the body is tested. Numerical-analytical solutions to both problems are constructed, the fields of characteristics of stress-strain states are illustrated, comments on states are made, and conclusions are done.