Lõplike elementide meetod

Lõplike elementide meetod ehk LEM (ingl k finite element method, FEM) on numbriline arvutusmeetod insener-tehniliste ja füüsikaliste probleemide lahendamiseks. Meetod hõlmab üldiste (osatuletistega) diferentsiaalvõrrandite ja/või nende süsteemi aproksimeerimist ehk lähendamist üle mingi (pideva) piirkonna (nt tasand), mis jaotatakse väiksemateks, lõplikeks elementideks (nn piirkonna diskreetimine). Numbriline lahend kogu piirkonna jaoks saadakse lõplike elementide summeerimisega mingi eeskirja järgi.

Ajalugu 

Esmakordselt kasutas piirkonna diskreetimist saksa matemaatik Karl Schellbach 1851. aastal, seda Plateau probleemi lahendamisel.  Lord Rayleigh', Walther Ritzi ja Boris Galjorkini hilisem arendustöö variatsiooniarvutuse ja kaalutud hälvete meetodites andis samuti märkimisväärse panuse LEMi teooriasse.
Hiljem on diskreetimist kasutanud Richard Courant, kes 1943. aastal uuris Saint-Venant' väändeprobleemi, kasutades ka Rayleigh'-Ritzi meetodit, kuigi tema töö jäi suurema tähelepanuta. Kaasaegse LEMi formuleerimisele aitas kaasa John Argyrise 1954. aastal avaldatud artikkel struktuurianalüüsi ja energiateoreemi teemadel. Argyrisest sõltumatult panustasid tänapäevase LEMi väljatöötamisse Ray William Clough (1960) ja M. J. Turner (1956), kellest viimane võttis kasutusele termini lõplik element.
Seoses aeronautikatööstuse ja arvutusjõudluse arenguga 1960. aastatel arendas NASA välja esimese LEMi tarkvara NASTRAN struktuurianalüüsi arvutuste jaoks. Edasistel aastatel kujunes LEM tööstusstandardiks, laienedes esialgsest struktuurianalüüsist soojusjuhtivuse, hüdromehaanika ja elektrodünaamika simulatsioonideni. Personaalarvutite tulekuga 1970. aastate lõpus ilmusid turule esimesed LEMi kommertstarkvarapaketid.

Meetodi põhimõte 

Üldjoontes viiakse LEM läbi järgnevalt:
 määratakse probleemi dimensionaalsus, vastavad võrrandid, nende kehtivus- ja määramispiirkonnad (defineerib elementide dimensionaalsuse);
 piirkond jaotatakse väiksemateks osadeks e lõplikeks elementideks (kolmnurk, ristkülik jne) e diskreeditakse, st koostatakse lõplike elementide võre (ingl k mesh), milles iga element on seotud teise elemendiga võre sõlmpunkti kaudu (vt parempoolset joonist); piirkonna lihtsamateks osadeks jaotamisel on mitu eelist:
 keeruka geomeetria täpne kirjeldus;
 erinevate materjaliomaduste arvestamine;
 lokaalsete efektide arvestamine.
 iga element, täpsemini võre sõlmpunkt, seostatakse
 vabadusastmete arvuga (määrab elemendi edasise käitumise);
 materjali, üldiselt välja omadustega (termodünaamilistes simulatsioonides nt soojusmahtuvusega);
 asjakohaste (füüsikaliste) harilike diferentsiaal- või algebraliste võrranditega (eeldatakse, et sõltumatud muutujad võre sõlmpunktides on teada või siis need aproksimeeritakse mingi polünoomi abil);
 kõikide elementide panused summeeritakse, saades kogu probleemi kirjeldava globaalse võrrandite süsteemi;
 saadud võrrandite süsteem lahendatakse numbriliselt (vaja algväärtusi, nt temperatuurijaotus);
 vajalikud füüsikalised suurused arvutatakse iga elemendi jaoks.
Ülaltoodud punktidest järeldub, et LEM sobib mittelineaarsete efektide uurimiseks. Samuti nähtub punktidest, et vaadeldav pidev süsteem lõpmatu arvu tundmatutega lähendatakse diskreetse süsteemi ja lõpliku arvu tundmatutega. Kui teatud piirkonnas vajatakse suuremat täpsust, siis seal saab konstrueerida tihedama elementide võrgustiku. Tegemist oleks optimaalse lahendusega, kuna kogu piirkonda pole mõtet väga tihedalt diskreetida: nõudmised arvuti muutmälu, salvestusseadmete ja protsessorijõudluse osas kasvaksid viimasel juhul drastiliselt.

Probleeminäide 

LEMi tehniliseks käsitluseks sobib ühedimensiooniline homogeenne Dirichlet' probleem

(1):

kus  on teadaolev funktsioon ja  otsitav lahend. Esmalt tuleb formuleerida ääreväärtusprobleemi nõrk formulatsioon e nõrk vorm, seejärel see diskreetida.

Korrutades esimest võrrandit valemis (1) skalaarselt suvalise funktsiooniga , mis rahuldab

(2):

saadakse

(3):

Võrrandi vasaku poole ositi integreerimine annab

milles viimase võrduse saamiseks kasutati valemit (2). Defineerides funktsionaali  järgnevalt:

saab valemi (3) kirja panna kui

(4):

Tegemist on nõrga vormiga probleemi (1) jaoks selles mõttes, et see seab nõrgemad tingimused lahendile  ja proovifunktsioonile . Kui ääreväärtusprobleem oleks defineeritud üle mingi suvalise piirkonna , siis vastavad testfunktsioonid moodustaksid Sobolevi ruumi . Testfunktsioonid rahuldavad , kus  on piirkonna  tükati sile raja. On võimalik näidata, et . Probleeminäite testfunktsioonide ruumiks on .

Diskreetimine 

Nõrga vormi (4) diskreetimine taandub lõpmatu-dimensioonilise lineaarse probleemi

asendamisele lõpliku-dimensioonilise probleemiga

(5):

kus  on lõpliku dimensiooniga tükiti polünoomide ruum.

Võttes intervalli , jagatakse see  osaks: . Ruumi  määratletakse kui

(6):

milles  tähistab pidevate funktsioonide hulka. Ülaltoodud definitsiooni kohaselt pole  diferentseeruv kogu määramispiirkonna ulatuses, st punktides  funktsiooni  tuletist ei leidu.

Baasi valik 

Ruumi  jaoks saab valida järgneva baasi:

Sellise definitsiooni kohaselt on igas intervallis  defineeritud baasifunktsioon  kolmnurkfunktsioon. Kehtib , kus  tähistab Kroneckeri sümbolit. Baas rahuldab järgmisi omadusi:

(7):

Maatriksi kujule viimine 

Kirjutades otsitava funktsiooni  baasifunktsioonide  kaudu

ja tähistades valemis , saab valemi (5) kirja panna kui

(8):

Defineerides vektorid  ja , kus , ning maatriksi , saab valemi (8) avaldada kui

(9):

Järelikult on ääreväärtusprobleem taandunud  leidmisele võrrandis (9). Omaduse (7) tõttu on maatriks  sümmeetriline, positiivselt määratud ja hõre, st enamus tema elementidest on nullid. Seda tüüpi võrrandite jaoks on välja töötatud mitmesuguseid algoritme, näiteks konjugeeritud gradiendi meetod ning LU- ja Cholesky dekompositsioon, mis arvutite peal efektiivselt realiseeruvad.

Üldistused 

LEMi on võimalik rakendada ka kõrgemaid (osa)tuletisi sisaldavatele ääreväärtusprobleemidele, näiteks Neumanni ja segaprobleemile, üle suvalise piirkonna . Diferentsiaalvõrrandi järk määrab baasifunktsioonide pidevuse tingimused, piirkonna dimensionaalsus aga elementide dimensiooni, näiteks  võib jaotada kolmnurkadeks, ruutudeks või muuks kumeraks hulknurgaks,  tetraeedriteks jne (ühedimensioonilises probleeminäites on piirkond sirglõikudeks jaotatud). Nõrk vorm formuleeritakse Greeni identsusi kasutades, nt homogeense Dirichlet' probleemi nõrgaks vormiks on

Baasifunktsioonid e interpolatsioonifunktsioonid on üldjuhul tükati polünoomid, nt Hermite'i ja Lagrange'i polünoomid, kusjuures nende väärtused peavad pidevuse tagamiseks erinevate elementide kokkupuuteäärel olema võrdsed. Nõrga vormi koostamist ja selle lähendamist lõpliku baasiga nimetatakse Galjorkini meetodiks. Ülesanne taandub lõpuks numbrilisele maatriksivõrrandi lahendamisele või Runge-Kutta meetodi rakendamisele.

Võrdlus lõplike vahede meetodiga 

Lõplike vahede meetod (LVM) on alternatiivne viis osatuletistega diferentsiaalvõrrandite (süsteemi) lahendamiseks. Põhilised erinevused LEMi ja LVMi vahel on:
 LEM suudab keeruka geomeetria ja ääreväärtusprobleemi lahendada, samas kui LVMi lahendamise piirkond on ristkülikukujuline (või ristkülikutest koosnev);
 LVMi on lihtsam implementeerida;
 LVM on teatud probleemide (nt Poissoni võrrandi) puhul LEMi erijuht;
 LEMi lähendus on üldjuhul täpsem vastavast LVMi omast (kuigi see sõltub üsna palju probleemi püstitusest; LVMiga sama täpsuse saavutamiseks vajab LEM palju vähem võre sõlmpunkte, järelikult implementeerimisel ka vähem mälu andmete töötlemiseks).
LEMi kasutatakse üldiselt struktuurianalüüsis, samas kui arvutuslikus hüdrodünaamikas (nt õhuvoolude simulatsioon lennukitiiva ümber) kasutatakse rohkem LVMi või lõpliku ruumala meetodit.

Rakendused 

LEM on populaarne mehaanika valdkonda kuuluvate insener-tehniliste probleemide simulatsioonides. LEM on integreeritud aeronautika-, ehitus-, laeva-, auto- ja isegi mööblitööstuse disainimis- ning tootmisprotsessi. Paljud LEMi tarkvarapaketid koosnevad komponentidest, millest igaüks käsitleb erinevaid füüsikavaldkondi, nt termodünaamika, elektromagnetism, hüdrodünaamika ja mehaanika. LEM aitab visualiseerida materjalide pingeid, nende jaotust ja konstruktsiooni vastupidavust, samuti on võimalik minimeerida erinevaid karakteristikuid, näiteks kaalu, materjali- ja tootmiskulu. LEMi tarkvara võimaldab simulatsiooni parametriseerida (nt konstruktsiooni pingete simulatsioon, kus parameetriks on erineva raskusega koormised), arvestada ajalist sõltuvust ja saada tulemus soovitud täpsuse piires.

LEM on märkimisväärselt parandanud tootearenduse standardeid ja metodoloogiat, seda eriti tööstusvaldkonnas. Näiteks on LEMi kasutuselevõtuga kontseptsioonist tootmisliinini jõudmise aeg oluliselt vähenenud. Seda tänu prototüüpide tarkvaralise testimise võimalusele, mis elimineerib mitteoptimaalsed või -töötavad prototüübid ning nendega seotud kulud, näiteks materjal ja aeg, tootearenduse algetappides.

Lisaks tööstuslikule rakendusele uuritakse LEMi abil füüsikalisi olukordi, mille tingimustes on raske eksperimenti läbi viia või seda seletav mudel vajaks kontrollimist.

Vaata ka 
 Galjorkini meetod
 Kompuuterfüüsika
 Lõplike vahede meetod
 Numbrilised meetodid
 Rayleigh'-Ritzi meetod
 Variatsiooniarvutus

Viited

Kirjandus

Eesti keeles 
 Lahe, Andres. Lõplike elementide meetod. Tallinn, 2008. 
 Kirs, Jüri. Sissejuhatus lõplike elementide meetodisse (I,II,III osa). Tallinn, 2001.
 Kasemägi, Heiki. Kompuuterfüüsika II. Tartu, 2009.

Võõrkeeltes 

 Ottosen, N.S.; Petersson, H. Introduction to the Finite Element Method. Prentice Hall, 1992. ISBN 978-0134738772
 Zienkiewicz, O.C; Taylor, R.L.; Zhu, J.Z. The Finite Element Method: Its Basis and Fundamentals (6th ed). Butterworth-Heinemann, 2006. ISBN 978-0750663205
 Segerlind, L.J. Applied Finite Element Analysis (2nd ed). John Wiley & Sons, 1984. ISBN 978-0-471-80662-2
 Lewis, R.W.; Nithiarasu, P.; Seetharamu, K. Fundamentals of the Finite Element Method for Heat and Fluid Flow. John Wiley & Sons, 2004. ISBN 978-0470847893

Välislingid 
 Scholarpedia – Finite element method