Progress in Civil and Structural Engineering Computing
Chapter 9 A.A. Becker
School of Mechanical, Materials, Manufacturing Engineering and Management, University of Nottingham, United Kingdom Keywords: boundary element method, contact, friction, inverse methods, creep. The stress analysis of contact problems is a major concern in many engineering applications such as ball bearings, gears, rollers, mechanical seals and pressure vessel attachments. The numerical modelling of practical contact problems requires special attention because the actual contact area between the contacting bodies is usually unknown in advance and if friction is present, the behaviour may be dependent on load history [1]. The two widely used numerical techniques in the analysis of contact problems are the Finite Element (FE) and the Boundary Element (BE) methods. The BE method, with its surface-only modelling capability in linear problems offers several advantages over the FE method in solving contact problems [2]. The contact stresses calculated by the BE methods are more accurate because there is a lesser degree of approximation imposed on the solution, and the surface tractions, which are fundamental in determining the contact pressures, are calculated to the same degree of accuracy as the displacements, which enables accurate coupling of the contact variables. In BE formulations, since the contact variables are directly coupled to satisfy equilibrium and compatibility relationships, there is no need to use any form of additional gap or interface elements across the contact interface. Another advantage of using the BE method in contact problems is the ease of re- meshing a local contact surface, rather than modifying a volume mesh. In this paper, a review of BE formulations for contact problems is presented, where the contact variables are coupled directly to produce a set of simultaneous linear equations with a unique solution. Several types of contact interface conditions are covered including heat conduction, frictionless and frictional stick-slip. The review also covers three-dimensional contact problems, independent meshing of contact surfaces, non-linear material behaviour and inverse techniques. The contact applications demonstrate the versatility of the BE technique in modelling industrial problems such as pressure vessels and three-dimensional geometries. The contact surfaces can be discretised using a matching node-on-node mesh or an independent mesh which allows the user to specify a coarse mesh on one contact surface and a fine mesh on the other surface [3]. Thermoelastic stresses can be easily incorporated into contact problems by coupling the temperatures and heat flux between the contact nodes. In three-dimensional contact stick-slip problems, robust iterative schemes are needed to identify the correct direction of slip along the contact plane [4]. The contact formulation is extended to model contact under creep conditions in which the contact solutions have to be marched in time to arrive at the final creep time [5]. The inverse BE method determines the surface contact stresses from stress information obtained inside the solution domain, and is shown to be an ideal companion to experimental stress measurement techniques, such as photoelasticity, because the interior domain can be modelled with un-connected points without the need for full inferior discretisation [6]. All of the contact applications demonstrate the versatility and suitability of the BE method for contact problems, particularly for linear three-dimensional problems with rapidly changing stress gradients and inverse problems. BE contact formulations enable the direct coupling of the contact variables on the surface through the direct incorporation of the equilibrium, compatibility and frictional slip relationships into the overall system of equations without the need for gap or interface elements. The BE method provides an alternative accurate technique for contact analysis. References
|