简单的三角形网格生成程序（Simple triangular mesh generator）.doc

简单的三角形网格生成程序（Simple triangular mesh generator） Title: a triangular mesh generation procedure is simple, for beginners reference Author: XURUI [3 28, 8:32] This is a simple triangle mesh generation program using advanced front preparation method, the reference for beginners. I see Zhou Hua software download area has been empty, will initiate, hope more people can contribute to the sharing of programs. If Zhou Hua saw, hoping to take this little program to download, because here will soon be flooded. The following is the source code Program bili Open (10, file=aa.map) Open (3, file=aa) Call coor0 Close (10) End Subroutine coor0 Implicit real*4 (A-H, o-z) Common/nxy/nx, NY, numnod, nume, numb, M Common /enod/ Nod1 (80000), NOD2 (80000), NOD3 (80000) Common /coor/ x (80000), y (80000) Integer n, Nod1, NOD2, NOD3 C write (* *)input = nx; ny; Open (1, file=aa.dat) Read (1 * NE) Do i=1, NE Read (1 *) x (I), y (I) Enddo Ns = 1 C = ne (nx+ny) *2 Nume = 0 Nums = 0 Call tricb (nume, NS, NE) Numnod = nume Call trimesh (nums, numnod, nume, 5, x, y) Call trinod (nume, NT, x, y, Nod1, NOD2, NOD3) Mtij=numnod Mtik= NT Write (10, mtij, mtik *) Write (10 (2f8.2) (x) (MTII), y (MTII), mtii=1, mtij) Write (10 (3i8) (Nod1) (MTII), NOD2 (MTII), NOD3 (MTII), mtii=1, mtik) Return End C.................................................................. C the following subroutines for generating triangular mesh C according to circle boundary nodal points C trimesh generates nodal points of mesh C trinod generates nodal No. of each triangle C tricb generates boundary edges according to C continuous nodal No. of circle boundary C tricbn generates boundary edges according to C discontinuous nodal No. of circle boundary C.................................................................. C trimesh (nums, numnod, numedge, MC, x, y) C trinod (numedge, NT, x, y, Nod1, NOD2, NOD3) C tricb (numedge, NS, NE) C tricbn (numedge, numnod, node) C...........................................................