duda elementos finitos

[vrnc]

Hola tengo un problema a resolver una barra sometida a torsion por el metodo de los elementos finitos.Tengo la barra discretizada en dos elementos uno lineal con dos nudos y el segundo elemento es subparametrico con aproximacion cuadratica en desplazamientos.(nudo 2 3 y 4) el nudo 1 esta empotrado y el 4 sometido a un momento torsor pl ¿como hallo las funciones de forma de los elementos??y la matriz de rigidez elemental?