Die filling is an important part of the powder handling process chain that greatly influences the characteristic structure and properties of the final part. Predictive modelling and simulation of the die-filling process can greatly contribute to the optimization of the part and the whole production procedure, e.g. by predicting the resulting powder compaction structure as a function of filling process parameters. The rheology of powders can be very difficult to model especially if heterogeneous agglomeration or time-dependent consolidation effects occur. We present a new discrete element contact force model that enables modelling complex powder flow characteristics including direct time-dependent consolidation effects and load history-dependent cohesion to describe the filling process of complex, difficult to handle powders. The model is demonstrated for simple flow and an industrial powder flow