Esope source files

   1 : C PARINV    SOURCE    CHAT      05/01/13    02:09:11     5004           
   2 : C PARINV     SOURCE    INSL       24/10/96    
   3 :       SUBROUTINE PARINV(S,SX,SY,SZ,SXY,SXZ,SYZ,NSTRS)
   4 : C             ****************************************************
   5 : C             *  CALCUL DES IVARIANT : I1, J2, J3 et COS(3TETA)  *
   6 : C             *           (  MODEL D'OTTOSEN   )                 *
   7 : C             ****************************************************
   8 :       IMPLICIT INTEGER(I-N)
   9 :       IMPLICIT REAL*8(A-H,O-Z)
  10 :       DIMENSION S1(6),S(NSTRS)
  11 : C
  12 : C      INITIALISATION
  13 : C
  13 :         CALL ZERO(S1,6,1)
  14 :         DO 11 I=1,NSTRS
  15 :           S1(I)=S(I)
  16 :    11   CONTINUE
  17 : C
  18 :         SX=0.D0
  19 :         SY=0.D0
  20 :         SZ=0.D0
  21 :         SXY=0.D0
  22 :         SXZ=0.D0
  23 :         SYZ=0.D0
  24 : C
  25 :       GO TO (10,1000,30,40,50,60),NSTRS
  26 :       GOTO 1000
  27 : C
  28 :    60 CONTINUE
  29 : C
  30 : C     COMPORTEMENT TRIDIMENSIONNEL
  31 : C
  32 :         SX =S1(1)
  33 :         SY =S1(2)
  34 :         SZ =S1(3)
  35 :         SXY=S1(4)
  36 :         SXZ=S1(5)
  37 :         SYZ=S1(6)
  38 :         GO TO 100
  39 : C
  40 :    50 CONTINUE
  41 : C
  42 : C     CONTRAINTES OU DEFORMATIONS PLANES
  43 : C        AVEC TOUS LES CISAILLEMENTS
  44 : C
  45 :         SX =S1(1)
  46 :         SY =S1(2)
  47 :         SXY=S1(3)
  48 :         SXZ=S1(4)
  49 :         SYZ=S1(5)
  50 :         GO TO 100
  51 : C
  52 :    40 CONTINUE
  53 : C
  54 : C       CAS AXISYMETRIQUE
  55 : C
  56 :         SX =S1(1)
  57 :         SY =S1(2)
  58 :         SZ =S1(3)
  59 :         SXY=S1(4)
  60 :         GO TO 100
  61 : C
  62 :    30 CONTINUE
  63 : C
  64 : C     CONTRAINTES OU DEFORMATIONS PLANES
  65 : C        SANS CISAILLEMENT TRANSVERSAL
  66 : C
  67 :         SX =S1(1)
  68 :         SY =S1(2)
  69 :         SXY=S1(3)
  70 :         GO TO 100
  71 : C
  72 :    10 CONTINUE
  73 : C
  74 : C     COMPORTEMENT UNIDIMENSIONNEL
  75 : C
  76 :         SX =S1(1)
  77 :         GO TO 100
  78 : C
  79 :  1000 WRITE(*,200) NSTRS
  80 :       STOP
  81 :   200 FORMAT(//,5X,' CAS NON DISPONIBLE DANS PARINV',/,
  82 :      * 20X,' NSTRS=',I2,//)
  83 :   100 CONTINUE
  84 : C
  85 :       RETURN
  86 :       END
  87 : 
  88 :  
  89 :