In this paper, mechanical behavior of unidirectional composites with flexible matrix was predicted by geometrical non-linear finite element analysis. Two typical idealized unit cells of square and hexagonal fiber arrays were modeled and these were subjected to different loadings. The stress-strain behavior of composites was predicted from which the effective properties were calculated. The hyperelasticity of polyurethane matrix was considered using Mooney-Rivlin model. In result, the stress-strain behavior of flexible composites shows non-linearity, especially it is remarkable under transverse normal and shear loading conditions. In this cases, there are great difference between square and hexagonal fiber array models.