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

ANSYS Mechanical

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

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

Рис. 2.45 Задача 2e. 2-ая форма потери устойчивости для SOLID45-модели. Критический момент Mизг = 9864 т·м

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

II-29

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

Задача 3 Устойчивость шарнирно опёртого стержня, загруженного равномерно распределённым изгибающим моментом

Источник

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

 

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

 

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

Тип задачи:

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

 

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

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

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

 

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

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

sliv_11.mac

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

Тестовый пример S12 – устойчивость изгибаемого прямого стержня (шарнирно опёртого длиной l по оси X, поперечное сечение с главными моментами инерции Iy и Iz и моментом инерции при кручении Ix), загруженного равномерно распределённым изгибающим моментом mизг. Эти примеры приведены в книге [Источник 1] на стр. 582 рис. 8.17. Там же даны теоретические результаты для критической нагрузки.

Требуется определить значение критического изгибающего момента mизг и формы потери устойчивости.

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

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

Геометрические характеристики

Длина стержня l = 2 м

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

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

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

В случае решения в стержневой постановке распределённый по длине изгибающий момент был заменён сосредоточенными моментами, приложенными в узлы. Сосредоточенные моменты задавались посредством создания жёстких рычагов длиной 0,01 м и приложения к их концам пары сил (см. рис…).

При решении задачи в объёмной постановке было рассмотрено 2 варианта закрепления балки (рис. 3.9, 3.13): на уровне половины высоты торцевого сечения (1-й случай) и по нижней грани торцевого сечения (2-й случай).

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

III-1

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

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

Тип КЭ

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

Размеры КЭ, м

Узлы×КЭ

1

BEAM188

126

0,1

21×20

2

SOLID45

3087

0,1×0,0167×0,0167

1029×960

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

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

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

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

Объёмная постановка

Стержневая постановка (см. рис…)

Край стержня x = 0 м

1-й случай см. рис.3.9

Ux = 0

2-й случай см. рис. 3.13

Uy = 0

 

Uz = 0

 

Rotx = 0

 

Край стержня x = 2 м

 

Uy = 0

 

Uz = 0

 

Rotx = 0

 

Нагрузки

При решении задачи в стержневой постановке к концам «жёстких» рычагов прикладывается система горизонтальных сил.

При решении задачи в объёмной постановке момент задан как пара равномерно распределённых по верхней и нижней граням стержня касательных давлений.

В обоих случаях приложенная нагрузка эквивалентна распределённому по длине изгибающему (в плоскости стержня) моменту mизг = 1 т·м/м.

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

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

I y = Iz = 0,8333 105 м4

Момент инерции при кручении

I x =1,4263 105 м4

Модуль сдвига

 

E

 

2 107

2

G =

 

=

 

= 7692307т/ м

 

2(1 +ν)

2 (1 0,3)

 

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

mиз = ± 22l2,3 EIzGIx = 2222,3 2 107 0.8333 105 7692307 1.4263 105 = 755,567 тм/ м

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

III-2

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

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

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

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

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

Таблица 3.1

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

Критический погонный изгибающий момент, mизг , т м / м

 

ANSYS

 

 

 

 

задачи

формы

 

 

 

Теория [1]

δ (%)

BEAM188

SOLID45

SOLID45*

 

3

1

756,727

723,615

657,227

 

755,567

0,154

2

1322,000

1350,000

1344,000

 

-

 

 

 

* – Закрепление по нижней грани балки (2 случай)

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

III-3

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

Решение задачи S11 в стержневой постановке, с применением КЭ BEAM188

Рис. 3.1 Изометрия BEAM188-модели балки с указанием закреплений и нагрузок

Рис. 3.2 Визуализация нагрузок (КЭ BEAM188)

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

III-4

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

Рис. 3.3 Изометрия BEAM188-модели балки. 3D-визуализация

Рис.3.4 1-ая форма потери устойчивости для BEAM188-модели. 3D-визуализация. Критический погонный момент mизг = 756,727 т·м/м

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

III-5

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

Рис. 3.5 1-ая форма потери устойчивости для BEAM188-модели. Критический момент mизг = 756,727 т·м/м

Рис.3.6 2-ая форма потери устойчивости для BEAM188-модели. 3D-визуализация. Критический погонный момент mизг = 1322,000 т·м/м

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

III-6

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

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

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

III-7

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

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

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

x = 0 м

x = 2 м

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

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

III-8

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

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

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

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

III-9