Приклади реальних моделей для пружного статичного аналізу¶
Моделі аналізу¶
У таблиці 9.2.1 наведено перелік перевірочних прикладів реальних моделей для пружного статичного аналізу. Геометрію моделей (за деякими винятками) показано на рис. 9.2.1–9.2.5. Для виконання прикладів з типами елементів 731 і 741 окремо потрібен розв’язувач прямого методу.
Таблиця 9.2.1. Перевірочні приклади реальних моделей для пружного статичного аналізу
| Назва випадку | Тип елемента | Перевірочна модель | Кількість вузлів | Кількість ступенів свободи |
|---|---|---|---|---|
| EX01A | 342 | Шатун (100 тис. вузлів) | 94,074 | 282,222 |
| EX01B | 342 | Шатун (330 тис. вузлів) | 331,142 | 993,426 |
| EX02 | 361 | Блок з отвором | 37,386 | 112,158 |
| EX03 | 342 | Лопатка турбіни | 10,095 | 30,285 |
| EX04 | 741 | Циліндрична оболонка | 10,100 | 60,600 |
| EX05A | 731 | Келих (груба сітка) | 7,240 | 43,440 |
| EX05B | 731 | Келих (середня сітка) | 48,803 | 292,818 |
| EX05C | 731 | Келих (дрібна сітка) | 100,602 | 603,612 |

Рис. 9.2.1. Шатун (EX01A)

Рис. 9.2.2. Блок з отвором (EX02)

Рис. 9.2.3. Лопатка турбіни (EX03, EX06)

Рис. 9.2.4. Циліндрична оболонка (EX04, EX09)

Рис. 9.2.5. Келих (EX05, EX10A)
Результати аналізу¶
Приклади результатів аналізу¶
Приклади результатів аналізу показано на рис. 9.2.6–9.2.9.

Рис. 9.2.6. Результат аналізу EX01A (напруження Мізеса та деформована форма, масштаб ×10)

Рис. 9.2.7. Результат аналізу EX02 (напруження Мізеса та деформована форма, масштаб ×100)

Рис. 9.2.8. Результат аналізу EX03 (деформована форма, масштаб ×10)

Рис. 9.2.9. Результат аналізу EX04 (деформована форма, масштаб ×100)
Результати перевірки точності та продуктивності аналізу на прикладі EX02¶
Аналіз було виконано за допомогою універсального комерційного програмного забезпечення на моделі, еквівалентній моделі блока з отвором із перевірочного прикладу EX02. На рис. 9.2.10 показано порівняння максимальних і мінімальних значень компонент напруження з результатами FrontISTR. Із рисунка видно, що компоненти напруження дуже добре узгоджуються.
Далі наведено результати дослідження впливу декомпозиції області на розподіл напружень. Декомпозицію області виконано методом RCB: по два поділи вздовж кожної з осей X, Y і Z, загалом на 8 областей. Стан розбиття показано на рис. 9.2.11. На рис. 9.2.12 показано розподіли напружень за результатами аналізу для однієї області та для поділу на 8 областей.

Рис. 9.2.10. Порівняння компонент напруження EX02 з універсальним програмним забезпеченням

Рис. 9.2.11. Результат поділу EX02 на 8 областей методом RCB

Рис. 9.2.12. Відмінності розподілу напруження Мізеса внаслідок декомпозиції області
З рис. 9.2.12 видно, що відмінностей між двома результатами немає і вони повністю збігаються.
Далі в таблиці 9.2.2 наведено порівняння часу виконання для різних налаштувань використаного розв’язувача HEC-MW. На рис. 9.2.13 також показано історію збіжності до отримання розв’язку.
Таблиця 9.2.2. Порівняння часу виконання розв’язувача HEC-MW
| Розв’язувач | Час виконання (s) |
|---|---|
| CGI | 38.79 |
| CGscale | 52.75 |
| BCGS | 60.79 |
| CG8 | 6.65 |

Рис. 9.2.13. Порівняння історії збіжності розв’язувача HEC-MW (порогове значення критерію збіжності: \(1.0 \times 10^{-8}\))
Порівняння часу обчислення за перевірочним прикладом EX01A¶
На перевірочному прикладі EX01A (шатун) досліджено прискорення обчислень завдяки декомпозиції області. Для обчислень використано 24-вузловий кластер із Xeon 2.8GHz. Результати показано на рис. 9.2.14. Із рисунка видно, що швидкість обчислень зростає пропорційно кількості областей.
Також досліджено відмінності часу обчислення залежно від обчислювального середовища. Результати наведено в таблиці 9.2.3.

Рис. 9.2.14. Ефект прискорення завдяки декомпозиції області
Таблиця 9.2.3. Порівняння часу обчислення на різних комп’ютерах (1 CPU)
| CPU | Частота (GHz) | OS | Час CPU (s) | Час розв’язувача (s) |
|---|---|---|---|---|
| Xeon | 2.8 | Linux | 850 | 817 |
| Pentium III | 0.866 | Win2000 | 2008 | 1980 |
| Pentium M | 0.760 | WinXP | 1096 | 1070 |
| Pentium 4 | 2.0 | WinXP | 802 | 785 |
| Pentium 4 | 2.8 | WinXP | 738 | 718 |
| Celeron | 0.700 | Win2000 | 2252 | 2215 |
| Pentium 4 | 2.4 | WinXP | 830 | 804 |