function K = TorsionalSpringElementAssemble(K,k,i,j)
% TorsionalSpringElementAssemble

 I = 2*i;
J = 2*j;

K(I,I) = K(I,I) + k(2,2);
K(I,J) = K(I,J) + k(2,4);
K(J,J) = K(J,J) + k(4,4);
K(J,I) = K(J,I) + k(4,2);
 
endfunction