Crater floor fractures are prominent post-cratering structural vestiges that are known from large impact craters on rocky celestial bodies. Two mechanisms have been proposed to explain the formation of crater floor fractures: emplacement of horizontal igneous sheets below crater floors and isostatic re-equilibration of crust underlying target rocks, i.e., crustal relaxation. Here, we use two-layer analogue experiments to model the deformation of lower and upper crust following crater formation, scaled to the physical conditions on Earth, to explore the structural and kinematic consequences of crustal relaxation. Specifically, the structural evolution of model upper crust was systematically analysed for various initial depths and diameters of crater floors, gleaned from previous numerical models for average continental crust. The analogue modelling results provide quantitative estimates of the duration, geometry and distribution of deformation zones in the upper crust and, for the first time, a quantitative relationship between the diameter, depth and fracture geometry of crater floors. The experiments also show that crater floor uplift is accomplished by long-wavelength subsidence of the crater periphery, which may operate on time scales of hundreds of thousands of years in nature. We conclude that patterns of natural crater floor fractures, including impact melt rock dikes known from the Sudbury and Vredefort impact structures, can be caused by long-term uplift of the crater floor, compensated by lateral crustal flow toward the crater centre.