Abstract
A method for automatic numerical generation of finite elements in the two-dimensional multiconnected region of an arbitrary shape is presented. The region of a complex shape is divided into some simply connected subregions and each subregion is mapped onto the square region in the transformed plane. The inverse functions of the mapping are determined by the numerical solutions to the boundary value problem of the simultaneous partial differential equations. After the mapped square regions are systematically divided into a set of triangular elements, each element is inversely mapped onto the physical plane by use of the functions. Finite element data for the global region are determined from those for the subregions.
