Light-weight, high-strength and high-stiffness composite materials are becoming increasingly important in high-performance structural design, such as in aerospace engineering. Since these materials exhibit degradation of both strength and stiffness under a hygrothermal environment, there is an urgent need to accurately estimate hygrothermal stresses in the composite materials. In this paper, a formulation in terms of a stress function is presented for a hygrothermoelastic problem in a multiply connected region, by deriving new integral conditions necessary for the assurance of the single-valuedness of rotation and displacements. A closed-from solution as an application of the proposed formulation is given for the hygrothemoelastic problem in a hollow circular cylinder made of T 300/5208 subjected to unaxisymmetric steady-state hygrothermal loading on the inner boundary surface. Numerical calculations of the distributions of temperature, moisture and hygrothermal stresses are carried out for the cases of various hygrometric and thermal loading conditions and various ratios of the outer radius to the inner one of the cylinder.