SUMMARYCorrect calculation of stresses at the interface of bonded or otherwise joined materials plays a significant role in many applications. It is therefore important that traction at the material interface is calculated as accurately as possible. This paper describes procedures that can be employed to achieve this goal by using centre‐based finite‐volume method. Total traction at the interface is calculated by decomposing it into normal and tangential components, both being calculated at each side of the interface, and applying the continuity assumption. The way in which the traction approximation is achieved depends on calculation of tangential gradient of displacement at the interface. To this end, three different methods are proposed and validated against problems with known solutions. It was shown that all methods can be successfully used to simulate problems with multi‐material domains, with the procedure based on finite area method being most accurate. Copyright © 2012 John Wiley & Sons, Ltd.