We present a fully isogeometric modeling and simulation method for geometrically exact, nonlinear 3D beams with spatially varying geometric and material distributions, both along the beam axis and through its cross-section. The approach is based on the modeling of 3D beams using the Cosserat rod theory and the numerical discretization using B-Spline and NURBS parameterizations in an isogeometric collocation method. Transversally varying material constitutions are represented using non-homogeneous, functionally graded beam cross-section definitions such as laminates and continuously graded cross-sections. Furthermore, to model the axial variation of material and geometry, we introduce the parameterization of cross-section properties as splin...