Esope source files

   1 : C DEFELA    SOURCE    AM        09/12/08    21:15:10     6582           
   2 :       SUBROUTINE DEFELA(DEP,EPSILO,NSTRS,SIGF,IFOUR,KERRE)
   3 : C_______________________________________________________________________
   4 : C
   5 : C CALCUL DE LA DEFORMATION ELASTIQUE EN FIN DE PAS
   6 : C_______________________________________________________________________
   7 : C                                                                            C   
   8 :       IMPLICIT INTEGER(I-N)
   9 :       IMPLICIT REAL*8(A-H,O-Z)
  10 :       DIMENSION DEP(NSTRS,NSTRS)
  11 :       DIMENSION EPSILO(NSTRS)
  12 :       DIMENSION SIGF(NSTRS)
  13 :       DIMENSION TEMP1(3,3),TEMP2(10,10)
  14 : C
  15 : C MATRICE DEP A INVERSER
  16 : C
  17 :       NDI =10
  17 :       CALL ZERO(TEMP2,NDI,NSTRS)
  18 : C
  19 : C
  20 :       IF (IFOUR.EQ.-2) THEN
  20 :         CALL ZERO(TEMP1,3,3)
  21 :         DO 10 I=1,3
  22 :           DO 10 J=1,3
  23 :             TEMP1(I,J)=DEP(I,J)
  24 :    10   CONTINUE     
  25 : C
  25 :         CALL INVMA1(TEMP1,3,3,ISING)
  26 :         IF (ISING.EQ.1) THEN
  27 :           WRITE(*,*) 'Matrice singuliere dans defela'
  28 :           KERRE = 99
  29 :           RETURN
  30 :         ENDIF
  31 : C
  32 :         TEMP2(1,1)=TEMP1(1,1)
  33 :         TEMP2(1,2)=TEMP1(1,2)
  34 :         TEMP2(2,1)=TEMP1(2,1)
  35 :         TEMP2(2,2)=TEMP1(2,2)
  36 :         TEMP2(3,3)=TEMP1(3,3)
  37 :       ELSE
  38 :         DO 30 I=1,NSTRS
  39 :           DO 30 J=1,NSTRS
  40 :             TEMP2(I,J)=DEP(I,J)
  41 :    30 CONTINUE     
  41 :         CALL INVMA1(TEMP2,NDI,NSTRS,ISING)
  42 :         IF (ISING.EQ.1) THEN
  43 :           WRITE(*,*) 'Matrice singuliere dans defela'
  44 :           KERRE = 99
  45 :           RETURN
  46 :         ENDIF
  47 :       ENDIF
  48 : C
  49 : C CALCUL DU VECTEUR DEFORMATION ELASTIQUE EN FIN DE PAS
  50 : C
  51 :       DO 50 I=1,NSTRS
  52 :         EPSILO(I)=0.D0
  53 :         DO 50 J=1,NSTRS
  54 :           EPSILO(I)=EPSILO(I)+TEMP2(I,J)*SIGF(J)
  55 :    50 CONTINUE     
  56 : C                                                                            
  57 :       END
  58 :  
  59 :