This paper investigates the maximum stress concentrations in a ยฎnite strip loaded by a bonded elastic rivet by using the complex variable method in conjunction with the least-square boundary collocation method (BCM). The rivet-load is modeled by a uniform distributed body force; and the resultant rivet-force is acting along the transverse direction. The accuracy of the BCM is checked by comparing the results to those of the ยฎnite element method for a speciยฎc ยฎnite geometry of a strip and by the exact solution for the case of an inยฎnite plane. Numerical results show that the maximum shear and hoop stresses at the interface decrease with increasing b/R, where b is half of the width of the strip and R is the radius of the rivet. The maximum shear stress at the interface increases with z = m 2 /m 1 (where m 1 and m 2 are the shear moduli of the strip and rivet respectively) while the maximum hoop stress decreases with z. For ze1, the maximum normal bond stress at the interface decreases initially to a local minimum before rising to a steady value as b/R further increases. As b/R increases, the angular location of maximum stress occurrence y max , which is measured from the direction of resultant rivet-force, increases from about 368 H 428 to 908 (the inยฎnite plane limit) for the shear bond stress, and jumps suddenly from a roughly constant value (508 H 558) to 08 (the inยฎnite plane limit) for the normal bond stress. Similar sudden shifts in the angular location of maximum stress are also observed in the hoop stress at the interface.