A computer algebra based finite element development environment

A finite element development environment based on the technical computing program Mathematica is described. The environment is used to automatically program standard element formulations and develop new elements with novel features. Source code can also be exported in a format compatible with commercial finite element program user-element facilities. The development environment is demonstrated for three mixed Petrov-Galerkin plane stress elements: a standard formulation, an advanced formulation incorporating rotational degrees of freedom and a standard formulation in which the stiffness matrix is integrated analytically, before being exported as ANSYS user elements. The results presented illustrate the accuracy of the standard mixed formulation element and the enhancement of performance when rotational degrees of freedom are added. Further, the analytically integrated element shows that computational requirements can be greatly reduced when analytical integration schemes are used in the formation.

[1]  Robert D. Cook,et al.  On improved hybrid finite elements with rotational degrees of freedom , 1989 .

[2]  D. Allman.  A compatible triangular element including vertex rotations for plane elasticity analysis , 1984 .

[3]  Thomas J. R. Hughes,et al.  A new family of stable elements for nearly incompressible elasticity based on a mixed Petrov-Galerkin finite element formulation , 1988 .

[4]  Amin Ghali,et al.  Hybrid hexahedral element for solids, plates, shells and beams by selective scaling , 1993 .

[5]  Stephen Wolfram,et al.  The Mathematica Book , 1996 .

[6]  T. Pian,et al.  Rational approach for assumed stress finite elements , 1984 .

[7]  E. Hinton,et al.  Spurious modes in two‐dimensional isoparametric elements , 1979 .

[8]  D. Mackenzie,et al.  Closed form integration of element stiffness matrices using a computer algebra system , 1995 .

[9]  Robert D. Cook,et al.  On the Allman triangle and a related quadrilateral element , 1986 .

[10]  Ian Smith,et al.  Properties of the 20-node brick , 1991 .

[11]  D. J. Kidger,et al.  Eigenvalues and eigenmodes of 8‐node brick elements , 1992 .

[12]  Theodore H. H. Pian,et al.  State-of-the-art development of hybrid/mixed finite element method , 1995 .

[13]  Nikolaos I. Ioakimidis.  Elementary applications of MATHEMATICA to the solution of elasticity problems by the finite element method , 1993 .

[14]  Shah M. Yunus,et al.  A study of different hybrid elements with and without rotational D.O.F. for plane stress/plane strain problems , 1988 .