Bolted spherical joints, due to their prominent merits in installation, have been widely used in modern spatial structures. Despite significant research, there is a lack of understanding of their flexural fracture behaviour, which is important for the catastrophe prevention of the whole structure. Given the recent development to fill this knowledge gap, it is the objective of this paper to experimentally investigate the flexural bending capacity of the overall fracture section featured by a heightened neutral axis and fracture behaviour related to variable crack depth in screw threads. Accordingly, two full-scale bolted spherical joints with different bolt diameters were evaluated under three-point bending. The fracture behaviour of bolted spherical joints is first revealed with respect to typical stress distribution and fracture mode. A new theoretical flexural bending capacity expression for the fracture section with a heightened neutral axis is proposed and validated. A numerical model is then developed to estimate the stress amplification and stress intensity factors related to the crack opening (mode-I) fracture for the screw threads of these joints. The model is validated against the theoretical solutions of the thread-tooth-root model. The maximum stress of the screw thread is shown to take place at the same location as the test bolted sphere, while its magnitude can be greatly reduced with an increased thread root radius and flank angle. Finally, different design variants related to threads that have influences on the SIFs are compared, and the moderate steepness of the flank thread has been found to be efficient in reducing the joint fracture. The research findings could thus be beneficial for further improving the fracture resistance of bolted spherical joints.