Télécharger eligeh.eso

Retour à la liste

Numérotation des lignes :

  1. C ELIGEH SOURCE PV 15/04/10 21:15:08 8474
  2. C
  3. SUBROUTINE ELIGEH(EPSO,G12,S1T,S2T,S3T,T1T,T3T,ALFA,EPPP,EPPM
  4. . ,EPPPF,EPPMF,XKTAN,TAUF)
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7. C----------------------------------------------------------------------
  8. C
  9. C LOI ELIGEHAUSEN 1D EXPRIMEE EN CONTRAINTE - GLISSEMENT
  10. C UTILISEE SUR DES ELEMENTS JOINTS 2D
  11. C ET POUR LES LOIS D'ANCRAGE (BARRE et FIBRE)
  12. C----------------------------------------------------------------------
  13. C Didier COMBESCURE - EMSI - MAI 2000
  14. C----------------------------------------------------------------------
  15. C
  16. C EPSO: Glissement total
  17. C G12: Module elastique (charge et décharge)
  18. C EPPP: Glissement plastique positif
  19. C EPPM: Glissement plastique négatif
  20. C EPPPF: Glissement plastique positif final
  21. C EPPMF: Glissement plastique négatif final
  22. C TAUF: Contrainte de cisaillement finale
  23. C KTAN: Module tangent
  24. C
  25. C S1T, S2T et S3T: Glissements caractéristiques
  26. C T1T et T3T: Contraintes de cisaillement caractéristiques
  27. C ALFA: Paramètre définissant la première partie de la courbe
  28. C
  29. C -----------------------------------
  30. C
  31. C VERIFICATION DES PARAMETRES
  32. C
  33. C
  34. -INC CCREEL
  35.  
  36. IF ((S1T.LT.0.D0).OR.(S2T.LT.0.D0).OR.(S3T.LT.0.D0).OR.
  37. . (T1T.LT.0.D0).OR.(T3T.LT.0.D0).OR.(S2T.LT.S1T).OR.
  38. . (S3T.LT.S2T).OR.(T1T.LT.T3T).OR.(ALFA.GT.1.D0).OR.
  39. . (ALFA.LT.0.D0)) THEN
  40. WRITE (*,*) 'Vos paramètres Eligehausen ne sont pas corrects'
  41. ENDIF
  42. C
  43. C
  44. C
  45. IF (EPSO.GE.EPPP) THEN
  46. C
  47. C CAS EPSO>EPPP : ELASTIQUE OU PLASTIQUE
  48. C
  49. C PREDICTION ELASTIQUE
  50. C
  51. TAUFELA = G12*(EPSO - EPPP)
  52. C
  53. C COURBE NON LINEAIRE
  54. C
  55. IF (EPSO.LT.S1T) THEN
  56. TAUFPLA = T1T*((EPSO/S1T)**ALFA)
  57. epso=max(xpetit*s1t,epso)
  58. XKTANP = ALFA*(T1T/S1T)*((EPSO/S1T)**(ALFA - 1.D0))
  59. ELSEIF (EPSO.LT.S2T) THEN
  60. TAUFPLA = T1T
  61. XKTANP = 0.D0
  62. ELSEIF (EPSO.LT.S3T) THEN
  63. TAUFPLA = T1T - ((T1T - T3T)*(EPSO - S2T)/(S3T - S2T))
  64. XKTANP = (T3T - T1T)/(S3T - S2T)
  65. ELSE
  66. TAUFPLA = T3T
  67. XKTANP = 0.D0
  68. ENDIF
  69. C
  70. C
  71. C
  72. IF (TAUFELA.GT.TAUFPLA) THEN
  73. TAUF = TAUFPLA
  74. EPPPF = EPSO - (TAUFPLA/G12)
  75. EPPMF = EPPM
  76. XKTAN = XKTANP
  77. ELSE
  78. TAUF = TAUFELA
  79. EPPPF = EPPP
  80. EPPMF = EPPM
  81. XKTAN = G12
  82. ENDIF
  83.  
  84. C
  85. ELSEIF (EPSO.GT.EPPM) THEN
  86. C
  87. C CAS EPPM<EPSO<EPPP : PHASE GLISSEMENT : TAUF + 0
  88. C
  89. TAUF = 0.D0
  90. XKTAN = 0.D0
  91. EPPPF = EPPP
  92. EPPMF = EPPM
  93. C
  94. ELSEIF (EPSO.LE.EPPM) THEN
  95. C
  96. C CAS EPSO<EPPM : ELASTIQUE OU PLASTIQUE
  97. C
  98. C PREDICTION ELASTIQUE
  99. C
  100. TAUFELA = G12*(EPSO - EPPM)
  101. C
  102. C COURBE NON LINEAIRE
  103. C
  104. IF (EPSO.GT.((-1.D0)*S1T)) THEN
  105. TAUFPLA = (-1.D0)*T1T*(((-1.D0)*EPSO/S1T)**ALFA)
  106. XKTANP = ALFA*(T1T/S1T)*(((-1.D0)*EPSO/S1T)**(ALFA - 1.D0))
  107. ELSEIF (EPSO.GT.((-1.D0)*S2T)) THEN
  108. TAUFPLA = (-1.D0)*T1T
  109. XKTANP = 0.D0
  110. ELSEIF (EPSO.GT.((-1.D0)*S3T)) THEN
  111. TAUFPLA = (-1.D0)*(T1T - ((T1T - T3T)*
  112. .(((-1.D0)*EPSO) - S2T)/(S3T - S2T)))
  113. XKTANP = (T3T - T1T)/(S3T - S2T)
  114. ELSE
  115. TAUFPLA = (-1.D0)*T3T
  116. XKTANP = 0.D0
  117. ENDIF
  118. C
  119. C
  120. C
  121. IF (TAUFELA.LT.TAUFPLA) THEN
  122. TAUF = TAUFPLA
  123. EPPMF = EPSO - (TAUFPLA/G12)
  124. EPPPF = EPPP
  125. XKTAN = XKTANP
  126. ELSE
  127. TAUF = TAUFELA
  128. EPPPF = EPPP
  129. EPPMF = EPPM
  130. XKTAN = G12
  131. ENDIF
  132.  
  133. ENDIF
  134. C
  135. C WRITE(*,*) XKTAN
  136. C
  137. C
  138. RETURN
  139. END
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  

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