Nearly Incompressible Materials

Nearly incompressible materials have a bulk to shear modulus ratio of several orders of magnitude. These materials can exhibit volumetric locking, that is, the computed displacements can be orders of magnitude smaller than expected. To overcome locking, Simcenter STAR-CCM+ uses a two-field approach, where both the displacement and the mean stress (or pressure) are independent variables.

With this approach, the stress-strain relationship is derived from a modified strain energy potential, which is a function of the displacement field u and the pressure p (Eqn. (4438)):

Figure 1. EQUATION_DISPLAY
Ψ˜(u,p)=Ψ(u)+Q(u,p)
(4546)

where Ψ(u) is the strain energy potential defined in Eqn. (4528) and:

Figure 2. EQUATION_DISPLAY
Q(u,p)=12kb(p¯p)2p¯=kb(J1)
(4547)

where p¯ is the internal pressure associated with the displacement field, kb is the bulk modulus, and J is the determinant of the deformation gradient (Eqn. (4431)).

The variation of the modified potential (Eqn. (4546)) produces a constraint equation relating p¯ and p. If p¯=p (that is, the constraint is exactly satisfied) then Ψ˜=Ψ.

The stress-strain relationship can then be written as:

Figure 3. EQUATION_DISPLAY
S˜=Ψ˜C
(4548)

with:

Figure 4. EQUATION_DISPLAY
S˜=S(p¯p)JC1
(4549)

where S˜ is the modified 2nd Piola-Kirchhoff stress, and C1 is the inverse of the right Cauchy Green strain tensor.