ANSYS Mechanical
.pdfvk.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 |