10 Dynamika konstrukcji, Budowlane, MES, MES. Łodygowski
[ Pobierz całość w formacie PDF ]
10. DYNAMIKA KONSTRUKCJI
1
10.
10. DYNAMIKA KONSTRUKCJI
10.1. Wprowadzenie
Ogólne równanie dynamiki zapisujemy w postaci:
M d
t
C d
t
Kd
t
=
P
t
(10.1)
Zapis powyższy oznacza, że równanie musi być spełnione w każdej chwili. Jest to odwołanie do
zasady d'Lamberta, która mówi, że dla układu będącego w ruchu równowaga musi być spełniona w każdej
chwili konkretnej przestrzeni czasowej.
Macierz
M
jest
macierzą masową, macierz
K
- macierzą sztywności. Macierz
P
t
jest macierzą
określającą przyłożone do układu obciążenia zewnętrzne. Natomiast C jest macierzą określającą tłumienie
układu. Macierz tą przyjmujemy najczęściej w postaci tzw. tłumienia proporcjonalnego (zależnego od
macierzy
K
i
M
) w postaci
C
=
1
M
2
K
(10.2)
Współczynniki
1
i
2
wyznaczamy na podstawie udziału poszczególnych postaci drgań własnych.
Jeśli założymy wartość tłumienia i obciążenia zewnętrzne równe zero, otrzymamy równanie
M d
t
Kd
t
=0
(10.3)
czyli problem drgań własnych układu. Idąc dalej, stosując podstawienie
d
t
=
d
0
⋅sin
t
(10.4)
i różniczkując dwukrotnie po czasie
t
d
t
=
d
0
⋅cos
t
d
t
=−
2
d
0
⋅sin
t
(10.5)
otrzymujemy wartości, które podstawiamy do równania (10.3) i dostajemy
K
−
2
M
d
0
=0
(10.6)
Czyli równanie, które definiuje nam uogólniony problem własny. Równanie to ma
n
równań
rzeczywistych w postaci par: wartości własnej i odpowiadającego jej wektora własnego.
J.Gieczewski, M.Kończl, A.Krzysztoń, D.Mejbaum, N.Roszak, M.Wojciechowski, J.Wojtkowiak
AlmaMater
10. DYNAMIKA KONSTRUKCJI
2
10.2. Metody całkowania
Jak wiemy równanie ruchu jest równaniem różniczkowym. Zastanówmy się zatem nad sposobami
jego rozwiązywania. Ze względu na pewne własności macierzy
K
,
M
i
C
w analizie ruchu ciała
przedstawionego przy pomocy elementów skończonych zasadniczo stosujemy dwie metody: metodę
całkowania bezpośredniego i metodę superpozycji modalnej.
10.2.1. Metody całkowania bezpośredniego
Metody całkowania bezpośredniego są metodami jawnymi. Polegają ona na tym, że równanie ruchu
jest całkowane krok po kroku. Równanie ma być spełnione tylko w wybranych chwilach „t”, a nie w całym
przedziale całkowania.
Zakładamy, że w chwili t=0 znane nam są przemieszczenia, prędkości i przyspieszenia (czyli znamy
d
0,
d
0,
d
0
). Rozpatrujemy określony przedział czasowy (0,T), który dzielimy na n równych przedziałów, w
których to poszukujemy naszych nieznanych wielkości. Rozważamy zatem chwile:
0,
t ,
2
t ,
...
, t , t
t ,
...
, T
(10.7)
Zadanie polega na zbudowaniu algorytmu, który pozwoli nam na obliczeniu poszukiwanych wartości
w danym kroku przy wykorzystaniu wyliczonych wartości z kroku poprzedniego. W taki sposób otrzymamy
wartości we wszystkich chwilach czasowych z przedziału
(0,T)
Pokażemy na przykładzie, tok postępowania przy rozwiązywaniu zadania z dynamiki przy pomocy
jednej z najbardziej efektywnych metod z grona metod całkowania bezpośredniego, a mianowicie
metodą
różnic centralnych
.
Zakładamy zmienność w czasie wektorów prędkości i przyspieszeń w postaci
d
t
≃
1
2
d
t
t
−
d
t
−
t
t
(10.8)
1
d
t
≃
t
2
d
t
t
−2
d
t
d
t
−
t
Jeśli podstawimy operatory różnicowe (10.8) do (10.1) otrzymamy
t
2
d
t
t
−2
d
t
d
t
−
t
M
2
d
t
t
−
d
t
−
t
t
C
Kd
t
=
P
t
(10.9)
Z równania (10.9) obliczamy poszukiwany stan przemieszczeń w chwili
t
t
czyli
d
t
t
.
J.Gieczewski, M.Kończl, A.Krzysztoń, D.Mejbaum, N.Roszak, M.Wojciechowski, J.Wojtkowiak
AlmaMater
1
10. DYNAMIKA KONSTRUKCJI
3
Należy zwrócić uwagę, że rozwiązanie to otrzymujemy na podstawie rozwiązania w chwili
t
. Stąd też
metodę tę zaliczamy do metod jawnych (explicite). Dużą zaletą tego sposobu rozwiązywania równania
(10.9) jest fakt, iż nie musimy odwracać macierzy sztywności.
Należy zwrócić uwagę, że obliczanie wyników w kolejnych chwilach z wykorzystaniem wyników
otrzymanych w chwilach poprzednich wymaga przyjęcia pewnej procedury startowej. Warto zaznaczyć, że
zakładamy tutaj, iż wektory
d
0,
d
0,
d
0
są znane w chwili początkowej czyli w chwili
t=0
. Stąd też
wykorzystując wzory (10.8) możemy wyznaczyć wektor przemieszczenia
d
w fikcyjnej chwili, która
poprzedzać będzie początek ruchu czyli dla chwili
t
−
t
:
d
t
−
t
=
d
0
−
t d
0
t
2
2
d
0,
(10.10)
Zaznaczmy, że metody jawne są tylko warunkowo stabilne, dlatego też wymagane jest zastosowanie
małych kroków całkowania
t
przy obliczaniu kolejnych wartości. Krok nie może być dowolnie duży,
lecz musi spełniać poniższą zależność
t
t
kr
=
T
n
,
(10.11)
gdzie
T
n
jest najmniejszym okresem drgań układu
Niespełnienie warunku (10.11) powoduje narastanie – akumulację błędów całkowania i zaokrągleń w
trakcie rozwiązywania równań ruchu.
Algorytm obliczeń dla metody całkowania jawnego:
•
Obliczamy macierze
K
,
C
,
M
•
Następnie obliczamy
d
0,
d
0,
d
0,
•
Określamy stałe
a
0
=
t
2
a
1
=
1
2
t
a
2
=2
a
0
a
3
=
1
(10.12)
a
2
•
Obliczamy
d
t
−
t
=
d
0
−
t d
0
t
2
2
d
0
•
Wyznaczamy
M
M
=
a
0
M
a
1
C
(10.13)
•
Triangularyzacja macierzy
M
przy pomocy wzoru
M
=
LDL
T
(10.14)
J.Gieczewski, M.Kończl, A.Krzysztoń, D.Mejbaum, N.Roszak, M.Wojciechowski, J.Wojtkowiak
AlmaMater
1
10. DYNAMIKA KONSTRUKCJI
4
•
Obliczenia dla każdego kroku:
- wektora obciążenia efektywnego
R
R
=
R
t
−
K
−
a
2
M
d
t
−
a
0
M
−
a
1
C
d
t
−
t
(10.15)
- rozwiązanie równania (10.9) dla chwili
t
t
L
T
DLd
t
t
=
R
(10.16)
-obliczenie wektorów prędkości i przyspieszeń:
d
t
≃
a
1
d
t
t
−
d
t
−
t
d
t
≃
a
0
d
t
t
−2
d
t
d
t
−
t
(10.17)
W przypadku braku tłumienia czyli gdy
C
=0
, równanie (10.9) upraszczamy do postaci
t
2
M
d
t
t
=
R
t
(10.18)
gdzie
R
=
R
−
K
−
t
2
M
2
d
−
−1
t
2
M
d
(10.19)
Jeśli w równaniu (10.18) macierz mas będzie diagonalna, wówczas rozwiązania otrzymujemy poprzez
wykonanie określonego wzorem (10.19) mnożenia
d
t
t
i
=
R
t
i
t
2
(10.20)
m
ii
oraz
R
t
i
będą oznaczać i-te składowe wektorów
d
t
t
i
R
t
, natomiast
m
ii
odnoszą się do i-tej składowej diagonalnej macierzy mas (należy jednakże spełnić założenie, że
m
ii
0
.
i
Zauważmy, że nie musimy znać macierzy globalnych (zarówno macierzy sztywności
K
jak i macierzy
mas
M
). Dzieje się tak dlatego, bo nie rozwiązujemy układu równań liniowych. Macierze
K
i
M
mogą być
określone na poziomie elementów.
Wśród metod całkowania jawnego możemy wymienić, oprócz metody różnic centralnych, między
innymi metodę Houbolta, Wilsona i Newmarka. Metody te, pod warunkiem przyjęcia pewnych wartości
współczynników charakterystycznych dla danej metody, należą do metod bezwarunkowo stabilnych.
J.Gieczewski, M.Kończl, A.Krzysztoń, D.Mejbaum, N.Roszak, M.Wojciechowski, J.Wojtkowiak
AlmaMater
1
gdzie
d
t
t
10. DYNAMIKA KONSTRUKCJI
5
10.2.2. Metody superpozycji modalnej
Równanie ruchu ma postać:
M d
t
C d
t
Kd
t
=
P
t
(10.21)
M d
t
t
C
˙
d
t
t
Kd
t
t
=
P
t
t
(10.22)
d
t
t
=
d
t
[
1 −
d
t
d
t
t
]
t
˙
(10.23)
d
t
t
=
d
t
d
t
t
2
[
2
−
d
t
d
t
t
]
ℵ
(10.24)
,
- parametry przyjmowane na podstawie rozwiązań dotyczących dokładności i stabilności
otrzymanych rozwiązań
=
1
6
=
1
2
Rozwiązując
ℵ
względem
d
t
t
otrzymamy:
d
t
t
=
1
t
2
[
d
t
t
−
d
t
−
d
t
t
−
t
2
2
−
d
t
]
(10.25)
d
t
t
=
t
d
t
t
−
d
t
1−
d
t
t
1−
2
d
t
(10.26)
1
t
2
[
d
t
t
−
d
t
−
d
t
t
−
t
2
2
−
d
t
]
⋅
M
=
P
t
t
(10.27)
Z równania tego obliczamy niewiadomy wektor przemieszczeń
d
t
t
i podstawiamy do
d
t
t
i
d
t
t
Jeśli liczba kroków i liczba stopni swobody układu jest duża, wówczas efektywność obliczeń
metodami całkowania bezpośredniego jest niesatysfakcjonująca. Należy wtedy posłużyć się innymi
metodami - metodami niejawnymi (implicite), do których można zaliczyć metodę superpozycji modalnej.
Należy tutaj dokonać przekształcenia równania równowagi (10.1) do postaci, która będzie wymagała
od nas mniejszego nakładu pracy.
J.Gieczewski, M.Kończl, A.Krzysztoń, D.Mejbaum, N.Roszak, M.Wojciechowski, J.Wojtkowiak
AlmaMater
[ Pobierz całość w formacie PDF ]