We investigate the configuration of gel sheets with centrosymmetric distribution of monomer concentration in this paper. The configuration energy of these gel sheets consists of the in-plane stretching energy and bending energy. The equilibrium shape equations are derived by variation principle. This provides a way to control the shape of gel sheets by the initial concentration and thickness. From the equilibrium shape equations, we know that the Gaussian curvature on boundary (K|C) of equilibrium shape is determined by the Poisson ratio [Formula: see text]. K|C is negative when [Formula: see text] but positive when [Formula: see text]. Specially, we derive two dome-like solutions from the equilibrium shape equations to compare with the experimental data. In these dome-like sheets, on the boundary part the Gaussian curvature is K < 0, which is different from the center part (K > 0). Furthermore, we deduce that the initial gel distribution of cylinder sheets is proportional to 1/r and find that N-isopropylacrylamide cylinder sheets cannot be formed without additional edges. Our theoretical results agree well with the experimental data [Klein et al., Science315, 1116 (2007)]. On the other hand, we predict a special type of gel sheets as minimal surface. Their residual stresses are constant and same along radial and circumference directions. For axisymmetric sheets, we give a criterion about the sign of Gaussian curvature K when thickness h is infinite small.