A meso-scale finite element (FE) modelling approach is developed for meso-scale hydraulic fracture of concrete in this study. The FE models, consisting of aggregates, mortar and pores, are first generated by both random aggregate generation and packing and direct conversion of real micro-CT images. Pore pressure cohesive interface elements are then inserted into solid FE meshes by a developed Python code. The modeling approach was validated by hydro-fracturing experiments of concrete cubes. The results show that the developed models can reliably simulate complicated hydraulic fracture process, and the contents of aggregates and pores have significant effects on the hydraulic fracture resistance and the crack patterns.