In this paper, a finite element modelling for pretwisted rod elements is presented. The kinematical relations are solved to find the shape functions. The spatial element developed can have variable cross section, and can rest on elastic distributed or point support. The loading can be applied uniformly distributed, or in point, in either a local or global coordinate system. Element matrices are generated. The accuracy and the efficiency are shown in an example.