Метод Обмеженого Елементу Пояснений: Від Мережі до Силових Напружень Моста

Від рівнянь, що керують, до алгебраїчної системи

Навантажений міст регулюється рівняннями лінійної еластичності: рівновага (∇·σ = 0), закон Гука, який пов’язує напруження зі деформацією (σ = C:ε), та відношення деформації до зміщення (ε = ½(∇u + (∇u)ᵀ)). Аналітичні розв'язки існують лише для найпростіших геометрій, тому FEM вирішує еквівалентну слабку (варіаційну) форму — інтегруючи умову рівноваги відносно віртуального зміщення замість розв’язання диференціального рівняння в точці. Ця переформулювання дозволяє комп'ютеру обробляти структури будь-якої складності, від окремого балону зварного мосту до повної плити моста.

Функції форми: наближення невідомого поля

У кожному елементі невідоме поле зсуву наближається як зважена сума функцій форми Nᵢ(x), де ваги – це нодальні значення:

u(x) ≈ Σᵢ Nᵢ(x) · uᵢ (в межах одного елемента) Основні властивості будь-якого допустимого набору функцій форми: Nᵢ(xⱼ) = δᵢⱼ (1 у своїй ноді, 0 у всіх інших нодах) Σᵢ Nᵢ = 1 (розділ unity – точне моделювання деформацій тіла) Для балону з напруженням функціями форми є лінійними вздовж довжини елемента; загальні 2D та 3D континууmnі елементи використовують трикутні, чотирикутні, тетраедричні або пірамідальні функції форми з зростаючим поліноміальним порядком. Функції форми вищого порядку (p-уточнення) покращують точність без зменшення розміру сітки та можуть досягати експоненційного збіжності для гладких задач.

u(x) ≈ Σᵢ Nᵢ(x) · uᵢ   (within one element)

Key properties of any valid shape function set:
  Nᵢ(xⱼ) = δᵢⱼ   (1 at its own node, 0 at every other node)
  Σᵢ Nᵢ = 1      (partition of unity — rigid-body motion is exact)

Жорсткість елемента та збирання

Підставляючи функціонування форми у слабку форму, отримують матрицю жорсткості елемента K_e = ∫ Bᵀ C B dV, де B містить похідні від функцій форми, а C – матриця еластичності, побудована з модуля Юнга E та коефіцієнта Пуассона ν. Для простих елементів у вигляді балок і стержнів цей інтеграл має точний замкнений вигляд; для загальних елементів continuum він обчислюється чисельно за допомогою квадратури Гаусса. Матриця жорсткості кожного елемента потім збирається в глобальну систему K·u = F – для конструкції з балонів N та 2 ступенями свободи на балон, K має розмірність 2N×2N і є роздільною, оскільки балон з'єднаний лише з сусідніми елементами.

Обмежені умови та розв’язування

Для моделі мосту необхідні два типи обмежень: Dirichlet (основні) умови визначають відомі деформації, наприклад, фіксована опора або роликовий кронштейн з u = 0, які застосовуються шляхом модифікації рядків та стовпців матриці жорсткості; Neumann (вихідні) умови визначають прикладені сили, що безпосередньо з’являються у векторі навантаження F. Для лінійно-еластичної конструкції K·u = F розв'язується один раз за допомогою прямого методу (LU або розкладання Чолевеса) або ітеративного методу (кон’югованого градієнта) — повторна ітерація Ньютона-Рафсона не потрібна, якщо навантаження на дах значні або елементи можуть пластично деформуватися.

Довіра до результату: збіжність та похибка

Метод Обмеженого Елементу генерує приблизне рішення, тому точність перевіряється шляхом уточнення моделі: h-уточнення зменшує розмір елементів (похибка деформації приблизно пропорційна h², похибка напружень – h для лінійних елементів); p-уточнення підвищує порядок поліноміальності замість цього. Оцінювачі помилок, такі як метод Zienkiewicz-Zhu, порівнюють згладжене відновлене поле напружень з початковим неперервним FEM полем напружень через межі елементів, виділяючи області – наприклад, з’єднання в балку, що несе концентроване навантаження – де сітку потрібно уточнити. У практичній інженерії зазвичай прагнуть досягти оціненої похибки менше 5-10% перед тим, як довіряти обчисленому прикладеному силі або зміщенню для проєктування.

Часті запитання

Спробуйте наживо

Усе, що вище, працює прямо у вашому браузері — відкрийте the simulation і змінюйте параметри під час роботи. Нічого не встановлюється, нічого не завантажується на сервер, уся модель живе в одній вкладці.