Based on the fibre orientation distribution of a glass fibre reinforced polymer (GFRP) an anisotropic material model is developed and implemented in the commercial Finite Element code LS-Dyna. The material model is coupled with the fibre orientation extracted from flow filling simulations performed in Moldex3D. Central mathematical models used are the Eshelby-Mori-Tanaka model for calculating a composites effective stiffness, Advani and Tucker’s orientation averaging method for considering the fibre orientation tensors, a Lagrange polynomial controlled radial return algorithm for the application of nonlinear hardening and finally Hashin’s damage initiation criterion to establish the initiation and evolution of damage. The material model is calibrated by adjusting Lagrange polynomial data points, until the simulation response corresponds to that of experimental tensile tests. Two different glass fibre reinforced composites are evaluated, and it is concluded that the model can capture the behaviour of the experimental curves with high accuracy. By validation against experimental results, it is shown that the material model performs well. The model is also compared to other advanced commercial material modelling software’s.