Next: Metoda konačnih elemenata
Up: Rubni problem
Previous: Metoda konačnih diferencija
  Sadržaj
  Indeks
Ritzova metoda
Iz varijacijskog principa slijedi da se umjesto rješavanja jednadžbe,
rješenje rubnog problema može dobiti rješavanjem problema
minimizacije odgovarajućeg funkcionala.
Ima raznih metoda koje koriste varijacijsku formulaciju. Jedna od njih
je Ritzova. Pogledajmo kako ona funkcionira kad se radi o rubnom
problemu
gdje je za svaki
Izaberemo linearno nezavisnih funkcija
koje
zadovoljavaju Dirichletov homogeni rubni uvjet. Dakle
Rješenje se pretpostavi u obliku
i neodređeni koeficijenti se odrede iz uvjeta da minimizira
pripadni funkcional energije (v. teorem 18)
Tako
dobiveni leži u vektorskom prostoru razapetom s funkcijama
tj. u vektorskom prostoru svih linearnih
kombinacija funkcija Rješenje ne mora biti
linearna kombinacija tih funkcija, pa u tom slučaju nije točno
već samo približno rješenje. No, grešku koju tom pretpostavkom
činimo, možemo umanjiti uzimanjem većeg
Prvi problem s kojim se susrećemo kod Ritzove metode je određivanje
funkcija koje se često zovu koordinatne funkcije. One
se obično biraju u skladu s problemom koji se rješava. Na pr. ako
iz fizikalnih razloga očekujemo periodičko rješenje, onda ćemo
takvima pretpostaviti i funkcije Inače se za koordinatne
funkcije mogu koristiti polinomi.
Nakon što smo izabrali funkcije
pretpostavljeno rješenje uvrstimo u funkcional
je derivabilna
funkcija od varijabli
pa jednadžbe
predstavljaju nužan uvjet za ekstrem funkcije u točki
Ovo je sustav od linearnih algebarskih
jednadžbi s nepoznanica. Ako stavimo
dobivamo sustav jednadžbi
On se može matrično zapisati
|
(3.27) |
gdje je
Matrica koja se inače zove matrica krutosti, je očito
simetrična, jer je
Ona je i pozitivno definitna. Doista,
Budući da je za svaki
ovaj integral se može poništavati samo tako da bude
odnosno
Odatle
a kako je
slijedi odnosno
Zbog linearne nezavisnosti koordinatnih
funkcija slijedi da je za svaki tj.
Dakle
za
pa je matrica pozitivno definitna.
Kao što smo vidjeli, takva
matrica se može dijagonalizirati i ima pozitivne vlastite
vrijednosti. To znači da je regularna. Zaista, neka su njezine
vlastite vrijednosti
Tada je
Odatle slijedi
da jednadžba (3.28) ima jedno i samo jedno rješenje.
Nedostatak ove metode je u tome što je matrica puna matrica, tj.
općenito je svaki njezin element različit od nule.
Primjer 3.20
Riješimo Ritzovom metodom sljedeći rubni problem
Rješenje. Za koordinatne funkcije uzmimo polinome. Budući da na lijevom
rubu imamo homogen Dirichletov uvjet, polinomi se moraju
poništavati u nuli. Neka su, dakle, koordinatne funkcije
Rješenje pretpostavljamo u obliku
Kad s
uđemo u funkcional, i njegove
derivacije po koeficijentima
izjednačimo s
nulom, dobivamo sustav jednadžbi
Odavde dobijemo za koeficijente
pa je rješenje
Točno rješenje je
Greška računata u točkama
iznosi
Next: Metoda konačnih elemenata
Up: Rubni problem
Previous: Metoda konačnih diferencija
  Sadržaj
  Indeks
2001-10-26