Contact Mechanics of Adhesive Beams. Part I: Moderate Indentation. Arxiv:1707.07543V1
Total Page:16
File Type:pdf, Size:1020Kb
Contact mechanics of adhesive beams. Part I: Moderate indentation. V. S. Punati∗, I. Sharma and P. Wahi Mechanics & Applied Mathematics Group, Department of Mechanical Engineering, Indian Institute of Technology Kanpur, Kanpur, Uttar Pradesh, India - 208016. Abstract We investigate the contact of a rigid cylindrical punch with an adhesive beam mounted on flexible end supports. Adhesion is modeled through an adhesive zone model. The resulting Fredholm integral equation of the first kind is solved by a Galerkin projection method in terms of Chebyshev polynomials. Results are re- ported for several combinations of adhesive strengths, beam thickness, and support flexibilities characterized through torsional and vertical translational spring stiff- nesses. Special attention is paid to the important extreme cases of clamped and simply supported beams. The popular Johnson-Kendall-Roberts (JKR) model for adhesion is obtained as a limit of the adhesive zone model. Finally, we compare our predictions with preliminary experiments and also demonstrate the utility of our approach in modeling complex structural adhesives. Keywords: contact mechanics; adhesive beams; integral transforms. 1 Introduction Research in patterned adhesives is often motivated by the structures of natural adhe- sives, such as those present in the feet of gekkos; see e.g. Hiller (1976), and Arul and Ghatak (2008). In conventional adhesives, such as thin, sticky tapes, only the top and bottom surfaces are active. However, multiple surfaces may be activated with appropri- ate patterning. With more surfaces participating in the adhesion process the adhesives show increased hysteresis and, so, better performance. One example of a patterned ad- hesive is the structural adhesive shown in Fig. 1(a), which was developed by Arul and arXiv:1707.07543v1 [cond-mat.soft] 24 Jul 2017 Ghatak (2008). Figure 1(b) shows a possible mechanical model of the structural adhesive of Fig. 1(a) that utilizes several interacting adhesive beams. This motivates the goal of this paper, which is to investigate the adhesive contact of a beam; see Fig. 2(a). Contact with a half-space has been well studied over the past century, and we refer the reader to Alexandrov and Pozharskii (2001), and also the texts of Galin and Gladwell (2008), Gladwell (1980), Johnson (1985) and Hills et al. (1993). At the same time, the contact of thin layers is an active area of research in view of applications to electronics ∗Electronic address: [email protected]; [email protected] 1 P P Layer Punch Punch kt1 kt1 ks1 ks1 Micro-channels kt2 kt2 ks2 ks2 Rigid Platform Rigid Platform (a) (b) Figure 1: (a) Structural adhesive designed by Arul and Ghatak (2008). (b) Mechanical model of the structural adhesive in (a) employing an interconnected stack of adhesive beams. The rigidity of the vertical walls is modeled through torsional (stiffness kt) and vertical translational (stiffness ks) springs, as shown. The system is indented by a rigid punch, pressed down by the force P . P x kt h kt k k s a a s y l l σ0 (a) a δc c P x ∞ h ∞ h kt kt k s a a ks y l l (b) Figure 2: (a) Indentation by a rigid cylindrical punch of an adhesive beam resting upon flexible supports. The flexible supports are modeled through torsional and vertical trans- lational springs with stiffnesses kt and ks, respectively. (b) Mathematical model of the indentation process shown in (a). The inset shows details of the adhesive zone active near the contact edges; see text for details. The vertical deflection is exaggerated for ease of representation. 2 and computer industry; see, e.g. Barthel and Perriot ( 2007), and Dalmeya et al. (2012). In contrast, the indentation of beams is much less studied. Seeking the solution to the indentation of a beam through a strength-of-materials approach overlooks the local contact mechanics. To probe the latter, it is necessary to formulate an elasticity problem with appropriate boundary conditions. This is typically a complex problem, and has prompted some alternate approaches to adhesionless contact, as discussed in the next paragraph. To the best of our knowledge there is no available work on the adhesive contact of beams. Two-dimensional adhesionless indentation of a beam has been studied in the past by Keer and Silva (1970), Keer and Miller (1983), Keer and Schonberg (1986), Sankar and Sun (1983), Sun and Sankar (1985), and Kim et al. (2014). Keer and Miller (1983) modeled the beam as a linear elastic layer of infinite extent with frictionless bottom and top surfaces. First, the elasticity problem was solved thorugh Fourier transforms, see e.g. Sneddon (1995, p. 395-414). Then, employing the Hankel transform (Gladwell 1980, p. 213) for the pressure distribution, and asymptotically matching the far-field displacements with those obtained from Euler-Bernoulli beam theory, a Fredholm integral equation of the second kind was obtained. This equation was solved numerically. Sankar and Sun (1983) employed Fourier series in their investigation of adhesionless contact with beams of finite length. Their results were in good agreement with those of Keer and Miller (1983). Recently, Kim et al. (2014) studied beam indentation through asymptotics. Finite element (FE) simulations were also carried out. The contact pa- rameters, i.e. contact area and the total load acting on the punch, obtained through asymptotics, matched results of FE simulations well. In all these studies, the interaction of the beam with the punch is non-adhesive. However, extending these methods to adhe- sive beams is difficult due to the presence of several iterated integral transforms. Here, we propose an alternative approach to the contact of a beam with a rigid punch, shown in Fig. 2(a), that models both non-adhesive and adhesive interactions in a straightforward manner. Adhesion is modeled through the introduction of an adhesive zone that extends beyond the contact zone, as shown in the inset of Fig. 2(b). An adhesive zone model allows us to investigate the effect of adhesion by admitting continuous variations in it's strength. The popular JKR (Johnson et al., 1971) and DMT (Derjaguin et al., 1975) models for adhesion are obtained as special cases. The non-adhesive case (`Hertzian contact') is found by setting the adhesion strength to zero. This paper is organized as follows. We start by formulating a mathematical model for the adhesive contact of the beam shown in Fig. 2(a). We will obtain a Fredholm integral equation of the first kind that relates the contact pressure distribution with the vertical displacement in the contact patch. Relevant conditions at the edge of the contact zone will also be derived. This integral equation is then solved numerically employing Galerkin projections in terms of Chebyshev polynomials. Theoretical results are compared with FE simulations whenever possible. We then explore the effect of various parameters in the problem, viz. flexibility of end supports, strength of adhesion, the beam's geometry, etc. We close with a comparison with preliminary experimental results, and an application of our solution to complex structural adhesives, such as the one shown in Fig. 1(a). 3 2 Mathematical model We begin by extending the beam of Fig. 2(a) beyond the supports to infinity, as shown in Fig. 2(b). The extension is done in a manner consistent with the kinematic and kinetic constraints imposed by the supports. Thus, the beam is extended linearly along its slope at the supports. The beam may now be thought of as a linear elastic layer of infinite length with thickness h. The beam is isotropic and homogeneous, with Young's modulus E and Poisson's ratio ν. The top and bottom surfaces of the beam are frictionless. 1 During indentation, a normal traction distribution Pc(x) acts on the top surface . − At this time, the vertical displacement of the bottom surface is vb(x), which is typically not known. When the contact area is less than the beam's thickness, we assume that vb(x) may be approximated by the displacement obtained from Euler-Bernoulli beam theory when a point load of magnitude P acts at the center of the top surface, as depicted in Fig. 20. The details of how vb(x) is calculated for a beam on flexible supports are provided in App. A. The governing equations for the horizontal (u) and vertical (v) displacements in the extended beam (elastic layer) of Fig. 2(b), assuming plane strain, are given by 2 (1 ν) @2u @2u 1 @2v − + + = 0 (2.1a) 1 2ν @x2 @y2 1 2ν @x@y − − @2v 2 (1 ν) @2v 1 @2u and + − + = 0; (2.1b) @x2 1 2ν @y2 1 2ν @x@y − − which reflect horizontal and vertical linear momentum balance, respectively; see, e.g. Timoshenko and Goodier (1970, p. 241) or Sadd (2005, p. 125). The boundary conditions may be taken to be σxy = 0; σyy = Pc(x) on the top surface, i.e. at y = 0, (2.2a) − and σxy = 0; v = vb(x) on the bottom surface, i.e. at y = h: (2.2b) We now follow Sneddon (1995, p. 402) to map the above problem into Fourier space by transforming the x coordinate. Solving for the vertical displacement in Fourier space and, then, taking the inverse Fourier transform yields the following integral equation for the vertical displacement of the top surface: 1 Z 2 2 ¯ sinh ξ h v (x; 0) = Pc (ξ) cos ξx dξ − πE∗ ξ (ξ h + sinh ξ h cosh ξ h) 0 1 1 Z sinh ξ h + ξ h cosh ξ h + v¯ (ξ) cos ξx dξ; (2.3) π b ξ h + sinh ξ h cosh ξ h 0 where E∗ = E= (1 ν2), − 1 1 Z Z ¯ Pc (ξ) = Pc (t) cos ξt dt andv ¯b (ξ) = vb (t) cos ξt dt: (2.4) − −∞ −∞ 1The negative sign is introduced in order to report compressive pressure during contact as positive.