This paper puts forward two different J(k)-integral-based methods, which can be used to perform mixed-mode fracture analysis of orthotropic functionally graded materials subjected to hygrothermal stresses. The first method requires the evaluation of both components of J(k)-integral, whereas the second method employs the first component.. 1 and the asymptotic crack tip displacement fields. Plane orthotropic hygrothermoelasticity is the basic theory behind the J(k)-integral formulation, which is carried out by assuming that all material properties are functions of the spatial coordinates. Developed procedures are implemented by means of the finite element method and integrated into a general purpose finite element analysis software. Temperature and specific moisture concentration fields needed in the fracture analyses are also computed through finite element analysis. Each of the developed methods is utilized in conjunction with the superposition technique to calculate the hygrothermal fracture parameters. An inclined crack located in a hygrothermally loaded orthotropic functionally graded layer is examined in parametric analyses. Comparisons of the results generated by the proposed methods do indicate that both methods lead to numerical results of high accuracy and that the developed form of the J(k)-integral is domain independent. Further results are presented so as to illustrate the influences of crack inclination angle, crack length, and crack location upon the modes I and II stress intensity factors.