Funkcja gęstości energii odkształcenia

Funkcja gęstości energii odkształcenia to matematyczny opis ilości energii zgromadzonej w jednostce objętości materiału w wyniku odkształcenia. W analizie MES stanowi ona fundament modelowania materiałów hipersprężystych, takich jak gumy, elastomery czy tkanki biologiczne, w których klasyczne prawo Hooke’a przestaje wystarczać.

Czym jest energia odkształcenia?

Gdy na obiekt działa siła lub gdy zostaje on podgrzany, dostarczona energia może powodować odkształcenie, zmianę temperatury lub inne efekty fizyczne. Energia zmagazynowana w odkształconym materiale w postaci sprężystego potencjału nazywana jest energią odkształcenia. Po usunięciu obciążenia może ona zostać oddana z powrotem — to ona przywraca odkształcony element do pierwotnego kształtu.

W rzeczywistych materiałach należy odróżnić rodzaje energii:

• energię odwracalną (sprężystą) — odzyskiwaną po odciążeniu,
• energię nieodwracalną — rozpraszaną w postaci ciepła, np. w plastyczności, tarciu, lepkości.

Funkcja gęstości energii odkształcenia opisuje wyłącznie część odwracalną.

Gęstość energii odkształcenia

Gęstość energii odkształcenia to energia odkształcenia przypadająca na jednostkę objętości materiału. Pozwala ona oderwać opis energetyczny od konkretnej geometrii i odnieść go bezpośrednio do właściwości materiału. W ujęciu lokalnym jest to wielkość polowa — przyjmuje różne wartości w różnych punktach odkształconego ciała.

Funkcja gęstości energii odkształcenia

Funkcja gęstości energii odkształcenia (oznaczana zwykle jako W) wyraża gęstość energii w postaci wzoru matematycznego, zależnego od stanu odkształcenia. Jej kluczową własnością jest to, że naprężenie w materiale otrzymuje się przez różniczkowanie W po odpowiedniej mierze odkształcenia. Innymi słowy, znajomość W wystarcza, by w pełni opisać sprężyste zachowanie materiału.

Funkcję W formułuje się najczęściej za pomocą:

• niezmienników tensora odkształcenia (I₁, I₂, I₃),
• współczynników wydłużenia głównego (λ₁, λ₂, λ₃),
• stałych materiałowych wyznaczanych eksperymentalnie.

Takie sformułowanie definiuje tzw. modele hipersprężyste, w których istnienie potencjału energii gwarantuje, że zachowanie materiału jest niezależne od drogi obciążenia.

Typowe modele hipersprężyste

W praktyce inżynierskiej stosuje się kilka dobrze ugruntowanych modeli funkcji gęstości energii odkształcenia:

• Neo-Hookean — najprostszy model jednoparametrowy, dobry dla małych i średnich odkształceń,
• Mooney-Rivlin — model dwu- lub wieloparametrowy, popularny dla gum,
• Ogden — oparty na współczynnikach wydłużenia, dobrze opisuje duże odkształcenia,
• Yeoh — trójparametrowy, sprawdza się przy dużych zakresach odkształcenia,
• Arruda-Boyce — model oparty na strukturze łańcuchowej polimerów,
• Polynomial / Reduced Polynomial — ogólne sformułowanie wielomianowe.

Wybór modelu zależy od zakresu odkształceń, dostępnych danych eksperymentalnych oraz typu zagadnienia.

Wyznaczanie parametrów materiałowych

Stałe materiałowe pojawiające się w funkcji W określa się eksperymentalnie. Standardowe testy parametrów materiałowych to:

• próba rozciągania jednoosiowego,
• próba dwuosiowa (biaxial),
• próba czystego ścinania (planar tension),
• próba ściskania.

Im więcej niezależnych prób, tym pewniejsze dopasowanie modelu. W programach MES dopasowanie parametrów odbywa się zazwyczaj automatycznie na podstawie wprowadzonych punktów krzywej eksperymentalnej.

Zastosowanie w MES

Funkcje gęstości energii odkształcenia są podstawą analizy konstrukcji wykonanych z materiałów hipersprężystych. Typowe zastosowania tej funkcji w MES obejmują:

• uszczelnienia gumowe i o-ringi,
• opony i elementy zawieszenia,
• wibroizolatory i poduszki silnika,
• elementy elastomerowe w przemyśle motoryzacyjnym i lotniczym,
• tkanki biologiczne w biomechanice,
• pianki i materiały komórkowe (z odpowiednimi modyfikacjami modelu).

Zalety i ograniczenia funkcji gęstości energii odkształcenia

Zalety:

• spójny opis dużych odwracalnych odkształceń,
• jeden potencjał W pozwala wyprowadzić zarówno naprężenia, jak i sztywność styczną,
• łatwa integracja z procedurami MES dla materiałów nieliniowych,
• dostępność wielu sprawdzonych modeli dla różnych klas materiałów.

Ograniczenia:

• opisuje tylko zachowanie odwracalne — nie uwzględnia plastyczności, pełzania ani uszkodzeń bez rozszerzeń modelu,
• wymaga starannego doboru modelu i jakościowych danych eksperymentalnych,
• niewłaściwy wybór modelu może prowadzić do niestabilności numerycznej,
• przy dużych odkształceniach konieczne bywa uwzględnienie nieściśliwości materiału.