Implementation of a C-1 triangular element based on the P-version of the finite element method

The implementation of a computer code CONE (for C(1) continuity) based on the p-version of the finite element method is described. A hierarchic family of triangular finite elements of degree p 5 is used. This family enforces C(1)-continuity across interelement boundaries, and the code is applicable to fourth order partial differential equations in two independent variables, in particular to the biharmonic equation. Applications to several benchmark problems in plate bending are presented. Sample results are examined and compared with theoretical predictions. In particular the analysis of the bending of a rhombic plate shows a significant improvement over othr published results.