Télécharger ricbet.eso

Retour à la liste

Numérotation des lignes :

  1. C RICBET SOURCE BP208322 17/03/01 21:18:06 9325
  2. SUBROUTINE RICBET(wrk52,wrk53,wrk54,nvari,iecou)
  3. C
  4. C====&===1=========2=========3=========4=========5=========6=========7==
  5. C Commentaires : Subroutine permettant de mettre en oeuvre le
  6. C modele RICBET (B. Richard) pour representer
  7. C le comportement 2D/3D d'un beton ordinaire sous
  8. C sollicitations cycliques
  9. C
  10. C Traits : - Endommagement scalaire
  11. C - Boucles hystŽrŽsis
  12. C - DŽformations permanentes
  13. C - Refermeture lineaire des fissures
  14. C - Effet unilatteral complet
  15. C - Non localite portant sur le seuil d'endommagement (Fd)
  16. C
  17. C Auteur : B. Richard
  18. C
  19. C Co-auteur : F. Ragueneau
  20. C
  21. C Date : 2010 - 2011
  22. C====&===1=========2=========3=========4=========5=========6=========7==
  23. C
  24. C-----DECLARATION GENERALE----------------------------------------------
  25. C
  26. IMPLICIT REAL*8(A-H,O-Z)
  27. REAL*8 DEPST2(4),DEPST3(6),VAR1(42)
  28. C
  29. C-----LISTE DES INCLUDES A CHARGER--------------------------------------
  30. C
  31. -INC CCOPTIO
  32. -INC DECHE
  33. C
  34. C-----DEFINITION DE L ISOTROPIE-----------------------------------------
  35. C
  36. CMATE = 'ISOTROPE'
  37. C
  38. C-----PARAMETERES-------------------------------------------------------
  39. C
  40.  
  41. C Dans le cas d'un calcul local :
  42. C --> NSUB = NSUB0
  43. C
  44. C Dans le cas d'un calculs non local :
  45. C --> calcul du seuil : NSUB = 1
  46. C --> calcul moyenne : NSUB = NSUB0
  47.  
  48.  
  49. C Choix de NSUB0 de maniere adaptative
  50. XTEMP = 0.0D0
  51.  
  52. DO I=1,(2*IDIM)
  53. XTEMP = MAX(ABS(DEPST(I)),XTEMP)
  54. ENDDO
  55.  
  56. IF (XTEMP.LE.1.0D-5) THEN
  57. NSUB0 = 10
  58. ELSEIF ((XTEMP.GT.1.0D-5).AND.(XTEMP.LE.1.0D-4)) THEN
  59. NSUB0 = 100
  60. ELSEIF ((XTEMP.GT.1.0D-4).AND.(XTEMP.LE.1.0D-3)) THEN
  61. NSUB0 = 1000
  62. ELSEIF (XTEMP.GT.1.0D-3) THEN
  63. NSUB0 = 1000
  64. ENDIF
  65.  
  66. IF (ISTEP.EQ.0) THEN
  67. NSUB = NSUB0
  68. ELSEIF (ISTEP.EQ.1) THEN
  69. NSUB = 1
  70. ELSEIF (ISTEP.EQ.2) THEN
  71. NSUB = NSUB0
  72. ENDIF
  73.  
  74. C
  75. C-----SWITCH SELON LA DIMENSION-----------------------------------------
  76. C
  77. IF (IDIM.EQ.2) THEN
  78.  
  79. PRINT*,'Pas operationnel... a travailler'
  80. STOP
  81.  
  82. C RICBET2D
  83. C-----------------------------------------------------------------------
  84. C
  85.  
  86. C On forme le vecteur des increments de deformations subdivisees
  87. DO K=1,4
  88. DEPST2(K)=DEPST(K)/NSUB
  89. ENDDO
  90.  
  91. C On retient les varibles internes de depart
  92. DO K=1,42
  93. VAR1(K)=VAR0(K)
  94. ENDDO
  95.  
  96. C On commence le substepping
  97. DO J=1,NSUB
  98.  
  99. CALL RICBET2D(XMAT,SIG0,SIGF,VAR1,VARF,DEPST2,ISTEP)
  100.  
  101. DO I=1,42
  102. VAR1(I)=VARF(I)
  103. ENDDO
  104.  
  105. ENDDO
  106.  
  107. ELSE
  108.  
  109. C RICBET3D
  110. C-----------------------------------------------------------------------
  111. C
  112.  
  113. C On forme le vecteur des increments de deformations subdivisees
  114. DO K=1,6
  115. DEPST3(K)=DEPST(K)/NSUB
  116. ENDDO
  117.  
  118. C On retient les varibles internes de depart
  119. DO K=1,42
  120. VAR1(K)=VAR0(K)
  121. ENDDO
  122.  
  123. C On commence le substepping
  124. DO J=1,NSUB
  125.  
  126. CALL RICBET3D(XMAT,SIG0,SIGF,VAR1,VARF,DEPST3,ISTEP)
  127.  
  128. C CALL RICBET3D2(XMAT,SIG0,SIGF,VAR1,VARF,DEPST3,ISTEP)
  129.  
  130. DO I=1,42
  131. VAR1(I)=VARF(I)
  132. ENDDO
  133.  
  134. ENDDO
  135.  
  136. ENDIF
  137.  
  138. RETURN
  139. END
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  
  147.  
  148.  
  149.  
  150.  
  151.  
  152.  
  153.  
  154.  
  155.  
  156.  
  157.  
  158.  
  159.  
  160.  
  161.  
  162.  
  163.  
  164.  

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