2418

Таблица 5

Низшие частоты собственных колебаний защемленной конической оболочки

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

Оболочка в виде сектора кругового тора с углом раствора 45имеет радиус поперечного сечения 7 10–2 м, радиус образующей 21 10–2 м и толщину 0.1 10–2 м. Материал оболочки имеет следующие характеристики: E = 2.1 1011 Н/м2, = 0.3, = 7.85 103 кг/м3. Одно крайнее поперечное сечение защемлено, другое – свободно. В таблице 10 представлены результаты чис-

ленного анализа.

Таблица 6

Низшие частоты собственных колебаний оболочки в виде сектора кругового тора, защемленной по одному крайнему поперечному сечению

На рис. 50 изображены формы собственных колебаний оболочки.

70

1

3

5

Рис. 52. Низшие собственные формы конической оболочки

71

1

3

5

Рис. 53. Низшие собственные формы сектора торообразной оболочки

72

Разработанный трапециевидный конечный элемент обладает следующими преимуществами:

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

-аналитические выражения для компонентов матриц инертности повышают точность и быстроту их вычисления, по сравнению с численным интегрированием;

-обеспечивает возможность моделирования оболочек с произвольным изменением толщин элементов вдоль меридианов и параллелей срединной поверхности вращения;

-алгоритм генерации конечноэлементной сетки для оболочек со срединными поверхностями вращения легко поддается автоматизации;

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

Моделирование массивных конструкций

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

При построении модели используется конечный элемент в виде пирамиды с треугольным основанием. В качестве узлов i, j, k, m приняты его вершины (рис.54).

Перемещения вдоль осей x,y,z распределяются по линейному закону:

73

Рис. 54. Конечный элемент в форме тетраэдра

иего узловые степени свободы

Ваппроксимирующем полиноме 12 независимых коэффициентов, соответствующих трем степеням свободы в каждом узле.

Для построения матрицы жесткости необходимо аппроксимацию полиномами (1) получить в явном виде, связав коэф-

перемещения u вдоль оси x зависят только от перемещений u1 , u2 , u3 , u4 ; перемещения вдоль оси y зависят только от

v1,v2 ,v3 ,v4 и т.д. (локальные номера узлов 1-4 соответствуют

глобальным номерам узлов i, j, k, m). Если в полиномы (1) подставить координаты узла i, то u = ui и т.д., т.е:

74

Из этого выражения можно найти коэффициенты аппроксимаций:

Аппроксимация перемещений u(x,y,z) вдоль оси x в явном виде примет вид:

Значения a2 ,b2 ,c2 , d2 и т.д. определяются с помощью кру-

говой перестановки узловых индексов. Выражения для v(x,y,z) и w(x,y,z) полностью аналогичны выражению для u(x,y,z).

75

конструкции и назначением для каждого элемента конструкции соответствующих характеристик. На основе этой модели формируются массивы координат узлов, расположенных в вершинах объемных элементов и матрица связей, в которую включаются номера узлов, принадлежащих данному конечному элементу. Задаются узловые нагрузки и граничные условия. В качестве примера моделируется защемленный одним своим крайним поперечным сечением брус и нагруженный двумя силами в плоскости свободного сечения (рис.55).

Рис. 55. Пример расчета

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

76

мент программно представляется в виде совокупности тетраэдральных конечных элемента (рис.56).

Рис. 56. Схема разбиения четырехгранной призмы на тетраэдры

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

77

конечно-элементной модели необходимо стремиться к такой форме конечных элементов, чтобы длины его ребер не отличались друг от друга более, чем в пять раз. Описанный процесс формирования исходных данных производится на предварительном этапе, который выполняет программа brus_d.f. На первом шаге по сформированным матрицам связей и координат узлов вычисляются матрицы жесткости и записываются в отдельный файл для дальнейших расчетов. Этот шаг выполняется программой s1.f . Для выполнения первого шага необходимо задать в файле DT1 количество элементов и узлов, механические характеристики материала элементов. На следующем шаге формируется глобальная матрица жесткости расчетной модели с учетом заданных граничных условий (программа S2.f). Граничные условия, нагрузки и размер формируемой фазы системы уравнений вводятся в файле DT2. Решение полученной системы линейных алгебраических уравнений проводится программой S3.f. Полученный массив узловых перемещений, позволяет затем вычислить напряжения в каждом конечном элементе (программа S4.f). На основании полученной информации можно сделать вывод о напряженно-деформи- рованном состоянии исследуемой конструкции и, при недостаточной прочности, изменить геометрические размеры или механические характеристики материала, проведя описанный расчет вновь. В рассмотренной модели последний узел нагружен двумя силами по 100 Н, узлы 1, 2,3, 4 лишены всех степеней свободы. Материал элементов модели имеет следующие

В приложении 4 представлены тексты программ этапов моделирования массивных конструкций.

Для рассматриваемого примера необходимые данные задаются в файлах:

DT0

44,132,10,60,4,2

0.,0.,0.

0.1,0.,0.

78

0.1,0.1,0.

0.,0.1,0.

0.,0.,0.25

0.1,0.,0.25

0.1,0.1,0.25

0.,0.1,0.25

0.,0.,0.5

0.1,0.,0.5

0.1,0.1,0.5

0.,0.1,0.5

0.,0.,0.75

0.1,0.,0.75

0.1,0.1,0.75

0.,0.1,0.75

0.,0.,1.00

0.1,0.,1.00

0.1,0.1,1.00

0.,0.1,1.00

0.,0.,1.25

0.1,0.,1.25

0.1,0.1,1.25

0.,0.1,1.25

0.,0.,1.50

0.1,0.,1.50

0.1,0.1,1.50

0.,0.1,1.50

0.,0.,1.75

0.1,0.,1.75

0.1,0.1,1.75

0.,0.1,1.75

0.,0.,2.00

0.1,0.,2.00

0.1,0.1,2.00

0.,0.1,2.00

0.,0.,2.25

0.1,0.,2.25

79