Esope source files

   1 : C ELASF1    SOURCE    AM        09/12/08    21:15:12     6582           
   2 :       SUBROUTINE ELASF1(DEP,AN,BN,CN,DEPST,SIGEF,
   3 :      $                   NSTRS,DSIGT,KERRE)
   4 : C
   5 : C ROUTINE APPELE DANS CAMECA
   6 : C_______________________________________________________________________
   7 : C
   8 : C SOUS PROGRAMME QUI CALCUL LES CONTRAINTES DE TEST FINALES
   9 : C  ELASTIQUES CORRIGE PAR LE FLUAGE TOTAL
  10 : C_______________________________________________________________________
  11 : C
  12 :       IMPLICIT INTEGER(I-N)
  13 :       IMPLICIT REAL*8(A-H,O-Z)
  14 :       DIMENSION AN(NSTRS),BN(NSTRS,NSTRS),CN(NSTRS,NSTRS)
  15 :       DIMENSION DEP(NSTRS,NSTRS),SIGEF(NSTRS),DSIGT(NSTRS)
  16 :       DIMENSION DEPST(NSTRS)
  17 : *
  18 :       DIMENSION EFLU(10,10),TEMP(10,10)
  19 :       DIMENSION DEFLUN(10),DEFELA(10)
  20 : 
  21 : C
  22 : C INITIALISATION DES VARIABLES
  23 : C
  24 :       NDI = 10
  24 :       CALL ZERO(DSIGT,NSTRS,1)
  24 :       CALL ZERO(EFLU,NDI,NSTRS)
  24 :       CALL ZERO(TEMP,NDI,NSTRS)
  24 :       CALL ZERO(DEFELA,NSTRS,1)
  24 :       CALL ZERO(DEFLUN,NSTRS,1)
  25 : C
  26 :       DO 5 I=1,NSTRS
  27 :         DO 5 J=1,NSTRS
  28 :           DO 5 K=1,NSTRS
  29 :             TEMP(I,J) = TEMP(I,J) + DEP(I,K) * CN(K,J)
  30 :   5   CONTINUE
  31 : C
  32 : C CONSTRUCTION DE LA MATRICE D ELASTICITE CORRIGE PAR LE FLUAGE
  33 : C
  34 : C  EQUATION : (2.4-1)
  35 : C
  36 :         DO 10 I=1,NSTRS
  37 :           TEMP(I,I) = 1.D0 + TEMP(I,I)
  38 :   10    CONTINUE
  39 : C
  39 :         CALL INVMA1(TEMP,NDI,NSTRS,ISING)
  40 :         IF (ISING.EQ.1) THEN
  41 :           WRITE(*,*) ' '
  42 :           WRITE(*,*) 'ERREUR DANS LE CALCUL DE LA MATRICE DE HOOKE'
  43 :           WRITE(*,*) 'CORRIGE PAR LE FLUAGE DANS ELASF1'
  44 :           WRITE(*,*) 'MATRICE NON INVERSIBLE'
  45 :           WRITE(*,*) ' '
  46 :           KERRE=99
  47 :           RETURN
  48 :         ENDIF
  49 : C
  50 : C CALCUL DE (1 + CN)-1 * E0
  51 : C
  52 :         DO 20 I=1,NSTRS
  53 :           DO 20 J=1,NSTRS
  54 :             DO 20 K=1,NSTRS
  55 :               EFLU(I,J) = EFLU(I,J) + TEMP(I,K)*DEP(K,J)
  56 :   20    CONTINUE
  57 : C
  58 : C        DO 30 I=1,NSTRS
  59 : C          DO 30 J=1,NSTRS
  60 : C            DEP(I,J) = EFLU(I,J)
  61 : C  30    CONTINUE
  62 : C
  63 : C CALCUL DES DEFORMATIONS DE FLUAGE PUIS ELASTIQUES
  64 : C
  65 :       DO 40 I=1,NSTRS
  66 :         DEFLUN(I) = DEFLUN(I) + AN(I)
  67 :         DO 40 J=1,NSTRS
  68 :           DEFLUN(I) = DEFLUN(I) + (BN(I,J)+CN(I,J)) * SIGEF(J)
  69 :   40  CONTINUE
  70 : C
  71 :       DO 50 I=1,NSTRS
  72 :         DEFELA(I) = DEPST(I) - DEFLUN(I)
  73 :   50  CONTINUE
  74 : C
  75 : C CALCUL DE L INCREMENT DE CONTRAINTES ELASTIQUES
  76 : C
  77 : C  EQUATION (2.4-2)
  78 : C
  79 :       DO 60 I=1,NSTRS
  80 :         DO 60 J=1,NSTRS
  81 :           DSIGT(I) = DSIGT(I) + EFLU(I,J) * DEFELA(J)
  82 :   60  CONTINUE
  83 : C
  84 :       END
  85 :  
  86 :