The linear complementarity approach has been utilized as a systematic and unified numerical process for determining the response of a rigid-plastic structure subjected to impulsive loading. However, the popular Lemke Algorithm for solving linear complementarity problems (LCP) encounters numerical instability issues whilst tracing the response of structures under extreme dynamic loading. This paper presents an efficient LCP approach with an enhanced initiation subroutine for resolving the numerical difficulties of the solver. The numerical response of the impulsively loaded structures is affected by the initial velocity profile, which if not found correctly can undermine the overall response. In the current study, the initial velocity profile is determined by a Linear Programming (LP) subroutine minimizing the energy function. An example of a uniform impulsively loaded simply supported beam is adduced to show the validity and accuracy of the proposed approach. The beam is approximated with bending hinges having infinite resistance to shear. Comparison of the numerical results to the available closed-form solution confirms the excellent performance of the approach. However, a subsequent investigation into a beam having the same support conditions and the applied loading, but with bending and shear deformation, results in numerical instability despite optimizing the initial velocity profile. Thus a more generic description of kinetics and kinematics is proposed that can further enhance the numerical efficiency of the LCP formulation. The ensuing numerical results are compared with the available close form solution to assess the accuracy and efficiency of the developed approach.