A simple algorithm for the stiffness matrix of triangular plate bending elements