A beam is subject to a constant axial load of magnitude \(P\), which is applied along the beam's neutral axis. Derive the element \(k_{31}\) of the local stiffness matrix for a beam element of length \(l\), including the effect of transverse displacement due to the axial load.