We address the issue of sample path simulation of Levy-driven continuous-time autoregressive moving average (CARMA) processes. Approximate discrete-time simulation schemes are constructed along with quantifiable error analysis for stable, second-order and non-negative CARMA processes, based upon the so-called series representation of infinitely divisible laws and associated Levy processes. We prove that under suitable conditions, the simulation scheme can be improved in terms of second-order structure, finite dimensional laws as well as sample path properties. The simulation procedure is often quite simple and allows one to conduct super-sampling without running the algorithm once again. The computational complexity of the proposed scheme is not affected much by the sampling scheme, such as sampling frequency and irregular spacing. Numerical results are presented throughout to illustrate the effectiveness of the proposed simulation scheme.