Изложен алгоритм построения градиентной теории термоупругости первого порядка (strain-gradient theory). Модель является обобщением прикладной модели межфазного слоя, и получена путем привлечения гипотезы Дюамеля-Неймана. Приводится математическая постановка модели, которая сводится к краевой задаче с разрешающим оператором повышенного порядка, как и для всех моделей градиентной теории упругости. Физические соотношения модели учитывают более высокие градиенты температурного поля, чем в классической теории. В работе приведены тестовые расчеты для слоистых композитных структур, на основе которых продемонстрированы особенности модели, связанные с возможностью учета локальной концентрации температурных напряжений в области границ сред.