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