postfr.f

关于网格自动生成delaunay算法的

***********************************************************************
*     Post Process for Fracture Mechanics
***********************************************************************
      subroutine postfr(mpoin,ng,x,mef,msrch,mtail,mnod,
     $     npoin,nef,nadrs,nodf,jac,isrch,ltail,iflnod,xi,bl,
     $     dldx,bmat,bmat2,phi,ntype,nifcls,aldm,dm,cm,kk,alcm,
     $     tu,disp,sts_node,eps_node,wn,jdata,rj)
      implicit real*8(a-h,o-z)
*
*==== dimensions for post process
*
      dimension x(2,mpoin+3)
      dimension tu(mpoin*ng)
*
*==== dimensions for J Integral
*
      dimension  xi(ng,mnod)
      dimension  nodf(3,mef),jac(0:2*mpoin+1,3)
      dimension  isrch(2,msrch), ltail(mtail)
      dimension  iflnod(mnod)
      dimension  bl(6,mnod),dldx(6,mnod,ng)
      dimension  bmat(3,mnod*ng),bmat2(3,mnod*ng)
      dimension  disp(mnod*ng)
      dimension  sts_node(4,mpoin),eps_node(4,mpoin)
      dimension  phi(2,mnod*ng)
      dimension  wn(mpoin)
      dimension  jdata(5),rj(10)
      do j= 1,jdata(1)
         ntip  = jdata(2)
         ndiv  = jdata(3)
         ndivr = jdata(4)
         ngauss= jdata(5)
         r= rj(j)
         call jinteg(mpoin,mef,msrch,mtail,mnod,npoin,nef,ng,nadrs,
     $        nodf,nifcls,aldm,dm,cm,kk,alcm,x,jac,isrch,ltail,
     $        iflnod,xi,bl,dldx,bmat,bmat2,phi,disp,tu,
     $        sts_node,eps_node,wn,ntip,ndiv,ndivr,ngauss,r,resj,tstar)
         write(1000+j,660) resj,tstar
 660     format(2e15.7)
      enddo
 100  format(a72)
 1100 return
      end