Higher-order zig-zag model for analysis of thick composite beams with including transverse normal stress and sublaminates approximations