2418

ЭВМ, не допускать неопределенных и нестандартных ситуаций, не требовать больших затрат машинного времени при реализации, обеспечивать алгоритмичность реализации задач механики на ЭВМ.

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

В первой схеме для конструирования дискретной модели используется семейство четырехугольных зон (квадратичные четырехугольники). Программа разработана для моделирования двумерных областей, которые могут быть представлены совокупностью прямоугольников и треугольников, границы которых описываются кривыми второго порядка. При описании треугольной зоны, две стороны квадратичного четырехугольника используются для задания одной стороны треугольной области.

Нумерация узлов в квадратичных четырехугольниках производится в соответствии с рис. 16. Безразмерные координаты узла 1 всегда равны = =-1.

Для каждого четырехугольника выполняются следующие процедуры.

1.Задается число строк и столбцов узлов в соответствии с требуемой степенью измельчения.

2.Производится проверка граничных узлов на предмет отсутствия присвоения им номера. Если такие узлы есть, то за ними сохраняются номера, присвоенные ранее.

20

Рис. 16

Узлы нумеруются последовательно, начиная от точки с координатами

= 1, =+1 и двигаясь слева направо (при изменении от –1 до +1) и сверху вниз (при изменении от +1 до –1).

Узлы, пронумерованные ранее, пропускаются. Номера граничных узлов сохраняются для последующего учета при формировании соседних зон.

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

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

Вкачестве примера рассмотрена плоская прямоугольная область с круговым вырезом (рис. 17). Для обеспечения достаточно полного описания границы кругового выреза выделено 12 четырехугольных зон (на рис. 17 обозначены римскими цифрами). Стороны четырехугольных зон, примыкающие к вырезу, разделены на два равных отрезка: это обеспечивается

21

заданием соответствующего числа строк и столбцов узлов в каждой зоне.

В результате такого расположения узлов вдоль периметра выреза получена равномерная сетка треугольных элементов. Исходные данные для рассматриваемой области размещены в файле «DT»:

e

22

63 13 57 41 11 12 13 40 55 56

73 13 55 40 13 14 15 39 53 54

85 13 53 39 15 16 17 18 19 38

95 3 51 52 53 38 19 20 21 37

10 5 3 49 50 51 37 21 22 23 36

11 5 7 27 35 49 36 23 24 25 26

12 3 7 29 34 47 48 49 35 27 28

Впервой строке обозначена решаемая проблема, в следующей указано число четырехугольных зон и количество граничных точек. Следующие 16 строк заполнены координатами x и y граничных точек. За координатами введены 12 строк, соответствующих числу четырехугольных зон, в которых указаны номера смежных зон. Данные каждой зоны считываются в следующем порядке: номер зоны, число строк и число столбцов, глобальные номера узлов, определяющие четырехугольник.

Вычисленные значения координат узлов и массив связей треугольных элементов сохраняются для проведения дискретного моделирования методом конечных элементов. Сетка конечных элементов для рассматриваемого примера представле-

23

на на рис. 18. Число узлов построенной сетки равно 290, число треугольных конечных элементов – 496.

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

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

Рис. 18

Тип элемента вводится перед началом вычислений: «1» – для треугольного элемента с узлами в вершинах, «2» – для треугольного элемента с 6 узлами; если ввести тип «0», то происходит переход к подпрограмме вывода полученных результантов. Затем задается число участков, на которые требуется разделить каждую сторону треугольной зоны. После этого осуществляется ввод номеров узлов и их координат, при этом номера угловых узлов надо подобрать таким образом, чтобы они соответствовали введенному числу деления сторон на элементы.

24

Если ввести отрицательный номер узла, то последующий ввод данных прекратится и произойдет переход к подпрограмме генерации сетки. После ввода номера последнего узла и его координат необходимо ввести число «0», при этом после генерации сетки осуществляется вывод результатов вычислений. При введении новых узлов используется общая система нумерации для всех треугольных зон. Вычисления координат узлов, расположенных вдоль сторон, осуществляется при помощи системы L – координат /3/.

Для угловых узлов треугольных зон вводятся функции формы:

N1 (2L1-1)L1; N2 (2L2-1)L2; N3 (2L3-1)L3 .

Для узлов, расположенных на серединах сторон, функции формы задаются в виде:

N4=4L1L2 ; N5=4L2L3 ; N6=4L3L1 .

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

X NiXi ; Y NiYi ; Z NiZi ( i=1,6) ,

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

В качестве примера, на рис. 19 изображена плоская область с нанесенной исходной сеткой триангуляции; сетка, измельченная с использованием описанного алгоритма, показана на рис. 20.

Третья схема построения сетки конечных элементов для тонкостенных конструкций производится на основе массивов координат узловых точек, полученных при помощи генератора узлов поверхности [6].

25

Использована локальная система координат , , связанная с исследуемой областью (рис. 21).

z

2

3

4

x

Рис. 21 Локальная и глобальная системы координат

характер изменения глобальных координат вдоль границы. Преобразование локальных координат в глобальную сис-

тему может быть определено, например, следующим образом

[6]:

функции для стороны n.

Для Y и Z могут быть составлены аналогичные соотношения. С помощью этих преобразований могут быть вычислены глобальные координаты внутренних узлов, если известны их

27

локальные координаты и граничные функции. Локальная ко-

Глобальные координаты, вычисленные в соответствии с описанным алгоритмом, используются для автоматизированного формирования сетки конечных элементов.

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

Сетка конечных элементов рассматриваемой оболочки представлена на рис. 22.

Рис. 22 Сетка конечных элементов

Постановка задачи предусматривает исследование собственных колебаний линейно-упругой тонкостенной оболочки, защемленной вдоль одной из сторон. Для решения неполной проблемы собственных значений использован метод итераций в подпространстве [7].

28

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

Некоторые из низших частот и форм собственных колебаний лопатки представлены на рис. 23–28.

Рис. 23 Низшая частота 1 = 526 с-1 собственных колебаний лопатки

Рис. 24 Низшая частота 2 = 1469 с-1 собственных колебаний лопатки

29