Many engineering devices and natural phenomena involve gels that swell under the constraint of hard materials. The constraint causes a field of stress in a gel, and often makes the swelling inhomogeneous even when the gel reaches a state of equilibrium. To analyze inhomogeneous swelling of a pH-sensitive gel, we implement a finite element method in the commercial software ABAQUS.  The program is attached here.  Contact Shenqiang Cai (shqcai@gmail.com) for a description of the program.