The popularity of pultruded glass-fibre reinforced polymer (pGFRP) profiles in construction is attributed to their lightness, strength, and durability. Despite their widespread use, a comprehensive understanding of their mechanics, particularly the phenomenon of "true" material crushing failure, remains elusive. Current methodologies for estimating the compressive resistance of pGFRP profiles rely on small-scale coupon test results [1] . However, there is compelling evidence of significant relative differences between full-section compressive strength obtained from tests and estimates derived from laminate testing. This study addresses this gap in knowledge by using computational methods to analyse the quasi-static crushing behaviour of pGFRP profiles. To that end, finite element models were developed, considering initial imperfections and geometrical and material non-linearities. In particular, the effect of end surface irregularities, measured experimentally, was investigated. A previously developed damage initiation and progression model was used [2,3], which considers the laminates as homogeneous materials, but with an additional feature – the coupling of delamination with in-plane compression – that was needed to predict more accurately the experimental results. Additionally, the effects of the variation of material properties at the web-flange junctions, including the resin-rich core and fillets, was also investigated. This study presents valuable insights into the stress and deformations induced by the above-mentioned geometric defects, with the aim of validation of the computational model providing a numerical tool that is able to accurately simulate the quasi-static crushing phenomenon.