top of page

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.

EQ2.JPG
EQ3.JPG
EQ1.JPG
Cook problem
Cook problem boundary conditions
Cook panel displacement

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.

Cook panel B-Bar displacement

B-Bar Method

Cook panel EAS displacement

Enhanced Assumed Strain

Cook panel mixed u/p displacement

Mixed u/p

Several simulations with different mesh configurations are conducted and vertical displacement of top edge node is plotted.

Cook panel displacement convergence

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. 

bottom of page