Gels are three-dimensional polymer networks with entrapped solvent (water etc.). They bear amazing features such as stimuli-responsive (temperature, PH, electric field etc.), high water content and biocompatibility and thus find a lot of applications. To understand the complex physics behind gel's swelling phenomenon, it is important to build up fundamental mechanical model and extend to complicated cases. In this dissertation, a coupled large deformation and diffusion model regarding gel's swelling behavior is presented. In this model, free-energy of the total gel is constituted by polymer stretching energy and polymer-solvent mixing energy. In-house nonlinear finite element code is implemented with fast computational capability. Complex phenomenon such as buckling and healing of cracked gel by swelling are studied. Due to the wide coverage of polymeric materials and solvents, solvent diffusion in gels not only follows Fickian diffusion law where concentration map is continuous but also follows non-Fickian diffusion law where concentration map shows high gradient. Phenomenological model with viscoelastic polymer constitutive and concentration dependent diffusivity is created. The model well captures this special diffusion phenomenon such as sharp diffusion front and distinctive swollen and unswollen region.