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