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