Télécharger cmaxwe.eso

Retour à la liste

cmaxwe
C CMAXWE    SOURCE    PV        17/12/08    21:16:02     9660           
       SUBROUTINE CMAXWE(WRK52,WRK53,WRK54,WR12,
     1      IB,IGAU,NCHAIN,DT,NWA)
C MAXWEL    SOURCE    DURAND    01/08/07    21:15:03     4179
c       SUBROUTINE MAXWEL(WRK0,WRK1,WRK5,WR12,MFR,
c     1      IB,IGAU,MELE,NCHAIN,KERRE,DT,CMATE,NWA,TEMP0)
*
*==============================================================
*      chaine de Maxwell
*==============================================================
*
* 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 PPARAM
-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
 
*
*        Formulation UNIDIMENSIONNELLE
*   ou   Formulation UNIDIRECTIONNELLE
*
      IF (MFR.EQ.27.OR.
     & ((MFR.EQ.1.OR.MFR.EQ.33).AND.( CMATE.EQ.'UNIDIREC'))) THEN
 
          CONS=1.
          IF(MFR.EQ.27) XSECT=xcarb(1)
          IF(MFR.EQ.1.OR.MFR.EQ.33) XSECT=1.D0
 
 
         CALL CMODFI(TPS1,TPS2,0,WR12,WRK52,WRK53,
     1                NWA,EFIC,TR,E1,E2)
           FLUA=EXP(-TR*DT)
           SM0(1)=SM0(1)/XSECT
           SM0(1)=SM0(1)*FLUA+CONS*EFIC*DEPST(1)
           SM0(1)=SM0(1)*XSECT
           SIGF(1)=SM0(1)
           SE1=SE1+E1
           SE2=SE2+E2
 
         CALL CMODFI(TPS1,TPS2,1,WR12,WRK52,WRK53,
     1                NWA,EFIC,TR,E1,E2)
           FLUA=EXP(-TR*DT)
           SM1(1)=SM1(1)/XSECT
           SM1(1)=SM1(1)*FLUA+CONS*EFIC*DEPST(1)
           SM1(1)=SM1(1)*XSECT
           SIGF(1)=SIGF(1)+SM1(1)
           SE1=SE1+E1
           SE2=SE2+E2
 
 
 
         CALL CMODFI(TPS1,TPS2,2,WR12,WRK52,WRK53,
     1                NWA,EFIC,TR,E1,E2)
           FLUA=EXP(-TR*DT)
           SM2(1)=SM2(1)/XSECT
           SM2(1)=SM2(1)*FLUA+CONS*EFIC*DEPST(1)
           SM2(1)=SM2(1)*XSECT
           SIGF(1)=SIGF(1)+SM2(1)
           SE1=SE1+E1
           SE2=SE2+E2
 
         CALL CMODFI(TPS1,TPS2,3,WR12,WRK52,WRK53,
     1                NWA,EFIC,TR,E1,E2)
           FLUA=EXP(-TR*DT)
           SM3(1)=SM3(1)/XSECT
           SM3(1)=SM3(1)*FLUA+CONS*EFIC*DEPST(1)
           SM3(1)=SM3(1)*XSECT
           SIGF(1)=SIGF(1)+SM3(1)
           SE1=SE1+E1
           SE2=SE2+E2
 
         CALL CMODFI(TPS1,TPS2,4,WR12,WRK52,WRK53,
     1                NWA,EFIC,TR,E1,E2)
           FLUA=EXP(-TR*DT)
           SM4(1)=SM4(1)/XSECT
           SM4(1)=SM4(1)*FLUA+CONS*EFIC*DEPST(1)
           SM4(1)=SM4(1)*XSECT
           SIGF(1)=SIGF(1)+SM4(1)
           SE1=SE1+E1
           SE2=SE2+E2
 
           IF (NCHAIN-1.GE.5) THEN
         CALL CMODFI(TPS1,TPS2,5,WR12,WRK52,WRK53,
     1                NWA,EFIC,TR,E1,E2)
             FLUA=EXP(-TR*DT)
             SM5(1)=SM5(1)/XSECT
             SM5(1)=SM5(1)*FLUA+CONS*EFIC*DEPST(1)
             SM5(1)=SM5(1)*XSECT
             SIGF(1)=SIGF(1)+SM5(1)
           SE1=SE1+E1
           SE2=SE2+E2
 
           ELSE IF (NCHAIN-1.GE.6) THEN
         CALL CMODFI(TPS1,TPS2,6,WR12,WRK52,WRK53,
     1                NWA,EFIC,TR,E1,E2)
             FLUA=EXP(-TR*DT)
             SM6(1)=SM6(1)/XSECT
             SM6(1)=SM6(1)*FLUA+CONS*EFIC*DEPST(1)
             SM6(1)=SM6(1)*XSECT
             SIGF(1)=SIGF(1)+SM6(1)
           SE1=SE1+E1
           SE2=SE2+E2
 
           ELSE IF (NCHAIN-1.GE.7) THEN
         CALL CMODFI(TPS1,TPS2,7,WR12,WRK52,WRK53,
     1                NWA,EFIC,TR,E1,E2)
             FLUA=EXP(-TR*DT)
             SM7(1)=SM7(1)/XSECT
             SM7(1)=SM7(1)*FLUA+CONS*EFIC*DEPST(1)
             SM7(1)=SM7(1)*XSECT
             SIGF(1)=SIGF(1)+SM7(1)
           SE1=SE1+E1
           SE2=SE2+E2
 
           ELSE IF (NCHAIN-1.GE.8) THEN
         CALL CMODFI(TPS1,TPS2,8,WR12,WRK52,WRK53,
     1                NWA,EFIC,TR,E1,E2)
             FLUA=EXP(-TR*DT)
             SM8(1)=SM8(1)/XSECT
             SM8(1)=SM8(1)*FLUA+CONS*EFIC*DEPST(1)
             SM8(1)=SM8(1)*XSECT
             SIGF(1)=SIGF(1)+SM8(1)
           SE1=SE1+E1
           SE2=SE2+E2
           ENDIF
 
           DEFP(1)=DEPST(1)-((SIGF(1))/(XSECT*SE2))
     1       +((SIG0(1))/(XSECT*SE1))
           EPINF(1)=EPIN0(1)-DEFP(1)
       DO I=1,NVARI
       VARF(I)=VAR0(I)
       END DO
       GOTO 999
 
*        Formulation TRIDIMENSIONNELLE MASSIVE
*
       ELSE IF (IFOUR.EQ.2.OR.IFOUR.EQ.1) THEN
         CONS=1.D0/((1.D0+ANU)*(1.D0-2.D0*ANU))
         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
     1            +CONS*EFIC*((UNANU)*DEPST(1)+ANU*(DEPST(2)+DEPST(3)))
             SM0(2)=SM0(2)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(2)+ANU*(DEPST(1)+DEPST(3)))
             SM0(3)=SM0(3)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(3)+ANU*(DEPST(1)+DEPST(2)))
             SM0(4)=SM0(4)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(4)
             SM0(5)=SM0(5)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(5)
             SM0(6)=SM0(6)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(6)
             SIGF(1)=SM0(1)
             SIGF(2)=SM0(2)
             SIGF(3)=SM0(3)
             SIGF(4)=SM0(4)
             SIGF(5)=SM0(5)
             SIGF(6)=SM0(6)
           ELSE IF (M.EQ.1) THEN
             SM1(1)=SM1(1)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(1)+ANU*(DEPST(2)+DEPST(3)))
             SM1(2)=SM1(2)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(2)+ANU*(DEPST(1)+DEPST(3)))
             SM1(3)=SM1(3)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(3)+ANU*(DEPST(1)+DEPST(2)))
             SM1(4)=SM1(4)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(4)
             SM1(5)=SM1(5)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(5)
             SM1(6)=SM1(6)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(6)
             SIGF(1)=SIGF(1)+SM1(1)
             SIGF(2)=SIGF(2)+SM1(2)
             SIGF(3)=SIGF(3)+SM1(3)
             SIGF(4)=SIGF(4)+SM1(4)
             SIGF(5)=SIGF(5)+SM1(5)
             SIGF(6)=SIGF(6)+SM1(6)
           ELSE IF (M.EQ.2) THEN
             SM2(1)=SM2(1)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(1)+ANU*(DEPST(2)+DEPST(3)))
             SM2(2)=SM2(2)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(2)+ANU*(DEPST(1)+DEPST(3)))
             SM2(3)=SM2(3)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(3)+ANU*(DEPST(1)+DEPST(2)))
             SM2(4)=SM2(4)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(4)
             SM2(5)=SM2(5)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(5)
             SM2(6)=SM2(6)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(6)
             SIGF(1)=SIGF(1)+SM2(1)
             SIGF(2)=SIGF(2)+SM2(2)
             SIGF(3)=SIGF(3)+SM2(3)
             SIGF(4)=SIGF(4)+SM2(4)
             SIGF(5)=SIGF(5)+SM2(5)
             SIGF(6)=SIGF(6)+SM2(6)
           ELSE IF (M.EQ.3) THEN
             SM3(1)=SM3(1)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(1)+ANU*(DEPST(2)+DEPST(3)))
             SM3(2)=SM3(2)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(2)+ANU*(DEPST(1)+DEPST(3)))
             SM3(3)=SM3(3)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(3)+ANU*(DEPST(1)+DEPST(2)))
             SM3(4)=SM3(4)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(4)
             SM3(5)=SM2(5)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(5)
             SM3(6)=SM2(6)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(6)
             SIGF(1)=SIGF(1)+SM3(1)
             SIGF(2)=SIGF(2)+SM3(2)
             SIGF(3)=SIGF(3)+SM3(3)
             SIGF(4)=SIGF(4)+SM3(4)
             SIGF(5)=SIGF(5)+SM3(5)
             SIGF(6)=SIGF(6)+SM3(6)
           ELSE IF (M.EQ.4) THEN
             SM4(1)=SM4(1)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(1)+ANU*(DEPST(2)+DEPST(3)))
             SM4(2)=SM4(2)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(2)+ANU*(DEPST(1)+DEPST(3)))
             SM4(3)=SM4(3)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(3)+ANU*(DEPST(1)+DEPST(2)))
             SM4(4)=SM4(4)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(4)
             SM4(5)=SM4(5)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(5)
             SM4(6)=SM4(6)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(6)
             SIGF(1)=SIGF(1)+SM4(1)
             SIGF(2)=SIGF(2)+SM4(2)
             SIGF(3)=SIGF(3)+SM4(3)
             SIGF(4)=SIGF(4)+SM4(4)
             SIGF(5)=SIGF(5)+SM4(5)
             SIGF(6)=SIGF(6)+SM4(6)
           ELSE IF (M.EQ.5) THEN
             SM5(1)=SM5(1)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(1)+ANU*(DEPST(2)+DEPST(3)))
             SM5(2)=SM5(2)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(2)+ANU*(DEPST(1)+DEPST(3)))
             SM5(3)=SM5(3)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(3)+ANU*(DEPST(1)+DEPST(2)))
             SM5(4)=SM5(4)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(4)
             SM5(5)=SM5(5)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(5)
             SM5(6)=SM5(6)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(6)
             SIGF(1)=SIGF(1)+SM5(1)
             SIGF(2)=SIGF(2)+SM5(2)
             SIGF(3)=SIGF(3)+SM5(3)
             SIGF(4)=SIGF(4)+SM5(4)
             SIGF(5)=SIGF(5)+SM5(5)
             SIGF(6)=SIGF(6)+SM5(6)
           ELSE IF (M.EQ.6) THEN
             SM6(1)=SM6(1)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(1)+ANU*(DEPST(2)+DEPST(3)))
             SM6(2)=SM6(2)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(2)+ANU*(DEPST(1)+DEPST(3)))
             SM6(3)=SM6(3)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(3)+ANU*(DEPST(1)+DEPST(2)))
             SM6(4)=SM6(4)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(4)
             SM6(5)=SM6(5)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(5)
             SM6(6)=SM6(6)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(6)
             SIGF(1)=SIGF(1)+SM6(1)
             SIGF(2)=SIGF(2)+SM6(2)
             SIGF(3)=SIGF(3)+SM6(3)
             SIGF(4)=SIGF(4)+SM6(4)
             SIGF(5)=SIGF(5)+SM6(5)
             SIGF(6)=SIGF(6)+SM6(6)
           ELSE IF (M.EQ.7) THEN
             SM7(1)=SM7(1)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(1)+ANU*(DEPST(2)+DEPST(3)))
             SM7(2)=SM7(2)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(2)+ANU*(DEPST(1)+DEPST(3)))
             SM7(3)=SM7(3)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(3)+ANU*(DEPST(1)+DEPST(2)))
             SM7(4)=SM7(4)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(4)
             SM7(5)=SM7(5)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(5)
             SM7(6)=SM7(6)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(6)
             SIGF(1)=SIGF(1)+SM7(1)
             SIGF(2)=SIGF(2)+SM7(2)
             SIGF(3)=SIGF(3)+SM7(3)
             SIGF(4)=SIGF(4)+SM7(4)
             SIGF(5)=SIGF(5)+SM7(5)
             SIGF(6)=SIGF(6)+SM7(6)
           ELSE IF (M.EQ.8) THEN
             SM8(1)=SM8(1)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(1)+ANU*(DEPST(2)+DEPST(3)))
             SM8(2)=SM8(2)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(2)+ANU*(DEPST(1)+DEPST(3)))
             SM8(3)=SM8(3)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(3)+ANU*(DEPST(1)+DEPST(2)))
             SM8(4)=SM8(4)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(4)
             SM8(5)=SM8(5)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(5)
             SM8(6)=SM8(6)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(6)
             SIGF(1)=SIGF(1)+SM8(1)
             SIGF(2)=SIGF(2)+SM8(2)
             SIGF(3)=SIGF(3)+SM8(3)
             SIGF(4)=SIGF(4)+SM8(4)
             SIGF(5)=SIGF(5)+SM8(5)
             SIGF(6)=SIGF(6)+SM8(6)
           END IF
10       CONTINUE
 
         DEFP(1)=DEPST(1)-((SIGF(1)-ANU*(SIGF(2)+SIGF(3)))/SE2
     1                   -(SIG0(1)-ANU*(SIG0(2)+SIG0(3)))/SE1)
         DEFP(2)=DEPST(2)-((SIGF(2)-ANU*(SIGF(1)+SIGF(3)))/SE2
     1                   -(SIG0(2)-ANU*(SIG0(1)+SIG0(3)))/SE1)
         DEFP(3)=DEPST(3)-((SIGF(3)-ANU*(SIGF(2)+SIGF(1)))/SE2
     1                   -(SIG0(3)-ANU*(SIG0(2)+SIG0(1)))/SE1)
         DEFP(4)=DEPST(4)-(1+ANU)*(SIGF(4)/SE2-SIG0(4)/SE1)
         DEFP(5)=DEPST(5)-(1+ANU)*(SIGF(5)/SE2-SIG0(5)/SE1)
         DEFP(6)=DEPST(6)-(1+ANU)*(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 ((IFOMOD.EQ.-1.AND.IFOUR.NE.-2).OR.
     &       (IFOMOD.EQ.0.OR.IFOMOD.EQ.1))  THEN
          CONS=1.D0/((1.D0+ANU)*(1.D0-2.D0*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
     1            +CONS*EFIC*((UNANU)*DEPST(1)+ANU*(DEPST(2)+DEPST(3)))
             SM0(2)=SM0(2)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(2)+ANU*(DEPST(1)+DEPST(3)))
             SM0(3)=SM0(3)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(3)+ANU*(DEPST(1)+DEPST(2)))
             SM0(4)=SM0(4)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(4)
             SIGF(1)=SM0(1)
             SIGF(2)=SM0(2)
             SIGF(3)=SM0(3)
             SIGF(4)=SM0(4)
           ELSE IF (M.EQ.1) THEN
             SM1(1)=SM1(1)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(1)+ANU*(DEPST(2)+DEPST(3)))
             SM1(2)=SM1(2)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(2)+ANU*(DEPST(1)+DEPST(3)))
             SM1(3)=SM1(3)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(3)+ANU*(DEPST(1)+DEPST(2)))
             SM1(4)=SM1(4)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(4)
             SIGF(1)=SIGF(1)+SM1(1)
             SIGF(2)=SIGF(2)+SM1(2)
             SIGF(3)=SIGF(3)+SM1(3)
             SIGF(4)=SIGF(4)+SM1(4)
           ELSE IF (M.EQ.2) THEN
             SM2(1)=SM2(1)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(1)+ANU*(DEPST(2)+DEPST(3)))
             SM2(2)=SM2(2)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(2)+ANU*(DEPST(1)+DEPST(3)))
             SM2(3)=SM2(3)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(3)+ANU*(DEPST(1)+DEPST(2)))
             SM2(4)=SM2(4)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(4)
             SIGF(1)=SIGF(1)+SM2(1)
             SIGF(2)=SIGF(2)+SM2(2)
             SIGF(3)=SIGF(3)+SM2(3)
             SIGF(4)=SIGF(4)+SM2(4)
           ELSE IF (M.EQ.3) THEN
             SM3(1)=SM3(1)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(1)+ANU*(DEPST(2)+DEPST(3)))
             SM3(2)=SM3(2)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(2)+ANU*(DEPST(1)+DEPST(3)))
             SM3(3)=SM3(3)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(3)+ANU*(DEPST(1)+DEPST(2)))
             SM3(4)=SM3(4)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(4)
             SIGF(1)=SIGF(1)+SM3(1)
             SIGF(2)=SIGF(2)+SM3(2)
             SIGF(3)=SIGF(3)+SM3(3)
             SIGF(4)=SIGF(4)+SM3(4)
           ELSE IF (M.EQ.4) THEN
             SM4(1)=SM4(1)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(1)+ANU*(DEPST(2)+DEPST(3)))
             SM4(2)=SM4(2)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(2)+ANU*(DEPST(1)+DEPST(3)))
             SM4(3)=SM4(3)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(3)+ANU*(DEPST(1)+DEPST(2)))
             SM4(4)=SM4(4)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(4)
             SIGF(1)=SIGF(1)+SM4(1)
             SIGF(2)=SIGF(2)+SM4(2)
             SIGF(3)=SIGF(3)+SM4(3)
             SIGF(4)=SIGF(4)+SM4(4)
           ELSE IF (M.EQ.5) THEN
             SM5(1)=SM5(1)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(1)+ANU*(DEPST(2)+DEPST(3)))
             SM5(2)=SM5(2)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(2)+ANU*(DEPST(1)+DEPST(3)))
             SM5(3)=SM5(3)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(3)+ANU*(DEPST(1)+DEPST(2)))
             SM5(4)=SM5(4)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(4)
             SIGF(1)=SIGF(1)+SM5(1)
             SIGF(2)=SIGF(2)+SM5(2)
             SIGF(3)=SIGF(3)+SM5(3)
             SIGF(4)=SIGF(4)+SM5(4)
           ELSE IF (M.EQ.6) THEN
             SM6(1)=SM6(1)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(1)+ANU*(DEPST(2)+DEPST(3)))
             SM6(2)=SM6(2)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(2)+ANU*(DEPST(1)+DEPST(3)))
             SM6(3)=SM6(3)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(3)+ANU*(DEPST(1)+DEPST(2)))
             SM6(4)=SM6(4)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(4)
             SIGF(1)=SIGF(1)+SM6(1)
             SIGF(2)=SIGF(2)+SM6(2)
             SIGF(3)=SIGF(3)+SM6(3)
             SIGF(4)=SIGF(4)+SM6(4)
           ELSE IF (M.EQ.7) THEN
             SM7(1)=SM7(1)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(1)+ANU*(DEPST(2)+DEPST(3)))
             SM7(2)=SM7(2)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(2)+ANU*(DEPST(1)+DEPST(3)))
             SM7(3)=SM7(3)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(3)+ANU*(DEPST(1)+DEPST(2)))
             SM7(4)=SM7(4)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(4)
             SIGF(1)=SIGF(1)+SM7(1)
             SIGF(2)=SIGF(2)+SM7(2)
             SIGF(3)=SIGF(3)+SM7(3)
             SIGF(4)=SIGF(4)+SM7(4)
           ELSE IF (M.EQ.8) THEN
             SM8(1)=SM8(1)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(1)+ANU*(DEPST(2)+DEPST(3)))
             SM8(2)=SM8(2)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(2)+ANU*(DEPST(1)+DEPST(3)))
             SM8(3)=SM8(3)*FLUA
     1            +CONS*EFIC*((UNANU)*DEPST(3)+ANU*(DEPST(1)+DEPST(2)))
             SM8(4)=SM8(4)*FLUA+CONS*EFIC*(1.D0-2.D0*ANU)*DEPST(4)
             SIGF(1)=SIGF(1)+SM8(1)
             SIGF(2)=SIGF(2)+SM8(2)
             SIGF(3)=SIGF(3)+SM8(3)
             SIGF(4)=SIGF(4)+SM8(4)
           END IF
20       CONTINUE
 
         DEFP(1)=DEPST(1)-((SIGF(1)-ANU*(SIGF(2)+SIGF(3)))/SE2
     1                   -(SIG0(1)-ANU*(SIG0(2)+SIG0(3)))/SE1)
         DEFP(2)=DEPST(2)-((SIGF(2)-ANU*(SIGF(1)+SIGF(3)))/SE2
     1                   -(SIG0(2)-ANU*(SIG0(1)+SIG0(3)))/SE1)
         DEFP(3)=DEPST(3)-((SIGF(3)-ANU*(SIGF(2)+SIGF(1)))/SE2
     1                   -(SIG0(3)-ANU*(SIG0(2)+SIG0(1)))/SE1)
         DEFP(4)=DEPST(4)-(1+ANU)*(SIGF(4)/SE2-SIG0(4)/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
*         CONS=1.D0/(UNANU*ANU)
         CONS=1.D0/(1.D0-(ANU*ANU))
         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
     1            +CONS*EFIC*(DEPST(1)+ANU*DEPST(2))
             SM0(2)=SM0(2)*FLUA
     1            +CONS*EFIC*(DEPST(2)+ANU*DEPST(1))
             SM0(3)=0.D0
             SM0(4)=SM0(4)*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPST(4)
             SIGF(1)=SM0(1)
             SIGF(2)=SM0(2)
             SIGF(3)=SM0(3)
             SIGF(4)=SM0(4)
           ELSE IF (M.EQ.1) THEN
             SM1(1)=SM1(1)*FLUA
     1            +CONS*EFIC*(DEPST(1)+ANU*DEPST(2))
             SM1(2)=SM1(2)*FLUA
     1            +CONS*EFIC*(DEPST(2)+ANU*DEPST(1))
             SM1(3)=0.D0
             SM1(4)=SM1(4)*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPST(4)
             SIGF(1)=SIGF(1)+SM1(1)
             SIGF(2)=SIGF(2)+SM1(2)
             SIGF(3)=SIGF(3)+SM1(3)
             SIGF(4)=SIGF(4)+SM1(4)
           ELSE IF (M.EQ.2) THEN
             SM2(1)=SM2(1)*FLUA
     1            +CONS*EFIC*(DEPST(1)+ANU*DEPST(2))
             SM2(2)=SM2(2)*FLUA
     1            +CONS*EFIC*(DEPST(2)+ANU*DEPST(1))
             SM2(3)=0.D0
             SM2(4)=SM2(4)*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPST(4)
             SIGF(1)=SIGF(1)+SM2(1)
             SIGF(2)=SIGF(2)+SM2(2)
             SIGF(3)=SIGF(3)+SM2(3)
             SIGF(4)=SIGF(4)+SM2(4)
           ELSE IF (M.EQ.3) THEN
             SM3(1)=SM3(1)*FLUA
     1            +CONS*EFIC*(DEPST(1)+ANU*DEPST(2))
             SM3(2)=SM3(2)*FLUA
     1            +CONS*EFIC*(DEPST(2)+ANU*DEPST(1))
             SM3(3)=0.D0
             SM3(4)=SM3(4)*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPST(4)
             SIGF(1)=SIGF(1)+SM3(1)
             SIGF(2)=SIGF(2)+SM3(2)
             SIGF(3)=SIGF(3)+SM3(3)
             SIGF(4)=SIGF(4)+SM3(4)
           ELSE IF (M.EQ.4) THEN
             SM4(1)=SM4(1)*FLUA
     1            +CONS*EFIC*(DEPST(1)+ANU*DEPST(2))
             SM4(2)=SM4(2)*FLUA
     1            +CONS*EFIC*(DEPST(2)+ANU*DEPST(1))
             SM4(3)=0.D0
             SM4(4)=SM4(4)*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPST(4)
             SIGF(1)=SIGF(1)+SM4(1)
             SIGF(2)=SIGF(2)+SM4(2)
             SIGF(3)=SIGF(3)+SM4(3)
             SIGF(4)=SIGF(4)+SM4(4)
           ELSE IF (M.EQ.5) THEN
             SM5(1)=SM5(1)*FLUA
     1            +CONS*EFIC*(DEPST(1)+ANU*DEPST(2))
             SM5(2)=SM5(2)*FLUA
     1            +CONS*EFIC*(DEPST(2)+ANU*DEPST(1))
             SM5(3)=0.D0
             SM5(4)=SM5(4)*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPST(4)
             SIGF(1)=SIGF(1)+SM5(1)
             SIGF(2)=SIGF(2)+SM5(2)
             SIGF(3)=SIGF(3)+SM5(3)
             SIGF(4)=SIGF(4)+SM5(4)
           ELSE IF (M.EQ.6) THEN
             SM6(1)=SM6(1)*FLUA
     1            +CONS*EFIC*(DEPST(1)+ANU*DEPST(2))
             SM6(2)=SM6(2)*FLUA
     1            +CONS*EFIC*(DEPST(2)+ANU*DEPST(1))
             SM6(3)=0.D0
             SM6(4)=SM6(4)*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPST(4)
             SIGF(1)=SIGF(1)+SM6(1)
             SIGF(2)=SIGF(2)+SM6(2)
             SIGF(3)=SIGF(3)+SM6(3)
             SIGF(4)=SIGF(4)+SM6(4)
           ELSE IF (M.EQ.7) THEN
             SM7(1)=SM7(1)*FLUA
     1            +CONS*EFIC*(DEPST(1)+ANU*DEPST(2))
             SM7(2)=SM7(2)*FLUA
     1            +CONS*EFIC*(DEPST(2)+ANU*DEPST(1))
             SM7(3)=0.D0
             SM7(4)=SM7(4)*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPST(4)
             SIGF(1)=SIGF(1)+SM7(1)
             SIGF(2)=SIGF(2)+SM7(2)
             SIGF(3)=SIGF(3)+SM7(3)
             SIGF(4)=SIGF(4)+SM7(4)
           ELSE IF (M.EQ.8) THEN
             SM8(1)=SM8(1)*FLUA
     1            +CONS*EFIC*(DEPST(1)+ANU*DEPST(2))
             SM8(2)=SM8(2)*FLUA
     1            +CONS*EFIC*(DEPST(2)+ANU*DEPST(1))
             SM8(3)=0.D0
             SM8(4)=SM8(4)*FLUA+CONS*EFIC*(UNANU)/2.D0*DEPST(4)
             SIGF(1)=SIGF(1)+SM8(1)
             SIGF(2)=SIGF(2)+SM8(2)
             SIGF(3)=SIGF(3)+SM8(3)
             SIGF(4)=SIGF(4)+SM8(4)
           END IF
30       CONTINUE
 
         DEFP(1)=DEPST(1)-((SIGF(1)-ANU*(SIGF(2)+SIGF(3)))/SE2
     1                   -(SIG0(1)-ANU*(SIG0(2)+SIG0(3)))/SE1)
         DEFP(2)=DEPST(2)-((SIGF(2)-ANU*(SIGF(1)+SIGF(3)))/SE2
     1                   -(SIG0(2)-ANU*(SIG0(1)+SIG0(3)))/SE1)
         DEFP(3)=DEPST(3)-((SIGF(3)-ANU*(SIGF(2)+SIGF(1)))/SE2
     1                   -(SIG0(3)-ANU*(SIG0(2)+SIG0(1)))/SE1)
         DEFP(4)=DEPST(4)-(1+ANU)*(SIGF(4)/SE2-SIG0(4)/SE1)
         DO I=1,NSTRS
           EPINF(I)=EPIN0(I)+DEFP(I)
         END DO
         GOTO 999
 
       ENDIF
       GOTO 999
 
 
999       RETURN
        END
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

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