A new computational method based on the J(k)-integral is put forward for the purpose of conducting fracture analysis of orthotropic materials subjected to hygrothermal stresses. By utilizing the constitutive relations of plane orthotropic hygrothermoelasticity, an alternative expression for the J(k)-integral is derived to replace the general limit definition. A numerical procedure is developed and integrated into a finite element analysis software to implement the proposed form of the J(k)-integral. Temperature and specific moisture concentration fields, which are required in fracture calculations, are also computed through finite element analysis. Numerical results are generated by considering an embedded crack in a polymer matrix fibrous composite laminate, that is subjected to steady-state hygrothermal loading. Comparisons of the mixed-mode stress intensity factors computed by the J(k)-integral based method to those evaluated via the displacement correlation technique demonstrate that, the proposed form of the J(k)-integral is domain independent and leads to numerical results of high accuracy. Presented parametric analyses illustrate the influences of the fiber volume fraction and the crack location on the modes I and II stress intensity factors, the energy release rate, and the T-stress.