Matrix representations for 3D strain-gradient elasticity