Politechnika Poznańska Zastosowanie algorytmów genetycznych do ...
Politechnika Poznańska Zastosowanie algorytmów genetycznych do ... Politechnika Poznańska Zastosowanie algorytmów genetycznych do ...
M T e (x, y) = T e j ∙ N e j (x, y), (4.11) j=1 gdzie M odpowiada liczbie węzłów przyjętego elementu e, T e j oznacza temperaturę w węźle j elementu e, N e j to funkcja kształtu elementu e. Stosując przedstawioną wcześniej metodę Galerkina, możemy zapisać: ∂ ∂Te λ ∂x ∂x + ∂ ∂Te λ ∂y ∂y + Q̇ N e i (x, y)dxdy = 0. (4.12) Ω e Korzystając z twierdzenia Greena, otrzymujemy następującą zależność: ∂ ∂Te λ ∂x Ω e ∂x N i e + ∂ ∂y λ ∂Te ∂y N i e dxdy = (4.13) = λ ∂Te ∂x dy − λ ∂Te ∂y dx N i e Γ e . Przekształcając lewą stronę równania otrzymujemy: ∂ ∂Te λ ∂x Ω e ∂x + ∂ ∂y λ ∂Te ∂y N i e dxdy = (4.14) = − λ ∂Te e e ∂N i ∂Te ∂N i ∂Te ∂Te + λ dxdy + λ dy − λ ∂x ∂x ∂y ∂y ∂x ∂y dx N i e , Ω e Γ e Biorąc pod uwagę, że Q̇ = λ ∂T ∂x n x + λ ∂T ∂y n y, (4.15) i podstawiając do równania Galerkina otrzymujemy: λ ∂Te ∂x Ω e ∂N i e ∂x + λ ∂Te ∂y ∂N i e dxdy = ∂y (4.16) = Q̇ N e i dxdy + λ ∂Te ∂x Ω e dy − λ ∂Te ∂y dx N i e Γ e . 41
Uwzględniając warunki brzegowe otrzymujemy: λ ∂Te ∂x Ω e ∂N i e ∂x + λ ∂Te ∂y ∂N i e dxdy = ∂y (4.17) = Q̇ N e i dxdy + Q̇ 2N e i ds + α(T m − T 1 ) N e i ds, Ω e Γe 2 Γe 3 gdzie Q̇ ds = λ ∂T ∂x n xds + λ ∂T ∂y n yds. (4.18) Wprowadzając do powyższego równania (4.17) zależność (4.11) otrzymujemy: M λ ∂N i e ∂x T ∂N j e j e ∂x + λ ∂N i e ∂y T ∂N j e j e ∂y dxdy = Ω e = Q̇ N i e dxdy Ω e j=1 + Q̇ 2N i e ds Γ 2 e M M j=1 − α T j e N j e Γ 3 e j=1 N e i ds + αT m N e i ds. Γ 3 e (4.19) Powyższe równanie można sprowadzić do algebraicznego układu równań postaci: Ka = f. (4.20) W klasycznej metodzie elementów skończonych macierz K zwykle nazywa się macierzą sztywności, zaś w przypadku zagadnień przewodnictwa ciepła spotyka się również określenie: macierz przewodności. Wektor a reprezentuje rozwiązanie, zaś f - wektor obciążeń. Poszczególne elementy równania macierzowego (4.19) przedstawiają się następująco: K = K e c + K e Γ3 , f = f e e q + f Γ2 + f e Γ3 , (4.21) (4.22) gdzie: 42
- Page 1 and 2: Politechnika Poznańska Wydział Bu
- Page 3 and 4: 5 Optymalizacja położenia źróde
- Page 5 and 6: Streszczenie Praca poświęcona jes
- Page 7 and 8: zegowych równań całkowych, co zm
- Page 9 and 10: wyliczane są współczynniki wagow
- Page 11 and 12: W pracy analizowane są zagadnienia
- Page 13 and 14: włókien, oraz przypadek odwrotny.
- Page 15 and 16: 2 Optymalizacja 2.1 Wprowadzenie Dl
- Page 17 and 18: W powyższym wzorze f oznacza maksy
- Page 19 and 20: z punktów w przestrzeni poszukiwa
- Page 21 and 22: osobniki co oznacza, że w wyniku r
- Page 23 and 24: manipulowanie wartościami poszczeg
- Page 25 and 26: twórcę programu znającego charak
- Page 27 and 28: metody jest liczność populacji ty
- Page 29 and 30: 3 Metoda rozwiązań podstawowych 3
- Page 31 and 32: a ij = Bφ i x j , j = 1, … , N,
- Page 33 and 34: Tabela 3.1 Przykładowe zbiory funk
- Page 35 and 36: najnowszych opublikowanych badań z
- Page 37 and 38: do obliczeń współrzędnych punkt
- Page 39 and 40: gdzie L i B oznaczają operator lin
- Page 41: 4.3 Rozwiązanie dwuwymiarowego ust
- Page 45 and 46: 5 Optymalizacja położenia źróde
- Page 47 and 48: Na brzegu Γ zadana jest stała tem
- Page 49 and 50: Ponieważ istnieje tak duża liczba
- Page 51 and 52: 5.5 Wyniki numeryczne Do obliczeń
- Page 53 and 54: Jak wspomniano wcześniej, w algory
- Page 55 and 56: W drugim cyklu obliczeń populacja
- Page 57 and 58: Rysunek 5.12 Konfiguracja 4 element
- Page 59 and 60: Zestawienie uzyskanych wyników dla
- Page 61 and 62: Rysunek 5.20 Rozwiązanie optymalne
- Page 63 and 64: W związku z tym, że położenie p
- Page 65 and 66: gdzie ∂ ∂n oznacza pochodną w
- Page 67 and 68: Tabela 6.1 Zestawienie wyników obl
- Page 69 and 70: Wykresy ilustrujące rozwiązanie i
- Page 71 and 72: Występująca na brzegu osobliwoś
- Page 73 and 74: odsunięta jest o 2 od brzegu obsza
- Page 75 and 76: 7 Wyznaczanie efektywnego współcz
- Page 77 and 78: 7.3 Równanie rządzące i warunki
- Page 79 and 80: T i = c ij ln r j 2 j dla i = M, F
- Page 81 and 82: 7.7 Wyniki numeryczne 7.7.1 Komórk
- Page 83 and 84: Rysunek 7.5 Układ włókien dla os
- Page 85 and 86: włókien. Parametry materiału prz
- Page 87 and 88: Rysunek 7.8 Układ włókien dla os
- Page 89 and 90: przedstawiono w tabeli 7.18. Wykorz
- Page 91 and 92: Rysunek 7.11 Układ włókien dla o
M<br />
T e (x, y) = T e j ∙ N e j (x, y), (4.11)<br />
j=1<br />
gdzie M odpowiada liczbie węzłów przyjętego elementu e, T e j oznacza temperaturę<br />
w węźle j elementu e, N e j to funkcja kształtu elementu e. Stosując przedstawioną<br />
wcześniej metodę Galerkina, możemy zapisać:<br />
∂ ∂Te<br />
λ<br />
∂x ∂x + ∂ ∂Te<br />
λ<br />
∂y ∂y + Q̇ N e i (x, y)dxdy = 0.<br />
(4.12)<br />
Ω e<br />
Korzystając z twierdzenia Greena, otrzymujemy następującą zależność:<br />
∂ ∂Te<br />
λ<br />
∂x<br />
Ω e<br />
∂x N i e + ∂ ∂y<br />
λ<br />
∂Te<br />
∂y N i e dxdy =<br />
(4.13)<br />
= λ ∂Te<br />
∂x<br />
dy − λ<br />
∂Te<br />
∂y dx N i e<br />
Γ e .<br />
Przekształcając lewą stronę równania otrzymujemy:<br />
∂ ∂Te<br />
λ<br />
∂x<br />
Ω e<br />
∂x + ∂ ∂y<br />
λ<br />
∂Te<br />
∂y N i e dxdy =<br />
(4.14)<br />
= − λ ∂Te e<br />
e<br />
∂N i ∂Te ∂N i ∂Te ∂Te<br />
+ λ dxdy + λ dy − λ<br />
∂x ∂x ∂y ∂y ∂x ∂y dx N i e ,<br />
Ω e Γ e<br />
Biorąc pod uwagę, że<br />
Q̇ = λ ∂T<br />
∂x n x + λ ∂T<br />
∂y n y, (4.15)<br />
i podstawiając <strong>do</strong> równania Galerkina otrzymujemy:<br />
λ ∂Te<br />
∂x<br />
Ω e<br />
∂N i<br />
e<br />
∂x<br />
+ λ<br />
∂Te<br />
∂y<br />
∂N i<br />
e<br />
dxdy =<br />
∂y<br />
(4.16)<br />
= Q̇ N e i dxdy + λ ∂Te<br />
∂x<br />
Ω e<br />
dy − λ<br />
∂Te<br />
∂y dx N i e<br />
Γ e .<br />
41