Neste trabalho um novo elemento de viga co-rotacionado é apresentado para a análise não-linear geométrica tridimensional, estática e dinâmica, de linhas marítimas de Materiais com Gradação Funcional (MGF). Assume-se que o módulo de elasticidade e a massa específica do material da viga variam ao longo da espessura da seção transversal tubular de acordo com uma lei de potência. Na discretização espacial das equações de equilíbrio, a linha marítima é representada por um elemento de viga de dois nós, com base nas hipóteses do modelo para vigas de Euler-Bernoulli, em que polinômios cúbicos de Hermite são utilizados na interpolação dos deslocamentos nodais e a cinemática do movimento é descrita através de grandezas referidas a um sistema coordenado local co-rotacionado. Consideram-se não linearidades geométricas envolvendo grandes deslocamentos e rotações, mas com pequenas deformações. Nas equações de movimento da linha marítima, são consideradas as seguintes influencias: do peso próprio, do empuxo, dos carregamentos hidrodinâmicos (devidos às ações de ondas, correntes e forças de massa adicional), dos deslocamentos prescritos (junto à fixação da embarcação), da ação de flutuadores e das forças de interação solo-estrutura. A integração temporal das equações de equilíbrio é realizada utilizando-se o algoritmo de discretização HHT (Hilbert-Hughes-Taylor) e a solução numérica obtida com a técnica iterativa de Newton Raphson. A metodologia numérica foi implementada e diversos exemplos são apresentados e discutidos enfatizando-se as diferenças de comportamento estrutural entre os modelos de viga com MGF e com material homogêneo. Resultados referentes a situações práticas da engenharia offshore são também tratados nos exemplos.