Pneumatic cushioning cylinders are commonly employed for vibration and shock control. A mathematical simulation model of a double acting pneumatic cushioning cylinder designed to absorb shock loads is presented which is based on the following assumptions; ideal equation of state isentropic flow through a port conservation of mass polytropic thermodynamics single degree of freedom piston dynamics and energy equivalent linear damping. These differential equation can be solved through numerical integration using the fourth order Runge-Kutta method. An experimental study was conducted to validate the results obtained by the numerical integra-tion technique. Simulated results show good agreement with experimental data. The computer simulation model presented here has been extremely useful not only in understanding the has been extremely useful not only in understanding the basic cushioning but also in evaluating different designs.