Esta pesquisa visa o desenvolvimento de uma ferramenta numérica capaz de simular a inclusão de elementos drenantes em solos, sejam eles representados por drenos subhorizontais ou poços, constituindo sistemas de drenagem amplamente utilizado em estabilidade de encostas, túneis, escavações e outros. A implementação foi gerada nos programas de fluxo SWMS_2D e SWMS_3D (Simunek e outros, 1994). A formulação proposta considera a equação de fluxo do elemento drenante e a estratégia numérica de sua inclusão na equação de fluxo do solo pelo método dos elementos finitos. Algumas análises numéricas foram realizadas visando a validação do algoritmo. Para os poços foram analisados os casos confinados e não confinados e em regime permanente e transiente, comparando as soluções numéricas obtidas com as soluções analíticas de Theis (1935) (Freeze, 1979), para aqüífero confinado, e de Neumann (1975), para aqüífero não confinado. Para os drenos subhorizontais, propõe - se aqui uma metodologia de análise, levando em consideração os parâmetros hidráulicos e geométricos de um elemento de dreno. Sua aplicação atual não requer muito rigor e este estudo vem com a proposta de ser uma ferramenta geotécnica na fase de projeto de uma obra.