Sources Esope | Cast3M

defela
C DEFELA    SOURCE    AM        09/12/08    21:15:10     6582
      SUBROUTINE DEFELA(DEP,EPSILO,NSTRS,SIGF,IFOUR,KERRE)
C_______________________________________________________________________
C
C CALCUL DE LA DEFORMATION ELASTIQUE EN FIN DE PAS
C_______________________________________________________________________
C                                                                            C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION DEP(NSTRS,NSTRS)
      DIMENSION EPSILO(NSTRS)
      DIMENSION SIGF(NSTRS)
      DIMENSION TEMP1(3,3),TEMP2(10,10)
C
C MATRICE DEP A INVERSER
C
      NDI =10
      CALL ZERO(TEMP2,NDI,NSTRS)
C
C
      IF (IFOUR.EQ.-2) THEN
        CALL ZERO(TEMP1,3,3)
        DO 10 I=1,3
          DO 10 J=1,3
            TEMP1(I,J)=DEP(I,J)
   10   CONTINUE
C
        CALL INVMA1(TEMP1,3,3,ISING)
        IF (ISING.EQ.1) THEN
          WRITE(*,*) 'Matrice singuliere dans defela'
          KERRE = 99
          RETURN
        ENDIF
C
        TEMP2(1,1)=TEMP1(1,1)
        TEMP2(1,2)=TEMP1(1,2)
        TEMP2(2,1)=TEMP1(2,1)
        TEMP2(2,2)=TEMP1(2,2)
        TEMP2(3,3)=TEMP1(3,3)
      ELSE
        DO 30 I=1,NSTRS
          DO 30 J=1,NSTRS
            TEMP2(I,J)=DEP(I,J)
   30 CONTINUE
        CALL INVMA1(TEMP2,NDI,NSTRS,ISING)
        IF (ISING.EQ.1) THEN
          WRITE(*,*) 'Matrice singuliere dans defela'
          KERRE = 99
          RETURN
        ENDIF
      ENDIF
C
C CALCUL DU VECTEUR DEFORMATION ELASTIQUE EN FIN DE PAS
C
      DO 50 I=1,NSTRS
        EPSILO(I)=0.D0
        DO 50 J=1,NSTRS
          EPSILO(I)=EPSILO(I)+TEMP2(I,J)*SIGF(J)
   50 CONTINUE
C
      END