Phase 3 : Préparation de l'appel à l'écoulement non-linéaire

1. Pas d'intervention dans cotype.eso : en effet les paramètres matériaux et les variables internes sont de type REAL*8, qui est le type par défaut.
2. La variable LOGVIS est initialisée à .true. dans comdef.eso
   
   

         if (  inplas .eq. 17 .or.
        2   (  inplas .ge. 19 .and. inplas .le. 25) .or.
        4      inplas .eq. 61 .or.
        4      inplas .eq. 63 .or.
        4      inplas .eq. 65 .or.
        1      inplas .eq. 53 .or.
        1      inplas .eq. 29 .or.
        2      inplas .eq. 44 .or.  inplas .eq. 45  .or.
   -->  2      INPLAS .EQ. 70 .OR.  inplas .eq. 74  .or.
        9      inplas .eq. 76 .or.
        9      inplas .eq. 77 .or.  inplas .eq.  82     ) then
            logvis=.true.
         endif
   

   
   
   La variable IND prend la valeur 1, le nombre de sous-pas et la précision des itérations interne sont initialisées
   
   

    
         IF ( INPLAS .EQ. 17 .OR. INPLAS .EQ. 19 .OR.
        2     INPLAS .EQ. 20 .OR. INPLAS .EQ. 61 .OR.
        3     INPLAS .EQ. 21 .OR. INPLAS .EQ. 65 .OR.
        4     INPLAS .EQ. 22 .OR. INPLAS .EQ. 23 .OR.
        5     INPLAS .EQ. 24 .OR. INPLAS .EQ. 70 .OR.
        6     INPLAS .EQ. 25 .OR. INPLAS .EQ. 29 .OR.
        7     INPLAS .EQ. 43 .OR. INPLAS .EQ. 44 .OR.
        8     INPLAS .EQ. 45 .OR. INPLAS .EQ. 53 .OR.
        9     INPLAS .EQ. 63 .OR. INPLAS .EQ. 66 .OR.
        a     INPLAS .EQ. 74 .OR. INPLAS .EQ. 76 .OR.
        b     INPLAS .EQ. 84 .OR. INPLAS .EQ. 85 .OR.
        c     INPLAS .EQ. 86 .OR. 
        d     INPLAS .EQ. 77 .OR. INPLAS .EQ. 82 .OR.
        e     INPLAS .EQ. 106 .OR. INPLAS .EQ. 102 .OR. 
        f     INPLAS .EQ. 107 .OR. INPLAS .EQ. 108) THEN
   *
   
            IND = 1
   *
          if (inplas.eq.44.or.inplas.eq.45) then
            MSOUPA = 2000
          else
            MSOUPA = 200
          endif
            PRECIS = 5.d-7
         ENDIF
   

   
   
   

3. Il est nécessaire d'intervenir dans coml7.eso. En effet, le modèle de KOCKS est intégré par une méthode classique de Runge-Kutta pour lequel un utilitaire est disponible : le sous-programme cconst.eso. Il faut donc appeler ce sous-programme dans la séquence suivante :

   
   
   

   *---------------------------------------------------------------------
   *     modeles de viscoplasticite integres par consti
   *---------------------------------------------------------------------
         IF ( INPLAS .EQ. 17 .OR.
        2     (INPLAS .GE. 19 .AND. INPLAS .LE. 25) .OR. 
        3      INPLAS .EQ. 61 .OR. INPLAS .EQ. 63  .OR.  
        4      INPLAS .EQ. 53 .OR. INPLAS .EQ. 102 .OR.     
        5      INPLAS .EQ. 44 .OR. INPLAS .EQ. 76  .OR.
        6      INPLAS .EQ. 45 .OR. INPLAS .EQ. 77  .OR.
        7      INPLAS .EQ. 84 .OR. INPLAS .EQ. 85  .OR.
        8      INPLAS .EQ. 86 .OR. INPLAS .EQ. 70  .OR.
        9      INPLAS .EQ. 107) THEN
   * 
    ...
   *
          CALL CCONST(wrk52,wrk53,wrk54,WRK7,WRK8,WRK9,
        1        NVARI,NSSINC,INV,IFORB,TETA1,TETA2,FI1,FI2,
        4     TLIFE,NCcor,IB,IGAU,NBPGAU,KERREU1,iecou,xecou)
   
   ...
   

   
   
   

4. Dans cet exemple, l'écoulement profite de la structure de cconst.eso ; cependant pour le modèle de KOCKS il faut appeler dans cconst.eso un nouveau sous-programme incre8.eso qui réalise le calcul d'un incrément de déformation et de variables internes. le sous-programme advan1.eso a lui aussi été légèrement modifié.