Hoppa till innehållet

Finita elementmetoden

Från Wikipedia
(Omdirigerad från FEM)
Ett exempel på ett beräkningsnät (i dagligt tal kallat en mesh) som kan användas för analys med finita elementmetoden. Man har här försökt modellera ett rör i ett magnetfält. De olika färgerna illustrerar att materialen är olika.
Beräkningsresultatet från finita elementmetoden visualiserat. Här har FEM använts för att beräkna magnetisk flödestäthet.

Finita Elementmetoden (FEM) är en numerisk metod för att lösa fysikaliska problem som beskrivs med partiella differentialekvationer med hjälp av datorer.

Teorierna för finita elementmetoden utvecklades redan i början av 1900-talet,[1][2] men det är först med tillgången till moderna datorer med stor beräkningskapacitet som metoden blivit praktiskt användbar. Sedan 1980-talet har FEM fått mycket bred användning inom vetenskap och teknologi eftersom den tillåter att fysikaliska problem som värmeledning, elektriska fält, hållfasthet med mera kan analyseras i detalj, även för komplexa geometrier och för anisotropa och olinjära materialegenskaper.

Inom till exempel mekanisk konstruktion används FEM bland annat för hållfasthetsanalys. FEM är integrerat i moderna CAD-system vilket tillåter konstruktörer att snabbt och realistiskt kontrollera hållfastheten på detaljer som ännu inte existerar annat än som datormodeller. Vid design av elektriska maskiner kan FEM användas för att i förväg i datormodeller beräkna elektriska fält och uppskatta förluster redan innan någon fysisk maskin har byggts.[3][4]

Vid lösandet av partiella differentialekvationer är den huvudsakliga utmaningen att skapa en ekvation som approximerar den partiella differentialekvationen men som är numeriskt stabil, vilket betyder att små fel i indata inte förstoras så att resultatet blir oanvändbart. Det finns flera sätt att göra detta och finita elementmetoden är bra för att lösa partiella differentialekvationer på problem med komplex geometri, till exempel mekaniska konstruktioner, exempelvis bilar, eller när vissa delar av området kräver större noggrannhet än andra. Vid väderprognoser är det viktigare att ha korrekt prognos över land än till havs, vilket kan uppnås genom att ha ett mer detaljerat beräkningsnät i områdena där särskilt hög noggrannhet önskas.

Metoden i korthet

[redigera | redigera wikitext]

Finita elementmetoden används för att lösa partiella differentialekvationer approximativt. Första steget i FE-metoden är att skriva om differentialekvationen i variationsform. Sedan delas problemdomänen upp i en stor mängd små delområden med nodpunkter längs områdesgränserna. Till varje nodpunkt definieras en funktion som är kontinuerlig över hela problemdomänen men nollskild endast i de delområden som omger noden. Ekvationens approximativa lösning beskrivs som en linjärkombination av alla

Funktionerna kallas för formfunktioner eller nodfunktioner. Viktningen kan nu bestämmas genom insättning i variationsproblemet.[5]

Finita elementmetoden har sitt ursprung i behoven att lösa komplexa elasticitets- och strukturanalysproblem inom flyg- och strukturmekanik. Dess utveckling började med arbeten av A. Hrennikoff (1941) och R. Courant (1942) . Fastän tillvägagångssättet var olika hade de en viktig sak gemensam: uppdelandet av ett kontinuerligt område till en uppsättning diskreta delområden. Hrennikoffs arbete delade upp området genom att använda galleranalogi medan Courant använde finita triangulära delområden för lösning av andra ordningens differentialekvationer som kommer från problemet med vridning av en cylinder.

Metoden fick en ordentlig matematisk grund 1973 i och med publiceringen av Strang och Fixs An Analysis of The Finite Element Method, och har sedan dess generaliserats till numerisk modellering av fysiska problem inom många ingenjörsområden, exempelvis elektromagnetism och strömningsmekanik.

Implicita och explicita finita elementmetoder

[redigera | redigera wikitext]

Det finns två principiellt olika lösningsmetoder för finita element: Implicit och explicit.

Den implicita metoden lämpar sig bäst för problem med små olinjäriteter och stora krav på noggrannhet. Lösningen är ovillkorligt stabil och steglängden påverkar endast lösningens noggrannhet. Denna metod har större krav på RAM men ger å andra sidan noggrannare resultat. Exempel på analyser som lämpar sig för implicit lösning är statiska spänningsanalyser.

Den explicita metoden lämpar sig bättre för problem med stora olinjäriteter (materialmodeller, stora deformationer, kontaktytor etcetera). Lösningens stabilitet är avhängig steglängden och ger felaktiga resultat om steglängden blir för stor. Fördelen med denna metod är att man lättare får konvergens vid stora olinjäriteter, enklare att anpassa för att fullt ut utnyttja datorer med många processorer, och är inte lika minneskrävande. Exempel på analyser som använder explicit lösare är krock-, kontakt- och explosionsanalyser.

Antag att vi vill lösa den tidsoberoende värmeledningsekvationen med homogen värmeutveckling i domänen och konstant temperatur vid ränderna. Detta problem beskrivs av Poissons ekvation och kan enkelt lösas med finita elementmetoden. Lösningen kan generaliseras till andra differentialekvationer och högre dimensioner.[5]

Problemets ekvation kan skrivas

med randvillkoren

och skrivs först om i variationsform. Detta görs genom att ekvationen multipliceras med en godtagbar funktion . En godtagbar funktion är en funktion som uppfyller problemets randvillkor . Därefter integreras båda sidor. Efter partiell integration har vi

Domänen delas upp i en mesh med 6 noder. Nodernas läge väljs till

Meshen och basfunktionen .

Till varje intern nod knyts en så kallad basfunktion (ett så kallat finit element)

Nu antas att lösningen till problemet kan skrivas

Insättning i variationsformuleringen ger

Denna ekvation gäller för alla godtagbara funktioner . I synnerhet gäller de för alla basfunktioner . Successiv insättning av dessa i ovanstående ekvation ger det algebraiska ekvationssystemet

Ekvationssystemet kan skrivas i matrisform som

där matrisen K har komponenterna

och

FE-metoden har sitt ursprung inom strukturanalys och därför kallas K traditionellt för systemets styvhetsmatris. Vektorn f kallas för systemets lastvektor.

Integralerna kan i detta fall lösas analytiskt

Finita elementlösningen tillsammans med analytisk lösning.

Lösningen till är

Den analytiska lösningen är

I detta fall är FEM-lösningen exakt i nodpunkterna.

Exempel på programvara för FE-analys

[redigera | redigera wikitext]

Nuförtiden ingår FEM som en naturlig komponent i alla större CAD-system. Dessa paket är dock oftast begränsade till enklare analyser. För allmänna FEM-beräkningar finns ett antal omfattande programpaket. Nedan följer en inkomplett lista[3]

Proprietära

[redigera | redigera wikitext]
  • MSC/NASTRAN - Omfattande paket för linjär och olinjär strukturanalys. Härstammar från NASA och flygindustrin. Vida använt med anor från FEM:s barndom.
  • HyperWorks - Heltäckande simuleringsvektyg för linjära- olinjära- och strömningsproblem samt stelkroppsdynamik. Specialiserat på optimering av strukturer.
  • LS-DYNA - Omfattande paket för olinjär explicit struktur- och strömningberäkningar. Härstammar från atombombsimuleringar. Vida använt för krocksimulering i bilindustrin.
  • ABAQUS - Omfattande paket för olinjär strukturanalys.
  • ANSYS - Omfattande paket. Har sitt ursprung inom kärnkraftindustrin.
  • COMSOL - Multiphysics simulering.
  • ADINA - Specialiserat på olinjär analys

Öppen källkod

[redigera | redigera wikitext]
  • Elmer - Multiphysics simulering.
  • Calculix - Allmän olinjär strukturanalys.
  • Code_Aster - Omfattande franskt program med rötter i kärnkraftsindustrin
  1. ^ K. Wieghardt (1906). ”Über einen Grenzübergans der Elastizitätslehre und seine Anwendung auf die Statik hochgradig statisch unbestimmter Fachwerke”. Verhandlungen des Vereins z. Beförderung des Gewerbefleisses. Abhandlungen 85: sid. 139-176. 
  2. ^ R. Courant (1943). ”Varitional Methods for the Solution of Problems of Equilibrium and Vibrations”. Bulletin of the American Mathematical Society 49: sid. 1-23. 
  3. ^ [a b] Sunnersjö, Staffan (1999). FEM i praktiken: en introduktion till finita elementmetodens praktiska tillämpning (2., [granskade och kompletterade] utg.). Stockholm: Sveriges verkstadsindustrier. Libris 8374937. ISBN 91-7548-541-9 
  4. ^ Zienkiewicz, O. C.; Taylor Robert Leroy (2005) (på engelska). Finite Element Method Set [Elektronisk resurs]. Butterworth-Heinemann. Libris 11954543 
  5. ^ [a b] Logg, Anders.; Mardal Kent-Andre., Wells Garth. (2012) (på engelska). Automated Solution of Differential Equations by the Finite Element Method [Elektronisk resurs : The FEniCS Book]. Lecture Notes in Computational Science and Engineering, 1439-7358 ; 84. Berlin, Heidelberg: Springer Berlin Heidelberg. Libris 12733043. ISBN 9783642230998. https://launchpadlibrarian.net/83776282/fenics-book-2011-10-27-final.pdf 

Externa länkar

[redigera | redigera wikitext]

Rensa denna sidas cache