Jtam.dvi JOURNAL OF THEORETICAL AND APPLIED MECHANICS 3, 39, 2001 CONTACT PROBLEMS WITH FRICTION, ADHESION AND WEAR IN ORTHOPAEDIC BIOMECHANICS. PART I – GENERAL DEVELOPMENTS Jerzy Rojek Józef Joachim Telega Institute of Fundamental Technological Research, Warsaw e-mail: jrojek@ippt.gov.pl, jtelega@ippt.gov.pl The present paper is a continuation of the contribution by Rojek and Telega (1999). An alternative adhesion law is used to the study of the bone-implant interface. Various problems related to the bone-implant interface are discussed. Key words: bone-implant interface, contact, adhesion, friction 1. Introduction Modelling of stress distribution and functioning of human (and animal) joints are inherently related to the study of complicated contact problems. One has to distinguish two major classes of such problems: the first class pertains to natural and diseased joints whilst the second class to artificial joints or joints after arthroplasty, cf Section 2 of this paper. The reader should be aware that understanding andmodelling of human joints is far from being complete. The aim of the present contribution is further development of the study of adhesion on bone-implant interface in human joints after arthroplasty, star- ted in Rojek and Telega (1999), cf also Bednarz et al. (2000). An important role played by adhesion on the bone-implant interface has been confirmed by inspection of retrieved implants, cf Rojek and Telega (1999), Laffargue et al. (1998). In the paper by Rojek and Telega (1999) the model of adhesion was essentially based on the idea due to Frémond (1987, 1988). Then the set of unilateral constraints is nonconvex. A modified approach was developed by Raous and his coworkers, cf Bretelle et al. (2000), Cocu and Roca (2000), 656 J.Rojek, J.J.Telega Cocu et al. (1997), Raous et al. (1999). Then one has to deal with two convex sets of constraints, one for a kinematic variable, and the second set for the density of adhesion β, cf Section 4.4. InRojek andTelega (1999), the influence of adhesion on the behaviour of a bone-femoral implantwas investigated. In thepresent paperwe shall study the tibial component of the condylar knee prosthesis, cf also Rojek et al. (1999), Bednarz et al. (2000). More realistic modelling of the interface bone-implant requires taking into account the influence of wear debris on the adhesion. This intriguing problem will be discussed in Part II (Rojek et al., 2001). However, the reader should be warned that a satisfactory model comprising friction, wear and adhesion and applicable to joints after arthroplasty is not available. The plan of the paper is as follows. In Section 2 and 3 various biome- chanical aspects of mainly knee joints after arthroplasty are briefly discussed. Theoretical modelling of adhesion is the subject of Section 4. 2. Knee joint: overview of selected biomechanical issues The aim of this section is to discuss concisely some basic problems related to knee joints, mainly after arthroplasty. The hip joint will also incidentally bementioned. Implants, and particularly the knee prostheses, have already a long histo- ry, sometimes even dramatic, see the comprehensive review papers by Jansons (1987), Prendergast (2000), Walker and Blunn (1997) and the relevant refe- rences cited in Rojek and Telega (1999). Knee models Hefzy and Grood (1988) and Hefzy and Cooke (1996) reviewed various static and dynamic models of the human knee joint, including the models developed to determine the forces in themuscles and ligaments during various activities, cf Bull and Amis (1998), Seireg and Arvikar (1973), Blankenvoort and Huiskes (1991), Toutoungi et al. (1997), Zavatsky and O’Connor (1994), Zatsiorsky (1998). Bao and Willems (1999) developed an involute-on-plane model of the knee. The parameters appearing in this model are the position of the involute basic circle and the radius of this circle. A method has been proposed to estimate these parameters. Contact problems with friction... 657 Determination of contact stresses Early results concerning the contact stresses between femur and tibia are due to Chand et al. (1976). These authors studied a natural knee joint, and the unilateral contact problem between the femur and tibia was solved as a quadratic programming problem. Moreover, the stress-freezing photoelastic technique was used to show an irregular stress distribution in the contact region. Blankenvoort et al. (1991) introduced the articular surface description which is based on a general theory for a thin elastic layer on a rigid foun- dation, cf also Blankenvoort and Huiskes (1996). Heegaard (1993) proposed a frictionless discretized model of patello- femoral joint. The general procedure is valid for a moderate slip between de- formable bodies or large slip contact with rigid obstacle, cf also Heegaard and Curnier (1993). Waldman and Bryant (1997) assessed the contribution of ultra-high mo- lecular weight polyethylene (UHMWPE) nonlinear viscoelastic properties in total knee replacement (TKR) from a developed mathematical rolling model. From such amodel, predictions of the effects of kinematic conditions on con- tact stress and rolling resistance were possible. A key assumption was made regarding the equivalence of rolling and sliding with adhesive friction in the presence of lubricants. The hysteretic nature of UHMWPE caused a signifi- cant portion of the rolling energy to become unrecoverable, which could be dissipated as heat. Localised temperature rises within the region of contact can result in increased creep, wear, and oxidation of UHMWPE. Thepaper in twoparts byLewis (1998) critically reviews the contact stress determination studies of cadaveric and artificial joints, with and without an associated TJR (total joint replacement), as well as physical systems that are deemed to represent idealized conditions at articular surfaces in TJRs. Analytical and numerical methods are typical for modern contact mechanics. The experimental methods discussed by this author include themethod using a pressure-sensitive film, resistive ink sensor, transducers and photoelasticity. Fuji film has been widely used in the studies of the contact mechanics of articular joints since 1980, cf Lewis (1998),Wu et al. (1998) and the references therein. The applications of Fuji Pressensor film range from the qualitative assessments of contact stress pattern to quantitative analysis of the contact pressure in articular joints using digital image techniques. Wu et al. (1998) studied the influence of inserting the Fuji Pressensor film into articular joints on the contactmechanics using the finite elementmethod.Thenumerical tests 658 J.Rojek, J.J.Telega showed that the contact pressure changes significantly when the Pressensor film is inserted into a joint. Harris et al. (1999) developed an easy, reproducible and, as the authors claim, reliable technique to evaluate TKA (total knee arthroplasty) systems and to determine the influence of the design parameters on the contact area. The system used is a computerized contact area and a pressuremeasurement system,K-Scan 4000 (Tekscan, SouthBoston, USA). This systemwas compa- redwith theFuji pressure-sensitive filmusing a custom testing jig designed for freedomofmovement. Contact areameasurements using this jig demonstrated the limitations of Fuji films in determining the contact areas. TheK-Scan sys- tem enabled the measurements of contact areas at different loads and flexion angles with the TKA system. Herzog et al. (1998) determined the in situ functional and material pro- perties of articular cartilage in an experimental model of joint injury and quantified the corresponding in situ joint contact mechanics. The cat ACL – transected knee was chosen as the experimental model (ACL – anterior cru- ciate ligament). Stewart et al. (1998) developed a simplifiedmethod to analyse the contact problems for the composite cushion knee joint replacement problems which can be readily incorporated into the design cycle. These include predictions of the contact radius and the maximum contact pressure as well as the detailed stress distributions within different layers and the interfaces. Experimental measurement of the contact area has also been conducted in order to verify the theoretical predictions. Parameters of knee design Sathasivam andWalker (1999) determined the effects of four design para- meterswhich vary to a large extent in the range of contemporary knee replace- ment designs. The parameters were the femoral frontal radius (30 or 70mm), the difference between the femoral and tibial frontal radii (2 or 10mm), the tibial sagittal radius (56 or 80mm), and the posterior-distal transition angle (−8or−20◦),which is the angle atwhich the small posterior arc of the sagittal profile transfer to the larger distal arc. Walker and Sathasivam (1999) claim that the ideal TKRwould be the one where the bearing surfaces are anatomical in shape and where the cruciate ligaments are preserved. However, for a total knee replacements without one or both cruciate ligaments, the ideal situation is where the bearing surfaces provide the same kinematic function as the cruciates, in particular, providing for the anterior and posterior translations of the femur on the tibia. In ad- Contact problems with friction... 659 dition, there should be a provision for internal-external rotation. Walker and Sathasivam (1999) investigated the feasibility of designing a knee replacement which would control and allow motion in this way. The scheme envisaged was to employ the usual type of partially conforming bearing surfaces for a fixed- bearing knee, or the closely-conforming surfaces of a mobile bearing knee, in combination with Guide Surfaces. A further requirement of the design was a sufficient femoral-tibial conformity for acceptable contact stresses. Effect of anisotropy on stress distribution Lewis et al. (1998) investigated the effect on the stresses in a construct comprising the tibial component and tibial insert of TKR anchored to the contiguous boneswith cement. A two-dimensionalmodelwith perfect bonding at the interfaces was studied. Two cases were considered. In the first case all the materials were assumed to have linear isotropic properties. Heterogeneity in the values of the Young and Poisson coefficients was taken into account. In the second case, the cortical bone was taken to have transversely isotropic elastic properties, with the representation of the elastic properties of all other materials remaining unchanged. The results obtained are not surprising, at least for us, since a strong effect of the representation of the elastic properties of material on the stresses was observed. Lewis et al. (1998) investigated also an elastomeric metacarpophalangeal joint implant in two cases. In the first case the elastomer was taken to have linear elastic properties whilst in the second case it was treated as a hyperelastic Mooney-Rivlin material. Again, the differences observed are obvious. Fatigue properties of HAPEX Ton That et al. (2000) investigated fatigue properties of HAPEX (hydroxyapatite-reinforced polyethylene composite) developed for bone repla- cement. The uniaxial tests comprised tension, compression and torsion to fa- ilure in order to determine the ultimate tensile strength. The biaxial fatigue tests of the samematerial were performed using various combinations of axial and torsional stresses. Micromovements Bynow it iswell established thatmicromovements of implants either stabi- lise, during thefirstpostoperative year, ormay contribute to aseptic loosening, cf the papers contained in Lewis and Galante (1985), Turner-Smith (1993). There are four regions of micromovement in relation to the development and 660 J.Rojek, J.J.Telega use of prostheses: (i) micromovement predicted by finite element and other numerical studies, (ii) strain-inducedmicromovement at the bone/implant in- terface in cadavers, (iii) strain-induced micromovement at the bone/implant interface, measured intra-operatively, post-operatively or both, and (iv) mi- cromovement observed in patients with the passage of time migration, cf the comments by Freeman in Turner-Smith (1993), pp. 93-98. Various measure- ment techniques ofmicromovement are also presented inTurner-Smith (1993), for instance, intra-operativemeasurement of tibial componentmicromovement in TKA. The paper by Berzins and Sumner (2000b) provides a useful review of some of the basic principles involved in performing in vitro tests of implant stability, cf also An andDraughn (2000). The technique described canmeasu- re both the inducible or recoverable motion and the nonrecoverable motion (often referred to as migration). 3. The bone-implant interface Micromovements in fact reflect the bone-implant interface behaviour. We observe that the cemented fixation of prosthesis is nothing else but the bone- implant contact where an implant consists of two phases: metal and UHMW- PE.Earlier ideas and results onmechanical, biological, chemical and electrical factors in tissue response to implants are summarized in theWorkshopReport edited by Lewis and Galante (1985), and in Boss et al. (1994). An important part of recent results on the bone-implant interface were summarized in the books edited by An and Draughn (2000), Helsen and Breme (1998). The bone-implant interface behaviour is far from being fully recognised and understood; also one lacks reliable phenomenological models. Anyway, adhesion, wear and friction play an important role in modelling the bone- implant interface. A model comprising all these phenomena would be rather complicated and is not yet available. Therefore, in the next section we will focus onmodelling the adhesion with friction. Primarily, however, we present some important ideas pertaining to the structure of the interface of interest. 3.1. Structure of interface Causes of tissue reactions around alloplastic materials are identified as cell mediated hypersensitivity to an implant component and tissue modification due to the presence of wear particles, debris and corrosion products of the Contact problems with friction... 661 prostheses used, cf Hildebrand andHornez (1998), Jallot et al. (1999b), Jones (1998), Thull (1998). Refractorymetals used for implants such as titanium, zirconium, niobium, tantalum and their alloys, as well as ceramic materials, are characterized by very low disintegration rates and justify the question of how the physioche- mical communication between material surfaces and the extracellular matrix actually occurs, cf Fig.1. Fig. 1. Overview of the complex reactions at the surface of passivatedmetal biomaterials. The initial adsorption pattern and conformational changes to constituents of the formed biofilm will be influenced by the physiochemical properties of the implant surface, after Thull (1998) In organic macromolecules, intramolecular and intermolecular bonds, to- gether with oxygen bridge bonds, may break down giving rise to structural or conformational changes, or both. Conformational changes may arise as a result of an exchange of charge carriers between the surface of the biomaterial and the biological macromolecules. In chronological order, the interaction bone-implant runs in the following steps: (i) The implant surface builds an interface according to the properties of the solid state and the surrounding liquid phase. (ii) Aprotein layer adsorbs and is structured in response to the physiochemi- cal properties of the surface in the equilibrium state. 662 J.Rojek, J.J.Telega (iii) Cell recognise the protein film and react. (iv) Tissue is structured according to the properties of the protein and cell layer on the surface. We observe that the equilibrium state of the material surface in the bio- logical environment is characterized by identical rates of adsorption and de- sorption of ions and constituents of, for instance, the extracellular matrix, cf Thull (1998). Themolecules making up an interface are different from those of the bulk phase because they are not surrounded by bulk phasemolecules. The bonding energy of the surface molecule is less than that associated with a bulk phase molecule, and the energy of the surfacemolecule is therefore higher than that of molecules in the bulk phase. The force in the surface that attempts to keep the surface area at aminimum is called the surface tension and is usually given in mJm−2. Indications in the literature show that the critical surface tension should be of the order of 60-120mJm−2 for good tissue adhesion, cf Fig.2. Fig. 2. Critical surface tension leading to adhesion or non-adhesion of the adjacent tissue, after Thull (1998) Remark 3.1. For a studyof thebehaviour of bone replacement-implant, such as hydroxyapatite, the reader is referred to Benhayoune et al. (1999), Jallot (1998), Jallot et al. (1999a), and the references therein. Contact problems with friction... 663 3.2. Mechanical testing of bone-implant interface The contact which will be studied in Section 4 is purely mechanical. Ho- wever, the material coefficients reflect macroscopic properties of interface. Let us now briefly discuss the relevant experimental techniques. Tests of interfacial bonding strengthhavebeenpursuedusingavariety of testmethods, implant geometries, and surface configurations. One type of test is a torque test, frequently involving screw-shaped implants in the rabbit tibia or femur, cf An and Draughan (2000), Berzins and Sumner (2000b), Nakamura and Nishoguchi (2000), and the references therein. Another typesof tests arepushoutandpullout tests, cfBerzins andSumner (2000a), Berzins et al. (1997). Themost common applications of pushout and pullout tests include testing for the effects of implantmaterial, surface texture, cross-sectional geometry, porosity, and surface composition in the context of cementless fixation by bone ingrowth or bone apposition to the implant.A few studies haveusedpushout tests to test either the cement-bone or cement-metal interfaces in cemented replacements. Berzins and Sumner (2000a) provided a practical guide to conducting these tests and the extensive list of citations. With both pushout and pullout tests, a load is applied to the implant via a device connected to the cross-head of the materials testing machine and a force-displacement curve is recorded, cf Fig.3. Fig. 3. A typical load-displacement curve from a pushout or pullout test, F – the maximum forces applied to the implant during the test, E′ – apparent shear stiffness (the slope of the load-displacement curve in its linear region); EA – energy adsorption to failure, after Berzins and Sumner (2000a) 664 J.Rojek, J.J.Telega Another methods for investigating interfacial strength is the tensile test, in which loads are applied in a direction normal to the tissue-biomaterial interface, cf Nakamura and Nishoguchi (2000). For a lot of other problems and methods related to mechanical testing of thebone-implant interface, the reader is referred toMcKoy et al. (2000), Davis et al. (2000), Dhert and Jansen (2000), Wang et al. (2000). According to Breme (1998), the microstructure of the implant interface is of great importance, cf also Jallot (1998), Jallot et al. (1999b). For instance, Young’s modulus value of the material must be lowered in order to decrease the stiffness of the implant in the direction of the bone. As can be seen from Fig.4a, a drastic change in Young’s modulus from the implant to the bone should be avoided because the difference in the elastic deformation can cau- se delamination stresses in the interface. In contrast, a gradual transition in Young’s modulus from the bone to the implant helps to avoid such a drastic change of stresses, cf Fig.4b. Fig. 4. Scheme of change in Young’s modulus at the bone-implant interface with a smooth and a porous implant surface, respectively, after Breme (1998) Contact problems with friction... 665 4. Modelling unilateral contact with adhesion and friction In the present section we shall introduce the description of adhesion with friction. Here we follow Raous et al. (1999). This approach constitutes a mo- dification of the theory primarily proposed by Frémond (1987, 1988), cf also Rojek and Telega (1999). 4.1. Adhesion: some basic notions Adhesionbetween two solidsmaybedue to ionic, covalent,metallic, hydro- gen or van derWaals forces; for more details the reader is referred toMaugis (1982), Maugis and Barquins (1980), Breme et al. (1998), Possart (1988). To cut these bonds and to separate the two solids 1 and 2 in contact on a unit area, energies γ1 and γ2 are needed to create the unit areas 1 and 2, whereas the excess energy γ12 (interfacial energy) is recovered. The quantity w= γ1+γ2−γ12 (4.1) is theDupré energy of adhesion or the thermodynamic work of adhesion. The separation never occurs as awhole, but by progression of a crack. Du- ring this propagation, the interface bonds are broken, elastic energy is released and irreversible work is dissipated at the crack tip. Work of adhesion is a useful quantity because it distinguishes the two states, contact and separation.Thiswork is done over a very small distance for vanderWaals forces, 99%of thework is achievedwhenthesurfaces arepulled1 nmapart. For other types of bonding, such as ionic and covalent, even smaller distances are involved. Thus, the precise shape of the force separation curve need not be known to understandmany phenomena. Indeed, the precise shape may even not be measurable because of the instability of the spontaneous jumping of smooth surfaces into contact. For other problems related to adhesion in biomechanics the reader is re- ferred to Rojek and Telega (1999), and the references therein. 4.2. Formulation of the initial–boundary value problem Let us consider two deformable bodies (Fig.5) occupying a domain Ω = Ω(1) ∪Ω(2) of R3 in their undeformed configuration. Let Γ(α) = ∂Ω(α) (α=1,2) consist of three disjoint parts: Γα0 , Γ α 1 and Γ α c such that Γ(α) =Γ α 0 ∪Γ α 1 ∪Γ α c (4.2) 666 J.Rojek, J.J.Telega Fig. 5. Two deformable bodies in contact Here the bar denotes the closure of set, and Γc = Γ (1) c = Γ (2) c is the contact surface of the two bodies Γc=Γ(1)∩Γ(2) (4.3) We assume that the considered bodies are initially in contact, that is Γc 6= ∅ (4.4) Moreover, the boundaries of the deformed bodies Γ(α), α = 1,2, possess a unique outward normal n(α) at each of their points. At the common part of the boundary, Γc, we have n (1) =−n(2) (4.5) Themotion of the two bodies is described by the following system of equ- ations σ (α) ij,j +f (α) i = ρü (α) i in Ω (α) u (α) i =u (α) i on Γ (α) 0 σ (α) ij n (α) j = g (α) i on Γ (α) 1 σ (α) ij = a (α) ijkleij(u (α)) (4.6) where α=1,2. All the static and kinematic quantities appearing in Eqs (4.6) depend on spatial variables and time t ∈ (0,T). Here and below there is no summation over α. The initial conditions for the displacement fields are specified by u(α)(0,x(α))=u (α) 0 (x (α)) (4.7) u̇(α)(0,x(α))=u (α) 1 (x (α)) x(α) ∈Ω(α) Contact problems with friction... 667 The elastic moduli aijkl are functions from L ∞(Ω) and satisfy the follo- wing condition ∃c1 ­ c0 > 0 ∀ǫ∈ E 3 s c0ǫijǫij ¬ aijkl(x)ǫijǫkl ¬ c1ǫijǫij (4.8) for almost every x∈Ω. Here E3s stands for the space of symmetric 3×3ma- trices.Obviously, the functions u(α),g(α),f(α),u (α) 0 , and u (α) 1 are prescribed. Our assumptions allow for the body to bemade of an anisotropic material. The set of equations (4.6) and (4.7) has to be supplemented with contact conditions on Γc. 4.3. Conditions for contact with adhesion We define the relative displacement [[u]] on the contact surface Γc as [[u]] =u(1)−u(2) (4.9) which can be decomposed into the normal and tangential components, un and uT , respectively un =unn uT = [[u]]−unn (4.10) where un = [[u]] ·n (4.11) and the unit normal vector n is taken as exterior to Ω1 n=n(1) =−n(2) (4.12) The interaction between the two bodies can be represented by the contact traction vectors, R(1) and R(2), which satisfy the following relations R (1) i =σ (1) ij n (1) j R (2) i =σ (2) ij n (2) j (4.13) By the Newton’s third law, we have R(1) =−R(2) (4.14) We take R=R(1) and similarly to (4.10) and (4.11), decompose R into the normal and tangential components, Rn and RT , respectively R=Rn+RT =Rnn+RT In the standard unilateral contact formulation no tensile normal contact forces are allowed Rn ¬ 0 (4.15) 668 J.Rojek, J.J.Telega In the present formulation contact forces can be either compressive or tensile (due to adhesion). The adhesion can generate reactions both in the normal and tangential directions. Possibility of decohesion (partial or total) is taken into account, the bond restitution will not be allowed, however. The formulation presented here employs a general framework for the study of contact problems with adhesion developed by Frémond (1987, 1988) and used in Rojek and Telega (1999), cf also Alves and Kikuchi (1998), Vena (1998). Some of the concepts presented by Raous et al. (1999) have also been adopted, cf also Cocu and Rocca (2000), Cocu et al. (1997), Bretelle et al. (2000).Themain idea consists in introducing the intensity of adhesion β(x, t), where x∈Γc and t stands for the time variable. The intensity of adhesion β is such that: (i) if β=1, all the bonds are active, (ii) if β=0, all the bonds are broken or the adhesion is absent, (iii) if 0<β < 1, a part β of the the bonds remains active, the remaining bonds are broken, the adhesion is partial. Initial adhesive bonds are given by β(0,x)=β0(x) x∈Γc (4.16) which supplements the set of equations (4.6) and (4.7). 4.4. Constitutive model of the contact interface Let us consider viscous contactwith adhesion and friction.We shall formu- late a constitutive model for such a contact interface allowing us to evaluate the contact reaction R. The contact traction will be split into reversible, Ren and ReT , and irreversible parts R i n and R i T (on Γc). As the state variables on the contact interface we take un,uT and β. Then the conjugate thermo- dynamic forces are reversible parts of the contact traction, Ren and R e T , as well as the thermodynamic force of decohesion Gβ. We define the generalized potential ϕ(β,un,uT) in the following form ϕ(β,un,uT)= kn 2 u2nβ 2+ kT 2 ‖uT‖ 2β2−wβ+ IP(β)+ IK(un) (4.17) where kn and kT are non-negative constants characterizing the interface stiffness, w stands for Dupré’s energy, IP and IK are indicator functions of the sets K and P K = {v | v­ 0}= R+ P = {γ | 0¬ γ¬ 1} (4.18) Contact problems with friction... 669 defined as, cf Rockafellar andWets (1998) IK(v) = { 0 if v∈K +∞ otherwise (4.19) IP(γ)= { 0 if γ ∈P +∞ otherwise Here R+ denotes the set of non-negative reals. The indicator functions impose appropriate constraints on un (impenetrability) and β. In Rojek and Telega (1999) only one indicator function is involved, nonconvex with respect to the couple (β,u), cf Frémond (1987, 1988). Here we follow a modified approach primarily used by Raous et al. (1999). The nonsmooth potential (4.17), so- metimes called a pseudo-potential (Panagiotopoulos, 1993), has a convenient feature, since it incorporates two convex indicator functions, IP and IK. The extended function ϕ is obviously nondifferentiable. Consequently, the state laws are written in the form of partial subdifferentials Ren ∈ ∂unϕ(β,un,uT) R e T ∈ ∂uTϕ(β,un,uT) (4.20) −Gβ ∈ ∂βϕ(β,un,uT) where, for instance, ∂unϕ denotes the subdifferential of the the potential ϕ with respect to the variable un. Subdifferentials (4.20) are defined as follows, cf Panagiotopoulos (1993), Rockafellar andWets (1998) ∂unϕ(β,un,uT)= {Qn |ϕ(β,vn,uT)−ϕ(β,un,uT)­Qn(vn−un) ∀vn} ∂uTϕ(β,un,uT)= {QT |ϕ(β,un,vT)−ϕ(β,un,uT)­QT · (vT −uT) ∀vT} (4.21) ∂βϕ(β,un,uT)= {Qβ |ϕ(η,un,vT)−ϕ(β,un,uT)­Qβ(η−β) ∀η} Explicit subdifferentiation (4.21) leads to the following relations, cf Raous et al. (1999) un ­ 0 −R e n+knunβ 2 ­ 0 (−Ren+knunβ 2)un =0 (4.22) ReT = kTuTβ 2 (4.23) 670 J.Rojek, J.J.Telega Gβ ­w if β=0 Gβ =w− (knu 2 n+kT‖uT‖ 2)β if β ∈]0,1[ Gβ ¬w− (knu 2 n+kT‖uT‖ 2) if β=1 (4.24) We observe that conditions (4.22) express generalized Signorini’s conditions for the contact with adhesion. To describe the irreversible parts Rin and R i T we introduce a positive function Φ, the so-called pseudo-potential of dissipation. Following Raous et al. (1999) we choose the following form for the extended function Φ Φ(Rn,un,β;u̇T , β̇)=µ|Rn−knunβ 2|‖u̇T‖+ b p+1 |β̇|p+1+ IIR−(β̇) (4.25) where µ is the Coulomb friction coefficient, whilst the parameters b and p characterize viscous properties of adhesive bonds and R− stands for the set of nonpositive reals. The indicator function IIR−(β̇) imposes the constraint β̇ ¬ 0, which means that adhesive bonds can only weaken and cannot be restituted. The first term on the right-hand side of Eq. (4.25) characterizes the dissipation due to friction, while the second term represents a viscous- type dissipation. The formof the friction dissipation inEq. (4.25) corresponds to the Coulomb law of friction; other friction models are possible. From the mathematical point of view, the space of functions with variation bounded in time is an appropriate space for βs, cf references in Rojek andTelega (1999). The complementary laws defining irreversible parts of contact traction Rin and RiT , and the thermodynamic force Gβ can be written as Rin =0 R i T ∈ ∂Φu̇T(Rn,un,β;u̇T , β̇) (4.26) Gβ ∈ ∂Φβ̇(Rn,un,β;u̇T , β̇) Differential inclusion (4.26)2 leads to the following condition ‖RiT‖¬µ|Rn−knunβ 2| (4.27) The partial subdifferential ∂Φ β̇ of the function Φ(Rn,un,β;u̇T , ·) is a sumof two subdifferentials, cf Rockafeller andWets (1998) ∂Φ β̇ (Rn,un,β;u̇T , β̇)= b|β̇| p∂|β̇|+∂IIR−(β̇) where ∂IIR−(β̇)=          0 if β̇ < 0 R − if β̇=0 ∅ otherwise Contact problems with friction... 671 Here ∅ denotes the empty set. Thus we finally get Gβ =−b|β̇| p if β̇ < 0 Gβ ∈ R − if β̇=0 (4.28) or β̇=− ( − G−β b )1/p since ∂|β̇|=          −1 if β̇ < 0 [−1,1] if β̇=0 +1 if β̇ > 0 Condition (4.27) describes two situations: (i) stick ‖RiT‖<µ|Rn−knunβ 2| ⇒ u̇T =0 (4.29) (ii) slip ‖RiT‖=µ|Rn−knunβ 2| ⇒ ∃λ­ 0 u̇T =−λR i T (4.30) The stick/slip conditions (4.29) and (4.30) can be written in the equivalent form φ¬ 0 λ­ 0 φλ=0 (4.31) u̇T +λR i T =0 where φ= ‖RiT‖−µ|Rn−knunβ 2| (4.32) Relation (4.28) can be combined with (4.24) to eliminate the thermodynamic force Gβ β̇=− { − 1 b [w− (knu 2 n+kT‖uT‖ 2)β]− }1/p if β ∈ [0,1[ β̇¬− { − 1 b [w− (knu 2 n+kT‖uT‖ 2)]− }1/p if β=1 (4.33) Summarizingwe can state that the complete set of contact conditions that must be taken into account consists of conditions for normal contact with adhesion (4.22), relations for tangential contact with adhesion and friction (4.23), (4.31), as well as the evolution law for the adhesion intensity β given by relations (4.33). 672 J.Rojek, J.J.Telega Remark 4.1. In the sense of the theory of hemivariational inequalities (Pa- nagiotopoulos, 1993), the interface laws involving adhesionare of nonmo- notone type, cf Mistakidis and Panagiotopoulos (1997), Tzaferopoulos andPanagiotopoulos (1994). In other words, the interface laws could be formulated in the form of a hemivariational inequality. Licht and Michaille (1997) studied a rather general case of bonding of two nonlinear elastic bodies, separated by a thin adhesive layer, cf also the references cited therein. Acknowledgment Theauthorswere supportedby theStateCommittee forScientificResearch(KBN, Poland) through the grant No. 8T11F01718. References 1. Alves M.K., Kikuchi N., 1998, Damage Theory Applied to the Analysis of Composite Materials Subjected to a Loss of Adherence Between its Constitu- ents,Eur. J. Mech., A/Solids, 17, 53-69 2. An Y.H., Draughn R.A., Edit., 2000, Mechanical Testing of Bone and the Bone-Implant Interface, CRCPress, Boca Raton 3. BaoH.,WillemsP.Y., 1999,On theKinematicModeling and theParameter Estimation of the HumanKnee Joint, J. Biomech. Eng., 121, 406-413 4. BednarzP.,RojekJ.,TelegaJ.J.,Kowalczewski J., 2000,ContactPro- blems inKnee Joint after Arthroplasty,Acta Bioeng. Biomech., 1, Supplement 1, 73-78 5. Benhayoune H., Jallot E., Frayssinet P., Balossier G., Bonhomme P., 1999,Etudede laResorptiondeL’Hydroxyapatite dense après Implementa- tion parMicroscopie Electronique à Balayage etMicroanalyse X, Innov. Tech. Biol. Med., 20, 181-186 6. Berzins A., Sumner D.R., 2000a, Implant Pushout and Pullout Tests, In: Y.H. An and R.A. Draughn, Edit.,Mechanical Testing of Bone and the Bone- Implant Interface, 463-476, CRCPress, Boca Raton 7. Berzins A., Sumner D.R., 2000b, In Vitro Measurements of Implant Stabi- lity, In: Y.H. An andR.A.Draughn, Edit.,Mechanical Testing of Bone and the Bone-Implant Interface, 515-526, CRCPress, Boca Raton Contact problems with friction... 673 8. Berzins A., Shah B., Weinans H., Sumner D.R., 1997, Nondestructive Measurements of Implant-Bone Interface ShearModulus andEffects of Implant Geometry in Pull-out Tests, J. Biomed. Mater. Res., 34, 337-340 9. BlankenvoortL.,HuiskesR., 1991,Ligament-Bone Interaction in aThree- Dimensional Model of the Knee, J. Biomech. Eng., 113, 263-269 10. Blankenvoort L., Huiskes R., 1996, Validation of a Three-Dimensional Model of the Knee, J. Biomechanics, 29, 955-961 11. Blankenvoort L., Kuiper J.H., Huiskes R., Grootenboer H.J., 1991, ArticularContact in aThree-DimensionalModel of theKnee,J. Biomechanics, 24, 1019-1031 12. Boss J.H., Shajrawi I., Mendes D.G., 1994, The Nature of Bone-Implant Interface,Med. Prog. Technol., 20, 119-142 13. Breme H.J., 1998, Surfaces, SurfaceModification and Tailoring, In: J.A. Hel- sen and H.J. Breme, Edit., Metals as Biomaterials, 153-175, John Wiley & Sons, Chichester 14. Breme H.J., Barbosa M.A., Rocha L.A., 1998, Adhesion to Ceramics, In: J.A. Helsen and H.J. Breme, Edit., Metals as Biomaterials, 219-264, John Wiley & Sons, Chichester 15. BretelleA.S., CocuM.,MonerieY., 2000,Formulation duContactAvec Adhérence en Elasticité Non Linéaire Entre Deux Solides Déformables, C.R. Acad. Sci. Paris, Série IIb, 328, 203-208 16. Bull A.M., Amis A.A., 1998, Knee JointMotion: Description andMeasure- ment,Proc. Inst. Mech. Engrs., 212, Part H, 357-372 17. Chand R., Haug E., Rim K., 1976, Stress in the Human Knee Joint, J. Biomechanics, 9, 417-422 18. Cocu M., Cangemi L., Raous M., 1997, ApproximationResults for a Class of Quasistatic Contact Problems Including Adhesion and Friction, In: Proc. IUTAM Symp. ‘Variations of Domain and Free Boundary Problems in Solid Mechanics’, Paris, April 22-25, 1997 19. Cocu M., Rocca R., 2000, Existence of Results for Unilateral Contact Pro- blemswithFrictionandAdhesion,Math.Modelling Numer.Anal.,34, 981-1001 20. Davis J.R., Lofthouse R., Jinnah R.H., 2000, In Vitro Testing of the Stability of Acetabular Components, In: Y.H. An and R.A. Draughn, Edit., Mechanical Testing of Bone and the Bone-Implant Interface, 527-540, CRC Press, Boca Raton 21. Dhert W., Jansen J., 2000, The Validity of a Single Pushout Test, In: Y.H. AnandR.A.Draughn,Edit.,Mechanical Testing ofBone and theBone-Implant Interface, 477-488, CRCPress, Boca Raton 674 J.Rojek, J.J.Telega 22. Frémond M., 1987, Adhérence des Solides, J. Méc. Théorique Appl., 6, 383- 407 23. Frémond M., 1988, Contact with Adhesion, In: J. Moreau and P.D. Pana- giotopoulos, Edit.,Nonsmooth Mechanics and Applications, 177-221, Springer- Verlag,Wien 24. Harris M.L., Morberg P., Bruce W.J.M., Walsh W.R., 1999, An Im- provedMethod forMeasuringTibiofemoralContactAreas inTotalKneeArth- roplasty: AComparison of K-Scan Sensor and Fuji Film, J. Biomechanics, 32, 951-958 25. Heegaard J.-H., 1993, Large Slip Contact in Biomechanics: Kinematics and Stress Analysis of the Patello-Femoral Joint, Thèse No. 1113, Ecole Polytech- nique Fédérale de Lausanne 26. Heegaard J.-H., Curnier A., 1993, AnAugmented LagrangianMethod for Discrete Large-Slip Contact Problems, Int. J. Num. Meth. Eng., 36, 569-593 27. Hefzy M.S., Cooke T.D.V., 1996, Review of Knee Models: 1996 Update, Appl. Mech. Reviews, 49, Part 2, S187-S193 28. Hefzy M.S., GroodE., 1988, Review ofKneeModels,Appl. Mech. Reviews, 41, 1-13 29. Helsen J.A., Breme H., Edit., 1998,Metals as Biomaterials, JohnWiley & Sons, Chichester 30. HerzogW., Diet S., Suter E.,Mayzus P., LeonardT.R.,MüllerC., Wu J.Z., Epstein M., 1998,Material and Functional Properties of Articular Cartilage and Patellofemoral ContactMechanics in an Experimental Model of Osteoarthritis, J. Biomechanics, 31, 1137-1145 31. Hildebrand H.F., Hornez J.-C., 1998, Biological Response and Biocompa- tibility, In: J.A.Helsen andH.J. Breme,Edit.,Metals as Biomaterials, 265-290, JohnWiley & Sons, Chichester 32. Jallot E., 1998, Correlation Between Hydroxyapatite Osseointegration and Young’sModulus,Medical Eng. Physics, 20, 697-701 33. JallotE., Irigaray J.L.,WeberG., Frayssinet P., 1999a, InVivoCha- racterization of the Interface Between Cortical Bone and Biphasic Calcium Phosphate by the PIXEMethod, Surface and Interface Annals, 27, 648-652 34. Jallot E., Benhayoune H., Kilian L., Balossier G., Bonhomme P., Oudadesse H., Irigaray J.L., 1999b, Ultrastructural Characterization of Soft Tissues Surrounding Implanted Hip Prosthesis by Complementary PIXE and TEMMethods.Nuclear Instruments Meth. Phys. Res.,B155, 315-321 35. Jansons H.A., 1987, Development of Endoprosthesis in the U.S.S.R. CRC Crictical Reviews in Biocompatibility, 3, 87-192 Contact problems with friction... 675 36. Jones D.B., 1998, Cells and Metals, In: J.A. Helsen and H.J. Breme, Edit., Metals as Biomaterials, 317-335, JohnWiley & Sons, Chichester 37. Laffargue P., Breme H.J., Helsen J.A., Hildebrand H.F., 1998, Re- trievalAnalysis, In: J.A.Helsen andH.J. Breme,Edit.,Metals as Biomaterials, 467-501, JohnWiley & Sons, Chichester 38. Lewis G., 1998, Contact Stress at Articular Surfaces in Total Joint Repla- cements. Part I: Experimental Methods, Bio-Medical Mat. Eng., 8, 91-110; Part II: Analytical and Numerical Methods, ibidem, 259-278 39. LewisG., JonathanV.,Kambhampati S., 1998,EffectofMaterialProperty Representation on Stress in Endoprostheses,Bio-Medical Mat. Eng., 8, 11-23 40. Lewis J.L., Galante J.O., Edit., 1985,The Bone-Implant Interface, Work- shop Report, American Academy of Orthopaedic Surgeons 41. Licht C., Michaille G., 1997, A Modelling of Elastic Adhesive Bonded Jo- ints,Adv. Math. Sci. Appl., 7, 711-740 42. MaugisD., 1982,Adherence ofSolids, In: J.Georges,Ed.,Microscopic Aspects of Adhesion and Lubrication, 221-252, Elsevier, Amsterdam 43. Maugis D., Barquins M., 1980, Fracture Mechanics and Adherence of Vi- scoplastic Solids, In: H.-L. Lee, Edit., Adhesion and Adsorption of Polymers, Part A, 203-277, PlenumPublishing Corporation, NewYork 44. McKoy B., An Y.H., Friedman R., 2000, Factors Affecting the Strength of the Bone-Implant Interface, In: Y.H. An andR.A. Draughn, Edit.,Mechanical Testing of Bone and the Bone-Implant Interface, 439-462, CRC Press, Boca Raton 45. Mistakidis E.S., Panagiotopoulos P.D., 1997, Numerical Treatment of Problems Involving Nonmonotone Boundary or Stress-Strain Laws, Comput. Struct., 64, 553-565 46. Nakamura T., Nishiguchi S., 2000, Tensile Testing of Bone-Implant Inter- face, In: Y.H. An and R.A. Draughn, Edit., Mechanical Testing of Bone and the Bone-Implant Interface, 489-497, CRCPress, Boca Raton 47. Panagiotopoulos P.D., 1993,Hemivariational Inequalities: Applications in Mechanics and Engineering, Springer-Verlag, Berlin 48. Possart W., 1998, Adhesion of Polymers, In: J.A. Helsen and H.J. Breme, Edit.,Metals as Biomaterials, 197-218, JohnWiley & Sons, Chichester 49. PrendergastP.J., 2000,BoneProstheses and Implants, In: S.C.Cowin,Ed., Bone Mechanics, (35-1)-(35-29), CRCPress, Boca Raton 50. Raous M., Cangemi L., Cocu M., 1999, A ConsistentModel Coupling Ad- hesion, Frictional Unilateral Contact, Comput. Meth. Appl. Mech. Eng., 177, 383-399 676 J.Rojek, J.J.Telega 51. Rockafellar R., Wets R.-B., 1998,Variational Analysis, Springer, Berlin 52. Rojek J., Telega J.J., 1999,Numerical Simulation ofBone-ImplantSystems Using a More Realistic Model of Contact Interfaces with Adhesion, J. Theor. Appl. Mech., 37, 659-686 53. Rojek J., Telega J.J., Bednarz P., 1999, Adhesion in Hip and Knee Im- plants: Modelling and Numerical Analysis, Acta Bioeng. Biomech., 1, Supple- ment 1, 377-380 54. Rojek J., Telega J.J., Stupkiewicz S., 2001,Contact Problemswith Fric- tion, Adhesion and Wear in Orthopaedic Biomechanics. Part II: Numerical Implementation and Application to Implanted Knee Joints. J. Theor. Appl. Mech., 39, 2001 55. Sathasivam S.,Walker P.S., 1999, TheConflictingRequirements of Laxity and Conformity in Total Knee Replacements, J. Biomechanics, 32, 239-247 56. Seireg A., Arvikar R.J., 1973, A Mathematical Model for Evaluation of Forces in Lower Extremities of theMusculo-Skeletal System, J. Biomechanics, 6, 313-326 57. Stewart T., Jim Z.M., Fisher J., 1998, Analysis of ContactMechanics for Composite Cushion Knee Joint Replacement, Proc. Inst. Mech. Engrs., 212, Part H, 1-10 58. Thull R., 1998, Tissue-Implant Interaction, In J.A. Helsen and H.J. Breme, Edit.,Metals as Biomaterials, 291-315, JohnWiley & Sons, Chichester 59. TonThatP.T., TannerK., BonfieldW., 2000,FatigueCharacterization of a Hydroxyapatite-Reinforced Polyethylene Composite. I. Uniaxial Fatigue, J. Biomed. Mat. Res., 51, 453-460; II. Biaxial Fatigue, ibidem, 461-468 60. Toutoungi D.E., ZavatskyA.B., O’Connor J.J., 1997, Parameter Sensi- tivity of aMathematical Model of the Anterior Cruciate Ligament,Proc. Inst. Mech. Engrs., 211, Part H, 235-246 61. Turner-SmithA., Ed., 1993,Micromovement in Orthopaedics, OxfordMedi- cal Engineering Series, 10, Clarendon Press, Oxford 62. TzaferopoulosM.A.,PanagiotopoulosP.D., 1994,ANumericalMethod for a Class of Hemivariational Inequalities,Comput. Mechanics, 15, 233-248 63. Vena P., 1998, Shape Optimization of a Femoral Head Endoprosthesis Allo- wing for Interface Behaviour, PhDThesis, Politecnico di Milano 64. Waldman S.D., Bryant J.T., 1997, Dynamic Contact Stress and Rolling Resistance Model for Total Knee Arthroplasties, J. Biomech. Eng., 119, 254- 260 65. Walker P.S., Blunn G.W., 1997, Biomechanical Principles of Total Knee Replacement Design, In: V.C. Mow andW.C. Hayes, Edit.,Basic Orthopaedic Biomechanics, 461-493, Lippincott-Raven Publishers, Philadelphia Contact problems with friction... 677 66. WalkerP.S., Sathasivam S., 1999,TheDesign ofGuide Surfaces for Fixed- Bearing andMobile-BearingKnee Replacements, J. Biomechanics, 32, 27-34 67. Wang X., Athanasiou K.A., Agrawal C.M., 2000, Fracture Toughness Tests of the Bone-Implant Interface, In: Y.H. An and R.A. Draughn, Edit., Mechanical Testing of Bone and the Bone-Implant Interface, 499-514, CRC Press, Boca Raton 68. Wu J.Z., Herzog W., Epstein M., 1998, Effects of Inserting a Pressensor Film intoArticular Joints on theActual ContactMechanics, J. Biomech. Eng., 120, 655-659 69. Zatsiorsky V.M., 1998,Kinematics of HumanMotion, HumanKinetics, Le- eds 70. Zavatsky A.B., O’Connor J.J., 1994, Three-Dimensional Geometrical Mo- dels of HumanKnee Ligaments,Proc. Inst. Mech. Engrs., 208, 229-240 Zagadnienia kontaktowe z tarciem, adhezją i zużyciem w biomechanice ortopedycznej. Część I – Rozważania ogólne Streszczenie Niniejsza praca stanowikontynuacjęwcześniejszej pracy autorów,por.Rojek iTe- lega (1999).Zastosowanoalternatywnymodel adhezji dla opsu interfazykość-implant. Przedyskutowano również szereg zagadnień związanych z tą interfazą. Manuscript received February 2, 2001; accepted for print May 23, 2001