Two techniques are employed to eliminate volumetric locking in this problem. The first involves refining the mesh in both corners at the bottom of the rubber model to reduce mesh distortion in these areas. The second introduces a small amount of compressibility into the rubber material model. Provided the amount of compressibility is small, the results obtained with a nearly incompressible material will be very similar to those obtained with an incompressible material; the presence of compressibility alleviates the volumetric locking.

Compressibility is introduced by setting the material constant to a nonzero value. The value is chosen so that the initial Poisson's ratio, , is close to 0.5. This can be verified using the equations given in “Hyperelastic behavior,” Section 10.5.1 of the ABAQUS Analysis User's Manual, for and (the initial shear and bulk moduli, respectively) for the polynomial form of the strain energy potential. For example, using the hyperelastic material coefficients obtained earlier from the test data (see “Results,” Section 8.6.8), 176051 and 4332.63; setting 5.E–7 yields 0.46.

A model that incorporates the above features is shown in Figure 8–38 (the input file for this model is not included in the appendix; however, it can be generated easily using ABAQUS/CAE).

The displaced shape associated with this model is shown in Figure 8–39. It is clear from this figure that the mesh distortion has been reduced significantly in the critical regions of the rubber model and that the volumetric locking has been eliminated.