Télécharger la plaquette Cast3M


fr en

Télécharger msheti.eso

Retour à la liste

Numérotation des lignes : 

msheti

  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.  
  51. -INC PPARAM
  52. -INC CCOPTIO
  53. C  Segment des materiaux
  54.        SEGMENT WRK0
  55.          REAL*8 XMAT(NCXMAT)
  56.        ENDSEGMENT
  57. C  Segment des contraintes
  58.        SEGMENT WRK1
  59.          REAL*8 DDHOOK(LHOOK,LHOOK),SIG0(NSTRS),DEPST(NSTRS)
  60.          REAL*8 SIGF(NSTRS),VAR0(NVARI),VARF(NVARI)
  61.          REAL*8 DEFP(NSTRS),XCAR(ICARA)
  62.        ENDSEGMENT
  63. C  Segment de la courbe
  64.        SEGMENT WRK2
  65.          REAL*8 TRAC(LTRAC)
  66.        ENDSEGMENT
  67. C
  68.        NVARI=VAR0(/1)
  69. C
  70. C  Lecture materiau
  71. C
  72.        YOUN   = XMAT(1)
  73.        XNU    = XMAT(2)
  74.        DTRANP = XMAT(11)
  75.        DTRANM = XMAT(12)
  76.        BETA   = XMAT(13)
  77.        NPELA  = INT(XMAT(14)+EPSILO)
  78. c+2003
  79.        TRFA   = XMAT(15)
  80.        DOCP   = XMAT(16)
  81.        DOCM   = XMAT(17)
  82. C+2003
  83. c+2004
  84.        EXPN   = XMAT(18)
  85. C+2004
  86. C
  87.        IF (IDIM.EQ.3) THEN
  88.         INRX = YOUN*XCAR(1)
  89.         INRY = YOUN*XCAR(2)
  90.         INRZ = YOUN*XCAR(3)
  91.         SECT = YOUN*XCAR(4)
  92.         SRDY = (YOUN/(2.D0*(1.D0+XNU)))*XCAR(5)
  93.         SECZ = (YOUN/(2.D0*(1.D0+XNU)))*XCAR(6)
  94. C
  95.         XDDEP = DEPST(3)
  96.         XFOR0 = SIG0(3)
  97. C
  98.        ELSE
  99. C
  100.         SECT = YOUN*XCAR(1)
  101.         INRZ = YOUN*XCAR(2)
  102.         SECZ = (YOUN/(2.D0*(1.D0+XNU)))*XCAR(3)
  103. C
  104.         XDDEP = DEPST(2)
  105.         XFOR0 = SIG0(2)
  106.        ENDIF
  107. C
  108. C Appel de la loi
  109. C
  110. C2003 CALL MSHEAR(XDDEP,XFOR0,XFORF,VAR0,VARF,NVARI,DDINL,
  111. C2003 & DTRANP,DTRANM,BETA,NPELA,
  112. C2004  CALL MSHEAR(XDDEP,XFOR0,XFORF,VAR0,VARF,NVARI,DDINL,
  113. C2004& DTRANP,DTRANM,BETA,NPELA,TRFA,DOCP,DOCM
  114.        CALL MSHEAR(XDDEP,XFOR0,XFORF,VAR0,VARF,NVARI,DDINL,
  115.      & DTRANP,DTRANM,BETA,NPELA,TRFA,DOCP,DOCM,EXPN,
  116.      & TRAC(IPOSFP),NCURFP,TRAC(IPOSKP),NCURKP,TRAC(IPOSLP),NCURLP,
  117.      & TRAC(IPOSFM),NCURFM,TRAC(IPOSKM),NCURKM,TRAC(IPOSLM),NCURLM,
  118.      & KERRE)
  119. C
  120. C==================================================================
  121. C
  122. C  On remplit les tableaux avant de sortir
  123. C
  124. C==================================================================
  125.        NSTRS=DEFP(/1)
  126.        DO IE1=1,NSTRS
  127.          DEFP(IE1)=XZER
  128.        ENDDO
  129.        IF (IDIM.EQ.3) THEN
  130.         SIGF(1) = SIG0(1) + SECT*DEPST(1)
  131.         SIGF(2) = SIG0(2) + SRDY*DEPST(2)
  132.         SIGF(3) = XFORF
  133.         DEFP(3) = DDINL
  134.         SIGF(4) = SIG0(4) + INRX*DEPST(4)
  135.         SIGF(5) = SIG0(5) + INRY*DEPST(5)
  136.         SIGF(6) = SIG0(6) + INRZ*DEPST(6)
  137.        ELSE
  138.         SIGF(1) = SIG0(1) + SECT*DEPST(1)
  139.         SIGF(2) = XFORF
  140.         DEFP(2) = DDINL
  141.         SIGF(3) = SIG0(3) + INRZ*DEPST(3)
  142.        ENDIF
  143. C
  144. C=====================================================
  145. C            FIN DE LA ROUTINE DU MODELE
  146. C======================================================
  147.       RETURN
  148.  
  149.       END
  150.  
  151.  
  152.  
  153.  
  154.  

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