Table Of ContentMarcus Wagner
Lineare und
nichtlineare FEM
Eine Einführung mit Anwendungen
in der Umformsimulation mit LS-DYNA®
Lineare und nichtlineare FEM
Marcus Wagner
Lineare und nichtlineare FEM
Eine Einführung mit Anwendungen in der
Umformsimulation mit LS-DYNA®
MarcusWagner
FakultätMaschinenbau
OTHRegensburg
Regensburg,Deutschland
ISBN978-3-658-17865-9 ISBN978-3-658-17866-6(eBook)
DOI10.1007/978-3-658-17866-6
DieDeutscheNationalbibliothekverzeichnetdiesePublikationinderDeutschenNationalbibliografie;detaillier-
tebibliografischeDatensindimInternetüberhttp://dnb.d-nb.deabrufbar.
SpringerVieweg
©SpringerFachmedienWiesbadenGmbH2017
DasWerkeinschließlichallerseinerTeileisturheberrechtlichgeschützt.JedeVerwertung,dienichtausdrücklich
vomUrheberrechtsgesetzzugelassenist,bedarfdervorherigenZustimmungdesVerlags.Dasgiltinsbesondere
fürVervielfältigungen,Bearbeitungen,Übersetzungen,MikroverfilmungenunddieEinspeicherungundVerar-
beitunginelektronischenSystemen.
DieWiedergabevonGebrauchsnamen,Handelsnamen,Warenbezeichnungenusw.indiesemWerkberechtigt
auchohnebesondereKennzeichnungnichtzuderAnnahme,dasssolcheNamenimSinnederWarenzeichen-
undMarkenschutz-Gesetzgebungalsfreizubetrachtenwärenunddahervonjedermannbenutztwerdendürften.
DerVerlag,dieAutorenunddieHerausgebergehendavonaus,dassdieAngabenundInformationenindiesem
WerkzumZeitpunktderVeröffentlichungvollständigundkorrektsind.WederderVerlagnochdieAutorenoder
dieHerausgeberübernehmen,ausdrücklichoderimplizit,GewährfürdenInhaltdesWerkes,etwaigeFehler
oderÄußerungen.DerVerlagbleibtimHinblickaufgeografischeZuordnungenundGebietsbezeichnungenin
veröffentlichtenKartenundInstitutionsadressenneutral.
Lektorat:ThomasZipsner
GedrucktaufsäurefreiemundchlorfreigebleichtemPapier.
SpringerViewegistTeilvonSpringerNature
DieeingetrageneGesellschaftistSpringerFachmedienWiesbadenGmbH
DieAnschriftderGesellschaftist:Abraham-Lincoln-Strasse46,65189Wiesbaden,Germany
Vorwort
Foreverycomplexproblemthereisananswer
thatisclear,simple–andwrong
ParaphrasenachH.L.Mencken,1917
Die Inhalte des Buches bauen auf den Vorlesungen über lineare und nichtlineare Finite-
Elemente-Methoden(FEM)imBachelor-undMasterstudiumanderOstbayerischenTech-
nischen Hochschule Regensburg auf, die ich seit einigen Jahren halte, sowie auf meiner
langjährigenTätigkeitalsBerechnungsingenieurinderBlechumformsimulation.
SowohlinderLehrealsauchinderPraxisbestehtdieSchwierigkeitzuvermitteln,dass
vorgeblicheinfachzubedienendeBerechnungsprogrammenurzusinnvollenundverläss-
lichenLösungenführen,wennderAnwenderversteht,washinterderBenutzerschnittstelle
abläuftundmansichdeswegenauchalsAnwendermitdentheoretischenHintergründen
beschäftigenmuss.Diestrifftumsomehrzu,dakommerzielleFE-Programmeeinenstetig
wachsendenFunktionsumfanghaben,mitimmerkomplexerenInhalten.
DasBuchzieltdaraufab,einenmöglichstabgeschlossenenÜberblicküberdielineare
undnichtlinearedynamischeBerechnungvonStrukturenmitFinite-Elemente-Methoden
zu bieten. Der Fokus liegt auf einer Erklärung der Zusammenhänge, die für Anwender
kommerziellerSoftwarenotwendigsind.Dabeiistesunumgänglichvielesverkürztdarzu-
stellenoderauchThemenfelderkomplettauszusparen.AufdieDarstellungtheoretischer
Inhaltekannnichtverzichtetwerden,dieFEMisteinnumerischesVerfahren.Siewerden
abernursoweitdargestellt,wiefürdenAnwendererforderlich.DerEinsatzdertheoreti-
schenInhalteinderBlechumformsimulationmitLS-DYNAzeigtfüreinenpraxisnahen
Fall,dassdiesesWissennotwendigist,auchfürdenreinenBenutzereinesFE-Programms.
Die Literatur in diesem Gebiet ist unüberschaubar umfangreich. Um den Lehrbuch-
charakter zu erhalten, wird nur ein kleiner Ausschnitt von möglichen Literaturstellen
angegeben.EswirdimmerdiekonkreteStelleineinerReferenzgenannt,sodassderLeser
vondortausauchdieweitereLiteraturerschließenkann.
DieKapitel1bis8umfasseneinegenerelleEinführungindielineareFEM.Daraufbau-
endieKapitel9bis13auf,diediewesentlichennichtlinearenAspektederFEMenthalten,
mit einem Fokus auf dynamische Effekte. Im letzten Kapitel 14 werden die theoreti-
schen Inhalte der nichtlinearen FEM am Beispiel einer Blechumformsimulation in Ihrer
Anwendung erläutert. Weiterhin wird eine kurze Einführung in das FE-Programmpaket
LS-DYNAimAnhangangegeben.
DasBuchrichtetsichanStudierendeausBachelor-undMasterstudiengängenderInge-
nieurwissenschaftensowieanIngenieureinderIndustrie,dieinihrerberuflichenPraxis
mitFinite-Element-BerechnungenzutunhabenundaneinemtieferenEinblickinteressiert
v
vi Vorwort
sind,alserüblicherweiseinGrundlagenkursenoderSchulungenvermitteltwird.Anwen-
derdesProgrammsLS-DYNAsindangesprochen,dadiesesProgrammgenutztwird,um
anpraktischenBeispielendiezuvoreingeführtentheoretischenInhaltezuverdeutlichen.
DasselbegiltfürBerechnerimBereichderBlechumformsimulation,dadienichtlinearen
InhaltedesBuchesandiesemBeispieldemonstriertwerden.
AlsVoraussetzungsolltemanGrundkenntnissederTechnischenMechanikmitbringen,
wiesieinjedemIngenieurstudiumvermitteltwerden.VorkenntnisseüberFinite-Elemente-
Methodensindnichtnotwendig.
Zusatzunterlagen (Beispielprogramme, Kommandodateien etc.) können auf der Ver-
lagsseite www.springer.com direkt beim Buch (im Download-Bereich unterhalb des In-
haltsverzeichnisses)heruntergeladenwerden.
AmGelingeneinessolchenBuchprojektessindvielePersonenbeteiligt,ohnediedies
sonichtmöglichwäre.MeinherzlicherDankgiltHerrnDr.-Ing.IngoHeinle,HerrnDr.-
Ing.BerndHochholdinger,FrauMagdaMartins-Wagner,M.Sc.undHerrnProf.Dr.-Ing.
habil. Kai Willner für das Korrekturlesen und die vielen Hinweise. Der Ostbayerischen
TechnischenHochschuleRegensburgseigedanktfürdiegroßzügigeUnterstützungwäh-
rendderErstellungdesBuches.DemSpringerViewegVerlag,vorallemFrauImkeZander
undHerrnThomasZipsnerausdemLektoratMaschinenbau,dankeichfürdieangenehme
ZusammenarbeitunddiewertvollenAnregungen.
Regensburg,Juli2017 MarcusWagner
Inhaltsverzeichnis
1 Einleitung.......................................................... 1
1.1 GrundlegendeProblemstellungundLösungsmethoden ................ 1
1.2 VorgehensweisebeiBerechnungsaufgaben .......................... 4
1.3 EinsatzgebietederFinite-Elemente-Methode......................... 4
Literaturverzeichnis .................................................. 6
2 EinführungindielineareFEM ....................................... 7
2.1 GrundgedankederFEMamBeispieldesStabs ....................... 8
2.1.1 MathematischeBeschreibungdesphysikalischenSystems ....... 8
2.1.2 DiskretisierunginfiniteElemente ........................... 11
2.1.3 BerechnungderElementmatrizen............................ 14
2.2 DiskretisierungeinesFachwerks ................................... 15
2.2.1 DiskretisierunginfiniteElemente ........................... 15
2.2.2 TransformationvonnatürlichenaufglobaleKoordinaten ........ 16
2.2.3 ZusammenbaudesGleichungssystems ....................... 19
2.2.4 EinbringenvonRandbedingungen ........................... 22
2.2.5 LösendesGleichungssystems............................... 22
2.3 Beispiel:StabmitveränderlichemQuerschnitt ....................... 24
2.4 Aufgaben ...................................................... 27
3 MechanischeGrößenderStrukturmechanik ........................... 29
3.1 FormulierungdesRandwertproblems............................... 29
3.2 DerSpannungszustand ........................................... 31
3.2.1 DerSpannungsdeviator .................................... 34
3.3 ZugeordneterVerzerrungszustand.................................. 35
3.4 Gleichgewichtsbedingungen....................................... 37
3.5 Voigt-Notation.................................................. 37
3.6 Verallgemeinerteslinear-elastischesMaterialgesetznachHooke......... 38
3.7 EbenerSpannungs-undVerzerrungszustand ......................... 39
Literaturverzeichnis .................................................. 40
vii
viii Inhaltsverzeichnis
4 MathematischeModellierungüberEnergieprinzipien ................... 41
4.1 DasPrinzipvomMinimumdesGesamtpotenzials .................... 41
4.1.1 Einführungsbeispiel ....................................... 42
4.1.2 GesamtpotenzialeinesStabs................................ 45
4.1.3 AllgemeinesPrinzipvomMinimumdesGesamtpotenzials....... 46
4.2 DasPrinzipdervirtuellenVerschiebung............................. 47
4.3 MethodedergewichtetenResiduenamBeispieldesStabs.............. 52
Literaturverzeichnis .................................................. 54
5 DiskretisierungmitfinitenElementen ................................. 55
5.1 DefinitiondesNäherungsansatzesfüreinElement .................... 56
5.1.1 DieFormfunktionsmatrix .................................. 56
5.1.2 NäherungsansatzfürDehnungenundSpannungen.............. 58
5.2 DiskretisierungdesPrinzipsvomMinimumdesGesamtpotenzials ...... 58
5.3 DiskretisierungdesPrinzipsdervirtuellenVerschiebung............... 60
5.4 AufbaudesGesamtgleichungssystems .............................. 62
5.4.1 EigenschaftenderGesamtsteifigkeitsmatrix ................... 64
5.5 EinbringenvonRandbedingungen.................................. 65
5.6 LösunglinearerGleichungssysteme ................................ 67
5.6.1 DirekteGleichungslöser ................................... 67
5.6.2 IterativeGleichungslöser................................... 70
5.6.3 Modellreduktionstechniken................................. 72
5.7 Aufgaben ...................................................... 74
Literaturverzeichnis .................................................. 74
6 Finite-Elemente-Klassen ............................................. 75
6.1 KlassifizierungvonElementen .................................... 75
6.2 DasisoparametrischeKonzept..................................... 76
6.3 EindimensionaleElemente........................................ 81
6.3.1 Zweiknotiges,linearesStabelement .......................... 81
6.3.2 Dreiknotiges,quadratischesStabelement...................... 83
6.3.3 Balkenelemente .......................................... 84
6.4 ZweidimensionaleElemente ...................................... 89
6.4.1 Scheibenelement.......................................... 90
6.4.2 Schalenelemente.......................................... 95
6.5 DreidimensionaleElemente .......................................100
6.5.1 Hexaederelemente ........................................100
6.5.2 Pentaederelemente ........................................101
6.5.3 Tetraederelemente ........................................102
6.6 Aufgaben ......................................................103
Literaturverzeichnis ..................................................103
7 MathematischeundnumerischeAspektederFEM......................105
7.1 MathematischeAnforderungenanfiniteElemente ....................105
7.1.1 BedingungenfürdieKonvergenzderLösung ..................105
7.1.2 VerfahrenzurReduktiondesDiskretisierungsfehlers............107
Inhaltsverzeichnis ix
7.2 NumerischeIntegration ..........................................110
7.2.1 Newton-Cotes-Quadratur...................................111
7.2.2 Gauß-Quadratur ..........................................114
7.2.3 MehrdimensionaleIntegrale ................................115
7.2.4 AnwendungshinweisezumIntegrationsverfahren...............118
7.3 Elementversteifung(Locking) .....................................120
7.3.1 BeschreibungdesLocking-Effekts...........................120
7.3.2 MaßnahmenzurVermeidungvonElementversteifung...........121
7.3.3 Null-Energie-Moden ......................................124
7.4 Praxis-HinweisezurModellierung .................................125
7.4.1 Vernetzungsmethoden .....................................126
7.4.2 AnforderungenandieElementauswahlundVernetzung .........127
7.4.3 AusnutzungvonSymmetrienbeiderVernetzung...............129
7.5 Aufgaben ......................................................130
Literaturverzeichnis ..................................................130
8 LinearezeitabhängigeFEM ..........................................131
8.1 HerleitungderdynamischenFEMübervirtuelleVerschiebungen........132
8.2 NumerischeModalanalyse ........................................133
8.2.1 ModaleTransformation ....................................137
8.2.2 ModaleReduktion ........................................138
8.2.3 NäherungsweiseBerechnungdesEigenwertproblems ...........139
8.2.4 AnwendungsgebietederModalanalyse .......................141
8.3 BerücksichtigungvonDissipationseffekten ..........................141
8.3.1 ProportionaleundmodaleDämpfung ........................143
8.4 Frequenzganganalyse ............................................145
Literaturverzeichnis ..................................................148
9 GeometrischeNichtlinearität ........................................149
9.1 EinführungzurgeometrischenNichtlinearität ........................150
9.2 KinematischeBeschreibung.......................................151
9.2.1 Konfigurationen ..........................................151
9.2.2 Deformations-undVerschiebungsgradient ....................153
9.3 BeispieleeindimensionalerVerzerrungs-undSpannungsmaße..........155
9.4 AllgemeinenichtlineareVerzerrungsmaße...........................158
9.5 ZeitlicheAbleitungenderDeformation..............................161
9.6 TransformationvonVolumen-undFlächenelementen .................163
9.7 SpannungsmaßebeinichtlinearerBetrachtung .......................165
9.8 EnergieprinzipieninnichtlinearerForm.............................166
9.8.1 Upgedatete-Lagrange-Formulierung .........................167
9.8.2 Totale-Lagrange-Formulierung..............................169
9.8.3 Diskretisierung ...........................................170
Literaturverzeichnis ..................................................172
x Inhaltsverzeichnis
10 MaterielleNichtlinearität ...........................................173
10.1 ÜbersichtüberkonstitutiveBeziehungen ............................174
10.2 Eindimensionales,zeitunabhängiges,elastoplastischesVerhalten ........176
10.2.1 MathematischeFormulierung ...............................177
10.3 MehrachsigedehnratenunabhängigeElastoplastizität ..................183
10.3.1 DieFließbedingung .......................................184
10.3.2 DieFließregel ............................................188
10.3.3 DasVerfestigungsgesetz ...................................189
10.3.4 DieKonsistenzbedingungunddieMaterialtangente.............191
10.3.5 BerücksichtigungderDehnratenabhängigkeit..................192
10.4 NumerischeUmsetzungder J -Plastizität............................193
2
Literaturverzeichnis ..................................................196
11 Kontaktmodellierung................................................197
11.1 GrundlegendeBegriffe ...........................................198
11.1.1 BedingungfürNormalkontakt ..............................199
11.1.2 BehandlungvontangentialemGleiten ........................201
11.2 VerfahrenzurKontaktdetektion....................................202
11.3 Kontaktformulierungen...........................................205
11.3.1 KinematischeZwangsbedingungen(Multi-Point-Constraint) .....206
11.3.2 Penalty-Verfahren.........................................206
11.3.3 Lagrange-Multiplikator-Verfahren ...........................209
11.3.4 Augmented-Lagrange-Verfahren.............................210
Literaturverzeichnis ..................................................210
12 GleichungslösungbeinichtlinearenstatischenProblemen................211
12.1 Newton-Raphson-Verfahren.......................................212
12.2 AnwendungdesVerfahrensaufdieFEM ............................213
12.2.1 Linearisierung............................................214
12.2.2 Inkrementell-iterativesVerfahren ............................217
12.2.3 KonvergenzdesNewton-Raphson-Verfahrens..................218
12.2.4 HinweiszurZeitabhängigkeit ...............................219
12.3 WeitereVerfahrenfürnichtlineareGleichungssysteme.................220
Literaturverzeichnis ..................................................222
13 ZeitintegrationvonnichtlinearendynamischenProblemen...............223
13.1 Einführung.....................................................224
13.2 ImpliziteZeitintegrationnachdemNewmark-(cid:12)-Verfahren .............227
13.3 ExpliziteZeitintegrationnachdemzentralenDifferenzenverfahren ......231
13.3.1 PraktischeUmsetzungdesVerfahrens ........................235
13.3.2 Punktmassenmatrix .......................................236
13.3.3 NutzungquadratischerAnsatzfunktioneninexplizitenVerfahren .237
13.3.4 Stabilitätskriterium........................................238
13.3.5 MaßnahmenzurReduktionderRechenzeit....................242
13.3.6 DynamischeRelaxation....................................243
13.4 GegenüberstellungderbeidenZeitintegrationsverfahren ...............244