(function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start': new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0], j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src= 'https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f); })(window,document,'script','dataLayer','GTM-5M6SH59');
8 minutes read
26 marca 2017

Analiza geometrycznie nieliniowa – jak to działa?

8 minutes read

Analiza geometrycznie nieliniowa jest niesamowicie przydatna w analizie konstrukcji. Pomaga sprawdzić, czy wprowadzony model jest poprawny i pozwala na dużą optymalizację rozwiązania. Dziś omówię podstawy tej analizy – nie teorię, lecz aspekty praktyczne.

Czego się dziś dowiesz

Wyjaśnię, co tak naprawdę robi analiza geometrycznie nieliniowa. Obiecuję, że nie będzie tu złożonej teorii. Zamiast tego będą proste przykłady z życia wzięte.

Oczywiście wybrałem takie przykłady, które w analizie liniowej i nieliniowej dają różne wyniki. Istnieją również przykłady, dla których oba podejścia dają praktycznie takie same rezultaty. Należy jednak pamiętać, że istnieje wiele takich przypadków, które wymagają podejścia nielinowego i właśnie tymi będziemy zajmować się dzisiaj.

Grunt to dobrze zacząć: Pranie!

Uwielbiam proste życiowe przykłady. Jeśli chodzi o analizę nieliniową, jej ideę najlepiej obrazuje sznurek z praniem.

Geometrically linear analysis

Spójrz na rysunek powyżej. Mokry sweter waży… 5 kg? Zakładając, że sznurek ma długość 5 m, możemy obliczyć moment zginający (0.0625 kNm). Wydaje się być mały, prawda? Weźmy jednak pod uwagę, że sznurek ma 5 mm średnicy (zakładając, że to naprawdę mocny sznurek… mój ma około 3…). Jego wskaźnik na zginanie wynosi więc ok. 0.0123 cm^3. Oznacza to, że gdyby sznurek chciał przenieść obciążenie od swetra za pomocą zginania, musiałby wytrzymać naprężenia 0.0625kNm/1.23e-8 m^3 = 5081MPa… czyli 20 razy większe, niż wytrzymuje stal.

Kiedy patrzę na mój sznurek od prania, raczej wątpię, by mógł przenieść takie naprężenia… a jednak ciągle działa.

Przytoczony tu przykład to typowe liniowe podejście do problemu. Być może zauważyłeś, że mój przykład trochę oszukiwał – nie narysowałem żadnego ugięcia. W rzeczywistości łatwo byłoby je obliczyć. Gdyby sznurek był zrobiony ze stali (to założenie nieco odległe od rzeczywistości), ugięcie wyniosłoby… 20.2 m.

Pomyśl tylko, jak bardzo biedny sznurek musiałby się rozciągnąć, żeby się tak odkształcić (pięciometrowy sznurek uginający się 20 m?!). To zagadnienie przybliża nas do idei analizy nielinowej, która stosuje teorię dużych ugięć.

Kiedy sznurek staje się nieliniowy

Sznurek nie może się wydłużyć o 20.2 m, jak to wyszło z podejścia linowego. Trochę się jednak wydłuży. Jak pewnie wiesz, to wydłużenie spowoduje powstanie siły rozciągającej w sznurku. Z kolei ta siła będzie stabiliozwać cały układ, zmniejszając ugięcie.

Geometrically nonlinear analysis - basic example

Obliczając ugięcie zgodnie z teorią dużych ugięć, otrzymamy wynik 57 mm (zakładając ciągle, że mamy stalowy sznurek). Dla nylonu ugięcie to wyniosłoby ok. 237 mm, natomiast siła osiowa wyniosłaby 0.26 kN (dla stali wynosi 1.07 kN).

Łatwo mogę uwierzyć, że mój sznurek jest w stanie przenieść w ten sposób jakieś 26 kg.

Jak działa analiza geometrycznie nieliniowa?

Geometrically nonlinear analysis - how does membrane state works

Spójrz na schemat powyżej. Załóżmy, że przecinamy sznurek na środku. Jest tam siła osiowa N (którą przed chwilą obliczyliśmy) ciągnąca sznurek „do środka”. W związku z tym w ścianie, w której jest umocowany sznurek, pojawia się pozioma reakcja również równa N. Obie siły tworzą parę sił, która może przenosić moment. Jest to pewne uproszczenie, dobrze jednak opisuje, co się tak naprawdę dzieje.

Łatwo można sprawdzić to, co napisałem. W obliczeniach liniowych otrzymaliśmy moment zginający 0.0625 kNm. W analizie nielinowej ugięcie 57 mm i siła osiowa 1.07 kN dają nam 0.057 m x 1.07 kN = 0.061 kNm. W przypadku nylonowego sznurka mamy natomiast 0.237 m x 0.26 kN = 0.0616 kNm. Oba wyniki nie zgadzają się idealnie, ale są dość bliskie, co po prostu pokazuje, jak należy myśleć o tym zagadnieniu.

Co należy zapamiętać odnośnie stanu cięgnowego:

  • Im większe odkształcenie, tym mniejsze siły osiowe (to dlatego w nylonowym sznurku spadła siła).
  • Teoretycznie wszystko dąży, by odkształcić się jak najmocniej w celu redukcji obciążeń.
  • Istnieje jednak ograniczenie tego, jak mocno element może się odkształcić. Materiał ma pewien moduł Younga, a przekrój poprzeczny ma określone pole i moment bezwładności. Dlatego odkształca się tyle, ile może, a następnie się zatrzymuje.
  • Jeśli odkształcenia są „małe”, belka (trudno nazwać to sznurkiem, ponieważ jest sztywna) będzie przenosić duże momenty zginające i małe siły osiowe. Dokładnie jak w analizie liniowej.
  • Jeśli odkształcenia są „duże”, moment zginający będzie mały, natomiast siły osiowe duże.
  • Zjawisko to jest oczywiście „ciągłe”, więc jeśli odkształcenia są „średnie”, występuje i trochę zginania, i trochę rozciągania.

Uwaga! Jest i haczyk!

W tym miejscu ktoś mógłby dojść do wniosku, że analiza nieliniowa ma tylko pozytywne aspekty, więc jej „ignorowanie” jest po stronie bezpiecznej. Niestety trzeba pamiętać:

W projektowaniu nic nigdy nie jest za darmo

Tutaj „płacimy” za ten miły efekt cięgnowy (który zwiększa nośność sznurka), otrzymując siłę rozciągającą w podporach. Chcesz sprawdzić, że to prawda? Spróbuj powiesić pranie na sznurku, który nie jest przymocowany do ściany, tylko leży swobodnie na stole. Dla typowej belki nie byłby to problem (stół wciąż przenosi siły pionowe). Ale wszyscy intuicyjnie wiemy, jak skończy się taki eksperyment.

Jest to też ostateczny dowód na to, że w podporach pojawia się wyrywanie – w końcu to ono właśnie sprawia, że system działa.

Charakter analizy geometrycznie nieliniowej

Geometryczna nieliniowość nie jest ani ściśle pozytywna, ani negatywna… po prostu jest!

Difference between linear and geometrically nonlinear analysis

Wyobraź sobie, że masz sznurek z praniem oparty na dwóch słupach. Jeśli przeanalizujesz ten problem w sposób liniowy, otrzymasz wyniki „złe” albo przynajmniej „konserwatywne”. Jest to spowodowane tym, że siły, które wyznaczysz w podejściu linowym, będą znacznie przeszacowane. Oznacza to, że zaprojektujesz sznurek bezpiecznie, ale bardzo konserwatywnie.

Z drugiej strony dla słupów podejście liniowe jest bardzo optymistyczne. Pokaże ono jedynie siły pionowe, całkowicie pomijając poziome. Oznacza to, że pominięte zostanie zginanie słupów.

Jak widać, geometryczna nieliniowość „pomaga” sznurkowi i „obciąża” słupy. Niestety nie da się wybrać, czy chce się jej używać, czy nie… te zjawiska po prostu tu są!

Tutaj potrzebne jest słowo komentarza. Zakładając, że wykonujemy linową statykę, projektujemy sznurek tak, by był w stanie przenieść moment zginający (obliczony jak dla belki). W takim przypadku nie mamy już do czynienia ze „sznurkiem”, a raczej z „belką”. To z kolei prowadzi do znacznego zmniejszenia deformacji i sił osiowych obciążających słupy. Być może w takim przypadku siły osiowe byłyby nawet pomijalnie małe, wciąż jednak by były.



Dlatego obliczenia statyczne są możliwe:

Jeśli odkształcenia w modelu są małe, również efekty geometrycznej nieliniowości będą niewielkie. Jedynym problemem jest to, że czasami trudno zgadnąć, czy deformacje są „wystarczająco małe”, by zignorować geometryczną nieliniowość.

W tym przykładzie warto jeszcze coś zauważyć. Jeśli słupy byłyby mocno obciążone siłą pionową, odkształcenia obserwowane w analizie nieliniowej spowodowałyby wzrost mimośrodu w słupach. To z kolei wywołałoby ich dodatkowe zginanie (nazywane często zginaniem drugiego rzędu), co w niektórych przypadkach może być bardzo istotne. To również jest efekt geometrycznej nieliniowości.

Podsumowując

Analiza geometrycznie nieliniowa jest bardzo przydatnym narzędziem w obliczeniach statycznych. Jeśli masz do czynienia z elementami, które się sporo odkształcają, jest to najlepsze podejście. Jeśli analizujesz przypadek, w którym podejście liniowe jest wystarczające, z analizy nieliniowej otrzymasz takie same wyniki jak z analizy liniowej. Jeśli jednak zadanie wymaga użycia analizy nielinowej, nie da się jej zastąpić liniową.

Analiza nieliniowa świetnie sprawdza się też w przypadkach zginania drugiego rzędu, gdzie odkształcenia modelu powodują zwiększenie sił prostopadłych wywołujących dodatkowe zginanie.

W obliczeniach statycznych popularny jest również dylemat dotyczący liniowej/nieliniowej analizy wyboczeniowej, mocno związana z zagadnieniami analizy geometrycznie linowej/nieliniowej. Temat ten opisałem w innym miejscu. Tutaj skupiłem się głownie na naprężeniach.

Darmowy kurs MES!

Jeśli udało mi się zainteresować Cię tematem koniecznie zapisz się na darmowy kurs o MES i stateczności. Możesz się do niego zapisać poniżej:

Author: Łukasz Skotny Ph.D.

Mam ponad 10 lat doświadczenia w praktycznym wykorzystaniu MES w projektowaniu (prowadzę własne biuro projektowe), a do tego przez dekadę byłem wykładowcą na Politechnice Wrocławskiej. Obecnie tutaj dzielę się swoją wiedzą z inżynierii i MES dzięki kursom oraz na blogu!

Read more

Join my FEA Newsletter

Get my 1h video Lecture on Nonlinear Material

    Your personal data administrator is Enterfea Łukasz Skotny, Skrzydlata 1/7, 54-129 Wrocław/POLAND, Email. By subscribing to the newsletter that includes marketing messages you consent to your personal data processing in accordance with this privacy policy

    Join the discussion

    Comments (18)

    Maciek - 2018-09-11 06:53:47

    i wszystko jasne, myślałem że się nie na studiach czegoś nie douczyłem :) Dziękuję za wyjaśnienie :)
    Pozdrawiam

    Reply
    Łukasz Skotny Ph.D. - 2018-09-11 07:21:02

    :)

    Reply
    Maciek - 2018-09-07 09:58:42

    Witam:)
    po raz kolejny gratuluję świetnego bloga i prostego wyjaśnienia skomplikowanej sprawy. Mam tylko jeden mały problem. Nie do końca rozumiem tylko czym jest ta magiczna teoria dużych ugięć. Albo stwierdzono, że na studiach mechanicznych tego uczyć nie trzeba albo pechowo musiałem ominąć jakiś wykład. Czy mógłbyś mi podać jakąś książkę, link... gdzie odkryję tą tajemnicę?

    Pozdrawiam

    Reply
    Łukasz Skotny Ph.D. - 2018-09-07 13:59:51

    Hej!

    Uj... z książką to bym nie rysykował :P Znaczy zawsze możesz spróbować Crisfielda (Nonlinear analysis of soils and structures)... ale to bolesne doświadczenie o ile mi wiadomo :)

    No więc ta teoria dużych ugięć to jest właśnie to, co sprawia że masz geometrycznie nieliniową analizę. Oczywiście tam są dodatkowe zagadnienia jak duży strain czy duże kąty, ale ogólnie jak to włączysz to masz nieliniową geometrię, a jak nie to masz liniową :)

    pozdr
    Ł

    Reply
    Maciek - 2018-09-11 06:18:14

    Chodzi mi w jaki sposób obliczyłeś ugięcie równe 57 mm. Analitycznie ( są na to jakieś konkretne wzory) , czy za pomocą programu ?

    Reply
    Łukasz Skotny Ph.D. - 2018-09-11 06:40:35

    Hej!

    Nie no, liczyłem to za pomocą programu oczywiście. :) Wydaje mi się że akurat taki temat jak ten sznurek dałoby się ogarnąć analitycznie (są krzywe ugięć kabli opisane), ale po co się tak męczyć? Wprowadzasz to do programu, robisz geometrycznie nieliniową analizę i... świętujesz zwycięstwo :)

    Pozdr
    Ł

    Reply
    Q - 2018-08-14 20:42:54

    Witam!

    Na wstępie gratuluje bloga, bo takie podejście pozwala na odczarowanie uniwersyteckiego podejścia do nauczania, które niestety miałem okazje przerobić i przeboleć. Całe szczęście na swojej drodze spotkałem kilku wykładowców którzy pokazali że można inaczej, a Ty Łukaszu robisz to samo i to jeszcze na blogu - chwała Ci za to! Jako że zaczynam projektować stal w MESie, a nie chcę robić tego po łebkach (na zasadzie wrzucania danych bez zrozumienia jak wpływają one na obliczenia i otrzymane wyniki) jest to pomoc nieodzowna i oszczędzająca wiele czasu i nerwów.

    Odnośnie wpisu - po przeczytaniu nasuwa mi się jedno pytanie, a mianowicie czy słuszne jest stwierdzenie, że gdybyśmy bardzo znacząco zmniejszyli działającą siłę (bez zmiany rozstawu podpór i przekroju sznurka) to sznurek będzie pracował jak belka - tj. będzie przenosił moment gnący bo odkształcenie od działającej siły będzie zmierzało do 0. I ogólniej - to czy element przenosi moment gnący i małe siły osiowe lub mały moment gnący i duże siły osiowe zależy od cech geometrycznych, cech materiałowych, warunków podparcia ale również od tego jakie siły na niego działają?

    Pozdrawiam,
    Q

    Reply
    Łukasz Skotny Ph.D. - 2018-08-17 14:45:16

    Hej!

    Dzięki za miłe słowa. Pytanie które stawiasz jest super - jedno z lepszych jakie dotej pory przeczytałem :)

    Tak, gdybyś mega zmniejszył obciążenia sznurek byłby belką... tylko że z ciężarem własnym niewiele zrobisz, więc jest "praktyczny limit" tego jak długo to działa. Po prostu w pewnym momencie już sam ciężar własny wystarcza... nawet jak nie ma w ogóle obciążenia, no i tutaj już zmniejszanie obciążenia nic nie da :p

    Geometryczna nieliniowość to złożona sprawa, i zależy od wielu rzeczy, na pewno od cech geometrycznych, również od materiału (choć raczej nie granicy plastyczności... no chyba że coś ci się zacznie uplastyczniać gdzieś po drodze :P). Zdecydowanie istotne są warunki podparcia (ostatecznie swobodnie odparta belka z przesuwem z jednej strony... nigdy się nie uciągli, bo nie ma jak!).

    Co do obciążeń to o ile mówimy o ich wartościach, to w zasadzie one nie odgrywają roli... chodzi jedynie o to że przy "małych" obciążeniach wyniki z analizy liniowej i nieliniowej są do siebie bardzo zbliżone :)

    pozdr
    Ł

    Reply
    Szymon - 2018-01-24 12:56:50

    Zacząłem od manuala zsoilu i skończyłem na tym blogu. Jest świetny ;)

    Pozdrawiam!

    Reply
    Łukasz Skotny Ph.D. - 2018-01-26 11:28:40

    Cześć Szymon!

    Cieszę się że blog Ci się podoba :)

    Pozdr
    Ł

    Reply
    Daria - 2018-01-09 23:03:21

    Witam,
    Super wytłumaczenie zagadnień, które studentom sprawiają wielkie trudności. Twoje podejście jest przykładem na to, że najtrudniejsze zagadnienia da się w prosty sposób wytłumaczyć nie używając przy tym stosu wzorów odwracających uwagę. Jak byłoby pięknie jakby na studiach Wykładowcy w tak przystępny sposób tłumaczyli nie proste zagadnienia....

    Reply
    Łukasz Skotny Ph.D. - 2018-01-10 09:24:34

    Cześć Daria!

    Bardzo mi miło to słyszeć - cieszę się że post Ci się spodobał.

    Z Uczelniami to długi temat... lekko nie jest to prawda, ale są przecież ludzi którzy uczą naprawdę świetnie. Czasem jednak trzeba się ich naszukać : )

    Pozdrawiam serdecznie
    Ł

    Reply
    Wojtek - 2017-11-29 12:32:50

    Również uważam że to super proste wytłumaczenie zasady działania bardziej skomplikowanego mechanizmu. Tak się powinien zaczynać wykład z analizy nieliniowej na studiach:)

    Reply
    Łukasz Skotny Ph.D. - 2017-11-29 15:52:22

    Cześć...

    ... gdzieś ty studiował! Ja nie miałem wykładów z nieliniowych analiz... (chyba :P)

    Cieszę się że Ci się podoba. Uważam że dosłownie wszystko da się tak wytłumaczyć :P Z resztą takich artykułów jest więcej - będę je jakoś sukcesywnie podsyłał :P

    Pozdr
    Ł

    Reply
    Piotr - 2017-10-28 09:56:00

    Super artykuł! Czapka z głowy :) Pozdrawiam.

    Reply
    Łukasz Skotny Ph.D. - 2017-10-28 10:35:34

    Cześć Piotr!

    Strasznie się cieszę że Ci się podobał - zawsze miło to usłyszeć :)

    Pozdr
    Ł

    Reply
    Piotr - 2017-09-29 07:38:07

    Bardzo fajne wytłumaczenie analizy nieliniowej :) Jednak można wytłumaczyć trudne przypadki prostymi przykładami:) Pozdrawiam

    Reply
    Łukasz Skotny Ph.D. - 2017-10-07 04:00:46

    Cześć Piotr :)

    Cieszę się że post przypadł Ci do gustu! Dzięki za miłe słowa :)

    Pozdr
    Ł

    Reply

    Zapisz się na mój Newsletter o MES (po angielsku)

    Co tydzień otrzymasz maila z anglojęzycznymi materiałami o MES

      Your personal data administrator is Enterfea Łukasz Skotny, Skrzydlata 1/7, 54-129 Wrocław/POLAND, Email. By subscribing to the newsletter that includes marketing messages you consent to your personal data processing in accordance with this privacy policy