Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

ANSYS Mechanical

.pdf
Скачиваний:
901
Добавлен:
31.08.2019
Размер:
61.81 Mб
Скачать

vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 3 (“исследовательские” задачи)

Решение задачи S11 в объёмной постановке, с применением КЭ SOLID45 (закрепление по нижней грани)

Рис. 3.12 Изометрия КЭ-модели балки (SOLID45) с указанием закреплений и нагрузок. 3D-визуализация

x = 0 м

x = 2 м

Рис. 3.13 Закрепления на торцах балки (SOLID45). 3D-визуализация

ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009

III-10

vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 3 (“исследовательские” задачи)

Рис 3.14 1-ая форма потери устойчивости для SOLID45-модели 2. Критический погонный момент mизг = 657,227 т·м/м

Рис 3.15 2-ая форма потери устойчивости для SOLID45-модели 2. Критический погонный момент mизг = 1344,000 т·м/м

ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009

III-11

vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 3 (“исследовательские” задачи)

Задача 4. Устойчивость П-образной рамы с абсолютно жёсткими стойками

Источник

1. А. В. Перельмутер, В. И. Сливкер “Устойчивость

 

равновесия конструкций и родственные проблемы”. (М.,

 

“СОФТ”, 2007, стр. 551-554)

Тип задачи:

Определение критических нагрузок и форм потери

 

устойчивости, задача на собственные значения (ANTYPE = 1)

Тип верифицируемых КЭ:

BEAM44 (3D 2-узловой балочный элемент)

 

BEAM188 (3D 2-узловой балочный элемент “Тимошенко)

 

SOLID45 (3-D 8-узловой объемный элемент)

Входные файлы:

sliv_12.mac

Постановка задачи

Рассматривается П-образная рама с абсолютно жёсткими стойками и ригелем конечной жёсткости, шарнирно опёртая, нагруженная сосредоточенными силами (см. рис.). Необходимо определить критические силы и формы потери устойчивости.

EI =

P

l

 

EI

8

EI =

P

8

h

Физические характеристики

Модуль упругости E = 2 107 т/м2 Коэффициент Пуассона ν = 0,3

Геометрические характеристики рамы (см. рисунок) l = 1 м

h = 1 м

Сечения горизонтального стержня – квадрат со стороной 0,1 м.

Описание КЭ-модели

Задача решалась 2-мя способами: в стержневой и объёмной постановках.

Для моделирования абсолютно жёстких стоек в случае стержневой модели использовались стержни по геометрическим характеристикам сечения идентичные таковым для горизонтального стержня. В случае решения задачи в объёмной постановке жёсткие стержни были приняты сечением 0,1×0,02 м.

Модуль упругости жёстких стержней превышал таковой для деформируемого стержня на 3 – 5 порядков.

Характерные размеры элементов, вычислительная размерность задачи (число степеней свободы) и количество узлов и элементов отображены в следующей таблице:

ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009

IV-1

vk.com/club152685050 | vk.com/id446425943

 

” задачи)

 

 

ANSYS Mechanical. Верификационный отчет. Том 3 (“исследовательские

 

 

 

 

 

 

 

 

 

Тип КЭ

Число степеней свободы

Характерные размеры

Узлы×КЭ

 

 

 

 

элементов, м

(количество)

 

 

BEAM44

72

 

0,2

12×11

 

 

BEAM188

72

 

0,2

12×11

 

4

SOLID45

825

Вертикальные стержни

315×136

 

 

 

0,2×0,01×0,05

 

 

 

 

 

 

 

 

 

 

Горизонтальные стержни

 

 

 

 

 

0,05×0,05×0,05

 

Для решения данной задачи применялись 2 типа КЭ:

BEAM188 – пространственный линейный элемент балки Тимошенко, имеющий 3 узла: 2 “содержательных” узла располагаются на оси элемента (по краям), третий является узлом ориентации и может быть общим для нескольких элементов;

SOLID45 – объемный элемент для моделирования трехмерного НДС. Определяется восемью узлами, каждый из которых имеет по три степени свободы и свойствами ортотропного материала.

Граничные условия

Нижние точки стоек

Ux = 0 Uy = 0

Нагрузки

Вертикальные сосредоточенные силы (см рис. …)

P = 1 т

Теоретическое решение

Моменты инерции стержня

I y = Iz =0,8333 105 м4

“Теоретическая” критическая сила вычисляются по формуле [Источник 1] (см. рис.

…)

 

 

 

 

Pcr =

6EI

=

6 2 107 0,8333 105

=1000т

lh

1 1

 

 

 

Методика расчёта

Определение собственных чисел (критических моментов) и форм потери устойчивости проводилось блочным методом Ланцоша. Вычислены 10 низших собственных чисел и соответствующих им форм потери устойчивости.

Результаты расчёта

Верифицируемым результатом расчёта являются низшие критические моменты и формы потери устойчивости. Ниже приведено сравнение результатов по ANSYS и данных [Источник 1].

 

 

Сопоставление результатов расчёта

Таблица 4.1

 

 

 

 

 

Критическая сила, Pcr , т

 

 

 

ANSYS

 

 

 

 

задачи

формы

 

 

Теория [1]

δ (%)

 

BEAM44

BEAM188

SOLID45

 

4

1

968,134

999,415

1071,000

1000

0,058

 

2

400333,000

400333,000

137806,000

-

 

 

 

 

 

ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009

IV-2

vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 3 (“исследовательские” задачи)

Решение задачи S12 в стержневой постановке, с применением элементов BEAM44

Рис. 4.1 КЭ-модель рамы с указанием закреплений, нагрузок и нумерацией узлов и элементов (элементы BEAM44)

Рис. 4.2 КЭ-модель рамы (элементы BEAM44). 3D-визуализация

ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009

IV-3

vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 3 (“исследовательские” задачи)

Рис. 4.3 1-ая форма потери устойчивости для BEAM44-модели. 3D-визуализация.

Критическая сила Pкр = 968,137 т

Рис. 4.4 1-ая форма потери устойчивости для BEAM44-модели.

Критическая сила Pкр = 968,137 т

ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009

IV-4

vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 3 (“исследовательские” задачи)

Рис. 4.5 2-ая форма потери устойчивости для BEAM44-модели. 3D-визуализация.

Критическая сила Pкр = 400333 т

Рис. 4.6 2-ая форма потери устойчивости для BEAM44-модели.

Критическая сила Pкр = 400333 т

ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009

IV-5

vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 3 (“исследовательские” задачи)

Решение задачи S12 в стержневой постановке, с применением элементов BEAM188

Рис. 4.7 КЭ-модель рамы с указанием закреплений, нагрузок и нумерацией узлов и элементов (элементы BEAM188)

Рис. 4.8 КЭ-модель рамы (элементы BEAM188). 3D-визуализация

ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009

IV-6

vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 3 (“исследовательские” задачи)

Рис. 4.9 1-ая форма потери устойчивости для BEAM188-модели. 3D-визуализация.

Критическая сила Pкр = 999,415 т

Рис. 4.10 1-ая форма потери устойчивости для BEAM188-модели.

Критическая сила Pкр = 999,415 т

ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009

IV-7

vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 3 (“исследовательские” задачи)

Рис. 4.11 2-ая форма потери устойчивости для BEAM188-модели. 3D-визуализация.

Критическая сила Pкр = 400333 т

Рис. 4.12 2-ая форма потери устойчивости для BEAM188-модели.

Критическая сила Pкр = 400333 т

ЗАО НИЦ СтаДиО, (www.stadyo.ru stadyo@stadyo.ru), НОЦ КМ МГСУ (niccm@mgsu.ru), 2009

IV-8