An All Quadrilateral Automesh Generation Technique and Explicit Integration Scheme for Finite Elements to Solve Some Elliptic Boundary Value Problems in Two Space

Abstract

This paper proposes the explicit integration scheme for a unique linear convex quadrilateral which can be obtained from an arbitrary linear triangle by joining the centroid to the midpoints of sides of the triangle. The explicit integration scheme proposed for these unique linear convex quadrilaterals is derived by using the standard transformations in two steps. We first map an arbitrary triangle into a standard right isosceles triangle by using an affine linear transformation from global (x, y) space into a local space (u, v). We discretise the standard right isosceles triangle in (u, v) space into three unique linear convex quadrilaterals. We have shown that any unique linear convex quadrilateral in (x, y) space can be mapped into one of the unique quadrilaterals in (u, v) space. We can always map these linear convex quadrilaterals into a 2-aquare in a natural space by an appropriate bilinear transformation. Using these two mappings, we have established an integral derivative product relation between the linear convex quadrilaterals in the (x, y) space interior to the arbitrary triangle and the linear convex quadrilaterals of the local space (u, v) interior to the standard right isosceles triangle. Further , we have shown that the product of global derivative integrals in (x, y) space can be expressed as a matrix triple product P ( X (2 * area of the arbitrary triangle in (x, y) space ) in which P is a geometric properties matrix and is the product of global derivative integrals in (u, v) space.We have shown that the explicit integration of the product of local derivative integrals in (u, v) space over the unique quadrilateral spanning vertices (1/3 ,1/3), (0, ½), (0, 0) , (1/2, 0) is now possible by use of symbolic processing capabilities in MATLAB which are based on Maple – V software package. In a similar manner we have found explicit integrals of ihe shape function and global derivative products as well as the product of shape functions.We use these integral values in computing stiffness matrix for some elliptic equations with constant coefficients.The proposed explicit integration scheme is a useful technique for boundary value problems governed by either a single equation or a system of second order partial differential equations. The physical applications of such problems are numerous in science, and engineering, the examples of single equations are the well known Laplace and Poisson equations with suitable boundary conditions. These problems are already studied in authors recent paper. We have demonstrated the proposed explicit integration scheme to solve the elliptic equations with constant coefficients subject to Dirichlet boundary conditions over simple polygonal domains like 1-square and 2-square.We have demonstrated the propsed scheme and technique to study five typical elliptic boundary value problems.Three FEM models with 2400,5400 and 7776 elements and 2481,5521 and 7921 nodes respectively are used.Tables of fem computed values,theoretical values at selected nodes and the finite element meshes are appended.The findings confirm with surface plots of exact solutions for the considered elliptic boundary value problems. Key words: Explicit Integration, Finite Element Method, Matlab Symbolic Mathematics, All Quadrilateral Mesh Generation Technique,Elliptic equations with constant coefficients ,Dirichlet Boundary Conditions ,Polygonal Domain,2-square,1-square