C CMAXGE    SOURCE    KICH      09/09/11    21:16:10     6505           
       SUBROUTINE CMAXGE(WRK52,WRK53,WRK54,WR12,
     1      IB,IGAU,NCHAIN,DT,NWA)
C MAXGEN    SOURCE    AM        00/12/13    21:37:45     4045           
*       SUBROUTINE MAXGEN(WRK0,WRK1,WRK5,WR12,MFR,
*     1      IB,IGAU,MELE,NCHAIN,KERRE,DT,CMATE,NWA,TEMP0)
*
*==============================================================
*      chaine de Maxwell pour les coques et plaques
*==============================================================
*
* entrees
*        XMAT(NCOMAT) = CARACTERISTIQUES MECANIQUES DU MATERIAU
*
*        SIG0(NSTRS)  = CONTR. AU DEBUT DU PAS D'INTEGRATION 
*        DEPST(NSTRS) = INCREMENT DES DEFORM. TOTALES 
*        VAR0(NVARI)  = VARIABLES INTERNES AU DEBUT DU PAS D'INTEGRATION
*        XCAR(ICARA)  = CARACT. GEOMETRIQUES DES ELEMENTS FINIS
*
*        EPINO(NSTRS) = DEFORMATION INELASTIQUE AU DEBUT DU PAS D'INTEGRATION
*
*   WRK12
*        EMi(2,NWA(i+1))    i=0 a 8
*                        = EVOLUTION DU MODUE D'YOUNG DE LA BRANCHE i
*        SMi(NSTRS)         i=1 a 8
*                         = VARIABLES INTERNES AU DEBUT DU PAS D'INTEGRATION
*
*   IB = NUMERO DE L'ELEMENT
*   IGAU = NUMERO DU POINT DE GAUSS
*   MELE = NUMERO DE L'ELEMENT FINI (TYPE)
*   NCHAIN = NOMBRE DE BRANCHE DU MODELE DE MAXWELL
*   DT = INCREMENT DE TEMPS
*   NWA = DIMENSION DES EVOLUTIONS EMi
*   TEMP0 = TEMPS COURANT
*

* sorties
*
*        SIGF(NSTRS)= CONTR. A LA FIN DU PAS D'INTEGRATION
*        DEFP(NSTRS)= INCREMENT DES DEFORM. PLASTIQUES A LA FIN DU PAS                 
*               D'INTEGRATION
*
*        EPINF(NSTRS) = DEFORMATION INELASTIQUE A LA FIN DU PAS D'INTEGRATION
*
*==============================================================
*
      IMPLICIT INTEGER(I-N)
       IMPLICIT REAL*8(A-H,O-Z)
-INC CCOPTIO
-INC DECHE
*
*
      SEGMENT WR12
        REAL*8  EM0(2,NWA(1)),EM1(2,NWA(2)),EM2(2,NWA(3))
        REAL*8  EM3(2,NWA(4)),EM4(2,NWA(5)),EM5(2,NWA(6))
        REAL*8  EM6(2,NWA(7)),EM7(2,NWA(8)),EM8(2,NWA(9))
        REAL*8  SM0(NSTRS),SM1(NSTRS),SM2(NSTRS),SM3(NSTRS)
        REAL*8  SM4(NSTRS),SM5(NSTRS),SM6(NSTRS),SM7(NSTRS)
        REAL*8  SM8(NSTRS)
      ENDSEGMENT
*

       DIMENSION NWA(9)
       NCXMAT=XMAT(/1)
       NVARI=VAR0(/1)
       NSTRS=SIG0(/1)
       ANU=XMAT(2)
       UNANU=1.D0-ANU
       TPS1=TEMP0
       TPS2=TEMP0+DT
       SE1=0.D0
       SE2=0.D0

*
         ep=xcarb(1)

      IF(CMATE.EQ.'ISOTROPE') THEN                                               
******************************************************            
*        Formulations  COQUES MINCES - ISOTROPE      *
******************************************************
       
***        TRIDIMENSIONNEL - FOURIER                                         
         IF (IFOUR.EQ.2.OR.IFOUR.EQ.1) THEN 
                       
         CONS=1.D0/(1.D0-(ANU*ANU))
         CONS1=1.D0/(1.D0+ANU)
         ep=xcarb(1)

         DO 10 M=0,NCHAIN-1 
         CALL CMODFI(TPS1,TPS2,M,WR12,WRK52,WRK53,
     1                NWA,EFIC,TR,E1,E2)
           SE1=SE1+E1
           SE2=SE2+E2
           FLUA=EXP(-TR*DT)
           IF (M.EQ.0) THEN   
             SM0(1)=SM0(1)*FLUA+CONS*EFIC*ep*(DEPST(1)+ANU*DEPST(2))
             SM0(2)=SM0(2)*FLUA+CONS*EFIC*ep*(ANU*DEPST(1)+DEPST(2))
             SM0(3)=SM0(3)*FLUA+CONS1*EFIC*(ep/2.D0)*DEPST(3)
             SM0(4)=SM0(4)*FLUA+
     1               CONS*EFIC*((ep**3)/12.D0)*(DEPST(4)+ANU*DEPST(5))
             SM0(5)=SM0(5)*FLUA+
     1               CONS*EFIC*((ep**3)/12.D0)*(ANU*DEPST(4)+DEPST(5))
             SM0(6)=SM0(6)*FLUA+CONS1*EFIC*(((ep**3)/24.D0)*DEPST(6))
             DO I=1,6
               SIGF(I)=SM0(I)
             ENDDO
           ELSE IF (M.EQ.1) THEN        
             SM1(1)=SM1(1)*FLUA+CONS*EFIC*ep*(DEPST(1)+ANU*DEPST(2))
             SM1(2)=SM1(2)*FLUA+CONS*EFIC*ep*(ANU*DEPST(1)+DEPST(2))
             SM1(3)=SM1(3)*FLUA+CONS1*EFIC*(ep/2.D0)*DEPST(3)
             SM1(4)=SM1(4)*FLUA+
     1               CONS*EFIC*((ep**3)/12.D0)*(DEPST(4)+ANU*DEPST(5))
             SM1(5)=SM1(5)*FLUA+
     1               CONS*EFIC*((ep**3)/12.D0)*(ANU*DEPST(4)+DEPST(5))
             SM1(6)=SM1(6)*FLUA+CONS1*EFIC*(((ep**3)/24.D0)*DEPST(6))
             DO I=1,6
               SIGF(I)=SIGF(I)+SM1(I)
             ENDDO
           ELSE IF (M.EQ.2) THEN
             SM2(1)=SM2(1)*FLUA+CONS*EFIC*ep*(DEPST(1)+ANU*DEPST(2))
             SM2(2)=SM2(2)*FLUA+CONS*EFIC*ep*(ANU*DEPST(1)+DEPST(2))
             SM2(3)=SM2(3)*FLUA+CONS1*EFIC*(ep/2.D0)*DEPST(3)
             SM2(4)=SM2(4)*FLUA+
     1               CONS*EFIC*((ep**3)/12.D0)*(DEPST(4)+ANU*DEPST(5))
             SM2(5)=SM2(5)*FLUA+
     1               CONS*EFIC*((ep**3)/12.D0)*(ANU*DEPST(4)+DEPST(5))
             SM2(6)=SM2(6)*FLUA+CONS1*EFIC*(((ep**3)/24.D0)*DEPST(6))
             DO I=1,6
               SIGF(I)=SIGF(I)+SM2(I)
             ENDDO
           ELSE IF (M.EQ.3) THEN
             SM3(1)=SM3(1)*FLUA+CONS*EFIC*ep*(DEPST(1)+ANU*DEPST(2))
             SM3(2)=SM3(2)*FLUA+CONS*EFIC*ep*(ANU*DEPST(1)+DEPST(2))
             SM3(3)=SM3(3)*FLUA+CONS1*EFIC*(ep/2.D0)*DEPST(3)
             SM3(4)=SM3(4)*FLUA+
     1               CONS*EFIC*((ep**3)/12.D0)*(DEPST(4)+ANU*DEPST(5))
             SM3(5)=SM3(5)*FLUA+
     1               CONS*EFIC*((ep**3)/12.D0)*(ANU*DEPST(4)+DEPST(5))
             SM3(6)=SM3(6)*FLUA+CONS1*EFIC*(((ep**3)/24.D0)*DEPST(6))
             DO I=1,6
               SIGF(I)=SIGF(I)+SM3(I)
             ENDDO
           ELSE IF (M.EQ.4) THEN
             SM4(1)=SM4(1)*FLUA+CONS*EFIC*ep*(DEPST(1)+ANU*DEPST(2))
             SM4(2)=SM4(2)*FLUA+CONS*EFIC*ep*(ANU*DEPST(1)+DEPST(2))
             SM4(3)=SM4(3)*FLUA+CONS1*EFIC*(ep/2.D0)*DEPST(3)
             SM4(4)=SM4(4)*FLUA+
     1               CONS*EFIC*((ep**3)/1.D02)*(DEPST(4)+ANU*DEPST(5))
             SM4(5)=SM4(5)*FLUA+
     1               CONS*EFIC*((ep**3)/12.D0)*(ANU*DEPST(4)+DEPST(5))
             SM4(6)=SM4(6)*FLUA+CONS1*EFIC*(((ep**3)/24.D0)*DEPST(6))
             DO I=1,6
               SIGF(I)=SIGF(I)+SM4(I)
             ENDDO
           ELSE IF (M.EQ.5) THEN
             SM5(1)=SM5(1)*FLUA+CONS*EFIC*ep*(DEPST(1)+ANU*DEPST(2))
             SM5(2)=SM5(2)*FLUA+CONS*EFIC*ep*(ANU*DEPST(1)+DEPST(2))
             SM5(3)=SM5(3)*FLUA+CONS1*EFIC*(ep/2.D0)*DEPST(3)
             SM5(4)=SM5(4)*FLUA+
     1               CONS*EFIC*((ep**3)/12.D0)*(DEPST(4)+ANU*DEPST(5))
             SM5(5)=SM5(5)*FLUA+
     1               CONS*EFIC*((ep**3)/12.D0)*(ANU*DEPST(4)+DEPST(5))
             SM5(6)=SM5(6)*FLUA+CONS1*EFIC*(((ep**3)/24.D0)*DEPST(6))
             DO I=1,6
               SIGF(I)=SIGF(I)+SM5(I)
             ENDDO
           ELSE IF (M.EQ.6) THEN
             SM6(1)=SM6(1)*FLUA+CONS*EFIC*ep*(DEPST(1)+ANU*DEPST(2))
             SM6(2)=SM6(2)*FLUA+CONS*EFIC*ep*(ANU*DEPST(1)+DEPST(2))
             SM6(3)=SM6(3)*FLUA+CONS1*EFIC*(ep/2.D0)*DEPST(3)
             SM6(4)=SM6(4)*FLUA+
     1               CONS*EFIC*((ep**3)/12.D0)*(DEPST(4)+ANU*DEPST(5))
             SM6(5)=SM6(5)*FLUA+
     1               CONS*EFIC*((ep**3)/12.D0)*(ANU*DEPST(4)+DEPST(5))
             SM6(6)=SM6(6)*FLUA+CONS1*EFIC*(((ep**3)/24.D0)*DEPST(6))
             DO I=1,6
               SIGF(I)=SIGF(I)+SM6(I)
             ENDDO
           ELSE IF (M.EQ.7) THEN
             SM7(1)=SM7(1)*FLUA+CONS*EFIC*ep*(DEPST(1)+ANU*DEPST(2))
             SM7(2)=SM7(2)*FLUA+CONS*EFIC*ep*(ANU*DEPST(1)+DEPST(2))
             SM7(3)=SM7(3)*FLUA+CONS1*EFIC*(ep/2.D0)*DEPST(3)
             SM7(4)=SM7(4)*FLUA+
     1               CONS*EFIC*((ep**3)/12.D0)*(DEPST(4)+ANU*DEPST(5))
             SM7(5)=SM7(5)*FLUA+
     1               CONS*EFIC*((ep**3)/12.D0)*(ANU*DEPST(4)+DEPST(5))
             SM7(6)=SM7(6)*FLUA+CONS1*EFIC*(((ep**3)/24.D0)*DEPST(6))
             DO I=1,6
               SIGF(I)=SIGF(I)+SM7(I)
             ENDDO
           ELSE IF (M.EQ.8) THEN
             SM8(1)=SM8(1)*FLUA+CONS*EFIC*ep*(DEPST(1)+ANU*DEPST(2))
             SM8(2)=SM8(2)*FLUA+CONS*EFIC*ep*(ANU*DEPST(1)+DEPST(2))
             SM8(3)=SM8(3)*FLUA+CONS1*EFIC*(ep/2.D0)*DEPST(3)
             SM8(4)=SM8(4)*FLUA+
     1               CONS*EFIC*((ep**3)/12.D0)*(DEPST(4)+ANU*DEPST(5))
             SM8(5)=SM8(5)*FLUA+
     1               CONS*EFIC*((ep**3)/12.D0)*(ANU*DEPST(4)+DEPST(5))
             SM8(6)=SM8(6)*FLUA+CONS1*EFIC*(((ep**3)/24.D0)*DEPST(6))
             DO I=1,6
               SIGF(I)=SIGF(I)+SM8(I)
             ENDDO
           END IF
10       CONTINUE  

         DEFP(1)=DEPST(1)-((SIGF(1)-ANU*SIGF(2))/(ep*SE2)
     1                   -(SIG0(1)-ANU*SIG0(2))/(ep*SE1))
         DEFP(2)=DEPST(2)-((SIGF(2)-ANU*SIGF(1))/(ep*SE2)
     1                   -(SIG0(2)-ANU*SIG0(1))/(ep*SE1)) 
         DEFP(3)=DEPST(3)-(2.D0*(1+ANU)/ep)*(SIGF(3)/SE2-SIG0(3)/SE1) 
         DEFP(4)=DEPST(4)-(12.D0/(ep**3))*((SIGF(4)-ANU*SIGF(5))/SE2
     1                     -(SIG0(4)-ANU*SIG0(5))/SE1)
         DEFP(5)=DEPST(5)-(12.D0/(ep**3))*((SIGF(5)-ANU*SIGF(4))/SE2
     1                     -(SIG0(5)-ANU*SIG0(4))/SE1)
         DEFP(6)=DEPST(6)
     1         -(24.D0*(1+ANU)/ep**3)*(SIGF(6)/SE2-SIG0(6)/SE1)    
         DO I=1,NSTRS
           EPINF(I)=EPIN0(I)+DEFP(I)  
         END DO
         GOTO 999

**        Formulation DEFORMATIONS PLANES /AXISYMETRIQUE
         ELSE IF(IFOUR.EQ.0.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN

         CONS=1.D0/(1.D0-(ANU*ANU))
         DO 20 M=0,NCHAIN-1 
         CALL CMODFI(TPS1,TPS2,M,WR12,WRK52,WRK53,
     1                NWA,EFIC,TR,E1,E2)
           SE1=SE1+E1
           SE2=SE2+E2
           FLUA=EXP(-TR*DT)
           IF (M.EQ.0) THEN         
             SM0(1)=SM0(1)*FLUA+CONS*EFIC*ep*(DEPST(1)+ANU*DEPST(2))
             SM0(2)=SM0(2)*FLUA+CONS*EFIC*ep*(DEPST(2)+ANU*DEPST(1))     
             SM0(3)=SM0(3)*FLUA
     1            +CONS*EFIC*((ep**3)/12.D0)*(DEPST(3)+ANU*DEPST(4))
             SM0(4)=SM0(4)*FLUA
     1            +CONS*EFIC*((ep**3)/12.D0)*(ANU*DEPST(3)+DEPST(4))
             DO I=1,4
               SIGF(I)=SM0(I)
             ENDDO
           ELSE IF (M.EQ.1) THEN
             SM1(1)=SM1(1)*FLUA+CONS*EFIC*ep*(DEPST(1)+ANU*DEPST(2))
             SM1(2)=SM1(2)*FLUA+CONS*EFIC*ep*(DEPST(2)+ANU*DEPST(1))     
             SM1(3)=SM1(3)*FLUA
     1            +CONS*EFIC*((ep**3)/12.D0)*(DEPST(3)+ANU*DEPST(4))
             SM1(4)=SM1(4)*FLUA
     1            +CONS*EFIC*((ep**3)/12.D0)*(ANU*DEPST(3)+DEPST(4))
             DO I=1,4
               SIGF(I)=SIGF(I)+SM1(I)
             ENDDO
           ELSE IF (M.EQ.2) THEN
             SM2(1)=SM2(1)*FLUA+CONS*EFIC*ep*(DEPST(1)+ANU*DEPST(2))
             SM2(2)=SM2(2)*FLUA+CONS*EFIC*ep*(DEPST(2)+ANU*DEPST(1))     
             SM2(3)=SM2(3)*FLUA
     1            +CONS*EFIC*((ep**3)/12.D0)*(DEPST(3)+ANU*DEPST(4))
             SM2(4)=SM2(4)*FLUA
     1            +CONS*EFIC*((ep**3)/12.D0)*(ANU*DEPST(3)+DEPST(4))
             DO I=1,4
               SIGF(I)=SIGF(I)+SM2(I)
             ENDDO
           ELSE IF (M.EQ.3) THEN
             SM3(1)=SM3(1)*FLUA+CONS*EFIC*ep*(DEPST(1)+ANU*DEPST(2))
             SM3(2)=SM3(2)*FLUA+CONS*EFIC*ep*(DEPST(2)+ANU*DEPST(1))     
             SM3(3)=SM3(3)*FLUA
     1            +CONS*EFIC*((ep**3)/12.D0)*(DEPST(3)+ANU*DEPST(4))
             SM3(4)=SM3(4)*FLUA
     1            +CONS*EFIC*((ep**3)/12.D0)*(ANU*DEPST(3)+DEPST(4))
             DO I=1,4
               SIGF(I)=SIGF(I)+SM3(I)
             ENDDO
           ELSE IF (M.EQ.4) THEN
             SM4(1)=SM4(1)*FLUA+CONS*EFIC*ep*(DEPST(1)+ANU*DEPST(2))
             SM4(2)=SM4(2)*FLUA+CONS*EFIC*ep*(DEPST(2)+ANU*DEPST(1))     
             SM4(3)=SM4(3)*FLUA
     1            +CONS*EFIC*((ep**3)/12.D0)*(DEPST(3)+ANU*DEPST(4))
             SM4(4)=SM4(4)*FLUA
     1            +CONS*EFIC*((ep**3)/12.D0)*(ANU*DEPST(3)+DEPST(4))
             DO I=1,4
               SIGF(I)=SIGF(I)+SM4(I)
             ENDDO
           ELSE IF (M.EQ.5) THEN
             SM5(1)=SM5(1)*FLUA+CONS*EFIC*ep*(DEPST(1)+ANU*DEPST(2))
             SM5(2)=SM5(2)*FLUA+CONS*EFIC*ep*(DEPST(2)+ANU*DEPST(1))     
             SM5(3)=SM5(3)*FLUA
     1            +CONS*EFIC*((ep**3)/12.D0)*(DEPST(3)+ANU*DEPST(4))
             SM5(4)=SM5(4)*FLUA
     1            +CONS*EFIC*((ep**3)/12.D0)*(ANU*DEPST(3)+DEPST(4))
             DO I=1,4
               SIGF(I)=SIGF(I)+SM5(I)
             ENDDO
           ELSE IF (M.EQ.6) THEN
             SM6(1)=SM6(1)*FLUA+CONS*EFIC*ep*(DEPST(1)+ANU*DEPST(2))
             SM6(2)=SM6(2)*FLUA+CONS*EFIC*ep*(DEPST(2)+ANU*DEPST(1))     
             SM6(3)=SM6(3)*FLUA
     1            +CONS*EFIC*((ep**3)/12.D0)*(DEPST(3)+ANU*DEPST(4))
             SM6(4)=SM6(4)*FLUA
     1            +CONS*EFIC*((ep**3)/12.D0)*(ANU*DEPST(3)+DEPST(4))
             DO I=1,4
               SIGF(I)=SIGF(I)+SM6(I)
             ENDDO
           ELSE IF (M.EQ.7) THEN
             SM7(1)=SM7(1)*FLUA+CONS*EFIC*ep*(DEPST(1)+ANU*DEPST(2))
             SM7(2)=SM7(2)*FLUA+CONS*EFIC*ep*(DEPST(2)+ANU*DEPST(1))     
             SM7(3)=SM7(3)*FLUA
     1            +CONS*EFIC*((ep**3)/12.D0)*(DEPST(3)+ANU*DEPST(4))
             SM7(4)=SM7(4)*FLUA
     1            +CONS*EFIC*((ep**3)/12.D0)*(ANU*DEPST(3)+DEPST(4))
             DO I=1,4
               SIGF(I)=SIGF(I)+SM7(I)
             ENDDO
           ELSE IF (M.EQ.8) THEN
             SM8(1)=SM8(1)*FLUA+CONS*EFIC*ep*(DEPST(1)+ANU*DEPST(2))
             SM8(2)=SM8(2)*FLUA+CONS*EFIC*ep*(DEPST(2)+ANU*DEPST(1))     
             SM8(3)=SM8(3)*FLUA
     1            +CONS*EFIC*((ep**3)/12.D0)*(DEPST(3)+ANU*DEPST(4))
             SM8(4)=SM8(4)*FLUA
     1            +CONS*EFIC*((ep**3)/12.D0)*(ANU*DEPST(3)+DEPST(4))
             DO I=1,4
               SIGF(I)=SIGF(I)+SM8(I)
             ENDDO
           END IF
20       CONTINUE

         DEFP(1)=DEPST(1)-((SIGF(1)-ANU*SIGF(2))/(SE2*ep)
     1                   -(SIG0(1)-ANU*SIG0(2))/(SE1*ep))
         DEFP(2)=DEPST(2)-((SIGF(2)-ANU*SIGF(1))/(SE2*ep)
     1                   -(SIG0(2)-ANU*SIG0(1))/(SE1*ep)) 
         DEFP(3)=DEPST(3)-(12/ep**3)*((SIGF(3)-ANU*SIGF(4))/SE2
     1                   -(SIG0(3)-ANU*SIG0(4))/SE1) 
         DEFP(4)=DEPST(4)-(12/ep**3)*((SIGF(4)-ANU*SIGF(3))/SE2
     1                   -(SIG0(4)-ANU*SIG0(3))/SE1) 
         DO I=1,NSTRS
           EPINF(I)=EPIN0(I)+DEFP(I)  
         END DO
         GOTO 999     

*        Formulation CONTRAINTES PLANES 
*
        ELSE IF(IFOMOD.EQ.-1.AND.IFOUR.EQ.-2) THEN
         er=(ep**3)/12.d0
         DO 30 M=0,NCHAIN-1 
         CALL CMODFI(TPS1,TPS2,M,WR12,WRK52,WRK53,
     1                NWA,EFIC,TR,E1,E2)
           SE1=SE1+E1
           SE2=SE2+E2
           FLUA=EXP(-TR*DT)
           IF (M.EQ.0) THEN
             SM0(1)=SM0(1)*FLUA+EFIC*ep*DEPST(1)
             SM0(3)=SM0(3)*FLUA+EFIC*er*DEPST(3)
             SIGF(1)=SM0(1)
             SIGF(3)=SM0(3)
           ELSE IF (M.EQ.1) THEN
             SM1(1)=SM1(1)*FLUA+EFIC*ep*DEPST(1)
             SM1(3)=SM1(3)*FLUA+EFIC*er*DEPST(3)
             SIGF(1)=SIGF(1)+SM1(1)
             SIGF(3)=SIGF(3)+SM1(3)
           ELSE IF (M.EQ.2) THEN
             SM2(1)=SM2(1)*FLUA+EFIC*ep*DEPST(1)
             SM2(3)=SM2(3)*FLUA+EFIC*er*DEPST(3)
             SIGF(1)=SIGF(1)+SM2(1)
             SIGF(3)=SIGF(3)+SM2(3)
           ELSE IF (M.EQ.3) THEN
             SM3(1)=SM3(1)*FLUA+EFIC*ep*DEPST(1)
             SM3(3)=SM3(3)*FLUA+EFIC*er*DEPST(3)
             SIGF(1)=SIGF(1)+SM3(1)
             SIGF(3)=SIGF(3)+SM3(3)
           ELSE IF (M.EQ.4) THEN
             SM4(1)=SM4(1)*FLUA+EFIC*ep*DEPST(1)
             SM4(3)=SM4(3)*FLUA+EFIC*er*DEPST(3)
             SIGF(1)=SIGF(1)+SM4(1)
             SIGF(3)=SIGF(3)+SM4(3)
           ELSE IF (M.EQ.5) THEN
             SM5(1)=SM5(1)*FLUA+EFIC*ep*DEPST(1)
             SM5(3)=SM5(3)*FLUA+EFIC*er*DEPST(3)
             SIGF(1)=SIGF(1)+SM5(1)
             SIGF(3)=SIGF(3)+SM5(3)
           ELSE IF (M.EQ.6) THEN
             SM6(1)=SM6(1)*FLUA+EFIC*ep*DEPST(1)
             SM6(3)=SM6(3)*FLUA+EFIC*er*DEPST(3)
             SIGF(1)=SIGF(1)+SM6(1)
             SIGF(3)=SIGF(3)+SM6(3)
           ELSE IF (M.EQ.7) THEN
             SM7(1)=SM7(1)*FLUA+EFIC*ep*DEPST(1)
             SM7(3)=SM7(3)*FLUA+EFIC*er*DEPST(3)
             SIGF(1)=SIGF(1)+SM7(1)
             SIGF(3)=SIGF(3)+SM7(3)
           ELSE IF (M.EQ.8) THEN
             SM8(1)=SM8(1)*FLUA+EFIC*ep*DEPST(1)
             SM8(3)=SM8(3)*FLUA+EFIC*er*DEPST(3)
             SIGF(1)=SIGF(1)+SM8(1)
             SIGF(3)=SIGF(3)+SM8(3)
           END IF
30       CONTINUE

         DEFP(1)=DEPST(1)-(1.d0/ep)*(SIGF(1)/SE2
     1                   -SIG0(1)/SE1)
         DEFP(2)=0.D0 
         DEFP(3)=DEPST(3)-(1.d0/er)*(SIGF(3)/SE2
     1                   -SIG0(3)/SE1)
         DEFP(4)=0.D0   
         DO I=1,NSTRS
           EPINF(I)=EPIN0(I)+DEFP(I)  
         END DO
         GOTO 999          
           END IF

*        Formulation  UNIDIRECTIONNELLES
*
        ELSE IF(CMATE.EQ.'UNIDIREC') THEN                                               
         DO 40 M=0,NCHAIN-1 
         CALL CMODFI(TPS1,TPS2,M,WR12,WRK52,WRK53,
     1                NWA,EFIC,TR,E1,E2)
           SE1=SE1+E1
           SE2=SE2+E2
           FLUA=EXP(-TR*DT)
           IF (M.EQ.0) THEN
             SM0(1)=SM0(1)*FLUA+EFIC*ep*DEPST(1)
             SIGF(1)=SM0(1)
           ELSE IF (M.EQ.1) THEN
             SM1(1)=SM1(1)*FLUA+EFIC*ep*DEPST(1)
             SIGF(1)=SIGF(1)+SM1(1)
           ELSE IF (M.EQ.2) THEN
             SM2(1)=SM2(1)*FLUA+EFIC*ep*DEPST(1)
             SIGF(1)=SIGF(1)+SM2(1)
           ELSE IF (M.EQ.3) THEN
             SM3(1)=SM3(1)*FLUA+EFIC*ep*DEPST(1)
             SIGF(1)=SIGF(1)+SM3(1)
           ELSE IF (M.EQ.4) THEN
             SM4(1)=SM4(1)*FLUA+EFIC*ep*DEPST(1)
             SIGF(1)=SIGF(1)+SM4(1)
           ELSE IF (M.EQ.5) THEN
             SM5(1)=SM5(1)*FLUA+EFIC*ep*DEPST(1)
             SIGF(1)=SIGF(1)+SM5(1)
           ELSE IF (M.EQ.6) THEN
             SM6(1)=SM6(1)*FLUA+EFIC*ep*DEPST(1)
             SIGF(1)=SIGF(1)+SM6(1)
           ELSE IF (M.EQ.7) THEN
             SM7(1)=SM7(1)*FLUA+EFIC*ep*DEPST(1)
             SIGF(1)=SIGF(1)+SM7(1)
           ELSE IF (M.EQ.8) THEN
             SM8(1)=SM8(1)*FLUA+EFIC*ep*DEPST(1)
             SIGF(1)=SIGF(1)+SM8(1)
           END IF
40       CONTINUE

         DEFP(1)=DEPST(1)-(1.d0/ep)*(SIGF(1)/SE2
     1                   -SIG0(1)/SE1)
         DO I=1,NSTRS
           EPINF(I)=EPIN0(I)+DEFP(I)  
         END DO
         GOTO 999          
       ENDIF
       GOTO 999        


 
999       RETURN
        END