Abstract: According to the complex plane regions in which there are discontinuous faces and multi-connected regions caused by various materials, faultages, cracks, holes or excavated regions, an algorithm of auto-arrangement of discrete nodes and its optimization are presented in this paper. In this algorithm an average step function is formed in the model of node arrangement according to the minimum step and radiation ratio, and every node as a parent node generates child nodes in terms of the step decided by the average step function. Whether the child nodes can live or not is decided by the existing nodes around them. The nodes adjacent to the edges of model are adjusted to meet the demand of computing. The algorithm has the advantages of perfect adaptability to regions with irregular shape, easy implementation and less time consumption. These auto-arranged discrete nodes are applicable to meshless method directly. The triangular and quadrangular meshes can also generated from these discrete nodes for finite element analysis.