Basalt columns, septarias, and mud cracks possess beautiful and intriguing crack patterns that are hard to predict because of the presence of cracks intersections and branches. The variational approach to brittle fracture provides a mathematically sound model based on minimization of the sum of bulk and fracture energies. It does not require any a priori assumption on fracture patterns and can therefore deal naturally with complex geometries. Here, we consider shrinkage cracks obtained during unidirectional drying of a colloidal suspension confined in a capillary tube. We focus on a portion of the tube where the cross-sectional shape cracks does not change as they propagate. We apply the variational approach to fracture to a tube cross-section and look for two-dimensional crack configurations minimizing the energy for a given loading level. We achieve qualitative and quantitative agreement between experiments and numerical simulations using a regularized energy (without any assumption on the cracks shape) or solutions obtained with traditional techniques (fixing the overall crack shape a priori). The results prove the efficiency of the variational approach when dealing with crack intersections and its ability to predict complex crack morphologies without any a priori assumption on their shape.