In order to avoid locking problems for incompressible models (such as the rubber models implemented), the weak form needs to be adjusted. This can be done through the Q1P0 family of elements, where the additional pressure variable is condensed out on the element level. However this is only valid for the discontinuous pressure field and higher order shape functions need to include the entire (u, p) function set and solve a partitioned system.
In order to avoid locking problems for incompressible models (such as the rubber models implemented), the weak form needs to be adjusted. This can be done through the Q1P0 family of elements, where the additional pressure variable is condensed out on the element level. However this is only valid for the discontinuous pressure field and higher order shape functions need to include the entire (u, p) function set and solve a partitioned system.
This is a tracking issue of the progress.