Neste trabalho apresenta-se uma metodologia pelo método dos elementos finitos para a análise estática-dinâmica bi-dimensional de linhas marítimas considerando-se o contato com o leito marinho. Para o cabo considera-se o elemento de dois nós do modelo de viga de Euler-Bernoulli com as grandezas estáticas e dinâmicas referidas a um sistema co- rotacionado. Considera-se a nãolinearidade geométrica associada ao movimento de corpo rígido da estrutura, supondo pequenas deformações. No movimento da linha, são consideradas as influências do peso próprio, empuxo, carregamentos hidrodinâmicos das correntes marinhas, deslocamentos prescritos no ponto superior da linha, forças de inercia e ação de flutuadores. Os efeitos sobre a linha associados ao contato com o fundo marinho são considerados empregando-se a técnica dos multiplicadores de Lagrange. Estes representam o carregamento (forças e momentos) requerido para garantir as restrições geométricas. No modelo não são considerados os efeitos das forças de atrito com o fundo marinho, cuja rigidez é considerada muito superior à do cabo. A matriz de rigidez usual do modelo de vigas é aumentada de forma a acomodar as novas incógnitas resultantes da imposição das condições de contato. Desta forma, condições de contato, separação repetidos (contato variável) e grandes movimentos relativos do cabo com o leito são contemplados pelo modelo. A solução das equações algébricas não-lineares, resultantes da integração temporal passo-a-passo de Newmark, é obtida com a técnica iterativa de Newton Raphson. A metodologia numérica foi implementada e resultados de alguns testes para situações da engenharia offshore são mostrados e comparados com os apresentados por outros autores independentes.