In recent years, the fundamental physics of spin-lattice (e.g., magnon-phonon) interaction has attracted significant experimental and theoretical interests given its potential paradigm-shifting impacts in areas like spin-thermoelectrics, spin-caloritronics, and spintronics. Modelling studies of the transport of magnons and phonons in magnetic crystals are very rare. In this paper, we use spin-lattice dynamics (SLD) simulations to model ferromagnetic crystalline iron, where the spin and lattice systems are coupled through the atomic position-dependent exchange function, and thus the interaction between magnons and phonons is naturally considered. We then present a method combining SLD simulations with spectral energy analysis to calculate the magnon and phonon harmonic (e.g., dispersion, specific heat, and group velocity) and anharmonic (e.g., scattering rate) properties, based on which their thermal conductivity values are calculated. This work represents an example of using SLD simulations to understand the transport properties involving coupled magnon and phonon dynamics.