Incompressible Materials in Finite Element Method
Incompressible and nearly incompressible behavior is often encountered in solid mechanics. Natural rubber is almost incompressible and the ratio of its bulk modulus to the shear modulus is in range of thousands.
Standard displacement formulation of elastic problems tends to produce inaccurate results as Poisson’s ratio reaches 0.5 or when the material becomes completely incompressible. The displacement response of the model is poor and elements lock due to volume preserving constraints imposed by incompressible materials. This phenomenon is known as volumetric locking.
Consider an elastic solid body with almost incompressible material. The stress strain relationship for such materials is given as
the pressure strain relationship is given as
where p is the hydrostatic pressure, k is the bulk modulus and is volumetric strain.
If we increase the bulk modulus to infinity, the volumetric strain reduces to zero and pressure cannot be calculated from the above equation. Hence displacement-based finite element method cannot predict the stresses within an incompressible material.
Advanced element formulations like B-Bar, Enhanced Assumed Strain and Mixed u/p formulations are used to address the problem of incompressible materials and associated volumetric locking response.
As an example, we consider the linear analysis of the Cook’s Membrane as shown in Fig 1. The Cook’s membrane problem is a tapered panel clamped at one end and subjected to an in-plane shear loading on the free end. The material properties are taken to be E = 250 and nu = 0:4999 such that a nearly incompressible response is obtained.
Cook's Membrane
Boundary Conditions
Pure Displacement Method
In standard element formulation, tip displacement of the membrane is 2.7 mm as against a target value of ~8 mm. This is due to volumetric locking which causes the membrane to lock and produce incorrect results.
Now we study the response of the membrane and compare the results of pure displacement formulation with B-Bar, Enhanced strain and mixed u/p formulations.
B-Bar Method
Enhanced Assumed Strain
Mixed u/p
Several simulations with different mesh configurations are conducted and vertical displacement of top edge node is plotted.
In B-Bar method, the vertical displacement converges slowly. It is able to model incompressible material behavior to some extent but still incapable to predict accurate displacements for v > 0.49 and avoid volumetric locking.
Enhanced Assumed Strain method and Mixed u/p formulation converges fast within few mesh configurations and predict accurate solutions for v > 0.49
note:Script file for above analysis is available in the SCIFESOL package.