Local Stiffness Matrix for each Element

I see some posts about output k matrix of full model only

Can you calculate element k matrix from it shape and properties?

Hello,

I want to avoid that.

But if there is no other way I could do this.

Hello,

I would like to extract the local stiffness elements as well to compare my FE implementation to Optistructs. Is there a method to do this? 

Thanks,

-Eric