This article presents a numerical technique for computing the biaxial yield surface of polymer-matrix composites with a given microstructure. Generalized Method of Cells in combination with an Improved Bodner-Partom Viscoplastic model is used to compute the inelastic deformation. The validation of presented model is proved by a fiber Bragg gratings (FBGs) strain test system through uniaxial testing under two different strain rate conditions. On this basis, the manufacturing process thermal residual stress and strain rate effect on the biaxial yield surface of composites are considered. The results show that the effect of thermal residual stress on the biaxial yield response is closely dependent on loading conditions. Moreover, biaxial yield strength tends to increase with the increasing strain rate.