No presente trabalho implementou-se um modelo de estado plano de tensões usando o método de elementos finitos híbrido com funções de aproximação baseada em soluções fundamentais (fem-hfs). a implementação foi utilizada para investigar os resultados de deslocamento, os efeitos da continuidade dos campos de tensão e a convergência do método. a continuidade inter-elementar é imposta mediante o uso de um princípio variacional modificado incorporando um campo de interface independente entre os elementos. a formulação pode envolver base equilibrada regular ou singular, no entanto as matrizes são, em ambos os casos, obtidas por integração em regiões regulares. as matrizes do domínio do elemento podem ser convertidas em integrais de contorno devido às características da base auto-equilibrada das soluções fundamentais. é realizada a discretização do problema para aplicação numérica via fem-hfs e o código computacional foi implementado em linguagem fortran. a malha é composta por elementos triangulares lineares de três nós para definir os resultados de deslocamentos no domínio de integração e elementos triangulares de até nove nós ao longo de cada aresta, em grade regular, para definir os deslocamentos inter-elementares. por fim, o método é testado em problemas padrão com solução de referência disponível: são os casos de placa de material linear em regime elástico sob carregamento distribuído uniforme, de flexão para pequenas deformações e deslocamentos, além de problema de placa com concentração de tensões. os resultados são comparados com a literatura e software ansys.