简单的三角形网格生成程序.docVIP

简单的三角形网格生成程序 标题： 一个简单的三角形网格生成程序，供初学者参考 作者:xurui [3月28日,8:32] 这是一个简单的采用advanced front方法编制的三角形网格生成程序，供初学者参考。我看周华的软件下载区一直空着，就抛砖引玉一下，希望有更多的人能贡献程序共享。若周华看到后，希望能把这个小小的程序放到下载区，因为在这里很快会被淹没。 下面是源码 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.................................................................. c nums ....... starting edge No. - 1 c numnod ..... number of nodes c numedge .... number of edges c mc ......... number of iteration for modifying nodal points c nt ......... number of triangles c ns ......... first nodal point c ne ......... last nodal point c.................................................................. c x,y ............... coordinate arrays o