Компьютерное моделирование плоских задач термоупругости: сравнительный анализ решений в связанной и несвязанной постановках

Проведено компьютерное моделирование плоской линейной задачи термоупругости в связанной и несвязанной постановках для неоднородной среды с отверстием. В случае несвязанной постановки сначала решалась независимая часть задачи термоупругости: уравнение нестационарной теплопроводности без учета члена, отвечающего за механическую мощность внутренних сил, при граничных условиях Дирихле. Находилось температурное поле, при котором из уравнений равновесия и соотношения Дюамеля-Неймана определялись статические термоупругие напряжения при граничных условиях в напряжениях. Дискретизация дифференциальных уравнений осуществлялась в соответствии с методом конечных элементов, базирующимся на построении векторно-матричной системы уравнений на основе слабой формы уравнений термомеханики при условии квазистатического деформирования. Метод конечных элементов реализован в коде PIONER, разработанном в ИГиЛ СО РАН. В качестве тестовой задачи рассматривалась задача охлаждения и деформирования полого цилиндра с заданными температурами и напряжениями на внутренней и внешней поверхностях. При решении задачи были использованы 8-узловые конечные элементы с неполной биквадратичной аппроксимацией геометрии и перемещений, позволяющие моделировать плоскую деформацию. Численные эксперименты показали, что для данного класса задач решения уравнений термоупругости в связанной и несвязанной постановках при определенных ограничениях, упрощающих задачу, а именно в отсутствие массовых сил, начальных напряжений, подвода тепла и конвекции на части поверхности, дают близкие по значениям температуры, напряжения, перемещения. Тем самым подтверждается сформировавшееся в научной среде мнение о том, что в линейной термоупругости член в уравнении теплопроводности, отвечающий за механическую мощность внутренних сил, не оказывает существенного влияния на результат решения термоупругой системы.