A NOVEL ISOPARAMETRIC FINITE-ELEMENT DISPLACEMENT FORMULATION FOR AXISYMMETRICAL ANALYSIS OF NEARLY INCOMPRESSIBLE MATERIALS

It is well accepted that severe numerical difficulties arise when using the conventional displacement method to analyse incompressible or nearly incompressible solids. These effects are caused by the kinematic constraints imposed on the nodal velocities by the constant volume condition. In elastic-plastic analysis, these effects are due to a conflict between the plastic flow rule and the finite element discretization. Although several methods have been proposed to cope with this problem, none has been based on the appropriate choice of displacement interpolation functions to minimize the constraints. The theoretical formulation of a new six-noded isoparametric displacement finite element, which is well suited for elastic-plastic analysis of axisymmetric constrained solids by using a rational displacement interpolation function, is presented in this paper. The proposed displacement interpolation function implies that the displacement in the axial direction and the product of the displacement in the radial direction and the radius should be treated as two independent basic variables. Alternatively, the proposed displacement interpolation function can also be implemented in a conventional displacement formulation simply by using a modified shape function matrix. The suitability of the proposed formulations is first studied theoretically by assessing the number of degrees of freedom per constraint and then verified by performing numerical experiments on typical boundary value problems which involve incompressible behaviour.