Télécharger ricbet.eso

Retour à la liste

Numérotation des lignes :

ricbet
  1. C RICBET SOURCE PV 22/01/23 21:15:05 11271
  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.  
  32. -INC PPARAM
  33. -INC CCOPTIO
  34. -INC DECHE
  35. C
  36. C-----DEFINITION DE L ISOTROPIE-----------------------------------------
  37. C
  38. CMATE = 'ISOTROPE'
  39. C
  40. C-----PARAMETERES-------------------------------------------------------
  41. C
  42.  
  43. C Dans le cas d'un calcul local :
  44. C --> NSUB = NSUB0
  45. C
  46. C Dans le cas d'un calculs non local :
  47. C --> calcul du seuil : NSUB = 1
  48. C --> calcul moyenne : NSUB = NSUB0
  49.  
  50.  
  51. C Choix de NSUB0 de maniere adaptative
  52. XTEMP = 0.0D0
  53.  
  54. DO I=1,(2*IDIM)
  55. XTEMP = MAX(ABS(DEPST(I)),XTEMP)
  56. ENDDO
  57.  
  58. IF (XTEMP.LE.1.0D-5) THEN
  59. NSUB0 = 10
  60. ELSEIF ((XTEMP.GT.1.0D-5).AND.(XTEMP.LE.1.0D-4)) THEN
  61. NSUB0 = 100
  62. ELSEIF ((XTEMP.GT.1.0D-4).AND.(XTEMP.LE.1.0D-3)) THEN
  63. NSUB0 = 1000
  64. ELSEIF (XTEMP.GT.1.0D-3) THEN
  65. NSUB0 = 1000
  66. ENDIF
  67.  
  68. IF (ISTEP.EQ.0) THEN
  69. NSUB = NSUB0
  70. ELSEIF (ISTEP.EQ.1) THEN
  71. NSUB = 1
  72. ELSEIF (ISTEP.EQ.2) THEN
  73. NSUB = NSUB0
  74. ENDIF
  75.  
  76. C
  77. C-----SWITCH SELON LA DIMENSION-----------------------------------------
  78. C
  79. IF (IDIM.EQ.2) THEN
  80.  
  81. PRINT*,'Pas operationnel... a travailler'
  82. STOP
  83.  
  84. C RICBET2D
  85. C-----------------------------------------------------------------------
  86. C
  87.  
  88. C On forme le vecteur des increments de deformations subdivisees
  89. DO K=1,4
  90. DEPST2(K)=DEPST(K)/NSUB
  91. ENDDO
  92.  
  93. C On retient les varibles internes de depart
  94. DO K=1,42
  95. VAR1(K)=VAR0(K)
  96. ENDDO
  97.  
  98. C On commence le substepping
  99. DO J=1,NSUB
  100.  
  101. ** CALL RICBET2D(XMAT,SIG0,SIGF,VAR1,VARF,DEPST2,ISTEP)
  102.  
  103. DO I=1,42
  104. VAR1(I)=VARF(I)
  105. ENDDO
  106.  
  107. ENDDO
  108.  
  109. ELSE
  110.  
  111. C RICBET3D
  112. C-----------------------------------------------------------------------
  113. C
  114.  
  115. C On forme le vecteur des increments de deformations subdivisees
  116. DO K=1,6
  117. DEPST3(K)=DEPST(K)/NSUB
  118. ENDDO
  119.  
  120. C On retient les varibles internes de depart
  121. DO K=1,42
  122. VAR1(K)=VAR0(K)
  123. ENDDO
  124.  
  125. C On commence le substepping
  126. DO J=1,NSUB
  127.  
  128. CALL RICBET3D(XMAT,SIG0,SIGF,VAR1,VARF,DEPST3,ISTEP)
  129.  
  130. C CALL RICBET3D2(XMAT,SIG0,SIGF,VAR1,VARF,DEPST3,ISTEP)
  131.  
  132. DO I=1,42
  133. VAR1(I)=VARF(I)
  134. ENDDO
  135.  
  136. ENDDO
  137.  
  138. ENDIF
  139.  
  140. RETURN
  141. END
  142.  
  143.  
  144.  
  145.  
  146.  
  147.  
  148.  
  149.  
  150.  
  151.  
  152.  
  153.  
  154.  
  155.  
  156.  
  157.  
  158.  
  159.  
  160.  
  161.  
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  
  168.  

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