Pultruded Fiber-reinforced plastic (FRP) shapes (beams and column) are thin-walled or moderately thick-walled open or closed sections consisting of assemblies of flat panels. Due to the high strength-to-stiffness ratio of composites and thin-walled sectional geometry of FRP shapes, buckling is the most likely mode of failure before material failure for FRP shapes. In this paper, a combined analytical and experimental approach is used to characterize the lateral buckling of pultruded FRP composite cantilever I-beams. An energy method based on nonlinear plate theory is developed, and it includes shear effects and bending-twisting coupling. Three types of buckling mode shape functions (exact transcendental function, polynomial function, and ha...