Перейти до змісту

Приклади реальних моделей для пружного статичного аналізу

Моделі аналізу

У таблиці 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

Шатун (EX01A)

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

Блок з отвором (EX02)

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

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

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

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

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

Келих (EX05, EX10A)

Рис. 9.2.5. Келих (EX05, EX10A)

Результати аналізу

Приклади результатів аналізу

Приклади результатів аналізу показано на рис. 9.2.6–9.2.9.

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

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

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

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

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

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

Результат аналізу EX04 (деформована форма, масштаб ×100)

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

Результати перевірки точності та продуктивності аналізу на прикладі EX02

Аналіз було виконано за допомогою універсального комерційного програмного забезпечення на моделі, еквівалентній моделі блока з отвором із перевірочного прикладу EX02. На рис. 9.2.10 показано порівняння максимальних і мінімальних значень компонент напруження з результатами FrontISTR. Із рисунка видно, що компоненти напруження дуже добре узгоджуються.

Далі наведено результати дослідження впливу декомпозиції області на розподіл напружень. Декомпозицію області виконано методом RCB: по два поділи вздовж кожної з осей X, Y і Z, загалом на 8 областей. Стан розбиття показано на рис. 9.2.11. На рис. 9.2.12 показано розподіли напружень за результатами аналізу для однієї області та для поділу на 8 областей.

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

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

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

Рис. 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

Порівняння історії збіжності розв’язувача HEC-MW (порогове значення критерію збіжності: 1.0×10⁻⁸)

Рис. 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