This research presents a simulation methodology for the Pedestrian Protection test using finite element analysis that enables the generation of robust windshield performance predictions across diverse laminate thicknesses, glass compositions, glass strengthening levels and designs. The approach outlined relies on the application of a nonlocal failure criterion based on critical energy and critical stress. Predictions generated by the model, including deceleration curves, HIC, and fracture patterns, have been validated by empirical test data demonstrating accurate results across a spectrum of windshield shapes, ply thicknesses, and glass compositions. Consequently, the outlined modeling approach may serve as a valuable tool for predicting and comparing the performance of different windshield designs and structures. Using the proposed modeling approach, we studied the performance of asymmetric laminates and unconventional windshield shapes.