Autor |
Analiza utraty stateczości |
kocyk
Użytkownik
Postów: 18
Data rejestracji: 05.06.13 |
Dodane dnia 09/19/2013 13:51 |
|
|
Witam, nie mogę znaleźć swojego starego tematu.
Wykonuję analizę utraty stateczności cienkościennych profili (kątownik oraz ceownik) dla różnych schematów podparcia w msc patran. Zależy mi na uzyskaniu ścieżek równowagi. wykorzystuję przy tym elementy quad4. Warunki brzegowe oraz obciążenia dodawane w środku ciężkości profili za pomocą MPC typu RBE2. Podczas analizy nieliniowej wykorzystuję metodę Newtona-Raphsona oraz arc-length.
Jednak w pliku .f06 otrzymuję taki komunikat:
^^^ USER INFORMATION MESSAGE 9005 (NLSTATIC)
^^^ THE SOLUTION FOR LOOPID= 622 IS SAVED FOR RESTART
*** USER WARNING MESSAGE 4698 (DCMPD)
STATISTICS FOR DECOMPOSITION OF MATRIX KLLRH.
THE FOLLOWING DEGREES OF FREEDOM HAVE FACTOR DIAGONAL RATIOS GREATER THAN
1.000000E+07 OR HAVE NEGATIVE TERMS ON THE FACTOR DIAGONAL.
USER INFORMATION:
THIS MESSAGE MAY BE IGNORED IF NO GRID POINT IDS OR HIGH RATIO MESSAGES APPEAR IN THE TABLE ON THE NEXT PAGE.
Czym to może być spowodowane?
Dziękuję za pomoc.
Edytowane przez kocyk dnia 09/19/2013 13:52 |
|
Autor |
RE: Analiza utraty stateczości |
dyku
Użytkownik
Postów: 124
Miejscowość: 3CITY
Data rejestracji: 13.05.08 |
Dodane dnia 09/20/2013 09:35 |
|
|
A jak masz podpartą belkę?
Może to jest wynik nie zablokowania któregoś stopnia swobody?
pozdr, |
|
Autor |
RE: Analiza utraty stateczości |
kocyk
Użytkownik
Postów: 18
Data rejestracji: 05.06.13 |
Dodane dnia 09/20/2013 14:16 |
|
|
Swoje analizy wykonuje dla schematu w ktorym jeden koniec jest zamocowany (odebrane wszystkie stopnie swobody), a drugi koniec swobodny (swoje analizy wykonuje rowniez dla innych schematow, ale na razie ograniczmy sie do tego). Katownik to 20x20 i symulacje przeprowadzam dla roznych wielkosci siatek w celu zobrazowania wplywu wielkosci siatki na otrzymane wyniki (siatka 20, 10, 5). Dodam, ze uwzgledniam nieliniowosc materialowa oraz geometryczna. Probowalem parametrem K6ROT zmieniac, wtedy tylko zmienial sie wezel w ktorym nastepuja te zbyt duzo przemieszczenia. Czy moze to wynikac z wielkosci siatki elememt skonczonych?
pozdrawiam:) |
|
Autor |
RE: Analiza utraty stateczości |
kojot
Użytkownik
Postów: 50
Miejscowość: Szczecin
Data rejestracji: 24.04.10 |
Dodane dnia 09/20/2013 20:09 |
|
|
W tutoraliach Ansysa czytałem że w nieliniowych analizach wyboczenia trzeba zadać małe obciążenie które faktycznie wyprowadzi układ z równowagi poza osiowym ( w rzeczywistości zawsze występuje jakaś nieosiowość obciążenia ). Może dopiero wtedy ciało zacznie sie wybaczać i problem zniknie. Nie znam się na Nastranie i niewiem czy faktycznie to jest jego źródłem.
Edytowane przez kojot dnia 09/20/2013 20:12 |
|
Autor |
RE: Analiza utraty stateczości |
kocyk
Użytkownik
Postów: 18
Data rejestracji: 05.06.13 |
Dodane dnia 09/21/2013 01:00 |
|
|
tzn model się wybacza, tylko właśnie występuje ten błąd o dużych translacjach lub obrotach w węźle i nie mam pojęcia czym to jest spowodowane... przy wyżej wspomnianym kątowniku natrafiam prawdopodobnie również na punkt bifurkacji, ponieważ na ścieżce równowagi dochodzi do siły granicznej, a następnie obciążenie oraz przemieszczenia maleją, ale to jest pewnie spowodowane złymi ustawieniami długości łuku.
Pozdrawiam. |
|
Autor |
RE: Analiza utraty stateczości |
Leonardo
Użytkownik
Postów: 90
Data rejestracji: 22.02.11 |
Dodane dnia 09/21/2013 08:29 |
|
|
Analiza wyboczeniowa (buckling) to teoretyczne wyznaczanie bifurkacyjnych punktów niestateczności konstrukcji perfekt. W jej wyniku otrzymuje się ciąg wartości własnych i odpowiadających im wektorów własnych. Kocyk natomiast wnioskując z opisu próbuje utworzyć ścieżkę równowagi. Aby to zadziałało powinna występować imperfekcja, czyli albo realna imperfekcja geometryczna albo zastępcze obciążenie poprzeczne. Bez tego nawet zastosowanie nieliniowości materiałowej i geometrycznej teoretycznie nie wystarczy. Chociaż tutaj pręt jest modelowany powłokami 4węzłowymi i w zależności od przyłożonego obciążenia zginanie może od razu wystąpić i wtedy zadziała.
Komunikat sugeruje za dużą różnicę wartości elementów macierzy sztywności, jakaś część zdaje się być geometrycznie zmienna, najprawdopodobniej problem leży w elementach rbe2. Parametr k6rot służy do sterowania 6 stopniem swobody w węzłach elementów powłokowych i związany jest z problemem obrotu węzła względem osi prostopadłej do płaszczyzny esa. Nie powinien mieć żadnego wpływu na problem, podobnie jak wielkość siatki. |
|
Autor |
RE: Analiza utraty stateczości |
kocyk
Użytkownik
Postów: 18
Data rejestracji: 05.06.13 |
Dodane dnia 09/26/2013 20:47 |
|
|
Witam,
Dalej borykam się z nieliniową analizą stateczności.
Mój badany kątownik bardzo dziwnie się zachowuje po przekroczeniu siły granicznej. Zamiast się przełamywać, powraca do swojej początkowej pozycji. czy to jest spowodowane wyłapaniem innej ścieżki równowagi...? Poniżej zamieszczam moje rysunki. Model, charakter pracy mojego kątownika, oraz przemieszczenia przy sile granicznej oraz dla ostatniej iteracji.
Proszę o pomoc.
Pozdrawiam.
Ł.K.
|
|
Autor |
RE: Analiza utraty stateczości |
piotr
Super Administrator
Postów: 1432
Data rejestracji: 17.05.07 |
Dodane dnia 09/27/2013 21:51 |
|
|
Czy mógłbyś to opisać trochę dokładniej. Nie bardzo rozumiem co się dzieje i co przedstawiają zamieszczone przez ciebie rysunki. Piszesz, że pręt wrócił do pierwotnej pozycji ale wykres temu przeczy (chyba)... |
|
Autor |
RE: Analiza utraty stateczości |
kocyk
Użytkownik
Postów: 18
Data rejestracji: 05.06.13 |
Dodane dnia 10/01/2013 13:33 |
|
|
Panie Piotrze na wykresie bardzo dobrze widać, że na początku wartość przemieszczenia rośnie, aż do uzyskania siły granicznej (w tej fazie kątownik się wybacza co przedstawione jest na powyższym rysunku 4), po przekroczeniu siły granicznej przemieszczenie maleje aż do punktu zero, co również widać na wykresie (czyli powraca, zamiast się przełamywać dalej, rysunek 5 przemieszczenia dla ostatniej wartości iteracji). Przemieszczenie powinno dalej rosnąć, dlatego wstawiłem czerwoną strzałkę na wykresie, aż do zniszczenia takiego kątownika. |
|
Autor |
RE: Analiza utraty stateczości |
piotr
Super Administrator
Postów: 1432
Data rejestracji: 17.05.07 |
Dodane dnia 10/01/2013 18:53 |
|
|
A nie jest to po prostu tak, że najpierw obciążenie rośnie, a potem maleje? Wskazywałby na to "histerezowy" wygląd wykresu. Czyli, że cykl obliczeniowy zawiera w sobie narastanie obciążenia do pewnej wartości a potem tej wartości zmniejszanie? |
|
Autor |
RE: Analiza utraty stateczości |
kocyk
Użytkownik
Postów: 18
Data rejestracji: 05.06.13 |
Dodane dnia 10/01/2013 19:23 |
|
|
Zmieniłem rodzaj przeliczania macierzy z automatycznej na obliczaną w każdym kroku iteracyjnym i muszę przyznać, że uzyskałem fajny wynik. Jak znajdę więcej czasu to postaram się zamieścić nowe rysunki.
Pozdrawiam
Ł.K. |
|
Autor |
RE: Analiza utraty stateczości |
kocyk
Użytkownik
Postów: 18
Data rejestracji: 05.06.13 |
Dodane dnia 10/01/2013 22:31 |
|
|
Przesyłam animację. Jednak jeszcze chyba coś nie gra w tej analizie....
http://www.youtube.com/watch?v=GbcOA9DE9yQ&feature=youtu.be |
|
Autor |
RE: Analiza utraty stateczości |
Leonardo
Użytkownik
Postów: 90
Data rejestracji: 22.02.11 |
Dodane dnia 10/02/2013 00:45 |
|
|
Zachowanie tego wspornika jest raczej absurdalne. Krzywa odkształcenia wykazuje podejrzaną podwójną krzywiznę widoczną w pobliżu wierzchołka. W początkowej fazie przyrostu obciążenia można zauważyć również cofnięcie się deformacji niczym "przeskoku" i ponowny wzrost ugięcia. Na koniec całość powraca do stanu wyjściowego.
Jeśli analiza prowadzona jest dla dużych odkształceń i z uwzględnieniem modelu sprężysto-plastycznego niemożliwa jest histereza, o której wspominał Piotr, gdyż po przekroczeniu granicy plastyczności powstały by już odkształcenia trwałe, które nie zniknęły by nawet przy usunięciu obciążenia.
Istotne może być wyjaśnienie sytuacji na rysunku 2. Siła pionowa 40kN zaczepiona jest do węzła, który znajduje się w środku ciężkości przekroju, środku ciężkości węzłów ujętych "pająkiem"? |
|
Autor |
RE: Analiza utraty stateczości |
kocyk
Użytkownik
Postów: 18
Data rejestracji: 05.06.13 |
Dodane dnia 10/02/2013 11:01 |
|
|
Zaczepiona jest w środku ciężkości przekroju.
Po przestawieniu rodzaju przeliczania macierzy otrzymałem poniższy wykres, który jest lepszy, ale dalej mnie niepokoi.
|
|