Esope source files

   1 : C MASHEJ    SOURCE    FANDEUR   12/07/18    21:15:50     7435           
   2 :       SUBROUTINE MASHEJ(wrk52,wrk53,WRK2,
   3 :      >                  NCURFP,NCURKP,NCURLP,NCURFM,NCURKM,NCURLM,
   4 :      >                  IPOSFP,IPOSKP,IPOSLP,IPOSFM,IPOSKM,IPOSLM)
   5 :       IMPLICIT INTEGER(I-N)
   6 :       IMPLICIT REAL*8(A-H,O-Z)
   7 : C
   8 : C======================================================================
   9 : C       MUR_SHEAR, Alessandra, Armelle  et Pierre, ELSA/ISPRA 
  10 : C                  05/2004
  11 : C======================================================================
  12 : C
  13 : C                    MODELE GLOBAL DE MUR EN CISAILLEMENT +
  14 : C                           COMPORTEMENT DE FLEXION EN TETE ET PIED
  15 : C            (Sur des elements macro  CIFL )
  16 : C            (pour appel de coml8)
  17 : C
  18 : C=======================================================================
  19 : C             CETTE ROUTINE EST APPELE DANS coml8
  20 : C
  21 : C WRK52     = Segment materiaux/contrainte/....
  22 : C WRK53     = Segment pratique (not used now)
  23 : C WRK2      = Segment courbe
  24 : C NCURFP    = longueur courbe de charge  +           (x>0,y>0)
  25 : C NCURKP    = longueur courbe de raideur +           (x>0,y>0)
  26 : C NCURLP    = longueur courbe de domaine elastique + (x>0,y>0)
  27 : C NCURFM    = longueur courbe de charge  -           (x<0,y<0)
  28 : C NCURKM    = longueur courbe de raideur -           (x<0,y>0)
  29 : C NCURLM    = longueur courbe de domaine elastique - (x<0,y>0)
  30 : C IPOS....  = pointeurs sur les courbes dans TRAC
  31 : C
  32 : C======================================================================
  33 : C XCARB     = Caracteristique de la section
  34 : C DEPST     = Increment de deformation axiale
  35 : C SIG0      = Contrainte initiale
  36 : C VAR0      = Variables internes initiales
  37 : C SIGF      = Contrainte finale
  38 : C VARF      = Variables internes finales
  39 : C DEFP      = Deformation plastique
  40 : C
  41 : C=======================================================================
  42 : C
  43 :       REAL*8 INRZ,MOMT
  44 :       PARAMETER (XZER=0.D0,UN=1.D0,EPSILO=1.D-16)
  45 : C
  46 : C=======================================================================
  47 : C      VARIABLES ET SEGMENTS NECESSAIRES
  48 : C=========================================================================
  49 : -INC CCOPTIO
  50 : -INC DECHE
  51 : C
  52 : C  Segment de la courbe
  53 :        SEGMENT WRK2
  54 :          REAL*8 TRAC(LTRAC)
  55 :        ENDSEGMENT
  56 : C
  57 :        NVARI=VAR0(/1)
  58 : C
  59 : C  Lecture materiau
  60 : C
  61 :        YOUN   = XMAT(1)
  62 :        XNU    = XMAT(2)
  63 :        DTRANP = XMAT(11)
  64 :        DTRANM = XMAT(12)
  65 :        BETA   = XMAT(13)
  66 :        NPELA  = INT(XMAT(14)+EPSILO)
  67 :        TRFA   = XMAT(15)
  68 :        DOCP   = XMAT(16)
  69 :        DOCM   = XMAT(17)
  70 :        EXPN   = XMAT(18)
  71 :        FNRF   = XMAT(19)
  72 :        FTDI   = XMAT(20)
  73 : C
  74 :        SECT   = XCARB(1)
  75 :        INRZ   = XCARB(2)
  76 :        BB     = 2*SQRT(3*INRZ/SECT)
  77 :        SS     = SECT/BB
  78 : C
  79 : C  PARTIE INFERIEURE DU MUR
  80 : C
  81 :        EPT   = epstf(1)
  82 :        COT   = epstf(2)
  83 :        IF(EPT*EPT.GE.(COT*BB/2)**2)THEN
  84 :          FNOT=0.5D0*YOUN*SECT*(EPT-ABS(EPT))
  85 :          MOMT=0.5D0*YOUN*INRZ*(UN-SIGN(UN,EPT))*COT
  86 : *        write(6,*)'1:',FNOT,MOMT
  87 :        ELSE
  88 :          FNOT=YOUN*SS/2*((BB*EPT)-(EPT*EPT/ABS(COT))-(ABS(COT)*BB*BB/4))
  89 :          MOMT=YOUN*SS/24*
  90 :      >       (((BB**3)*COT)+
  91 :      >        (SIGN(UN,COT)*((4*(EPT**3)/(COT**2))-(3*BB*BB*EPT))))
  92 : *        write(6,*)'2:',FNOT,MOMT
  93 :        ENDIF
  94 :        SIGF(1) = FNOT
  95 :        SIGF(2) = MOMT 
  96 : C
  97 : C  PARTIE SUPERIEURE DU MUR
  98 : C
  99 :        EPT   = epstf(4)
 100 :        COT   = epstf(5)
 101 :        IF(EPT*EPT.GE.(COT*BB/2)**2)THEN
 102 :          FNOT=YOUN*SECT/2*(EPT-ABS(EPT))
 103 :          MOMT=YOUN*INRZ/2*(UN-SIGN(UN,EPT))*COT
 104 :        ELSE
 105 :          FNOT=YOUN*SS/2*((BB*EPT)-(EPT*EPT/ABS(COT))-(ABS(COT)*BB*BB/4))
 106 :          MOMT=YOUN*SS/24*
 107 :      >       (((BB**3)*COT)+
 108 :      >        (SIGN(UN,COT)*((4*EPT**3/COT**2)-(3*BB*BB*EPT))))
 109 :        ENDIF
 110 :        SIGF(4) = FNOT
 111 :        SIGF(5) = MOMT 
 112 : C
 113 : C  PARTIE MEDIANE DU MUR
 114 : C
 115 :        XDDEP = DEPST(3)
 116 :        XFOR0 = SIG0(3)
 117 : C
 118 : C  Appel de la loi de shear
 119 : C
 119 :        CALL MSHEAR(XDDEP,XFOR0,XFORF,VAR0,VARF,NVARI,DDINL,
 120 :      & DTRANP,DTRANM,BETA,NPELA,TRFA,DOCP,DOCM,EXPN,
 121 :      & TRAC(IPOSFP),NCURFP,TRAC(IPOSKP),NCURKP,TRAC(IPOSLP),NCURLP,
 122 :      & TRAC(IPOSFM),NCURFM,TRAC(IPOSKM),NCURKM,TRAC(IPOSLM),NCURLM,
 123 :      & KERRE)
 124 : C
 125 : C  RENORMALISATION
 126 : C
 127 :        SIGF(3) = XFORF*SQRT((SECT*FTDI+(ABS(SIGF(1)+SIGF(4))/2))/
 128 :      >                      (SECT*FTDI+          FNRF          )  )
 129 : C      write(6,*)SIGF(3)/XFORF,FNRF,ABS(SIGF(1)+SIGF(4))/2
 130 : C
 131 : C==================================================================
 132 : C
 133 : C  On remplit les tableaux avant de sortir
 134 : C
 135 : C==================================================================
 136 :        NSTRS=DEFP(/1)
 137 :        DEFP(1)=DEPST(1)-((SIGF(1)-SIG0(1))/(YOUN*SECT))
 138 :        DEFP(2)=DEPST(2)-((SIGF(2)-SIG0(2))/(YOUN*INRZ))
 139 :        DEFP(3)=DDINL
 140 :        DEFP(4)=DEPST(4)-((SIGF(4)-SIG0(4))/(YOUN*SECT))
 141 :        DEFP(5)=DEPST(5)-((SIGF(5)-SIG0(5))/(YOUN*INRZ))
 142 : C
 143 : C=====================================================
 144 : C            FIN DE LA ROUTINE DU MODELE
 145 : C======================================================
 146 :       RETURN
 147 :       END
 148 : 
 149 :