next up previous contents
Next: Metoda konačnih elemenata Up: Rubni problem Previous: Metoda konačnih diferencija   Contents


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

$\displaystyle (p(x)\,u'(x))' +f(x) = 0,\hspace{1cm}u(0) = 0,\;\;u'(\ell{}) = 0,$

gdje je $ p(x)>0$ za svaki $ x.$

Izaberemo $ n$ linearno nezavisnih funkcija $ v_i, i=1,2,\ldots, n$ koje zadovoljavaju Dirichletov homogeni rubni uvjet. Rješenje se pretpostavi u obliku

$\displaystyle u_n(x) = \sum_{i=1}^n c_i\,v_i(x),$

i neodređeni koeficijenti se odrede iz uvjeta da $ u_n$ minimizira funkcional

$\displaystyle F(u) = \frac{1}{2}\,\int_0^{\ell}\,p\,{u'}^2\,dx - \int_0^{\ell}\,f\,u\,dx.$

Tako dobiveni $ u_n$ leži u vektorskom prostoru razapetom s funkcijama $ v_1,v_2,\ldots,v_n,$ tj. u vektorskom prostoru svih linearnih kombinacija funkcija $ v_1,v_2,\ldots,v_n.$ Rješenje ne mora biti linearna kombinacija tih funkcija, pa u tom slučaju $ u_n$ nije točno već samo približno rješenje. No, grešku koju tom pretpostavkom činimo, možemo umanjiti uzimanjem većeg $ n.$

Prvi problem s kojim se susrećemo kod Ritzove metode je određivanje funkcija $ v_i,$ 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 $ v_i.$ Inače se za koordinatne funkcije mogu koristiti polinomi.

Nakon što smo izabrali funkcije $ v_i,$ pretpostavljeno rješenje uvrstimo u funkcional

$\displaystyle F(u_n) = \frac{1}{2}\,\int_0^{\ell}\,p\,\left(\sum_{i=1}^n
c_i\,v'_i\right)^2\,dx - \int_0^{\ell}\,f\,\sum_{i=1}^n c_i\,v_i\,dx$

$\displaystyle = \frac{1}{2}\,\sum_{i=1}^n\sum_{j=1}^n
c_i\,c_j\,\int_0^{\ell}\,p\,{v'}_i\, v'_j\,dx - \sum_{i=1}^n c_i\,\int_0^{\ell}\,
f\,v_i\,dx$

$\displaystyle = \Phi(c_1,c_2,\ldots,c_n).$

$ \Phi$ je derivabilna funkcija od $ n$ varijabli $ c_1,c_2,\ldots,c_n,$ pa jednadžbe

$\displaystyle \frac{\partial \Phi}{\partial c_i} = \sum_{j=1}^n c_j\,\int_0^{\ell}\,
p\,v'_i\, v'_j\,dx - \int_0^{\ell}\,f\,v_i\,dx = 0,$

za $ i=1,2,\ldots,n,$ predstavljaju nužan uvjet za ekstrem funkcije $ \Phi$ u točki $ (c_1,c_2,\ldots,c_n).$ Ovo je sustav od $ n$ linearnih algebarskih jednadžbi s $ n$ nepoznanica. Ako stavimo

$\displaystyle K_{i\,j} = \int_0^{\ell}\,p\,v'_i\,v'_j\,dx,\hspace{1cm}b_i = \int_0^{\ell}\,
f\,v_i\,dx,$

dobivamo sustav jednadžbi

$\displaystyle \sum_{j=1}^n K_{i\,j}\,c_j = b_i,\qquad i=1,2,\ldots,n.$

On se može matrično zapisati

$\displaystyle K\,\boldsymbol{c} = \boldsymbol{b},$ (3.23)

gdje je $ K = [K_{ij}],\; \boldsymbol{c}=[c_j],\;
\boldsymbol{b}=[b_i].$

Matrica $ K,$ koja se inače zove matrica krutosti, je očito simetrična, jer je

$\displaystyle K_{i\,j} = \int_0^{\ell}\,p\,v'_i\,v'_j\,dx = \int_0^{\ell}\,p\,v'_j\,v'_i\,dx = K_{ji}.$

Ona je i pozitivno definitna. Doista,

$\displaystyle K\,\boldsymbol{c}\cdot{}\boldsymbol{c} = \left[\sum_{j=1}^n
K_{ij}\,c_j\right]\cdot{}[c_i] = \sum_{i=1}^n \sum_{j=1}^n
K_{ij}\,c_j\,c_i =$

$\displaystyle \sum_{i=1}^n \sum_{j=1}^n \int_0^{\ell}\,
p\,v'_i\,v'_j\,dx\,c_j\...
...i\,c_j\,v'_j\,dx = \int_0^{\ell}\,p\,\left(\sum_{i=1}^n
c_i\,v'_i\right)^2\,dx.$

Budući da je $ p(x)>0,$ za svaki $ x,$ ovaj integral se može poništavati samo tako da bude

$\displaystyle \left(\sum_{i=1}^n c_i\,v'_i\right)^2 = 0,$

odnosno

$\displaystyle \sum_{i=1}^n c_i\,v'_i = 0,$

a to se može dogoditi jedino u slučaju da je $ c_i=0$ za svaki $ i,$ tj. $ \boldsymbol{c}=\boldsymbol{0}.$ Dakle $ K\boldsymbol{c}\cdot\boldsymbol{c}>0,$ za $ \boldsymbol{c}\neq\boldsymbol{0},$ pa je matrica $ K$ 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 $ \lambda{}_1,\lambda{}_2,\ldots,\lambda{}_n.$ Tada je

$\displaystyle \det{}K = \det{}(X^{-1}[\lambda{}_i\,\delta_{ij}]\,X) =
\det{}X^{...
...i\,\delta_{ij}]\,\det{}X =
\lambda{}_1\,\lambda{}_2\,\cdots{}\,\lambda{}_n > 0.$

Odatle slijedi da jednadžba (3.23) ima jedno i samo jedno rješenje.

Nedostak ove metode je u tome što je matrica $ K$ puna matrica, tj. općenito je svaki njezin element različit od nule.


next up previous contents
Next: Metoda konačnih elemenata Up: Rubni problem Previous: Metoda konačnih diferencija   Contents
Salih Suljagic
1999-01-27