Esope source files

   1 : C MSHETI    SOURCE    CHAT      05/01/13    01:53:08     5004           
   2 : C MSHETI    SOURCE    PP        02/01/10    23:59:59     ????
   3 : C======================================================================
   4 : C       MUR_SHEAR, Elisa, Armelle et Pierre, ELSA/ISPRA 01/2002
   5 : C======================================================================
   6 : C
   7 : C                    MODELE GLOBAL DE MUR EN CISAILLEMENT
   8 : C            (Sur des elements de poutre TIMO - Effort tranchant/Cisail.)
   9 : C            (pour appel de ecou60)
  10 : C
  11 :       SUBROUTINE MSHETI(WRK0,WRK1,WRK2,
  12 :      >                  NCURFP,NCURKP,NCURLP,NCURFM,NCURKM,NCURLM,
  13 :      >                  IPOSFP,IPOSKP,IPOSLP,IPOSFM,IPOSKM,IPOSLM,
  14 :      >                  KERRE)
  15 :       IMPLICIT INTEGER(I-N)
  16 :       IMPLICIT REAL*8(A-H,O-Z)
  17 : C
  18 : C=======================================================================
  19 : C             CETTE ROUTINE EST APPELE DANS ECOU60
  20 : C
  21 : C
  22 : C WRK0      = Segment materiaux
  23 : C WRK1      = Segment contraintes
  24 : C WRK2      = Segment courbe
  25 : C NCURFP    = longueur courbe de charge  +           (x>0,y>0)
  26 : C NCURKP    = longueur courbe de raideur +           (x>0,y>0)
  27 : C NCURLP    = longueur courbe de domaine elastique + (x>0,y>0)
  28 : C NCURFM    = longueur courbe de charge  -           (x<0,y<0)
  29 : C NCURKM    = longueur courbe de raideur -           (x<0,y>0)
  30 : C NCURLM    = longueur courbe de domaine elastique - (x<0,y>0)
  31 : C IPOS....  = pointeurs sur les courbes dans TRAC
  32 : C
  33 : C======================================================================
  34 : C XCAR      = Caracteristique de la section
  35 : C DEPST     = Increment de deformation axiale
  36 : C SIG0      = Contrainte initiale
  37 : C VAR0      = Variables internes initiales
  38 : C SIGF      = Contrainte finale
  39 : C VARF      = Variables internes finales
  40 : C DEFP      = Deformation plastique
  41 : C
  42 : C=======================================================================
  43 : C
  44 :       REAL*8 INRX,INRY,INRZ
  45 :       PARAMETER (XZER=0.D0,UN=1.D0,EPSILO=1.D-16)
  46 : C
  47 : C=======================================================================
  48 : C      VARIABLES ET SEGMENTS NECESSAIRES
  49 : C=========================================================================
  50 : -INC CCOPTIO
  51 : C  Segment des materiaux
  52 :        SEGMENT WRK0
  53 :          REAL*8 XMAT(NCXMAT)
  54 :        ENDSEGMENT
  55 : C  Segment des contraintes
  56 :        SEGMENT WRK1
  57 :          REAL*8 DDHOOK(LHOOK,LHOOK),SIG0(NSTRS),DEPST(NSTRS)
  58 :          REAL*8 SIGF(NSTRS),VAR0(NVARI),VARF(NVARI)
  59 :          REAL*8 DEFP(NSTRS),XCAR(ICARA)
  60 :        ENDSEGMENT
  61 : C  Segment de la courbe
  62 :        SEGMENT WRK2
  63 :          REAL*8 TRAC(LTRAC)
  64 :        ENDSEGMENT
  65 : C
  66 :        NVARI=VAR0(/1)
  67 : C
  68 : C  Lecture materiau
  69 : C
  70 :        YOUN   = XMAT(1)
  71 :        XNU    = XMAT(2)
  72 :        DTRANP = XMAT(11)
  73 :        DTRANM = XMAT(12)
  74 :        BETA   = XMAT(13)
  75 :        NPELA  = INT(XMAT(14)+EPSILO)
  76 : c+2003
  77 :        TRFA   = XMAT(15)
  78 :        DOCP   = XMAT(16)
  79 :        DOCM   = XMAT(17)
  80 : C+2003
  81 : c+2004
  82 :        EXPN   = XMAT(18)
  83 : C+2004
  84 : C
  85 :        IF (IDIM.EQ.3) THEN
  86 :         INRX = YOUN*XCAR(1)
  87 :         INRY = YOUN*XCAR(2)
  88 :         INRZ = YOUN*XCAR(3)
  89 :         SECT = YOUN*XCAR(4)
  90 :         SRDY = (YOUN/(2.D0*(1.D0+XNU)))*XCAR(5)
  91 :         SECZ = (YOUN/(2.D0*(1.D0+XNU)))*XCAR(6)
  92 : C
  93 :         XDDEP = DEPST(3)
  94 :         XFOR0 = SIG0(3)
  95 : C
  96 :        ELSE
  97 : C
  98 :         SECT = YOUN*XCAR(1)
  99 :         INRZ = YOUN*XCAR(2)
 100 :         SECZ = (YOUN/(2.D0*(1.D0+XNU)))*XCAR(3)
 101 : C
 102 :         XDDEP = DEPST(2)
 103 :         XFOR0 = SIG0(2)
 104 :        ENDIF
 105 : C
 106 : C Appel de la loi
 107 : C
 107 : C2003 CALL MSHEAR(XDDEP,XFOR0,XFORF,VAR0,VARF,NVARI,DDINL,
 108 : C2003 & DTRANP,DTRANM,BETA,NPELA,
 108 : C2004  CALL MSHEAR(XDDEP,XFOR0,XFORF,VAR0,VARF,NVARI,DDINL,
 109 : C2004& DTRANP,DTRANM,BETA,NPELA,TRFA,DOCP,DOCM
 109 :        CALL MSHEAR(XDDEP,XFOR0,XFORF,VAR0,VARF,NVARI,DDINL,
 110 :      & DTRANP,DTRANM,BETA,NPELA,TRFA,DOCP,DOCM,EXPN,
 111 :      & TRAC(IPOSFP),NCURFP,TRAC(IPOSKP),NCURKP,TRAC(IPOSLP),NCURLP,
 112 :      & TRAC(IPOSFM),NCURFM,TRAC(IPOSKM),NCURKM,TRAC(IPOSLM),NCURLM,
 113 :      & KERRE)
 114 : C
 115 : C==================================================================
 116 : C
 117 : C  On remplit les tableaux avant de sortir
 118 : C
 119 : C==================================================================
 120 :        NSTRS=DEFP(/1)
 121 :        DO IE1=1,NSTRS
 122 :          DEFP(IE1)=XZER
 123 :        ENDDO
 124 :        IF (IDIM.EQ.3) THEN
 125 :         SIGF(1) = SIG0(1) + SECT*DEPST(1)
 126 :         SIGF(2) = SIG0(2) + SRDY*DEPST(2)
 127 :         SIGF(3) = XFORF
 128 :         DEFP(3) = DDINL
 129 :         SIGF(4) = SIG0(4) + INRX*DEPST(4)
 130 :         SIGF(5) = SIG0(5) + INRY*DEPST(5)
 131 :         SIGF(6) = SIG0(6) + INRZ*DEPST(6)
 132 :        ELSE
 133 :         SIGF(1) = SIG0(1) + SECT*DEPST(1)
 134 :         SIGF(2) = XFORF
 135 :         DEFP(2) = DDINL
 136 :         SIGF(3) = SIG0(3) + INRZ*DEPST(3)
 137 :        ENDIF
 138 : C
 139 : C=====================================================
 140 : C            FIN DE LA ROUTINE DU MODELE
 141 : C======================================================
 142 :       RETURN
 143 : 
 144 :       END
 145 : 
 146 : 
 147 :  
 148 :