Esope source files

   1 : C CHRP3D    SOURCE    CHAT      07/10/22    21:15:32     5921           
   2 :       SUBROUTINE CHRP3D(PPP,SIEPIN,SIEPOO,IOP)
   3 :       IMPLICIT REAL*8(A-H,o-Z)
   4 :       IMPLICIT INTEGER (I-N)
   5 : *
   6 : *     ------------------------------------------------------------------
   7 : *
   8 : *         CHGT DE REPERE POUR LES CONTRAINTES ET LES DEFORMATIONS
   9 : *           ECRITES SOUS FORME VECTORIELLE EN 3D
  10 : *                                               H. BUNG      02-93
  11 : *     ------------------------------------------------------------------
  12 : *
  13 : *     PPP      : MATRICE DE PASSAGE
  14 : *     SIEPIN   : VECTEUR DE CONTRAINTES OU DE DEFORMATIONS EN ENTREE
  15 : *     SIEPOO   : VECTEUR DE CONTRAINTES OU DE DEFORMATIONS EN SORTIE
  16 : *     IOP            ENTREE     SORTIE 
  17 : *          = 0       SIG_G      SIG_L
  18 : *          = 1       SIG_L      SIG_G
  19 : *          = 2       EPS_G      EPS_L
  20 : *          = 3       EPS_L      EPS_G
  21 : *
  22 : *--    RAPPEL  
  23 : *       [SIG] = [ S_11, S_22, S_33,   S_12,   S_23,   S_13 ]
  24 : *       [EPS] = [ E_11, E_22, E_33, 2*E_12, 2*E_23, 2*E_13 ]
  25 : *
  26 : *
  27 : *      IMPLICIT NONE
  28 : *
  29 : *
  30 : *---   VARIABLES GLOBALES
  31 : *      REAL *8 PPP(3,3),SIEPIN(*),SIEPOO(*)
  32 : *      INTEGER IOP
  33 :        DIMENSION PPP(3,*),SIEPIN(*),SIEPOO(*)
  34 : *
  35 : *---   VARIABLES LOCALES
  36 :       
  37 :       DIMENSION  AIN(9),AOO(9)
  38 :       DIMENSION LLM(9),LLV(6),LDEF(6)
  39 : *
  40 :       DATA LLM/1,4,6, 4,2,5, 6,5,3/
  41 :       DATA LLV/1,5,9, 2,6,3/
  42 :       DATA LDEF/2,3, 4,6, 7,8/
  43 : *
  44 : *
  45 : *----    PASSAGE VERS LA FORME MATRICIELLE
  46 :       DO IK=1,9
  47 :          AIN(IK) = SIEPIN(LLM(IK))
  48 :       END DO
  49 : *
  50 :       IF(IOP.GE.2) THEN
  51 : *---     TENSEUR DE DEFORMATIONS 
  52 :          DO IK=1,6
  53 :             IKK      = LDEF(IK)
  54 :             AIN(IKK) = 0.5D0*AIN(IKK)
  55 :          ENDDO
  56 :       ENDIF
  57 : *
  58 :       KOP=MOD(IOP,2) 
  59 :       IF(KOP.EQ.0)THEN
  60 : *          ENTREE : GLOBAL --> SORTIE : LOCAL
  60 :          CALL DR3GLM1(PPP,AIN,AOO)
  61 :       ENDIF
  62 :       IF(KOP.EQ.1) THEN
  63 : *          ENTREE : LOCAL --> SORTIE : GLOBAL
  63 :          CALL DR3GLM2(PPP,AOO,AIN)
  64 :       ENDIF
  65 :       IF(KOP.NE.0.AND.KOP.NE.1)THEN
  66 : *          STOP
  67 :          STOP 'CHR3D_1'
  68 :       ENDIF
  69 : C      SELECT CASE(KOP)
  70 : C      CASE (0) 
  71 : C*          ENTREE : GLOBAL --> SORTIE : LOCAL
  71 : C         CALL DR3GLM1(PPP,AIN,AOO)
  72 : C      CASE (1)
  73 : C*          ENTREE : LOCAL --> SORTIE : GLOBAL
  73 : C         CALL DR3GLM2(PPP,AOO,AIN)
  74 : C      CASE DEFAULT
  75 : C*          STOP
  76 : C         STOP 'CHR3D_1'
  77 : C      END SELECT
  78 : *
  79 : *----    PASSAGE VERS LA FORME VECTORIELLE
  80 :       DO IK=1,6
  81 :          SIEPOO(IK) = AOO(LLV(IK))
  82 :       END DO
  83 : *
  84 : *----    POUR LES DEFORMATIONS
  85 :       IF(IOP.GE.2) THEN
  86 :          SIEPOO(4) = 2.D0*SIEPOO(4)
  87 :          SIEPOO(5) = 2.D0*SIEPOO(5)
  88 :          SIEPOO(6) = 2.D0*SIEPOO(6)
  89 :       ENDIF
  90 : *
  91 :       END
  92 : 
  93 :