Modeling the axial-torsional response of metallic strands accounting for the deformability of the internal contact surfaces: Derivation of the symmetric stiffness matrix