No presente trabalho se estuda um elemento finito subparamétrico que aproxima o campo de deslocamentos com funções spline, implementando um programa que pode ser utilizado para calculo estático, dinâmico e de instabilidade de estruturas compostas de placas, vigas de paredes finas, vigas caixão e em geral em elementos alongados (pontes e perfis metálicos). O grau de liberdade de rotação perpendicular ao plano do elemento é introduzido na formulação para possibilitar uma análise tridimensional. Apresenta-se um método que serve como base para determinar a constante de rigidez correspondente. Nos exemplos apresentados avalia-se a precisão obtida utilizando pouco número de divisões longitudinais do continuo, vantagem que justifica o uso desses elementos em estudos de pré-projeto ou otimização de estruturas. Comparam-se os resultados com soluções teóricas ou resultados de outros programas estruturais, permitindo apreciar as possibilidades e limitações da modelagem usando elementos finitos com funções spline. As diferenças observadas, que surgem principalmente em placas espessas, são explicadas pela aproximação da deformação de cisalhamento encontrada na literatura para os elementos utilizados na comparação. Mostra-se também exemplos de instabilidade analisados em três dimensões que permitem considerar diferentes condições de apoio e discutir os resultados de fórmulas conhecidas.