## 1 Introduction  
This paper presents an alternative solution to overcome the drawbacks of the boundary element method (BEM). The conventional BEM solves the continuum mechanics problems by using the reciprocal work theorem also known as Betti's theorem, the Navier-Cauchy equations, and the divergence theorem. The divergence theorem simplifies the continuum mechanics equations solution by modeling only at the boundaries [1]. Such modeling consists of dividing the boundary of the body under analysis into a discrete set of functions, denominated \"elements\" with nodes, and describing the displacements and tractions by polynomial functions. Since BEM is a consequence of the Betti's theorem, there are two sets of displacements and tractions. One set is known in advance through entities called traction and displacement kernels, which are valid for any equilibrium geometry. The other set belongs to the body under analysis which is the problem to be solved. Naturally, to obtain a unique solution, boundary conditions such as prescribed tractions and displacements must be applied. The calculated traction and displacement kernels depend on the normal of the surface, the geometric variable, and the interpolation points that describe the boundary. However, the discretization process produces a lack of continuity on the boundary. An alternative to solve this problem is to use the Non-Uniform Rational B-Splines functions (NURBS) from the CAD model instead of using interpolation functions. This concept was called Isogeometric Analysis (IGA), and it was proposed by [PERSON] et al. [2] in 2005. Initially, the IGA was used with the FEM framework, but eventually, it was combined with BEM and with the finite cell method (FD). [PERSON] et al. [3, 4] proposed the Isogeometric Boundary Element Method (IGA-BEM) in 2012 and established the guidelines for combining IGA with BEM, e.g., how the elements are generatedusing the NURBS knot vector and the NURBS basis functions. They also used the Greville abscissae scheme to get the collocation points and the displacement and traction kernels. Finally, they applied the method to solve the L-plate, the hole-within-an-infinite-plate, and the L-shaped wedge problems. [PERSON] and [PERSON] [5] proposed an IGA-BEM combined with the fast multipole method (FMM), to decrease the control points of the NURBS. They proved that this combination could achieve the same accuracy as the conventional BEM but with fewer degrees of freedom and less computing time. [PERSON] et al. [6] presented an extended IGA-BEM using a partition-of-unity method on NURBS functions; they applied it to solve a two-dimensional Helmholtz problem. Likewise, [PERSON] et al. [7] refined the IGA-BEM, using T-Splines instead of NURBS, to solve linear elastostatic problems; they showed the performance of this method with a patch test and a propeller analysis. Moreover, [PERSON] et al. [8] continued the work of [3,4,7], and they used the T-Splines with IGA-BEM, to optimize the shape of three-dimensional elastic bodies; they chose the control points as design variables to modify the geometry. The IGA-BEM combination was not only used to reduce the degrees of freedom, changing basis functions, or solving linear problems. For example, [PERSON] and [PERSON] [9] developed a method to calculate the singular integrals on a 3D potential problem and [PERSON] et al. [10] used IGA-BEM to analyze 3D Stock flows. [PERSON] et al. [11] used a geometric algorithm to propagate fractures, based on the fatigue Paris law, and IGA-BEM to analyze them. [PERSON] and [PERSON] [12] used IGA-BEM to solve the acoustic scattering problem of a submarine. The applications of the IGA-BEM in a wide variety of problems have given good results. Nevertheless, the solution depends on calculating the kernels at the interpolation points, restricting the solution domain to the geometry discretization domain. In this work, a methodology based on IGA-BEM is proposed.