Stress relaxation and creep of polymer gels in solvent under various deformation modes such as uniaxial, strip-biaxial, and equibiaxial were theoretically investigated. The magnitudes of relaxed stress and the creep at equilibrium under each deformation mode were derived by a thermodynamic consideration of gel system. Combining a constitutive equation of gel with the equation of motion of polymer network, the stress and strain under each deformation mode have been formulated. The theory proposed here was applied to the rectangular gels under various deformations to calculate the stress relaxation and creep behavior of polymer networks in solvent. Two methods different in treatment of swelling under the constant strain or stress were employed for the calculation: one is based on the assumption that the swelling proceeds isotropically, and the other considers the anisotropic swelling process. The results obtained by the two methods mainly differ in the diffusion mode determining the swelling behavior. The possibility of undershoot of relative strain in load-free direction in the creep is also expected.