Recent developments have demonstrated that the behaviour of segmental grey cast iron (GCI) tunnel linings can be simulated accurately with sophisticated 3D numerical models where the geometry of the segments and joints is considered explicitly. Such a 3D model is, however, impractical for geotechnical numerical analyses where the tunnel lining is usually simulated with simpler, structural elements. In this paper, a novel tunnel joint model specifically developed for simulating the behaviour of GCI longitudinal joints with beam elements in geotechnical finite element analysis is introduced. Unlike existing models, the model considers the contribution of the bolts to the joint response and the nonlinear decay of the rotational stiffness after joint opening as well as the dependency on the compression level. After calibration of the model parameters, the model was utilised in a series of parametric analyses replicating the elastic continuum model using 2D beam elements. The response of the segmental GCI lining was essentially equivalent to that predicted by an advanced 3D model which validated the new modelling approach. The development of the joint beam model enables the nonlinear behaviour of segmental GCI tunnel linings to be accounted for in routine geotechnical numerical analyses.