phasekinetic2.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	phasekinetic2 (nel0, time, tempel, tempo, tempmin, ae1, ae3, bs, ms, fcfer, fcper, fcbai, fgrain, frac1, frac2, frac3, frac4, frac5, x2, x3, x4, x5, gfac_f, phi_f, psi_f, cr_f, cf, gfac_p, phi_p, psi_p, cr_p, cp, gfac_b, phi_b, psi_b, cr_b, cb, phi_m, psi_m, n_m, fgfer, fgper, fgbai, qr2, qr3, qr4, kper, kbain, alpha, xeq2, xeq4, xgama, totfrac, timestep, nicool, index, theaccfact)

Function/Subroutine Documentation

phasekinetic2()

subroutine phasekinetic2	(	integer	nel0,
			time,
			tempel,
			tempo,
			tempmin,
			ae1,
			ae3,
			bs,
			ms,
			fcfer,
			fcper,
			fcbai,
			fgrain,
			frac1,
			frac2,
			frac3,
			frac4,
			frac5,
			x2,
			x3,
			x4,
			x5,
			gfac_f,
			phi_f,
			psi_f,
			cr_f,
			cf,
			gfac_p,
			phi_p,
			psi_p,
			cr_p,
			cp,
			gfac_b,
			phi_b,
			psi_b,
			cr_b,
			cb,
			phi_m,
			psi_m,
			n_m,
			fgfer,
			fgper,
			fgbai,
			qr2,
			qr3,
			qr4,
			kper,
			kbain,
			alpha,
			xeq2,
			xeq4,
			xgama,
			totfrac,
			timestep,
		integer	nicool,
		integer, dimension(nicool)	index,
		intent(in)	theaccfact )

Definition at line 29 of file phasekinetic2.F.

     
#include      "implicit_f.inc"
C-----------------------------------------------
C   I N P U T   A r g u m e n t s
C-----------------------------------------------
        INTEGER NEL0,NICOOL,INDEX(NICOOL)
        my_real ,intent(in) :: theaccfact
        my_real
     .        tempel(nel0),tempo(nel0),tempmin(nel0),timestep
        my_real
     .        time,alpha,tref,ae1,ae3,bs,ms,gsize,nu,fcfer,fcper,fcbai,
     .        fgrain,qr2,qr3,qr4,kper,kbain,xeq2,xeq4,yx,
     .        gfac_f,phi_f,psi_f,cr_f,cf,gfac_p,phi_p,psi_p,cr_p,cp,
     .        gfac_b,phi_b,psi_b,cr_b,cb,phi_m,psi_m,n_m,fgfer,fgper,fgbai
C-----------------------------------------------
        my_real
     .        frac1(nel0),frac2(nel0),frac3(nel0),frac4(nel0),frac5(nel0)
        my_real
     .        x2(nel0),x3(nel0),x4(nel0),x5(nel0),totfrac(nel0),xgama(nel0)
        my_real
     .        ftemp,ux,vx,udot,vdot,f,const,fdot,x2old(nel0),
     .        x3old(nel0),x4old(nel0),x5old(nel0),gx,gdot,dti,term1,term2 ,term3
        INTEGER I,K,J 
C-----------------------------------------------
      ! SUM OF FRAC1+FRAC2+FRAC3+FRAC4+FRAC5=1--------
       dti= one/max(timestep*theaccfact,em10)
      DO j=1,nicool ! NEL0
 
       i=index(j)
       ! CHECK TEMPERATURE FOR CORRESPONDING PHASE CHANGE
       IF(totfrac(i)<one.AND.tempel(i)<tempmin(i))THEN
       IF (tempel(i)<ae3)THEN 
        IF(tempel(i)>ae1)THEN
          IF(x2(i)<0.999)THEN
           ! FERRITE FORMATION X2
           !X2OLD(I)= FRAC2(I)/XEQ2
           x2old(i)=x2(i) 
           IF(x2(i)==zero)x2(i)=em10
           ftemp = exp(-qr2/tempel(i))*abs(tempel(i)-ae3)**3
           const = ftemp * fgfer * cf
           DO k=1,4
              ux = x2(i)     **(phi_f*(one-x2(i)))
              vx = (one-x2(i))**(psi_f*x2(i)) 
              gx = exp(cr_f *x2(i)*x2(i))
              IF (gx<zero) gx=one
              f =(x2(i)-x2old(i))*dti - const * ux*vx /gx
              udot = phi_f* ux * ((one-x2(i))/max(em10,x2(i))-log(max(em10,x2(i))))                        
              vdot = psi_f* vx * (log(one-x2(i))-x2(i)/(one-x2(i)))
              gdot = two * cr_f * x2(i) * gx
              fdot = dti-const*((udot*vx+vdot*ux)*gx-ux*vx*gdot)/gx/gx
              x2(i)= max(em20,x2(i)-f/fdot)
           ENDDO
           frac2(i)=x2(i)*xeq2
           frac1(i)=one-frac2(i)-frac3(i)-frac4(i)-frac5(i)
          ENDIF
c---------------------------------          
        ELSEIF(tempel(i)>bs)THEN 
          
         IF(x3(i)<0.999)THEN
            ! PEARLITE X3                                       
            x3old(i)= x3(i)  !FRAC3(I)/(ONE-XEQ2) !!
            IF(x3(i)==zero)x3(i)=em10
            ftemp=6.17*exp(-qr3/tempel(i))*abs(tempel(i)-ae1)**3
            const=ftemp* fgper *cp !!! adjusted
            DO k=1,4
              ux = x3(i)     **(phi_p*(one-x3(i)))
              vx = (one-x3(i))**(psi_p*x3(i)) 
              gx = exp(cr_p *x3(i)*x3(i))
              IF (gx<zero) gx=one
              f =(x3(i)-x3old(i))*dti - const * ux*vx /gx
              udot = phi_p* ux * ((one-x3(i))/max(em10,x3(i))-log(max(em10,x3(i))))                        
              vdot = psi_p* vx * (log(one-x3(i))-x3(i)/(one-x3(i)))
              gdot = two * cr_p * x3(i) * gx
              fdot = dti-const*((udot*vx+vdot*ux)*gx-ux*vx*gdot)/gx/gx
              x3(i)= max(em20,x3(i)-f/fdot)
            ENDDO       
c            FRAC3(I)=X3(I)*(ONE-X2(I))
            frac3(i)=x3(i)*(one-xeq2)
            frac1(i)=one-frac2(i)-frac3(i)-frac4(i)-frac5(i)
            ! CHECK IF FERRITE FORMATION
          ENDIF
c---------------------------------                    
        ELSEIF(tempel(i)>ms)THEN
          IF(x4(i)<0.999)THEN
           ! BAINITE  FORMATION 
            x4old(i) =  x4(i)  
            IF(x4(i) == zero) x4(i) = x3(i)  
            ! INITIALISE AU TAUX DE PERLITE
            ftemp=exp(-qr4/tempel(i)) *(tempel(i)-bs)**2
            const=ftemp* fgbai *cb
            DO k=1,4
              ux = x4(i)     **(phi_b*(one-x4(i)))
              vx = (one-x4(i))**(psi_b*x4(i)) 
              gx = exp(cr_b *x4(i)*x4(i))
              IF (gx<zero) gx=one
              f =(x4(i)-x4old(i))*dti - const * ux*vx /gx
              udot = phi_b* ux * ( (one-x4(i))/max(em10,x4(i))-log(max(em10,x4(i))) )                        
              vdot = psi_b* vx * (log(one-x4(i))-x4(i)/(one-x4(i)))
              gdot = two * cr_b * x4(i) * gx
              fdot = dti-const*((udot*vx+vdot*ux)*gx-ux*vx*gdot)/gx/gx
              x4(i)= max(em20,x4(i)-f/fdot)
           ENDDO
            frac4(i)=x4(i)*(one-x2(i))
            frac1(i)=one-frac2(i)-frac3(i)-frac4(i)-frac5(i)
          ENDIF
            
        ELSE           
          ! MARTENSITE  FORMATION X5
          IF (frac5(i)==zero)xgama(i)= frac1(i)
          x5old(i) =  x5(i)  
          term1 = alpha*(ms-tempel(i))**n_m
          term2 = max(em20,x5(i))     ** phi_m
          IF(psi_m >=zero) THEN
             term3 = (one - x5(i)) ** psi_m
          ELSE
             term3 = (one - x5(i)) ** (psi_m*two -xgama(i) )
          ENDIF
          x5(i) = x5old(i) + (term1 * term2 *term3)*(tempo(i)-tempel(i))
          frac5(i)=x5(i)*xgama(i)
          frac1(i)=one-frac2(i)-frac3(i)-frac4(i)-frac5(i)
         ENDIF
 
 
       ENDIF   
       ENDIF   
      ENDDO    
     
      RETURN