This work aims to discuss the behavior of pultruded glass fiber-polymer composites related to issues not fully understood, comprising the damage and failure mechanisms, as well as the multi-crack composites response. Experimental and numerical investigations are conducted mainly on web-flange junctions (WFJ) extracted from varied profiles and bridge decks systems. Factors that favor the composites delamination, such as the influence of local defects/imperfections and the consequent modification on the junctions stresses trajectories are addressed. Theoretical equations were proposed to determine the WFJs rotational stiffness, resulting in maximum differences of 24 percent and 38 percent for I-sections and channels, respectively. To complement the scarcity of fracture parameter data related to the material, adaptations of classical fracture mechanic experiments for both Modes I and II were proposed and conducted. Nine data reduction methods were analyzed and compared with finite element models. The Modified Beam Theory method presented the best results for crack propagation in Mode I, whereas for Mode II, the Corrected Beam Theory using effective crack length and the Experimental Compliance Methods presented the closest results to numerical models. Finally, a novel approach for the multi-crack assessment is proposed. Compliance and R-curves are obtained based on the sum of the crack lengths and the strain energy release rate for each crack is evaluated separately. The results of this novel approach were compared with numerical analyses developed with the use of cohesive elements and a good agreement was found in terms of load vs. displacement response.