Télécharger msheti.eso

Retour à la liste

Numérotation des lignes :

  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
  108. C2003 CALL MSHEAR(XDDEP,XFOR0,XFORF,VAR0,VARF,NVARI,DDINL,
  109. C2003 & DTRANP,DTRANM,BETA,NPELA,
  110. C2004 CALL MSHEAR(XDDEP,XFOR0,XFORF,VAR0,VARF,NVARI,DDINL,
  111. C2004& DTRANP,DTRANM,BETA,NPELA,TRFA,DOCP,DOCM
  112. CALL MSHEAR(XDDEP,XFOR0,XFORF,VAR0,VARF,NVARI,DDINL,
  113. & DTRANP,DTRANM,BETA,NPELA,TRFA,DOCP,DOCM,EXPN,
  114. & TRAC(IPOSFP),NCURFP,TRAC(IPOSKP),NCURKP,TRAC(IPOSLP),NCURLP,
  115. & TRAC(IPOSFM),NCURFM,TRAC(IPOSKM),NCURKM,TRAC(IPOSLM),NCURLM,
  116. & KERRE)
  117. C
  118. C==================================================================
  119. C
  120. C On remplit les tableaux avant de sortir
  121. C
  122. C==================================================================
  123. NSTRS=DEFP(/1)
  124. DO IE1=1,NSTRS
  125. DEFP(IE1)=XZER
  126. ENDDO
  127. IF (IDIM.EQ.3) THEN
  128. SIGF(1) = SIG0(1) + SECT*DEPST(1)
  129. SIGF(2) = SIG0(2) + SRDY*DEPST(2)
  130. SIGF(3) = XFORF
  131. DEFP(3) = DDINL
  132. SIGF(4) = SIG0(4) + INRX*DEPST(4)
  133. SIGF(5) = SIG0(5) + INRY*DEPST(5)
  134. SIGF(6) = SIG0(6) + INRZ*DEPST(6)
  135. ELSE
  136. SIGF(1) = SIG0(1) + SECT*DEPST(1)
  137. SIGF(2) = XFORF
  138. DEFP(2) = DDINL
  139. SIGF(3) = SIG0(3) + INRZ*DEPST(3)
  140. ENDIF
  141. C
  142. C=====================================================
  143. C FIN DE LA ROUTINE DU MODELE
  144. C======================================================
  145. RETURN
  146.  
  147. END
  148.  
  149.  
  150.  
  151.  
  152.  

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