A series of numerical material tests based on homogenization theory are carried out to characterize the anisotropic elastoplastic and hardening behavior of thermoplastic resin reinforced in one direction by carbon fibers (UD-CFRTP), and a new constitutive model is proposed to describe the characteristic material behavior. First, macroscopic stress-strain curves are obtained from numerical material tests (NMTs) conducted on a unit cell model, i.e., a representative volume element of periodic microstructures, to which uniaxial and combined stress states are applied under various conditions of stress ratios and stress paths. Using the obtained NMTs results, the initial shape of the yield surface and its evolution are investigated in detail. The results suggest that even if a resin has isotopic hardening characteristics, UD-CFRTP exhibits complex anisotropic post-yielding or, equivalently, hardening behavior. Then, a constitutive model that can describe such material behavior is originally formulated, and a step-wise identification strategy is presented to determine the material parameters. Finally, by carrying out verification analyses, we confirmed that the macroscopic stress-strain curves represented by the proposed model can be fitted to those obtained by NMTs with sufficient accuracy for practical use and that the performance of the proposed model is much better than that of a conventional constitutive law.