Осесимметричная задача теории упругости рассмотрена в работе Безухова [4].
Расчет методом конечных элементов
Расчет МКЭ выполняется в программных пакетах, например, ANSYS.
Первоначально строится твердотельная трехмерная модель, выполняется ее корректировка по результатам гидродинамического расчета численным методом (методом конечных объемов) в программной пакете МКЭ. Программный пакет выполняет междисциплинарные расчеты, при проектировании насосы осуществляется сквозное проектирование от расчета течения жидкости в проточной части и оптимизации по энергетическому параметру до прочностных расчетов с учетом влияния температуры.
Под расчетом по методу конечных элементов понимается вычислительный процесс на компьютере, состоящий из [5,с.6]:
– описания конечных элементов, численного интегрирования для вычисления элементов матриц,
– объединение матриц отдельных конечных элементов в общую матрицу ансамбля элементов,
– численное решение системы уравнений равновесия.
Решение уравнений равновесия для статических и динамических задач занимает основные затраты машинного времени на вычисления. Инженер-расчетчик может контролировать ход вычисления.
При расчете МКЭ оболочек (т.е. корпусов аппаратов) предполагается [5,с.73], связь конечных элементов в узловых точек (которых конечное число), перемещения узловых точек определяют перемещения конечных элементов (поля конечных элементов). За счет этого используя принцип возможных перемещений можно составить уравнения равновесия для совокупности всех конечных элементов.
Расчет трехмерной задачи
Пример формы трехмерного конечного элемента:
Перемещения тетраэдрического элемента определяется перемещением 12 компонентами перемещений его узлов [6,с.107]:
Компонентами u, v, w определяется вектор перемещений точки.
Матрица деформаций [6,с.108]:
Матрица тепловых деформаций [14,с.109] (θε – средняя температура элемента):
Матрица упругости [14,с.109]:
Матрица напряжений [14,с.109] ({σ0 – аддитивный член}):
Объединяя тетраэдрические элементы, можно разбивать пространство на «кирпичики». В этом случае повышается наглядность разбиения.
Зенкевич указывает [6,с.115] о расчете сосуда высокого давления МКЭ с использованием конечных элементов в виде «кирпичиков». В приводимом примере расчета выполнялся для 10000 степеней свободы. И Зенкевич указывает на то, что при применении более сложных конечных элементов расчет упрощается за счет уменьшения степеней свободы. Но использование сложных элементов не даст преимуществ в сокращении времени подготовки расчета, если процесс разбиения автоматизирован [6,с.169]. В настоящее время в программных пакетах МКЭ используется автоматизированное построение расчетной сетки. При этом при применении сложных элементов сокращается время вычислений, однако ширина матрицы увеличивается и сокращение времени может не происходить. Увеличение размеров конечных элементов приводит к ухудшению аппроксимации конструкции.
Расчет осесимметричной задачи
По данным [6,с.89] , [7,с.229] для решения осесимметричной задачи может быть использован подход плоской задачи. В этом случае треугольный симплекс-элемент вращением образует треугольный тор [7,с.229]. Такой тор показан на рисунке в работе О. Зенкевича [6,с.87]
Объемное тело 3D-модели представляет собой объем, по которому берется интеграл таких треугольных элементов. Отличие осесимметричной задачи от плоской состоит в том, что при деформации оболочки в радиальном направлении вызывает деформацию в окружном направлении. И в рассмотрение должна быть введена четвертая компонента деформации и напряжения по сравнению со случаем плоской задачи [6,с.88]. В плоской задаче компоненты напряжения в направлении, перпендикулярном к координатной плоскости, равны нулю. Трехмерный симплекс-элемент рассматривается аналогично двумерному конечному элементу [7,с.226].
Векторы напряжений и деформаций и матрица упругости по данным [7,с.229]:
Вектор начальной деформации от теплового воздействия [7,с.230]:
Напряжения вычисляются по закону Гука [7,с.233]:
или через узловые перемещения после подстановки