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