During the past few decades, several studies have been conducted to understand the behaviour of powders in vibrated beds. This paper introduces a technique of incorporating the agglomeration and deagglomeration phenomena in the simulation of vibrated fine powders. Two-dimensional direct simulations are performed using 300 spheres 2.99 mm in diameter in a trapezoidal container vibrated vertically at an amplitude of 2.5 mm and 20 Hz frequency as preliminary conditions. Under non-cohesive conditions, the results are in agreement with those found in the Literature. As a preliminary effort to predict the behaviour of cohesive fine powders under vibrated conditions, agglomeration and deagglomeration processes are modelled as the formation and destruction, respectively, of interparticle bonds during particle collisions. Two parameters used to model agglomeration and deagglomeration are the ease of cohesion and cohesivity of the powder. Dependencies of these parameters on certain physical properties of cohesive powders have been suggested. Simulation results reveal two aggregate populations, one with uniform size aggregates and another population with multi-sized aggregates. The former aggregates were more prevalent in weakly cohesive powders while the latter in highly cohesive powders. Interesting macroscopic bed behaviour such as alternating cycles of agglomeration and deagglomeration were also observed. Further work is needed in which the aerodynamic forces are taken into account and cohesion mechanisms at the particle surface are modelled.