Nonconforming boundary elements and finite elements for interface and contact problems with friction [Elektronische Ressource] : hp-version for mortar, penalty and Nitsche's methods / von Alexey Chernov

Alexey Chernov

Sujets

• Physik

Informations

• Publié par: gottfried_wilhelm_leibniz_universitat_hannover
• Publié le: 01 janvier 2006
• Nombre de lectures: 20
• Langue: English
• Poids de l'ouvrage: 6 Mo

Extrait

Nonconforming boundary elements
and finite elements for interface
and contact problems with friction –
hp-version for mortar,
penalty and Nitsche’s methods
Von der Fakulta¨t fu¨r Mathematik und Physik
der Universitat Hannover¨
zur Erlangung des Grades
Doktor der Naturwissenschaften
Dr. rer. nat.
genehmigte Dissertation
von
Dipl.-Math. Alexey Chernov
geboren am 13.08.1981 in Moskau
2006Referent: Prof. Dr. E. P. Stephan, Universita¨t Hannover
Korreferent: Prof. Dr. P. Hansbo, Chalmers University of Technology, Go¨teborg
Korreferent: Prof. Dr. J. Gwinner, Universitat der Bundeswehr, Munchen¨ ¨
Tag der Promotion: 27.06.2006Abstract
Theobjectiveofthisthesisistheconstruction,analysisandimplementationofhighorder
FE,BEandFE/BEcoupling methods forinterface andfrictionalcontact problems with
nonmatching discretizations, which have a wide industrial application.
A new hp-Nitsche’s FE/BE coupling method for interface problems is designed and
analysed. The method is proven to be consistent and stable, independent of the dis-
cretization parameters. A priori error analysis shows that the method is optimal in h
and is suboptimal in p on quasiuniform meshes.
The question of unique solvability is addressed for the one-body contact problem with
Tresca’s friction. Constructing a chain of equivalent formulations the frictional contact
problem is approximated with a sequence of frictionless contact problems (Uzawa algo-
rithm). Conditions for the convergence of the algorithm are obtained. An hp-penalty
BEMforone-bodyfrictionlesscontact isdeveloped. Astheapriorierroranalysisshows,
1−ǫthe penalty parameter ε must be chosen proportional to (h/p) , for optimal conver-n
gence rate (in the energy norm) of the discrete penalty solution to the solution of the
original variationalinequality formulation. A residual based a posteriori error estimator
for the h-version of penalty FE and BE for one-body contact with Tresca friction is
investigated. The error estimator, motivated so far with heuristical arguments and only
for FEM, is shown to be reliable and efficient for both FEM and BEM.
The two-body elastoplastic contact problem with Coulomb’s law of friction is solved
with the FE/BE coupling and pure BE methods. The incremental loading procedure
with Newton iterations on each loading step is used. Linearization of the frictional con-
tact and plasticity terms as well as a description of the solution procedure are given in
detail. Theresidualaposteriorierrorestimate, obtainedforone-bodyfrictionalcontact,
is generalized to this two-body frictional contact problem. A novel hp-mortar method
for two-body contact with Tresca’s friction is designed and analysed for a variational
inequality formulation. The contact constraints are imposed on the discrete global set
of Gauss-Lobatto points. The nonmatched meshes are connected in terms of the hp-
1/4mortar projection. The a priori error analysis shows the convergence rate O((h/p) )
in the energy norm under additional assumptions on the discretization parameters. A
Dirichlet-to-NeumannalgorithmandanUzawa algorithmareused tosolve theproblem.
A heuristically motivated error indicator is used to perform an hp automatic refine-
ment procedure. The h-version of the constructed method is extended onto two-body
elastoplastic frictional contact problems and is compared to the results provided by the
penalty method.
The theoretical results are supported by numerical benchmark computations.
Key words. frictional contact, interface problems, finite elements, boundary elements,
FE/BEcoupling,hp-methods,apriorierror,aposteriorierror,mortar,penalty,Nitsche’s
method
3Zusammenfassung
DasZieldieser DissertationistdieKonstruktion,AnalyseundImplementierung vonFE,
BE und FE/BE Kopplungsverfahren fur Interface- und Reibungskontaktprobleme mit¨
unpassenden Diskretisierungen.
Eine neue hp-Nitsche FE/BE Kopplungsmethode fur Interface-Probleme wird konstru-¨
iert und analysiert. Es wird bewiesen, dass das Verfahren konsistent und stabil ist,
unabha¨ngig von den Diskretisierungsparametern. Die durchgefu¨hrte a priori Fehlerana-
lyse zeigt, dass die Methode optimal in h und suboptimal in p auf den quasiuniformen
Gittern ist.
DieFragedereindeutigenLosbarkeitwirdfurdasEin-Korper-KontaktproblemmitTres-¨ ¨ ¨
ca Reibung untersucht. Die Losung des Ausgangsproblems mit Reibung wird durch eine¨
FolgevonreibungslosenProblemenapproximiert(UzawaAlgorithmus).DieKonvergenz-
bedingungen fur den Algorithmus werden hergeleitet. Eine hp-Penalty BE Methode fur¨ ¨
das Ein-Ko¨rper-Kontaktproblem wird entwickelt. Wie die a priori Fehleranalyse zeigt,
1−ǫmuß der Penalty-Parameter ε proportional zu (h/p) , gew¨ahlt werden, um die opti-n
male Konvergenzordnung (in der Energienorm) der diskreten Penalty-Lo¨sung gegen die
exakte L¨osung der variationellen Ungleichung zu erreichen. Als na¨chstes wird ein resi-
dueller a posteriori Fehlerscha¨tzer fu¨r die h-Versionen von FEM und BEM untersucht.
FurdenFehlerschatzer, derbisher nurmitheuristischen Argumentenmotiviert undaus-¨ ¨
schließlich fur FEM benutzt wurde, wird bewiesen, dass er zuverlassig und effizient ist.¨ ¨
Ferner werden die Zwei-Ko¨rper-Kontaktprobleme mit Coulomb’scher Reibung fu¨r die
h-Versionen von FE/BE und reinem BE Verfahren betrachtet. Die inkrementelle Last-
aufbringung mit dem Newton-Verfahren in jedem Iterationsschritt wird eingesetzt. Die
Linearisierung der Reibungskontaktterme und der Plastizita¨tsterme sowie die Beschrei-
bung der L¨osungsprozedur werden detailliert angegeben. Die residuelle a posteriori Feh-
lerabschatzung, die im Falle des Ein-Korper-Reibungskontaktproblems gewonnen wur-¨ ¨
de, wird auf ein Zwei-Korper-Reibungskontaktproblem verallgemeinert. Eine neue hp-¨
Mortar Methode fur das Zwei-Korper-Kontaktproblem mit Tresca Reibung wird kon-¨ ¨
struiert und die variationelle Ungleichung analysiert. Die Kontaktbedingungen sind auf
der diskreten globalen Menge der Gauß-Lobatto Knoten definiert. Nichtpassende Git-
ter sind durch die hp-Mortarprojektion verbunden. Die a priori Fehleranalyse zeigt die
1/4KonvergenzordnungO((h/p) )inderEnergienormunterzusa¨tzlichen Bedingungenfu¨r
die Diskretisierungsparametern. Dirichlet-zu-Neumann Verfahren und Uzawa Verfahren
werden als Lo¨sungsprozedur benutzt. Ein Fehlerindikator wird heuristisch begru¨ndet
und in einer automatischen hp Gitterverfeinerungsprozedur eingesetzt. Die h-Version
des obigen Verfahrens wird auf Zwei-Korper elastoplastische Kontaktprobleme mit Rei-¨
bung generalisiert und mit den Ergebnissen des Penalty-Verfahrens verglichen.
Die theoretische Ergebnisse werden durch die Benchmark-Rechnungen unterstu¨tzt.
Schlagworte: Reibungskontakt, Interface-Probleme, finite Elemente, Randelemente,
FE/BE Kopplung, hp-Methoden, a priori, a posteriori, Mortar, Penalty, Nitsche Ver-
fahren
4Acknowledgements
It is a great pleasure for me to thank my advisor, Prof. Dr. Ernst P. Stephan, for
inspiring me to work in this area. I am very grateful to him for his intensive guidance
of my work, for helpful discussions and important remarks. Also I would like to thank
PD Dr. Matthias Maischak for his support and help in my theoretical investigations
as well as in the programming. His software package maiprogs became a basis for
implementation of numerical experiments, presented in the thesis.
Iamverygratefultomyco-advisorandmyco-refereeProf. Dr. PeterHansbo(Chalmers
University ofTechnology, Gothenburg, Sweden) forhishospitalityandnumerous discus-
sions during my visit in Gothenburg. I would like to thank my co-referee Prof. Dr.
Joachim Gwinner (Universitat der Bundeswehr, Munchen, Germany) for his support¨ ¨
and readiness to examine this thesis in a very short period of time. I would also like to
thank Prof. Dr. Patrick Hild (Laboratoire de Math´ematiques, Universit´e de Franche-
Comt´e, Besanc¸on, France) for his comments and sending me several important papers.
Furthermore, I would like to give my thanks to the members of our working group
in the Institute of Applied Mathematics, University of Hanover, for the friendly and
stimulating atmosphere. My special thanks go to Dipl.-Math. Sergey Geyn for his
cooperation, particularly for performing several numerical examples from Section 4.2.7.
Also I am very grateful to all members of the DFG Graduiertenkolleg 615 for their
collaboration and interest in my work.
Finally, I would like to thank wholeheartedly my family and give a very special thanks
to my wife Natalia for her strong support and patience during my work on this thesis.
This thesis was supported by the DFG Graduiertenkolleg 615 ”Interaction of Modeling,
Computation Methods and Software Concepts for Scientific-Technological Problems”.
Hanover, May 2006 Alexey Chernov
56Contents
Introduction 13
1 Foundations 19
1.1 Sobolev spaces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
1.2 Boundary integral operators for elliptic problems . . . . . . . . . . . . . 22
1.3 Symmetricboundaryelementformulationformixedboundaryvalueprob-
lems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
1.4 Discretization of the St