Developing models that capture the spatial characteristics of wireless channels is vital for studying the performance of wireless communication systems particularly that exploit spatial processing receivers. In this paper, we propose an extension of the 2D model presented in [1] to the 3D case. We derive the joint cumulative distribution function (CDF) and the joint probability density function (pdf) for both Direction-Of-Arrival (DOA) and Direction-Of-Departure (DOD) for a single interaction multipath component. To illustrate the behavior of the proposed 3D model, extensive numerical simulations are presented.