Télécharger la plaquette Cast3M


fr en

Télécharger t4m.eso

Retour à la liste

Numérotation des lignes : 

t4m

  1. C T4M       SOURCE    JK148537  23/03/20    21:15:07     11638          
  2. C t4m    SOURCE   G-M ROUX  CHAUFFAGE+REFROIDISSEMENT V6
  3.       SUBROUTINE t4m(iwrk52,iwrk53,iele,igau)
  4.       IMPLICIT INTEGER(I-N)
  5.       IMPLICIT REAL*8(A-H,O-Z)
  6. C
  7.  
  8. -INC PPARAM
  9. -INC CCOPTIO
  10. -INC CCHAMP
  11. -INC DECHE
  12.       REAL*8 TINC0,TINC1
  13.  
  14.  
  15.        TINC0 = 0.D0
  16.        TINC1 = 0.D0
  17.        wrk52 = iwrk52
  18.        wrk53 = iwrk53
  19. C>>>>>>Initialisations des champs
  20.        T0 = ture0(1)
  21.        T1 = turef(1)
  22.        Tm=(T0 + T1) / 2.
  23.        TPS0=TEMP0
  24.        TPS1=TEMPF
  25. cJK148537 28/02/23
  26. c       PHA10 = valma0(1)
  27. c       PHA20 = valma0(2)
  28. c       PHA30 = valma0(3)
  29.        PHA10 = rhas0(1)
  30.        PHA20 = rhas0(2)
  31.        PHA30 = rhas0(3)
  32.        TINC0 = valma0(1)
  33.        PHA40 = rhas0(4)
  34.        rhasF(1)=PHA10
  35.        rhasF(2)=PHA20
  36.        rhasF(3)=PHA30
  37.        XMATF(1)=TINC0
  38.          XMATF(2)=0.
  39.        rhasF(4)=PHA40
  40.        XMATF(12)=valma0(12)
  41. C>>>>>>Initialisation des constantes
  42.        PARR_Ar=valma0(9)
  43.        PARR_Br=valma0(10)
  44.        PARR_Cr=valma0(11)
  45.        AC1=valma0(2)
  46.        PARJMA_k=valma0(3)
  47.        PARJMA_m=valma0(4)
  48.        PARKM_ms=valma0(7)
  49.        PARKM_ks=valma0(8)
  50.        PARHB_k=valma0(5)
  51.        PARHB_e=valma0(6)
  52.        PARHB_n=valma0(13)
  53.        TF1=valma0(14)
  54.        TF2=valma0(15)
  55. C
  56. C>>>>>germination de la phase austénitique
  57.        IF (Tm.GT.AC1.AND.TINC0.LT.1.) THEN
  58.          t_inc1= Tm - AC1
  59.          t_inc2 = (PARR_Cr / t_inc1)
  60.          t_inc3 = EXP(t_inc2)
  61.          t_inc4 = PARR_Br - Tm
  62.          IF (t_inc4.LE.0.) THEN
  63.            XMATF(1)=1.
  64.          ELSE
  65.            t_inc=t_inc3 * t_inc4 * PARR_Ar
  66.            TINC1 = TINC0 + ((TPS1 - TPS0) / t_inc)
  67.            XMATF(1)=TINC1
  68.          ENDIF
  69.          IF (TINC1.GT.1) THEN
  70.            XMATF(1) = 1.
  71.          ENDIF
  72.        ENDIF
  73. C      ->si la germination est terminé, tinc=1
  74.          IF (Tm.GE.AC1.AND.TINC0.EQ.1.OR.PHA20.GT.0) THEN
  75.            XMATF(1) = 1.
  76.          ENDIF
  77. C      ->si T<Ms et yg=0
  78.          IF (Tm.LT.PARKM_ms.AND.PHA20.EQ.0) THEN
  79.            XMATF(1) = 0.
  80.          ENDIF
  81. C      ->si T>Assat
  82.          IF (Tm.GT.PARR_Br) THEN
  83.            XMATF(1) = 1.
  84.          ENDIF
  85. C
  86. C>>>>>>croissance de la phase austénitique (méthode de Heun)
  87.        IF (T1.GT.AC1.AND.XMATF(1).EQ.1.AND.PHA20.LT.1) THEN
  88.          IF (T1.LE.1000) THEN
  89. C      ->on passes les caractéristiques thermiques en Kelvin
  90.          AC1=AC1 + 273.
  91.          T1b=T1 + 273.
  92.          T0b=T0 + 273.
  93. C      ->calcul de Teq(t=ti)
  94.          Teq_11 = AC1 * PARJMA_k
  95.          Teq_12 = 1. - PHA20
  96.          Teq_13 = LOG(Teq_12)
  97.          Teq_14 = Teq_11 + ((-1. * Teq_13) ** (1. / PARJMA_m))
  98.          Teq_1 = (Teq_14 / PARJMA_k)
  99. C      ->calcul de dyA/dt(t=ti)
  100.          dPHA2_1_1 = (-1. * PARHB_e) / (8.31 * T0b)
  101.          dPHA2_1_2 = PARHB_k * EXP(dPHA2_1_1)
  102.           dPHA2_1_3=ABS(T0b - Teq_1)
  103.          dPHA2_1 = dPHA2_1_2 * (dPHA2_1_3 ** PARHB_n) * (1. - PHA20)
  104. C      ->calcul de yA(t=ti)
  105.          PHA2_1=PHA20 + dPHA2_1 * (TPS1 - TPS0)
  106.          IF (PHA2_1.GT.9.999d-01) THEN
  107.            PHA21=1.
  108.          ELSE
  109. C          ->calcul de Teq(t=ti+1)
  110.            Teq_21 = AC1 * PARJMA_k
  111.            Teq_22 = 1. - PHA2_1
  112.            Teq_23 = LOG(Teq_22)
  113.            Teq_24 = Teq_21 + ((-1. * Teq_23) ** (1. / PARJMA_m))
  114.            Teq_2 = (Teq_24 / PARJMA_K)
  115. C          ->calcul de dyA/dt(t=ti+1)
  116.            dPHA2_2_1 = (-1. * PARHB_e) / (8.31 * T1b)
  117.            dPHA2_2_2 = PARHB_k * EXP(dPHA2_2_1)
  118.              dPHA2_2_3=ABS(T1b - Teq_2)
  119.            dPHA2_2 = dPHA2_2_2 * (dPHA2_2_3 ** PARHB_n) * (1. - PHA2_1)
  120. C          ->calcul de yA(t)(ti+1)
  121.            PHA21=PHA20 + ( (dPHA2_1 + dPHA2_2) / 2.)*(TPS1 - TPS0)
  122.            IF (PHA21.GT.9.9d-01) THEN
  123.              PHA21=1.
  124.            ENDIF
  125.          ENDIF
  126.          IF (PHA30.GT.0.AND.PHA10.GT.0.) THEN
  127. C          ->calcul des nouvelles phases si martensite+martensite revenue
  128.            PHA11=PHA10-(PHA21-PHA20)/2.
  129.            PHA31=PHA30-(PHA21-PHA20)/2.
  130.            IF (PHA11.GT.0.) THEN
  131.              PHA11=0.
  132.              PHA31=1 - PHA21
  133.            ENDIF
  134.            IF (PHA31.GT.0.) THEN
  135.              PHA31=0.
  136.              PHA11=1 - PHA21
  137.            ENDIF
  138.            rhasF(2)=PHA21
  139.            rhasF(1)=PHA11
  140.            rhasF(3)=PHA31
  141.            XMATF(12)=PHA21
  142.          ELSE
  143.            IF (PHA10.GT.0.) THEN
  144. C            ->calcul des nouvelles phases si martensite revenue
  145.              rhasF(2)=PHA21
  146.              rhasF(1)=1. - PHA21
  147.              rhasF(3)=0.
  148.              XMATF(12)=PHA21
  149.            ELSE
  150. C          ->calcul des nouvelles phases si martensite
  151.              rhasF(2)=PHA21
  152.              rhasF(3)=1. - PHA21
  153.              rhasF(1)=0.
  154.              XMATF(12)=PHA21
  155.            ENDIF
  156.          ENDIF
  157.        ENDIF
  158.        ENDIF
  159.        IF (T1.GT.1000.AND.T1.LT.TF1) THEN
  160.          rhasF(2)=1.
  161.          rhasF(3)=0.
  162.          rhasF(1)=0.
  163.          XMATF(12)=1.
  164.          ENDIF
  165. C>>>>>transformation liquide-solide(loi linéaire)
  166.         IF (T1.GE.TF1.AND.T1.LE.TF2) THEN
  167.           rhasF(4)=(T1 - TF1) * (1.d0 / (TF2 - TF1))
  168.           rhasF(2)=1.d0 - rhasF(4)
  169.             rhasF(1)=0.
  170.             rhasF(3)=0.
  171.             XMATF(12)=1.
  172.         ENDIF
  173.         IF (T1.GE.TF2) THEN
  174.           rhasF(4)=1.
  175.           rhasF(2)=0.
  176.             rhasF(1)=0.
  177.             XMATF(12)=1.
  178.             XMATF(1)=1.
  179.             rhasF(3)=0.
  180.         ENDIF
  181.         IF (T1.LT.TF1.AND.PHA40.GT.0.) THEN
  182.           rhasF(4)=0.
  183.           rhasF(2)=1.
  184.             rhasF(1)=0.
  185.             XMATF(12)=1.
  186.             rhasF(3)=0.
  187.         ENDIF
  188. C
  189. C>>>>>transformation martensitique
  190.        IF (T1.LT.T0.AND.PHA20.GT.0.AND.T1.LE.PARKM_ms) THEN
  191.          PHA31_1 = -1. * PARKM_ks *(PARKM_ms - T1)
  192.          PHA31_2 = EXP(PHA31_1)
  193.          PHA31_3=XMATF(12) * (1-PHA31_2)
  194.          PHA11=PHA10
  195.          PHA21=XMATF(12) - PHA31_3
  196.          PHA31=1 - (PHA11 + PHA21)
  197.          IF (PHA21.LT.1.d-03) THEN
  198.            PHA21=0.
  199.            PHA11=PHA10
  200.            PHA31=1 - PHA11
  201.              rhasF(4)=0.
  202.          ENDIF
  203. C        ->sorties
  204.          rhasF(1)=PHA11
  205.          rhasF(2)=PHA21
  206.          rhasF(3)=PHA31
  207.        ENDIF
  208. C
  209. C>>>>>transformation martensite revenue
  210.          rhasF(5)=0. 
  211. C
  212. c       if(iele.eq.1.and.igau.eq.3) write(6,*) (rhasF(jj),jj= 1,5)
  213.        do 80 jj = 2,15
  214.          if (jj.eq.12) goto 80
  215.         xmatf(jj) = valma0(jj)
  216.   80   continue
  217. C
  218.        RETURN
  219.        END
  220.  
  221.  
  222.  
  223.  
  224.  
  225.  
  226.  
  227.  
  228.  
  229.  
  230.  
  231.  
  232.  
  233.  
  234.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales