Télécharger la plaquette Cast3M


fr en

Télécharger parinv.eso

Retour à la liste

Numérotation des lignes : 

parinv

  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
  14.         CALL ZERO(S1,6,1)
  15.         DO 11 I=1,NSTRS
  16.           S1(I)=S(I)
  17.    11   CONTINUE
  18. C
  19.         SX=0.D0
  20.         SY=0.D0
  21.         SZ=0.D0
  22.         SXY=0.D0
  23.         SXZ=0.D0
  24.         SYZ=0.D0
  25. C
  26.       GO TO (10,1000,30,40,50,60),NSTRS
  27.       GOTO 1000
  28. C
  29.    60 CONTINUE
  30. C
  31. C     COMPORTEMENT TRIDIMENSIONNEL
  32. C
  33.         SX =S1(1)
  34.         SY =S1(2)
  35.         SZ =S1(3)
  36.         SXY=S1(4)
  37.         SXZ=S1(5)
  38.         SYZ=S1(6)
  39.         GO TO 100
  40. C
  41.    50 CONTINUE
  42. C
  43. C     CONTRAINTES OU DEFORMATIONS PLANES
  44. C        AVEC TOUS LES CISAILLEMENTS
  45. C
  46.         SX =S1(1)
  47.         SY =S1(2)
  48.         SXY=S1(3)
  49.         SXZ=S1(4)
  50.         SYZ=S1(5)
  51.         GO TO 100
  52. C
  53.    40 CONTINUE
  54. C
  55. C       CAS AXISYMETRIQUE
  56. C
  57.         SX =S1(1)
  58.         SY =S1(2)
  59.         SZ =S1(3)
  60.         SXY=S1(4)
  61.         GO TO 100
  62. C
  63.    30 CONTINUE
  64. C
  65. C     CONTRAINTES OU DEFORMATIONS PLANES
  66. C        SANS CISAILLEMENT TRANSVERSAL
  67. C
  68.         SX =S1(1)
  69.         SY =S1(2)
  70.         SXY=S1(3)
  71.         GO TO 100
  72. C
  73.    10 CONTINUE
  74. C
  75. C     COMPORTEMENT UNIDIMENSIONNEL
  76. C
  77.         SX =S1(1)
  78.         GO TO 100
  79. C
  80.  1000 WRITE(*,200) NSTRS
  81.       STOP
  82.   200 FORMAT(//,5X,' CAS NON DISPONIBLE DANS PARINV',/,
  83.      * 20X,' NSTRS=',I2,//)
  84.   100 CONTINUE
  85. C
  86.       RETURN
  87.       END
  88.  
  89.  
  90.  
  91.  

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