Bone fractures are among the most common and potentially serious injuries to the skeleton, femoral shaft fractures being especially severe. Thanks to recent advances in the area of in silico analysis, several approximations of the bone healing process have been achieved. In this context, the objective of this work was to simulate the initial phase of callus formation in long bones, without a pre-meshed domain in the 3D space. A finite element approach was computationally implemented to obtain the values of the cell concentrations along the whole domain and evaluate the areas where the biological quantities reached the thresholds necessary to trigger callus growth. A voxel model was used to obtain the 3D domain of the bone fragments and callus. A mesh growth algorithm controlled the addition of new elements to the domain at each step of the iterative procedure until complete callus formation. The implemented approach is able to reproduce the generation of the primary callus, which corresponds to the initial phase of fracture healing, independently of the fracture type and complexity, even in the case of several bone fragments. The proposed approach can be applied to the most complex bone fractures such as oblique, severely comminuted or spiral-type fractures, whose simulation remains hardly possible by means of the different existing approaches available to date.