MES - warunki brzegowe

rafal snipped-for-privacy@op.pl pisze:

Witam,

Dokładnie jak opisałeś. W macierzy sztywności w wierszach, które odpowiadają węzłom na brzegu wpisujesz na głównej przekątnej 1 a poza główną przekątną 0. Natomiast w wektorze prawych stron wartość temperatury na brzegu. Uwaga: Sprawdz w jakich jednostkach masz współczynnik przewodności termicznej lambda. Bo jezeli w [W/mK], a nie w [W/mC] to wtedy musisz podac temperature w skali absolutnej, czyli T=293K.

Pozdrawiam Szymon

Dzieki. Mam jeszcze kilka pytan.

1.W macierzy sztywno[ci i wektorze obciazeD mam caBki krzywoliniowe po brzegu (przy wyprowadzaniu wzorów na elementy macierzowe mam uwzglednione warunki brzegowe). Czy wpisywanie warto[ci na przekatnej w takim wypadku nadal stosowac ?

1. Mam na kwadracie 25 wezBów (16 elementów kwadratowych). Po rozwiazaniu równania macierzowego mam 25 wspóBczynników. Mam równie| 4 funkcje ksztaBtu dla elementu w wspóBrzdnych referencyjnych. Jak teraz mo|na wyko|ystac znalezione wspóBczynniki do zbudowania szukanej funkcji rozkBadu temperatury ?

3.Przy transformacjach zachowujących kształt (kwadrat=kwadrat). Jakobian wynosi J=S/4 gdzie S to pole kwadratu w wspóBrzdnych globalnych. Równanie macierzowe: S*T=F mam podane wzory w wspóBrzdnych globalnych Hij=caBka_po_obszarze(kx*pochodna_po_x(f_ksztaBtu_i) +ky*pochodna_po_y(f_ksztaBtu_j)) + +caBka_po_brzegu(alfa*f_ksztaBtu_i*f_ksztaBtu_j) Fi=-caBka_po_brzegu((alfa*T0-q)*f_ksztaBtu_i)- caBka_po_obszarze(Q*f_ksztaBtu_i)

Czy macierze lokalne w wspóBrzdnych referencyjnych beda tak wygldac ?: Hlok_ij=(tu jakobian sie skróci)caBka_po_obszarze(to_samo)

+J*caBka_po_brzegu(to_samo) Flok_i=-J*caBka_po_brzegu(to_samo)-J*caBka_po_obszarze(to_samo)