We propose a model order reduction (MOR) approach for networks containing simple and complex components. Simple components are modeled by linear ODE (and/or DAE) systems, while complex components are modeled by nonlinear PDE (and/or PDAE) systems. These systems are coupled through the network topology using the Kirchhoff laws. As application we consider MOR for electrical networks, where semiconductors form the complex components which are modeled by the transient drift-diffusion equations (DDEs). We sketch how proper orthogonal decomposition (POD) combined with discrete empirical interpolation (DEIM) and passivity-preserving balanced truncation methods for electrical circuits (PABTEC) can be used to reduce the dimension of the model. Furthermore we investigate residual-based sampling to construct reduced order models which are valid over a certain parameter range.