The paper is concerned with the numerical solution of a thermoelastic rolling contact problem with wear. The friction between the bodies is governed by Coulomb law. A frictional heat generation and heat transfer across the contact surface as well as Archard's law of wear in contact zone are assumed. The friction coefficient is assumed to depend on temperature. In the paper quasistatic approach to solve this contact problem is employed. This approach is based on the assumption that for the observer moving with the rolling body the displacement of the supporting foundation is independent on time. The original thermoelastic contact problem described by the hyperbolic inequality governing the displacement and the parabolic equation governing the heat flow is transformed into elliptic inequality and elliptic equation, respectively. In order to solve numerically this system we decouple it into mechanical and thermal parts. Finite element method is used as a discretization method. Numerical examples showing the influence of the temperature dependent friction coefficient on the temperature distribution and the length of the contact zone are provided.
Contact Problem FormulationConsider deformations of an elastic strip lying on a rigid foundation (see Fig. 1). The strip occupies domain Ω ∈ R 2 with the boundary Γ. A wheel rolls along the upper surface Γ C of the strip with velocity V . Denote by u = (u 1 , u 2 ), u = u(x, t), x ∈ Ω, t ∈ (0, T ), T > 0, a displacement of the strip and by θ = θ(x, t) the absolute temperature of the strip. The displacement u and the temperature θ of the strip satisfiy the system of equations [2,3,5,9,10] with the following initial conditions:∂n (x, t) = q on Γ × (0, T ). Functionsū 0 ,ū 1 ,θ, q are given. Constants ρ, α,κ, c p denote mass density, thermal expansion, thermal conductivity, heat capacity coefficient of the strip material, respectively. λ and γ are Lame coefficients, Γ 0 = Γ \ Γ C , n = (n 1 , n 2 ) is outward normal versor to the boundary Γ. The operators A, B and D are defined as [1,7] A denotes a transpose of A. By F = (F 1 , F 2 ) we denote surface traction vector on the boundary Γ. This vector is a priori unknown and is given by conditions of contact and friction. These conditions on Γ C × (0, T ) take a form [1,2,4,6,7,9,10] Friction coefficient µ = µ(θ) is dependent on temperature θ. g r denotes the gap between bodies in contact. w = w(x, t) denotes an internal state variable to model the distance between the bodies due to wear [9] and satisfies the Archard law
Friction Coefficient Dependent on TemperatureUsing results of experimental and analytical investigation in [8], we employ the following Lipschitz continuous temperature dependence of friction coefficient (see Fig. 2):where µ 0 , 0 < µ F < µ 0 and 0 < θ 0 < θ F are given positive constants.