Lambdoid craniosynostosis (LC) is a rare non-syndromic craniosynostosis characterised by fusion of the lambdoid sutures at the back of the head. Surgical correction including the spring assisted cranioplasty is the only option to correct the asymmetry at the skull in LC. However, the aesthetic outcome from spring assisted cranioplasty may remain suboptimal. The aim of this study is to develop a parametric finite element (FE) model of the LC skulls that could be used in the future to optimise spring surgery. The skull geometries from three different LC patients who underwent spring correction were reconstructed from the pre-operative computed tomography (CT) in Simpleware ScanIP. Initially, the skull growth between the pre-operative CT imaging and surgical intervention was simulated using MSC Marc. The osteotomies and spring implantation were performed to simulate the skull expansion due to the spring forces and skull growth between surgery and post-operative CT imaging in MSC Marc. Surface deviation between the FE models and post-operative skull models reconstructed from CT images changed between ± 5 mm over the skull geometries. Replicating spring assisted cranioplasty in LC patients allow to tune the parameters for surgical planning, which may help to improve outcomes in LC surgeries in the future.