В работе показано, что численно-аналитический метод граничных элементов является эффективным инструментом решения трехмерных задач теории упругости для тел с включениями, которые характеризуются как сложной геометрической формой, так и условиями скользящего контакта с окружающей средой. В частности, он позволяет непосредственное включение таких условий в гранично-интегральную формулировку задач, приводя к уменьшению количества решаемых уравнений. На примере гладкого контакта конечного цилиндрического волокна с матрицей при всестороннем сжатии установлено качественные и количественные различия в распределениях параметров напряженно-деформированного состояния структуры сравнительно со случаем идеального контакта. Этот фактор следует учитывать при моделировании эффективных механических свойств гранулированных и волокнистых композитов, допускающих нарушения сцепления наполнителя с матрицей.