Energy absorbed during damage depends on the toughness of solids. Material with high toughness can dissipate more energy for the same damage compared to a lower-toughness material. Hence the use of high-toughness material can add to the energy dissipation mechanism of structures. Bio-inspired composites that combine a soft phase with the majority of a hard phase referred to as brick and mortar (BM) structure, show an enhanced fracture toughness compared to the constituent materials. Recent experimental studies show that topologically interlocked composites possess even superior properties than conventional bio-inspired composites with BM-type structures. In this study, we examined how different interlocking angles influence the mechanical behavior of topologically interlocked composites. Specifically, we developed open-source codes by implementing a phase-field model (PFM) using a Julia-based new finite element (FE) library, Gridap. These codes serve as a numerical aid to investigate the influence of interlocking angles on fracture toughness. Gridap enables a compact and user-friendly implementation, requiring minimal memory usage and offering users significant flexibility in writing weak forms of partial differential equations. In addition, Gridap provides a very convenient way to impose displacement continuity at the material interface.