An example of nonuniqueness for the continuous static unilateral contact model with Coulomb friction

Similar documents
A mixed finite element method and solution multiplicity for Coulomb frictional contact

A uniqueness criterion for the Signorini problem with Coulomb friction

On finite element uniqueness studies for Coulomb s frictional contact model

On Mixed Methods for Signorini Problems

The Mortar Finite Element Method for Contact Problems

ON THE ANALYSIS OF A VISCOPLASTIC CONTACT PROBLEM WITH TIME DEPENDENT TRESCA S FRIC- TION LAW

An error estimate for the Signorini problem with Coulomb friction approximated by finite elements

An optimal control problem governed by implicit evolution quasi-variational inequalities

A domain decomposition algorithm for contact problems with Coulomb s friction

Solution of contact problems with friction

Some improvements of Xfem for cracked domains

DYNAMIC CONTACT WITH SIGNORINI S CONDITION AND SLIP RATE DEPENDENT FRICTION

Fundamentals of Linear Elasticity

PROPAGATION OF A MODE-I CRACK UNDER THE IRWIN AND KHRISTIANOVICH BARENBLATT CRITERIA

Legendre-Fenchel duality in elasticity

Expression of Dirichlet boundary conditions in terms of the strain tensor in linearized elasticity

EXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON

Second-gradient theory : application to Cahn-Hilliard fluids

Existence of minimizers for the pure displacement problem in nonlinear elasticity

Transactions on Engineering Sciences vol 14, 1997 WIT Press, ISSN

LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)

Série Mathématiques. MIRCEA SOFONEA Quasistatic processes for elastic-viscoplastic materials with internal state variables

ABHELSINKI UNIVERSITY OF TECHNOLOGY

VARIATIONAL ANALYSIS OF A ELASTIC-VISCOPLASTIC CONTACT PROBLEM WITH FRICTION AND ADHESION

Weak Convergence Methods for Energy Minimization

New constructions of domain decomposition methods for systems of PDEs

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1

COMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Computational Methods in Elasticity - Michel SALAUN

Measurement of deformation. Measurement of elastic force. Constitutive law. Finite element method

Existence and uniqueness of the weak solution for a contact problem

arxiv: v1 [math.na] 27 Jan 2016

Analytical formulation of Modified Upper Bound theorem

MATHEMATICAL MODELLING OF HUMAN LIMB IN MATLAB

Spectrum and Exact Controllability of a Hybrid System of Elasticity.

Mechanics PhD Preliminary Spring 2017

quantitative information on the error caused by using the solution of the linear problem to describe the response of the elastic material on a corner

Two-scale Dirichlet-Neumann preconditioners for boundary refinements

LIMIT LOAD OF A MASONRY ARCH BRIDGE BASED ON FINITE ELEMENT FRICTIONAL CONTACT ANALYSIS

PREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO

A new algorithm for solving 3D contact problems with Coulomb friction

The Finite Element Method for Computational Structural Mechanics

On the uniform Poincaré inequality

Intrinsic finite element modeling of a linear membrane shell problem

Generalized Korn s inequality and conformal Killing vectors

Instabilities and Dynamic Rupture in a Frictional Interface

A FINITE ELEMENT FORMULATION FOR NONLINEAR 3D CONTACT PROBLEMS

ON THE ASYMPTOTIC BEHAVIOR OF ELLIPTIC PROBLEMS IN PERIODICALLY PERFORATED DOMAINS WITH MIXED-TYPE BOUNDARY CONDITIONS

Existence and Uniqueness of the Weak Solution for a Contact Problem

An Elastic Contact Problem with Normal Compliance and Memory Term

Application of Regularized Digital Images Correlation Analysis to the Tensile Test of NiTi Shape Memory Alloy

MAE4700/5700 Finite Element Analysis for Mechanical and Aerospace Design

ON SOME ELLIPTIC PROBLEMS IN UNBOUNDED DOMAINS

Granular Micro-Structure and Avalanche Precursors

Topology Optimization of Compliant Mechanism with Geometrical Advantage

Constitutive models: Incremental plasticity Drücker s postulate

PEAT SEISMOLOGY Lecture 2: Continuum mechanics

Zanaboni Theorem and Saint-Venant s Principle

2D OR 3D FRICTIONAL CONTACT ALGORITHMS AND APPLICATIONS ln A LARGE DEFORMA TION CONTEXT

Homogenization limit for electrical conduction in biological tissues in the radio-frequency range

LINEARIZED ELASTICITY AND KORN S INEQUALITY REVISITED

The Rotating Inhomogeneous Elastic Cylinders of. Variable-Thickness and Density

EXISTENCE OF SOLUTIONS FOR A CLASS OF VARIATIONAL INEQUALITIES

Glowinski Pironneau method for the 3D ω-ψ equations

Prediction of geometric dimensions for cold forgings using the finite element method

ON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR TERM

Partitioned Formulation for Solving 3D Frictional Contact Problems with BEM using Localized Lagrange Multipliers

Duality method in limit analysis problem of non-linear elasticity

Generalized Budan-Fourier theorem and virtual roots

Debonding process in composites using BEM

Intrinsic finite element modeling of a linear membrane shell problem

Optimized Surface Acoustic Waves Devices With FreeFem++ Using an Original FEM/BEM Numerical Model

FRICTIONAL CONTACT PROBLEM IN DYNAMIC ELECTROELASTICITY

SIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED

Determination of thin elastic inclusions from boundary measurements.

Optimal shape and position of the support for the internal exact control of a string

A Domain Decomposition Algorithm for Contact Problems: Analysis and Implementation

Interface and contact problems

Solution of contact problems in linear elasticity using a feasible interior point algorithm for nonlinear complementarity problems

arxiv: v1 [math.ap] 10 May 2013

Chapter 3. Load and Stress Analysis

Generalized Budan Fourier theorem and virtual roots

Vector hysteresis models

Basic concepts to start Mechanics of Materials

c 2014 Society for Industrial and Applied Mathematics

Transactions on Engineering Sciences vol 6, 1994 WIT Press, ISSN

A DUALITY ALGORITHM FOR THE OBSTACLE PROBLEM

Domain optimisation using Trefftz functions application to free boundaries

HANDBUCH DER PHYSIK HERAUSGEGEBEN VON S. FLÜGGE. BAND VIa/2 FESTKÖRPERMECHANIK II BANDHERAUSGEBER C.TRUESDELL MIT 25 FIGUREN

Shock enhancement of cellular structures under impact loading: Part II Analysis

Arbitrary Normal and Tangential Loading Sequences for Circular Hertzian Contact

Stator/Rotor Interface Analysis for Piezoelectric Motors

Generalized Budan-Fourier theorem and virtual roots

Optimal Slope of Dramix Type Fibers in Reinforced Concrete

Models for dynamic fracture based on Griffith s criterion

QUASISTATIC INDENTATION OF A RUBBER COVERED ROLL BY A RIGID ROLL - THE BOUNDARY ELEMENT SOLUTION

NUMERICAL ANALYSIS OF CERTAIN CONTACT PROBLEMS IN ELASTICITY WITH NON-CLASSICAL FRICTION LAWS

Finite Element Method in Geotechnical Engineering

The effects of a discontinues weight for a problem with a critical nonlinearity

Finite Elements for Elastic Shell Models in

Mathematical Background

Transcription:

An example of nonuniqueness for the continuous static unilateral contact model with Coulomb friction Un exemple de non-unicité pour le modèle continu statique de contact unilatéral avec frottement de Coulomb Patrick Hild a a Laboratoire de Mathématiques de Besançon, Université de Franche-Comté / CNRS UMR 6623, 16 route de Gray, 25030 Besançon, France. Abstract The aim of this note is to propose an example of nonuniqueness for the continuous static unilateral contact model with Coulomb friction in linear elasticity. To cite this article: P. Hild, C. R. Acad. Sci. Paris, Ser. I 336 (2003). Résumé Le but de cette note est de proposer un exemple de non-unicité pour le modèle continu statique de contact unilatéral avec frottement de Coulomb en élasticité linéaire. Pour citer cet article : P. Hild, C. R. Acad. Sci. Paris, Ser. I 336 (2003). 1. Introduction The Coulomb model [2] is the most common law of friction used in solid mechanics in order to describe slipping or sticking bodies on a contact surface. This law is very often coupled with the unilateral contact model which takes into account the possible separation of the body from the surface. In the simple case of elastostatics, the variational formulation of the unilateral contact problem with Coulomb friction (see [3,]) was followed some years later by an existence result in the case of an infinitely long strip with small friction (see [8]). These results were generalized with more classical geometries and the bounds ensuring existence were improved, particularly in references [7] and [5]. Nevertheless this (simple) model shows Email address: patrick.hild@math.univ-fcomte.fr (Patrick Hild). Preprint submitted to Elsevier Science 3rd June 2007

numerous mathematical difficulties so that there does not exist, to our knowledge neither uniqueness results nor nonuniqueness or nonexistence examples. The aim of this note is to propose a simple example of nonuniqueness for the continuous unilateral contact model with Coulomb friction for a linear elastic body lying on a rigid foundation. This example admits at least two solutions provided the friction coefficient is greater than a critical value. Moreover these two solutions (one which sparates from the foundation and another one corresponding to stick on the foundation) do not depend on the friction coefficient. 2. Problem set-up Let us consider the deformation of an elastic body occupying, in the initial unconstrained configuration a domain Ω in R 2. The boundary Ω of Ω consists of Γ D, Γ N and Γ C where the measure of Γ D does not vanish. The body Ω is submitted to given displacements U on Γ D, it is subjected to surface traction forces F on Γ N and the body forces are denoted f. In the initial configuration, the part Γ C is considered as the candidate contact surface on a rigid foundation which means that the contact zone cannot enlarge during the deformation process. The contact is assumed to be frictional and the stick, slip and separation zones on Γ C are not known in advance. We denote by µ > 0 the given friction coefficient on Γ C. The unit outward normal and tangent vectors on Ω are n = (n x, n y ) and t = ( n y, n x ) respectively. The unilateral contact problem with Coulomb s friction law consists of finding the displacement field u : Ω R 2 satisfying (1) (6): div σ(u) + f = 0 in Ω, σ(u) = C ε(u) in Ω, (2) u = U on Γ D, (3) σ(u)n = F on Γ N. () The notation σ(u) : Ω S 2 represents the stress tensor field lying in S 2, the space of second order symmetric tensors on R 2. The linearized strain tensor field is ε(u) = ( u + T u)/2 and C is the fourth order symmetric and elliptic tensor of linear elasticity. Afterwards we choose the following notation for any displacement field u and for any density of surface forces σ(u)n defined on Ω: u = u n n + u t t and σ(u)n = σ n (u)n + σ t (u)t. On Γ C, the three conditions representing unilateral contact are as follows u n 0, σ n (u) 0, σ n (u) u n = 0, (5) and the Coulomb friction law on Γ C is described by the following conditions: u t = 0 = σ t (u) µ σ n (u), u t 0 = σ t (u) = µ σ n (u) u t u t. (6) Remark. Let us mention that the true Coulomb friction law involves the tangential contact velocities and not the tangential displacements. However, a problem analogous to the one discussed here is obtained by time discretization of the quasi-static frictional contact evolution problem. In this case (see [1]) u, f 2 (1)

and F stand for u((i + 1) t), f((i + 1) t) and F((i + 1) t) respectively and u t has to be replaced by u t ((i + 1) t) u t (i t), where t denotes the time step. For simplicity and without any loss of generality only the static case described above will be considered in the following. The variational formulation of problem (1) (6) consists of finding u K satisfying (see [3,]): a(u, v u) µσ n (u)( v t u t ) dγ L(v u), v K, (7) Γ C where a(u, v) = Ω (Cε(u)) : ε(v) dω, L(v) = Ω f.v dω + F.v dγ, Γ N for any u and v in the Sobolev space (H 1 (Ω)) 2. In these definitions the notations and : represent the canonical inner products in R 2 and S 2 respectively. In (7), the set K stands for the convex of admissible displacement fields: } K = {v (H 1 (Ω)) 2 ; v = U on Γ D, v n 0 on Γ C. As far as we know there only exist existence results in the case of small friction coefficients (see [8,7,5]) for problem (1) (6) and there are neither uniqueness results (unless the loads f,f and U are equal to zero) nor nonuniqueness or nonexistence examples available. 3. A nonuniqueness example We consider the triangle Ω of vertexes A = (0, 0), B = (1, 0) and C = (3/, 1/) and we define Γ D = [B, C], Γ N = [A, C], Γ C = [A, B]. The body Ω lies on the rigid foundation, the half-space delimited by the straight line (A, B). We suppose that the body Ω is governed by Hooke s law concerning homogeneous isotropic materials so that (2) becomes σ(u) = Eν (1 2ν)(1 + ν) tr(ε(u))i + E 1 + ν ε(u) where I represents the identity matrix, tr is the trace operator, E and ν denote Young s modulus and Poisson ratio, respectively. The chosen material characteristics are ν = 1/5 and E = 1 (the choice of E is only made for the sake of simplicity and any choice of a postive E would lead to the same kind of nonuniqueness example). Let (x = (1, 0), y = (0, 1)) stand for the canonical basis of R 2. We suppose that the volume forces f = (f x, f y ) = (0, 0) are absent and that the surface forces denoted F = (F x, F y ) are such that F x = 35 10 α, F y = 0, 8 where α > 0. On Γ D, the presribed displacements U = (U x, U y ) are given by U x = 6α(x 1), U y = 3 α(x 1). where α < /3 (to avoid some penetration of Γ D in the rigid foundation). Set 0 < α < /3 and introduce two linear displacement fields u = (u x, u y ) and u = (u x, u y ) in Ω: u x = (7x + y 7)α, u y = ( x 7 y + 1 ) α, (8) 3

u x = 6yα, u y = 3 yα. (9) The displacement field u moves points A and C to ( 7α, α) and (3/ 3α/2, 1/ 3α/16) respectively whereas position of point B remains unchanged. When considering u the points A and B are stuck and the new position of point C becomes (3/ 3α/2, 1/ 3α/16). In the next proposition we show that the two displacement fields u and u are solutions of the frictional contact problem(1) (6) if the friction coefficient µ is large enough. Proposition. Let be given Ω, Γ D, Γ N, Γ C, E, ν, f, F, U, α as previously. For any µ 3 there exist at least two solutions (given by (8) and (9)) of the Coulomb frictional contact problem (1) (6). Proof. Using the constitutive relation (2), one easily obtains: 175 σ(u) = α 2 0, 0 0 5 σ(u) = α 2 5 2 5. 2 5 6 As a result the equilibrium equations div σ(u) = 0 and div σ(u) = 0 are satisfied in Ω. On Γ D, the two fields u and u coincide with the prescribed displacement field U and the stress vectors σ(u)n and σ(u)n are equal to F on Γ N (since n = ( 1/ 10, 3/ 10)). It remains to verify the fulfillment of the frictional contact conditions for both fields. We begin with u. On Γ C, n = (0, 1), t = (1, 0) and we get σ n (u) = 0, σ t (u) = 0, u n = α(x 1), u t = 7α(x 1). Since α > 0 it follows that u n 0. This together σ n (u) = σ t (u) = 0 implies that u satisfies conditions (5) (6) for any positive µ. In fact u is a solution which separates from the rigid foundation except at point B. We now consider the displacement field u. The displacements, the normal and tangential stresses on Γ C are: σ n (u) = 5 6 α, σ t(u) = 5 2 α, u n = 0, u t = 0. Obviously σ n (u) < 0 and σ t (u) µ σ n (u) when µ 3. As a consequence, the displacement field u satisfies conditions (5) (6) for any µ 3. Note that this solution is stuck on the rigid foundation. This concludes the proof of the proposition. Remarks. 1. The values µ 3 in the proposition correspond to quite large friction coefficients from an practical point of view. Nevertheless such a choice allows to exhibit a simple non-uniqueness example. In fact it is possible to obtain examples of non-uniqueness with µ 1 + ɛ for any ɛ > 0. This result is obtained together with the theoretical framework considering the general setting in which such nonunique isolated solutions occur. This work is in preparation. 2. The are probably some other types of nonuniqueness cases for problem (1) (6): a different approach from the one presented in this note consists of searching sufficient conditions of nonuniqueness involving infinitely many solutions (which all remain in slipping contact) for critical (eigen)values of the friction coefficient (see [6]). 3. The existence of examples with nonunique solutions to (1) (6) for arbitrary small friction coefficients is an open question. References [1] M. Cocu, E. Pratt and M. Raous, Formulation and approximation of quasistatic frictional contact, Int. J. Engng. Sci. 22 (1996) 783 798.

[2] C. A. de Coulomb, Théorie des machines simples, en ayant égard au frottement de leurs parties et la raideur des cordages. Pièce qui a remporté le Prix double de L Académie des Sciences pour l année 1781, Mémoire des Savants Etrangers, X (1785), 163-332. Reprinted by Bachelier, Paris 1809. [3] G. Duvaut, Problèmes unilatéraux en mécanique des milieux continus. Actes du Congrès International des Mathématiciens (Nice, 1970), Tome 3, pp. 71 77. Gauthier-Villars, Paris, 1971. [] G. Duvaut and J.L. Lions, Les inéquations en mécanique et en physique, Dunod, Paris, 1972. [5] C. Eck and J. Jarušek, Existence results for the static contact problem with Coulomb friction, Math. Mod. Meth. Appl. Sci. 8 (1998) 5 68. [6] R. Hassani, P. Hild and I. Ionescu, Sufficient condition of non-uniqueness for the Coulomb friction problem, Math. Meth. Appl. Sci. 27 (3) (200), to be published. [7] J. Jarušek, Contact problems with bounded friction. Coercive case, Czechoslovak. Math. J. 33 (1983) 237 261. [8] J. Nečas, J. Jarušek and J. Haslinger, On the solution of the variational inequality to the Signorini problem with small friction, Boll. Unione Mat. Ital. 17-B(5) (1980) 796 811. 5