Lõplike elementide meetod

Allikas: Vikipeedia

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[muuda | redigeeri lähteteksti]

Esmakordselt kasutas piirkonna diskreetimist saksa matemaatik Karl Schellbach 1851. aastal, seda Plateau probleemi lahendamisel. [1] 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.[2] 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.[3] 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. [4]

Meetodi põhimõte[muuda | redigeeri lähteteksti]

Lõplike elementide võrgustik
Lõplike elementide 2D-võrgustik
LEM 2D-lahend
LEMi 2D-lahend

Üldjoontes viiakse LEM läbi järgnevalt:[5]

  1. määratakse probleemi dimensionaalsus, vastavad võrrandid, nende kehtivus- ja määramispiirkonnad (defineerib elementide dimensionaalsuse);
  2. 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.
  3. 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);
  4. kõikide elementide panused summeeritakse, saades kogu probleemi kirjeldava globaalse võrrandite süsteemi;
  5. saadud võrrandite süsteem lahendatakse numbriliselt (vaja algväärtusi, nt temperatuurijaotus);
  6. 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[muuda | redigeeri lähteteksti]

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

(1):u''(x)=f(x), x\in(0,1),\qquad u(0)=u(1)=0,

kus f(x) on teadaolev funktsioon ja u(x) 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 v(x), mis rahuldab

(2):v(0)=v(1)=0,

saadakse

(3):\int_0^1 v(x)u''(x) dx=\int_0^1 v(x)f(x) dx.

Võrrandi vasaku poole ositi integreerimine annab

\int_0^1 v(x)u''(x) dx=v(x)u'(x)|_0^1-\int_0^1 v'(x)u'(x)dx=-\int_0^1 v'(x)u'(x)dx,

milles viimase võrduse saamiseks kasutati valemit (2). Defineerides funktsionaali \phi(u,\cdot) järgnevalt:

\phi(u,v)=\int_0^1 u'(x)v'(x) dx,

saab valemi (3) kirja panna kui

(4):\int_0^1 v(x)f(x) dx=-\phi(u,v).

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

Diskreetimine[muuda | redigeeri lähteteksti]

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

\mbox{ leia } u\in H^1(0,1)\mbox{ nii, et }\forall v\in H^1(0,1)\mbox{ korral }-\phi(u,v)=\int fv

asendamisele lõpliku-dimensioonilise probleemiga

(5):\mbox{ leia } u\in V\subset H^1(0,1)\mbox{ nii, et }\forall v\in V\mbox{ korral }-\phi(u,v)=\int fv

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

Võttes intervalli (0,1), jagatakse see n+1 osaks: 0=x_0<x_1<\ldots<x_n<x_{n+1}=1. Ruumi V määratletakse kui

(6):V=\{v:[0,1]\rightarrow\mathbb{R} \big| v\in C^0, v|_{[x_k,x_{k+1}]}\mbox{ on lineaarne }\forall k=0,\ldots,n\mbox{ ja } v(x_0)=v(x_{n+1})=0\},

milles C^0 tähistab pidevate funktsioonide hulka. Ülaltoodud definitsiooni kohaselt pole v diferentseeruv kogu määramispiirkonna ulatuses, st punktides \{x_k\}_{k=0}^{n+1} funktsiooni v tuletist ei leidu.

Baasi valik[muuda | redigeeri lähteteksti]

Baasifunktsioonid vk (sinised) ning lineaarkombinatsioon neist (punane)
Baasifunktsioonid vk (sinised) ning lineaarkombinatsioon neist (punane)

Ruumi V jaoks saab valida järgneva baasi:

v_k(x)=\begin{cases}
\frac{x-x_{k-1}}{x_k-x_{k-1}}, & \mbox{ kui } x\in[x_{k-1},x_k],\\
\frac{x_{k+1}-x}{x_{k+1}-x_k}, & \mbox{ kui } x\in[x_k,x_{k+1}],\qquad\forall k=1,\ldots,n\\
0 & \mbox{ muul juhul }
\end{cases}.

Sellise definitsiooni kohaselt on igas intervallis [x_{k-1},x_{k+1}] defineeritud baasifunktsioon v_k kolmnurkfunktsioon. Kehtib v_i(x_j)=\delta_{ij}, kus \delta_{mn} tähistab Kroneckeri sümbolit. Baas rahuldab järgmisi omadusi:

(7):\int_0^1 v_jv_k dx=0,\qquad\phi(v_j,v_k)=0,\qquad\mbox{ kui kehtib }|j-k|>1.

Maatriksi kujule viimine[muuda | redigeeri lähteteksti]

Kirjutades otsitava funktsiooni u(x) baasifunktsioonide v_k(x) kaudu

u(x)=\sum_{k=1}^nu_kv_k(x)

ja tähistades valemis v(x)=v_j(x), saab valemi (5) kirja panna kui

(8):-\sum_{k=1}^nu_k\phi(v_k,v_j)=\sum_{k=1}^n\int f(x)v_k(x) dx,\qquad\forall j=1,\ldots,n.

Defineerides vektorid \mathbf{u}=(u_1,\ldots,u_n)^T ja \mathbf{b}=(b_1,\ldots,b_n)^T, kus b_i=\int f(x)v_i(x) dx, ning maatriksi L=L_{ij}=\phi(v_i,v_j), saab valemi (8) avaldada kui

(9):-L\mathbf{u}=\mathbf{b}.

Järelikult on ääreväärtusprobleem taandunud \mathbf{u} leidmisele võrrandis (9). Omaduse (7) tõttu on maatriks L 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[muuda | redigeeri lähteteksti]

LEMi on võimalik rakendada ka kõrgemaid (osa)tuletisi sisaldavatele ääreväärtusprobleemidele, näiteks Neumanni ja segaprobleemile, üle suvalise piirkonna \Omega. Diferentsiaalvõrrandi järk määrab baasifunktsioonide pidevuse tingimused, piirkonna dimensionaalsus aga elementide dimensiooni, näiteks \Omega\subset\mathbb{R}^2 võib jaotada kolmnurkadeks, ruutudeks või muuks kumeraks hulknurgaks, \Omega\subset\mathbb{R}^3 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

\int_\Omega fv dV=-\int_\Omega \nabla u \cdot \nabla v dV = -\phi(u,v).

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[muuda | redigeeri lähteteksti]

Lõplike vahede meetod (LVM) on alternatiivne viis osatuletistega diferentsiaalvõrrandite (süsteemi) lahendamiseks. Põhilised erinevused LEMi ja LVMi vahel on[6]:

  • 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[muuda | redigeeri lähteteksti]

LEM on populaarne mehaanika valdkonda kuuluvate insener-tehniliste probleemide simulatsioonides. LEM on integreeritud aeronautika-, ehitus-,[7] laeva-,[8] auto- ja isegi mööblitööstuse[9] 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.[10] 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[muuda | redigeeri lähteteksti]

Viited[muuda | redigeeri lähteteksti]

  1. Kurrer, Karl-Eugen (2008). The History of the Theory of Structures: From Arch Analysis to Computational Mechanics. Ernst & Sohn, 655. ISBN 978-3433018385. 
  2. Argyris, J.H.; Kelsey, S. (1954). "Energy Theorems and Structural Analysis". Aircraft Engineering 26: 347—356. 
  3. Gupta, K.K.; Meek, J.L. (1996). "A Brief History of the Finite Element Method". International Journal for Numerical Methods in Engineering 39: 3761—3774. 
  4. "Geschichte der FEM". www.cadfem.de. Vaadatud 21. sept 2013. (saksa k)
  5. Kasemägi, Heiki (2009). Kompuuterfüüsika II. 
  6. "FEM-FDM" (PDF). www.delcam-ural.ru. Vaadatud 24. sept 2013. lk 3 (ingl k)
  7. Mahendran, M. (2007). "Applications of Finite Element Analysis in Structural Engineering" 38—46. eprints.qut.edu.au (ingl k)
  8. Piala, Peter; Kalina, Tomas (2010). "Finite Element Method in Ship's Industry". 
  9. Koc; Kizilkaya; Erdinler; Korkut (2011). "The use of finite element method in the furniture industry". African Journal of Business Management 5(3): 855—865. doi:10.5897/AJBM10.551. 
  10. Petrin, H.; Wiesler, B. (20.—22. sept 1999). "Application of FEM — Tools in the Engine Development Process". web.mscsoftware.com lk 3 (ingl k)

Kirjandus[muuda | redigeeri lähteteksti]

Eesti keeles[muuda | redigeeri lähteteksti]

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

Võõrkeeltes[muuda | redigeeri lähteteksti]

  • 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[muuda | redigeeri lähteteksti]