Numerical Methods in Contact Mechanics E-Book

Vladislav A. Yastrebov

• Herausgeber: John Wiley & Sons
• Kategorie: Fachliteratur
• Sprache: Englisch

Seitenzahl: 429

Veröffentlichungsjahr: 2013

Beschreibung

Computational contact mechanics is a broad topic which brings together algorithmic, geometrical, optimization and numerical aspects for a robust, fast and accurate treatment of contact problems. This book covers all the basic ingredients of contact and computational contact mechanics: from efficient contact detection algorithms and classical optimization methods to new developments in contact kinematics and resolution schemes for both sequential and parallel computer architectures. The book is self-contained and intended for people working on the implementation and improvement of contact algorithms in a finite element software. Using a new tensor algebra, the authors introduce some original notions in contact kinematics and extend the classical formulation of contact elements. Some classical and new resolution methods for contact problems and associated ready-to-implement expressions are provided. Contents: 1. Introduction to Computational Contact. 2. Geometry in Contact Mechanics. 3. Contact Detection. 4. Formulation of Contact Problems. 5. Numerical Procedures. 6. Numerical Examples. About the Authors Vladislav A. Yastrebov is a postdoctoral-fellow in Computational Solid Mechanics at MINES ParisTech in France. His work in computational contact mechanics was recognized by the CSMA award and by the Prix Paul Caseau of the French Academy of Technology and Electricité de France.

Leseprobe

Contents

Foreword

Preface

Notations

Chapter 1. Introduction to Computational Contact

1.1. Historical remark

1.2. Basics of the numerical treatment of contact problems

Chapter 2. Geometry in Contact Mechanics

2.1. Introduction

2.2. Interaction between contacting surfaces

2.3. Variations of geometrical quantities

2.4. Numerical validation

2.5. Discretized geometry

2.6. Enrichment of contact geometry

Chapter 3. Contact Detection

3.1. Introduction

3.2. All-to-all detection

3.3.Bucket sort detection

3.4. Case of unknown master–slave

3.5. Parallel contact detection

3.6. Conclusion

Chapter 4. Formulation of Contact Problems

4.1. Contact of a deformable solid with a rigid plane

4.2. Contact of a deformable solid with an arbitrary rigid surface

4.3. Contact between deformable solids

4.4. Variational equality and resolution methods

4.5. Penalty method

4.6. Method of Lagrange multipliers

4.7. Augmented Lagrangian Method

Chapter 5. Numerical Procedures

5.1. Newton’s method

5.2. Return mapping algorithm

5.3. Finite element method

5.4. Residual vectors and tangent matrices for contact elements

5.5. Method of partial Dirichlet–Neumann boundary conditions

5.6. Technical details

Chapter 6. Numerical Examples

6.1. Two dimensional problems

6.2. Three-dimensional problems

Appendix 1. Vectors, Tensors and s-Structures

A1.1. Fundamentals

A1.2. Vector space basis

A1.3. Sub-basis, vector function of v-scalar argument

A1.4. Tensors

A1.6. S-structures

A1.7. Reducedformof s-structures

Appendix 2. Variations of Geometrical Quantities

A2.1. First-order variations

A2.2. Second-order variations

Bibliography

Index

To Alexandra, Andrey and Daniel

First published 2013 in Great Britain and the United States by ISTE Ltd and John Wiley & Sons, Inc.

Apart from any fair dealing for the purposes of research or private study, or criticism or review, as permitted under the Copyright, Designs and Patents Act 1988, this publication may only be reproduced, stored or transmitted, in any form or by any means, with the prior permission in writing of the publishers, or in the case of reprographic reproduction in accordance with the terms and licenses issued by the CLA. Enquiries concerning reproduction outside these terms should be sent to the publishers at the undermentioned address:

ISTE Ltd

27-37 St George’s Road

London SW19 4EU

UK

www.iste.co.uk

John Wiley & Sons, Inc.

111 River Street

Hoboken, NJ 07030

USA

www.wiley.com

© ISTE Ltd 2013

The rights of Vladislav A. Yastrebov to be identified as the author of this work have been asserted by him in accordance with the Copyright, Designs and Patents Act 1988.

Library of Congress Control Number: 2012954082

British Library Cataloguing-in-Publication Data

A CIP record for this book is available from the British Library

ISBN: 978-1-84821-519-1

Foreword

The area of contact mechanics has become a vivid research field since the modeling of engineering problems has become a lot more sophisticated. This is due to the available computing power that has led to more refined models including contact constraints. This book by Vladislav Yastrebov is related to this still emerging area of contact mechanics. It starts from the basic principles and geometrical relations of contact mechanics and then moves to the essential issue, the detection of contact constraints. Here, the book includes valuable help for those who want to implement contact algorithms since the detection procedures are described in detail, including many exceptions. The formulation of contact problems, again, provides many insights into the complex contact behavior and gives a complete overview with respect to frictionless and frictional contact. This is also true for the chapter that describes the numerical procedures for contact problems. These are discussed in detail and well presented such that the reader will understand the different approaches that can be applied to solve nonlinear contact.

The book will be useful as an introduction to contact mechanics and related algorithms for graduate students who have the necessary background in mathematics and continuum mechanics. However, the book is also a reliable and comprehensive source for researchers who are interested in implementing algorithms and discretization schemes for the solution of nonlinear contact problems. Last but not least, design engineers from industry can use this book as background information for contact analyses related to, e.g. forming, forging and other problems that involve contact and friction.

Prof. Dr. Ing. Peter WRIGGERS December 2012 Institute of Continuum Mechanics Leibniz Universität Hannover

Preface

Nowadays, contact and friction are particularly important for our civilization. Think, for example, about car brakes, wheel–rail contact, assembled pieces in engines and turbines, bearings and gears in mechanical devices and electromechanical contacts. Numerical simulations permit us to study and improve these complex systems involving contact, friction and wear. This book answers the question of what is behind these simulations and uncovers for readers the underlying machinery of the finite element analysis in contact mechanics.

Regardless of the prevalence of contact and friction, these phenomena are hard to study experimentally because of their multiscale/multiphysical nature and the inaccessibility of contact interfaces to direct observations. Likewise, these problems are challenging for numerical treatment due to the particularity of contact and friction conditions and complexity of involved algorithms. Moreover, the mathematics and non-trivial notions introduced in this branch of computational mechanics are hard to comprehend for beginners. Thus, the first motivation of this book is to introduce new people to this field. The second motivation is to expose all involved components of computational contact in its integrity and interconnection: geometry, detection and resolution. And finally, I would like to expose some original developments in computational contact mechanics.

I address this book to students, engineers and researchers who solve contact problems by means of the finite element method. Also, I am aiming at developers wishing to implement or improve contact algorithms in their commercial or in-house finite element software. To make the book accessible to people unfamiliar with basics of the computational contact, I shall introduce all terms and notions and give many examples. For all developments in contact geometry I used a new tensor algebra, so some effort are needed from the reader to “get used to it”. But I believe that for readers familiar with programming, this novelty should not present a difficulty, because the main concept is transparent – array of arrays.

Contact algorithms are rich in details that are seemingly negligible but are crucial for the robustness and accuracy of the code. So, based on our experience, I made an attempt to expose most of them. Furthermore, as the implementation of contact algorithms is delicate (both for contact detection and resolution steps), it requires an extended validation and testing. For that purpose, I standardized and exposed many tests from the literature and suggested some new ones.

I hope that this book will introduce new people to the field of computational contact mechanics and that the ideas expressed here will engender the development of new methods and approaches to make the simulation of contact more reliable and accurate.

This book would not be possible without the help and encouragement of Georges Cailletaud, Frédéric Feyel and the financial support of CNRS and SNECMA, which I gratefully acknowledge. I express my thanks to my dear wife Alexandra, my sons, Andrey and Daniel, my parents and my brother for their love, patience and comprehension. I am grateful to my colleagues Djamel Missoum-Benziane and Nikolay Osipov for their constant support, help and friendship. I also acknowledge André Pineau, Jacques Besson and Samuel Forest for creating a stimulating scientific atmosphere. I thank also Liliane Locicero, Konaly Sar, Odile Adam, and Anne Piant for their permanent administrative help, empathy and friendly attention.

Vladislav A. YASTREBOVCentre des MatériauxMINES ParisTechEvryDecember 2012

Notations

Vectors and tensors

– Scalar (zero-order tensor) – small Latin and Greek letters:

a, α, b, . . .

– Vector (first-order tensor) – underlined small bold Latin and Greek letters:

c,β,d, . . .

– Second-order tensor – capital bold Latin letters underlined twice:

– Higher order tensor – capital bold Latin letters underlined twice with upper left index of order:

V-Vectors and V-tensors

– V-scalar (“vector of scalars”) – small Latin and Greek letters underlined by a wave:

– V-vector (“vector of vectors”) – small Latin and Greek letters underlined by a line and a wave:

– V-tensor (“vector of tensors”) – capital bold Latin letters underlined by a double line and a wave:

T-Vectors and T-tensors

– T-scalar (“tensor of scalars”) – capital bold Latin letter underlined by a double wave:

– T-vector (“tensor of vectors”) – small Latin and Greek letters underlined by a line and a double wave:

– T-tensor (“tensor of tensors”) – capital bold Latin letters underlined by a double line and a double wave:

Vector and tensor operations

– : Euclidean norm of a vector;

– det : determinant of a tensor;

– : unit tensor;

– : unit t-scalar;

– : trace of a tensor;

– : inverse of a tensor;

– : transpose of a tensor;

– : scalar or dot product;

– : vector or cross product;

– : tensor product;

– : tensor contraction.

Other operations

– : full time derivative;

– : first variations;

– : full first variations;

– : second variation;

– : full second variation;

– : gradient;

– : divergence;

– : rotor.

Miscellaneous

Abbreviations

– PM: penalty method;

– LMM: Lagrange multiplier method;

– ALM: augmented Lagrangian method;

– FEM: finite element method;

– FEA: finite element analysis;

– CAD: computer-aided design;

– NTN: node-to-node;

– NTS: node-to-segment discretizations;

– MPC: multi-point constraints;

– PDN: partial dirichlet–neumann;

– SDMR: single detection multiple resolution;

– MDMR: multiple detection multiple resolution.

Remark. Macaulay brackets, dist (., .) and θ (.) functions.

Throughout the book, we use the notation of Macaulay brackets.

Theθfunction is a similar notation widely used in both engineering and mathematical literature:

or a more general dist(, ) function:

where dist(x,∂Ω) is a somehow defined distance from pointxto the closure of the setΩ. For example, in the simplest case

All these functions are equivalent for the considered case and interchangeable, so the reader is invited to interpret the Macaulay brackets as one of the abovementioned functions to which he/she is more accustomed.

Chapter 1

Introduction to Computational Contact

From a mechanical point of view, at macroscopic scale, contact is the notion for the interaction between bodies coming in to touch and exchanging load and energy (heat and electric charge). The physics of contact – rich and complex – is hard to study because of its multiscale and multiphysical nature and also because of the contact zone’s inaccessibility for direct observations1. Both experimental and numerical investigations of the contact experience difficulties. Tribology is an experimental science that describes and characterizes the contact, adhesion, friction, wear and lubrification as well as the involved mechanical, physical and chemical effects at different scales. Mathematics formalizes these descriptions by some measurable quantities (the coefficient of friction, the real contact area, the heat transfer coefficient, etc.). On the basis of observations, we suggest some models for the evolution of these quantities and then integrate them in complete numerical models to study particular systems at a common basis. The simulations thus imply strong simplifications2 that may be crucial for their validity. On the other hand, even the simplest models appear complicated for numerical simulations both from mathematical and programming points of view. All together, the oversimplification of the phenomenon and the imperfect implementation of established models may easily lead to incorrect results. The aim of this book is to resolve some of these problems, provide a consistent basis for the numerical treatment of contact problems at all stages, avoid unnecessary simplifications and enhance existing numerical models. The book focuses on the mechanics of contact and its numerical treatment by the finite element method (FEM). The underlying physics and mathematics are covered only partly and superficially3.

We suppose that the reader is unfamiliar with this field, so from the very beginning we provide the reader with the vocabulary and all the basic notions of computational contact mechanics. Many notions are very technical and complex, so we explain them in many ways throughout the book4 and often enhance these explanation with figures and examples. The material covers three subjects: geometry, detection and resolution. However, these parts are hard to comprehend independently. To ensure the readability of the book, this chapter introduces the subject, provides the global context of the numerical treatment and briefly discusses every component.

Contact problems in mechanics of deformable solids can be singled out into a particular class. First, they imply a discontinuity: the contact occurs at the interface between two separate continuous bodies. Second, the contact constraints at this interface cannot be replaced by ordinary boundary conditions imposed on both the contacting surfaces. Third, the contact interface itself cannot be simply considered as an internal surface. In an idealized case, the contact interface is a zero thickness layer, which sustains only compressive stress in the direction orthogonal to the contact interface (Figure 1.1(a)); any stretching leads to the vanishing of the contact interface (Figure 1.1(b)). In the case of frictionless contact, the contact interface contrary to an ordinary internal surface, does not sustain any tangential efforts, which allows two solids to slide relative to each other (Figure 1.1(a)). In the case of frictional contact, things become more complex. For example, in case of the classic Coulomb’s friction law, in stick state, the contact interface under pressure is similar to an internal interface – no separation, no sliding – locally both surfaces remain glued to each other (Figure 1.1(c)). However, if locally we reach a critical shear stress, the surfaces start to slip relatively to each other (). It follows from this simple representation that the contribution of the contact interface to the energy of the system is always zero except in the case of frictional slip.