Télécharger concyc1.eso

Retour à la liste

Numérotation des lignes :

  1. C CONCYC1 SOURCE PV 17/12/08 21:16:57 9660
  2. SUBROUTINE CONCYC1(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 CONCYC (M. Vassaux & 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 hysteretiques
  12. C - Deformations 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 : M. Vassaux
  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 PAS DE SUBSTEPPING, tant que l'on ne sait pas gérer ISTEP depuis le script...
  77. NSUB = 1
  78.  
  79. C RICBET3D
  80. C-----------------------------------------------------------------------
  81. C
  82.  
  83. C On forme le vecteur des increments de deformations subdivisees
  84. C DO K=1,6
  85. C DEPST3(K)=DEPST(K)/NSUB
  86. C ENDDO
  87.  
  88. C On retient les varibles internes de depart
  89. C DO K=1,35
  90. C VAR1(K)=VAR0(K)
  91. C ENDDO
  92.  
  93. C On commence le substepping
  94. C DO J=1,NSUB
  95.  
  96. CALL CONCYC2(XMAT,SIG0,SIGF,VAR0,VARF,DEPST,ISTEP)
  97.  
  98.  
  99. C DO I=1,35
  100. C VAR1(I)=VARF(I)
  101. C ENDDO
  102. C
  103. C ENDDO
  104.  
  105.  
  106. RETURN
  107. END
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  
  134.  

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