We present computationally efficient models and approaches and utilize them to investigate the dynamics of microbeams under mechanical shock. We explore using a hybrid approach utilizing a beam model combined with the shock spectrum of a spring-mass-damper model. We conclude that this approach is computationally efficient and yields accurate results in both quasi-static and dynamic loading conditions. We utilize a reduced-order model based on the nonlinear Euler-Bernoulli beam model. We demonstrate that this model is capable of capturing accurately the dynamic behavior of microbeams under shock pulses of various amplitudes (low-g and high-g), in various damping conditions, structural boundaries (clamped-clamped and clamped-free), and can capture both linear and nonlinear behavior. We investigate high-g loading cases. We report significant increase in the computational cost of simulations when using traditional nonlinear finite-element models because of the activation of higher-order modes. We demonstrate that the developed reduced-order model can be very efficient in such cases.