This work aims to detail the development of a three-dimensional numerical code for simulating fracturing and fragmentation processes in geomaterials. The code is based on the combined finite-discrete element method, being able to simulate the transition from continuum to discontinuum. At first, a brief introduction to the combined finite-discrete element method is presented, and its main implementations are highlighted. In addition, the basic concepts of the cohesive zone models are shown, emphasizing the PPR (Park-Paulino-Roesler) potential based cohesive model. Then, the details of the code development are presented. The development of the code was divided into four parts: modelling of the continuum, transition from continuum to discontinuum, detection and interaction between contacts, and the solution of the equilibrium equations. The implementations made on the code are verified upon the simulation of three types of laboratory tests commonly employed in rock mechanics (Brazilian test, uniaxial compressive strength test, fracture toughness test), in which the parameters used and the results for comparison were obtained from the literature. At last, the results exhibiting the fracture patterns for each type of test and its force-displacement, or stress-strain, curves are displayed, as well as the final considerations and suggestions for future works.