Este trabalho tem como objetivo apresentar detalhes do desenvolvimento de um código numérico tridimensional para simulação de processos de fraturamento e fragmentação em geomateriais. O código se baseia no método combinado dos elementos finitos e discretos, sendo capaz de simular a transição do contínuo para o descontínuo. A princípio é apresentada uma breve introdução ao método combinado dos elementos finitos e discretos, e as suas principais implementações são destacadas. Além disso, são apresentados os conceitos básicos dos modelos de zona coesiva, ressaltando o modelo baseado em potencial PPR (Park-Paulino-Roesler). Em seguida, são apresentados os detalhes do desenvolvimento do código. O desenvolvimento do código foi dividido em quatro partes: modelagem do contínuo, transição entre o contínuo e o descontínuo, detecção e interação entre contatos e a resolução das equações de equilíbrio. As implementações feitas no código são verificadas mediante a simulação de três tipos de ensaios de laboratório comumente empregados na mecânica das rochas (ensaio brasileiro, ensaio de compressão simples e ensaio de tenacidade à fratura), cujos parâmetros utilizados e os resultados para comparação foram obtidos da literatura. Por fim, são apresentados os resultados exibindo os padrões de fraturamento para cada tipo de ensaio e suas curvas força-deslocamento, ou tensão-deformação, bem como as considerações finais e sugestões para trabalhos futuros.