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

ANSYS Mechanical

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

vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 2 (примеры из Verification Manual)

Формы колебаний для КЭ-модели (вариант BEAM44)

7 форма

8 форма

9 форма

ЗАО НИЦ СтаДиО (www.stadyo.ru stadyo@stadyo.ru), 2009

I-11

vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 2 (примеры из Verification Manual)

 

 

 

 

 

 

Таблица 1.1

 

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

Тип КЭ

 

№ частоты

Частота, Гц

δ (%)

Примечание

 

 

ANSYS

Источник

 

 

 

 

(NAFEMS)

 

 

 

 

 

 

 

 

Изгиб 1

 

 

 

1,2

42,730

42,650

0,19

 

 

 

3

71,453

71,200

0,35

Кручение 1

 

BEAM188

 

4

124,872

125,000

-0,10

Осевое смещение 1

 

 

5,6

148,544

148,150

0,27

Изгиб 2

 

 

 

 

 

 

7

212,598

213,610

-0,47

Кручение 2

 

 

 

8,9

282,239

283,470

-0,43

Изгиб 3

 

 

 

1,2

42,634

42,650

-0,04

Изгиб 1

 

 

 

3

71,713

71,200

0,72

Кручение 1

 

BEAM189

 

4

125,000

125,000

0,00

Осевое смещение 1

 

 

5,6

148,195

148,150

0,03

Изгиб 2

 

 

 

 

 

 

7

215,080

213,610

0,69

Кручение 2

 

 

 

8,9

284,095

283,470

0,22

Изгиб 3

 

 

 

1,2

42,639

42,650

-0,03

Изгиб 1

 

 

 

3

71,600

71,200

0,56

Кручение 1

 

BEAM44

 

4

125,129

125,000

0,10

Осевое смещение 1

 

 

5,6

148,566

148,150

0,28

Изгиб 2

 

 

 

 

 

 

7

216,569

213,610

1,39

Кручение 2

 

 

 

8,9

287,175

283,470

1,31

Изгиб 3

 

Максимальная по абсолютной величине погрешность δ (расхождение с [Источник]):

для BEAM188 δ = 0,47%

для BEAM189 δ = 0,72%

для BEAM44 δ = 1,39%

ЗАО НИЦ СтаДиО (www.stadyo.ru stadyo@stadyo.ru), 2009

I-12

vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 2 (примеры из Verification Manual)

Пример 2 (VMP09-T52). Определение собственных частот и форм колебаний шарнирно опертой квадратной плиты

Источник

Reference NAFEMS FVB Test 52

Тип задачи:

Определение собственных частот и форм колебаний

 

(частичная проблема собственных значений)

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

SHELL181 (3-D 4-узловые оболочечные элементы)

 

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

 

SOLID186 (3-D 20-узловые объемные элементы)

 

SOLID187 (3-D 10-узловые объемные элементы)

 

SHELL281 (3-D 8-узловые оболочечные элементы)

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

vmp09-t52-181.mac

 

vmp09-t52-185.mac

 

vmp09-t52-186.mac

 

vmp09-t52-187.mac

 

vmp09-t52-281.mac

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

Рассматривается квадратная шарнирно опертая плита “средней” толщины. Требуется определить первые 10 собственных частот и форм колебаний с использованием различных типов оболочечных и объемных КЭ.

Рис. 2.1 Расчетная схема плиты

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

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

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

Длина стороны плиты L = 10 м Толщина W = 1 м

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

Квадратная плита расположена в плоскости X-Y.Локальные координатные оси элементов сонаправлены глобальным осям.

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

SHELL181 – пространственно-оболочечный конечный элемент, имеющий четыре узла. В данном элементе реализована теория оболочек Миндлина-Рейсснера.

SOLID185 – объемный (трехмерный) гексаэдрический восьмиузловой конечный элемент.

ЗАО НИЦ СтаДиО (www.stadyo.ru stadyo@stadyo.ru), 2009

II-1

vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 2 (примеры из Verification Manual)

SOLID186 – объемный (трехмерный) гексаэдрический квадратичный элемент, имеющий двадцать узлов. Элемент имеет квадратичное представление перемещений и в состоянии использовать нерегулярную форму сетки.

SOLID187 – объемный (трехмерный) тетраэдрический квадратичный элемент. Определяется десятью узлами.

SHELL281 – пространственно-оболочечный восьмиузловой элемент. В данном элементе реализована теория оболочек Миндлина-Рейсснера.

При расчёте с применением элементов SHELL181 и SHELL281 была включена опция полного интегрирования с несовместными формами.

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

 

 

Число

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

Узлы×Элементы

Тип КЭ

степеней

(длина×ширина×высота), м

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

 

 

свободы

 

 

 

 

1

SHELL181

2646

0,5×0,5

441×400

2

SOLID185

3993

1,0×1,0×0,1

1331×1000

3

SOLID186

465

2,5×2,5×1,0

155×16

4

SOLID187

1047

2,5 (размер грани)

349×152

5

SHELL281

7686

0,5×0,5

1281×400

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

При использовании элементов SHELL181, SHELL281 и SOLID185 запрещены вертикальные перемещения узлов, находящихся на боковых гранях плиты.

x = 0 м; x = 10 м; y = 0 м; y = 10 м.

Uz = 0

x = 0 м; x = 10 м; y = 0 м; y = 10 м. z = 0-1 м.

Uz = 0

При использовании элементов SOLID186 и SOLID187 запрещены вертикальные перемещения узлов находящихся на нижних ребрах боковых граней плиты.

x =0 м; x = 10 м; y = 0 м; y = 10 м. z = 0м. Uz = 0

ЗАО НИЦ СтаДиО (www.stadyo.ru stadyo@stadyo.ru), 2009

II-2

vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 2 (примеры из Verification Manual)

Рис. 2.2 Изометрия КЭ-модели плиты с указанием закреплений (элементы SHELL181). 3D-визуализация

ЗАО НИЦ СтаДиО (www.stadyo.ru stadyo@stadyo.ru), 2009

II-3

vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 2 (примеры из Verification Manual)

Рис. 2.3 Изометрия КЭ-модели плиты с указанием закреплений (элементы SOLID185). 3D-визуализация

Рис. 2.4 Изометрия КЭ-модели плиты с указанием закреплений (элементы SOLID186). 3D-визуализация

ЗАО НИЦ СтаДиО (www.stadyo.ru stadyo@stadyo.ru), 2009

II-4

vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 2 (примеры из Verification Manual)

Рис. 2.5 Изометрия КЭ-модели плиты с указанием закреплений (элементы SOLID187). 3D-визуализация

Рис. 2.6 Изометрия КЭ-модели плиты с указанием закреплений (элементы SHELL281). 3D-визуализация

ЗАО НИЦ СтаДиО (www.stadyo.ru stadyo@stadyo.ru), 2009

II-5

vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 2 (примеры из Verification Manual)

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

Результатом расчёта являются первые (низшие) 10 собственных частот и форм колебаний пластины, из которых первые 3 нулевые соответствуют движению как жесткого тела в плоскости пластины, а 7 следующих – колебаниям как упругого тела. Ниже в таблице 2.1 приведено сравнение результатов по ANSYS и данных NAFEMS.

Собственные формы упругих колебаний для КЭ-модели (элемент SOLID185)

4 форма

5 форма

7 форма

8 форма

9 форма

10 форма

ЗАО НИЦ СтаДиО (www.stadyo.ru stadyo@stadyo.ru), 2009

II-6

vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 2 (примеры из Verification Manual)

Собственные формы колебаний для КЭ-модели (элемент SHELL281)

4 форма

5 форма

7 форма

8 форма

9 форма

10 форма

ЗАО НИЦ СтаДиО (www.stadyo.ru stadyo@stadyo.ru), 2009

II-7

vk.com/club152685050ANSYS Mechanical. |Верификационныйvk.com/id446425943отчет. Том 2 (примеры из Verification Manual)

 

 

 

 

 

 

Таблица 2.1

 

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

Тип КЭ

 

№ частоты

Частота, Гц

δ (%)

Примечание

 

 

ANSYS

NAFEMS

 

 

 

 

 

 

 

 

 

4

44,396

45,897

3,271

 

 

SHELL181

 

5,6

108,221

109,440

1,114

 

 

 

7

164,503

167,890

2,017

 

 

 

 

8

194,097

193,590

0,262

 

 

 

 

9,10

206,734

206,190

0,264

 

 

 

 

4

48,925

45,897

6,597

 

 

SOLID185

 

5,6

123,931

109,440

13,241

 

 

 

7

184,648

167,890

9,981

 

 

 

 

8

195,655

193,590

1,066

 

 

 

 

9,10

209,403

206,190

1,558

 

 

 

 

4

43,837

45,897

4,488

 

 

 

 

5,6

106,751

109,440

2,457

 

 

SOLID186

 

7

161,176

167,890

3,999

 

 

 

8

193,796

193,590

0,107

 

 

 

 

 

 

 

 

9

205,226

206,190

0,467

 

 

 

 

10

205,971

206,190

0,106

 

 

 

 

4

45,387

45,897

1,111

 

 

 

 

5

112,311

109,440

2,623

 

 

SOLID187

 

6

112,386

109,440

2,692

 

 

 

7

179,239

167,890

6,760

 

 

 

 

8

193,940

193,590

0,181

 

 

 

 

9

206,600

206,190

0,199

 

 

 

 

10

206,610

206,190

0,204

 

 

 

 

4

44,110

45,897

3,894

 

 

 

 

5,6

106,982

109,440

2,246

 

 

SHELL281

 

7

162,337

167,890

3,308

 

 

 

8

193,609

193,590

0,010

 

 

 

 

 

 

 

 

9

202,216

206,190

1,927

 

 

 

 

10

202,570

206,190

1,756

 

 

Максимальная по абсолютной величине погрешность δ:

для SHELL181 δ = 3,270%

для SOLID185 δ = 13,241%

для SOLID186 δ = 4,488%

для SOLID187 δ = 6,760%

для SHELL281 δ = 3,894%

ЗАО НИЦ СтаДиО (www.stadyo.ru stadyo@stadyo.ru), 2009

II-8