The geometries of many problems of practical interest are created from circular or elliptic arcs. Arc boundary elements can represent these boundaries exactly, and consequently, errors caused by representing such geometries using polynomial shape functions can be removed. To fully utilize the geometry of circular boundary, the non-singular boundary integral equations (BIEs) and a general nonlinear transformation technique available for arc elements are introduced to remove or damp out the singular or nearly singular properties of the integral kernels. Several benchmark 2D elastostatic problems demonstrate that the present algorithm can effectively handle singular and nearly singular integrals occurring in the boundary element method (BEM) for boundary layer effect and thin-walled structural problems. Owing to the employment of exact geometrical representation, only a small number of elements need to be divided along the boundary and high accuracy can be achieved without increasing other more computational efforts.