The Barcelona Basic Model (BBM) is an elasto-plastic model used to describe the coupled thermo-hydro-mechanical behaviors of unsaturated soil. BBM is frequently adopted to model the unique swelling behavior of bentonite, which is generally considered as the buffer material between the host rock and the canister containing high-level radioactive waste in deep geological repositories, under the changing thermal, hydraulic, mechanical and chemical conditions during the lifetime of repository. Therefore, a variety of the continuum-based numerical methods tried to add the BBM for modelling the multi barrier systems of geological repository and succeeded to describe the elasto-plastic deformation of bentonite. However, to demonstrate the entire barrier systems the host rock should be modelled simultaneously with the buffer materials, and the continuum-based methods may be limited in their ability to reflect the fracture networks in the host rock which could be the major flow channels of groundwater. This research applies BBM in 3DEC, a three-dimensional block-based discrete element method, and validates the model by comparing the change of specific volume and mean effective stress during three numerical test cases. Discontinuum-based numerical methods with BBM can be extended to describe the coupled thermo-hydro-mechanical processes of multi-barrier systems in geological repositories, with a focus on the interaction between the host rock and bentonite.