Abstract This paper shows and discusses a generic implementation of the global-local analysis toward generalized finite element method (GFEMgl). This implementation, performed into an academic computational platform, follows the object-oriented approach presented by the authors in a previous work for the standard version of GFEM in which the shape functions of finite elements are hierarchically enriched by analytical functions, according to the problem behavior. In global-local GFEM, however, the enrichment functions are constructed numerically from the solution of a local problem. This strategy allows the use of a coarse mesh even when the problem produces complex stress distributions. On the other hand, a local problem is defined where the stress field presents high gradients and it is discretized using a large number of elements. The results of the local problem are used to enrich the global problem which improves the approximate solution. The great advantage is allowing a well-refined description of the local problem, when necessary, avoiding an overburden for the computation of the global solution. Details of the implementation are presented and important aspects of using this strategy are highlighted in the numerical examples.